function [x,n,err,vx]=secanti(fun,x0,x1,toll,nmax)
  %% [x,n,err]=secanti(fun,x0,x1,toll,nmax)
  n=0;				
  err='';
  x=x0;
  fx0=eval(fun);
  x=x1;
  vx=[x0,x];
  do
    fx=eval(fun);
    if fx-fx0 == 0
      err='Divisione per zero'
      return;
    end
    inc=(x-x0)*fx/(fx-fx0);
    x0=x;
    fx0=fx;
    x=x-inc;
    vx=[vx,x];
    n=n+1;
  until abs(inc) < toll | n==nmax
  if n==nmax
    err='No converge';
  end