Для четырехпланетной задачи построена осредненная численно-аналитическая теория движения второго порядка по массам планет. В работе рассматриваются результаты численного интегрирования осредненных уравнений движения для системы Солнце-Юпитер-Сатурн-Уран-Нептун на интервале времени 10 млрд лет. Интегрирование проводилось методами Рунге-Кутты и Эверхарта высоких порядков. Движение планет имеет почти периодический характер. Эксцентриситеты и наклоны орбит планет остаются малыми. Короткопериодические возмущения сохраняются малыми на всем интервале интегрирования. Сделаны выводы о резонансных свойствах движения. Даны оценки точности численного интегрирования.