We investigate the long-time stability of the Sun-Jupiter-Saturn-Uranus system by considering a planar secular model, that can be regarded as a major refinement of the approach first introduced by Lagrange. Indeed, concerning the planetary orbital revolutions, we improve the classical circular approximation by replacing it with a solution that is invariant up to order two in the masses; therefore, we investigate the stability of the secular system for rather small values of the eccentricities. First, we explicitly construct a Kolmogorov normal form, so as to find an invariant KAM torus which approximates very well the secular orbits. Finally, we adapt the approach that is at basis of the analytic part of the Nekhoroshevs theorem, so as to show that there is a neighborhood of that torus for which the estimated stability time is larger than the lifetime of the Solar System. The size of such a neighborhood, compared with the uncertainties of the astronomical observations, is about ten times smaller.