Aunque las leyes de Kepler pueden ser demostradas empíricamente aplicando las leyes de Newton a la dinámica de dos partículas atraídas por interacción gravitatoria, una fórmula explícita para el movimiento como función del tiempo permanece indefinida. Este documento propone una solución cuasi-analítica para abordar este desafío. Aproxima la dinámica real de los cuerpos celestes con un grado de precisión satisfactorio y un costo computacional mínimo. Este problema está estrechamente relacionado con la ecuación de Kepler, ya que resolver las ecuaciones de movimiento como función del tiempo también proporciona una solución a la ecuación de Kepler. Los resultados se presentan para cada planeta del sistema solar, incluido Plutón, y la solución se compara con las órbitas reales.