На примере задачи движения спутника реализуется метод приближенного аналитического вычисления предельных циклов трёхмерных систем обыкновенных дифференциальных уравнений. Оценки точности полученных периодических решений даются в терминах относительной невязки. Строятся траектории последовательных аналитических решений, и демонстрируется их сходимость к искомому устойчивому предельному циклу, причём геометрическая форма данного цикла оказывается близкой к эллипсу.