METU Departments of Mathematics

Runge-Kutta methods for Hamiltonian systems in non-standard symplectic two-form

B. Karasozen

Department of Mathematics, Middle East Technical University 06531 Ankara - Turkey

Runge-Kutta methods are applied to Hamiltonian systems on Poisson manifolds with a non-standard symplectic two-form. It has been shown that the Gauss Legendre Runge-Kutta (GLRK) methods and combination of the partitioned Runge-Kutta methods of Lobatto IIIA and IIIb type are symplectic up to the second order in terms of the step size. Numerical results on Lotka-Volterra and Kermack-McKendrick epidemic disease model reveals that the application of the symplectic Runge-Kutta methods preserves the integral invariants of the underlying system for long-time computations.\\

C.R. Categories: G 1.7.1

Keywords: Hamiltonian equations, Poisson manifolds, symplectic two-forms, symplectic Runge-Kutta methods, Lotka-Volterra equations, epidemic models