A first-order, semi-analytical method for the long-term motion of resonant satellites is introduced. The method provides long-term solutions, valid for nearly all eccentricities and inclinations, and for all commensurability ratios. The method allows the inclusion of all zonal and tesseral harmonics of a nonspherical planet. We present here an application of the method to a synchronous satellite including J2 and J22 harmonics. Global, long-term solutions for this problem are given for arbitrary values of eccentricity, argument of perigee and inclination.
