# Program 23
# This program initiates the use of RSquadrature.uniform, and
# obtains the resulting weights and abscissae.  See Example 7.3.1

p <- 3
Nr <- 3
Nq <- 4

#lower limits for parameters in first column, upper limits in second column
limits <- matrix(c(-1,-3,-3,1,3,3),p,2)

set.seed(12345)
out <- RSquadrature.uniform(p, limits, Nr, Nq)
wts <- out$w
abscissae <- out$a
cbind(wts,abscissae)