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