pub fn _1_e(x: f64) -> (Value, Result)
This routine computes the first-order Debye function D_1(x) = (1/x) \int_0^x dt (t/(e^t - 1)).