수다닷컴

질문 요약 :가우스-조단 소거법에서 좀 더 큰 행렬을 받고 싶어요.
질문 내용 : 가우스-조단 소거법으로 REF랑 RREF를 만드는데요.
행렬이 4*4를 벗어난 큰 행렬은 값이 아예 안나오거든요.
다른 5*5 라든지 더 큰 행렬의 REF와 RREF가 나오게 하려면 여기서 어떻게 손을 봐야할까요??
머리가 안 돌아가네요ㅜㅜ 도와주세요ㅜㅜㅜㅜ#include stdio.h
#include stdlib.h
#include math.h
#define MAX 5
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); // 연산 결과 출력함수
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);
for(i=0; i m; i++)
{
A-a[i] = (double *)malloc( sizeof(double) * n);
}
}
void destroy(matrix_t *A)
{
int i;
for(i=0; i A-m; i++)
{
free(A-a[i]);
}
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)
{
int i, j;
for(i = 0; i A-m; i++)
{
for(j = 0; j A-n; j++)
{
printf( %.3f\t, A-a[i][j]);
}
puts();
}
puts();
}
void gauss(matrix_t *A)
{
int i, j, k;
double pivot;
printf(\n 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;
print(A);
for(k=0; k A-m; k++)
{
if(i == k) continue;
pivot = A-a[k][i];
for(j=i+1; j A-n; j++)
{
A-a[k][j] -= A-a[i][j] * pivot;
}
A-a[k][i] = 0;
}
print(A);
}
return;
}
int main()
{
int i = 0, j = 0;
int n;
matrix_t A;
FILE *fptr;
fptr = fopen(matrix.txt,rt);
load(&A,fptr);
print(&A);
gauss(&A);
for(i=0; i A.m; i++)
{
printf(X%d = %f\n,i+1,A.a[i][A.n-1]);
}
destroy(&A);
fclose(fptr);
}

신청하기