雅克比迭代法求解线性方程组的程序

13623121332812832215xxxxxxxxx雅克比迭代法求解线性方程组的程序：#includestdio.h#includemath.h#includeconio.h#defineN50#defineeps1e-8voidmain(){doublex[3]={0,0,0},y[3],a[3][3]={5,2,1,2,8,-3,1,-3,-6},b[3]={8,21,1},s1,s2;inti,k=0,j;do{s2=0;k=k+1;for(i=0;i3;i++){s1=0;for(j=0;j3;j++){if(i==j)continue;s1+=a[i][j]*x[j];}y[i]=(b[i]-s1)/a[i][i];s2+=fabs(x[i]-y[i]);}if(s2eps)printf(迭代次数K=%d\n,k);if(kN){printf(不收敛！);exit();}for(i=0;i3;i++)x[i]=y[i];}while(s2eps);for(i=0;i3;i++)printf(x[%d]=%lf\n,i,x[i]);}