pub fn dc_prefactor(theta_z: f64) -> f64
Precompute the x-independent DC prefactor: (4α/3π) θ_z² × I₄^pl / (1 + 14.16 θ_z)
The only x-dependent part remaining is H_dc(x) = exp(-2x) × polynomial.