An upper bound estimate and stability for the global error of numerical integration using double exponential transformation

Abstract

The double exponential formula was introduced for calculating definite integrals with singular point oscillation functions and Fourier integral. The double exponential transformation is not only useful for numerical computations but it is also used in different methods of Sinc theory. In this paper we give an upper bound estimate for the error of double exponential transformation. By improving integral estimates having singular final points, in theorem 1 we prove that the method is convergent and the rate of convergence is O(h 2 ) where h is a step size. Our main tool in the proof is DE formula in Sinc theory. The advantage of our method is that the time and space complexity is drastically reduced. Furthermore, we discovered upper bound error in DE formula independent of N truncated number, as a matter of fact we proved stability. Numerical tests are presented to verify the theoretical predictions and confirm the convergence of the numerical solution. MSC: 65D30, 65D32.