Composite integrators for bi-Hamiltonian systems

Symmetric composition methods are applied to dynamical systems in bi-Hamiltonian form; to Lotka-Volterra equations, to the completely integrable Lorenz equation and to the periodic Toda lattice. The numerical results obtained show that the Hamiltonians are preserved with high accuracy and the periodicity of the solutions are retained.