A way to generate a model of secular perturbations for planetary orbits

The reciprocal distance between two material points that rotate around a central body in nonintersecting orbits is expanded and the results are presented. The expansion is obtained accurate to the tenth order with respect to small parameters: the eccentricities and sine of the orbital inclination angle. The result is the basis of the averaging operation of the perturbation function in the system of eight major planets in the solar system, and of the numerical integration of the averaged equations of motion. The averaged Hamiltonian contains the terms whose period of variation is greater than 200 years. Forty eight equation of first order are numerically integrated with increments of 100 years for two intervals from the beginning of the Christian era: 25 million years forward and 25 million years backward over time. To present the results of calculation, the website (URL: http://vadimchazov.narod.ru/secequat.htm) was developed, where the initial codes, executable program modules, the results of calculations presented in graphical form, text files with initial conditions, tables for expanding the reciprocal distance between two material points, and the tables with the results of expansion of the perturbation function for eight major planets of the solar system are presented.