************80
! F1_ABSCISSAS_AB computes Fejer type 1 abscissas for the interval [A,B].
Licensing:
This code is distributed under the GNU LGPL license.
Modified:
29 December 2007
Author:
John Burkardt
Reference:
Philip Davis, Philip Rabinowitz,
Methods of Numerical Integration,
Second Edition,
Dover, 2007,
ISBN: 0486453391,
LC: QA299.3.D28.
Walter Gautschi,
Numerical Quadrature in the Presence of a Singularity,
SIAM Journal on Numerical Analysis,
Volume 4, Number 3, 1967, pages 357-362.
Joerg Waldvogel,
Fast Construction of the Fejer and Clenshaw-Curtis Quadrature Rules,
BIT Numerical Mathematics,
Volume 43, Number 1, 2003, pages 1-18.
Parameters:
Input, real ( kind = 8 ) A, B, the endpoints of the interval.
Input, integer ( kind = 4 ) N, the order of the rule.
Output, real ( kind = 8 ) X(N), the abscissas.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=8) | :: | a | ||||
| real(kind=8) | :: | b | ||||
| integer(kind=4) | :: | n | ||||
| real(kind=8) | :: | x(n) |