R/integrate.q.R

Defines functions integrate.q

# Quick integration:
integrate.q <-
function(x, y){

  n <- length(x)
  
# Shift to find lower and upper vectors: 1-(n - 1) and 2-n
  lx <- x[-n]; ux <- x[-1]
  ly <- y[-n]; uy <- y[-1]
 
  sum( (ux - lx) * (ly + uy) / 2 )

} # END integrate.q

Try the catIrt package in your browser

Any scripts or data that you put into this service are public.

catIrt documentation built on May 28, 2022, 1:09 a.m.