Visibility of Luna and Lunar Crescent

mcela=sin(!lar*PI/180.0)'sin(mdcliaprPvl80.0)1si.(mall'PV180.0);
nft la=merdcos(ndellarPvl80.0)'@s(mlt'PUl80 0));
i(fah6(nela)<= 1.0)
I goroxy(45,20)i@ut<<"
meeta=acos(m€ela)' 180.(VPIi
mdehFnddiap+@t'1cos{mera*Pvl 80 0);
(@t'sii(6ela'rvl80.oycos(mdehe'Pl/180 0)yls 0i
malFsin(plarrPvl80.0)'sin(ndehd'Pl/180.0)+cos(plal'PVl80 0)*co(hdcllai
?v180 0)rcos{mhdll 5.0'PYl80 0)i
nralt=ain(mah)'l 80 o/PI;
nalt+=(aco(eMj(enej+hik/1000 0))x 180.0/PI);
]
)
cour<<'.eb geater lhan I i
{ hdeluFndellapi
; a,m-si1(m,lelrar'Pl lSOo)_s,n{plat'Pl l80ofs'n(moh'Pl l800r:
nDzm=mdnv(cos(olaf PVI 80 0lrcos(mah'PVI 80 0)li
nd=lcodmaa)r 180 o/PI;
i(nha-12.0)
nazm=160.0:
i(mla<12-0)
mm=1600_n@;
nrmidia=(n6<Vmdsr)' 180 o/PI;
..r,r,="."u.otrtio"o'pr rtOOttsrn(eps'lon'Pl l80 0)rpos(l 0_
po$,lsi;(epcilonrPl,I80 0)'rn(slong'Pl/l80 0' 2 0) 0 t)))
ii;i'iriio",tr.o.po",rsinrcp'it.n'Pvl80 0)rs'n{slon8'PVl 80 0) 2 0) 0 5)li
iii=.co(iii)'l 80.0/PI;
mele hor=90 o.Dhil.mdai
aii*rir.:-.*.iai*z zf oqplatrPl/l 80 0ycos{l)hil'ddt'Pvl80 0):
rda; l=l0 t'.os(plal'PVl80 0,/.osohilcmda'PVl80 0)i
217