determines the (n+1)/2 nonnegative points x(i) and
the corresponding weights w(i) of the n-point
gauss-legendre integration rule, normalized to the
interval [-1,1]. the x(i) appear in descending order.
Note
This routine is from 'methods of numerical integration', P.J. Davis and P. Rabinowitz, page 369.
.. Scalar Arguments ..
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer | :: | N | ||||
double precision | :: | X(*) | ||||
double precision | :: | W(*) |