broman (version 0.72-4)

simp: Numerical integration

Description

Perform numerical integration by Simpson's rule or the trapezoidal rule.

Usage

simp(f, a, b, tol = 0.00000001, max.step = 1000, ...)

Arguments

f

The integrand; must be a vectorized function.

a

Lower limit of integration.

b

Upper limit of integration.

tol

Tolerance for choosing the number of grid points.

max.step

Log base 2 of the total number of grid points.

...

Other arguments passed to the integrand, f.

Value

The integral of f from a to b.

Details

Iterately doubles the number of grid points for the numerical integral, stopping when the integral decreases by less than tol.

See Also

stats::integrate()

Examples

# NOT RUN {
f <- function(x) x*x*(1-x)*sin(x*x)
I1 <- trap(f,0,2)
I2 <- simp(f,0,2)

# }