ClearAll[h,v,m,ω,x,ψ,e,n]

(*Setting Universal Constants*)
h=1;

Manipulate[
Plot[

{(1/2)*m*(ω^2)*(x^2),(((((m*ω)/(Pi*h))))^(1/4))*(1/(Sqrt[(2^n)*Factorial[n]]))*HermiteH[n,Sqrt[((m*ω)/h)]*x]*Exp[-((m*ω)/(2*h))*(x^2)]+ (n+(1/2))*h*ω,
(n+(1/2))*h*ω },

{x,-3,3}, PlotRange -> {0,7}, PlotLabel -> "Wave Function"], 

{{n,0,"Energy Level"},0,5,1},
{{m,5,"Mass"},3,7}, 
{{ω,1,"Angular Velocity"},0.5,1}]

Manipulate[
Plot[

{(1/2)*m*(ω^2)*(x^2),((((((((m*ω)/(Pi*h))))^(1/4))*(1/(Sqrt[(2^n)*Factorial[n]]))*HermiteH[n,Sqrt[((m*ω)/h)]*x]*Exp[-((m*ω)/(2*h))*(x^2)]))^2)+ (n+(1/2))*h*ω,
(n+(1/2))*h*ω },

{x,-3,3}, PlotRange -> {0,7}, PlotLabel -> "Probability Wave Function"], 

{{n,0,"Energy Level"},0,5,1},
{{m,5,"Mass"},3,7}, 
{{ω,1,"Angular Velocity"},0.5,1}]