有没有大神帮看下鄙人写的高斯完全主元素消元法代码哪里出了问题,真心感谢啊！

Question

/*矩阵A用于存放线性方程组的增广矩阵，向量X表示线性方程组的解
#include
#include
int main()
{
double m,p,A[10][10],X[10];
int n,i,j,k,q,l,c,w;
char a;
printf("请输入方程的阶数（小于等于8）：\n");
scanf("%d",&n);
for(i=1;i<=n;i++)//方便起见从A[1][1]开始存入数据
{
printf("请输入增广矩阵第%d行：\n",i);
for(j=1;j<=n+1;j++)
scanf("%lf",&A[i][j]);
}//for
for(i=1;i {
m=fabs(A[i][i]);j=i;l=i; /*m表示矩阵第i列中的最大值，j用来标记最大值所在的 行数,l用来标记最大值所在的列数*/
for(k=i;k for(c=i;c if(fabs(A[k][c])>m)
{m=fabs(A[k][c]);j=k;l=c;}
for(q=i;q<=n+1;q++)//交换第i行和第j行
{
p=A[i][q];
A[i][q]=A[j][q];
A[j][q]=p;
}//for
for(w=i;w<=n;w++)//交换第i列与第c列
{
p=A[w][i];
A[w][i]=A[w][c];
A[w][c]=p;
}
//*************************消去过程*********
for(k=i+1;k<=n;k++)
{
m=-A[k][i]/A[i][i];
for(j=i;j<=n+1;j++)
A[k][j]=A[k][j]+A[i][j]*m;
}//for
}//for
//*************回代过程*************
X[n]=A[n][n+1]/A[n][n];
for(i=n-1;i>=1;i--)
{
p=0.0;
for(j=i+1;j<=n;j++)
p=p+A[i][j]*X[j];
X[i]=(A[i][n+1]-p)/A[i][i];
}//for
printf("线性方程组的解为：\n");
for(i=1;i<=n;i++)
printf("x[%d]=%lf\n",i,X[i]);
a=getchar();
return 0;
}

Answer 1

http://download.csdn.net/detail/henyjone/6743183