jacobiff.m

来自「数值分析各种公式的解法」· M 代码 · 共 15 行

function x=jacobiff(a,b,eps)
n=length(a);
d=diag(diag(a));
l=tril(a,-1);
u=triu(a,1);
x0=diag(zeros(n));
k=1
x=-inv(d)*(l+u)*x0+inv(d)*b;
while norm(x-x0)>eps
    k=k+1;
    p=x;
    x=-inv(d)*(l+u)*x+inv(d)*b
    x0=p;
end