■円周率の計算(その40)

  π=2{AGM(1,1/√2)}^2/(1−Σ2^ncn^2)

  AGM(1,1/√2)=0.8472130847

より,ガウスは

  AGM(1,1/√2)=π/2・1/K(1/√2)

であることを見抜き,その証明を与えた.

===================================

【1】ガウスの公式による計算

100 PFILE$="scrn:"

110 OPEN PFILE$ FOR OUTPUT AS #1

120 A=1:X=1

130 B=1/SQR(2)

140 C=1/4

150 I=1

160 '

170 WHILE I<=20

180 Y=A

190 A=(A+B)/2

200 B=SQR(B*Y)

210 C=C-X*(A-Y)*(A-Y)

220 X=X*2

230 PRINT #1,I,(A+B)*(A+B)/C/4

240 I=I+1

250 WEND

260 CLOSE #1

270 END

1 3.14058

2 3.14159

3 3.14159

4 3.14159

5 3.14159

6 3.14159

7 3.14159

8 3.14159

9 3.14159

10 3.14159

11 3.14159

12 3.14159

13 3.14159

14 3.14159

15 3.14159

16 3.14159

17 3.14159

18 3.14159

19 3.14159

20 3.14159

===================================