Home » Archimedes archive » Micro User » MU 1990-08.adf » Aug\BBC/Eclipse

Aug\BBC/Eclipse

This website contains an archive of files for the Acorn Electron, BBC Micro, Acorn Archimedes, Commodore 16 and Commodore 64 computers, which Dominic Ford has rescued from his private collection of floppy disks and cassettes.

Some of these files were originally commercial releases in the 1980s and 1990s, but they are now widely available online. I assume that copyright over them is no longer being asserted. If you own the copyright and would like files to be removed, please contact me.

Tape/disk: Home » Archimedes archive » Micro User » MU 1990-08.adf
Filename: Aug\BBC/Eclipse
Read OK: ✔
File size: 169B bytes
Load address: FFFFFB43
Exec address: 6D421948
Duplicates

There is 1 duplicate copy of this file in the archive:

File contents
   10REM ** Lunar Eclipse Finder **
   20REM ** (c) S Hatch December 1989 **
   30*TV255,1
   40MODE7:VDU14:VDU23;8202;0;0;0;:FOR C%=0 TO 1:PRINTTAB(7,C%)CHR$141;CHR$131;"Lunar Eclipse Finder":NEXT C%:DIM month$(12):FOR C%=1 TO 12:READ month$(C%):NEXT C%
   50DATA January,February,March,April,May,June,July,August,September,October,November,December
   60REPEAT:PRINTTAB(11,8)"Which year? ":INPUTTAB(23,8)""Y%:Y$=STR$(Y%)
   70MODE1:PRINT''"Searching.....";:M%=1:D=1:PROCdays
   80REPEAT:PROCFIND:UNTIL yZ=Y%+1:*FX21
   90MODE7:PRINT'''"No more lunar eclipses in ";Y$''"<Q> to quit, any other key for more";:g%=GET OR 32:UNTIL g%=113
  100CLS:END
  110:
  120DEFPROCdays
  130m%=M%:y%=Y%
  140IFNOT (m%=1 OR m%=2) THEN y%=y%+1:m%=m%-12
  150y%=y%-1:m%=m%+12:A%=y% DIV 100:B%=2-A%+A% DIV 4:C%=INT(y%*365.25):DD%=INT(30.6001*(m%+1)):DAYS=B%+C%+DD%+D-723244
  160ENDPROC
  170:
  180DEFPROCFIND
  190PROCSUN:PROCMOON:difference=LAMBDASUN-LAMBDAMOON:check=ABS(difference):IF check<PI/2 THEN PROCdate:DAYS=DAYS+15:ENDPROC:ELSE PROCFULL
  200DAYS=DAYS+30
  210ENDPROC
  220:
  230DEFPROCSUN
  240N=(DAYS/365.2422)*2*PI:N=FNrt(N):M=N-.065674349:IF M<=0 THEN M=M+2*PI
  250E=FNKepler(M):LAMBDASUN=E+4.932237687:LAMBDASUN=FNrt(LAMBDASUN)
  260ENDPROC
  270:
  280DEFPROCMOON
  290l=1.13403578+DAYS*.2299715042:l=FNrt(l):Mm=l-.0019443683*DAYS-6.0978848:Mm=FNrt(Mm):EN=2.652035286-.00092422*DAYS:EN=FNrt(EN):Ev=.0222337493*SIN(2*(l-LAMBDASUN)-Mm):Ae=.0032428218*SIN(M):A3=.0064577182*SIN(M):Mtickm=Mm+Ev-Ae-A3
  300Ec=.1097567753*SIN(Mtickm):A4=.0037350046*SIN(2*Mtickm):ltick=l+Ev+Ec-Ae+A4:V=.0114895025*SIN(2*(ltick-LAMBDASUN)):ldoubletick=ltick+V:ENtick=EN-.0027925268*SIN(M):y=SIN(ldoubletick-ENtick)*.995970321:x=COS(ldoubletick-ENtick)
  310tan=ATN(y/x):tan=FNrtquad(tan):LAMBDAMOON=ENtick+tan:LAMBDAMOON=FNrt(LAMBDAMOON):BETAMOON=ASN(SIN(ldoubletick-ENtick)*.089683442)
  320ENDPROC
  330:
  340DEFFNrt(angle)
  350IF angle<0 THEN angle=INT(ABS(angle)/2*PI)*2*PI+angle+2*PI
  360=2*PI*(angle/(2*PI)-INT(angle/(2*PI)))
  370:
  380DEFFNKepler(ANGLE)
  390kl=ANGLE:REPEAT:delta=kl-.016718*SIN(kl)-ANGLE:correction=delta/(1-.016718*COSkl):kl=kl-correction:UNTIL ABS(delta)<=.00000001
  400=2*ATN(TAN(kl/2)*1.016860112)
  410:
  420DEFFNrtquad(tan)
  430IF y>=0 AND x>=0 THEN =tan
  440IF x<0 THEN =tan+PI
  450IF y<0 AND x>0 THEN =tan+2*PI
  460:
  470DEFPROCdate
  480DAYSS=DAYS+2444238.5:IZ=INT(DAYSS+.5):FZ=DAYSS+.5-INT(DAYSS+.5):IF IZ>2299160 THEN AZ=INT((IZ-1867216.25)/36524.25):ELSE AZ=IZ
  490BZ=IZ+1+AZ-INT(AZ/4):CZ=BZ+1524:DZ=INT((CZ-122.1)/365.25):EZ=INT(365.25*DZ):GZ=INT((CZ-EZ)/30.6001):dZ=CZ-EZ+FZ-INT(30.6001*GZ):ddZ=INT(dZ):hrs=DAYS-INT(DAYS):hr=INT(24*hrs):min=INT(60*(24*hrs-hr))
  500IF GZ<13.5 THEN mZ=GZ-1:ELSE mZ=GZ-13
  510IF mZ>2.5 THEN yZ=DZ-4716:ELSE yZ=DZ-4715
  520min$=STR$(min):IF LEN(min$)<2 min$="0"+min$
  530hr$=STR$(hr):IF LEN(hr$)<2 hr$="0"+hr$
  540date$=STR$(ddZ):month$=month$(mZ):year$=STR$(yZ)
  550digit$=RIGHT$(date$,1):IF digit$="1" suffix$="st" ELSEIF digit$="2" suffix$="nd" ELSEIF digit$="3" suffix$="rd" ELSE suffix$="th"
  560IF LEFT$(date$,1)="1" suffix$="th"
  570ENDPROC
  580:
  590DEFPROCFULL
  600REPEAT:catchup=8.8837259E-3+1.0471976E-3*COS(Mtickm):IF LAMBDASUN<PI THEN difference=difference+PI ELSE difference=difference-PI:IF ABS(difference)>PI AND SGN(difference)=-1 THEN difference=2*PI+difference
  610DAYS=DAYS+(difference/(24*catchup)):PROCSUN:PROCMOON:difference=LAMBDASUN-LAMBDAMOON
  620UNTIL ABS(ABS(difference)-PI)< 0.000001:PROCdate
  630test=ABS(ABS(ldoubletick-ENtick)-PI):IF test<.2138028334 OR test>2.92778982 THEN VDU7:PROC_diagram
  640ENDPROC
  650:
  660DEFPROC_diagram
  670rhotick=.99698599/(1+.0549*COS(Mtickm+Ec)):pi=.016592845/rhotick
  680deltalambda=.009599311+.00104719755*COS(Mtickm):deltabeta=8.726646E-4*COS(ldoubletick-ENtick)
  690Su=pi-.00471238898:thetam=.009042551/rhotick:rate=deltalambda-7.16783E-4:factor=256/rate
  700a=1+deltabeta^2/rate^2:b=2*BETAMOON*deltabeta/rate:H=(Su+thetam/2)^2:c=BETAMOON^2-H
  710root=b^2-4*a*c:IF SGN(root)=-1 PROC_pen:ENDPROC
  720X2=(-b+SQR(root))/(2*a):X1=(-b-SQR(root))/(2*a):IF SGN(X2)=-1 X=X1:X1=X2:X2=X
  730starthrs=X1/rate:endhrs=X2/rate:mid=24*hrs:starttime$=FN_DDtoDM(mid+starthrs):endtime$=FN_DDtoDM(mid+endhrs):mid$=hr$+":"+min$+" GMT"
  740CLS:CLG:VDU20
  750PRINT"Eclipse particulars are:"''TAB(5)"Date: ";date$;suffix$;" ";month$;" ";year$''
  760PRINTTAB(5)"Eclipse begins: ";starttime$'TAB(5)"Mid-eclipse:    ";mid$;'TAB(5)"Eclipse ends:   ";endtime$
  770*FX21
  780PRINT'''"Any key for the diagram";:G%=GET:VDU19,2,4,0,0,0:COLOUR130:GCOL0,0
  790CLS:MOVE0,512:DRAW1280,512:MOVE0,550:VDU5:PRINT"Ecliptic":MOVE1018,550:PRINTdate$;":";STR$(mZ);":";RIGHT$(year$,2):VDU4
  800r=Su*factor:PROC_circle(640,512,r)
  810lh=(2.5*deltabeta+BETAMOON)*factor+512:MOVE0,lh:DRAW1280,(BETAMOON-2.5*deltabeta)*factor+512:MOVE0,lh:DRAW50,(2.3*deltabeta+BETAMOON)*factor+530:MOVE0,lh:DRAW50,(BETAMOON+2.3*deltabeta)*factor+500
  820GCOL0,3
  830r=thetam/2*factor:PROC_circle(640-X1*factor,(BETAMOON+starthrs*deltabeta)*factor+512,r):PRINTTAB(33,15+12*SGN(BETAMOON))LEFT$(starttime$,5)
  840y=BETAMOON*factor+512:GCOL0,1:COLOUR1:PROC_circle(640,y,r):PRINTTAB(17,15+12*SGN(BETAMOON))LEFT$(mid$,5):GCOL0,3:COLOUR3
  850PROC_circle(640-X2*factor,(BETAMOON+endhrs*deltabeta)*factor+512,r):PRINTTAB(3,15+12*SGN(BETAMOON))LEFT$(endtime$,5):*FX21
  860PRINTTAB(10,31)"Any key to continue...";:g%=GET:PRINTTAB(10,31)"Search resumed........";
  870ENDPROC
  880:
  890DEFPROC_pen
  900mid=24*hrs:mid$=hr$+":"+min$+" GMT":CLS:CLG:VDU20
  910 PRINT'''"This is a penumbral eclipse:"''TAB(5)"Date: ";date$;suffix$;" ";month$;" ";year$''TAB(5)"Mid-eclipse:    ";mid$;:*FX21
  920PRINT'''"Any key for the diagram";:G%=GET:VDU19,129,132,0,0,0,19,2,0,0,0,0:COLOUR 129:GCOL0,0
  930CLS:MOVE0,512:DRAW1280,512:MOVE0,550:VDU5:PRINT"Ecliptic":MOVE1018,550:PRINTdate$;":";STR$(mZ);":";RIGHT$(year$,2):VDU4:r=Su*factor:PROC_circle(640,512,r)
  940lh=(2.5*deltabeta+BETAMOON)*factor+512:MOVE0,lh:DRAW1280,(BETAMOON-2.5*deltabeta)*factor+512:MOVE0,lh:DRAW50,(2.3*deltabeta+BETAMOON)*factor+530:MOVE0,lh:DRAW50,(BETAMOON+2.3*deltabeta)*factor+500:GCOL0,3
  950y=BETAMOON*factor+512:r=thetam/2*factor:PROC_circle(640,y,r):PRINTTAB(17,15+12*SGN(BETAMOON))LEFT$(mid$,5):*FX21
  960PRINTTAB(10,31)"Any key to continue...";:g%=GET:PRINTTAB(10,31)"Search resumed........";
  970ENDPROC
  980:
  990DEFFN_DDtoDM(time)
 1000LOCAL hrs,hrs$,mins,mins$
 1010hrs=INT(time):hrs$=STR$(hrs):IF LEN(hrs$)<2 hrs$="0"+hrs$
 1020mins=INT((time-hrs)*60+.5):mins$=STR$(mins):IF LEN(mins$)<2 mins$="0"+mins$
 1030=hrs$+":"+mins$+" GMT"
 1040:
 1050DEFPROC_circle(ox,oy,radius):LOCAL x,y
 1060MOVE ox-radius,oy:DRAW ox+radius,oy:FOR y=1 TO radius STEP 4
 1070x=SQR((radius+y)*(radius-y))
 1080oxpx=ox+x:oxmx=ox-x:oypy=oy+y:oymy=oy-y
 1090MOVE oxmx,oypy:DRAW oxpx,oypy
 1100MOVE oxmx,oymy:DRAW oxpx,oymy
 1110NEXT
 1120ENDPROC


 � ** Lunar Eclipse Finder **
