fn='exp(x)-3*x^2';
x0=-2
x1=0
x=x1;
n=0;
while abs(eval(fn)) > 1.e-6
  x=x0;f0=eval(fn);
  x=x1;f1=eval(fn);
  x2=x1-f1*(x1-x0)/(f1-f0);
  [x0,x1,x2]
  x0=x1;x1=x2;
  x=x1;
  n=n+1;
  f=eval(fn);
  [n,f]
end
