jacobi.m

来自「这是一些matlab函数：欧拉公式一、格斯比函数等等 它是一些已知的m文件 可以」· M 代码 · 共 28 行

function tx=jacobi(A,b,imax,x0,tol)
del=10^-10;
tx=[x0]; n=length(x0);
for i=1:n
    dg=A(i,i);
    if abs(dg)<del
        disp('diagonal element is too small');
        return
    end
end
for k=1:imax
    for i=1:n
        sm=b(i);
        for j=1:n
            if j~=i
                sm=sm-A(i,j)*x0(j);
            end
        end
        x(i) =sm/A(i,i);
    end
    tx=[tx;x];
    if norm(x-x0)<tol
        return
    else
        x0=x;
    end
end