powered by
gaussint(f, x, w) gauss.legendre(f, m = 5) gauss.laguerre(f, m = 5) gauss.hermite(f, m = 5)
gaussint
adaptint
giniquintile
mcint
midpt
revolution-solid
romberg
simp38
simp
trap
w = c(1, 1) x = c(-1 / sqrt(3), 1 / sqrt(3)) f <- function(x) { x^3 + x + 1 } gaussint(f, x, w)
Run the code above in your browser using DataLab