func I1(x float64) float64 { return float64(C.gsl_sf_bessel_I1(C.double(x))) }
// First-order modified Bessel function of the first kind (I1) func ModifiedBesselFirstKindFirst(x float64) float64 { return float64(C.gsl_sf_bessel_I1(C.double(x))) }