A Parametric Six-Step Method for Second-Order IVPs with Oscillating Solutions

In this paper, we develop an explicit symmetric six-step method for the numerical solution of second-order initial value problems (IVPs) with oscillating solutions. The proposed method is phase-fitted and incorporates a free coefficient as a parameter to optimize its performance. By exploring a wide range of values for this parameter, we computationally determine the periodicity interval. The objective of this procedure is to identify the range of the parameter’s values for which the method remains stable. Based on the output from the periodicity interval analysis, we then aim to define the optimal values for the parameter by numerically solving three initial value problems. The results guided us in identifying these optimal values and confirm the high efficiency of the new method. The method’s efficiency is further validated for the chosen optimal parameter value for specific oscillatory problems, where it is compared with well-known phase-fitted methods.