Simpson integration

Computes the integral of fdataobj$data with respect to fdataobj$argvals using simpson or trapezoid rule integration.

Usage

int.simpson(fdataobj, method = NULL)

int.simpson2(x, y, equi = TRUE, method = NULL)

Arguments

fdataobj

fdata objtect.

method

Method for numerical integration, see details.

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).

Details

Posible values for method are:

• "TRAPZ": Trapezoid rule integration.

• "CSR": Composite Simpson's rule integration.

• "ESR": Extended Simpson's rule integration.

If method=NULL (default), the value of par.fda.usc$int.method is used.

See also

See also integrate.

Author

Manuel Febrero-Bande, Manuel Oviedo de la Fuente manuel.oviedo@udc.es

Examples

if (FALSE) { # \dontrun{
x<-seq(0,2*pi,length=1001)
fx<-fdata(sin(x)/sqrt(pi),x)
fx0<-fdata(rep(0,length(x)),x)
int.simpson(fx0)
int.simpson(fx)
} # }