Visibility of Luna and Lunar Crescent

)
nmFsinc(so 0/60 of?V l 80 0)-
sin(;lar'Pvl80 o)rsin(sd€l$prPVl 80 0);
dem=o(iL'Pvl80 0)t@s(sdel6P'Pvl80 0)i
iflfab(nud/d.M)<= 1.0)
I
hN=@s(.w/d@)rl800/Pl;
mot=sdphaFplonS/l 5 o-gnstzdoi
nontmot'h4Yl5 0;
ntwetuot+hlsJl5,0;
i(dot<o 0) nnot+=24 0i
i(l@1>24-0) nnor*24 0;
i(mon?<o 0) mooe+=24 01
i(oon>z4 0) mone'=24 0;
i(nt*o<0.o) mtwo+=24.0i
i(olwo>24.0)nt{o-=24 0;
slsrgmsldo+mon tl 00273?90935:
none=mone+d€ltarlin(lyed-1620.0)l/3600 0i
jd=juliedate(y.e,6odhdat,mo.e1240);
F(jd-2451545-0Y365250,0a
su._coodc):
@ph=srsr+plotrg/l 5.oelphap;
altsFbincid(Pl51Pvl80.0)isinGdellap'Pt/l 80 0)
+cos(! lat'Pl/ 1 80.0)'cos(sdellap'Pv l 80 0)
'@s(qph' 15.01PUI80 0))' I 80 o/Pli
ddlan{(all.$50.0/60.o/(cos(plat'PI/180.0)'
.os(sdehap'PI/l 80.0Xsin(€phx l5-04PU180.0))yl 5.0;
ulsFhone+delrmi
)
ulsFursFdellallin(ly.aFl 620.0)V3600.0i
uts!+=(aco(e@y(eoaj+hite/1000 0))' I 80.0/PD/15.0)i
rac=drsd24.0;
,
llsR$r+4nlrm;
i(hsr<o.o) hs.F24-0i
i(lts!>24.0) Itsr:24.0;
fo(lcnt=l Jcnt<4jcil+)
{
jd-j'iiedardyeu,nonrhdate,f.e);
jdFjd+d€llarlini(yed-1620.0)l(3600.0'24-0);
nutarion0;
stjd=juliedlle(]dJnonthdat,0.0)i
st=Grjd-245 t 545.0y3652t.0;
sid-1imc(sr)i
t=(jde,245t545.0y36525.0;
225