雅可比迭代法的MATLAB程序(共2页).doc

精选优质文档-倾情为你奉上雅可比迭代法的MATLAB程序：Functionx,k,index=Jacobi(A,b,ep,it-max)% 求线性方程组的雅可比法;% A为方程组的系数矩阵;% b为方程组的右端项;% x为方程组的解;% ep为精度要求,缺省值为le-5;% it_max为最大迭代次数,缺省值为100;% k为迭代次数;% index 为指标变量,index=0表示计算失败,index=1表示计算成功;if nargin4 it_max=100;endif nargin3ep=le-5;endn=length(A);k=0;x=zeros(n,1);y=zeros(n,1);index=1;while k=it_max for i=1:n if abs (A(i,i)le-10 index=0; return;endy(i)=(b(i)-A(i,1:n)*x(1:n)+A(i,i)*x(i)/A(i,i