雅可比迭代法解线性方程组(共1页).doc

《雅可比迭代法解线性方程组(共1页).doc》由会员分享，可在线阅读，更多相关《雅可比迭代法解线性方程组(共1页).doc（1页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

精选优质文档-倾情为你奉上雅可比迭代法解线性方程组#includestdio.h #includemath.h#includealloc.hJacobidiedai(n,a,b,x) int n; double *a,*b,*x; int i,j; double *x0,m=0,eps; x0=malloc(n*sizeof(double); for(i=0;in;i+) x0i=xi; while(1) for(i=0;in;i+) for(j=0;jn;j+) m+=ai*n+j*x0j; m-=ai*n+i*x0i; xi=(bi-m)/ai*n+i; m=0; for(i=1;in;i+) eps=fabs(x0-x00); if(epsfabs(xi-x0i) eps=fabs(xi-x0i); if(eps1e-3) return 1; for(i=0;in;i+) x0i=xi; free(x0);main() int i; double a33=5,2,1,-1,4,2,2,-3,10,b3=-12,20,3,x3=0,0,0; Jacobidiedai(3,a,b,x); for(i=0;i3;i+) printf(x%d=%fn,i,xi);专心-专注-专业