您的当前位置:首页用高斯用列主元消元法求解下面的方程组

用高斯用列主元消元法求解下面的方程组

来源:锐游网


用高斯用列主元消元法求解下面的方程组

x1x2x34x425x4x3x12x412342x1x2x311x432x1x27x3x40

#include

int main(void)

{

float A[4][4]={{1,-1,1,4 },

{5,-4,3,2 },

{2,1,1,11},

{2,-1,7,-1},};

float b[4]={2,4,3,0};

float x[4]={0};

float Aik,S;

int i,j,k;

int size=4;

//消去过程

for(k=0;k{

if(!A[k][k])

return -1;

for(i=k+1;i{

Aik=A[i][k]/A[k][k];

for(j=k;j{

A[i][j]=A[i][j]-Aik*A[k][j];

}

b[i]=b[i]-Aik*b[k];

}

}

//消去的结果

for(i=0;i//回代过程

x[size-1]=b[size-1]/A[size-1][size-1];

for(k=size-2;k>=0;k--)

{

S=b[k];

for(j=k+1;j{

S=S-A[k][j]*x[j];

}

x[k]=S/A[k][k];

}

//solution

printf(\"The solution x[]=\\n\");

for(i=0;iprintf(\"%f \

return 0;

}

因篇幅问题不能全部显示,请点此查看更多更全内容

Top