我希望(1/y)*(2/(1+(log(y))^2))从0到1 进行整合.Wolfram alpha告诉我这应该是pi.但是当我在R中进行蒙特卡罗整合时,我在尝试超过10次后仍然获得3.00和2.99.这就是我所做的:
y=runif(10^6)
f=(1/y)*(2/(1+(log(y))^2))
mean(f)
我将确切的函数复制到wolfram alpha中以检查积分应该是pi
我试图通过检查它的平均值并绘制历史图来检查我的y是否正确分布,它似乎没问题.我的电脑有问题吗？
编辑:也许其他人可以复制我的代码并自己运行,以确认它不是我的电脑表现.