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