StarInfo/Amery/!CtrlLimit/!RunImage

Home » Archimedes archive » Acorn User » AU 1997-01 B.adf » Regulars » StarInfo/Amery/!CtrlLimit/!RunImage

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 » Acorn User » AU 1997-01 B.adf » Regulars
Filename:	StarInfo/Amery/!CtrlLimit/!RunImage
Read OK:	✔
File size:	06D3 bytes
Load address:	0000
Exec address:	0000

File contents

   10REM Central Limit Theorum
   20REM LEN1995 Jonathan Amery
   30REM Started 09/07/1995
   40REM Current Finished --/--/199-
   50:
   60ON ERROR REPORT:PRINTERL:END
   70VDU 23,180,0,0,127,204,204,120,0,0
   80DIM Data%(640)
   90ScaleFactor=16
  100Sz%=16
  110MODE 12
  120IF INKEY(15)<>13 THEN
  130PRINT "Rectangular Distribution"
  140PROCINPUT ("Sample size:",Sz%)
  150PRINT "Press any key to start"
  160M=0:n=0:Tot=0:Totq=0
  170G=GET
  180CLS
  190ELSE
  200QQ%=OPENIN"<Obey$Dir>.!Data"
  210INPUT#QQ%,ScaleFactor,Tot,Totq,M,n
  220ScaleFactor+=ScaleFactor
  230FOR N=0 TO 640
  240INPUT#QQ%,Data%(N)
  250NEXT
  260CLOSE#QQ%
  270PROCReplot
  280ENDIF
  290REPEATUNTILINKEY(30)<>13
  300REPEAT
  310T=0
  320FOR N=1 TO Sz%
  330n+=1:PRINT TAB(54,0);"Number of tests:";n
  340t=RND(640)-1
  350T+=t
  360NEXT
  370T=T/Sz%
  380T=INT(T+.5)
  390Tot+=T
  400Totq+=T^2
  410M+=1:PRINT TAB(0,0);"n:";M;TAB(30,0);"Sample Average:";T;"   "
  420PRINTTAB(0,1);"�:";Tot/M;TAB(68,1);"�:";SQR((Totq/M)-(Tot/M)^2)
  430Data%(T)+=1
  440IF Data%(T)>(256/ScaleFactor) THEN PROCReplot ELSE GCOL 4:LINE (T*2),0,(T*2),(Data%(T)*4*ScaleFactor)
  450QQQ=INKEY(0)
  460UNTIL QQQ=32 OR QQQ=13
  470IF QQQ=13 THEN S=TRUE ELSE S=FALSE
  480REPEATUNTIL INKEY$(50)<>" " AND INKEY(50)<>13
  490IF S THEN
  500QQ%=OPENOUT"<Obey$Dir>.!Data"
  510PRINT#QQ%,ScaleFactor,Tot,Totq,M,n
  520FOR N=0 TO 640
  530PRINT#QQ%,Data%(N)
  540NEXT
  550CLOSE#QQ%
  560ENDIF
  570V=(Totq/M)-(Tot/M)^2
  580E=Tot/M
  590Q=1/(SQR(2*PI*V))
  600GCOL 3
  610Z=M
  620MOVE 0,0
  630FOR N=0 TO 640
  640DRAW N*2,(Q*EXP((-1*((N-E)^2))/(2*V)))*Z*ScaleFactor*4
  650NEXT
  660REPEATUNTIL INKEY$(0)=" "
  670REPEATUNTIL INKEY$(50)<>" "
  680END
  690:
  700DEFPROCINPUT(S$,RETURN V)
  710LOCAL ERROR
  720PRINTS$;
  730OFF
  740ON ERROR LOCAL ERROR
  750REPEAT
  760IF V<>0 THEN PRINT STR$(V);"            ";
  770IF V<>0 THEN
  780FOR N=1 TO LEN (STR$(V)+"            ")
  790PRINT CHR$(8);
  800NEXT
  810ENDIF
  820IF V=0 THEN PRINT" ";CHR$(8);
  830G$=GET$
  840IF ASC(G$)=13 THEN UNTILTRUE:RESTORE ERROR:PRINT:ENDPROC
  850IF (G$<"0" OR G$>"9") AND ASC(G$)<>127 AND ASC(G$)<>8 THEN UNTIL FALSE
  860IF ASC(G$)<>127 AND ASC(G$)<>8 THEN V=V*10+EVAL(G$):UNTIL FALSE
  870V=V DIV 10
  880UNTIL FALSE
  890:
  900DEFPROCReplot
  910ScaleFactor=ScaleFactor/2
  920CLS
  930FOR T=0 TO 640
  940GCOL 4:LINE (T*2),0,(T*2),(Data%(T)*4*ScaleFactor)
  950NEXT
  960ENDPROC


