c编的sor迭代法解线性方程组的程序(3页).doc

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

1、-c编的sor迭代法解线性方程组的程序-第 3 页c编的sor迭代法解线性方程组的程序2010-12-15 20:33#include #include double norm(double *x,double *y,int n) int i=0; double s=0; for(i=0;in;i+) s=s+fabs(xi-yi)*fabs(xi-yi); return sqrt(s); void Jacobi(double A,double b,double x0,double x,double w,int n1)int n=100;/最大迭代次数double ee=0.0001,s=0;

2、/精度int i,j,k;for(i=1;i=n;i+) for(j=0;jn1;j+) xj=x0j; for(j=0;jn1;j+) s=0; for(k=0;kn1;k+) s=s+Aj*n1+k*x0k; x0j=xj+w*(bj-s)/Aj*n1+j; if (norm(x,x0,n1)ee) break; int main() double a33=8,-1,1,2,10,-1,1,1,-5;double b3= 1,4,3;double x03= 0,0,0.0;double x3= 0,0,0,w=1.005;int n1=3,i;Jacobi(a,b,x0,x,w,n1);for (i=0;in1; i+)printf(x(%d)=%6.4fn,i+1,xi);system(pause); return(0);