A Fortran 90 separable Hamiltonian system solver Original Research Article

We discuss a prototype Fortran 90 separable Hamiltonian system solver and present a template illustrating its use. The solver permits a wide choice of symplectic and non-symplectic integrators, and fixed and error-controlled integration step sizes. It is designed to permit insertion of new formulas as they become available.