Выполнено осреднение уравнений движения двупланетной задачи до второго порядка по возмущающим массам в гелиоцентрической системе отсчета. Использована комплексная форма канонических элементов Пуанкаре. Разложение возмущающей функции построено при помощи коэффициентов Лапласа. Осреднение выполнялось методом Хори Депри. Результаты настоящей работы легко обобщаются на случай исследования движения произвольного числа N планет.