(*Morse oscillator with mass  mu, the potential (depth DD) for  J=0 and for  J=1*)
DD=0.010
mu=400. 
alpha=1.
a=((Sqrt[2 mu DD])/alpha)
beta=(alpha/(Sqrt[8 mu DD]))
kwant=2 alpha (Sqrt[(DD/(2 mu))])
kforce=2 (alpha^2) DD
BB=0.002
vMor=DD (E^(-alpha (x-2))) ((E^(-alpha (x-2)))-2)
Show[
Plot[vMor,{x,0.,7.},PlotPoints->101,
PlotStyle->{GrayLevel[0.],Thickness[0.005]},AxesLabel->{R,E},LabelStyle->Directive[Bold,Larger]],
Plot[{vMor+BB*J*(J+1)/x^2,J={1,2}},{x,0.,7.},PlotStyle->Directive[Green,Thickness[0.005]]]
]
Plot[vMor,{x,0.,7.},PlotPoints->101,
PlotStyle->{GrayLevel[0.],Thickness[0.005]},AxesLabel->{R,E},LabelStyle->Directive[Bold,Larger]]
Plot[vMor+0.008*Exp[-(x-4)^2],{x,0.,7.},PlotPoints->101,
PlotStyle->{GrayLevel[0.],Thickness[0.005]},AxesLabel->{R,E},LabelStyle->Directive[Bold,Larger]]

Null