Implicit symmetric and symplectic exponentially fitted modified Runge–Kutta–Nyström methods for solving oscillatory problems

Symplectic exponentially fitted RK and RKN methods have been extensively studied by many researchers. Due to their good property, they have been applied to many problems such as pendulum, Morse oscillator, harmonic oscillator, Lennard–Jones oscillator, Kepler’s orbit problem, and so on. In this paper, we construct an implicit symmetric and symplectic exponentially fitted modified Runge–Kutta–Nyström (ISSEFMRKN) method. The new integrator integrates exactly differential systems whose solutions can be expressed as linear combinations of functions from the set [Formula: see text] , [Formula: see text] , or equivalently [Formula: see text] when [Formula: see text] , [Formula: see text] . When [Formula: see text] approaches zero, the ISSEFMRKN method reduces to the classical symplectic, symmetric RKN integrator. Numerical experiments are accompanied to show the efficiency and competence of the new method compared with some efficient codes in the literature.