pub fn delta_rho_over_rho(x_grid: &[f64], delta_n: &[f64]) -> f64
Compute fractional energy in distortion: Δρ/ρ = ∫x³ Δn dx / G₃