La fuerza gravitacional ejercida sobre una masa m originada en un centro de masa M está dada por la siguiente fórmula:

transformando la ecuación anterior en una ecuación diferencial obtenemos:

la resolución de dicha ecuación diferencial nos proporciona los elementos matemáticos para poder elaborar un programa que nos permita calcular la posición de los planetas del sistema solar en el cielo terrestre para una fecha determinada.

Si bien la resolución antes mencionada puede resultar un poco compleja para aquellos que no conocen Cálculo Diferencial, he decidido incluir aquí un link para ella, de forma que aquellas personas interesadas puedan seguirla paso a paso.   Resolución

Basado en el desarrollo anterior he colocado aquí un programa en QBASIC de Microsoft (que acompaña el DOS de las versiones antiguas del Windows), llamado PLANETAS.BAS
La interface gráfica del programa es bastante simple, aparece un menu con el siguiente aspecto:

La opción [P]osicion retorna las coordenadas celestes (Ascención Recta, Declinación y Distancia) para la fecha solicitada. La pantalla tiene el siguiente aspecto:

La opción [G]rafico retorna un gráfico con las trayectorias de los planetas en la esferas celeste para un intervalo determinado a partir de la fecha solicitada. La pantalla tiene el siguiente aspecto:

La opción [O]rbitas retorna un gráfico con las trayectorias de los planetas en sus órbitas, vistos desde el sector norte del sistema solar (también aquí para un intervalo determinado a partir de la fecha solicitada). La pantalla tiene el siguiente aspecto:

Finalmente a opción [Q]uit nos permite salir del programa y volver al modo editor del QBASIC

	Para los interesados, el programa puede ser bajado aquí: (versión en portugués)	PLANETAS.BAS (18 KB)
	Recuerden que necesitarán abrir este programa dentro del editor de QBASIC y después apretar la tecla F5 para correrlo. Para aquéllos que no disponen de QBASIC en sus máquinas, aquí lo tienen disponible (QBASIC.EXE) junto con el archivo de ayuda (QBASIC.HLP) y el archivo del ícono (QBasic.ico):	Qbasic.zip (288 KB)
SkyGlobe DOS versión 3.6 Skyglobe for Windows version 2.0	Últimas versiones actualizadas! El Skyglobe for Windows realmente funciona!!!	skyglobewin.zip (1199 KB)