gseidel.m

来自「这是一些matlab函数：欧拉公式一、格斯比函数等等 它是一些已知的m文件 可以」· M 代码 · 共 29 行

function tx=gseidel(A,b,imax,x0,tol)
del=10^-10;
tx=[x0];n=length(x0);
for i=1:n
    dg=A(i,i);
    if abs(dg)<del
        disp('diagonal element is too small');
        return
    end
end
for k=1:imax
    x=x0;
    for i=1:n
        sm=b(i);
        for j=1:n
            if j~=i
                sm=sm-A(i,j)*x(j);
            end
        end
        x(i)=sm/A(i,i);
    end
    tx=[tx,x];
    if norm(x-x0)<tol
        return
    else
        x0=x;
    end
end