Deutsch: Modifikation der Anziehungskraft durch Einführung einer Abhängigkeit des Typs ${\displaystyle r/(r^{2}+\epsilon ^{2})^{3/2}}$ nach Aarseth S.J., Dynamical evolution of clusters in galaxies I, in: Monthly Notices of the Royal Astronomical Society Band 126, S.223 ff, 1963, hier demonstriert anhand der Simulation eines Zweikörpersystems mit dem Hermite-Polynome-Verfahren und dynamischen Zeitschritten über einen Zeitraum von 10 Jahren. ${\displaystyle \epsilon }$ ist auf 0 (blaue Kurve), 0.1 AE (rote Kurve) und 0.5 AE (grüne Kurve gesetzt. Der Körper im Zentrum weist eine Sonnenmasse, das umlaufende Objekt eine Erdmasse auf. Die Startpositionen sind (0/0) und (1 AE / 0), die Anfangsgeschwindigkeiten (0/0) und (0 / 35 km/s). Die Modifikation beruht primär auf einer solchen der potentiellen Energie gemäß einer Abhängigkeit der Form ${\displaystyle 1/{\sqrt {r^{2}+\epsilon ^{2}}}}$