R

时间:2018-03-19 08:19:20

标签: r algorithm simulation montecarlo

我想使用命中和未命中方法计算以下积分。 I=∫x^3dx lower= 0upper =1

我知道如何解决它,但我无法在R中找到正确的代码来计算它并生成 - 例如100000随机 - 然后将它们绘制成如下:

https://i.stack.imgur.com/r0CmK.png

谢谢。

1 个答案:

答案 0 :(得分:2)

1。从所需长度的均匀分布生成2个向量

$( document ).ready(function() {
   $('#btnUpdate').attr('disabled', true);  
});

2. 积分的近似值是(x,y)点低于要整合的函数的情况数:

l = 10000
x = runif(l)
y = runif(l)

3. 对于绘图,您只需绘制点,根据曲线的上方或下方更改颜色,并使用sum(y<x^3)/l 添加函数:

curve()