7 #include <gsl/gsl_math.h> 
    8 #include <gsl/gsl_sf_hyperg.h> 
    9 #include <gsl/gsl_sf_bessel.h> 
   15   double confluentHypergeometric1F1(
const double a, 
const double b, 
const double x)
 
   17     return gsl_sf_hyperg_1F1(a, b, x);
 
   20   double besselI1(
const double x)
 
   22     return gsl_sf_bessel_I1(x);
 
   25   bool rfe(
const double a, 
const double b, 
double tol)
 
   27     return gsl_fcmp(a, b, tol) == 0;