mulstablepoint.m

来自「源代码包括插值」· M 代码 · 共 25 行

function [r,n]=mulStablePoint(F,x0,eps)
%非线性方程组：f
%初始解：a
%解的精度：eps
%求得的一组解：r
%迭代步数：n

if nargin==2
    eps=1.0e-6;
end

x0 = transpose(x0);
n=1;
tol=1;
while tol>eps
    r= subs(F,findsym(F),x0);                           %迭代公式
    tol=norm(r-x0);                    %注意矩阵的误差求法，norm为矩阵的欧几里德范数
    n=n+1;
    x0=r;
    if(n>100000)                        %迭代步数控制
        disp('迭代步数太多，可能不收敛！');
        return;
    end
end