(*y = Sqrt[m ω/hbar ] x *)

psi0[y_] = E^( -y^2/2)
psi1[y_] = y E^(-y^2/2) /2
psi2[y_] = (4 y^2 - 2) E^(-y^2/2) /8

Plot[{psi0[y](*, Abs[psi0[y]]*)}, {y, -10, 10}, PlotRange ->Full]

Plot[{psi1[y](*, Abs[psi1[y]]*)}, {y, -10, 10}, PlotRange ->Full]

Plot[{psi2[y](*, Abs[psi2[y]]*)}, {y, -10, 10}, PlotRange ->Full]

SetDirectory["C:\\Users\\Peeter\\cygwin_home\\physicsplay\\notes\\blogit"]