モンテカルロ法を使った円周率を求める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"となる。