simp {broman} | R Documentation |
Perform numerical integration by Simpson's rule or the trapezoidal rule.
simp(f, a, b, tol = 0.00000001, max.step = 1000, ...)
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, |
Iterately doubles the number of grid points for the numerical
integral, stopping when the integral decreases by less than
tol
.
The integral of f
from a
to b
.
f <- function(x) x*x*(1-x)*sin(x*x)
I1 <- trap(f,0,2)
I2 <- simp(f,0,2)