x <- 1:100 y <- dexp(1:100) plot(x, y) #覦覯2 library("pracma") trapz(x,y) #覦覯1 library("flux") auc(x,y)