  Some explicit algorithms for higher order symplectic integration of a large
class of Hamilton's equations have recently been discussed by Mushtaq \emph{et.
al}. Here we present a Python program for automatic numerical implementation of
these algorithms for a given Hamiltonian, both for double precision and
multiprecision computations. We provide examples of how to use this program,
and illustrate behaviour of both the code generator and the generated solver
module(s).
