We analyze the secular interactions of two coplanar planets based on a high order (order 12) expansion of the perturbative potential in powers of the eccentricities. The model depends on only two parameters (the ratio of semi-major axis and the mass ratio of the planets) and can be reduced to a one degree of freedom system, allowing for an exhaustive parametric analysis. Following Pauwels (1984) we map the phase space on a sphere, avoiding in this way the artificial singularities introduced by other mappings. We show that the twelve order expansion is able to describe correctly most of the exosolar planetary systems discovered so far, even if the eccentricities of these planets are considerably larger than the eccentricities of our own solar system. The expansion is even able to reproduce, at moderate eccentricities, the secular resonances discovered numerically by Michtchenko and Malhotra (2004) at moderate to large eccentricities.To search for other articles by the author(s) go to: http://adsabs.harvard.edu/abstract_service.html