Xn=0;%Xn+1, the iterates root X=10;%Xn a=0;%actual value f=0;%Xn+1=f(Xn), the function of the iterates of the root, X = alpha for i=1:1:100 f=X.^(1/3); Xn = f; [i,Xn,f,a-Xn] X=Xn; end