Вычисление периодических колебаний спутника

Рассматривается обыкновенное дифференциальное уравнение второго порядка, описывающее колебания спутника в плоскости эллиптической орбиты. Оно содержит два параметра: $e$ и $\mu$. При $0\leq e<1$ оно регулярно, а при $e=1$ имеет особенность Для $e=0,0.1,0.5,0.9,0.99$ и $0.999$ вычислены пять семейств симметричных (нечетных) периодических решений при $|\mu|\leq20$. Вычислены соответствующие значения следа, характеризующие устойчивость. При $e>0.9$ использована регуляризация с помощью эксцентрической аномалии. Результаты представлены графиками. Из них видно, что при $e\to1$ эти семейства стремятся к определенным предельным положениям.