Visibility of Luna and Lunar Crescent

mooi coord(i),
llii'i--.Ji. i:s'pv r ro.or-'i"olat'Pvr 80 o
''srd
mder@P'Pr't I 80 0):
dem<os{plaf Pl/l 80 Of ds{mdelcp'Pv l 80 0):
ifffabslnw/deM)'= I o)
{ hass=a@s(num/domrrl80o/PI:
'
Dnor{alph.Fplong/15.0_gnseeroi
mone{nol-has/15 0;
mtwr.fuol+h6J15.0;
i(rndol<o 0) nnol+=24 0;
i(mtoP24 0) mnot-=24 0i
i(mone<o 0) nonar=24.0;
i(mone24 0) nonc'=24 0i
i(mrivo<o 0) ntqcF24.oi
i{mt{o>24 0)mt*o-24 0:
sts!=sstzeGrmlm' 1.00273?90935i
mlwoattu+deltarlin(lve{-1620 0)l/3600 0;
jd=jdiandat (ye&,mo hdlt€,mtwoz4 0);
r=(jd-24s1545.0)/36s25.0;
eph=qrs+plone/l 5.0-malPhoP,
dlr t -6rn(simelar'Pv180.0)'5in(ndelup'Pvl80 0)
_
c;.{pisr Pvl80 0fLo)(nddlap'Pul80 (
'cos(eph'15 o'Pvl80 0))'l80 cvPl
d.linn=(ali5e0.l25y(cos(plat'P/180 0)
'cos(ndelap'Pvl80 0)'sin(@ph'15 o'Ptl800))y15 0:
!ts=mtwo+dellaml
I
uhs+(acos(cruj(eMl+h,td l000.0J f IE0 0/Pl/l5 0)i
)
ltns=utns+aniin;
i(xms<0.o) ltmsts24.o;
i(1lme24.0) llme=24.0;
for(lcnFl Jcnr<4Jcnt++)
{
jd-juliedate(y@,monthdat€,fiac)i
jd.=jd+d.lhrlin(y.d-I620.0)l(3600 0'24 0);
stdjuliddaG8q,monthd.lc,0.0);
st=Gtjd.2451545.0Y16525.0i
t=(ide-2451545.0Y3652s.0;
nllrsin(o.125'PV180.o)-sin(llatrPVl80.0)'sin(ndctaprPVl80 0);
226