\\ Julien est en l'honneur de Jules Scaliger, le père du mathématicien \\ qui a introduit cette échelle de temps en 1583. Le calendrier \\ grégorien date de 1582. \\ Calcul de la date en jours juliens JD(J,M,A,h,m,s) = { local(AT,BT,JDT) ; if (M<3 , A = A-1; M = M+12); AT = floor(A/100); BT = 2-AT+floor(AT/4); JDT = floor(365.25*(A+4716)) +floor(30.6001*(M+1)) +J +(h+(m+s/60)/60)/24 +BT -1524.5; return(JDT); } \\ Calcul de la date en jours juliens rapportée à l'époque 2000 \\ Utile pour l'utilisation des théories planétaires actuelles JD2000(J,M,A,h,m,s) = JD(J,M,A,h,m,s)-JD(1,1,2000,12,0,0); \\ Calcul de la date calendaire à partir du JD \\ La valeur retounée est le vecteur [J,M,A,h,m,s,JS] \\ (JS est le jour de la semaine) CD(jd) = { local(a,A,B,C,D,E,F,G,h,m,s,Z); Z = floor(jd+0.5); F = jd+0.5-Z; if (Z>=2299161, a = floor((Z-1867216.25)/36524.25); A = Z+1+a-floor(a/4), A = Z; ); B = A+1524; C = floor((B-122.1)/365.25); D = floor(365.25*C); E = floor((B-D)/30.6001); G = B-D-floor(30.6001*E)+F; h = (G-floor(G))*24; m = (h-floor(h))*60; s = (m-floor(m))*60; if (E<14,E=E-1,E=E-13); if (floor(E)>2,C=C-4716,C=C-4715); return([floor(G),E,C,floor(h),floor(m),floor(s),(Z+1)%7]); }