gausseidel.m

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

function x=gausseidel(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+l)*b
while norm(x-x0)>eps
    k=k+1;
    p=x;
    x=-inv(d+l)*u*x+inv(d+l)*b
    x0=p;
end