雅克比迭代法.doc

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

1、雅克比迭代法上机题目：用雅克比迭代法解线性方程组上机程序：#include#includemain(){doubleA[3][3]={{5,2,1},{-1,4,2},{2,-3,10}},b[3]={-12,20,3};//输入系数矩阵A和右端向量bdoublen=3,tol=1.0e-3,x[3]={-3,1,1};//输入方程大小n，误差限tol，和初始向量xdoubley[3];//记录每次迭代产生的新的近似解doublemaxerr;//maxerr记录相邻两个数值解想x与y的差的绝对值

2、的最大的那个分量。intk,i,j;printf("此线性方程组的增广矩阵如下：");//下面循环是输出增广矩阵（A,b）for(k=0;k

3、x[i]);}printf(")");//输出初始迭代向量完printf("误差限是%f",tol);//输出误差限printf("Jacobi迭代解序列X^(k)max

4、x^(k+1)-x^(k)

5、");printf("x^%d=",k=0);for(i=0;i

6、x[j];}y[i]=(b[i]-T)/A[i][i];}//Jacobi迭代完//求相邻两个数值解想x与y的差的绝对值的最大的那个分量maxerr=fabs(y[0]-x[0]);for(j=1;j

7、[i]);printf("%f",maxerr);printf("");}while(maxerr>tol);}运行结果：