void bisection_min ( double f ( double x ), double a_init, double b_init, 
  double x_tol, double &a, double &m, double &b, int &it, int &nf );