Approximate integrals of motion

We determine approximate numerical integrals of motion of 2D symmetric Hamiltonian systems. We detail for a few gravitational potentials the conditions under which quasi-integrals are obtained as polynomial series. We show that each of these potentials has a wide range of regular orbits that are accurately modelled with a unique approximate integral of motion.