Desiderei sapere come è possibile calcolare la fase della luna di un qualsiasi giorno di un qualsiasi anno.

Per calcolare
il moto della Luna con precisione sono necessarie centinaia di termini
dipendenti dai
Ecco la formula:

i = 180° – D – 6°,289 sin M’ + 2°,100
sin M – 1°,274 sin (2D–M’) – 0°,658 sin
2D – 0°,214 sin 2M’ – 0°,110 sin D

Dove gli
angoli D, M e M’ si ricavano dalle formule seguenti:

Elongazione
media della Luna
D = 297,8502042 + 445.267,1115168 T – 0,0016300 T² + T³/545.868
– T^4/113.065.000

Anomalia
media del Sole
M = 357,5291092 + 35.999,0502909 T – 0,0001536 T² + T³/24.490.000

Anomalia
media della Luna
M’ = 134,9634114
+ 477.198,8676313 T + 0,0089770 T² + T³/69.699 – T^4/14.712.000