Home » Archimedes archive » Acorn User » AU 1997-11 B.adf » Regulars » StarInfo/vanMourik/PI
StarInfo/vanMourik/PI
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-11 B.adf » Regulars |
| Filename: | StarInfo/vanMourik/PI |
| Read OK: | ✔ |
| File size: | 02C0 bytes |
| Load address: | 0000 |
| Exec address: | 0000 |
File contents
10REM >PI 20REM JEAN VAN MOURIK 30PRINT"ITERATED CALCULATION OF PI WITH A VERY QUICK CONVERCENCE." 40PRINT"THIS SERIES FORMULA IS BY DR. PETER BORWEIN. (SCIENCE NEWS 28 OCT. 1995)":PRINT:PRINT"S=1:I=0:SUM=0":PRINT"REPEAT:V=8*I:T=S*(4/(V+1)-2/(V+4)-1/(V+5)-1/(V+6)):SUM=SUM+T:I=I+1:S=S/16:UNTIL0" 50S=1:I=0:SM=0:PRINT"PI = ";PI:PRINT 60 70REPEAT 80V=8*I:T=S*(4/(V+1)-2/(V+4)-1/(V+5)-1/(V+6)):SM=SM+T:I+=1:S=S/16 90PRINT"I= ";I;" SUM= ";SM 100REPEAT:MOUSE A,B,C:UNTILC:REPEAT:MOUSE A,B,C:UNTILC=0 110UNTILI=9 120 130PRINT:PRINT"A SIMPLER FORMULA AND LESS QUICK TO CONVERGE;":PRINT"I=0:SUM=0:C%=4:N%=1":PRINT"REPEAT:SUM=SUM+C%/N%:N%=N%+2:C%=-C%:UNTIL0":PRINT 140I=0:SM=0:C%=4:N%=1 145 150REPEAT 160D=C%/N%:SM+=D 170PRINT"N%= ";N%;" ";D;" + SUM= ";SM 175N%+=2:C%=-C%:REPEAT:MOUSE A,B,C:UNTILC 180UNTIL0
� >PI � JEAN VAN MOURIK @�"ITERATED CALCULATION OF PI WITH A VERY QUICK CONVERCENCE." (��"THIS SERIES FORMULA IS BY DR. PETER BORWEIN. (SCIENCE NEWS 28 OCT. 1995)":�:�"S=1:I=0:SUM=0":�"REPEAT:V=8*I:T=S*(4/(V+1)-2/(V+4)-1/(V+5)-1/(V+6)):SUM=SUM+T:I=I+1:S=S/16:UNTIL0" 2S=1:I=0:SM=0:�"PI = ";�:� < F� PCV=8*I:T=S*(4/(V+1)-2/(V+4)-1/(V+5)-1/(V+6)):SM=SM+T:I+=1:S=S/16 Z�"I= ";I;" SUM= ";SM d!�:ȗ A,B,C:�C:�:ȗ A,B,C:�C=0 n�I=9 x �}�:�"A SIMPLER FORMULA AND LESS QUICK TO CONVERGE;":�"I=0:SUM=0:C%=4:N%=1":�"REPEAT:SUM=SUM+C%/N%:N%=N%+2:C%=-C%:UNTIL0":� �I=0:SM=0:C%=4:N%=1 � �� �D=C%/N%:SM+=D �&�"N%= ";N%;" ";D;" + SUM= ";SM �N%+=2:C%=-C%:�:ȗ A,B,C:�C ��0 �
00000000 0d 00 0a 09 f4 20 3e 50 49 0d 00 14 15 f4 20 4a |..... >PI..... J| 00000010 45 41 4e 20 56 41 4e 20 4d 4f 55 52 49 4b 0d 00 |EAN VAN MOURIK..| 00000020 1e 40 f1 22 49 54 45 52 41 54 45 44 20 43 41 4c |.@."ITERATED CAL| 00000030 43 55 4c 41 54 49 4f 4e 20 4f 46 20 50 49 20 57 |CULATION OF PI W| 00000040 49 54 48 20 41 20 56 45 52 59 20 51 55 49 43 4b |ITH A VERY QUICK| 00000050 20 43 4f 4e 56 45 52 43 45 4e 43 45 2e 22 0d 00 | CONVERCENCE."..| 00000060 28 b7 f1 22 54 48 49 53 20 53 45 52 49 45 53 20 |(.."THIS SERIES | 00000070 46 4f 52 4d 55 4c 41 20 49 53 20 42 59 20 44 52 |FORMULA IS BY DR| 00000080 2e 20 50 45 54 45 52 20 42 4f 52 57 45 49 4e 2e |. PETER BORWEIN.| 00000090 20 20 28 53 43 49 45 4e 43 45 20 4e 45 57 53 20 | (SCIENCE NEWS | 000000a0 32 38 20 4f 43 54 2e 20 31 39 39 35 29 22 3a f1 |28 OCT. 1995)":.| 000000b0 3a f1 22 53 3d 31 3a 49 3d 30 3a 53 55 4d 3d 30 |:."S=1:I=0:SUM=0| 000000c0 22 3a f1 22 52 45 50 45 41 54 3a 56 3d 38 2a 49 |":."REPEAT:V=8*I| 000000d0 3a 54 3d 53 2a 28 34 2f 28 56 2b 31 29 2d 32 2f |:T=S*(4/(V+1)-2/| 000000e0 28 56 2b 34 29 2d 31 2f 28 56 2b 35 29 2d 31 2f |(V+4)-1/(V+5)-1/| 000000f0 28 56 2b 36 29 29 3a 53 55 4d 3d 53 55 4d 2b 54 |(V+6)):SUM=SUM+T| 00000100 3a 49 3d 49 2b 31 3a 53 3d 53 2f 31 36 3a 55 4e |:I=I+1:S=S/16:UN| 00000110 54 49 4c 30 22 0d 00 32 1d 53 3d 31 3a 49 3d 30 |TIL0"..2.S=1:I=0| 00000120 3a 53 4d 3d 30 3a f1 22 50 49 20 3d 20 22 3b af |:SM=0:."PI = ";.| 00000130 3a f1 0d 00 3c 05 20 0d 00 46 05 f5 0d 00 50 43 |:...<. ..F....PC| 00000140 56 3d 38 2a 49 3a 54 3d 53 2a 28 34 2f 28 56 2b |V=8*I:T=S*(4/(V+| 00000150 31 29 2d 32 2f 28 56 2b 34 29 2d 31 2f 28 56 2b |1)-2/(V+4)-1/(V+| 00000160 35 29 2d 31 2f 28 56 2b 36 29 29 3a 53 4d 3d 53 |5)-1/(V+6)):SM=S| 00000170 4d 2b 54 3a 49 2b 3d 31 3a 53 3d 53 2f 31 36 0d |M+T:I+=1:S=S/16.| 00000180 00 5a 19 f1 22 49 3d 20 22 3b 49 3b 22 20 20 53 |.Z.."I= ";I;" S| 00000190 55 4d 3d 20 22 3b 53 4d 0d 00 64 21 f5 3a c8 97 |UM= ";SM..d!.:..| 000001a0 20 41 2c 42 2c 43 3a fd 43 3a f5 3a c8 97 20 41 | A,B,C:.C:.:.. A| 000001b0 2c 42 2c 43 3a fd 43 3d 30 0d 00 6e 08 fd 49 3d |,B,C:.C=0..n..I=| 000001c0 39 0d 00 78 05 20 0d 00 82 7d f1 3a f1 22 41 20 |9..x. ...}.:."A | 000001d0 53 49 4d 50 4c 45 52 20 46 4f 52 4d 55 4c 41 20 |SIMPLER FORMULA | 000001e0 41 4e 44 20 4c 45 53 53 20 51 55 49 43 4b 20 54 |AND LESS QUICK T| 000001f0 4f 20 43 4f 4e 56 45 52 47 45 3b 22 3a f1 22 49 |O CONVERGE;":."I| 00000200 3d 30 3a 53 55 4d 3d 30 3a 43 25 3d 34 3a 4e 25 |=0:SUM=0:C%=4:N%| 00000210 3d 31 22 3a f1 22 52 45 50 45 41 54 3a 53 55 4d |=1":."REPEAT:SUM| 00000220 3d 53 55 4d 2b 43 25 2f 4e 25 3a 4e 25 3d 4e 25 |=SUM+C%/N%:N%=N%| 00000230 2b 32 3a 43 25 3d 2d 43 25 3a 55 4e 54 49 4c 30 |+2:C%=-C%:UNTIL0| 00000240 22 3a f1 0d 00 8c 16 49 3d 30 3a 53 4d 3d 30 3a |":.....I=0:SM=0:| 00000250 43 25 3d 34 3a 4e 25 3d 31 0d 00 91 05 20 0d 00 |C%=4:N%=1.... ..| 00000260 96 05 f5 0d 00 a0 11 44 3d 43 25 2f 4e 25 3a 53 |.......D=C%/N%:S| 00000270 4d 2b 3d 44 0d 00 aa 26 f1 22 4e 25 3d 20 22 3b |M+=D...&."N%= ";| 00000280 4e 25 3b 22 20 20 20 20 20 22 3b 44 3b 22 20 2b |N%;" ";D;" +| 00000290 20 53 55 4d 3d 20 22 3b 53 4d 0d 00 af 1e 4e 25 | SUM= ";SM....N%| 000002a0 2b 3d 32 3a 43 25 3d 2d 43 25 3a f5 3a c8 97 20 |+=2:C%=-C%:.:.. | 000002b0 41 2c 42 2c 43 3a fd 43 0d 00 b4 06 fd 30 0d ff |A,B,C:.C.....0..| 000002c0
.