pub fn normalize(psi: &mut [Complex], dx: f64)
Normalize a 1D wave function so that integral |psi|^2 dx = 1.