数值积分的算法

https://www.r-bloggers.com/lang/chinese/841

1、用integrate函数直接求出数值积分
integrate(function(x) sin(x),0,pi)

2、蒙特卡罗仿真
x <- runif(100000,min=0,max=pi)
y <- runif(100000)
pi*sum(y<sin(x))/100000

解释：x的范围是从0到π，y的范围是从sinπ/2到sin0（runif默认即min=0，max=1），求解过程即：sum（y<sin(x)）*（π-0）/10000即高乘宽除以试验次数。

3、