pub fn qho_wavefunction(n: u32, x: f64, omega: f64, mass: f64, hbar: f64) -> f64
QHO wave function psi_n(x) = N_n * H_n(xi) * exp(-xi^2/2) where xi = sqrt(m*omega/hbar) * x