An adaptive optimized Nystrom method for second-order IVPs

This research work deals with the development, analysis, and implementation of an adaptive optimized one-step Nystrom method for solving second-order initial value problems of ODEs and time-dependent partial differential equations. The new method is developed through a collocation technique with a new approach for selecting the collocation points. An embedding-like procedure is used to estimate the error of the proposed optimized method. The current approach has been used to compute efficiently approximate solutions to general second-order IVPs. The numerical experiments demonstrate that the introduced error estimation and step-size control strategy presented in this manuscript have produced a good performance compared to some of the other existing numerical methods.