pub fn delta_n_over_n(x_grid: &[f64], delta_n: &[f64]) -> f64
Compute fractional photon number change: ΔN/N = ∫x² Δn dx / G₂