Ce problème est aussi examiné dans la page Force centrale dans laquelle on détermine la trajectoire uniquement par intégration numérique des équations du système formé par les deux masses.
Ici, on commence par calculer les élément de la trajectoire dans le cas elliptique à partir des conditions initiales :
x = xo, y = 0, vx = 0, vy = vyo.
Ensuite, on vérifie par intégration numérique que cette trajectoire est suivie par la masse mobile.
Comme la force d'attraction entre les deux masses est f = − k.m / r², l'énergie potentielle est U = − k.m / r.
La conservation de l'énergie de ce système isolé donne :

v² / 2 − k / r = E (énergie mécanique par unité de masse) Avec votre cours de mécanique, retrouvez les relations suivantes qui permettent de calculer a et b (demi-axes) de l'ellipse et T (période de révolution). v et r sont relatifs au même point de la trajectoire.

La relation (1) met en évidence une vitesse critique v²_c = 2.k / r. Pour les vitesses initiales supérieures, le mouvement est hyperbolique. Cette relation montre aussi que pour v² = k / r, la trajectoire est circulaire.

---

Utilisation :
Avec la souris, glissez le point bleu situé à l'extrémité du vecteur vitesse initiale (en rouge) afin de modifier les conditions initiales.
Le programme affiche les paramètres de la trajectoire.
Les unités sont arbitraires et la constante k est égale à 10. (4.π.²/ k est donc voisin de 4).

La case à cocher "Aires" permet de visualiser ou non la loi des aires.
En mode normal, le programme affiche les "rayons vecteurs" joignant le mobile aux foyers de l'orbite et en rouge le vecteur vitesse qui est tangent à la trajectoire.
Lors de la modification des conditions initiales, le programme affiche le graphe 2.ln(T) − ln(4) en fonction de 3.ln(a) afin de visualiser la 3° loi de Kepler : lors des modifications, le point représentatif se déplace sur la première diagonale du graphe.

Les valeurs trop faibles de r et de v sont automatiquement corrigées.