计算方法——雅克比迭代法(C程序)

#includestdio.h#includemath.h#defineL1000#defineM3#defineN4voidmain(){inti,j,k,l=0;floatm,e,a[M][N],x[M],y[M],temp[M];printf(\n请输入方程组的增广矩阵:\n);for(i=0;i=M-1;i++)for(j=0;j=N-1;j++)scanf(%f,&a[i][j]);printf(\n请输入初始向量:\n);for(i=0;i=M-1;i++)scanf(%f,&x[i]);printf(\n请输入误差范围:\n);scanf(%f,&e);for(i=0;i=M-1;i++){m=1/a[i][i];for(j=0;j=N-1;j++)a[i][j]=a[i][j]*m;for(j=0;j=N-2;j++)a[i][j]=-a[i][j];a[i][i]=0;}printf(\n迭代过程如下:\n);for(k=1;k=L;k++){for(i=0;i=M-1;i++){temp[i]=x[i];y[i]=a[i][0]*x[0];for(j=1;j=N-2;j++)y[i]=y[i]+a[i][j]*x[j];y[i]=y[i]+a[i][N-1];}for(i=0;i=M-1;i++){x[i]=y[i];printf(%12f,x[i]);if(fabs(temp[i]-x[i])=e)l=l+1;}printf(\n);if(l==3)break;}printf(\n方程组的解为:\n);for(i=0;i=M-1;i++)printf(%12f,x[i]);printf(\n\n);}