Expansion of the Hamiltonian of the two-planetary problem into the Poisson series in all elements: Application of the Poisson series processor

This paper is the third in a series of articles devoted to one of the basic problems of celestial mechanics: the study of the evolution of solar-type planetary systems. In the previous papers a brief review of the history and current state of the problem was given; the plan of the study was outlined; the Jacobi coordinates and the related osculating elements were introduced; the form of the Poisson expansion of the Hamiltonian in all elements was given; and the expansion coefficients for the Hamiltonian of the two-planetary Sun-Jupiter-Saturn problem were obtained (though with impure accuracy) by a simple algorithm that is reduced to the calculation of multiple integrals of elementary functions. In the present paper the expansion of the Hamiltonian of the two-planetary Sun-Jupiter-Saturn problem into the Poisson series in all elements is constructed with the help of the PSP Poisson series processor, which is capable of required accuracy.