#includevoid main()
{
float array[4][5], x, s; //定义二维数组和浮点变量//
int i, j, m, k; //定义整型变量//
printf("please input coefficient matrix array:\n"); //输入增广矩阵array//
for(i=0; i<4; i++)
{
for(j=0; j<5; j++)
{
scanf("%f", &array[i][j]); //输入一个4*5阶矩阵//
}
printf("\n");
}
for(i=0; i<4; i++) //输出输入的矩阵//
{
for(j=0; j<5; j++)
{
printf("%8.4f",array[i][j]);
}
printf("\n");
}
printf("\n\n\n\n\n\n");
for(k=0;k<=3;k++) //在相应的列中选主元//
{
m = k;
for(i=k+1;i<=3;i++)
{
x = fabs(array[m][k]);
if(fabs(array[i][k] > x))
m = i;
}
printf("%d\n\n\n\n",m);
if(array[m][k]==0)
{
printf("ERROR");
return;
}
else //两行进行比较交换两行//
{
for(j=k;j<=4;j++)
{
x = array[k][j];
array[k][j] = array[m][j];
array[m][j] = x;
}
for(i=0;i<4;i++)
{
for(j=0;j<5;j++)
{
printf("%8.4f ",array[i][j]);
}
printf("\n");
}
printf("\n\n\n\n\n\n");
for(i=k+1;i<=3;i++)
{
s = array[i][k]/array[k][k]; //消元过程//
for(j=0;j<=4;j++)
{
array[i][j] = array[i][j] - s * array[k][j];
}
}
}
}
for(i=0;i<4;i++) //输出变换后的矩阵//
{
for(j=0;j<5;j++)
{
printf("%8.4f", array[i][j]);
}
printf("\n");
}
array[3][4] = array[3][4]/array[3][3]; //回代过程//
for(i=2;i>=0;i--)
{
x = 0;
for(j=i+1;j<4;j++)
{
x = x + array[i][j] * array[j][4];
}
array[i][4]=(array[i][4] - x) / array[i][i];
}
printf("x1 = %f", array[0][4]); //输出结果//
printf("\n");
printf("x2 = %f", array[1][4]);
printf("\n");
printf("x3 = %f", array[2][4]);
printf("\n");
printf("x4 = %f", array[3][4]);
printf("\n");
}
验证程序的正确性,输入已知解的方程增广矩阵,运行结果与方程解一致。
图为输入的增广矩阵
图为运行结果方程的解
图为随机输入的增广矩阵运行的结果