(*Oscylator Morse' a .  Czastka ma mase mu,  
      a potecjal odpowiada energii
wiazania rownej DD .  Najpierw porownujemy krzywa  energii potencjalnej
dla oscylatora harmonicznego z taka krzywa dla oscylatora Morse (o tej samej
 
              stalej silowej kforce) .  
          Potem wykreslamy poziomy energetyczne
dla oscylatora Morse' a dla czastki o masie mu,  
      a nastepnie dla czastki
o masie dwa razy wiekszej .  Pozniej wykreslamy poziomy,  
      gdy 
studnia nieco sie wyplyca *) 
DD = 0.01
mu = 1840.  
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
z = (E^(-alpha x))
VMorse = DD z (z - 2)
VHooke = -DD + ((kforce (x^2))/2)
Plot[{VMorse, VHooke}, {x, -2. , 5. }]
Del = a - 0.5
E0 = -DD + kwant 0.5 - kwant 0.25 beta
E1 = -DD + kwant 1.5 - kwant (1.5^2) beta
E2 = -DD + kwant 2.5 - kwant (2.5^2) beta
E3 = -DD + kwant 3.5 - kwant (3.5^2) beta
E4 = -DD + kwant 4.5 - kwant (4.5^2) beta
E5 = -DD + kwant 5.5 - kwant (5.5^2) beta
Plot[{VMorse, E0, E1, E2, E3, E4, E5}, {x, -2. , 5. }]
DD = 0.01
mu = 1840.  2
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
z = (E^(-alpha x))
VMorse = DD z (z - 2)
VHooke = -DD + ((kforce (x^2))/2)
Del = a - 0.5
E0 = -DD + kwant 0.5 - kwant 0.25 beta
E1 = -DD + kwant 1.5 - kwant (1.5^2) beta
E2 = -DD + kwant 2.5 - kwant (2.5^2) beta
E3 = -DD + kwant 3.5 - kwant (3.5^2) beta
E4 = -DD + kwant 4.5 - kwant (4.5^2) beta
E5 = -DD + kwant 5.5 - kwant (5.5^2) beta
E6 = -DD + kwant 6.5 - kwant (6.5^2) beta
E7 = -DD + kwant 7.5 - kwant (7.5^2) beta
E8 = -DD + kwant 8.5 - kwant (8.5^2) beta
Plot[{VMorse, E0, E1, E2, E3, E4, E5, E6, E7, E8}, {x, -2. , 
        5. }]
DD = 0.008
mu = 1840. 
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
z = (E^(-alpha x))
VMorse = DD z (z - 2)
VHooke = -DD + ((kforce (x^2))/2)
Del = a - 0.5
E0 = -DD + kwant 0.5 - kwant 0.25 beta
E1 = -DD + kwant 1.5 - kwant (1.5^2) beta
E2 = -DD + kwant 2.5 - kwant (2.5^2) beta
E3 = -DD + kwant 3.5 - kwant (3.5^2) beta
E4 = -DD + kwant 4.5 - kwant (4.5^2) beta
Plot[{VMorse, E0, E1, E2, E3, E4}, {x, -2. , 5. }]

Null