Function dc_heating_integral

Source
pub fn dc_heating_integral(
    x_grid: &[f64],
    delta_n: &[f64],
    theta_z: f64,
    theta_e: f64,
) -> f64
Expand 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.