Abstract
A family of explicit, eighth‐order, four‐step methods for the numerical solution of
y′′=f(x,y) is considered. This family is derived through an interpolatory approach after using three stages (ie, function evaluations) per step. In the present work, we alter three of the coefficients of such a method in order to become phase fitted. We conclude with numerical tests over a set of problems justifying our effort of dealing with the new methods.