fn='exp(x)-3*x^2';
xl=0
xu=2
x=.5*(xu+xl)
while abs(eval(fn)) > 1.e-6
  x=xl;fl=eval(fn);
  x=xu;fu=eval(fn);
  xr=xl-fl*(xu-xl)/(fu-fl);
  [xl,xu,xr]
  x=xr;fr=eval(fn);
  if fl*fr <= 0 
   xu=xr;
  elseif fr*fu <= 0
   xl=xr;
  else
   disp(' Error - root not bisected ')
   break
  end
   x=xr;
   eval(fn)   
end
