Mapping orbits with low station keeping costs for constellations of satellites based on the integral over the time of the perturbing forces