  We propose Poisson integrators for the numerical integration of separable
Poisson systems. We analyze three situations in which the Poisson systems are
separated in three ways and the Poisson integrators can be constructed by using
the splitting method. Numerical results show that the Poisson integrators
outperform the higher order non-Poisson integrators in phase orbit tracking,
long-term energy conservation and efficiency.
