Print "Enter date
Input M
Input D
Input Y
Gosub DJ
M=(J+4.867)/29.53058
M=2*(M-int M)-1
about=abs M
Print "Moon lit
Print about
Print "Headed for a
If M>0 Goto full
Print "new moon.
End 
Label full
Print "full moon.
End 
Label DJ
J=int (365.2422Y+30.44*(M-1)+D+1)
m₃=1
If M<3 Goto m₃
m₃=0
Label m₃
N=M-2+12*m₃
Z=Y-m₃
E=int (Z/100)
Z=Z-100E
W=int (2.61N-.2)+D+Z+int (Z/4)+int (E/4)-2E
W=W-7*int (W/7)
X=J-7*int (J/7)
xw=1
If X<W Goto xw
xw=0
Label xw
J=J-X+W-7*xw+1721061
Return