Families of third and fourth algebraic order trigonometrically fitted symplectic methods for the numerical integration of Hamiltonian systems

The numerical integration of Hamiltonian systems by symplectic and trigonometrically fitted (TF) symplectic method is considered in this work. We construct new trigonometrically fitted symplectic methods of third and fourth order. We apply our new methods as well as other existing methods to the numerical integration of the harmonic oscillator, the 2D harmonic oscillator with an integer frequency ratio and an orbit problem studied by Stiefel and Bettis.