1 REM  ESTACIONES    ------- SAVE ESTACION ------
18 CLS:PRINT "*  SOLSTICIO/EQUINOCCIO  *
20 PRINT :INPUT "ENTRE EL AÑO COMO YYYY";Y
30 PRINT
80 JD#= 1721139.2855# +365.2421376#*Y  +.067919*(Y/1000)^2 
     -.0027879*(Y/1000)^3
90 T1#=JD#:GOSUB 350
100 DY=INT(GD):HR=INT((GD-DY)*24):MN=((GD-DY)*24-HR)*60
110 PRINT Y;" EQUINOCCIO DE MARZO     = ";DY;" DIA ";HR;" HRS ";MN;"MIN  GMT"
120 PRINT
130 JD# =1721233.2486# +365.2417284#*Y -.053018*(Y/1000)^2 
           +.009332*(Y/1000)^3
140 T2#=JD#:GOSUB 350
150 DY=INT(GD):HR=INT((GD-DY)*24):MN=((GD-DY)*24-HR)*60
160 PRINT Y;" SOLSTICIO DE JUNIO      = ";DY;" DIA ";HR;" HRS ";MN;
        " MIN  GMT"
170 PRINT
180 JD# =1721325.6978# +365.2425055#*Y -.126689*(Y/1000)^2 
       +.0019401*(Y/1000)^3
190 T3#=JD#:GOSUB 350
200 DY=INT(GD):HR=INT((GD-DY)*24):MN=((GD-DY)*24-HR)*60
210 PRINT Y;" EQUINOCCIO DE SEPTIEMBRE= ";DY;" DIA ";HR;" HRS ";
       MN;" MIN  GMT"
220 PRINT
230 JD#=1721414.392# +365.2428898#*Y -.010965*(Y/1000)^2 
        -8.488501E-03*(Y/1000)^3
240 T4#=JD#:GOSUB 350
250 DY=INT(GD):HR=INT((GD-DY)*24):MN=((GD-DY)*24-HR)*60
260 PRINT Y;" SOLSTICIO DE DICIEMBRE  =";DY;" DIA ";HR;" HRS ";MN;
       " MIN  GMT"
270 PRINT
280 PRINT "DURACION DE LAS ESTACIONES :"
285 PRINT
290 PRINT "DURACION DE LA PRIMAVERA =";:
       PRINT USING "###.####";T2#-T1#;:PRINT"DIAS":PRINT
300 PRINT "DURACION DEL VERANO=";:
       PRINT USING "###.####";T3#-T2#;:PRINT"DIAS"
305 PRINT
310 PRINT "DURACION DEL OTOÑO=";:
       PRINT USING "###.####";T4#-T3#;:PRINT"DIAS"
315 PRINT
320 PRINT "DURACION DEL INVIERNO=";:
       PRINT USING "###.####"; T1#+365.2422#-T4#;:PRINT"DIAS"
325 PRINT :INPUT "OTRO CALCULO ";A$
327 IF A$="s" OR A$="S" THEN 1
330 END
340 REM "*  DIA JULIANO A FECHA CALENDARIO GREGORIANO  *
350 JZ#=JD#+.5:JI=INT(JZ#):JF=JZ#-JI
360 IF JI<2299161# THEN JA=JI
370 IF JI=>2299161# THEN AL=INT((JI-1867216.25#)/36524.25#):JA=JI+1+AL-INT(AL/4)
380 JB =JA+1524 :JC=INT((JB-122.1)/365.25): 
        KD=INT(365.25*JC):JE=INT((JB-KD)/30.6001)
390 GD=JB-KD-INT(30.6001*JE)+JF
400 IF JE<13.5 THEN GM=JE-1 ELSE GM=JE-13
410 IF GM >2.5 THEN GY=JC-4716 ELSE GY =JC-4715
420 RETURN

2)Aplica el programa a 1248 . En dicho año observa :

2a.- Por no haber hecho todavía la reforma del calendario y ser el año juliano de duración un poco mayor que el año real, las fechas de comienzo de las estaciones se iban retrasando.

2b.-En este año las duraciones de Primavera y Verano eran iguales, mientras las de otoño e invierno también. Explicar a que se debe la variación de la duración de las estaciones y que causó esta coincidencia.