function x=Steffensen(fun,x0,ep)
%Steffensen法求方程之根
%求根函数：fun
%与根相近点：x0
%精度设置：ep
x=x0;
x1=x0;
while 1
x1=x;
y=fun(x)+x;
z=fun(y)+y;
x=x-(y-x)*(y-x)/(z-2*y+x);
if abs(x-x1)<ep
    break;
end
end