First order planetary perturbations with elliptic functions

The differential equations of planetary theory are solved analytically to first order for the two-dimensional case, using only Jacobian elliptic functions and the elliptic integrals of the first and second kind. This choice of functions leads to several new features potentially of importance for planetary theory. The first of these is that the solutions do not require the expansion of the reciprocal of the distance between two planets, even for those variables which depend on two angular arguments. A second result is that the solution is free from small divisors with the exception of two special resonances. In fact, not only are the solutions for resonant orbits free from small divisors, the perturbations for all variables are expressible in closed form. A subset of the resonant orbits maintains this form and in addition has the remarkable feature that the first order perturbations are purely periodic; they contain no secular terms. A solution for the 13 resonance case is given as an example.