Visibility of Luna and Lunar Crescent

iisha<o.0)
sdelb-sinGlarrPvl80-0f cos(epsilon'PI/180 0)l
sdelta=ldclls+cos(sla1+ PV l 80.O)'sin(€ps ilon'Pv l 80 0)'si.Glons'PI/ 1 80 0) i
d€lta=6inGd€lla)' I 80.o/Pl
salphap=elpha+.lcladpha
elr- sin$laf P! I 80.o sinrsdelrap'PI/ I 80.0 ) rosulat'PL I 80 0r.on sdehtp'PV
I 80.0)rcos(sba*15.0'Pvl " 80.0):
salt=asinGalt)' I 80.0/PI;
i(sw=l)
{
;Ffra l .oz(r&(eh+l o.v(sali+5.1 l ))'Pvl80 0)fe8l 0(271 0+teoF))!(pr€s
/l0l0.o);
l
I
s.cb-nn{olarrPul80 0Is'n(!d.l!lp'PU180 0,'sinleli'Pv!80 0):
tela=$etd(cos(cdelbp'PVl 80 0/'(olleht PLI 80 0r,:
seela=acos(seela);l 80 0/Pl;
sdcltrsde liap+srcfi'cos{seela * P l/ 1 80 0);
slDr=sh!-(sefi*sin(seetai Pl/ I 80.o)ho(sdelta/u/180 0)/ I 5'0i
sdehFsdelk: shFsha, l
idh=sinLplat'PUl80O)'iirsdelt,r.Pt/lSOOFcos(pldfPLlS00rrcoil3dcltat'PI/
180.0)'cos(shd'15.0'Pvl80 0):
satFdin(s.lr)' 180 o/PIi
*lr+=r,c6ermar4erur+tnle/1000 0D' 180 o/PL/l 5 0
sa6=s,n(rd.ld. pU I 80 O Fsid ptalr PI/ I 80.0)! sintsallr Pl/ | 80 0):
sm=szn/Gos(platrPvl80 o)'co(elt'Pl/180 0))l
sad=acos(saar)1 180 o/PI;
it(shdF=o 0)
ehe
s@=180 0i
itlshaF=12o)
el* i(shd<I20)
s*nidi.=Gn(ydsf l 80.0/PIi
)
void moon c@d(do!ble)
{
fDk2=foDcnf elP2000 ur'."r):
sumv=o.o,suntp=0 0lsumvtp=0 0isunvppp-0 0:
fo(i=0ii<2l8ii+)
I r(mfl rDu2.'%[ 0/"1r 0/0lf q.lf %lf
" 0lf ,& \,&!lpo'&!lp l.&alp2 &alpl &'ly'l
sMFsuov+v.sin((atpof atpl rr+alp2r1.r/10000.0
+alp3iPo*(i,3Y1000000 0
2t4