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