Home » CEEFAX disks » telesoftware12.adl » 04-02-89/Alginon

04-02-89/Alginon

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 » CEEFAX disks » telesoftware12.adl
Filename: 04-02-89/Alginon
Read OK: ✔
File size: 1554 bytes
Load address: 1B00
Exec address: FFFF8023
Duplicates

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

File contents
>L.

10000DEF PROCfile(open,name$)

10010IF NOT open THEN PTR#E%=0:PRINT#E%,total:CLOSE#0:ENDPROC

10020CLOSE#0

10030E%=OPENUP(name$)

10040IF E%=0 CLOSE#0:E%=OPENOUT(name$):CLOSE#0:E%=OPENUP(name$):total=1:PRINT#E%,total,E%,E%:ENDPROC

10050IF E%<>0 INPUT#E%,total

10060ENDPROC

10070DEF PROCdata(pair,x,y,write,alterX$,alterY$)

10080LOCAL Z%,xn

10090IF pair<1 THEN PROCfile(FALSE,""):PRINT "ERROR PAIR < 1":END

10100IF pair>total THEN PTR#E%=12*total-6:INPUT#E%,Z%,Z%:FOR Z%=total+1 TO pair:PRINT#E%,Z%,Z%:NEXT

10110PTR#E%=12*pair-6

10120IF write THEN PRINT#E%,x,y ELSE INPUT#E%,X,Y

10130IF pair>total THEN total=pair

10140IF NOT write xn=EVAL(alterX$):Y=EVAL(alterY$):X=xn

10150ENDPROC

10160DEF PROCgstart

10170CLG

10180PROCpos(TRUE,0,0,0,0,0,0)

10190PROCsca(TRUE,0,0,0,0,0,0)

10200CLOSE#0

10210ON ERROR PROCerror:GOTO L%

10220L%=0

10230ENDPROC

10240DEF PROCgstop

10250CLOSE#0

10260ENDPROC

10270DEF PROCpos(default,Xtmi,Xtma,xs,Ytmi,Ytma,ys)

10280IF default THEN Xtmin=15:Xtmax=1279:xsl=15:Ytmin=15:Ytmax=1023:ysl=15

10290IF NOT default THEN Xtmin=Xtmi:Xtmax=Xtma:xsl=xs:Ytmin=Ytmi:Ytmax=Ytma:ysl=ys

10300Pry=Ytmax-Ytmin

10310Prx=Xtmax-Xtmin

10320ENDPROC

10330DEF PROCsca(default,Xsmi,Xsma,sx,Ysmi,Ysma,sy)

10340IF default THEN Xsmin=-1:Xsmax=1:stx=0.1:Ysmin=-1:Ysmax=1:sty=0.1

10350IF NOT default THEN Xsmin=Xsmi:Xsmax=Xsma:stx=sx:Ysmin=Ysmi:Ysmax=Ysma:sty=sy

10360PROCscalevars

10370ENDPROC

10380DEF PROCAutoScales(name$,alterX$,alterY$)

10390LOCAL RangeX,RangeY,Pwx%,Pwy%,X%

10400PROCfile(TRUE,name$)

10410PROCdata(1,0,0,FALSE,alterX$,alterY$)

10420Xsmax=X:Xsmin=X

10430Ysmin=Y:Ysmax=Y

10440FOR X%=2 TO total

10450PROCdata(X%,0,0,FALSE,alterX$,alterY$)

10460IF X>Xsmax THEN Xsmax=X

10470IF X<Xsmin THEN Xsmin=X

10480IF Y>Ysmax THEN Ysmax=Y

10490IF Y<Ysmin THEN Ysmin=Y

10500NEXT X%

10510RangeX=Xsmax-Xsmin

10520RangeY=Ysmax-Ysmin

10530Pwx%=FNpower(RangeX)

10540Pwy%=FNpower(RangeY)

10550Xsmax=(INT(Xsmax/10^Pwx%+1))*10^Pwx%

10560Xsmin=(INT(Xsmin/10^Pwx%))*10^Pwx%

10570Ysmax=(INT(Ysmax/10^Pwy%+1))*10^Pwy%

10580Ysmin=(INT(Ysmin/10^Pwy%))*10^Pwy%

10590stx=FNstepsize(Xsmax-Xsmin)

10600sty=FNstepsize(Ysmax-Ysmin)

10610PROCscalevars

10620ENDPROC

10630DEF FNpower(number)

10640=INT(LOG(number))

10650DEF FNstepsize(range)

10660LOCAL power% 

10670power%=FNpower(range)

10680range=range/10^power%

10690IF range<=3 THEN =2*10^(power%-1)

10700IF range>3 AND range<7 THEN =5*10^(power%-1)

10710=1*10^power%

10720DEF PROCscalevars

10730Sry=Ysmax-Ysmin

10740Y0=Ytmax-Pry*Ysmax/Sry

10750Srx=Xsmax-Xsmin

10760X0=Xtmax-Prx*Xsmax/Srx

10770Xsmin2=Xsmin

10780ENDPROC

10790DEF PROCDrawfile(name$,pen,alterX$,alterY$,crosses)

10800LOCAL PEN,X%,P

10810PROCfile(TRUE,name$)

10820FOR X%=1 TO total

10830PROCdata(X%,0,0,FALSE,alterX$,alterY$)

10840Y=FNy(Y):X=FNx(X)

10850IF X%=1 THEN PEN=0 ELSE PEN=pen 

10860IF NOT crosses THEN PROCdraw(X,Y,PEN):GOTO10950

10870P=0

10880PROCdraw(X-10,Y,P)

10890P=pen

10900PROCdraw(X+10,Y,P)

10910P=0

10920PROCdraw(X,Y-10,P)

10930P=pen

10940PROCdraw(X,Y+10,P)

10950NEXT X%

10960PROCfile(FALSE,name$)

10970ENDPROC

10980DEF PROCdraw(x,y,pen)

10990IF pen=0 MOVE x,y ELSE GCOL0,pen:DRAWx,y 

11000ENDPROC

11010DEF FNx(X)

11020=(Prx/Srx)*X+X0

11030DEF FNy(Y)

11040=(Pry/Sry)*Y+Y0

11050DEF PROCscale(pen)

11060LOCAL X,X2,Y,Y2,P

11070IF Y0>Ytmax THEN Y2=Ytmax:GOTO11120

11080IF Y0<Ytmin THEN Y2=Ytmin:GOTO11120

11090Y2=Y0

11100X=Xtmin:Y=Y2:P=0:PROCdraw(X,Y,P)

11110X=Xtmax:P=pen:PROCdraw(X,Y,P)

11120IF X0>Xtmax THEN X2=Xtmax:GOTO11170

11130IF X0<Xtmin THEN X2=Xtmin:GOTO11170

11140X2=X0

11150X=X2:Y=Ytmin:P=0:PROCdraw(X,Y,P)

11160Y=Ytmax:P=pen:PROCdraw(X,Y,P)

11170FOR X=Xsmax TO Xsmin STEP -1*stx

11180x1=FNx(X):Y=Y2:P=0

11190PROCdraw(x1,Y,P):IF X=0 THEN P=0 ELSE P=pen

