Runge–Kutta Pairs of Orders 5(4) Trained to Best Address Keplerian Type Orbits

The derivation of Runge–Kutta pairs of orders five and four that effectively uses six stages per step is considered. The coefficients provided by such a method are 27 and have to satisfy a system of 25 nonlinear equations. Traditionally, various solutions have been tried. Each of these solutions makes use of some simplified assumptions and offers different families of methods. Here, we make use of the most celebrated family to appear in the literature, where we may use as the last stage the first function evaluation from the next step (FSAL property). The family under consideration has the advantage of being solved explicitly. Actually, we arrive at a subsystem where all the coefficients are found with respect to five free parameters. These free parameters are adjusted (trained) in order to deliver a pair that outperforms other similar pairs of orders 5(4) in Keplerian type orbits, e.g., Kepler, perturbed Kepler, Arenstorf orbit or Pleiades. The training uses differential evolution technique. The finally proposed pair has a remarkable performance and offers on average more than a digit of accuracy in a variety of orbits.