This is documentation for an obsolete product.

 8.4 The Lunation and LunationNumber Functions Lunation allows the dates of new moons to be addressed sequentially. Numbering new moons. The thousandth new moon since 1900 January 1 occurs on 1980 November 8. In[29]:=Lunation[1000] Out[29]= This confirms that the previous new moon date is the thousandth new moon. In[30]:=LunationNumber[%] Out[30]= Lunation[n], with n as an integer, gives the date of a new moon. If n is a half-integer, Lunation instead returns the date of a full moon. Thus Lunation[999.5] is the thousandth full moon since 1900 January 1. To produce a sequence of full moon dates, use LunationNumber to determine the starting value of n and then add 0.5 to get a full moon. Add any integer to get later full moons. An expression like Table[Lunation[LunationNumber[{1993,1,1}] - 0.5 + i], {i, 1, 13}], therefore, returns the dates of the 13 full moons that occur in the year 1993. It is also possible to determine the precise rising and setting times of the Moon with the functions MoonRise and MoonSet. Atmospheric refraction is taken into account in the same way as with the SunRise and SunSet functions. You can use the option Refract -> False to suppress refraction. MoonCalendar is a function that produces a text table showing the dates and zodiac positions of all the full and new moons in a given year. It also shows the dates when the Moon is at its halfway waxing and waning phases. Rising and setting times of the Moon. The Moon rises at 08:50 local time in Melbourne on 1993 November 17. In[31]:=MoonRise[{1993,11,17}] Out[31]= This shows that the first full moon of 1993 occurs on January 9 in the direction of the zodiac constellation of Gemini. In[32]:=MoonCalendar[1993] Out[32]=