2017年6月11日日曜日

Excel VBAで数学 7

 Excel VBAを使った簡単な数学プログラミングの続きです。

 今回は円周率の計算です。
 値は知っていても、実際に計算したことがある人は、多くないのではないでしょうか。
 今回は、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 件のコメント:

コメントを投稿