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