Meine Antwort war auch eher eine auf Torpedos Antwort.
Jelly hat es besser formuliert als ich : wie gesagt, die Formel an sich ist einfach, aber je mehr Sonderfälle man einbaut, um so schlimmer wird sie. Ich würde das Problem deshalb immer nur auf 2 Massen reduzieren und Elipsen, Bahnneigung usw. weg lassen. Z.B. so : der Mond dreht sich nur um die Erde und wird von Sonne usw. nicht beeinflußt, genauso wenig wie ein Satellit von Sonne oder Mond.
Ein Swing By Manöver zu simulieren, au Backe.
Alleine schon den Winkel der Ablenkung zu berechnen um den die Sonde am Mars umgeleitet und beschleunigt wird. Dann spielt die Zeit noch eine gewaltige Rolle, denn wo befindet sich die Erde überhaupt, falls die Sonde per swing by vom Mars zurückkommt ? Verfehlt sie die Erde, dann gute Nacht, bzw. guten Tag Sonne !
Solange kein Modell für so was steht, brauchst du auch keinen Source.