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