2 yudacanhe123 yudacanhe123 于 2016.03.15 17:37 提问

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

/*矩阵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;
}

1个回答

caozhy
caozhy   Ds   Rxr 2016.03.15 18:04
Csdn user default icon
上传中...
上传图片
插入图片
准确详细的回答,更有利于被提问者采纳,从而获得C币。复制、灌水、广告等回答会被删除,是时候展现真正的技术了!