An analytic solution for the J2 perturbed equatorial orbit

An analytic solution for the J2 perturbed equatorial orbit is obtained in terms of elliptic functions and integrals. The necessary equations for computing the position and elocity vectors, and the time are given in terms of known functions. The perturbed periapsis and apoapsis distances are determined from the roots of a characteristic cubic.