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

《雅克比迭代法求解线性方程组的程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、雅克比迭代法求解线性方程组的程序：#include"stdio.h"#include"math.h"#include"conio.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;i<3;i++){s1=0;for(j=0;j<3;j++){if(i==j)continue;s1+=a[i][j]*x[j];}y[i]=(b[i]-s1)/a[i][i];

2、s2+=fabs(x[i]-y[i]);}if(s2N){printf("不收敛！");exit();}for(i=0;i<3;i++)x[i]=y[i];}while(s2>eps);for(i=0;i<3;i++)printf("x[%d]=%lf",i,x[i]);}