NuMFor 9f2ab49 (2024-04-08)
Numerical (Modern) Fortran. Library for Simple Numerical computing
qextr Interface Reference

Subroutine qextr carries out the Epsilon extrapolation algorithm. More...

Detailed Description

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.

Parameters
[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]resultthe 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.

The documentation for this interface was generated from the following file: