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;