Numerical integration of periodic orbits in the elliptic restricted three-body problem

The problem of finding periodic orbits in the circular restricted three-body problem has been very extensively studied in celestial mechanics. It is well known that continuous families of periodic orbits exist, for which the period varies in a continuous way. However, all the applications which are found in the solar system correspond to cases with non-zero eccentricities, and the elliptic restricted three-body problem is thus a better approximation than the circular one. For instance, for the motion of a satellite in the Earth-Moon system, as a first approximation, we may assume that the Moon moves around the Earth in circular motion; but as a much better approximation, we can also assume that the Moon moves in an elliptic orbit around the Earth.