� Central Limit Theorum
� �1995 Jonathan Amery
� Started 09/07/1995
(!� Current Finished --/--/199-
2:
<� � �:�:�
F$� 23,180,0,0,127,204,204,120,0,0
P� Data%(640)
ZScaleFactor=16
d
Sz%=16
n� 12
x� �(15)<>13 �
� � "Rectangular Distribution"
��INPUT ("Sample size:",Sz%)
�� "Press any key to start"
�M=0:n=0:Tot=0:Totq=0
�G=�
��
��
�QQ%=�"<Obey$Dir>.!Data"
�"�#QQ%,ScaleFactor,Tot,Totq,M,n
�ScaleFactor+=ScaleFactor
�� N=0 � 640
��#QQ%,Data%(N)
��
	�#QQ%
�Replot
�
"���(30)<>13
,�
6T=0
@� N=1 � Sz%
J&n+=1:� �54,0);"Number of tests:";n
Tt=�(640)-1
^T+=t
h�
rT=T/Sz%
|
T=�(T+.5)
�
Tot+=T
�
Totq+=T^2
�8M+=1:� �0,0);"n:";M;�30,0);"Sample Average:";T;"   "
�7�0,1);"�:";Tot/M;�68,1);"�:";�((Totq/M)-(Tot/M)^2)
�Data%(T)+=1
�Y� Data%(T)>(256/ScaleFactor) � �Replot � � 4:� (T*2),0,(T*2),(Data%(T)*4*ScaleFactor)
�QQQ=�(0)
�� QQQ=32 � QQQ=13
�� QQQ=13 � S=� � S=�
��� �(50)<>" " � �(50)<>13
�	� S �
�QQ%=�"<Obey$Dir>.!Data"
�"�#QQ%,ScaleFactor,Tot,Totq,M,n
� N=0 � 640
�#QQ%,Data%(N)
�
&	�#QQ%
0�
:V=(Totq/M)-(Tot/M)^2
DE=Tot/M
NQ=1/(�(2*�*V))
X� 3
bZ=M
l	� 0,0
v� N=0 � 640
�5� N*2,(Q*�((-1*((N-E)^2))/(2*V)))*Z*ScaleFactor*4
��
��� �(0)=" "
��� �(50)<>" "
��
�:
���INPUT(S$,� V)
�� �
��S$;
��
�� � � �
��
�#� V<>0 � � �(V);"            ";
� V<>0 �
#� N=1 � � (�(V)+"            ")
� �(8);
 �
*�
4� V=0 � �" ";�(8);
>G$=�
H� �(G$)=13 � ��:� �:�:�
R5� (G$<"0" � G$>"9") � �(G$)<>127 � �(G$)<>8 � � �
\.� �(G$)<>127 � �(G$)<>8 � V=V*10+�(G$):� �
fV=V � 10
p� �
z:
���Replot
�ScaleFactor=ScaleFactor/2
��
�� T=0 � 640
�0� 4:� (T*2),0,(T*2),(Data%(T)*4*ScaleFactor)
��
��
�

