蒙特卡落积分求圆周率π

一、方法

        假设我们想要估计圆周率π。我们知道半径为r的圆的面积为πr²，同时也等于下面定义的积分I。因此，圆周率π=I/(r²)。半径是已知量，而关键就在于如何求得圆的面积I。让我们通过蒙特卡洛积分求得I的近似值。在下面的公式中，指示函数用于判断一个二维点是否在圆内，如果二维点在圆内则指示函数的值为1否则为0.

        为了更加直观地理解这个方法，请看下面的图2.19。其中上面公式(2.102)中的符号S为下图中的所有点的个数。

二、参考 

        Machine Learning: A Probabilistic Perspective (The MIT Press)