在R软件中,我们如何计算定积分呢?
下面讲解两个例子。
R中的积分命令是:
"integrate"
例1:对常见密度函数积分:
integrate(dnorm, -1.96, 1.96)
integrate(dnorm, -Inf, Inf)
例2:对自定义函数积分:
## a slowly-convergent integral
integrand <- function(x) {1/((x+1)*sqrt(x))}
integrate(integrand, lower = 0, upper = Inf)