pub fn Shi_e(x: f64) -> Result<Result, Value>
This routine computes the integral Shi(x) = \int_0^x dt \sinh(t)/t.