Extrasolar systems with planets on eccentric orbits close to or in mean-motion resonances are common. The classical low-order resonant Hamiltonian expansion is unfit to describe the long-term evolution of these systems. We extend the Laplace-Lagrange secular approximation for coplanar systems with two planets by including (near-)resonant harmonics, and realize an expansion at high order in the eccentricities of the resonant Hamiltonian both at orders one and two in the masses. We show that the expansion at first order in the masses gives a qualitative good approximation of the dynamics of resonant extrasolar systems with moderate eccentricities, while the second order is needed to reproduce more accurately their orbital evolutions. The resonant approach is also required to correct the secular frequencies of the motion given by the Laplace-Lagrange secular theory in the vicinity of a mean-motion resonance. The dynamical evolutions of four (near-)resonant extrasolar systems are discussed, namely GJ 876 (2:1 resonance), HD 60532 (3:1), HD 108874 and GJ 3293 (close to 4:1).

We study the secular dynamics of extrasolar planetary systems by extending the Lagrange-Laplace theory to high order and by including the relativistic effects. We investigate the long-term evolution of the planetary eccentricities via normal form and we find an excellent agreement with direct numerical integrations. Finally we set up a simple analytic criterion that allows to evaluate the impact of the relativistic effects in the long-time evolution.

Mercury is the unique known planet that is situated in a 3:2 spin-orbit resonance nowadays. Observations and models converge to the same conclusion: the planet is presently deeply trapped in the resonance and situated at the Cassini state $1$, or very close to it. We investigate the complete non-linear stability of this equilibrium, with respect to several physical parameters, in the framework of Birkhoff normal form and Nekhoroshev stability theory. We use the same approach adopted for the 1:1 spin-orbit case with a peculiar attention to the role of Mercury's non negligible eccentricity. The selected parameters are the polar moment of inertia, the Mercury's inclination and eccentricity and the precession rates of the perihelion and node. Our study produces a bound to both the latitudinal and longitudinal librations (of 0.1 radians) for a long but finite time (greatly exceeding the age of the solar system). This is the so-called effective stability time. Our conclusion is that Mercury, placed inside the 3:2 spin-orbit resonance, occupies a very stable position in the space of these physical parameters, but not the most stable possible one.

We study the secular evolution of several exoplanetary systems by extending the Laplace-Lagrange theory to order two in the masses. Using an expansion of the Hamiltonian in the Poincar\'e canonical variables, we determine the fundamental frequencies of the motion and compute analytically the long-term evolution of the keplerian elements. Our study clearly shows that, for systems close to a mean-motion resonance, the second order approximation describes their secular evolution more accurately than the usually adopted first order one. Moreover, this approach takes into account the influence of the mean anomalies on the secular dynamics. Finally, we set up a simple criterion that is useful to discriminate between three different categories of planetary systems: (i) {\it secular} systems (HD 11964, HD 74156, HD 134987, HD 163607, HD 12661 and HD 147018); (ii) systems {\it near a mean-motion resonance} (HD 11506, HD 177830, HD 9446, HD 169830 and $\upsilon$ Andromedae); (iii) systems {\it really close to} or {\it in a mean-motion resonance} (HD 108874, HD 128311 and HD 183263).