The family of Numerov-type methods that effectively uses seven stages per step is considered. All the coefficients of the methods belonging to this family can be expressed analytically with respect to four free parameters. These coefficients are trained through a differential evolution technique in order to perform best in a wide range of Keplerian-type orbits. Then it is observed with extended numerical tests that a certain method behaves extremely well in a variety of orbits (e.g., Kepler, perturbed Kepler, Arenstorf, Pleiades) for various steplengths used by the methods and for various intervals of integration.
