Xn=0;%Xn+1 X=10;%Xn Xo=9;%Xn-1 F=X.^5;%F(Xn) Fo=Xo.^5;%F(Xn-1) a=0;%real root for i=0:1:200 Xn = X-F*((X-Xo)/(F-Fo)); i; Xn; D = Xn-X; Xo=X; X=Xn; F = X.^5;%F(Xn) but we use X because X = Xn Fo=Xo.^5;%F(X) but we use Xo because Xo = X if abs(a-F)<.00001 [i,Xn,F,D] break end end