pub fn jl(l: i32, x: f64) -> f64
This routine computes the regular spherical Bessel function of order l, j_l(x), for l >= 0 and x >= 0.