Schroeding equation : harmonic oscillator

From https://demonstrations.wolfram.com/QuantizedSolutionsOfThe1DSchroedingerEquationForAHarmonicOsc/

EnergyLevel[e_]:=Line[{{-Sqrt[e],e},{Sqrt[e],e}}];
shoLevelPlot[e_]:=Plot[x^2,{x,-2.75,2.75},Epilog->Flatten[{Red,Thick,Line[{{-3.35,e},{3.35,e}}],Dashed,Black,Thickness[0.004],Table[EnergyLevel[w],{w,0.5,7.5,1}]},1],PlotStyle->{Black,Thickness[0.005]},AxesLabel->{"x"},PlotLabel->"harmonic oscillator potential",ImageSize->{275,275}]

shoLevelPlot[3.5]

Solve Schoedinger eq .

-h/2m Psi[x]’’+1/2 m om^2 x^2 Psi[x] = energy Psi[x]

with h=1, m=1, om=1

?Manipulate

Manipulate[Plot[a*x^2+b,{x,-3,3},PlotRange->{{-3,3},{-20,20}}],{a,-1,1},{b,{-2,0,2}}]

sol[energy_,pf_]:=NDSolve[{-0.5y''[x]+0.5 (x^2)y[x]-energy y[x]==0,y[1.]==1.,y[-1.]==pf 1.},y,{x,-5.9,5.9}]

Manipulate[
Grid[{{Plot[Evaluate[y[x]/.sol[energy,pf]],{x,-5.5,5.5},PlotRange->{{-6,6},{-10,10}},Frame->True,PlotStyle->Thickness[.007],AxesLabel->{"x",""},PlotLabel->"wavefunction ψ(x)",ImageSize->{275,275},ImagePadding->{{25,25},{25,25}}],shoLevelPlot[energy]}}],{{energy,0.5},0.,7.5,.00001},
{{energy,0.5,"allowed energy values (even parity)"},{0.5,2.5,4.5,6.5},Enabled->pf==1},
{{energy,0.5,"allowed energy values (odd parity)"},{1.5,3.5,5.5,7.5},Enabled->pf==-1},
{{pf,1,"parity"},{1->"even",-1->"odd"}}]