NuMFor 9f2ab49 (2024-04-08)
Numerical (Modern) Fortran. Library for Simple Numerical computing
|
Subroutine qextr carries out the Epsilon extrapolation algorithm. More...
Subroutine qextr carries out the Epsilon extrapolation algorithm.
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 given. The condensed epsilon table is computed. Only those elements needed for the computation of the next diagonal are preserved.
[in] | n | (integer) The entry of EPSTAB which contains the new element in the first column of the epsilon table. |
[in] | epstab | (real, array) the two lower diagonals of the triangular epsilon table. The elements are numbered starting at the right-hand corner of the triangle. |
[out] | result | the estimated value of the integral |
[out] | res3la | (real, array dimension(3)) the last three results |
[in,out] | nres | (integer) the number of calls to the routine. This should be zero on the first call, and is automatically updated before return. |