seidel.m

求解非线性不动点方程组x=G(x)的非线性塞德尔迭代方法

function x0=seidel(x0)
%求解非线性不动点方程组x=G(x)，给定初值x0
%非线性塞德尔迭代

max1=1000;
eps=5e-6;
N=length(x0);

for k=1:max1
    x1=x0;
    for j=1:N
        A=G(x1);
        x1(j)=A(j);
    end
    err=abs(norm(x1-x0));
    x0=x1;
    if err<eps
       break;
    end
end        
ss=x1