gseid.m

来自「matlab下雅克比在数值计算中有很好的应用!」· M 代码 · 共 22 行

function x=gseid(A,b,p,delta,max1)
N=length(b);
for k=1:max1
for j=1:N
    if j==1
        x(1)=(b(1)-A(1,2:N)*p(2:N))/A(1,1);
    elseif j==N
        x(N)=(b(N)-A(N,1:N-1)*(x(1:N-1))')/A(N,N);
    else 
        x(j)=(b(j)-A(j,1:j-1)*(x(1:j-1))'-A(j,j+1:N)*p(j+1:N))/A(j,j);
    end
end
err=abs(norm(x'-p));
p=x';
if (err<delta)
    break
end
end
x=x';
err,k