Visibility of Luna and Lunar Crescent

+alp3'pow(!4)/100000000.0)rPI/180.0);
)
fo{i=oji?44;i+f)
{ fed(fpa2,"o4lf%lf%lf',&v,&rip0,&alpl);
smlT=s@vpf v'sitr((alpot.lp l'l)'Pvl 80.0);
)
fo(i=oii<154;i+r)
I fsc f(tptt2, r6lf %lfo/o1t',&v,&alp0,&alPi );
sumlpp=sumvpfi v's n((alpo+llPlil)iPI/l 80 0)i
)
fo(i=0ii<25;i++)
I
lsconi(fpu2, 9'lf ?'lf eolf'.&v.&slP0.&alpI)i
s!mwpp-smvppF ! | sin(alpo-alp I i t)t PVI 80 0li
)
mton{=218.11665+481267.88134't-13.2581'l/10000.0
+ 1 .856'po*(!3u 000000.0-1 .514'pos{tjy1000000{0 0
+sum!4Gur?+swvppl+sMvppp't'i/l 0000.0)/1000 0i
nrois=(nloD936o.o-(long (mlo.s,/360.0)))1360 0;
sumu=o,0tsumup=o.0rsunupp{,0isunuPpp=o 01
sunr=o.0;sumrFo.0isumrtp=o.0;suntppp-o.0i
lb(i=0ri<l88ri++)
r iscanfi fok2.'/'lf t'll"/.lf %lf %lf
' %lf ,&v.&.1p0,&alp l ,&alp2 &alpl.&'lp4 ):
!mu+mu+vt sin(alPo-aipl 'r falp2'r't/l0000 0'alpJ'pov'1t.J y' t000000 0
+alpl.pov(!4y100000000-0)'Ptl 80 0)i
)
fodi=0;i<64ii++)
{ f$anf( lol12,i%lf'l"lf e"lf .& v,&tlpo.&olP I )
suNup=sumufi! rsml( tlpo+alp l'r )r Pl/ l 80 0)
)
r fscanll tDk2. '/"la
%lf "/.1f ,&v,Aalpo.&alp | ):
sumuip=smupp+vtin(alPo+alpI'0'PI/l 80 0li
)
I l*anfi fDr2."%lf %ll c,6lf ,&t,&dpo.&dp | )
sunu;pp'sumuppF vrnn(.lpo'llpl'0'Pl/l 80 0):
I
;hFsumu+GuouCsunupplt+smupPp1irVl0000.0yl000 0;
fo(i-0;i<l54ii++)
I
Gcsnf(fplr2, %lf "/.lf %lr9"lf 90lt
9"lf .&v.&alD0.&alDt,&.102,&alp3,&dlp4)i
snrsmrr'cosriltpO alpt'l+;lp2'rt/lO0o0o' Jpt'po*rt l/10000000
+alp3xpowt,4yl00000000.0)'Pv180.0)r
2t5