モンテカルロ法によるπの推定 - Qiita

モンテカルロ法によってπを推定する。 方法 正方形の中に、x,y座標ともに一様な乱数を点としてプロットする。そのうち、円の中にある点の数を数え上げることで、円の中に存在する点の確率が計算できる。 この確率から面積を求めることができる。 具体的には、x,y座標ともに(-1,-1)と(1,1)となる正方形とそこに密接する円を想定する。 import matplotlib.pyplot as plt import matplotlib.patches as patches import numpy as np fig,ax = plt.subplots(figsize=(5, 5)) rect = patches.Rectangle((-1,-1),2,2,linewidth=1, fill=False, color='blue') ax.add_patch(rect) ax.add_artis