jacobi.m

来自「关于Euler、jacobi、rungkutta4数值算法的matlab实现」· M 代码 · 共 15 行

function x=jacobi(A,b,x0,e,N)

if nargin<5,N=10;end
if nargin<4,e=1e-4;end
if nargin<3,x0=zeros(n,1);end
x=x0;
x0=x+2*e;
k=0;A1=tril(A);iA1=inv(A1);
while norm(x0-x,inf)>e&k<N,
    k=k+1;
    %x0=x;
    x=-iA1*(A-A1)*x+iA1*b;
    disp(x');
end
if k==N,warning('已达到迭代次数上限');end