jacobi.m

来自「该程序包包含大量的matlab程序」· 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