pub fn dc_heating_integral(
x_grid: &[f64],
delta_n: &[f64],
theta_z: f64,
theta_e: f64,
) -> f64Expand description
DC emission rate integrated over frequency (for electron temperature equation).
H_DC = (1/(4 G₃ θ_z)) ∫ [1 − n(e^{x_e} − 1)] K_DC(x)/x³ · x³ dx
This enters the electron temperature evolution as a cooling term.