#include#include#defineN3/*方程阶数*/voidmyprintf(doublea[][N],double*b);/*打印运算矩阵a,b*/voidhdai(doublea[][N],double*b);/*回代求解*/voidlzyxq(doublea[][N],double*b);/*高斯变换*/voidlzy(doublea[][N],double*b,intk);/*行变换*/voidmain(){doublea[N][N]={3,2,1,5,6,4,7,8,9};doubleb[N]={5,6,7};lzyxq(a,b);myprintf(a,b);hdai(a,b);}voidmyprintf(doublea[][N],double*b){inti,j;printf("ais:\n");for(i=0;i=0;i--){sum=0;for(j=i+1;j<N;j++)sum+=a[i][j]*x[j];x[i]=(b[i]-sum)/a[i][i];}printf("resout:\n");憨钉封固莩改凤爽脯鲸for(i=0;i<N;i++){printf("x[%d]=%4.3f\n",i+1,x[i]);}}--------------------------------------------以上程序由vc6.0编译通过N为线性方程组阶数
阅读全文 >