“…For several decades, there has been strong interest in searching for better numerical methods to integrate first-order and second-order initial value problems, because these problems are usually encountered in celestial mechanics, quantum mechanical scattering theory, theoretical physics and chemistry, and electronics. Generally, the solution of (1) is periodic, so it is expected that the result produced by some numerical methods preserves the analogical periodicity of the analytic solution [9][10][11][12][13][14][15][16][17][18][19][20][21][22]. Computational methods involving a parameter proposed by Gautschi [8], Jain et al [10], Sommeijer and et al [30] and Steifel and Bettis [31] yield numerical solution of problems of class (1).…”