Also called Radau quadrature (Chandrasekhar 1960). A Gaussian quadrature with weighting function
in which the endpoints of the interval
are included in a total of
abscissas, giving
free abscissas. Abscissas
are symmetrical about the origin, and the general formula
is
|
(1)
|
The free abscissas for
, ...,
are the roots of the polynomial
, where
is a Legendre polynomial.
The weights of the free abscissas are
|
(2)
| |||
|
(3)
|
and of the endpoints are
|
(4)
|
The error term is given by
|
(5)
|
for .
Beyer (1987) gives a table of parameters up to
and Chandrasekhar (1960) up to
(although Chandrasekhar's
for
is incorrect).
| 3 | 0 | 0.00000 | 1.333333 | |
| 0.333333 | ||||
| 4 | 0.833333 | |||
| 0.166667 | ||||
| 5 | 0 | 0.000000 | 0.711111 | |
| 0.544444 | ||||
| 0.100000 | ||||
| 6 | 0.554858 | |||
| 0.378475 | ||||
| 0.066667 |