B. Karasozen

A reversible Hamiltonian in Birkhoff-Gustavson normal form is solved using symplectic and symmetric Runge-Kutta methods with constant and variable step sizes. The numerical methods are compared with respect to Hamiltonian preservation and the accuracy of solutions.