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