Pythonを使って円周率を求める #モンテカルロ法 - Qiita

初めに この記事では ラムダさんの動画を見てそのまま、深夜テンションで書いたモンテカルロ法を用いて円周率を求めるプログラムを解説する記事です。深夜テンションなのであまり期待しないで下さい！ コードが酷い import matplotlib.pyplot as plt import random Xlist=[] Ylist=[] Number_list=[] pi=0 pi_list=[] for i in range(15000): Xlist.append(random.uniform(0,1)) Ylist.append(random.uniform(0,1)) if (i+1)%10 == 0: Number_list.append(i+1) pi_list.append((pi/(i+1))*4) if Xlist[i]**2+Ylist[i]**2 <= 1 and (i+1)