Poniżej prezentuję niedawno (1993) odkryte choreografię w problemie 3 ciał. Obecnie tego typu rozwiązań znamy znacznie więcej.

G=1;                       (* Gravitational constant *)
m=1;                       (* Mass of each body *)

(*One full circuit of the 8 takes≈ 6.32591398 time‑units*)
T=6.32591398;

threeBodyFigure8=NBodySimulation["InverseSquare",
(*----------3‑body figure‑8 initial data----------*)
{
<|"Mass"->1.03m,"Position"->{-0.97000436,0.24308753,0.0},"Velocity"->{0.466203685,0.43236573,0.}|>,
<|"Mass"->0.99m,"Position"->{0.97000436,-0.24308753,0.},"Velocity"->{0.466203685,0.43236573,0.}|>,
<|"Mass"->m,"Position"->{0.,0.,0.},"Velocity"->{-0.93240737,-0.86473146,0.}|>
}
,4 T    (*simulate four periods*)
];

ParametricPlot3D[Evaluate[threeBodyFigure8[All,"Position",t]],{t,0,4T}]

Manipulate[Show[
Graphics3D[Point[Evaluate[threeBodyFigure8[All,"Position",t]]]],
PlotRange->1.3{{-1,1},{-1,1},{-1,1}}
]
,{t,0,4T}]

srodekMasy=threeBodyFigure8[All,"Mass"].threeBodyFigure8[All,"Position",t]/Total[threeBodyFigure8[All,"Mass"]];

ParametricPlot3D[srodekMasy,{t,0,4T}]

ParametricPlot3D[Evaluate[threeBodyFigure8[All,"Position",t]-threeBodyFigure8[3,"Position",t]],{t,0,4T}]