steffensen.m

这是非线性方程（组）的数值解法

function x=Steffensen(f,x0,eps)
% Steffensen迭代方法，为加速迭代收敛的方法
% 比Picard的收敛要快
if nargin==2
    eps=1e-5;
end
y1=f(x0);
z1=f(y1);
x1=x0-((y1-x0)^2)/(z1-2*y1+x0);
while abs((x1-x0))>eps
    x0=x1;
    y1=f(x0);
    z1=f(y1);
    x1=x0-((y1-x0)^2)/(z1-2*y1+x0);
end
x=x1;
end