pub fn i1_scaled(x: f64) -> f64
This routine computes the scaled regular modified spherical Bessel function of first order, \exp(-|x|) i_1(x).