format long
g=@(x)(x^3-2*x+1)
f=@(x)x^3-3*x+1
tol=10^(-5)
xv=.5 % x iniziale
maxiter=50
for n=1:maxiter 
    xn=g(xv);
    [n, xn]
    pause
    if abs(xn-xv)<tol, break, else xv=xn; end
end
if n>=maxiter, disp('Troppe iterazioni'), end