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