Estimate the integral of one-dimensional function using the trapezoidal rule.
vectors of x and y values. In the normal case, these
vectors are both supplied, and of equal length. There are also two
special cases. First, if
Flag indicating the desired return value (see “Value”).
Optional numbers indicating the range of the integration.
These values may be used to restrict the range of integration, or to
extend it; in either case,
type="A" (the default), a single value is returned,
containing the estimate of the integral of
type="dA", a numeric vector of the same length as
x, of which
the first element is zer0, the second element is the integral between
x, etc. If
type="cA", the result is the
cumulative sum (as in
cumsum) of the values that would be
type="dA". See “Examples”.
There is no handling of
1 2 3 4 5 6 7 8 9 10
x <- seq(0, 1, length.out=10) # try larger length.out to see if area approaches 2 y <- 2*x + 3*x^2 A <- integrateTrapezoid(x, y) dA <- integrateTrapezoid(x, y, "dA") cA <- integrateTrapezoid(x, y, "cA") print(A) print(sum(dA)) print(tail(cA, 1)) print(integrateTrapezoid(diff(x[1:2]), y)) print(integrateTrapezoid(y))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.