pub fn E1_e(x: f64) -> (Value, Result)
This routine computes the exponential integral E_1(x),
E_1(x) := \Re \int_1^\infty dt \exp(-xt)/t.