DQELG(3)

SYNOPSIS
SUBROUTINE DQELG (N, EPSTAB, RESULT, ABSERR, RES3LA, NRES)

NAME
dqelg - The routine determines the limit of a given sequence of approximations, by means of the Epsilon algorithm of P.Wynn. An estimate of the absolute error is also given. The condensed Epsilon table is computed. Only those elements needed for the computation of the next diagonal are preserved.

LIBRARY
SLATEC

TYPE
DOUBLE PRECISION (QELG-S, DQELG-D)

KEYWORDS
CONVERGENCE ACCELERATION, EPSILON ALGORITHM, EXTRAPOLATION

AUTHOR
Piessens, Robert Applied Mathematics and Programming Division K. U. Leuven de Doncker, Elise Applied Mathematics and Programming Division K. U. Leuven

DESCRIPTION
Epsilon algorithm Standard fortran subroutine Double precision version PARAMETERS N - Integer EPSTAB(N) contains the new element in the first column of the epsilon table. EPSTAB - Double precision Vector of dimension 52 containing the elements of the two lower diagonals of the triangular epsilon table. The elements are numbered starting at the right-hand corner of the triangle. RESULT - Double precision Resulting approximation to the integral ABSERR - Double precision Estimate of the absolute error computed from RESULT and the 3 previous results RES3LA - Double precision Vector of dimension 3 containing the last 3 results NRES - Integer Number of calls to the routine (should be zero at first call)

SEE ALSO
DQAGIE, DQAGOE, DQAGPE, DQAGSE

ROUTINES CALLED
D1MACH

REVISION HISTORY (YYMMDD)
800101 DATE WRITTEN
890531 Changed all specific intrinsics to generic. (WRB)
890531 REVISION DATE from Version 3.2
891214 Prologue converted to Version 4.0 format. (BAB)
900328 Added TYPE section. (WRB)