A system of averaged equations of planetary motion around a central star is constructed. An astrocentric coordinate system is used. The two-planet problem is considered, but all constructions are easily generalized to an arbitrary number N of planets. The motion is investigated in modified (complex) Poincarécanonical elements. The averaging is performed by the Hori–Deprit method over the fast mean longitudes to the second order relative to the planetary masses. An expansion of the disturbing function is constructed using the Laplace coefficients. Some terms of the expansion of the disturbing function and the first terms of the expansion of the averaged Hamiltonian are given. The results of this paper can be used to investigate the evolution of orbits with moderate eccentricities and inclinations in various planetary systems.