  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.
