Abstract
A new method and its computer program for the double exponential formula for numerical integration is presented. The present method avoids underflow encountered in the computation of weights of the formula and sums the integral up to a sufficiently large number of terms, so that it does not deteriorate the accuracy of the double exponential formula even for integrands with strong end-point singularity.
