jacobi.m

来自「本文件提供了几个数值分析中的Matlab小程序。水平有限」· M 代码 · 共 22 行

function[x,k]=jacobi(a,b,x0,e)
n=length(a);
for i=1:n;
    v(i)=1/a(i,i);
end
    d=diag(v);
    bj=eye(n)-d*a;
    g=d*b;
    x=bj*x0+g;
    k=1;
    while norm(x-x0,'inf')/norm(x,'inf')>e
        x0=x;
        k=k+1;
        x=bj*x0+g;
        if k==1000
            break;
        end 
    end
    f=norm(x-x0,'inf')/norm(x,'inf')