fda.usc (version 0.9.4)

int.simpson: Simpson integration

Description

Computes the integral of y with respect to x using simpson or trapezoid rule integration.

Usage

int.simpson(fdataobj,equi=TRUE,method="CSR")
int.simpson2(x, y, equi = TRUE, method = "CSR")

Arguments

fdataobj

fdata objtect.

x

Sorted vector of x-axis values: argvals.

y

Vector of y-axis values.

equi

=TRUE, the observed points on each curve are equally spaced (by default).

method

Method for numerical integration, see details.

Details

If method="CSR", composite Simpson's rule integration is used. If method="ESR", extended Simpson's rule integration is used. If method="TRAPZ", Trapezoid rule integration is used. If x is a fdata class object it computes the integral of x$data with Respect to x$argvals.

See Also

See also integrate.

Examples

x<-seq(0,2*pi,length=1001)
fx<-fdata(sin(x)/sqrt(pi),x)
fx0<-fdata(rep(0,length(x)),x)
int.simpson(fx0$data,fx$data)
int.simpson(fx)
int.simpson2(fx$argvals,fx$data^2)
int.simpson2(fx$argvals,fx$data^2,method="TRAPZ")
int.simpson2(fx$argvals,fx$data^2,method="ESR")