diedai.rar_K._矩阵设初始值

#include<stdio.h> #include<math.h> #include<string.h> #define M 3 #define N 4 void main() { static float a[M][N]={{5,-1,-1,3}, {1,-10,20,11},{-10,1,-5,-14}}; float x[M]={0,0,0}; double r; float t,q,e=0.0001; int i=0,j,k,lm=100; for(i;i<M;i++) { for(j=0;j<N;j++) printf("%15f",a[i][j]); printf("\n"); } for(k=0;k<lm;k++) { r=0; for(i=0;i<M;i++) { q=0; t=x[i]; for(j=0;j<M;j++) { if(j!=i) q=q+a[i][j]*x[j]; x[i]=(a[i][N-1]-q)/a[i][i]; } if(fabs(x[i]-t)>r) r=fabs(x[i]-t); } if(r<e)break; } if(k==lm)printf("no"); else for(i=0;i<M;i++) printf("x(%d)=%15.7f\n",i,x[i]); }