Canonical difference scheme for the Hamiltonian equation

In this paper we construct canonical difference schemes of any order accuracy based on Padé approximation for Linear canonical systems with constant coefficients. For non-linear Hamiltonian equations we will use an infinitesimally canonical transformation to construct canonical schemes of any order accuracy.