(*the ground state wave function for the Morse oscillator*)
alpha=0.2
m=1
Den=3
a=Sqrt[2*m*Den]/alpha
psi1=Exp[-a*Exp[-alpha*x]]*(2a*Exp[-alpha*x]^(a-0.5))
Plot[psi1,{x,-5,5}]
Hpsi1=-hkresl^2/(2*m)*psi1''+Den*Exp[-alpha*x]*(Exp[-alpha*x]-2)*psi1