R语言 蒙特卡洛方法（Monte Carlo）

蒙特卡洛方法是指通过产生大量的随机数，来用于数值统计计算以获得问题的近似解。一个经典的例子是用蒙特卡洛模拟计算圆周率π，如下图所示，在正方形内产生足够多的随机点，这些点落在圆内的概率就是圆面积和正方形面积的比值，所以落在圆内的点的个数比上所有的点的个数（在正方形内的点）就等于落在圆内的概率。据此即可计算出π = 4*圆面积/正方形面积，下面用代码说明。

> calpi <- function(n){
   
+   x <- matrix(runif(