pub fn E2(x: f64) -> f64
This routine computes the second-order exponential integral E_2(x),
E_2(x) := \Re \int_1^\infty dt \exp(-xt)/t^2.