GaussElimination_MCP

#include<stdio.h>
#include<math.h>
int main()
{
int i,j,k,p,n;
float m,temp,temp1,A[10][10],X[10];
printf("请输入方程组未知个数n:");
scanf("%d",&n);
printf("请输入增广矩阵A相关信息：\n");
for(i=1;i<=n;++i)
for(j=1;j<=n+1;++j)
scanf("%f",&A[i][j]);
for(i=1;i<=n-1;++i)//化成上三角型式
{

temp=fabs(A[i][i])*1.0;//最大元
p=i;
for(j=i+1;j<=n;++j)//找出最大元
{
if(fabs(A[j][i])>temp)
{
temp=A[j][i];
p=j;
}
}
printf("%f\n",temp);
if(temp==0)
{
printf("数据输入有误或不满足Cauuss列主元法的规则!\n");
return 0;
}
if(p!=i)//如果最大元不是i行i列，则交换i，p两行的数据
{
for(j=i;j<=n+1;++j)
{


temp=A[i][j];
A[i][j]=A[p][j];
A[p][j]=temp;
}
}
for(j=i+1;j<=n;++j)//消去和列主元同列的数据，使其变为0
{
m=A[j][i]/A[i][i];
for(k=i;k<=n+1;++k)//列
{
A[j][k]=A[j][k]-A[i][k]*m;
}
}
for(k=1;k<=n;++k)
{
    for(j=1;j<=n+1;++j)\
    printf("%f ",A[k][j]);
printf("\n");
}
}


if(A[n][n]==0)//判断是否满足列主元是否为零
{
printf("数据输入有误或不满足Cauuss列主元法的规则!\n");
return 0;
}
X[n]=A[n][n+1]/A[n][n];
for(i=n-1;i>=1;--i)
{
temp=0;
for(j=i+1;j<=n;++j)//累加计算
{
temp=temp+X[j]*A[i][j];
}
X[i]=(A[i][n+1]-temp)/A[i][i];
}
printf("该方程的解为：\n");
for(i=1;i<=n;i++)
{
printf("%f ",X[i]);
}
printf("\n");
return 0;


}