11200Y=Y2-xsl:PROCdraw(x1,Y,P)

11210NEXT X

11220FOR Y=Ysmin TO Ysmax STEP sty

11230y1=FNy(Y):X=X2:P=0

11240PROCdraw(X,y1,P):IF Y=0 THEN P=0 ELSE P=pen

11250X=X2-ysl:PROCdraw(X,y1,P)

11260NEXT Y

11270ENDPROC

11280DEF PROCplot(pen,equ$)

11290LOCAL X,y1,y2,P

11300y2=2

11310FOR X=Xsmin2 TO Xsmax STEP Srx/Prx*10

11320Y=FNy(EVAL(equ$)) 

11330IF Y>Ytmax OR Y<Ytmin THEN y1=FALSE:y2=2 ELSE y1=TRUE:y2=y2-1

11340IF y1=FALSE THEN P=0:Y=Ytmin

11350IFy2>0 P=0 ELSE P=pen

11360PROCdraw(FNx(X),Y,P)

11370NEXT X:P=0:PROCdraw(FNx(X),Y,P)

11380ENDPROC

11390DEF PROCstring(x,y,string$)

11400MOVE FNx(x),FNy(y)

11410VDU5

11420PRINT string$

11430VDU4

11440ENDPROC

11450DEF PROCerror

11460IF ERR<>18 AND ERR<>20 THEN 11510

11470IF ERL<>11030 OR L%=0 THEN 11510

11480Xsmin2=X+(Srx/Prx)

11490PROCplot(pen,equ$)

11500ENDPROC

11510CLS:PRINT '"EXECUTION HALTED, there is an error" 

11520@%=0:PRINT'"at line "ERL" It is:-"'

11530REPORT

11540PRINT ''"Please record how you 'crashed' the" 

11550PRINT "program and contact the author."

11560PROCgstop:STOP

11570DEF PROCbestfit(name$,alterX$,alterY$)

11580LOCAL sigX,sigY,sigXY,sigX2,sigY2,top,base,X%

11590PROCfile(TRUE,name$)

11600PROCdata(1,0,0,FALSE,alterX$,alterY$)

11620FOR X%=1 TO total

11630PROCdata(X%,0,0,FALSE,alterX$,alterY$)

11680sigX=sigX+X

11690sigY=sigY+Y

11700sigXY=sigXY+X*Y

11710sigX2=sigX2+X^2

11720sigY2=sigY2+Y^2

11730NEXT X%

11740PROCfile(FALSE,name$)

11750m=(sigXY-(sigX*sigY/total))/(sigX2-sigX^2/total)

11760c=sigY/total-m*sigX/total

11761top=total*sigXY-sigX*sigY

11762base=SQR((total*sigX2-sigX^2)*(total*sigY2-sigY^2)) 

11763r=top/base

11770ENDPROC

>*SPOOL


