picard.m

来自「这是非线性方程（组）的数值解法」· M 代码 · 共 14 行

function x=Picard(f,x0,eps)
% f(x)=x-sin(x)-1/2
% x=F(x)=sin(x)+1/2 ,如果迭代构造不好可能就发散了
% f就是F(x),x0就是初值
if nargin==2
    eps=1e-5;
end
x1=f(x0);
while abs(x1-x0)>eps && abs(f(x1))>eps
    x0=x1;
    x1=f(x0);
end
x=x1;
end