A Short Tale of Long Tail Integration
For scientists and engineers computing Fourier/Laplace-type integrals numerically, this provides a low-cost way to improve accuracy.
The paper presents a simple end-point correction to approximate tail integrals of oscillatory functions, significantly reducing truncation error with negligible extra computation.
Integration of the form $\int_a^\infty {f(x)w(x)dx} $, where $w(x)$ is either $\sin (ω{\kern 1pt} x)$ or $\cos (ω{\kern 1pt} 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.