Мат.модель движения Земли вокруг Солнца и солнечных часов
 

Есть модель движения Земли вокруг Солнца в виде системы уравненийX'' = -mu*X/|r^3|Y'' = -mu*Y/|r^3|Z'' = -mu*Z/|r^3|Интегрированием по времени с начальными уловиями для точки Афелия, где расстояние до Cолнца = 1,52*10^11 ; скорость = 2,92*10^4 можно решать эту систему для любого момента времени. И переводя координаты в нужную систему координат строить вектор солнечной тени на поверхности Земли. И все отлично считается и работает, но есть незадача, если рисовать суточную тень на 22 число каждого месяца, то насколько я знаю на Южном полюсе и на Северном тень должна совпадать в феврале (ю.п.) и апреле (с.п.), январе (ю.п.) и мае (с.п.). А у меня не совпадает и есть подозрение, что возможно в перигелий нужно было задаться новыми условиями или еще что-то я упустил. Не знаю даже.Подскажите, плиз, где копать