# script_quad_glgr
# Autor
# 02/05/2020
# Quadratura de Gauss-Laguerre para calcular integral do tipo
#   \int_{0}^{\infty}f(x)e^{-x}dx
rm(list = ls())

# pacote
library(statmod)

# função f(x)
f <- function(x){
  out <- x+1
  return(out)
}

# pontos e pesos da quadratura de Gauss-Laguerre
n <- 5
lgr <- gauss.quad(n,'laguerre')
# pontos
x <- lgr$nodes
# pesos
w <- lgr$weights

# fórmula da quadratura
soma <- 0
for(i in 1:n){
  soma <- soma+(w[i]*f(x[i]))
}

# Aproximação da integral
soma

# ou 
cat('O valor aproximado da integral é:', soma, '\n')