In[]:=
a=1.5;​​Ω=0.8;​​ξe=0.1;​​ξg=0.5;​​ξa=0.0;​​rα=Sqrt[0.25];​​Ma=
1
ξg
ξg
2
rα
-
2
ξe
+2ξe*ξg
;​​qppa=
h''[t]
α''[t]
;​​Ka=
1
ξe
ξe
2
ξe
+
2
rα
/
2
Ω
;​​Fexta=
0
1/a
0
ξa/a
;​​i=1;​​Va=0;​​While[Va<1.8,​​Stato=Inverse[Ma].((-
2
Ω
*Ka)-(
2
Va
*Fexta));​​Eigenvalues[Stato];​​a2=%[[1]];​​b2=%%[[2]];​​Print[a2];​​λ[i]={a2,b2};​​i=i+1;​​Va=Va+0.01;]​​