00000000  0d 00 0a 1b f4 20 43 65  6e 74 72 61 6c 20 4c 69  |..... Central Li|
00000010  6d 69 74 20 54 68 65 6f  72 75 6d 0d 00 14 1a f4  |mit Theorum.....|
00000020  20 a9 31 39 39 35 20 4a  6f 6e 61 74 68 61 6e 20  | .1995 Jonathan |
00000030  41 6d 65 72 79 0d 00 1e  18 f4 20 53 74 61 72 74  |Amery..... Start|
00000040  65 64 20 30 39 2f 30 37  2f 31 39 39 35 0d 00 28  |ed 09/07/1995..(|
00000050  21 f4 20 43 75 72 72 65  6e 74 20 46 69 6e 69 73  |!. Current Finis|
00000060  68 65 64 20 2d 2d 2f 2d  2d 2f 31 39 39 2d 0d 00  |hed --/--/199-..|
00000070  32 05 3a 0d 00 3c 0e ee  20 85 20 f6 3a f1 9e 3a  |2.:..<.. . .:..:|
00000080  e0 0d 00 46 24 ef 20 32  33 2c 31 38 30 2c 30 2c  |...F$. 23,180,0,|
00000090  30 2c 31 32 37 2c 32 30  34 2c 32 30 34 2c 31 32  |0,127,204,204,12|
000000a0  30 2c 30 2c 30 0d 00 50  10 de 20 44 61 74 61 25  |0,0,0..P.. Data%|
000000b0  28 36 34 30 29 0d 00 5a  12 53 63 61 6c 65 46 61  |(640)..Z.ScaleFa|
000000c0  63 74 6f 72 3d 31 36 0d  00 64 0a 53 7a 25 3d 31  |ctor=16..d.Sz%=1|
000000d0  36 0d 00 6e 08 eb 20 31  32 0d 00 78 11 e7 20 a6  |6..n.. 12..x.. .|
000000e0  28 31 35 29 3c 3e 31 33  20 8c 0d 00 82 20 f1 20  |(15)<>13 .... . |
000000f0  22 52 65 63 74 61 6e 67  75 6c 61 72 20 44 69 73  |"Rectangular Dis|
00000100  74 72 69 62 75 74 69 6f  6e 22 0d 00 8c 1f f2 49  |tribution".....I|
00000110  4e 50 55 54 20 28 22 53  61 6d 70 6c 65 20 73 69  |NPUT ("Sample si|
00000120  7a 65 3a 22 2c 53 7a 25  29 0d 00 96 1e f1 20 22  |ze:",Sz%)..... "|
00000130  50 72 65 73 73 20 61 6e  79 20 6b 65 79 20 74 6f  |Press any key to|
00000140  20 73 74 61 72 74 22 0d  00 a0 18 4d 3d 30 3a 6e  | start"....M=0:n|
00000150  3d 30 3a 54 6f 74 3d 30  3a 54 6f 74 71 3d 30 0d  |=0:Tot=0:Totq=0.|
00000160  00 aa 07 47 3d a5 0d 00  b4 05 db 0d 00 be 05 cc  |...G=...........|
00000170  0d 00 c8 1b 51 51 25 3d  8e 22 3c 4f 62 65 79 24  |....QQ%=."<Obey$|
00000180  44 69 72 3e 2e 21 44 61  74 61 22 0d 00 d2 22 e8  |Dir>.!Data"...".|
00000190  23 51 51 25 2c 53 63 61  6c 65 46 61 63 74 6f 72  |#QQ%,ScaleFactor|
000001a0  2c 54 6f 74 2c 54 6f 74  71 2c 4d 2c 6e 0d 00 dc  |,Tot,Totq,M,n...|
000001b0  1c 53 63 61 6c 65 46 61  63 74 6f 72 2b 3d 53 63  |.ScaleFactor+=Sc|
000001c0  61 6c 65 46 61 63 74 6f  72 0d 00 e6 0f e3 20 4e  |aleFactor..... N|
000001d0  3d 30 20 b8 20 36 34 30  0d 00 f0 12 e8 23 51 51  |=0 . 640.....#QQ|
000001e0  25 2c 44 61 74 61 25 28  4e 29 0d 00 fa 05 ed 0d  |%,Data%(N)......|
000001f0  01 04 09 d9 23 51 51 25  0d 01 0e 0b f2 52 65 70  |....#QQ%.....Rep|
00000200  6c 6f 74 0d 01 18 05 cd  0d 01 22 0f f5 fd a6 28  |lot......."....(|
00000210  33 30 29 3c 3e 31 33 0d  01 2c 05 f5 0d 01 36 07  |30)<>13..,....6.|
00000220  54 3d 30 0d 01 40 0f e3  20 4e 3d 31 20 b8 20 53  |T=0..@.. N=1 . S|
00000230  7a 25 0d 01 4a 26 6e 2b  3d 31 3a f1 20 8a 35 34  |z%..J&n+=1:. .54|
00000240  2c 30 29 3b 22 4e 75 6d  62 65 72 20 6f 66 20 74  |,0);"Number of t|
00000250  65 73 74 73 3a 22 3b 6e  0d 01 54 0e 74 3d b3 28  |ests:";n..T.t=.(|
00000260  36 34 30 29 2d 31 0d 01  5e 08 54 2b 3d 74 0d 01  |640)-1..^.T+=t..|
00000270  68 05 ed 0d 01 72 0b 54  3d 54 2f 53 7a 25 0d 01  |h....r.T=T/Sz%..|
00000280  7c 0d 54 3d a8 28 54 2b  2e 35 29 0d 01 86 0a 54  ||.T=.(T+.5)....T|
00000290  6f 74 2b 3d 54 0d 01 90  0d 54 6f 74 71 2b 3d 54  |ot+=T....Totq+=T|
000002a0  5e 32 0d 01 9a 38 4d 2b  3d 31 3a f1 20 8a 30 2c  |^2...8M+=1:. .0,|
000002b0  30 29 3b 22 6e 3a 22 3b  4d 3b 8a 33 30 2c 30 29  |0);"n:";M;.30,0)|
000002c0  3b 22 53 61 6d 70 6c 65  20 41 76 65 72 61 67 65  |;"Sample Average|
000002d0  3a 22 3b 54 3b 22 20 20  20 22 0d 01 a4 37 f1 8a  |:";T;"   "...7..|
000002e0  30 2c 31 29 3b 22 b5 3a  22 3b 54 6f 74 2f 4d 3b  |0,1);".:";Tot/M;|
000002f0  8a 36 38 2c 31 29 3b 22  b4 3a 22 3b b6 28 28 54  |.68,1);".:";.((T|
00000300  6f 74 71 2f 4d 29 2d 28  54 6f 74 2f 4d 29 5e 32  |otq/M)-(Tot/M)^2|
00000310  29 0d 01 ae 0f 44 61 74  61 25 28 54 29 2b 3d 31  |)....Data%(T)+=1|
00000320  0d 01 b8 59 e7 20 44 61  74 61 25 28 54 29 3e 28  |...Y. Data%(T)>(|
00000330  32 35 36 2f 53 63 61 6c  65 46 61 63 74 6f 72 29  |256/ScaleFactor)|
00000340  20 8c 20 f2 52 65 70 6c  6f 74 20 8b 20 e6 20 34  | . .Replot . . 4|
00000350  3a 86 20 28 54 2a 32 29  2c 30 2c 28 54 2a 32 29  |:. (T*2),0,(T*2)|
00000360  2c 28 44 61 74 61 25 28  54 29 2a 34 2a 53 63 61  |,(Data%(T)*4*Sca|
00000370  6c 65 46 61 63 74 6f 72  29 0d 01 c2 0c 51 51 51  |leFactor)....QQQ|
00000380  3d a6 28 30 29 0d 01 cc  15 fd 20 51 51 51 3d 33  |=.(0)..... QQQ=3|
00000390  32 20 84 20 51 51 51 3d  31 33 0d 01 d6 18 e7 20  |2 . QQQ=13..... |
000003a0  51 51 51 3d 31 33 20 8c  20 53 3d b9 20 8b 20 53  |QQQ=13 . S=. . S|
000003b0  3d a3 0d 01 e0 1d f5 fd  20 bf 28 35 30 29 3c 3e  |=....... .(50)<>|
000003c0  22 20 22 20 80 20 a6 28  35 30 29 3c 3e 31 33 0d  |" " . .(50)<>13.|
000003d0  01 ea 09 e7 20 53 20 8c  0d 01 f4 1b 51 51 25 3d  |.... S .....QQ%=|
000003e0  ae 22 3c 4f 62 65 79 24  44 69 72 3e 2e 21 44 61  |."<Obey$Dir>.!Da|
000003f0  74 61 22 0d 01 fe 22 f1  23 51 51 25 2c 53 63 61  |ta"...".#QQ%,Sca|
00000400  6c 65 46 61 63 74 6f 72  2c 54 6f 74 2c 54 6f 74  |leFactor,Tot,Tot|
00000410  71 2c 4d 2c 6e 0d 02 08  0f e3 20 4e 3d 30 20 b8  |q,M,n..... N=0 .|
00000420  20 36 34 30 0d 02 12 12  f1 23 51 51 25 2c 44 61  | 640.....#QQ%,Da|
00000430  74 61 25 28 4e 29 0d 02  1c 05 ed 0d 02 26 09 d9  |ta%(N).......&..|
00000440  23 51 51 25 0d 02 30 05  cd 0d 02 3a 18 56 3d 28  |#QQ%..0....:.V=(|
00000450  54 6f 74 71 2f 4d 29 2d  28 54 6f 74 2f 4d 29 5e  |Totq/M)-(Tot/M)^|
00000460  32 0d 02 44 0b 45 3d 54  6f 74 2f 4d 0d 02 4e 12  |2..D.E=Tot/M..N.|
00000470  51 3d 31 2f 28 b6 28 32  2a af 2a 56 29 29 0d 02  |Q=1/(.(2*.*V))..|
00000480  58 07 e6 20 33 0d 02 62  07 5a 3d 4d 0d 02 6c 09  |X.. 3..b.Z=M..l.|
00000490  ec 20 30 2c 30 0d 02 76  0f e3 20 4e 3d 30 20 b8  |. 0,0..v.. N=0 .|
000004a0  20 36 34 30 0d 02 80 35  df 20 4e 2a 32 2c 28 51  | 640...5. N*2,(Q|
000004b0  2a a1 28 28 2d 31 2a 28  28 4e 2d 45 29 5e 32 29  |*.((-1*((N-E)^2)|
000004c0  29 2f 28 32 2a 56 29 29  29 2a 5a 2a 53 63 61 6c  |)/(2*V)))*Z*Scal|
000004d0  65 46 61 63 74 6f 72 2a  34 0d 02 8a 05 ed 0d 02  |eFactor*4.......|
000004e0  94 0f f5 fd 20 bf 28 30  29 3d 22 20 22 0d 02 9e  |.... .(0)=" "...|
000004f0  11 f5 fd 20 bf 28 35 30  29 3c 3e 22 20 22 0d 02  |... .(50)<>" "..|
00000500  a8 05 e0 0d 02 b2 05 3a  0d 02 bc 13 dd f2 49 4e  |.......:......IN|
00000510  50 55 54 28 53 24 2c f8  20 56 29 0d 02 c6 07 ea  |PUT(S$,. V).....|
00000520  20 85 0d 02 d0 08 f1 53  24 3b 0d 02 da 05 87 0d  | ......S$;......|
00000530  02 e4 0b ee 20 85 20 ea  20 85 0d 02 ee 05 f5 0d  |.... . . .......|
00000540  02 f8 23 e7 20 56 3c 3e  30 20 8c 20 f1 20 c3 28  |..#. V<>0 . . .(|
00000550  56 29 3b 22 20 20 20 20  20 20 20 20 20 20 20 20  |V);"            |
00000560  22 3b 0d 03 02 0c e7 20  56 3c 3e 30 20 8c 0d 03  |";..... V<>0 ...|
00000570  0c 23 e3 20 4e 3d 31 20  b8 20 a9 20 28 c3 28 56  |.#. N=1 . . (.(V|
00000580  29 2b 22 20 20 20 20 20  20 20 20 20 20 20 20 22  |)+"            "|
00000590  29 0d 03 16 0b f1 20 bd  28 38 29 3b 0d 03 20 05  |)..... .(8);.. .|
000005a0  ed 0d 03 2a 05 cd 0d 03  34 16 e7 20 56 3d 30 20  |...*....4.. V=0 |
000005b0  8c 20 f1 22 20 22 3b bd  28 38 29 3b 0d 03 3e 08  |. ." ";.(8);..>.|
000005c0  47 24 3d be 0d 03 48 1b  e7 20 97 28 47 24 29 3d  |G$=...H.. .(G$)=|
000005d0  31 33 20 8c 20 fd b9 3a  f7 20 85 3a f1 3a e1 0d  |13 . ..:. .:.:..|
000005e0  03 52 35 e7 20 28 47 24  3c 22 30 22 20 84 20 47  |.R5. (G$<"0" . G|
000005f0  24 3e 22 39 22 29 20 80  20 97 28 47 24 29 3c 3e  |$>"9") . .(G$)<>|
00000600  31 32 37 20 80 20 97 28  47 24 29 3c 3e 38 20 8c  |127 . .(G$)<>8 .|
00000610  20 fd 20 a3 0d 03 5c 2e  e7 20 97 28 47 24 29 3c  | . ...\.. .(G$)<|
00000620  3e 31 32 37 20 80 20 97  28 47 24 29 3c 3e 38 20  |>127 . .(G$)<>8 |
00000630  8c 20 56 3d 56 2a 31 30  2b a0 28 47 24 29 3a fd  |. V=V*10+.(G$):.|
00000640  20 a3 0d 03 66 0c 56 3d  56 20 81 20 31 30 0d 03  | ...f.V=V . 10..|
00000650  70 07 fd 20 a3 0d 03 7a  05 3a 0d 03 84 0c dd f2  |p.. ...z.:......|
00000660  52 65 70 6c 6f 74 0d 03  8e 1d 53 63 61 6c 65 46  |Replot....ScaleF|
00000670  61 63 74 6f 72 3d 53 63  61 6c 65 46 61 63 74 6f  |actor=ScaleFacto|
00000680  72 2f 32 0d 03 98 05 db  0d 03 a2 0f e3 20 54 3d  |r/2.......... T=|
00000690  30 20 b8 20 36 34 30 0d  03 ac 30 e6 20 34 3a 86  |0 . 640...0. 4:.|
000006a0  20 28 54 2a 32 29 2c 30  2c 28 54 2a 32 29 2c 28  | (T*2),0,(T*2),(|
000006b0  44 61 74 61 25 28 54 29  2a 34 2a 53 63 61 6c 65  |Data%(T)*4*Scale|
000006c0  46 61 63 74 6f 72 29 0d  03 b6 05 ed 0d 03 c0 05  |Factor).........|
000006d0  e1 0d ff                                          |...|
000006d3

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