pst-ode et 3 corps animés d’un mouvement périodique sur une trajectoire en forme de huit

Alexander Grahn a récemment mis en ligne un package dédié à la résolution numérique des équations différentielles ordinaires : pst-ode.
http://tug.ctan.org/tex-archive/graphics/pstricks/contrib/pst-ode
Il utilise une méthode un peu plus performante que celles utilisées par Dominique Rodriguez dans pst-eqdf et pstricks-add.
Ceci est un premier test qui ne fait pas le tour des possibilités du package et auquel Alexander pourra certainement apporter des améliorations.
 On trouve beaucoup de renseignements sur ce mouvement en forme de 8 ainsi que des études théoriques sur les chorégraphies de N corps soumis uniquement aux forces gravitationnelles sur la page de Carles Simó :
http://melusine.eu.org/syracuse/swf/1-nbody
Par ailleurs Jean-Michel Sarlat a réalisé de magnifiques animations au format flash, très fluides, de ces chorégraphies, d’après les données calculées par Carles Simó :
http://melusine.eu.org/syracuse/swf/1-nbody/
Ces chorégraphies peuvent, par ailleurs, être très simplement animées avec pst-animate, peut-être prochainement ?

Voici une image de l'animation :

Cette animation incluse dans le fichier pdf et son fichier source (.tex), peuvent être téléchargés sans le sous-dossier Ncorps de  :

 gravitation