- #1
Dustinsfl
- 2,281
- 5
The L4 position is stable in the Earth Moon and I perturbing a satellite by km in the x direction to see the trajectory over the course of the year. However, the satellite isn't moving. Can anyone see if there is something wrong? I gave the satellite no initial velocity.
Code:
In[2587]:= ClearAll["Global`*"]
me = 5.974*10^(24);
mm = 7.348*10^(22);
G = 6.67259*10^(-20);
re = 6378;
rm = 1737;
r12 = 384400;
In[2594]:= \[Mu] = G*(me + mm);
\[Pi]1 = me/(me + mm);
\[Pi]2 = mm/(me + mm);
M = me + mm;
\[CapitalOmega] = Sqrt[\[Mu]/r12^3];
\[Mu]1 = G*me;
\[Mu]2 = G*mm;
In[2601]:= xl4 = 384400/2 - 4671
yl4 = Sqrt[3]/2*384400 // N
Out[2601]= 187529
Out[2602]= 332900.
In[2612]:= r0 = {xl4+1, yl4, 0}
v0 = {0, 0, 0};
Out[2612]= {187529, 332900., 0}
{187529, 332900.16521473817`, 0}
In[2614]:=
s = NDSolve[{x1''[t] -
2*\[CapitalOmega]*x2'[t] - \[CapitalOmega]^2*
x1[t] == -\[Mu]1/(Sqrt[(x1[t] + \[Pi]2*r12)^2 +
x2[t]^2])^3*(x1[t] + \[Pi]2*
r12) - \[Mu]2/(Sqrt[(x1[t] - \[Pi]1*r12)^2 +
x2[t]^2])^3*(x1[t] - \[Pi]1*r12),
x2''[t] -
2*\[CapitalOmega]*x1'[t] - \[CapitalOmega]^2*
x2[t] == -\[Mu]1/(Sqrt[(x1[t] + \[Pi]2*r12)^2 + x2[t]^2])^3*
x2[t] - \[Mu]2/(Sqrt[(x1[t] - \[Pi]1*r12)^2 + x2[t]^2])^3*
x2[t],
x3''[t] == -\[Mu]1/(Sqrt[(x1[t] + \[Pi]2*r12)^2 + x2[t]^2])^3*
x3[t] - \[Mu]2/(Sqrt[(x1[t] - \[Pi]1*r12)^2 + x2[t]^2])^3*
x3[t], x1[0] == r0[[1]], x1'[0] == v0[[1]], x2[0] == r0[[2]],
x2'[0] == v0[[2]], x3[0] == r0[[3]], x3'[0] == v0[[3]]}, {x1, x2,
x3}, {t, 0, 24*3600*365}];
In[2617]:= ParametricPlot3D[
Evaluate[{x1[t], x2[t], x3[t]} /. s], {t, 0, 200000},
PlotStyle -> {Red, Thick}]
Last edited: