解方程c语言程序,C语言程序解线性方程组

匿名用户

1级

2011-06-04 回答

我今天写了一个,参考一下吧:

//作品：多元一次方程组的计算

//作者：与你看日出

//日期：2009年4月25日 星期六

//说明：输出值只能是小数(最多六位)，如x=1.876546

//Han： 初始设的最多未知数的个数，运行程序后只能比它的个数小

//hang：计算中途中自己改变的未知数的个数

//JUZHEN：初始的示例矩阵

#include

#include

#define Han 200//(可自设)多元一次方程组有n行n+1列(多的一列是等号右边的值)，给出行数就能确定矩阵，

#define JUZHEN {1,1,1,1,5},{1,2,-1,4,-2},{2,-3,-1,-5,-2},{3,1,2,11,0}//示例一个

main()

{

int i,j,k,m,n,t,cf,hang=4;

float temp;

float AA[Han][Han+1]={JUZHEN};//定义所要计算的数组

do//判断是否重试

{

for(i=0;i

{

printf("\n");

for(j=0;j

{

printf("%g\t",AA[i][j]);

}

}

printf("\n");

printf("是否自己输入？是：1；否：0");

scanf("%d",&t);

if(t==1)//判断是否自己输入数组

{

printf("输入未知数的个数");

scanf("%d",&hang);

for(i=0;i

{

for(j=0;j

{

printf("第%d行第%d列的数为：",i+1,j+1);

scanf("%f",&AA[i][j]);

}

}

}

for(k=0;k

{

while(AA[k][k]==0)//如果第K行K列的那个数为0,则加和重组一行。

{

for(m=k+1;m

for(n=k;n

{

AA[k][n]+=AA[m][n];

}

}

for(i=k;i

{

temp=AA[i][k];

for(j=k;j

{

AA[i][j]/=temp;

}

}

for(i=k+1;i

{

for(j=0;j

{

AA[i][j]-=AA[k][j];

}

}

}

for(k=hang-2;k>=0;k--)//这个大循环将数组的右上角转化为0

{

for(i=k+1;i

{

AA[k][hang+1-1]-=AA[k][i]*AA[i][hang+1-1];

AA[k][i]=0;

}

}

for(i=0;i

{

printf("\n");

for(j=0;j

{

printf("%g\t",AA[i][j]);

}

}

printf("\n未知数的值为：\n");

for(i=0;i

{

printf("x(%d)=\t%g\n",i+1,AA[i][hang+1-1]);

}

printf("\n");

printf("是否再试一次？是：1；否：0");

scanf("%d",&cf);

}

while(cf==1);//判断是否重试

}