수다닷컴

질문 제목 : 가우스-조단 소거법으로 연립방정식 해 구하기질문 요약 :디버그 에러 수정질문 내용 :컴파일 에러는 없으나 실행 후 디버그 에러가 납니다.
#include stdio.h
#include stdlib.h
#include time.h
#include math.h
#define max 1000
void printout(double**, double [], int);
void gauss(double**, double [], int);
void back(double**, double [], int);
int main(void)
{
double **a, *b;
int i, j, n;
srand(time(null));
printf(nxn 행렬의 n을 입력하세요 : );
scanf(%d, &n);
if (n = 2 && n = max)
{
a = (double**) malloc(sizeof(double)*n);
for (i = 0; i n; i++)
a[i] = (double*) malloc(sizeof(double)*n);
b = (double*) malloc(sizeof(double)*n);
for (i = 0; i n; i++)
{
for (j = 0; j n; j++)
{
a[i][j] = rand() % (21 -10) / 1.0;
}
}
for (i = 0; i n; i++)
b[i] = rand() % (21 -10) / 1.0;
printout(a, b, n);
gauss(a, b, n);
for (i = 0; i n; i++)
free(a[i]);
free(a);
free(b);
}
return 0;
}
void printout(double **a, double b[], int n)
{
int i, j;
for (i = 0; i n; i++)
{
for (j = 0; j n; j++)
{
printf(%+.2lf\t, a[i][j]);
}
printf( x%d , i+1);
printf( %+.2lf\n, b[i]);
}
printf(\n);
}
void gauss(double **a, double b[], int n)
{
int i, j, k, l;
enum
{
false, true
} error;
double mult, temp;
printf(------------------------------------------------\n);
error = false;
k = 0;
while (k n && error != true)
{
l = k;
for (j = k+1; j n; j++)
{
if (fabs(a[j][k]) fabs(a[l][k]))
{
l = j;
}
for (j = k; j n; j++)
{
temp = a[k][j];
a[k][j] = a[l][j];
a[l][j] = temp;
}
temp = b[k];
b[k] = b[l];
b[l] = temp;
if (a[k][k] != 0)
{
for (j = k+1; j n; j++)
{
mult = -1 * (a[j][k] / a[k][k]);
for (i = k; i n; i++)
{
if (mult != 0)
a[j][i] += (mult * a[k][i]);
}
b[j] += (mult * b[k]);
printout(a, b, n);
}
}
else
error = true;
}
k += 1;
}
if (error == true)
printf(오류입니다.\n);
printout(a, b, n);
printf(------------------------------------------------\n);
back(a, b, n);
}
void back(double **a, double b[], int n)
{
int i, j, k;
double sum;
double x[] = {0.0};
printf(a[n][n] = %.2lf\n\n, a[n-1][n-1]);
if (fabs(a[n-1][n-1]) == 0 || fabs(a[n-1][n-1]) 1/1000000)
{
printf(비가역행렬입니다.\n);
}
else
{
x[n-1] = b[n-1] / a[n-1][n-1];
for (i = n-2; i -1; i--)
&n
{
sum = 0.0;
for (j = n-1; j i; j--)
{
sum += (a[i][j] * x[j]);
}
if (fabs(a[n-1][n-1]) == 0 || fabs(a[n-1][n-1]) 1/1000000)
{
printf(비가역행렬입니다.\n);
}
else
{
x[i] = (b[i] - sum) / a[i][i];
}
}
printf(선형시스템의 해\n);
for (k = 0; k n; k++)
{
printf(x%d = , k+1);
printf(%.2lf\t, x[k]);
}
printf(\n\n);
}
}

신청하기