Prueba de la precisión de una simulación del sistema solar

He escrito un programa Java para modelar el sistema solar y necesito probar su precisión. ¿Cuáles serían algunas buenas maneras de hacer eso? He comprobado la conservación de la energía y el momento angular. Estoy usando el algoritmo de velocidad verlet, que parece haber conservado perfectamente ambos.

Me preguntaba si había alguna manera de verificar si estaban de vuelta donde comenzaron después de un período orbital, pero no pude encontrar la forma de programar una verificación. Realmente me gustaría una forma de medir el efecto de cambiar el paso de tiempo, pero no sé cuál sería una buena medida. ¿Hay algo más que pueda hacer?

¿A qué detalle de las interacciones va su programa? ¿Tiene en cuenta la interacción de los planetas? ¿Utiliza GR o mecánica newtoniana? Podría intentar comprobar si se cumplen las leyes de Kepler.
@MartinUeding Solo mecánica newtoniana, los planetas interactúan pero solo se modelan como masas puntuales.

Respuestas (1)

Una forma simple y genérica de probar tales programas es evolucionar hacia atrás en el tiempo: elija algunas condiciones iniciales, evolucione por un tiempo T , luego invertir todas las velocidades, evolucionar nuevamente por un tiempo T , y verifique qué tan cerca está de las condiciones iniciales.