eval ("function y = f (x), y = x^3-3*x-1; end;")
  
  y=bisection('f',1,2,50,.000001)