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

精选优质文档-倾情为你奉上雅可比迭代法解线性方程组#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);