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
swnydick/catIrt documentation built on June 13, 2022, 9:30 p.m.