수다닷컴

8x1 + x2 -x3 =82x1 + x2 +9x3 =12x1 - 7x2 +2x3 =-4이식을 가우스 자델 반복법으로 구해야되는데요밑에는 프로그램샘플입니다.. 여기서 멀 고쳐야 위에식을 돌릴수잇을까요?ㅠㅜ파일도 올려놓겟습니다.#include stdio.h
#include stdlib.h
#include math.h#define MAX_EQ20float calc_e(float[], float[], int);print_equation(float a[][MAX_EQ+1], float b[], int n)
{
int i, j;for(i=1;i=n;i++) {
for(j=1;j=n;j++) {
printf(%8.3f*x%d, a[i][j], j);
if(j == n) printf( = );
else printf( + );
}
printf(%-8.3f\n, b[i]);
}
}print_solution(int p, float x[], int n)
{
int i;
printf([Iteration %02d] , p);
for(i=1;i=n;i++) {
printf(x%d = %-10.5f, i, (float)x[i]);
if(i == n) printf(\n);
else printf(, );
}
}read_input(char *file, float a[][MAX_EQ+1], float b[], float x[], int*n)
{
int i, j;
FILE *fp;if((fp = fopen(file, r)) == (FILE *)0) return -1;fscanf(fp, %d, (char *)n);
for(i=1;i=(int)*n;i++) {
for(j=1;j=(int)*n;j++) {
fscanf(fp, %f, (char *)&a[i][j]);
}
fscanf(fp, %f, (char *)&b[i]);
}for(i=1;i=(int)*n;i++)
fscanf(fp, %f, (char *)&x[i]);return 1;
}float calc_e(float xo[], float xe[], int n)
{
int i;
float s;for(i=1, s=0.0; i = n; i++)
s += fabs(xo[i]-xe[i]);return( (float)(s/(float)n) );
}main(int argc, char *argv[])
{
int i, j, p, n;
float a[MAX_EQ+1][MAX_EQ+1], b[MAX_EQ+1];
float xo[MAX_EQ+1], xe[MAX_EQ+1], e;if(argc 3) {
printf(Usage: %s input-file e\n, argv[0]);
exit(0);
}
if(read_input(argv[1], a, b, xe, &n) 0) {
printf(cannot read the input file: %s\n, argv[0]);
exit(-1);
}e = (float)atof(argv[2]);print_equation(a, b, n);
printf(epsilon = %8.3f\n\n, e);print_solution(0, xe, n);for(p=1;1;p++) {
for(i=1;i=n;i++) {
if(p%2) {
xo[i] = 0.0;
for(j=1;j=n;j++) {
if(j = (i-1))
xo[i] -= (a[i][j]*xo[j]);
else if(j = (i+1))
xo[i] -= (a[i][j]*xe[j]);
else
xo[i] += b[i];
}
xo[i] /= a[i][i];
} else {

xe[i] = 0.0;
for(j=1;j=n;j++) {
if(j = (i-1))
xe[i] -= (a[i][j]*xe[j]);
else if(j = (i+1))
xe[i] -= (a[i][j]*xo[j]);
else
xe[i] += (b[i]);
}
xe[i] /= a[i][i];
}if(p%2)
print_solution(p, xo, n);
else
print_solution(p, xe, n);

if(calc_e(xo, xe, n) e) break;
}
}
}

신청하기