かまりくんの遊び場

電子工作とか

モンテカルロ法を使った円周率を求めるpythonプログラム

import random
gco = 0
#nは求めたい桁数
n = input("求めたい桁数を入力しようね!")
for i in range(pow(10,n)):
    x = random.random()
    y = random.random()
    a = pow(x,2)+pow(y,2)
    i = i - 1
    if(a<=1):
        gco = gco + 1

pi = 4*gco/(pow(10,n))
print(pi)

さて、このプログラムの解説に入ろう。このプログラムで一番大事なのはファイル名である。まずpythonでプログラムを書いているので拡張子は".py"となり、そして円周率はπ(pi)である。つまりファイル名は"pi.py"となる。