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