m1 = 0.035453/(6.022*10^23)

μ=m1 m1/(m1+m1)

(OverBar[ν])=56000.0

c=2.99*10^8

ℏ=(6.626*10^-34) / (2 Pi)

k = (c (OverBar[ν]) 2Pi)^2μ

α = ((k μ)/ ℏ^2)^0.5

NormConst[v_]:= 1/(2^v Factorial[v])^(1/2)(α/Pi)^(1/4)

Psi[x_,v_]:=NormConst[v] HermiteH[v,α^(1/2)x]Exp[-α x^2/2]

xvar1 = -1*10^-10

xvar2 = -xvar1

Plot[Psi[x,1],{x,xvar1,xvar2}, PlotRange -> Full]