pub fn g_bb(x: f64) -> f64
Blackbody derivative: G_bb(x) = x e^x / (e^x - 1)² This is -x ∂n_pl/∂x = x²/(4T) ∂B_ν/∂T normalized