pub fn Ci_e(x: f64) -> (Value, Result)
This routine computes the Cosine integral Ci(x) = -\int_x^\infty dt \cos(t)/t for x > 0.