2733db63-e93b-4ac1-9e29-d01618afbd8b 117 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

Arenstorf orbit

The Arenstorf orbits are closed trajectories of the restricted three-body problem. Two bodies of masses m and (1-m) moving in a circular rotation, and a third body of negligible mass moving in the same plane (such as a satellite-earth-moon system.)

t y μ D μp 1 μ - : y1 y 3 el : y2 y 4 el : r1 y 1 el μ + (2^y 2 el 2^+ : r1 r1 r1 sqrt * : r2 y 1 el μp - (2^y 2 el 2^+ : r2 r2 r2 sqrt * : y3 y 1 el 2 y 4 el * + μp r1 / y 1 el μ + (* - μ r2 / y 1 el μp - (* - : y4 y 2 el 1 μp r1 / - μ r2 / - (* 2 y 3 el * - : y1 y2 y3 y4 stack 101 line :

4 loops: Mass parameter: m = 0.012277471 Coordinates: y1(0) = 0.994 y2(0) = 0 y'1(0) = 0 y'2(0) = -2.00158510637908252240537862224 Period: tper = 17.0652165601579625588917206249

Y 0.994 00 2.00158510637908252240537862224 - 41 mat :

μ 0.012277471 :

t0 0 :

n 200 :

tmax 17.0652165601579625588917206249 :

res Y t0 tmax n t y μ D dn_GearsBDF :

res Y t0 tmax n t y μ D dn_ImplicitRK5 :

res Y t0 tmax n t y μ D dn_AdamsMoulton :

res Y t0 tmax n t y μ D dn_ExplicitRK45 :

T res 1 col :

Y11 res 2 col :

Y12 res 3 col :

Y11 Y12 augment

3 loops: Mass parameter: m = 0.012277471 Coordinates: y1(0) = 0.994 y2(0) = 0 y'1(0) = 0 y'2(0) = -2.0317326295573368357302057924 Period: tper = 11.124340337266085134999734047

Y 0.994 00 2.0317326295573368357302057924 - 41 mat :

μ 0.012277471 :

t0 0 :

n 200 :

tmax 11.124340337266085134999734047 :

res Y t0 tmax n t y μ D dn_GearsBDF :

res Y t0 tmax n t y μ D dn_ImplicitRK5 :

res Y t0 tmax n t y μ D dn_AdamsMoulton :

res Y t0 tmax n t y μ D dn_ExplicitRK45 :

T res 1 col :

Y11 res 2 col :

Y12 res 3 col :

Y11 Y12 augment

2 loops: Mass parameter: m = 0.012277471 Coordinates: y1(0) = 1.2 y2(0) = 0 y'1(0) = 0 y'2(0) = -1.049357510 Period: tper = 6.192169331

Y 1.2 00 1.049357510 - 41 mat :

t0 0 :

n 200 :

tmax 6.192169331 :

μ 0.012277471 :

res Y t0 tmax n t y μ D dn_GearsBDF :

res Y t0 tmax n t y μ D dn_ImplicitRK5 :

res Y t0 tmax n t y μ D dn_AdamsMoulton :

res Y t0 tmax n t y μ D dn_ExplicitRK45 :

T res 1 col :

Y11 res 2 col :

Y12 res 3 col :

Y11 Y12 augment