steffensen.m

steffenson迭代法Matlab源码

function X=steffensen(f,x0,n)
i=0;
X=x0;
y=feval(f,X);
z=feval(f,y);
x0=X-(y-X).^2/(z-2.*y+X);
while abs(X-x0)>0.1^(n+1)
   X=x0;
   y=feval(f,X);
   z=feval(f,y);
   x0=X-(y-X).^2/(z-2.*y+X);
   i=i+1;
   if i>10000
       X='error - result divergence!';
       break
   end
end