void golden_section ( double f ( double x ), int n, double x_tol, double &a, 
  double &b, int &it );