pub fn expectation_x(psi: &WaveFunction1D) -> f64
Expectation value of position: = integral psi* x psi dx.