Use
ode78
and
ode89
to solve a celestial mechanics problem that requires high accuracy in each step from the ODE solver for a successful integration. Both
ode45
and
ode113
are unable to solve the problem using the default error tolerances. Even with more stringent error thresholds,
ode89
performs best on the problem due to the high accuracy of the Runge-Kutta formulas it uses in each step.