jacobi.m

来自「都是大家需要的程序,大家要下的快下,大家要多多支持」· M 代码 · 共 12 行

function y=jacobi(a,b,x0)
D=diag(diag(a));
U=-triu(a,1);
L=-tril(a,-1);
B=D\(L+U);
f=D\b;
y=B*x0+f;n=1;
while norm(y-x0)>=1.0e-6
    x0=y;
    y=B*x0+f;n=n+1;
end
n