Phase‐fitted,six‐step methods for solving x′′ = f(t,x)

An explicit, six‐step method of sixth order is presented and tuned for the numerical solution of x′′ = f(t,x). This method is explicit, hybrid, and uses two function evaluations (stages) per step. Its coefficients are varied and depend on the step size. This variance comes from the demand of the method to nullify the phase errors produced when solving the standard simple oscillator. The first and second derivative of this error vanish also. Numerical tests in a set of relevant problems illustrate the efficiency of the newly derived method.