powered by
This function calculates nodes and weights for Gaussian quadrature. See help("gauss.quad") from package statmod.
help("gauss.quad")
statmod