Visibility of Luna and Lunar Crescent

{ fs@f(fpt'%lf yolf y.lf,&sx,&bx,&c);
rmr=l€mF4.cd(bx+c.1);
)
slong+{tmp.po*{!i):
)
slonrslong/l 00000000.0;
slong=jon8. | 80.0/PIi
slong=Glon g/360.Odouble(longGlons/360.0)));
slongl=360.0 slone+-180.0
i(slone<0.0)
dong+=360.0;
i(slons>350.0)
slong=slong-360.0:
// L.firurL OF EARTH
fo(i=0:i<2;i+)
ro(=oj<adlilj+)
I [sco(fpr,"%lf yolf %lr,&ax,&bx,&c);
temp-aeosrdrco(bx+c0t
)
slai+=Cemprpov{t,t);
)
slat=*IaU100000000.0j
slatl-(t80.0/PI);
fo(i=0ri<sli+)
)
to{=0i<dsrflil;j+)
{ fs.n{&rr,"o/,lf %lf.4lf'.&a,&bx3c);
r.mpaettp+u'cos{bx+cn);
)
sdst=sdsl+renp*pow(,i);
sdsi=sdst/100000000.0r
dst=sdsr. 149597E70.0:
slo.g=slong+dclt!si{20.4898/sdrry3600.0i
elph!=sin(slons.Pvl80.0).6s{.Filon+Pyl 80.0);
slph!-slph!-tdGlal'Pvl 80.0).ein(.plilon.Pvl 80.0);
salpha=slphd(@s(slong.Prl 80-0));
lalpha=aran(salPha). I 80.0/PI:
i(Glong>90.0)&&Glorg<2?0.0)) sglphaFl8o.oi
if(slong>2?0.0) salpha+=150.0i
elpha-salphe/15.0;
i(s.lph*24.0)
sha=lastcur-salph!;
salpha=elph.'24.0:
2t3