Me parecía raro que el error sea prácticamente el mismo al inicio y al final de cada órbita. Viéndolo con más detall, se puede ver que el error final es un poco más grande que el error inicial y que la solución numérica no es un circunferencia perfecta, si no que es un poco más elíptica. Me imagino que con pasos más grandes esto se ve claramente. En ese caso, parece que los gráficos están bien. Saludos!
