Eighth order, Numerov-like schemes with coefficients tailored for superior performance on ODE systems with oscillatory solutions

Second order Ordinary Differential Equations (ODE) were considered. Numerov-like techniques employing effectively seven stages per step and sharing eighth algebraic order were under examination for numerically solving them. The coefficients of these methods were contingent on four independent parameters. To tackle issues with oscillatory solutions, we typically aimed to fulfill specific criteria such as minimizing phase-lag, expanding the periodicity interval, or even neutralizing amplification errors. These latter attributes stemmed from a test problem mimicking an ideal trigonometric trajectory. Here, we suggested training the coefficients of the chosen method family across a broad spectrum of pertinent problems. Following this training using the differential evolution method, we identified a particular method that surpassed others in this category across an even broader array of oscillatory problems.