Computing accurate Poincaré maps

We present a numerical method particularly suited for computing Poincaré maps for systems of ordinary differential equations. The method is a generalization of a stopping procedure described by Hénon [Physica D 5 (1982) 412], and it applies to a wide family of systems.