Вот что приводит Meeus:
E = L0 - 0.0057183d - alpha + Dpsi*cos(eps)
где L0-средняя долгота солнца. Используя VSOP87, получаем:
L0 (градусы)=
= 280.4664567 +
+ 360007.6982779*tau +
+ 0.03032028 * tau^2 +
+ tau^3/49931 -
- tau^4/15300 -
- tau^5/2000000
tau = (JD-2451545.0)/365250
(внимание! tau в этой формуле в тысячелетиях, а не столетиях)
alpha - видимое прямое восх. солнца, вычисленное с учетом аберрации и нутации.
Dpsi*cos(eps)
Dpsi-нутация по долготе, eps-наклонение эклиптики.
Это точная формула (определение ур. времени). Ссылаясь на Smart, Meeus приводит и приближенную формулу:
E =
y*sin(2*L0) -
- 2*e*sin(M) +
+ 4*e*y*sin(M)*cos(2*L0) -
- 0.5*y^2*sin(4*L0) -
- 1.25*e^2*sin(2M)
где y=[tg(eps/2)]^2
L0-ср. долгота солнца (см. выше)
e-эксцентриситет земной орбиты
M-ср. аномалия солнца
eps= 23d26'21.448" -
- 46.8150" *T -
- 0.00059 *T^2 +
+ 0.001813 *T^3
e = 0.016708634 -
- 0.000042037*T -
- 0.0000001267*T^2
M (градусы) =
357.52911 +
+ 35999.05029*T -
- 0.0001537*T^2