A short tale of long tail integration

  • Authors:
  • Xiaolin Luo;Pavel V. Shevchenko

  • Affiliations:
  • CSIRO Mathematics, Informatics and Statistics, Sydney, Australia 1670;CSIRO Mathematics, Informatics and Statistics, Sydney, Australia 1670

  • Venue:
  • Numerical Algorithms
  • Year:
  • 2011

Quantified Score

Hi-index 0.00

Visualization

Abstract

Integration of the form $\int_a^\infty {f(x)w(x)dx} $ , where w(x) is either sin(驴x) or cos(驴x), is widely encountered in many engineering and scientific applications, such as those involving Fourier or Laplace transforms. Often such integrals are approximated by a numerical integration over a finite domain (a, b), leaving a truncation error equal to the tail integration $\int_b^\infty {f(x)w(x)dx} $ in addition to the discretization error. This paper describes a very simple, perhaps the simplest, end-point correction to approximate the tail integration, which significantly reduces the truncation error and thus increases the overall accuracy of the numerical integration, with virtually no extra computational effort. Higher order correction terms and error estimates for the end-point correction formula are also derived. The effectiveness of this one-point correction formula is demonstrated through several examples.