pub fn ei(x: f64) -> f64
Exponential integral Ei(x) for x > 0.
Uses a series expansion for small x and an asymptotic series for large x.