%� ** (c) S Hatch December 1989 **
*TV255,1
(}�7:�14:�23;8202;0;0;0;:� C%=0 � 1:�7,C%)�141;�131;"Lunar Eclipse Finder":� C%:� month$(12):� C%=1 � 12:� month$(C%):� C%
2[� January,February,March,April,May,June,July,August,September,October,November,December
<0�:�11,8)"Which year? ":�23,8)""Y%:Y$=�(Y%)
F*�1:�''"Searching.....";:M%=1:D=1:�days
P�:�FIND:� yZ=Y%+1:*FX21
Ze�7:�'''"No more lunar eclipses in ";Y$''"<Q> to quit, any other key for more";:g%=� � 32:� g%=113
d�:�
n:
x
��days
�m%=M%:y%=Y%
�'� (m%=1 � m%=2) � y%=y%+1:m%=m%-12
�my%=y%-1:m%=m%+12:A%=y% � 100:B%=2-A%+A% � 4:C%=�(y%*365.25):DD%=�(30.6001*(m%+1)):DAYS=B%+C%+DD%+D-723244
��
�:
�
��FIND
�m�SUN:�MOON:difference=LAMBDASUN-LAMBDAMOON:check=�(difference):� check<�/2 � �date:DAYS=DAYS+15:�:� �FULL
�DAYS=DAYS+30
��
�:
�	��SUN
�BN=(DAYS/365.2422)*2*�:N=�rt(N):M=N-.065674349:� M<=0 � M=M+2*�
�AE=�Kepler(M):LAMBDASUN=E+4.932237687:LAMBDASUN=�rt(LAMBDASUN)
�
:

��MOON
"�l=1.13403578+DAYS*.2299715042:l=�rt(l):Mm=l-.0019443683*DAYS-6.0978848:Mm=�rt(Mm):EN=2.652035286-.00092422*DAYS:EN=�rt(EN):Ev=.0222337493*�(2*(l-LAMBDASUN)-Mm):Ae=.0032428218*�(M):A3=.0064577182*�(M):Mtickm=Mm+Ev-Ae-A3
,�Ec=.1097567753*�(Mtickm):A4=.0037350046*�(2*Mtickm):ltick=l+Ev+Ec-Ae+A4:V=.0114895025*�(2*(ltick-LAMBDASUN)):ldoubletick=ltick+V:ENtick=EN-.0027925268*�(M):y=�(ldoubletick-ENtick)*.995970321:x=�(ldoubletick-ENtick)
6}tan=�(y/x):tan=�rtquad(tan):LAMBDAMOON=ENtick+tan:LAMBDAMOON=�rt(LAMBDAMOON):BETAMOON=�(�(ldoubletick-ENtick)*.089683442)
@�
J:
Tݤrt(angle)
^3� angle<0 � angle=�(�(angle)/2*�)*2*�+angle+2*�
h%=2*�*(angle/(2*�)-�(angle/(2*�)))
r:
|ݤKepler(ANGLE)
�skl=ANGLE:�:delta=kl-.016718*�(kl)-ANGLE:correction=delta/(1-.016718*�kl):kl=kl-correction:� �(delta)<=.00000001
�=2*�(�(kl/2)*1.016860112)
�:
�ݤrtquad(tan)
�� y>=0 � x>=0 � =tan
�� x<0 � =tan+�
�� y<0 � x>0 � =tan+2*�
�:
�
��date
�uDAYSS=DAYS+2444238.5:IZ=�(DAYSS+.5):FZ=DAYSS+.5-�(DAYSS+.5):� IZ>2299160 � AZ=�((IZ-1867216.25)/36524.25):� AZ=IZ
�BZ=IZ+1+AZ-�(AZ/4):CZ=BZ+1524:DZ=�((CZ-122.1)/365.25):EZ=�(365.25*DZ):GZ=�((CZ-EZ)/30.6001):dZ=CZ-EZ+FZ-�(30.6001*GZ):ddZ=�(dZ):hrs=DAYS-�(DAYS):hr=�(24*hrs):min=�(60*(24*hrs-hr))
�"� GZ<13.5 � mZ=GZ-1:� mZ=GZ-13
�&� mZ>2.5 � yZ=DZ-4716:� yZ=DZ-4715
)min$=�(min):� �(min$)<2 min$="0"+min$
$hr$=�(hr):� �(hr$)<2 hr$="0"+hr$
.date$=�(ddZ):month$=month$(mZ):year$=�(yZ)
&sdigit$=�date$,1):� digit$="1" suffix$="st" �� digit$="2" suffix$="nd" �� digit$="3" suffix$="rd" � suffix$="th"
0 � �date$,1)="1" suffix$="th"
:�
D:
N
��FULL
X��:catchup=8.8837259E-3+1.0471976E-3*�(Mtickm):� LAMBDASUN<� � difference=difference+� � difference=difference-�:� �(difference)>� � �(difference)=-1 � difference=2*�+difference
bRDAYS=DAYS+(difference/(24*catchup)):�SUN:�MOON:difference=LAMBDASUN-LAMBDAMOON
l(� �(�(difference)-�)< 0.000001:�date
vWtest=�(�(ldoubletick-ENtick)-�):� test<.2138028334 � test>2.92778982 � �7:�_diagram
��
�:
���_diagram
�Brhotick=.99698599/(1+.0549*�(Mtickm+Ec)):pi=.016592845/rhotick
�]deltalambda=.009599311+.00104719755*�(Mtickm):deltabeta=8.726646E-4*�(ldoubletick-ENtick)
�\Su=pi-.00471238898:thetam=.009042551/rhotick:rate=deltalambda-7.16783E-4:factor=256/rate
�Wa=1+deltabeta^2/rate^2:b=2*BETAMOON*deltabeta/rate:H=(Su+thetam/2)^2:c=BETAMOON^2-H
�'root=b^2-4*a*c:� �(root)=-1 �_pen:�
�JX2=(-b+�(root))/(2*a):X1=(-b-�(root))/(2*a):� �(X2)=-1 X=X1:X1=X2:X2=X
ڇstarthrs=X1/rate:endhrs=X2/rate:mid=24*hrs:starttime$=�_DDtoDM(mid+starthrs):endtime$=�_DDtoDM(mid+endhrs):mid$=hr$+":"+min$+" GMT"
��:�:�20
�Q�"Eclipse particulars are:"''�5)"Date: ";date$;suffix$;" ";month$;" ";year$''
�`�5)"Eclipse begins: ";starttime$'�5)"Mid-eclipse:    ";mid$;'�5)"Eclipse ends:   ";endtime$
	*FX21
?�'''"Any key for the diagram";:G%=�:�19,2,4,0,0,0:�130:�0,0
X�:�0,512:�1280,512:�0,550:�5:�"Ecliptic":�1018,550:�date$;":";�(mZ);":";�year$,2):�4
 #r=Su*factor:�_circle(640,512,r)
