r语言用重要抽样蒙特卡罗积分估计e^-x/(1+x^2)在[0,1]上的积分并计算方差的代码，不写注释

# 定义被积函数
f <- function(x) exp(-x)/(1 + x^2)

# 定义重要抽样函数
w <- function(x) 1/(2*sqrt(x))

# 定义抽样个数
n <- 10000

# 生成服从w分布的随机数
x <- runif(n)
y <- w(x)

# 计算被积函数在随机数处的值
fx <- f(x)

# 计算积分估计值和方差
integral <- sum(fx*y)/n
variance <- sum((fx*y)^2)/n - integral^2

# 输出结果
cat("Integral estimate: ", integral, "\n")
cat("Variance: ", variance, "\n")