pub fn I1_e(x: f64) -> (Value, Result)
This routine computes the regular modified cylindrical Bessel function of first order, I_1(x).