수다닷컴

질문내용 :

가우스-조단 소거법으로 구하는 프로그램 짜고있는데요.
5행 7열인 행렬을 하면 답이 안나와요ㅠㅠ

제가넣은 input값
5 7
4 5 36 0 12 40 119
4 10 36 -12 6 44 94
8 10 72 4 18 92 226
-4 -5 -36 0 -6 -40 -95
-12 -15 -108 -8 -24 -144 -333잘못된 아웃풋... 맨 마지막 문단이 최종 값으로 나온거에요..

이게정답

1, 0, 9, 0, 0, 0, 9;
0, 1, 0, 0, 0, 8, 7;
0, 0, 0, 1, 0, 3, 3;
0, 0, 0, 0, 1, 0, 4;
0, 0, 0, 0, 0, 0, 0;

위처럼 나오려면 어딜 어떻게 수정해야할까요?ㅠㅠㅠㅠ
질문 내용 :

#include stdio.h
#include stdlib.h
#include math.h

typedef struct matrix_t
{
int m, n;
double **a;
}
matrix_t;

void create(matrix_t *A, int m, int n);
void destroy(matrix_t *A);
void print(matrix_t *A, const char *title); // 연산 결과 출력함수
void gauss(matrix_t *A); // 가우스 요르단

void create(matrix_t *A, int m, int n)
{
int i;

A-m = m;
A-n = n;
A-a = (double **)malloc( sizeof(double*) * m);
A-a[0] = (double *)malloc( sizeof(double) * m * n);
for(i=1; i m; i++)
{
A-a[i] = A-a[i-1] + n;
}
}

void destroy(matrix_t *A)
{
int i;

if(A-a) {
if(A-a[0]) {
free(A-a[0]);
}
free(A-a);
}
A-m = A-n = 0;
A-a = 0;
}

void load(matrix_t *A, FILE *fptr)
{
int m, n;
int i, j;

fscanf(fptr,%d %d,&m,&n);
create(A,m,n);
for(i = 0; i A-m; i++)
{
for(j = 0; j A-n; j++)
{
fscanf(fptr,%lf, &A-a[i][j]);
}
}
}

void print(matrix_t *A, const char *title)
{
int i, j;

if(title) puts(title);
for(i = 0; i A-m; i++)
{
for(j = 0; j A-n; j++)
{
printf( % 6.2f\t, A-a[i][j]);
}
puts();
}
puts();
}

void gauss(matrix_t *A)
{
int i, j, k;
double pivot;
// char title[10];

printf(***************************************************** \n);
printf(Gauss-Jordan Elimination \n);
printf(----------------------------------------------------- \n\n);

for(i=0; i A-m; i++)
{
pivot = A-a[i][i];
if(pivot == 0)
{
return;
}
for(j=i+1; j A-n; j++)
{
A-a[i][j] /= pivot;
}
A-a[i][i] = 1;

// sprintf(title,%d = ,i);
print(A,0);

for(k=0; k A-m; k++)
{
if(i == k) continue;
pivot = A-a[k][i];
 sp; for(j=i+1; j A-n; j++)
{
A-a[k][j] -= A-a[i][j] * pivot;
}
A-a[k][i] = 0;
}
print(A,0);
}
printf(***************************************************** \n\n);
return;
}

int main()
{
int i;
matrix_t A = {0,};
FILE *fptr;

fptr = fopen(input2.txt,rt);
load(&A,fptr);
fclose(fptr);

// print(&A,A = );
print(&A,\nOriginal expression \n);
gauss(&A);

puts(Solution );
printf(\n);
for(i=0; i A.m; i++)
{
printf( X%d = %8.2f\n,i+1,A.a[i][A.m]);
}

destroy(&A);
}

신청하기