高斯列主元消元法的VC语言代码

【这是自己编写的程序，请自己阅读着使用。

】
/*shuzhi1*/
#include <stdio.h>
#define C 3
void main()
{
float A[C][C];
float b[C];
float X[C];
float t,r,s;
int i,j,k,m,n;
printf("Please input A:\n");
for(i=0;i<=C-1;i++)
{
for(j=0;j<=C-1;j++)
scanf("%f",&A[i][j]);
}
printf("The data in the computer is:\nA=\n");
for(i=0;i<=C-1;i++)
{
for(j=0;j<=C-1;j++)
printf("%4.4f ",A[i][j]);
printf("\n");
} /*A数据的输入与检验*/
printf("Is there any mistake? Press 1 for yes,Press 0 for no! ");
scanf("%d",&k);
for(;k==1;)
{
printf("Please input the adress(i,j) of the error data: ");
scanf("%d",&i);
scanf("%d",&j);
printf("Please input the right data: ");
scanf("%f",&t);
A[i-1][j-1]=t;
printf("The data in the computer is:\nA=\n");
for(i=0;i<=C-1;i++)
{
for(j=0;j<=C-1;j++)
printf("%4.4f ",A[i][j]);
printf("\n");
}
printf("Is there any mistake? Press 1 for yes,Press 0 for no! ");
scanf("%d",&k);
} /*修改出错数据*/
printf("Please input b:\n");
for(i=0;i<=C-1;i++)
scanf("%f",&b[i]);
printf("The data in the computer is:\nb=\n");
for(i=0;i<=C-1;i++)
printf("%4.4f ",b[i]); /*b数据的输入与检验*/
printf("\nIs there any mistake? Press 1 for yes,Press 0 for no! ");
scanf("%d",&k);
for(;k==1;)
{
printf("Please input the adress(i) of the error data: ");
scanf("%d",&i);
printf("Please input the right data: ");
scanf("%f",&t);
b[i-1]=t;
for(i=0;i<=C-1;i++)
printf("%4.4f ",b[i]);
printf("\nIs there any mistake? Press 1 for yes,Press 0 for no! ");
scanf("%d",&k);
} /*修改出错数据*/
/*选择主元*/
for(i=0;i<=C-2;i++)
{
for(j=i+1;j<=C-1;j++)
{
if(A[i][i]<=0)
r=-A[i][i];
else
r=A[i][i];
if(A[j][i]<=0)
s=-A[j][i];
else
s=A[j][i];
if (r<s)
{
for(k=0;k<=C-1;k++)
{
t=0;
t=A[j][k];
A[j][k]=A[i][k];
A[i][k]=t; /*交换A[i]与A[j]*/
}
t=0;
t=b[j];
b[j]=b[i];
b[i]=t; /*交换b*/
}
}
printf("After %d change,the the A and b are:\n\n",i+1);
for(m=0;m<=C-1;m++)
{
for(n=0;n<=C-1;n++)
printf("%4.4f ",A[m][n]);
printf("b[%d]=%4.4f ",m+1,b[m]);
printf("\n");
}
printf("\n"); /*打印变化后的数组*/
for(m=i+1;m<=C-1;m++)
{
t=(-A[m][i]);
for(n=0;n<=C-1;n++)
A[m][n]=A[m][n]+(A[i][n]/A[i][i])*t;
b[m]=b[m]+b[i]/A[i][i]*t;
}
for(m=0;m<=C-1;m++)
{
for(n=0;n<=C-1;n++)
printf("%4.4f ",A[m][n]);
printf("b[%d]=%4.4f ",m+1,b[m]);
printf("\n");
}
printf("\n");
} /*主元消零并打印当前数组*/
/*迭代求X[C]*/
for(i=0;i<=C-1;i++)
X[i]=0;
for(i=C-1;i>=0;i--)
{
t=0;
for(j=i+1;j<=C-1;j++)
t=t+A[i][j]*X[j];
X[i]=(b[i]-t)/A[i][i];
}
printf("The answer of 'Ax=b' is:\n");
for(i=0;i<=C-1;i++)
printf("X[%d]=%4.4f ",i,X[i]);
printf("\n");
/*程序结束*/
}。