Potential energy for diatomic molecules

b = 1 ;

lo = 2 ;

potener[x_] := (1 - Exp[-b (x - lo)])^2 - 1

fMorse = Plot[{potener[z]}, {z, 0, 5}, Frame→True, PlotStyle→ {Thickness[.02]}] ;

[Graphics:../HTMLFiles/MATH5_TUT1_611.gif]

Harmonic oscillator  potential energy

k = D[potener[x], {x, 2}] <br />

2 ^(4 - 2 x) - 2 ^(2 - x) (1 - ^(2 - x))

kk = k/.x→lo

2

[Graphics:../HTMLFiles/MATH5_TUT1_617.gif]

Display["MHO.JPG",fHarmonic,"JPEG"]

fHarmonic2 = Show[fHarmonic,Graphics[{Text["Morse",{3.7,-.5}],Text["Harmonic",{3.6,-0.1}]}]];

[Graphics:../HTMLFiles/MATH5_TUT1_618.gif]

Display["MHO.JPG",fHarmonic2,"JPEG"]

Created by Mathematica  (September 7, 2006) Valid XHTML 1.1!