read("J2000-sol.gp"); \\ v - u où v est l'anomale vraie et u l'anomalie excentrique VmoinsU(u,e) = asin((e-(1-sqrt(1-e*e))*cos(u))*sin(u)/(1-e*cos(u))); \\ Calcul de t connaissant l pour un passage au périgée caractérisé \\ par tO et lO. tpourl(l,e,lO,tO) = { local(u,v,i); v = l-lO; u = v; for(i=1,20,u = v-VmoinsU(u,e)); return((u-e*sin(u))*365.25/2/Pi+tO); } \\ Instant du début de la saison, A l'année et l (un multiple de Pi/2) \\ la longitude caractérisant la saison. SAISON(A,l) = { local(jd,tO,lO,e); jd = JD2000(1,1,A,0,0,0); tO = TempsPerigeeSoleilM(A); lO = LaEvSoleil((jd+tO)/36525); e = ExcentriciteSoleil((jd+tO)/36525); return(CD(tpourl(l,e,lO,tO)+JD(1,1,A,0,0,0))); } printemps(A) = SAISON(A,2*Pi); ete(A) = SAISON(A,2.5*Pi); automne(A) = SAISON(A,3*Pi); hiver(A) = SAISON(A,3.5*Pi);