Returns two vectors which hold the needed quadrature nodes and weights.
Arguments
n
number of nodes (and weights).
type
quadrature type -- currently Gauss-Legendre, "GL", and Radau, "Ra", are supported.
x1
lower limit of the integration interval.
x2
upper limit of the integration interval.
Author
Sven Knoth
Details
A more detailed description will follow soon. The algorithm for the Gauss-Legendre quadrature was delivered by
Knut Petras to me, while the one for the Radau quadrature was taken from John Burkardt.
References
H. Brass and K. Petras (2011),
Quadrature Theory. The Theory of Numerical Integration on a Compact Interval,
Mathematical Surveys and Monographs, American Mathematical Society.