수다닷컴

질문 제목 : 선형 방정식 계산 프로그램을 짜는 것 좀 도와주세요질문 요약 :gauss와 main를 어떻게 바꿔야 할지 모르겠습니다.질문 내용 :

파일(*.txt)로 입력을 받고 풀이과정을 역대입법이 아닌 가우스-조단 소거법을 이용하여 선형방정식을 푸는건데

예를 들면

x1 - 3 x2 = -3
2 x1 + x2 = 8

와 같은 방정식을 아래와 같이 파일로 입력받고

2
1 3
2 1
-3
8

풀이과정을 출력할 때에는

1 -3 -3
2 1 8

2 18
1 -3 -3

21 8
2 -6 -6

2 1 8
0 7 14

2 1 8
0 1 2

2 0 6
0 1 2

1 0 3
0 1 2

와 같이 바꿔야 하는데 머리에 쥐날 듯 하여서...
도와주세요... 흑

#include stdio.h
#include stdlib.h
#include math.h
#define MAX 5
typedef struct _matrix{double a[MAX][MAX];}matrix;
typedef struct _column{double c[MAX];}column;
typedef enum{E_FALSE, E_TRUE}E_BOOL;
void printout(matrix a, column c, int n);// 연산 결과 출력함수
void gauss(matrix a, column c, int n);// 가우스 소거법 계산함수
void backsub(matrix a, column c, int n);// 역대입법 계산함수
void printout(matrix a, column c, int n)
{
int i, j;
for(i = 0; i n; i++)
{
for(j = 0; j n; j++)
{
printf( %.3f\t, a.a[i][j]);
}

printf( x%d , i+1);
printf( = %.3f \n, c.c[i]);
}
printf(\n);
}
void gauss(matrix a, column c, int n)
{
E_BOOL error;
int i = 0, j = 0, k = 0, l = 0;
double multi = 0, temp = 0;
printf(\n Gauss-Jordan Elimination 풀이 과정 : \n);
printf( ***************************************************** \n\n);
error = E_FALSE;
k = 0;

while(k n && error != E_TRUE)
{
l = k;
for(j = k+1; j n; j++)
{
if(fabs(a.a[j][k]) fabs(a.a[l][k]))
{
l = j;
}
for(j = k; j n; j++)
{
temp = a.a[k][j];
a.a[k][j] = a.a[l][j];
a.a[l][j] = temp;
}
temp = c.c[k];
c.c[k] = c.c[l];
c.c[l] = temp;
if(a.a[k][k] != 0)
{
for(j = k + 1; j n; j++)
{
multi = -1 * (a.a[j][k]) / a.a[k][k];

for(i = k; i n; i++)
{
if(multi != 0)
a.a[j][i] = a.a[j][i] + multi * a.a[k][i];
}
c.c[j] = c.c[j] + multi * c.c[k];
&;printout(a, c, n);
}
}
else
error = E_TRUE;
}
k = k + 1;
}
if(error == E_TRUE)
printf(Trap condition......);
printout(a, c, n);
printf( ***************************************************** \n\n);
backsub(a, c, n);
}
void backsub(matrix a, column c, int n)
{
int i = 0, j = 0, k = 0;
float sum;
column mat = {0};
printf( a[n, n] = %.3f\n\n, a.a[n-1][n-1]);
if((fabs(a.a[n-1][n-1]) == 0) || (fabs(a.a[n-1][n-1]) 1/1000000))
printf( This matrix is singular, does not have unique solution. \n\n);
else
{
mat.c[n-1] = c.c[n-1] / (a.a[n-1][n-1]);
for(i = n - 2; i -1; i--)
{
sum = 0.0f;
for(j = n - 1; j i; j--)
sum = sum + a.a[i][j] * mat.c[j];
if(fabs(a.a[n-1][n-1]) == 0 || fabs(a.a[n-1][n-1]) 1/1000000)
printf( This matrix is singular, does not have unique solution. \n\n);
else
mat.c[i] = (c.c[i] - sum) / a.a[i][i];
}
printf( 입력한 선형방정식의 해답 \n);
for(k = 0; k n; k++)
{
printf( x%d = , k+1);
printf(%.3f \t, mat.c[k]);
}
printf(\n\n);
}
}
void main()
{
int i = 0, j = 0;
int n;
matrix a;
column c;
FILE *fptr = stdin;
printf( ** 가우스 - 조단 소거법을 이용한 선형방정식 계산 프로그램 ** \n\n);
printf( 선형방정식의 최대 차수를 입력하세요 : );
scanf(%d, &n);
printf(\n);
printf( 선형방정식의 수식을 입력하세요. \n);

for(i = 0; i n; i++)
{
for(j = 0; j n; j++)
{
printf( %d 번째 선형방정식 x%d 의 값 : , i+1, j+1);
scanf(%lf, &a.a[i][j]);
}
}
printf(\n 선형방정식의 결과 값을 입력하세요. \n);

for(i = 0; i n; i++)
{
printf( %d 번째 선형방정식의 결과 값 : , i+1);
scanf(%lf, &c.c[i]);
}
printf(\n);
gauss(a, c, n);
}

신청하기