fixpt.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 24 行

function[p,k,err,P]=fixpt(f1021,p0,tol,max1)
%  g是给定的迭代函数。
%  p0是给定的初始值。
%  tol是给定的误差界。
%  max1是所允许的最大迭代次数。
%  k是所进行的迭代次数加1。
%  p是不动点的近似值。
%  err是误差。
%  P={p1,p2,…,pn}
P(1)=p0;
for k=2:max1
    P(k)=feval('f1021',P(k-1));
        k,err=abs(P(k)-P(k-1))
  p=P(k);
    if(err<tol),
           break;
end
    if k==max1
           disp('maximum number of iterations exceeded');
end

end
P=P