pub fn j1_e(x: f64) -> (Value, Result)
This routine computes the regular spherical Bessel function of first order, j_1(x) = (\sin(x)/x - \cos(x))/x.