pub fn il_scaled(l: i32, x: f64) -> f64
This routine computes the scaled regular modified spherical Bessel function of order l, \exp(-|x|) i_l(x)