Visibility of Luna and Lunar Crescent

]
fo(i=o;i<l l4ii+)
{ fs.d(fpL2,'%li %lf %lf'3v.&rlp0,&alpl );
suotFsuDr/*rcos((alpolalPl rl)'?tl 80 0).
)
fo(i=oii<68;i+)
I
fscdf(fp$2,'/"lf o/.lf %lf'.&v.&alp0,&alpI):
vmDp=suffpPrvtosl alp0+tlp l rt l'PI/ l 80 0),
l
fo(i=oii<9:i+)
r fsdfi fDk2."!/olf./olf %lf'.& v.&tlPo.&.lp l)i
)unrppp runrtpp{'coe(dlpo'alpl't'Pl 180 oli
)
md{-li5O0O.5 7+sunr sudrsr suffpp'l_sumrPpptl't
| 0000 0:
nla{=(ml5r/l60.o{lonc (hla/160 0)))'160 0i
'"i.""',-O OOOrolla-O.O0OOt019'sin(225 O'47?lo8 o'lf PUI80 0)):
mhr+-={.d o0o0l754'sin(l83 1+481202 0r0'Pl/! 80 0))i
ndst+=(o o708rcos((225.0+a?7198 e'rl'PVl 80 0))
maipha=sin(mton8'PUlSOofcskFilon'Pvl800)i - -
nrn;ha- malDhard(mlal'Pt/l 80 0t'sin(cPsrlon'PUl80 0)'
mal;ha=ialpha(@(mlong'u/180.0));
maloha=.iu(maLlha)+ I 80 o/PI
-ri,,ir.no'go oraai.r.ng< u o o I I mdlpha- _1800:
iftmbn;27o or nalPha+-1600:
if(mh*24.0)
frba_-24 0;
inmhi<o.or dha+=24 0i
m;?lk=$n(mltr I Pv l80 ofcosl€psrlorr PV lao 0): --^ ^ -
-i"iu=na.r. -'Lmr"t'pvrtOO)tsin('psilontPVlSO O)rsininlonsrPl/1800):
nrdelh=asin(mdella)'l 80 o/Pl:
rulphlp=malpha+dell4lPhai
mha=hh.nelraatPha;
-lii 'iil""iJ" er'i do ol'"'t.aerl!pr Pr/ r 80 o,
"os(l)rar'Pv
r80 0)'co(mderbP'
Pl/l8o.o)'cos(mnar lt 0rPyl80 0l,
malr=6in(malt)1 180 o/Pli
i(sr=l)
;rcfr( 1.oz(bn(nal(+ l0 la mah s lt)fPl/l80oDft283 042730 remfl))'(tr
drcfFmefr/60 0;
216