Les pasamos el link de una demostración hecha con MATHEMATICA que demuestra las trayectorias de cada masa del problema de Kepler:

http://demonstrations.wolfram.com/TheCelestialTwoBodyProblem/.

La evolución temporal la obtienen resolviendo la ecuación de Kepler  en forma aproximada.

Un control permite elegir el cociente de masas de 1 a 50. El control de excentricidad permite órbitas desde circular (e=0) hasta elíptica alargada (e=0.9). La llamada longitud de apoapsis permite  achicar o agrandar la órbita, y finalmente el angulo de apoapsis define la orientacion de la elipse en el plano de la órbita.

Si tienen instalado un MATHEMATICA deberán bajar la notebook. Si no, pueden instalar el player de MATHEMATICA desde:

http://demonstrations.wolfram.com/download-cdf-player.html

Tomen nota de los cambios en la velocidad a medida que los cuerpos se acercan, aumenta la velocidad pues el potencial se hace mas negativo.

Disfrútenlo.