%
% Il metodo delle secanti
%
format long
f=@(x)x^2-1   
tol=10^(-10)
x0=3
f0=f(x0)
x1=2
f1=f(x1) % x0 e x1 sono le "stime" iniziali
n=0;
while(abs(f1)>tol & n< 20)
    n=n+1;
    xn=x0-f0*(x1-x0)/(f1-f0)
    x0=x1; f0=f1;
    x1=xn; f1=f(x1);
end
[n,x1]