function fluna()  {
//(Approssimazione alla 3a cifra)
pig=4*Math.atan(1)
ar=pig/180
eps=23.4392911*ar
giulia(gg,mm,aaaa,hh,mese)
Ts=(JD-2415020)/36525     //Periodo secolare
L1=270.434+481267.883*Ts  //Longitudine media della Luna
M=358.475+35999.049*Ts   //Anomalia media del Sole
M1=296.104+477198.849*Ts+.009*Ts*Ts   //Anomalia media della Luna
D=350.737+445267.114*Ts         //Elongazione della Luna
F=11.250+483202.025*Ts  //Distanza media della Luna dal suo nodo ascendente
E=1-.002*Ts
L1=L1-fint(L1/360)*360
M1=M1*ar
D=D*ar
F=F*ar
M=M*ar
i=180-D/ar-6.2*sin(M1)+2.1*sin(M)-1.2*sin(2*D-M1)
i=i*ar
K=(1+cos(i))/2            //Frazione illuminata del disco  Lunare

//Calcolo delle coordinate eclittiche
//e aggiunta dei termini periodici

//lambda
lam=L1+6.288*sin(M1)+1.274*sin(2*D-M1)+.658*sin(2*M)+.213*sin(2*M1)
lam=lam-.185*sin(M)*E-.114*sin(2*F)+.058*sin(2*D-2*M1)+.057*sin(2*D-M-M1)*E
lam=lam+.053*sin(2*D+M1)+.045*sin(2*D-M)*E+.041*sin(M1-M)*E-.034*sin(D)
lam=lam-.03*sin(M+M1)+.015*sin(2*D-2*F)
//beta
beta=5.128*sin(F)+.280*sin(M1+F)+.277*sin(M1-F)+.173*sin(2*D-F)
beta=beta+.055*sin(2*D+F-M1)+.046*sin(2*D-F-M1)+.032*sin(2*D+F)+.017*sin(2*M1+F)
beta=beta+.009*sin(2*D+M1-F)+.009*sin(2*M1-F)
lam=lam*ar
beta=beta*ar
//tanga=((sin(lam)*cos(eps)-Math.tan(beta)*sin(eps))/cos(lam))
//alflu=Math.atan(tanga)
//alert("lam="+lam+",beta="+beta)
//seno=sin(beta)*cos(eps)+cos(beta)*cos(eps)*sin(lam)
//declu=Math.asin(seno)
//alflu=alflu*12/pig
//if ((cos(lam))<0) alflu=alflu+12
//declu=declu/ar
//lam=lam/ar
//beta=beta/ar
//if alflu<0 then alflu=alflu+24
//ar1=alflu:dec1=declu
//ar=ar1*pig/12
//gosub zeta
//gosub trasfo
//dec=90-dec
//px=-dec*cos(ar):py=dec*sin(ar)
//px=px*scala:py=py*scala
}