R中的“命中与未命中”积分方法

我想使用命中和未命中方法计算以下积分。

I=∫x^3dx

带

下=0

和

上=1

我知道如何解决它，但我无法在

R

中找到正确的代码来计算它并生成-例如100000 random-然后像这样绘制它们：

---

谢谢。

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

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

2.积分的近似值是（x，y）点低于要积分的函数的情况数：

sum(y<x^3)/l

---

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

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

2.积分的近似值是（x，y）点低于要积分的函数的情况数：

sum(y<x^3)/l

---

你知道我怎么用对偶变量计算同一个函数吗？不知道，但是你知道我怎么用对偶变量计算同一个函数吗？不知道，但是你知道