今回は円周率の計算です。
値は知っていても、実際に計算したことがある人は、多くないのではないでしょうか。
今回は、Machinの公式というのを使っています。計算機がない時代に手計算をしていた人々のことを思うと、感慨深いです。
ソースコードは、ご自由にご利用ください。ただし、趣味のプログラムなので、保証はありません。
Option Explicit
'piの計算 (Machinの公式)
Public Sub Macro()
Dim i As Long
Dim m As Long
Dim p As Double
Dim temp As Double
m = 100
temp = 16 / 5
For i = 0 To 3 * m + 2
p = p + temp / (2 * i + 1)
temp = -temp / 5 / 5
Next i
temp = 4 / 239
For i = 0 To m
p = p - temp / (2 * i + 1)
temp = -temp / 239 / 239
Next i
Sheet1.Cells(1, 1) = p
End Sub
0 件のコメント:
コメントを投稿