00000000  3e 4c 2e 0a 0d 31 30 30  30 30 44 45 46 20 50 52  |>L...10000DEF PR|
00000010  4f 43 66 69 6c 65 28 6f  70 65 6e 2c 6e 61 6d 65  |OCfile(open,name|
00000020  24 29 0a 0d 31 30 30 31  30 49 46 20 4e 4f 54 20  |$)..10010IF NOT |
00000030  6f 70 65 6e 20 54 48 45  4e 20 50 54 52 23 45 25  |open THEN PTR#E%|
00000040  3d 30 3a 50 52 49 4e 54  23 45 25 2c 74 6f 74 61  |=0:PRINT#E%,tota|
00000050  6c 3a 43 4c 4f 53 45 23  30 3a 45 4e 44 50 52 4f  |l:CLOSE#0:ENDPRO|
00000060  43 0a 0d 31 30 30 32 30  43 4c 4f 53 45 23 30 0a  |C..10020CLOSE#0.|
00000070  0d 31 30 30 33 30 45 25  3d 4f 50 45 4e 55 50 28  |.10030E%=OPENUP(|
00000080  6e 61 6d 65 24 29 0a 0d  31 30 30 34 30 49 46 20  |name$)..10040IF |
00000090  45 25 3d 30 20 43 4c 4f  53 45 23 30 3a 45 25 3d  |E%=0 CLOSE#0:E%=|
000000a0  4f 50 45 4e 4f 55 54 28  6e 61 6d 65 24 29 3a 43  |OPENOUT(name$):C|
000000b0  4c 4f 53 45 23 30 3a 45  25 3d 4f 50 45 4e 55 50  |LOSE#0:E%=OPENUP|
000000c0  28 6e 61 6d 65 24 29 3a  74 6f 74 61 6c 3d 31 3a  |(name$):total=1:|
000000d0  50 52 49 4e 54 23 45 25  2c 74 6f 74 61 6c 2c 45  |PRINT#E%,total,E|
000000e0  25 2c 45 25 3a 45 4e 44  50 52 4f 43 0a 0d 31 30  |%,E%:ENDPROC..10|
000000f0  30 35 30 49 46 20 45 25  3c 3e 30 20 49 4e 50 55  |050IF E%<>0 INPU|
00000100  54 23 45 25 2c 74 6f 74  61 6c 0a 0d 31 30 30 36  |T#E%,total..1006|
00000110  30 45 4e 44 50 52 4f 43  0a 0d 31 30 30 37 30 44  |0ENDPROC..10070D|
00000120  45 46 20 50 52 4f 43 64  61 74 61 28 70 61 69 72  |EF PROCdata(pair|
00000130  2c 78 2c 79 2c 77 72 69  74 65 2c 61 6c 74 65 72  |,x,y,write,alter|
00000140  58 24 2c 61 6c 74 65 72  59 24 29 0a 0d 31 30 30  |X$,alterY$)..100|
00000150  38 30 4c 4f 43 41 4c 20  5a 25 2c 78 6e 0a 0d 31  |80LOCAL Z%,xn..1|
00000160  30 30 39 30 49 46 20 70  61 69 72 3c 31 20 54 48  |0090IF pair<1 TH|
00000170  45 4e 20 50 52 4f 43 66  69 6c 65 28 46 41 4c 53  |EN PROCfile(FALS|
00000180  45 2c 22 22 29 3a 50 52  49 4e 54 20 22 45 52 52  |E,""):PRINT "ERR|
00000190  4f 52 20 50 41 49 52 20  3c 20 31 22 3a 45 4e 44  |OR PAIR < 1":END|
000001a0  0a 0d 31 30 31 30 30 49  46 20 70 61 69 72 3e 74  |..10100IF pair>t|
000001b0  6f 74 61 6c 20 54 48 45  4e 20 50 54 52 23 45 25  |otal THEN PTR#E%|
000001c0  3d 31 32 2a 74 6f 74 61  6c 2d 36 3a 49 4e 50 55  |=12*total-6:INPU|
000001d0  54 23 45 25 2c 5a 25 2c  5a 25 3a 46 4f 52 20 5a  |T#E%,Z%,Z%:FOR Z|
000001e0  25 3d 74 6f 74 61 6c 2b  31 20 54 4f 20 70 61 69  |%=total+1 TO pai|
000001f0  72 3a 50 52 49 4e 54 23  45 25 2c 5a 25 2c 5a 25  |r:PRINT#E%,Z%,Z%|
00000200  3a 4e 45 58 54 0a 0d 31  30 31 31 30 50 54 52 23  |:NEXT..10110PTR#|
00000210  45 25 3d 31 32 2a 70 61  69 72 2d 36 0a 0d 31 30  |E%=12*pair-6..10|
00000220  31 32 30 49 46 20 77 72  69 74 65 20 54 48 45 4e  |120IF write THEN|
00000230  20 50 52 49 4e 54 23 45  25 2c 78 2c 79 20 45 4c  | PRINT#E%,x,y EL|
00000240  53 45 20 49 4e 50 55 54  23 45 25 2c 58 2c 59 0a  |SE INPUT#E%,X,Y.|
00000250  0d 31 30 31 33 30 49 46  20 70 61 69 72 3e 74 6f  |.10130IF pair>to|
00000260  74 61 6c 20 54 48 45 4e  20 74 6f 74 61 6c 3d 70  |tal THEN total=p|
00000270  61 69 72 0a 0d 31 30 31  34 30 49 46 20 4e 4f 54  |air..10140IF NOT|
00000280  20 77 72 69 74 65 20 78  6e 3d 45 56 41 4c 28 61  | write xn=EVAL(a|
00000290  6c 74 65 72 58 24 29 3a  59 3d 45 56 41 4c 28 61  |lterX$):Y=EVAL(a|
000002a0  6c 74 65 72 59 24 29 3a  58 3d 78 6e 0a 0d 31 30  |lterY$):X=xn..10|
000002b0  31 35 30 45 4e 44 50 52  4f 43 0a 0d 31 30 31 36  |150ENDPROC..1016|
000002c0  30 44 45 46 20 50 52 4f  43 67 73 74 61 72 74 0a  |0DEF PROCgstart.|
000002d0  0d 31 30 31 37 30 43 4c  47 0a 0d 31 30 31 38 30  |.10170CLG..10180|
000002e0  50 52 4f 43 70 6f 73 28  54 52 55 45 2c 30 2c 30  |PROCpos(TRUE,0,0|
000002f0  2c 30 2c 30 2c 30 2c 30  29 0a 0d 31 30 31 39 30  |,0,0,0,0)..10190|
00000300  50 52 4f 43 73 63 61 28  54 52 55 45 2c 30 2c 30  |PROCsca(TRUE,0,0|
00000310  2c 30 2c 30 2c 30 2c 30  29 0a 0d 31 30 32 30 30  |,0,0,0,0)..10200|
00000320  43 4c 4f 53 45 23 30 0a  0d 31 30 32 31 30 4f 4e  |CLOSE#0..10210ON|
00000330  20 45 52 52 4f 52 20 50  52 4f 43 65 72 72 6f 72  | ERROR PROCerror|
00000340  3a 47 4f 54 4f 20 4c 25  0a 0d 31 30 32 32 30 4c  |:GOTO L%..10220L|
00000350  25 3d 30 0a 0d 31 30 32  33 30 45 4e 44 50 52 4f  |%=0..10230ENDPRO|
00000360  43 0a 0d 31 30 32 34 30  44 45 46 20 50 52 4f 43  |C..10240DEF PROC|
00000370  67 73 74 6f 70 0a 0d 31  30 32 35 30 43 4c 4f 53  |gstop..10250CLOS|
00000380  45 23 30 0a 0d 31 30 32  36 30 45 4e 44 50 52 4f  |E#0..10260ENDPRO|
00000390  43 0a 0d 31 30 32 37 30  44 45 46 20 50 52 4f 43  |C..10270DEF PROC|
000003a0  70 6f 73 28 64 65 66 61  75 6c 74 2c 58 74 6d 69  |pos(default,Xtmi|
000003b0  2c 58 74 6d 61 2c 78 73  2c 59 74 6d 69 2c 59 74  |,Xtma,xs,Ytmi,Yt|
000003c0  6d 61 2c 79 73 29 0a 0d  31 30 32 38 30 49 46 20  |ma,ys)..10280IF |
000003d0  64 65 66 61 75 6c 74 20  54 48 45 4e 20 58 74 6d  |default THEN Xtm|
000003e0  69 6e 3d 31 35 3a 58 74  6d 61 78 3d 31 32 37 39  |in=15:Xtmax=1279|
000003f0  3a 78 73 6c 3d 31 35 3a  59 74 6d 69 6e 3d 31 35  |:xsl=15:Ytmin=15|
00000400  3a 59 74 6d 61 78 3d 31  30 32 33 3a 79 73 6c 3d  |:Ytmax=1023:ysl=|
00000410  31 35 0a 0d 31 30 32 39  30 49 46 20 4e 4f 54 20  |15..10290IF NOT |
00000420  64 65 66 61 75 6c 74 20  54 48 45 4e 20 58 74 6d  |default THEN Xtm|
00000430  69 6e 3d 58 74 6d 69 3a  58 74 6d 61 78 3d 58 74  |in=Xtmi:Xtmax=Xt|
00000440  6d 61 3a 78 73 6c 3d 78  73 3a 59 74 6d 69 6e 3d  |ma:xsl=xs:Ytmin=|
00000450  59 74 6d 69 3a 59 74 6d  61 78 3d 59 74 6d 61 3a  |Ytmi:Ytmax=Ytma:|
00000460  79 73 6c 3d 79 73 0a 0d  31 30 33 30 30 50 72 79  |ysl=ys..10300Pry|
00000470  3d 59 74 6d 61 78 2d 59  74 6d 69 6e 0a 0d 31 30  |=Ytmax-Ytmin..10|
00000480  33 31 30 50 72 78 3d 58  74 6d 61 78 2d 58 74 6d  |310Prx=Xtmax-Xtm|
00000490  69 6e 0a 0d 31 30 33 32  30 45 4e 44 50 52 4f 43  |in..10320ENDPROC|
000004a0  0a 0d 31 30 33 33 30 44  45 46 20 50 52 4f 43 73  |..10330DEF PROCs|
000004b0  63 61 28 64 65 66 61 75  6c 74 2c 58 73 6d 69 2c  |ca(default,Xsmi,|
000004c0  58 73 6d 61 2c 73 78 2c  59 73 6d 69 2c 59 73 6d  |Xsma,sx,Ysmi,Ysm|
000004d0  61 2c 73 79 29 0a 0d 31  30 33 34 30 49 46 20 64  |a,sy)..10340IF d|
000004e0  65 66 61 75 6c 74 20 54  48 45 4e 20 58 73 6d 69  |efault THEN Xsmi|
000004f0  6e 3d 2d 31 3a 58 73 6d  61 78 3d 31 3a 73 74 78  |n=-1:Xsmax=1:stx|
00000500  3d 30 2e 31 3a 59 73 6d  69 6e 3d 2d 31 3a 59 73  |=0.1:Ysmin=-1:Ys|
00000510  6d 61 78 3d 31 3a 73 74  79 3d 30 2e 31 0a 0d 31  |max=1:sty=0.1..1|
00000520  30 33 35 30 49 46 20 4e  4f 54 20 64 65 66 61 75  |0350IF NOT defau|
00000530  6c 74 20 54 48 45 4e 20  58 73 6d 69 6e 3d 58 73  |lt THEN Xsmin=Xs|
00000540  6d 69 3a 58 73 6d 61 78  3d 58 73 6d 61 3a 73 74  |mi:Xsmax=Xsma:st|
00000550  78 3d 73 78 3a 59 73 6d  69 6e 3d 59 73 6d 69 3a  |x=sx:Ysmin=Ysmi:|
00000560  59 73 6d 61 78 3d 59 73  6d 61 3a 73 74 79 3d 73  |Ysmax=Ysma:sty=s|
00000570  79 0a 0d 31 30 33 36 30  50 52 4f 43 73 63 61 6c  |y..10360PROCscal|
00000580  65 76 61 72 73 0a 0d 31  30 33 37 30 45 4e 44 50  |evars..10370ENDP|
00000590  52 4f 43 0a 0d 31 30 33  38 30 44 45 46 20 50 52  |ROC..10380DEF PR|
000005a0  4f 43 41 75 74 6f 53 63  61 6c 65 73 28 6e 61 6d  |OCAutoScales(nam|
000005b0  65 24 2c 61 6c 74 65 72  58 24 2c 61 6c 74 65 72  |e$,alterX$,alter|
000005c0  59 24 29 0a 0d 31 30 33  39 30 4c 4f 43 41 4c 20  |Y$)..10390LOCAL |
000005d0  52 61 6e 67 65 58 2c 52  61 6e 67 65 59 2c 50 77  |RangeX,RangeY,Pw|
000005e0  78 25 2c 50 77 79 25 2c  58 25 0a 0d 31 30 34 30  |x%,Pwy%,X%..1040|
000005f0  30 50 52 4f 43 66 69 6c  65 28 54 52 55 45 2c 6e  |0PROCfile(TRUE,n|
00000600  61 6d 65 24 29 0a 0d 31  30 34 31 30 50 52 4f 43  |ame$)..10410PROC|
00000610  64 61 74 61 28 31 2c 30  2c 30 2c 46 41 4c 53 45  |data(1,0,0,FALSE|
00000620  2c 61 6c 74 65 72 58 24  2c 61 6c 74 65 72 59 24  |,alterX$,alterY$|
00000630  29 0a 0d 31 30 34 32 30  58 73 6d 61 78 3d 58 3a  |)..10420Xsmax=X:|
00000640  58 73 6d 69 6e 3d 58 0a  0d 31 30 34 33 30 59 73  |Xsmin=X..10430Ys|
00000650  6d 69 6e 3d 59 3a 59 73  6d 61 78 3d 59 0a 0d 31  |min=Y:Ysmax=Y..1|
00000660  30 34 34 30 46 4f 52 20  58 25 3d 32 20 54 4f 20  |0440FOR X%=2 TO |
00000670  74 6f 74 61 6c 0a 0d 31  30 34 35 30 50 52 4f 43  |total..10450PROC|
00000680  64 61 74 61 28 58 25 2c  30 2c 30 2c 46 41 4c 53  |data(X%,0,0,FALS|
00000690  45 2c 61 6c 74 65 72 58  24 2c 61 6c 74 65 72 59  |E,alterX$,alterY|
000006a0  24 29 0a 0d 31 30 34 36  30 49 46 20 58 3e 58 73  |$)..10460IF X>Xs|
000006b0  6d 61 78 20 54 48 45 4e  20 58 73 6d 61 78 3d 58  |max THEN Xsmax=X|
000006c0  0a 0d 31 30 34 37 30 49  46 20 58 3c 58 73 6d 69  |..10470IF X<Xsmi|
000006d0  6e 20 54 48 45 4e 20 58  73 6d 69 6e 3d 58 0a 0d  |n THEN Xsmin=X..|
000006e0  31 30 34 38 30 49 46 20  59 3e 59 73 6d 61 78 20  |10480IF Y>Ysmax |
000006f0  54 48 45 4e 20 59 73 6d  61 78 3d 59 0a 0d 31 30  |THEN Ysmax=Y..10|
00000700  34 39 30 49 46 20 59 3c  59 73 6d 69 6e 20 54 48  |490IF Y<Ysmin TH|
00000710  45 4e 20 59 73 6d 69 6e  3d 59 0a 0d 31 30 35 30  |EN Ysmin=Y..1050|
00000720  30 4e 45 58 54 20 58 25  0a 0d 31 30 35 31 30 52  |0NEXT X%..10510R|
00000730  61 6e 67 65 58 3d 58 73  6d 61 78 2d 58 73 6d 69  |angeX=Xsmax-Xsmi|
00000740  6e 0a 0d 31 30 35 32 30  52 61 6e 67 65 59 3d 59  |n..10520RangeY=Y|
00000750  73 6d 61 78 2d 59 73 6d  69 6e 0a 0d 31 30 35 33  |smax-Ysmin..1053|
00000760  30 50 77 78 25 3d 46 4e  70 6f 77 65 72 28 52 61  |0Pwx%=FNpower(Ra|
00000770  6e 67 65 58 29 0a 0d 31  30 35 34 30 50 77 79 25  |ngeX)..10540Pwy%|
00000780  3d 46 4e 70 6f 77 65 72  28 52 61 6e 67 65 59 29  |=FNpower(RangeY)|
00000790  0a 0d 31 30 35 35 30 58  73 6d 61 78 3d 28 49 4e  |..10550Xsmax=(IN|
000007a0  54 28 58 73 6d 61 78 2f  31 30 5e 50 77 78 25 2b  |T(Xsmax/10^Pwx%+|
000007b0  31 29 29 2a 31 30 5e 50  77 78 25 0a 0d 31 30 35  |1))*10^Pwx%..105|
000007c0  36 30 58 73 6d 69 6e 3d  28 49 4e 54 28 58 73 6d  |60Xsmin=(INT(Xsm|
000007d0  69 6e 2f 31 30 5e 50 77  78 25 29 29 2a 31 30 5e  |in/10^Pwx%))*10^|
000007e0  50 77 78 25 0a 0d 31 30  35 37 30 59 73 6d 61 78  |Pwx%..10570Ysmax|
000007f0  3d 28 49 4e 54 28 59 73  6d 61 78 2f 31 30 5e 50  |=(INT(Ysmax/10^P|
00000800  77 79 25 2b 31 29 29 2a  31 30 5e 50 77 79 25 0a  |wy%+1))*10^Pwy%.|
00000810  0d 31 30 35 38 30 59 73  6d 69 6e 3d 28 49 4e 54  |.10580Ysmin=(INT|
00000820  28 59 73 6d 69 6e 2f 31  30 5e 50 77 79 25 29 29  |(Ysmin/10^Pwy%))|
00000830  2a 31 30 5e 50 77 79 25  0a 0d 31 30 35 39 30 73  |*10^Pwy%..10590s|
00000840  74 78 3d 46 4e 73 74 65  70 73 69 7a 65 28 58 73  |tx=FNstepsize(Xs|
00000850  6d 61 78 2d 58 73 6d 69  6e 29 0a 0d 31 30 36 30  |max-Xsmin)..1060|
00000860  30 73 74 79 3d 46 4e 73  74 65 70 73 69 7a 65 28  |0sty=FNstepsize(|
00000870  59 73 6d 61 78 2d 59 73  6d 69 6e 29 0a 0d 31 30  |Ysmax-Ysmin)..10|
00000880  36 31 30 50 52 4f 43 73  63 61 6c 65 76 61 72 73  |610PROCscalevars|
00000890  0a 0d 31 30 36 32 30 45  4e 44 50 52 4f 43 0a 0d  |..10620ENDPROC..|
000008a0  31 30 36 33 30 44 45 46  20 46 4e 70 6f 77 65 72  |10630DEF FNpower|
000008b0  28 6e 75 6d 62 65 72 29  0a 0d 31 30 36 34 30 3d  |(number)..10640=|
000008c0  49 4e 54 28 4c 4f 47 28  6e 75 6d 62 65 72 29 29  |INT(LOG(number))|
000008d0  0a 0d 31 30 36 35 30 44  45 46 20 46 4e 73 74 65  |..10650DEF FNste|
000008e0  70 73 69 7a 65 28 72 61  6e 67 65 29 0a 0d 31 30  |psize(range)..10|
000008f0  36 36 30 4c 4f 43 41 4c  20 70 6f 77 65 72 25 20  |660LOCAL power% |
00000900  0a 0d 31 30 36 37 30 70  6f 77 65 72 25 3d 46 4e  |..10670power%=FN|
00000910  70 6f 77 65 72 28 72 61  6e 67 65 29 0a 0d 31 30  |power(range)..10|
00000920  36 38 30 72 61 6e 67 65  3d 72 61 6e 67 65 2f 31  |680range=range/1|
00000930  30 5e 70 6f 77 65 72 25  0a 0d 31 30 36 39 30 49  |0^power%..10690I|
00000940  46 20 72 61 6e 67 65 3c  3d 33 20 54 48 45 4e 20  |F range<=3 THEN |
00000950  3d 32 2a 31 30 5e 28 70  6f 77 65 72 25 2d 31 29  |=2*10^(power%-1)|
00000960  0a 0d 31 30 37 30 30 49  46 20 72 61 6e 67 65 3e  |..10700IF range>|
00000970  33 20 41 4e 44 20 72 61  6e 67 65 3c 37 20 54 48  |3 AND range<7 TH|
00000980  45 4e 20 3d 35 2a 31 30  5e 28 70 6f 77 65 72 25  |EN =5*10^(power%|
00000990  2d 31 29 0a 0d 31 30 37  31 30 3d 31 2a 31 30 5e  |-1)..10710=1*10^|
000009a0  70 6f 77 65 72 25 0a 0d  31 30 37 32 30 44 45 46  |power%..10720DEF|
000009b0  20 50 52 4f 43 73 63 61  6c 65 76 61 72 73 0a 0d  | PROCscalevars..|
000009c0  31 30 37 33 30 53 72 79  3d 59 73 6d 61 78 2d 59  |10730Sry=Ysmax-Y|
000009d0  73 6d 69 6e 0a 0d 31 30  37 34 30 59 30 3d 59 74  |smin..10740Y0=Yt|
000009e0  6d 61 78 2d 50 72 79 2a  59 73 6d 61 78 2f 53 72  |max-Pry*Ysmax/Sr|
000009f0  79 0a 0d 31 30 37 35 30  53 72 78 3d 58 73 6d 61  |y..10750Srx=Xsma|
00000a00  78 2d 58 73 6d 69 6e 0a  0d 31 30 37 36 30 58 30  |x-Xsmin..10760X0|
00000a10  3d 58 74 6d 61 78 2d 50  72 78 2a 58 73 6d 61 78  |=Xtmax-Prx*Xsmax|
00000a20  2f 53 72 78 0a 0d 31 30  37 37 30 58 73 6d 69 6e  |/Srx..10770Xsmin|
00000a30  32 3d 58 73 6d 69 6e 0a  0d 31 30 37 38 30 45 4e  |2=Xsmin..10780EN|
00000a40  44 50 52 4f 43 0a 0d 31  30 37 39 30 44 45 46 20  |DPROC..10790DEF |
00000a50  50 52 4f 43 44 72 61 77  66 69 6c 65 28 6e 61 6d  |PROCDrawfile(nam|
00000a60  65 24 2c 70 65 6e 2c 61  6c 74 65 72 58 24 2c 61  |e$,pen,alterX$,a|
00000a70  6c 74 65 72 59 24 2c 63  72 6f 73 73 65 73 29 0a  |lterY$,crosses).|
00000a80  0d 31 30 38 30 30 4c 4f  43 41 4c 20 50 45 4e 2c  |.10800LOCAL PEN,|
00000a90  58 25 2c 50 0a 0d 31 30  38 31 30 50 52 4f 43 66  |X%,P..10810PROCf|
00000aa0  69 6c 65 28 54 52 55 45  2c 6e 61 6d 65 24 29 0a  |ile(TRUE,name$).|
00000ab0  0d 31 30 38 32 30 46 4f  52 20 58 25 3d 31 20 54  |.10820FOR X%=1 T|
00000ac0  4f 20 74 6f 74 61 6c 0a  0d 31 30 38 33 30 50 52  |O total..10830PR|
00000ad0  4f 43 64 61 74 61 28 58  25 2c 30 2c 30 2c 46 41  |OCdata(X%,0,0,FA|
00000ae0  4c 53 45 2c 61 6c 74 65  72 58 24 2c 61 6c 74 65  |LSE,alterX$,alte|
00000af0  72 59 24 29 0a 0d 31 30  38 34 30 59 3d 46 4e 79  |rY$)..10840Y=FNy|
00000b00  28 59 29 3a 58 3d 46 4e  78 28 58 29 0a 0d 31 30  |(Y):X=FNx(X)..10|
00000b10  38 35 30 49 46 20 58 25  3d 31 20 54 48 45 4e 20  |850IF X%=1 THEN |
00000b20  50 45 4e 3d 30 20 45 4c  53 45 20 50 45 4e 3d 70  |PEN=0 ELSE PEN=p|
00000b30  65 6e 20 0a 0d 31 30 38  36 30 49 46 20 4e 4f 54  |en ..10860IF NOT|
00000b40  20 63 72 6f 73 73 65 73  20 54 48 45 4e 20 50 52  | crosses THEN PR|
00000b50  4f 43 64 72 61 77 28 58  2c 59 2c 50 45 4e 29 3a  |OCdraw(X,Y,PEN):|
00000b60  47 4f 54 4f 31 30 39 35  30 0a 0d 31 30 38 37 30  |GOTO10950..10870|
00000b70  50 3d 30 0a 0d 31 30 38  38 30 50 52 4f 43 64 72  |P=0..10880PROCdr|
00000b80  61 77 28 58 2d 31 30 2c  59 2c 50 29 0a 0d 31 30  |aw(X-10,Y,P)..10|
00000b90  38 39 30 50 3d 70 65 6e  0a 0d 31 30 39 30 30 50  |890P=pen..10900P|
00000ba0  52 4f 43 64 72 61 77 28  58 2b 31 30 2c 59 2c 50  |ROCdraw(X+10,Y,P|
00000bb0  29 0a 0d 31 30 39 31 30  50 3d 30 0a 0d 31 30 39  |)..10910P=0..109|
00000bc0  32 30 50 52 4f 43 64 72  61 77 28 58 2c 59 2d 31  |20PROCdraw(X,Y-1|
00000bd0  30 2c 50 29 0a 0d 31 30  39 33 30 50 3d 70 65 6e  |0,P)..10930P=pen|
00000be0  0a 0d 31 30 39 34 30 50  52 4f 43 64 72 61 77 28  |..10940PROCdraw(|
00000bf0  58 2c 59 2b 31 30 2c 50  29 0a 0d 31 30 39 35 30  |X,Y+10,P)..10950|
00000c00  4e 45 58 54 20 58 25 0a  0d 31 30 39 36 30 50 52  |NEXT X%..10960PR|
00000c10  4f 43 66 69 6c 65 28 46  41 4c 53 45 2c 6e 61 6d  |OCfile(FALSE,nam|
00000c20  65 24 29 0a 0d 31 30 39  37 30 45 4e 44 50 52 4f  |e$)..10970ENDPRO|
00000c30  43 0a 0d 31 30 39 38 30  44 45 46 20 50 52 4f 43  |C..10980DEF PROC|
00000c40  64 72 61 77 28 78 2c 79  2c 70 65 6e 29 0a 0d 31  |draw(x,y,pen)..1|
00000c50  30 39 39 30 49 46 20 70  65 6e 3d 30 20 4d 4f 56  |0990IF pen=0 MOV|
00000c60  45 20 78 2c 79 20 45 4c  53 45 20 47 43 4f 4c 30  |E x,y ELSE GCOL0|
00000c70  2c 70 65 6e 3a 44 52 41  57 78 2c 79 20 0a 0d 31  |,pen:DRAWx,y ..1|
00000c80  31 30 30 30 45 4e 44 50  52 4f 43 0a 0d 31 31 30  |1000ENDPROC..110|
00000c90  31 30 44 45 46 20 46 4e  78 28 58 29 0a 0d 31 31  |10DEF FNx(X)..11|
00000ca0  30 32 30 3d 28 50 72 78  2f 53 72 78 29 2a 58 2b  |020=(Prx/Srx)*X+|
00000cb0  58 30 0a 0d 31 31 30 33  30 44 45 46 20 46 4e 79  |X0..11030DEF FNy|
00000cc0  28 59 29 0a 0d 31 31 30  34 30 3d 28 50 72 79 2f  |(Y)..11040=(Pry/|
00000cd0  53 72 79 29 2a 59 2b 59  30 0a 0d 31 31 30 35 30  |Sry)*Y+Y0..11050|
00000ce0  44 45 46 20 50 52 4f 43  73 63 61 6c 65 28 70 65  |DEF PROCscale(pe|
00000cf0  6e 29 0a 0d 31 31 30 36  30 4c 4f 43 41 4c 20 58  |n)..11060LOCAL X|
00000d00  2c 58 32 2c 59 2c 59 32  2c 50 0a 0d 31 31 30 37  |,X2,Y,Y2,P..1107|
00000d10  30 49 46 20 59 30 3e 59  74 6d 61 78 20 54 48 45  |0IF Y0>Ytmax THE|
00000d20  4e 20 59 32 3d 59 74 6d  61 78 3a 47 4f 54 4f 31  |N Y2=Ytmax:GOTO1|
00000d30  31 31 32 30 0a 0d 31 31  30 38 30 49 46 20 59 30  |1120..11080IF Y0|
00000d40  3c 59 74 6d 69 6e 20 54  48 45 4e 20 59 32 3d 59  |<Ytmin THEN Y2=Y|
00000d50  74 6d 69 6e 3a 47 4f 54  4f 31 31 31 32 30 0a 0d  |tmin:GOTO11120..|
00000d60  31 31 30 39 30 59 32 3d  59 30 0a 0d 31 31 31 30  |11090Y2=Y0..1110|
00000d70  30 58 3d 58 74 6d 69 6e  3a 59 3d 59 32 3a 50 3d  |0X=Xtmin:Y=Y2:P=|
00000d80  30 3a 50 52 4f 43 64 72  61 77 28 58 2c 59 2c 50  |0:PROCdraw(X,Y,P|
00000d90  29 0a 0d 31 31 31 31 30  58 3d 58 74 6d 61 78 3a  |)..11110X=Xtmax:|
00000da0  50 3d 70 65 6e 3a 50 52  4f 43 64 72 61 77 28 58  |P=pen:PROCdraw(X|
00000db0  2c 59 2c 50 29 0a 0d 31  31 31 32 30 49 46 20 58  |,Y,P)..11120IF X|
00000dc0  30 3e 58 74 6d 61 78 20  54 48 45 4e 20 58 32 3d  |0>Xtmax THEN X2=|
00000dd0  58 74 6d 61 78 3a 47 4f  54 4f 31 31 31 37 30 0a  |Xtmax:GOTO11170.|
00000de0  0d 31 31 31 33 30 49 46  20 58 30 3c 58 74 6d 69  |.11130IF X0<Xtmi|
00000df0  6e 20 54 48 45 4e 20 58  32 3d 58 74 6d 69 6e 3a  |n THEN X2=Xtmin:|
00000e00  47 4f 54 4f 31 31 31 37  30 0a 0d 31 31 31 34 30  |GOTO11170..11140|
00000e10  58 32 3d 58 30 0a 0d 31  31 31 35 30 58 3d 58 32  |X2=X0..11150X=X2|
00000e20  3a 59 3d 59 74 6d 69 6e  3a 50 3d 30 3a 50 52 4f  |:Y=Ytmin:P=0:PRO|
00000e30  43 64 72 61 77 28 58 2c  59 2c 50 29 0a 0d 31 31  |Cdraw(X,Y,P)..11|
00000e40  31 36 30 59 3d 59 74 6d  61 78 3a 50 3d 70 65 6e  |160Y=Ytmax:P=pen|
00000e50  3a 50 52 4f 43 64 72 61  77 28 58 2c 59 2c 50 29  |:PROCdraw(X,Y,P)|
00000e60  0a 0d 31 31 31 37 30 46  4f 52 20 58 3d 58 73 6d  |..11170FOR X=Xsm|
00000e70  61 78 20 54 4f 20 58 73  6d 69 6e 20 53 54 45 50  |ax TO Xsmin STEP|
00000e80  20 2d 31 2a 73 74 78 0a  0d 31 31 31 38 30 78 31  | -1*stx..11180x1|
00000e90  3d 46 4e 78 28 58 29 3a  59 3d 59 32 3a 50 3d 30  |=FNx(X):Y=Y2:P=0|
00000ea0  0a 0d 31 31 31 39 30 50  52 4f 43 64 72 61 77 28  |..11190PROCdraw(|
00000eb0  78 31 2c 59 2c 50 29 3a  49 46 20 58 3d 30 20 54  |x1,Y,P):IF X=0 T|
00000ec0  48 45 4e 20 50 3d 30 20  45 4c 53 45 20 50 3d 70  |HEN P=0 ELSE P=p|
00000ed0  65 6e 0a 0d 31 31 32 30  30 59 3d 59 32 2d 78 73  |en..11200Y=Y2-xs|
00000ee0  6c 3a 50 52 4f 43 64 72  61 77 28 78 31 2c 59 2c  |l:PROCdraw(x1,Y,|
00000ef0  50 29 0a 0d 31 31 32 31  30 4e 45 58 54 20 58 0a  |P)..11210NEXT X.|
00000f00  0d 31 31 32 32 30 46 4f  52 20 59 3d 59 73 6d 69  |.11220FOR Y=Ysmi|
00000f10  6e 20 54 4f 20 59 73 6d  61 78 20 53 54 45 50 20  |n TO Ysmax STEP |
00000f20  73 74 79 0a 0d 31 31 32  33 30 79 31 3d 46 4e 79  |sty..11230y1=FNy|
00000f30  28 59 29 3a 58 3d 58 32  3a 50 3d 30 0a 0d 31 31  |(Y):X=X2:P=0..11|
00000f40  32 34 30 50 52 4f 43 64  72 61 77 28 58 2c 79 31  |240PROCdraw(X,y1|
00000f50  2c 50 29 3a 49 46 20 59  3d 30 20 54 48 45 4e 20  |,P):IF Y=0 THEN |
00000f60  50 3d 30 20 45 4c 53 45  20 50 3d 70 65 6e 0a 0d  |P=0 ELSE P=pen..|
00000f70  31 31 32 35 30 58 3d 58  32 2d 79 73 6c 3a 50 52  |11250X=X2-ysl:PR|
00000f80  4f 43 64 72 61 77 28 58  2c 79 31 2c 50 29 0a 0d  |OCdraw(X,y1,P)..|
00000f90  31 31 32 36 30 4e 45 58  54 20 59 0a 0d 31 31 32  |11260NEXT Y..112|
00000fa0  37 30 45 4e 44 50 52 4f  43 0a 0d 31 31 32 38 30  |70ENDPROC..11280|
00000fb0  44 45 46 20 50 52 4f 43  70 6c 6f 74 28 70 65 6e  |DEF PROCplot(pen|
00000fc0  2c 65 71 75 24 29 0a 0d  31 31 32 39 30 4c 4f 43  |,equ$)..11290LOC|
00000fd0  41 4c 20 58 2c 79 31 2c  79 32 2c 50 0a 0d 31 31  |AL X,y1,y2,P..11|
00000fe0  33 30 30 79 32 3d 32 0a  0d 31 31 33 31 30 46 4f  |300y2=2..11310FO|
00000ff0  52 20 58 3d 58 73 6d 69  6e 32 20 54 4f 20 58 73  |R X=Xsmin2 TO Xs|
00001000  6d 61 78 20 53 54 45 50  20 53 72 78 2f 50 72 78  |max STEP Srx/Prx|
00001010  2a 31 30 0a 0d 31 31 33  32 30 59 3d 46 4e 79 28  |*10..11320Y=FNy(|
00001020  45 56 41 4c 28 65 71 75  24 29 29 20 0a 0d 31 31  |EVAL(equ$)) ..11|
00001030  33 33 30 49 46 20 59 3e  59 74 6d 61 78 20 4f 52  |330IF Y>Ytmax OR|
00001040  20 59 3c 59 74 6d 69 6e  20 54 48 45 4e 20 79 31  | Y<Ytmin THEN y1|
00001050  3d 46 41 4c 53 45 3a 79  32 3d 32 20 45 4c 53 45  |=FALSE:y2=2 ELSE|
00001060  20 79 31 3d 54 52 55 45  3a 79 32 3d 79 32 2d 31  | y1=TRUE:y2=y2-1|
00001070  0a 0d 31 31 33 34 30 49  46 20 79 31 3d 46 41 4c  |..11340IF y1=FAL|
00001080  53 45 20 54 48 45 4e 20  50 3d 30 3a 59 3d 59 74  |SE THEN P=0:Y=Yt|
00001090  6d 69 6e 0a 0d 31 31 33  35 30 49 46 79 32 3e 30  |min..11350IFy2>0|
000010a0  20 50 3d 30 20 45 4c 53  45 20 50 3d 70 65 6e 0a  | P=0 ELSE P=pen.|
000010b0  0d 31 31 33 36 30 50 52  4f 43 64 72 61 77 28 46  |.11360PROCdraw(F|
000010c0  4e 78 28 58 29 2c 59 2c  50 29 0a 0d 31 31 33 37  |Nx(X),Y,P)..1137|
000010d0  30 4e 45 58 54 20 58 3a  50 3d 30 3a 50 52 4f 43  |0NEXT X:P=0:PROC|
000010e0  64 72 61 77 28 46 4e 78  28 58 29 2c 59 2c 50 29  |draw(FNx(X),Y,P)|
000010f0  0a 0d 31 31 33 38 30 45  4e 44 50 52 4f 43 0a 0d  |..11380ENDPROC..|
00001100  31 31 33 39 30 44 45 46  20 50 52 4f 43 73 74 72  |11390DEF PROCstr|
00001110  69 6e 67 28 78 2c 79 2c  73 74 72 69 6e 67 24 29  |ing(x,y,string$)|
00001120  0a 0d 31 31 34 30 30 4d  4f 56 45 20 46 4e 78 28  |..11400MOVE FNx(|
00001130  78 29 2c 46 4e 79 28 79  29 0a 0d 31 31 34 31 30  |x),FNy(y)..11410|
00001140  56 44 55 35 0a 0d 31 31  34 32 30 50 52 49 4e 54  |VDU5..11420PRINT|
00001150  20 73 74 72 69 6e 67 24  0a 0d 31 31 34 33 30 56  | string$..11430V|
00001160  44 55 34 0a 0d 31 31 34  34 30 45 4e 44 50 52 4f  |DU4..11440ENDPRO|
00001170  43 0a 0d 31 31 34 35 30  44 45 46 20 50 52 4f 43  |C..11450DEF PROC|
00001180  65 72 72 6f 72 0a 0d 31  31 34 36 30 49 46 20 45  |error..11460IF E|
00001190  52 52 3c 3e 31 38 20 41  4e 44 20 45 52 52 3c 3e  |RR<>18 AND ERR<>|
000011a0  32 30 20 54 48 45 4e 20  31 31 35 31 30 0a 0d 31  |20 THEN 11510..1|
000011b0  31 34 37 30 49 46 20 45  52 4c 3c 3e 31 31 30 33  |1470IF ERL<>1103|
000011c0  30 20 4f 52 20 4c 25 3d  30 20 54 48 45 4e 20 31  |0 OR L%=0 THEN 1|
000011d0  31 35 31 30 0a 0d 31 31  34 38 30 58 73 6d 69 6e  |1510..11480Xsmin|
000011e0  32 3d 58 2b 28 53 72 78  2f 50 72 78 29 0a 0d 31  |2=X+(Srx/Prx)..1|
000011f0  31 34 39 30 50 52 4f 43  70 6c 6f 74 28 70 65 6e  |1490PROCplot(pen|
00001200  2c 65 71 75 24 29 0a 0d  31 31 35 30 30 45 4e 44  |,equ$)..11500END|
00001210  50 52 4f 43 0a 0d 31 31  35 31 30 43 4c 53 3a 50  |PROC..11510CLS:P|
00001220  52 49 4e 54 20 27 22 45  58 45 43 55 54 49 4f 4e  |RINT '"EXECUTION|
00001230  20 48 41 4c 54 45 44 2c  20 74 68 65 72 65 20 69  | HALTED, there i|
00001240  73 20 61 6e 20 65 72 72  6f 72 22 20 0a 0d 31 31  |s an error" ..11|
00001250  35 32 30 40 25 3d 30 3a  50 52 49 4e 54 27 22 61  |520@%=0:PRINT'"a|
00001260  74 20 6c 69 6e 65 20 22  45 52 4c 22 20 49 74 20  |t line "ERL" It |
00001270  69 73 3a 2d 22 27 0a 0d  31 31 35 33 30 52 45 50  |is:-"'..11530REP|
00001280  4f 52 54 0a 0d 31 31 35  34 30 50 52 49 4e 54 20  |ORT..11540PRINT |
00001290  27 27 22 50 6c 65 61 73  65 20 72 65 63 6f 72 64  |''"Please record|
000012a0  20 68 6f 77 20 79 6f 75  20 27 63 72 61 73 68 65  | how you 'crashe|
000012b0  64 27 20 74 68 65 22 20  0a 0d 31 31 35 35 30 50  |d' the" ..11550P|
000012c0  52 49 4e 54 20 22 70 72  6f 67 72 61 6d 20 61 6e  |RINT "program an|
000012d0  64 20 63 6f 6e 74 61 63  74 20 74 68 65 20 61 75  |d contact the au|
000012e0  74 68 6f 72 2e 22 0a 0d  31 31 35 36 30 50 52 4f  |thor."..11560PRO|
000012f0  43 67 73 74 6f 70 3a 53  54 4f 50 0a 0d 31 31 35  |Cgstop:STOP..115|
00001300  37 30 44 45 46 20 50 52  4f 43 62 65 73 74 66 69  |70DEF PROCbestfi|
00001310  74 28 6e 61 6d 65 24 2c  61 6c 74 65 72 58 24 2c  |t(name$,alterX$,|
00001320  61 6c 74 65 72 59 24 29  0a 0d 31 31 35 38 30 4c  |alterY$)..11580L|
00001330  4f 43 41 4c 20 73 69 67  58 2c 73 69 67 59 2c 73  |OCAL sigX,sigY,s|
00001340  69 67 58 59 2c 73 69 67  58 32 2c 73 69 67 59 32  |igXY,sigX2,sigY2|
00001350  2c 74 6f 70 2c 62 61 73  65 2c 58 25 0a 0d 31 31  |,top,base,X%..11|
00001360  35 39 30 50 52 4f 43 66  69 6c 65 28 54 52 55 45  |590PROCfile(TRUE|
00001370  2c 6e 61 6d 65 24 29 0a  0d 31 31 36 30 30 50 52  |,name$)..11600PR|
00001380  4f 43 64 61 74 61 28 31  2c 30 2c 30 2c 46 41 4c  |OCdata(1,0,0,FAL|
00001390  53 45 2c 61 6c 74 65 72  58 24 2c 61 6c 74 65 72  |SE,alterX$,alter|
000013a0  59 24 29 0a 0d 31 31 36  32 30 46 4f 52 20 58 25  |Y$)..11620FOR X%|
000013b0  3d 31 20 54 4f 20 74 6f  74 61 6c 0a 0d 31 31 36  |=1 TO total..116|
000013c0  33 30 50 52 4f 43 64 61  74 61 28 58 25 2c 30 2c  |30PROCdata(X%,0,|
000013d0  30 2c 46 41 4c 53 45 2c  61 6c 74 65 72 58 24 2c  |0,FALSE,alterX$,|
000013e0  61 6c 74 65 72 59 24 29  0a 0d 31 31 36 38 30 73  |alterY$)..11680s|
000013f0  69 67 58 3d 73 69 67 58  2b 58 0a 0d 31 31 36 39  |igX=sigX+X..1169|
00001400  30 73 69 67 59 3d 73 69  67 59 2b 59 0a 0d 31 31  |0sigY=sigY+Y..11|
00001410  37 30 30 73 69 67 58 59  3d 73 69 67 58 59 2b 58  |700sigXY=sigXY+X|
00001420  2a 59 0a 0d 31 31 37 31  30 73 69 67 58 32 3d 73  |*Y..11710sigX2=s|
00001430  69 67 58 32 2b 58 5e 32  0a 0d 31 31 37 32 30 73  |igX2+X^2..11720s|
00001440  69 67 59 32 3d 73 69 67  59 32 2b 59 5e 32 0a 0d  |igY2=sigY2+Y^2..|
00001450  31 31 37 33 30 4e 45 58  54 20 58 25 0a 0d 31 31  |11730NEXT X%..11|
00001460  37 34 30 50 52 4f 43 66  69 6c 65 28 46 41 4c 53  |740PROCfile(FALS|
00001470  45 2c 6e 61 6d 65 24 29  0a 0d 31 31 37 35 30 6d  |E,name$)..11750m|
00001480  3d 28 73 69 67 58 59 2d  28 73 69 67 58 2a 73 69  |=(sigXY-(sigX*si|
00001490  67 59 2f 74 6f 74 61 6c  29 29 2f 28 73 69 67 58  |gY/total))/(sigX|
000014a0  32 2d 73 69 67 58 5e 32  2f 74 6f 74 61 6c 29 0a  |2-sigX^2/total).|
000014b0  0d 31 31 37 36 30 63 3d  73 69 67 59 2f 74 6f 74  |.11760c=sigY/tot|
000014c0  61 6c 2d 6d 2a 73 69 67  58 2f 74 6f 74 61 6c 0a  |al-m*sigX/total.|
000014d0  0d 31 31 37 36 31 74 6f  70 3d 74 6f 74 61 6c 2a  |.11761top=total*|
000014e0  73 69 67 58 59 2d 73 69  67 58 2a 73 69 67 59 0a  |sigXY-sigX*sigY.|
000014f0  0d 31 31 37 36 32 62 61  73 65 3d 53 51 52 28 28  |.11762base=SQR((|
00001500  74 6f 74 61 6c 2a 73 69  67 58 32 2d 73 69 67 58  |total*sigX2-sigX|
00001510  5e 32 29 2a 28 74 6f 74  61 6c 2a 73 69 67 59 32  |^2)*(total*sigY2|
00001520  2d 73 69 67 59 5e 32 29  29 20 0a 0d 31 31 37 36  |-sigY^2)) ..1176|
00001530  33 72 3d 74 6f 70 2f 62  61 73 65 0a 0d 31 31 37  |3r=top/base..117|
00001540  37 30 45 4e 44 50 52 4f  43 0a 0d 3e 2a 53 50 4f  |70ENDPROC..>*SPO|
00001550  4f 4c 0a 0d                                       |OL..|
00001554
04-02-89/Alginon.m0
04-02-89/Alginon.m1
04-02-89/Alginon.m2
04-02-89/Alginon.m4
04-02-89/Alginon.m5

© Dominic Ford 2021–2024.

Our privacy policy is here.
Last updated: 28 Jul 2024, 23:31 UTC
Website designed by .