Движение вязкоупругой планеты

Необходимо выполнить следующие задания, с теорией
проверить равенства (2.8), вычислить скалярные произведения в правых частях системы уравнений (2.5) и получить замкнутую систему обыкновенных дифференциальных уравнений относительно переменных Андуайе описывающую вращательное движение вязкоупругой планеты на эллиптической орбите.
На следующем шаге планируется на базе этой системы уравнений получить усредненную систему уравнений в окрестности резонанса 1:1.