雅可比迭代法MATLAB程序.docVIP

function [y,n]=jacobi(A,b,x0,ep) if nargin==3 ep=1.0e-6; else if nargin 3 error return end end D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,10); B=D\(L+U); f=D\b; y=B*x0+f; n=1; while norm(y-x0)=ep x0=y; y=B*x0+f; n=n+1; end