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