*�lh=(2.5*deltabeta+BETAMOON)*factor+512:�0,lh:�1280,(BETAMOON-2.5*deltabeta)*factor+512:�0,lh:�50,(2.3*deltabeta+BETAMOON)*factor+530:�0,lh:�50,(BETAMOON+2.3*deltabeta)*factor+500
4�0,3
>~r=thetam/2*factor:�_circle(640-X1*factor,(BETAMOON+starthrs*deltabeta)*factor+512,r):�33,15+12*�(BETAMOON))�starttime$,5)
H[y=BETAMOON*factor+512:�0,1:�1:�_circle(640,y,r):�17,15+12*�(BETAMOON))�mid$,5):�0,3:�3
Rm�_circle(640-X2*factor,(BETAMOON+endhrs*deltabeta)*factor+512,r):�3,15+12*�(BETAMOON))�endtime$,5):*FX21
\L�10,31)"Any key to continue...";:g%=�:�10,31)"Search resumed........";
f�
p:
z
��_pen
�/mid=24*hrs:mid$=hr$+":"+min$+" GMT":�:�:�20
�z �'''"This is a penumbral eclipse:"''�5)"Date: ";date$;suffix$;" ";month$;" ";year$''�5)"Mid-eclipse:    ";mid$;:*FX21
�Q�'''"Any key for the diagram";:G%=�:�19,129,132,0,0,0,19,2,0,0,0,0:� 129:�0,0
�x�:�0,512:�1280,512:�0,550:�5:�"Ecliptic":�1018,550:�date$;":";�(mZ);":";�year$,2):�4:r=Su*factor:�_circle(640,512,r)
��lh=(2.5*deltabeta+BETAMOON)*factor+512:�0,lh:�1280,(BETAMOON-2.5*deltabeta)*factor+512:�0,lh:�50,(2.3*deltabeta+BETAMOON)*factor+530:�0,lh:�50,(BETAMOON+2.3*deltabeta)*factor+500:�0,3
�cy=BETAMOON*factor+512:r=thetam/2*factor:�_circle(640,y,r):�17,15+12*�(BETAMOON))�mid$,5):*FX21
�L�10,31)"Any key to continue...";:g%=�:�10,31)"Search resumed........";
��
�:
�ݤ_DDtoDM(time)
�� hrs,hrs$,mins,mins$
�5hrs=�(time):hrs$=�(hrs):� �(hrs$)<2 hrs$="0"+hrs$
�Gmins=�((time-hrs)*60+.5):mins$=�(mins):� �(mins$)<2 mins$="0"+mins$
=hrs$+":"+mins$+" GMT"
:
!��_circle(ox,oy,radius):� x,y
$4� ox-radius,oy:� ox+radius,oy:� y=1 � radius � 4
.x=�((radius+y)*(radius-y))
8+oxpx=ox+x:oxmx=ox-x:oypy=oy+y:oymy=oy-y
B� oxmx,oypy:� oxpx,oypy
L� oxmx,oymy:� oxpx,oymy
V�
`�
�
00000000  0d 00 0a 20 f4 20 2a 2a  20 4c 75 6e 61 72 20 45  |... . ** Lunar E|
00000010  63 6c 69 70 73 65 20 46  69 6e 64 65 72 20 2a 2a  |clipse Finder **|
00000020  0d 00 14 25 f4 20 2a 2a  20 28 63 29 20 53 20 48  |...%. ** (c) S H|
00000030  61 74 63 68 20 44 65 63  65 6d 62 65 72 20 31 39  |atch December 19|
00000040  38 39 20 2a 2a 0d 00 1e  0c 2a 54 56 32 35 35 2c  |89 **....*TV255,|
00000050  31 0d 00 28 7d eb 37 3a  ef 31 34 3a ef 32 33 3b  |1..(}.7:.14:.23;|
00000060  38 32 30 32 3b 30 3b 30  3b 30 3b 3a e3 20 43 25  |8202;0;0;0;:. C%|
00000070  3d 30 20 b8 20 31 3a f1  8a 37 2c 43 25 29 bd 31  |=0 . 1:..7,C%).1|
00000080  34 31 3b bd 31 33 31 3b  22 4c 75 6e 61 72 20 45  |41;.131;"Lunar E|
00000090  63 6c 69 70 73 65 20 46  69 6e 64 65 72 22 3a ed  |clipse Finder":.|
000000a0  20 43 25 3a de 20 6d 6f  6e 74 68 24 28 31 32 29  | C%:. month$(12)|
000000b0  3a e3 20 43 25 3d 31 20  b8 20 31 32 3a f3 20 6d  |:. C%=1 . 12:. m|
000000c0  6f 6e 74 68 24 28 43 25  29 3a ed 20 43 25 0d 00  |onth$(C%):. C%..|
000000d0  32 5b dc 20 4a 61 6e 75  61 72 79 2c 46 65 62 72  |2[. January,Febr|
000000e0  75 61 72 79 2c 4d 61 72  63 68 2c 41 70 72 69 6c  |uary,March,April|
000000f0  2c 4d 61 79 2c 4a 75 6e  65 2c 4a 75 6c 79 2c 41  |,May,June,July,A|
00000100  75 67 75 73 74 2c 53 65  70 74 65 6d 62 65 72 2c  |ugust,September,|
00000110  4f 63 74 6f 62 65 72 2c  4e 6f 76 65 6d 62 65 72  |October,November|
00000120  2c 44 65 63 65 6d 62 65  72 0d 00 3c 30 f5 3a f1  |,December..<0.:.|
00000130  8a 31 31 2c 38 29 22 57  68 69 63 68 20 79 65 61  |.11,8)"Which yea|
00000140  72 3f 20 22 3a e8 8a 32  33 2c 38 29 22 22 59 25  |r? ":..23,8)""Y%|
00000150  3a 59 24 3d c3 28 59 25  29 0d 00 46 2a eb 31 3a  |:Y$=.(Y%)..F*.1:|
00000160  f1 27 27 22 53 65 61 72  63 68 69 6e 67 2e 2e 2e  |.''"Searching...|
00000170  2e 2e 22 3b 3a 4d 25 3d  31 3a 44 3d 31 3a f2 64  |..";:M%=1:D=1:.d|
00000180  61 79 73 0d 00 50 1b f5  3a f2 46 49 4e 44 3a fd  |ays..P..:.FIND:.|
00000190  20 79 5a 3d 59 25 2b 31  3a 2a 46 58 32 31 0d 00  | yZ=Y%+1:*FX21..|
000001a0  5a 65 eb 37 3a f1 27 27  27 22 4e 6f 20 6d 6f 72  |Ze.7:.'''"No mor|
000001b0  65 20 6c 75 6e 61 72 20  65 63 6c 69 70 73 65 73  |e lunar eclipses|
000001c0  20 69 6e 20 22 3b 59 24  27 27 22 3c 51 3e 20 74  | in ";Y$''"<Q> t|
000001d0  6f 20 71 75 69 74 2c 20  61 6e 79 20 6f 74 68 65  |o quit, any othe|
000001e0  72 20 6b 65 79 20 66 6f  72 20 6d 6f 72 65 22 3b  |r key for more";|
000001f0  3a 67 25 3d a5 20 84 20  33 32 3a fd 20 67 25 3d  |:g%=. . 32:. g%=|
00000200  31 31 33 0d 00 64 07 db  3a e0 0d 00 6e 05 3a 0d  |113..d..:...n.:.|
00000210  00 78 0a dd f2 64 61 79  73 0d 00 82 0f 6d 25 3d  |.x...days....m%=|
00000220  4d 25 3a 79 25 3d 59 25  0d 00 8c 27 e7 ac 20 28  |M%:y%=Y%...'.. (|
00000230  6d 25 3d 31 20 84 20 6d  25 3d 32 29 20 8c 20 79  |m%=1 . m%=2) . y|
00000240  25 3d 79 25 2b 31 3a 6d  25 3d 6d 25 2d 31 32 0d  |%=y%+1:m%=m%-12.|
00000250  00 96 6d 79 25 3d 79 25  2d 31 3a 6d 25 3d 6d 25  |..my%=y%-1:m%=m%|
00000260  2b 31 32 3a 41 25 3d 79  25 20 81 20 31 30 30 3a  |+12:A%=y% . 100:|
00000270  42 25 3d 32 2d 41 25 2b  41 25 20 81 20 34 3a 43  |B%=2-A%+A% . 4:C|
00000280  25 3d a8 28 79 25 2a 33  36 35 2e 32 35 29 3a 44  |%=.(y%*365.25):D|
00000290  44 25 3d a8 28 33 30 2e  36 30 30 31 2a 28 6d 25  |D%=.(30.6001*(m%|
000002a0  2b 31 29 29 3a 44 41 59  53 3d 42 25 2b 43 25 2b  |+1)):DAYS=B%+C%+|
000002b0  44 44 25 2b 44 2d 37 32  33 32 34 34 0d 00 a0 05  |DD%+D-723244....|
000002c0  e1 0d 00 aa 05 3a 0d 00  b4 0a dd f2 46 49 4e 44  |.....:......FIND|
000002d0  0d 00 be 6d f2 53 55 4e  3a f2 4d 4f 4f 4e 3a 64  |...m.SUN:.MOON:d|
000002e0  69 66 66 65 72 65 6e 63  65 3d 4c 41 4d 42 44 41  |ifference=LAMBDA|
000002f0  53 55 4e 2d 4c 41 4d 42  44 41 4d 4f 4f 4e 3a 63  |SUN-LAMBDAMOON:c|
00000300  68 65 63 6b 3d 94 28 64  69 66 66 65 72 65 6e 63  |heck=.(differenc|
00000310  65 29 3a e7 20 63 68 65  63 6b 3c af 2f 32 20 8c  |e):. check<./2 .|
00000320  20 f2 64 61 74 65 3a 44  41 59 53 3d 44 41 59 53  | .date:DAYS=DAYS|
00000330  2b 31 35 3a e1 3a 8b 20  f2 46 55 4c 4c 0d 00 c8  |+15:.:. .FULL...|
00000340  10 44 41 59 53 3d 44 41  59 53 2b 33 30 0d 00 d2  |.DAYS=DAYS+30...|
00000350  05 e1 0d 00 dc 05 3a 0d  00 e6 09 dd f2 53 55 4e  |......:......SUN|
00000360  0d 00 f0 42 4e 3d 28 44  41 59 53 2f 33 36 35 2e  |...BN=(DAYS/365.|
00000370  32 34 32 32 29 2a 32 2a  af 3a 4e 3d a4 72 74 28  |2422)*2*.:N=.rt(|
00000380  4e 29 3a 4d 3d 4e 2d 2e  30 36 35 36 37 34 33 34  |N):M=N-.06567434|
00000390  39 3a e7 20 4d 3c 3d 30  20 8c 20 4d 3d 4d 2b 32  |9:. M<=0 . M=M+2|
000003a0  2a af 0d 00 fa 41 45 3d  a4 4b 65 70 6c 65 72 28  |*....AE=.Kepler(|
000003b0  4d 29 3a 4c 41 4d 42 44  41 53 55 4e 3d 45 2b 34  |M):LAMBDASUN=E+4|
000003c0  2e 39 33 32 32 33 37 36  38 37 3a 4c 41 4d 42 44  |.932237687:LAMBD|
000003d0  41 53 55 4e 3d a4 72 74  28 4c 41 4d 42 44 41 53  |ASUN=.rt(LAMBDAS|
000003e0  55 4e 29 0d 01 04 05 e1  0d 01 0e 05 3a 0d 01 18  |UN).........:...|
000003f0  0a dd f2 4d 4f 4f 4e 0d  01 22 de 6c 3d 31 2e 31  |...MOON..".l=1.1|
00000400  33 34 30 33 35 37 38 2b  44 41 59 53 2a 2e 32 32  |3403578+DAYS*.22|
00000410  39 39 37 31 35 30 34 32  3a 6c 3d a4 72 74 28 6c  |99715042:l=.rt(l|
00000420  29 3a 4d 6d 3d 6c 2d 2e  30 30 31 39 34 34 33 36  |):Mm=l-.00194436|
00000430  38 33 2a 44 41 59 53 2d  36 2e 30 39 37 38 38 34  |83*DAYS-6.097884|
00000440  38 3a 4d 6d 3d a4 72 74  28 4d 6d 29 3a 45 4e 3d  |8:Mm=.rt(Mm):EN=|
00000450  32 2e 36 35 32 30 33 35  32 38 36 2d 2e 30 30 30  |2.652035286-.000|
00000460  39 32 34 32 32 2a 44 41  59 53 3a 45 4e 3d a4 72  |92422*DAYS:EN=.r|
00000470  74 28 45 4e 29 3a 45 76  3d 2e 30 32 32 32 33 33  |t(EN):Ev=.022233|
00000480  37 34 39 33 2a b5 28 32  2a 28 6c 2d 4c 41 4d 42  |7493*.(2*(l-LAMB|
00000490  44 41 53 55 4e 29 2d 4d  6d 29 3a 41 65 3d 2e 30  |DASUN)-Mm):Ae=.0|
000004a0  30 33 32 34 32 38 32 31  38 2a b5 28 4d 29 3a 41  |032428218*.(M):A|
000004b0  33 3d 2e 30 30 36 34 35  37 37 31 38 32 2a b5 28  |3=.0064577182*.(|
000004c0  4d 29 3a 4d 74 69 63 6b  6d 3d 4d 6d 2b 45 76 2d  |M):Mtickm=Mm+Ev-|
000004d0  41 65 2d 41 33 0d 01 2c  da 45 63 3d 2e 31 30 39  |Ae-A3..,.Ec=.109|
000004e0  37 35 36 37 37 35 33 2a  b5 28 4d 74 69 63 6b 6d  |7567753*.(Mtickm|
000004f0  29 3a 41 34 3d 2e 30 30  33 37 33 35 30 30 34 36  |):A4=.0037350046|
00000500  2a b5 28 32 2a 4d 74 69  63 6b 6d 29 3a 6c 74 69  |*.(2*Mtickm):lti|
00000510  63 6b 3d 6c 2b 45 76 2b  45 63 2d 41 65 2b 41 34  |ck=l+Ev+Ec-Ae+A4|
00000520  3a 56 3d 2e 30 31 31 34  38 39 35 30 32 35 2a b5  |:V=.0114895025*.|
00000530  28 32 2a 28 6c 74 69 63  6b 2d 4c 41 4d 42 44 41  |(2*(ltick-LAMBDA|
00000540  53 55 4e 29 29 3a 6c 64  6f 75 62 6c 65 74 69 63  |SUN)):ldoubletic|
00000550  6b 3d 6c 74 69 63 6b 2b  56 3a 45 4e 74 69 63 6b  |k=ltick+V:ENtick|
00000560  3d 45 4e 2d 2e 30 30 32  37 39 32 35 32 36 38 2a  |=EN-.0027925268*|
00000570  b5 28 4d 29 3a 79 3d b5  28 6c 64 6f 75 62 6c 65  |.(M):y=.(ldouble|
00000580  74 69 63 6b 2d 45 4e 74  69 63 6b 29 2a 2e 39 39  |tick-ENtick)*.99|
00000590  35 39 37 30 33 32 31 3a  78 3d 9b 28 6c 64 6f 75  |5970321:x=.(ldou|
000005a0  62 6c 65 74 69 63 6b 2d  45 4e 74 69 63 6b 29 0d  |bletick-ENtick).|
000005b0  01 36 7d 74 61 6e 3d 99  28 79 2f 78 29 3a 74 61  |.6}tan=.(y/x):ta|
000005c0  6e 3d a4 72 74 71 75 61  64 28 74 61 6e 29 3a 4c  |n=.rtquad(tan):L|
000005d0  41 4d 42 44 41 4d 4f 4f  4e 3d 45 4e 74 69 63 6b  |AMBDAMOON=ENtick|
000005e0  2b 74 61 6e 3a 4c 41 4d  42 44 41 4d 4f 4f 4e 3d  |+tan:LAMBDAMOON=|
000005f0  a4 72 74 28 4c 41 4d 42  44 41 4d 4f 4f 4e 29 3a  |.rt(LAMBDAMOON):|
00000600  42 45 54 41 4d 4f 4f 4e  3d 98 28 b5 28 6c 64 6f  |BETAMOON=.(.(ldo|
00000610  75 62 6c 65 74 69 63 6b  2d 45 4e 74 69 63 6b 29  |ubletick-ENtick)|
00000620  2a 2e 30 38 39 36 38 33  34 34 32 29 0d 01 40 05  |*.089683442)..@.|
00000630  e1 0d 01 4a 05 3a 0d 01  54 0f dd a4 72 74 28 61  |...J.:..T...rt(a|
00000640  6e 67 6c 65 29 0d 01 5e  33 e7 20 61 6e 67 6c 65  |ngle)..^3. angle|
00000650  3c 30 20 8c 20 61 6e 67  6c 65 3d a8 28 94 28 61  |<0 . angle=.(.(a|
00000660  6e 67 6c 65 29 2f 32 2a  af 29 2a 32 2a af 2b 61  |ngle)/2*.)*2*.+a|
00000670  6e 67 6c 65 2b 32 2a af  0d 01 68 25 3d 32 2a af  |ngle+2*...h%=2*.|
00000680  2a 28 61 6e 67 6c 65 2f  28 32 2a af 29 2d a8 28  |*(angle/(2*.)-.(|
00000690  61 6e 67 6c 65 2f 28 32  2a af 29 29 29 0d 01 72  |angle/(2*.)))..r|
000006a0  05 3a 0d 01 7c 13 dd a4  4b 65 70 6c 65 72 28 41  |.:..|...Kepler(A|
000006b0  4e 47 4c 45 29 0d 01 86  73 6b 6c 3d 41 4e 47 4c  |NGLE)...skl=ANGL|
000006c0  45 3a f5 3a 64 65 6c 74  61 3d 6b 6c 2d 2e 30 31  |E:.:delta=kl-.01|
000006d0  36 37 31 38 2a b5 28 6b  6c 29 2d 41 4e 47 4c 45  |6718*.(kl)-ANGLE|
000006e0  3a 63 6f 72 72 65 63 74  69 6f 6e 3d 64 65 6c 74  |:correction=delt|
000006f0  61 2f 28 31 2d 2e 30 31  36 37 31 38 2a 9b 6b 6c  |a/(1-.016718*.kl|
00000700  29 3a 6b 6c 3d 6b 6c 2d  63 6f 72 72 65 63 74 69  |):kl=kl-correcti|
00000710  6f 6e 3a fd 20 94 28 64  65 6c 74 61 29 3c 3d 2e  |on:. .(delta)<=.|
00000720  30 30 30 30 30 30 30 31  0d 01 90 1d 3d 32 2a 99  |00000001....=2*.|
00000730  28 b7 28 6b 6c 2f 32 29  2a 31 2e 30 31 36 38 36  |(.(kl/2)*1.01686|
00000740  30 31 31 32 29 0d 01 9a  05 3a 0d 01 a4 11 dd a4  |0112)....:......|
00000750  72 74 71 75 61 64 28 74  61 6e 29 0d 01 ae 18 e7  |rtquad(tan).....|
00000760  20 79 3e 3d 30 20 80 20  78 3e 3d 30 20 8c 20 3d  | y>=0 . x>=0 . =|
00000770  74 61 6e 0d 01 b8 12 e7  20 78 3c 30 20 8c 20 3d  |tan..... x<0 . =|
00000780  74 61 6e 2b af 0d 01 c2  1a e7 20 79 3c 30 20 80  |tan+...... y<0 .|
00000790  20 78 3e 30 20 8c 20 3d  74 61 6e 2b 32 2a af 0d  | x>0 . =tan+2*..|
000007a0  01 cc 05 3a 0d 01 d6 0a  dd f2 64 61 74 65 0d 01  |...:......date..|
000007b0  e0 75 44 41 59 53 53 3d  44 41 59 53 2b 32 34 34  |.uDAYSS=DAYS+244|
000007c0  34 32 33 38 2e 35 3a 49  5a 3d a8 28 44 41 59 53  |4238.5:IZ=.(DAYS|
000007d0  53 2b 2e 35 29 3a 46 5a  3d 44 41 59 53 53 2b 2e  |S+.5):FZ=DAYSS+.|
000007e0  35 2d a8 28 44 41 59 53  53 2b 2e 35 29 3a e7 20  |5-.(DAYSS+.5):. |
000007f0  49 5a 3e 32 32 39 39 31  36 30 20 8c 20 41 5a 3d  |IZ>2299160 . AZ=|
00000800  a8 28 28 49 5a 2d 31 38  36 37 32 31 36 2e 32 35  |.((IZ-1867216.25|
00000810  29 2f 33 36 35 32 34 2e  32 35 29 3a 8b 20 41 5a  |)/36524.25):. AZ|
00000820  3d 49 5a 0d 01 ea b7 42  5a 3d 49 5a 2b 31 2b 41  |=IZ....BZ=IZ+1+A|
00000830  5a 2d a8 28 41 5a 2f 34  29 3a 43 5a 3d 42 5a 2b  |Z-.(AZ/4):CZ=BZ+|
00000840  31 35 32 34 3a 44 5a 3d  a8 28 28 43 5a 2d 31 32  |1524:DZ=.((CZ-12|
00000850  32 2e 31 29 2f 33 36 35  2e 32 35 29 3a 45 5a 3d  |2.1)/365.25):EZ=|
00000860  a8 28 33 36 35 2e 32 35  2a 44 5a 29 3a 47 5a 3d  |.(365.25*DZ):GZ=|
00000870  a8 28 28 43 5a 2d 45 5a  29 2f 33 30 2e 36 30 30  |.((CZ-EZ)/30.600|
00000880  31 29 3a 64 5a 3d 43 5a  2d 45 5a 2b 46 5a 2d a8  |1):dZ=CZ-EZ+FZ-.|
00000890  28 33 30 2e 36 30 30 31  2a 47 5a 29 3a 64 64 5a  |(30.6001*GZ):ddZ|
000008a0  3d a8 28 64 5a 29 3a 68  72 73 3d 44 41 59 53 2d  |=.(dZ):hrs=DAYS-|
000008b0  a8 28 44 41 59 53 29 3a  68 72 3d a8 28 32 34 2a  |.(DAYS):hr=.(24*|
000008c0  68 72 73 29 3a 6d 69 6e  3d a8 28 36 30 2a 28 32  |hrs):min=.(60*(2|
000008d0  34 2a 68 72 73 2d 68 72  29 29 0d 01 f4 22 e7 20  |4*hrs-hr))...". |
000008e0  47 5a 3c 31 33 2e 35 20  8c 20 6d 5a 3d 47 5a 2d  |GZ<13.5 . mZ=GZ-|
000008f0  31 3a 8b 20 6d 5a 3d 47  5a 2d 31 33 0d 01 fe 26  |1:. mZ=GZ-13...&|
00000900  e7 20 6d 5a 3e 32 2e 35  20 8c 20 79 5a 3d 44 5a  |. mZ>2.5 . yZ=DZ|
00000910  2d 34 37 31 36 3a 8b 20  79 5a 3d 44 5a 2d 34 37  |-4716:. yZ=DZ-47|
00000920  31 35 0d 02 08 29 6d 69  6e 24 3d c3 28 6d 69 6e  |15...)min$=.(min|
00000930  29 3a e7 20 a9 28 6d 69  6e 24 29 3c 32 20 6d 69  |):. .(min$)<2 mi|
00000940  6e 24 3d 22 30 22 2b 6d  69 6e 24 0d 02 12 24 68  |n$="0"+min$...$h|
00000950  72 24 3d c3 28 68 72 29  3a e7 20 a9 28 68 72 24  |r$=.(hr):. .(hr$|
00000960  29 3c 32 20 68 72 24 3d  22 30 22 2b 68 72 24 0d  |)<2 hr$="0"+hr$.|
00000970  02 1c 2e 64 61 74 65 24  3d c3 28 64 64 5a 29 3a  |...date$=.(ddZ):|
00000980  6d 6f 6e 74 68 24 3d 6d  6f 6e 74 68 24 28 6d 5a  |month$=month$(mZ|
00000990  29 3a 79 65 61 72 24 3d  c3 28 79 5a 29 0d 02 26  |):year$=.(yZ)..&|
000009a0  73 64 69 67 69 74 24 3d  c2 64 61 74 65 24 2c 31  |sdigit$=.date$,1|
000009b0  29 3a e7 20 64 69 67 69  74 24 3d 22 31 22 20 73  |):. digit$="1" s|
000009c0  75 66 66 69 78 24 3d 22  73 74 22 20 8b e7 20 64  |uffix$="st" .. d|
000009d0  69 67 69 74 24 3d 22 32  22 20 73 75 66 66 69 78  |igit$="2" suffix|
000009e0  24 3d 22 6e 64 22 20 8b  e7 20 64 69 67 69 74 24  |$="nd" .. digit$|
000009f0  3d 22 33 22 20 73 75 66  66 69 78 24 3d 22 72 64  |="3" suffix$="rd|
00000a00  22 20 8b 20 73 75 66 66  69 78 24 3d 22 74 68 22  |" . suffix$="th"|
00000a10  0d 02 30 20 e7 20 c0 64  61 74 65 24 2c 31 29 3d  |..0 . .date$,1)=|
00000a20  22 31 22 20 73 75 66 66  69 78 24 3d 22 74 68 22  |"1" suffix$="th"|
00000a30  0d 02 3a 05 e1 0d 02 44  05 3a 0d 02 4e 0a dd f2  |..:....D.:..N...|
00000a40  46 55 4c 4c 0d 02 58 b4  f5 3a 63 61 74 63 68 75  |FULL..X..:catchu|
00000a50  70 3d 38 2e 38 38 33 37  32 35 39 45 2d 33 2b 31  |p=8.8837259E-3+1|
00000a60  2e 30 34 37 31 39 37 36  45 2d 33 2a 9b 28 4d 74  |.0471976E-3*.(Mt|
00000a70  69 63 6b 6d 29 3a e7 20  4c 41 4d 42 44 41 53 55  |ickm):. LAMBDASU|
00000a80  4e 3c af 20 8c 20 64 69  66 66 65 72 65 6e 63 65  |N<. . difference|
00000a90  3d 64 69 66 66 65 72 65  6e 63 65 2b af 20 8b 20  |=difference+. . |
00000aa0  64 69 66 66 65 72 65 6e  63 65 3d 64 69 66 66 65  |difference=diffe|
00000ab0  72 65 6e 63 65 2d af 3a  e7 20 94 28 64 69 66 66  |rence-.:. .(diff|
00000ac0  65 72 65 6e 63 65 29 3e  af 20 80 20 b4 28 64 69  |erence)>. . .(di|
00000ad0  66 66 65 72 65 6e 63 65  29 3d 2d 31 20 8c 20 64  |fference)=-1 . d|
00000ae0  69 66 66 65 72 65 6e 63  65 3d 32 2a af 2b 64 69  |ifference=2*.+di|
00000af0  66 66 65 72 65 6e 63 65  0d 02 62 52 44 41 59 53  |fference..bRDAYS|
00000b00  3d 44 41 59 53 2b 28 64  69 66 66 65 72 65 6e 63  |=DAYS+(differenc|
00000b10  65 2f 28 32 34 2a 63 61  74 63 68 75 70 29 29 3a  |e/(24*catchup)):|
00000b20  f2 53 55 4e 3a f2 4d 4f  4f 4e 3a 64 69 66 66 65  |.SUN:.MOON:diffe|
00000b30  72 65 6e 63 65 3d 4c 41  4d 42 44 41 53 55 4e 2d  |rence=LAMBDASUN-|
00000b40  4c 41 4d 42 44 41 4d 4f  4f 4e 0d 02 6c 28 fd 20  |LAMBDAMOON..l(. |
00000b50  94 28 94 28 64 69 66 66  65 72 65 6e 63 65 29 2d  |.(.(difference)-|
00000b60  af 29 3c 20 30 2e 30 30  30 30 30 31 3a f2 64 61  |.)< 0.000001:.da|
00000b70  74 65 0d 02 76 57 74 65  73 74 3d 94 28 94 28 6c  |te..vWtest=.(.(l|
00000b80  64 6f 75 62 6c 65 74 69  63 6b 2d 45 4e 74 69 63  |doubletick-ENtic|
00000b90  6b 29 2d af 29 3a e7 20  74 65 73 74 3c 2e 32 31  |k)-.):. test<.21|
00000ba0  33 38 30 32 38 33 33 34  20 84 20 74 65 73 74 3e  |38028334 . test>|
00000bb0  32 2e 39 32 37 37 38 39  38 32 20 8c 20 ef 37 3a  |2.92778982 . .7:|
00000bc0  f2 5f 64 69 61 67 72 61  6d 0d 02 80 05 e1 0d 02  |._diagram.......|
00000bd0  8a 05 3a 0d 02 94 0e dd  f2 5f 64 69 61 67 72 61  |..:......_diagra|
00000be0  6d 0d 02 9e 42 72 68 6f  74 69 63 6b 3d 2e 39 39  |m...Brhotick=.99|
00000bf0  36 39 38 35 39 39 2f 28  31 2b 2e 30 35 34 39 2a  |698599/(1+.0549*|
00000c00  9b 28 4d 74 69 63 6b 6d  2b 45 63 29 29 3a 70 69  |.(Mtickm+Ec)):pi|
00000c10  3d 2e 30 31 36 35 39 32  38 34 35 2f 72 68 6f 74  |=.016592845/rhot|
00000c20  69 63 6b 0d 02 a8 5d 64  65 6c 74 61 6c 61 6d 62  |ick...]deltalamb|
00000c30  64 61 3d 2e 30 30 39 35  39 39 33 31 31 2b 2e 30  |da=.009599311+.0|
00000c40  30 31 30 34 37 31 39 37  35 35 2a 9b 28 4d 74 69  |0104719755*.(Mti|
00000c50  63 6b 6d 29 3a 64 65 6c  74 61 62 65 74 61 3d 38  |ckm):deltabeta=8|
00000c60  2e 37 32 36 36 34 36 45  2d 34 2a 9b 28 6c 64 6f  |.726646E-4*.(ldo|
00000c70  75 62 6c 65 74 69 63 6b  2d 45 4e 74 69 63 6b 29  |ubletick-ENtick)|
00000c80  0d 02 b2 5c 53 75 3d 70  69 2d 2e 30 30 34 37 31  |...\Su=pi-.00471|
00000c90  32 33 38 38 39 38 3a 74  68 65 74 61 6d 3d 2e 30  |238898:thetam=.0|
00000ca0  30 39 30 34 32 35 35 31  2f 72 68 6f 74 69 63 6b  |09042551/rhotick|
00000cb0  3a 72 61 74 65 3d 64 65  6c 74 61 6c 61 6d 62 64  |:rate=deltalambd|
00000cc0  61 2d 37 2e 31 36 37 38  33 45 2d 34 3a 66 61 63  |a-7.16783E-4:fac|
00000cd0  74 6f 72 3d 32 35 36 2f  72 61 74 65 0d 02 bc 57  |tor=256/rate...W|
00000ce0  61 3d 31 2b 64 65 6c 74  61 62 65 74 61 5e 32 2f  |a=1+deltabeta^2/|
00000cf0  72 61 74 65 5e 32 3a 62  3d 32 2a 42 45 54 41 4d  |rate^2:b=2*BETAM|
00000d00  4f 4f 4e 2a 64 65 6c 74  61 62 65 74 61 2f 72 61  |OON*deltabeta/ra|
00000d10  74 65 3a 48 3d 28 53 75  2b 74 68 65 74 61 6d 2f  |te:H=(Su+thetam/|
00000d20  32 29 5e 32 3a 63 3d 42  45 54 41 4d 4f 4f 4e 5e  |2)^2:c=BETAMOON^|
00000d30  32 2d 48 0d 02 c6 27 72  6f 6f 74 3d 62 5e 32 2d  |2-H...'root=b^2-|
00000d40  34 2a 61 2a 63 3a e7 20  b4 28 72 6f 6f 74 29 3d  |4*a*c:. .(root)=|
00000d50  2d 31 20 f2 5f 70 65 6e  3a e1 0d 02 d0 4a 58 32  |-1 ._pen:....JX2|
00000d60  3d 28 2d 62 2b b6 28 72  6f 6f 74 29 29 2f 28 32  |=(-b+.(root))/(2|
00000d70  2a 61 29 3a 58 31 3d 28  2d 62 2d b6 28 72 6f 6f  |*a):X1=(-b-.(roo|
00000d80  74 29 29 2f 28 32 2a 61  29 3a e7 20 b4 28 58 32  |t))/(2*a):. .(X2|
00000d90  29 3d 2d 31 20 58 3d 58  31 3a 58 31 3d 58 32 3a  |)=-1 X=X1:X1=X2:|
00000da0  58 32 3d 58 0d 02 da 87  73 74 61 72 74 68 72 73  |X2=X....starthrs|
00000db0  3d 58 31 2f 72 61 74 65  3a 65 6e 64 68 72 73 3d  |=X1/rate:endhrs=|
00000dc0  58 32 2f 72 61 74 65 3a  6d 69 64 3d 32 34 2a 68  |X2/rate:mid=24*h|
00000dd0  72 73 3a 73 74 61 72 74  74 69 6d 65 24 3d a4 5f  |rs:starttime$=._|
00000de0  44 44 74 6f 44 4d 28 6d  69 64 2b 73 74 61 72 74  |DDtoDM(mid+start|
00000df0  68 72 73 29 3a 65 6e 64  74 69 6d 65 24 3d a4 5f  |hrs):endtime$=._|
00000e00  44 44 74 6f 44 4d 28 6d  69 64 2b 65 6e 64 68 72  |DDtoDM(mid+endhr|
00000e10  73 29 3a 6d 69 64 24 3d  68 72 24 2b 22 3a 22 2b  |s):mid$=hr$+":"+|
00000e20  6d 69 6e 24 2b 22 20 47  4d 54 22 0d 02 e4 0b db  |min$+" GMT".....|
00000e30  3a da 3a ef 32 30 0d 02  ee 51 f1 22 45 63 6c 69  |:.:.20...Q."Ecli|
00000e40  70 73 65 20 70 61 72 74  69 63 75 6c 61 72 73 20  |pse particulars |
00000e50  61 72 65 3a 22 27 27 8a  35 29 22 44 61 74 65 3a  |are:"''.5)"Date:|
00000e60  20 22 3b 64 61 74 65 24  3b 73 75 66 66 69 78 24  | ";date$;suffix$|
00000e70  3b 22 20 22 3b 6d 6f 6e  74 68 24 3b 22 20 22 3b  |;" ";month$;" ";|
00000e80  79 65 61 72 24 27 27 0d  02 f8 60 f1 8a 35 29 22  |year$''...`..5)"|
00000e90  45 63 6c 69 70 73 65 20  62 65 67 69 6e 73 3a 20  |Eclipse begins: |
00000ea0  22 3b 73 74 61 72 74 74  69 6d 65 24 27 8a 35 29  |";starttime$'.5)|
00000eb0  22 4d 69 64 2d 65 63 6c  69 70 73 65 3a 20 20 20  |"Mid-eclipse:   |
00000ec0  20 22 3b 6d 69 64 24 3b  27 8a 35 29 22 45 63 6c  | ";mid$;'.5)"Ecl|
00000ed0  69 70 73 65 20 65 6e 64  73 3a 20 20 20 22 3b 65  |ipse ends:   ";e|
00000ee0  6e 64 74 69 6d 65 24 0d  03 02 09 2a 46 58 32 31  |ndtime$....*FX21|
00000ef0  0d 03 0c 3f f1 27 27 27  22 41 6e 79 20 6b 65 79  |...?.'''"Any key|
00000f00  20 66 6f 72 20 74 68 65  20 64 69 61 67 72 61 6d  | for the diagram|
00000f10  22 3b 3a 47 25 3d a5 3a  ef 31 39 2c 32 2c 34 2c  |";:G%=.:.19,2,4,|
00000f20  30 2c 30 2c 30 3a fb 31  33 30 3a e6 30 2c 30 0d  |0,0,0:.130:.0,0.|
00000f30  03 16 58 db 3a ec 30 2c  35 31 32 3a df 31 32 38  |..X.:.0,512:.128|
00000f40  30 2c 35 31 32 3a ec 30  2c 35 35 30 3a ef 35 3a  |0,512:.0,550:.5:|
00000f50  f1 22 45 63 6c 69 70 74  69 63 22 3a ec 31 30 31  |."Ecliptic":.101|
00000f60  38 2c 35 35 30 3a f1 64  61 74 65 24 3b 22 3a 22  |8,550:.date$;":"|
00000f70  3b c3 28 6d 5a 29 3b 22  3a 22 3b c2 79 65 61 72  |;.(mZ);":";.year|
00000f80  24 2c 32 29 3a ef 34 0d  03 20 23 72 3d 53 75 2a  |$,2):.4.. #r=Su*|
00000f90  66 61 63 74 6f 72 3a f2  5f 63 69 72 63 6c 65 28  |factor:._circle(|
00000fa0  36 34 30 2c 35 31 32 2c  72 29 0d 03 2a b6 6c 68  |640,512,r)..*.lh|
00000fb0  3d 28 32 2e 35 2a 64 65  6c 74 61 62 65 74 61 2b  |=(2.5*deltabeta+|
00000fc0  42 45 54 41 4d 4f 4f 4e  29 2a 66 61 63 74 6f 72  |BETAMOON)*factor|
00000fd0  2b 35 31 32 3a ec 30 2c  6c 68 3a df 31 32 38 30  |+512:.0,lh:.1280|
00000fe0  2c 28 42 45 54 41 4d 4f  4f 4e 2d 32 2e 35 2a 64  |,(BETAMOON-2.5*d|
00000ff0  65 6c 74 61 62 65 74 61  29 2a 66 61 63 74 6f 72  |eltabeta)*factor|
00001000  2b 35 31 32 3a ec 30 2c  6c 68 3a df 35 30 2c 28  |+512:.0,lh:.50,(|
00001010  32 2e 33 2a 64 65 6c 74  61 62 65 74 61 2b 42 45  |2.3*deltabeta+BE|
00001020  54 41 4d 4f 4f 4e 29 2a  66 61 63 74 6f 72 2b 35  |TAMOON)*factor+5|
00001030  33 30 3a ec 30 2c 6c 68  3a df 35 30 2c 28 42 45  |30:.0,lh:.50,(BE|
00001040  54 41 4d 4f 4f 4e 2b 32  2e 33 2a 64 65 6c 74 61  |TAMOON+2.3*delta|
00001050  62 65 74 61 29 2a 66 61  63 74 6f 72 2b 35 30 30  |beta)*factor+500|
00001060  0d 03 34 08 e6 30 2c 33  0d 03 3e 7e 72 3d 74 68  |..4..0,3..>~r=th|
00001070  65 74 61 6d 2f 32 2a 66  61 63 74 6f 72 3a f2 5f  |etam/2*factor:._|
00001080  63 69 72 63 6c 65 28 36  34 30 2d 58 31 2a 66 61  |circle(640-X1*fa|
00001090  63 74 6f 72 2c 28 42 45  54 41 4d 4f 4f 4e 2b 73  |ctor,(BETAMOON+s|
000010a0  74 61 72 74 68 72 73 2a  64 65 6c 74 61 62 65 74  |tarthrs*deltabet|
000010b0  61 29 2a 66 61 63 74 6f  72 2b 35 31 32 2c 72 29  |a)*factor+512,r)|
000010c0  3a f1 8a 33 33 2c 31 35  2b 31 32 2a b4 28 42 45  |:..33,15+12*.(BE|
000010d0  54 41 4d 4f 4f 4e 29 29  c0 73 74 61 72 74 74 69  |TAMOON)).startti|
000010e0  6d 65 24 2c 35 29 0d 03  48 5b 79 3d 42 45 54 41  |me$,5)..H[y=BETA|
000010f0  4d 4f 4f 4e 2a 66 61 63  74 6f 72 2b 35 31 32 3a  |MOON*factor+512:|
00001100  e6 30 2c 31 3a fb 31 3a  f2 5f 63 69 72 63 6c 65  |.0,1:.1:._circle|
00001110  28 36 34 30 2c 79 2c 72  29 3a f1 8a 31 37 2c 31  |(640,y,r):..17,1|
00001120  35 2b 31 32 2a b4 28 42  45 54 41 4d 4f 4f 4e 29  |5+12*.(BETAMOON)|
00001130  29 c0 6d 69 64 24 2c 35  29 3a e6 30 2c 33 3a fb  |).mid$,5):.0,3:.|
00001140  33 0d 03 52 6d f2 5f 63  69 72 63 6c 65 28 36 34  |3..Rm._circle(64|
00001150  30 2d 58 32 2a 66 61 63  74 6f 72 2c 28 42 45 54  |0-X2*factor,(BET|
00001160  41 4d 4f 4f 4e 2b 65 6e  64 68 72 73 2a 64 65 6c  |AMOON+endhrs*del|
00001170  74 61 62 65 74 61 29 2a  66 61 63 74 6f 72 2b 35  |tabeta)*factor+5|
00001180  31 32 2c 72 29 3a f1 8a  33 2c 31 35 2b 31 32 2a  |12,r):..3,15+12*|
00001190  b4 28 42 45 54 41 4d 4f  4f 4e 29 29 c0 65 6e 64  |.(BETAMOON)).end|
000011a0  74 69 6d 65 24 2c 35 29  3a 2a 46 58 32 31 0d 03  |time$,5):*FX21..|
000011b0  5c 4c f1 8a 31 30 2c 33  31 29 22 41 6e 79 20 6b  |\L..10,31)"Any k|
000011c0  65 79 20 74 6f 20 63 6f  6e 74 69 6e 75 65 2e 2e  |ey to continue..|
000011d0  2e 22 3b 3a 67 25 3d a5  3a f1 8a 31 30 2c 33 31  |.";:g%=.:..10,31|
000011e0  29 22 53 65 61 72 63 68  20 72 65 73 75 6d 65 64  |)"Search resumed|
000011f0  2e 2e 2e 2e 2e 2e 2e 2e  22 3b 0d 03 66 05 e1 0d  |........";..f...|
00001200  03 70 05 3a 0d 03 7a 0a  dd f2 5f 70 65 6e 0d 03  |.p.:..z..._pen..|
00001210  84 2f 6d 69 64 3d 32 34  2a 68 72 73 3a 6d 69 64  |./mid=24*hrs:mid|
00001220  24 3d 68 72 24 2b 22 3a  22 2b 6d 69 6e 24 2b 22  |$=hr$+":"+min$+"|
00001230  20 47 4d 54 22 3a db 3a  da 3a ef 32 30 0d 03 8e  | GMT":.:.:.20...|
00001240  7a 20 f1 27 27 27 22 54  68 69 73 20 69 73 20 61  |z .'''"This is a|
00001250  20 70 65 6e 75 6d 62 72  61 6c 20 65 63 6c 69 70  | penumbral eclip|
00001260  73 65 3a 22 27 27 8a 35  29 22 44 61 74 65 3a 20  |se:"''.5)"Date: |
00001270  22 3b 64 61 74 65 24 3b  73 75 66 66 69 78 24 3b  |";date$;suffix$;|
00001280  22 20 22 3b 6d 6f 6e 74  68 24 3b 22 20 22 3b 79  |" ";month$;" ";y|
00001290  65 61 72 24 27 27 8a 35  29 22 4d 69 64 2d 65 63  |ear$''.5)"Mid-ec|
000012a0  6c 69 70 73 65 3a 20 20  20 20 22 3b 6d 69 64 24  |lipse:    ";mid$|
000012b0  3b 3a 2a 46 58 32 31 0d  03 98 51 f1 27 27 27 22  |;:*FX21...Q.'''"|
000012c0  41 6e 79 20 6b 65 79 20  66 6f 72 20 74 68 65 20  |Any key for the |
000012d0  64 69 61 67 72 61 6d 22  3b 3a 47 25 3d a5 3a ef  |diagram";:G%=.:.|
000012e0  31 39 2c 31 32 39 2c 31  33 32 2c 30 2c 30 2c 30  |19,129,132,0,0,0|
000012f0  2c 31 39 2c 32 2c 30 2c  30 2c 30 2c 30 3a fb 20  |,19,2,0,0,0,0:. |
00001300  31 32 39 3a e6 30 2c 30  0d 03 a2 78 db 3a ec 30  |129:.0,0...x.:.0|
00001310  2c 35 31 32 3a df 31 32  38 30 2c 35 31 32 3a ec  |,512:.1280,512:.|
00001320  30 2c 35 35 30 3a ef 35  3a f1 22 45 63 6c 69 70  |0,550:.5:."Eclip|
00001330  74 69 63 22 3a ec 31 30  31 38 2c 35 35 30 3a f1  |tic":.1018,550:.|
00001340  64 61 74 65 24 3b 22 3a  22 3b c3 28 6d 5a 29 3b  |date$;":";.(mZ);|
00001350  22 3a 22 3b c2 79 65 61  72 24 2c 32 29 3a ef 34  |":";.year$,2):.4|
00001360  3a 72 3d 53 75 2a 66 61  63 74 6f 72 3a f2 5f 63  |:r=Su*factor:._c|
00001370  69 72 63 6c 65 28 36 34  30 2c 35 31 32 2c 72 29  |ircle(640,512,r)|
00001380  0d 03 ac bb 6c 68 3d 28  32 2e 35 2a 64 65 6c 74  |....lh=(2.5*delt|
00001390  61 62 65 74 61 2b 42 45  54 41 4d 4f 4f 4e 29 2a  |abeta+BETAMOON)*|
000013a0  66 61 63 74 6f 72 2b 35  31 32 3a ec 30 2c 6c 68  |factor+512:.0,lh|
000013b0  3a df 31 32 38 30 2c 28  42 45 54 41 4d 4f 4f 4e  |:.1280,(BETAMOON|
000013c0  2d 32 2e 35 2a 64 65 6c  74 61 62 65 74 61 29 2a  |-2.5*deltabeta)*|
000013d0  66 61 63 74 6f 72 2b 35  31 32 3a ec 30 2c 6c 68  |factor+512:.0,lh|
000013e0  3a df 35 30 2c 28 32 2e  33 2a 64 65 6c 74 61 62  |:.50,(2.3*deltab|
000013f0  65 74 61 2b 42 45 54 41  4d 4f 4f 4e 29 2a 66 61  |eta+BETAMOON)*fa|
00001400  63 74 6f 72 2b 35 33 30  3a ec 30 2c 6c 68 3a df  |ctor+530:.0,lh:.|
00001410  35 30 2c 28 42 45 54 41  4d 4f 4f 4e 2b 32 2e 33  |50,(BETAMOON+2.3|
00001420  2a 64 65 6c 74 61 62 65  74 61 29 2a 66 61 63 74  |*deltabeta)*fact|
00001430  6f 72 2b 35 30 30 3a e6  30 2c 33 0d 03 b6 63 79  |or+500:.0,3...cy|
00001440  3d 42 45 54 41 4d 4f 4f  4e 2a 66 61 63 74 6f 72  |=BETAMOON*factor|
00001450  2b 35 31 32 3a 72 3d 74  68 65 74 61 6d 2f 32 2a  |+512:r=thetam/2*|
00001460  66 61 63 74 6f 72 3a f2  5f 63 69 72 63 6c 65 28  |factor:._circle(|
00001470  36 34 30 2c 79 2c 72 29  3a f1 8a 31 37 2c 31 35  |640,y,r):..17,15|
00001480  2b 31 32 2a b4 28 42 45  54 41 4d 4f 4f 4e 29 29  |+12*.(BETAMOON))|
00001490  c0 6d 69 64 24 2c 35 29  3a 2a 46 58 32 31 0d 03  |.mid$,5):*FX21..|
000014a0  c0 4c f1 8a 31 30 2c 33  31 29 22 41 6e 79 20 6b  |.L..10,31)"Any k|
000014b0  65 79 20 74 6f 20 63 6f  6e 74 69 6e 75 65 2e 2e  |ey to continue..|
000014c0  2e 22 3b 3a 67 25 3d a5  3a f1 8a 31 30 2c 33 31  |.";:g%=.:..10,31|
000014d0  29 22 53 65 61 72 63 68  20 72 65 73 75 6d 65 64  |)"Search resumed|
000014e0  2e 2e 2e 2e 2e 2e 2e 2e  22 3b 0d 03 ca 05 e1 0d  |........";......|
000014f0  03 d4 05 3a 0d 03 de 13  dd a4 5f 44 44 74 6f 44  |...:......_DDtoD|
00001500  4d 28 74 69 6d 65 29 0d  03 e8 19 ea 20 68 72 73  |M(time)..... hrs|
00001510  2c 68 72 73 24 2c 6d 69  6e 73 2c 6d 69 6e 73 24  |,hrs$,mins,mins$|
00001520  0d 03 f2 35 68 72 73 3d  a8 28 74 69 6d 65 29 3a  |...5hrs=.(time):|
00001530  68 72 73 24 3d c3 28 68  72 73 29 3a e7 20 a9 28  |hrs$=.(hrs):. .(|
00001540  68 72 73 24 29 3c 32 20  68 72 73 24 3d 22 30 22  |hrs$)<2 hrs$="0"|
00001550  2b 68 72 73 24 0d 03 fc  47 6d 69 6e 73 3d a8 28  |+hrs$...Gmins=.(|
00001560  28 74 69 6d 65 2d 68 72  73 29 2a 36 30 2b 2e 35  |(time-hrs)*60+.5|
00001570  29 3a 6d 69 6e 73 24 3d  c3 28 6d 69 6e 73 29 3a  |):mins$=.(mins):|
00001580  e7 20 a9 28 6d 69 6e 73  24 29 3c 32 20 6d 69 6e  |. .(mins$)<2 min|
00001590  73 24 3d 22 30 22 2b 6d  69 6e 73 24 0d 04 06 1a  |s$="0"+mins$....|
000015a0  3d 68 72 73 24 2b 22 3a  22 2b 6d 69 6e 73 24 2b  |=hrs$+":"+mins$+|
000015b0  22 20 47 4d 54 22 0d 04  10 05 3a 0d 04 1a 21 dd  |" GMT"....:...!.|
000015c0  f2 5f 63 69 72 63 6c 65  28 6f 78 2c 6f 79 2c 72  |._circle(ox,oy,r|
000015d0  61 64 69 75 73 29 3a ea  20 78 2c 79 0d 04 24 34  |adius):. x,y..$4|
000015e0  ec 20 6f 78 2d 72 61 64  69 75 73 2c 6f 79 3a df  |. ox-radius,oy:.|
000015f0  20 6f 78 2b 72 61 64 69  75 73 2c 6f 79 3a e3 20  | ox+radius,oy:. |
00001600  79 3d 31 20 b8 20 72 61  64 69 75 73 20 88 20 34  |y=1 . radius . 4|
00001610  0d 04 2e 1e 78 3d b6 28  28 72 61 64 69 75 73 2b  |....x=.((radius+|
00001620  79 29 2a 28 72 61 64 69  75 73 2d 79 29 29 0d 04  |y)*(radius-y))..|
00001630  38 2b 6f 78 70 78 3d 6f  78 2b 78 3a 6f 78 6d 78  |8+oxpx=ox+x:oxmx|
00001640  3d 6f 78 2d 78 3a 6f 79  70 79 3d 6f 79 2b 79 3a  |=ox-x:oypy=oy+y:|
00001650  6f 79 6d 79 3d 6f 79 2d  79 0d 04 42 1b ec 20 6f  |oymy=oy-y..B.. o|
00001660  78 6d 78 2c 6f 79 70 79  3a df 20 6f 78 70 78 2c  |xmx,oypy:. oxpx,|
00001670  6f 79 70 79 0d 04 4c 1b  ec 20 6f 78 6d 78 2c 6f  |oypy..L.. oxmx,o|
00001680  79 6d 79 3a df 20 6f 78  70 78 2c 6f 79 6d 79 0d  |ymy:. oxpx,oymy.|
00001690  04 56 05 ed 0d 04 60 05  e1 0d ff                 |.V....`....|
0000169b

© Dominic Ford 2021–2025.

Our privacy policy is here.
Last updated: 01 May 2025, 21:49 UTC
Website designed by .