void bisection ( double &a, double &b, double tol, double f ( double x ),
  int &it );