C++用高斯消元法求解方程组

Answer 1

#include<stdio.h>

int main()

{    

int a[3][4]={{1,2,3,1},{2,7,5,6},{1,4,9,-3}};

int i,j;

//第一行求解，现在第一行第一个参数就为1，所以不用作处理，直接消其它行第一列参数

for (i=1;i<3;i++)

{

for (j=3;j>=0;j--)

{

a[i][j]=a[i][j]-a[i][0]*a[0][j];

}

}

//第二行求解，先将第二行参数归一化

a[1][3]=a[1][3]/a[1][1];

a[1][2]=a[1][2]/a[1][1];

a[1][1]=1;

//消其它行第二列参数

//先消第一行

for (j=3;j>0;j--)

{

a[0][j]=a[0][j]-a[0][1]*a[1][j];

}

//消第三行

for (j=3;j>0;j--)

{

a[2][j]=a[2][j]-a[2][1]*a[1][j];

}

//求解第三行，先将第三行第参数归一化

a[2][3]=a[2][3]/a[2][2];

a[2][2]=1;

//消其它行第三列参数

for (i=0;i<2;i++)

{

for (j=3;j>1;j--)

{

a[i][j]=a[i][j]-a[i][2]*a[2][j];

}

}

printf("方程组解为:\nx1=%d\nx2=%d\nx3=%d\n",a[0][3],a[1][3],a[2][3]);

return 0;

}

结果：

Answer 2