Helper for Efficient Integration
integrateVector(integrand, upper, ...)(function)
to be integrated.
(numeric)
upper limits of integration.
additional arguments to be passed to integrand.
This function returns for each upper limit the estimates of the integral.
integrand <- function(x) x^2
upper <- c(0, 1, 0.4, 2, 5, 2, 0.3, 0.4, 1)
integrateVector(integrand, upper = upper)
#> [1] 0.00000000 0.33333333 0.02133333 2.66666667 41.66666667 2.66666667
#> [7] 0.00900000 0.02133333 0.33333333