r语言蒙特卡洛求积分

使用 Monte Carlo 方法求解积分的基本思路是随机生成一组样本点，并根据这些样本点的函数值来计算积分的近似值。在 R 语言中，可以通过以下代码来实现：

# 定义被积函数
f <- function(x) {
  return(x^2)
}

# 定义积分范围
a <- 0
b <- 1

# 生成样本点
n <- 10000
x <- runif(n, a, b)
y <- f(x)

# 计算积分的近似值
integral <- (b-a) * mean(y)
print(integral)

以上代码中，我们定义了被积函数 f(x)，并指定了积分范围 [a, b]。然后，我们随机生成 n 个样本点，并根据这些样本点的函数值来计算积分的近似值。