Tanh-sinh quadrature
From Wikipedia, the free encyclopedia
Tanh-sinh quadrature is a method for numerical integration introduced by Hidetosi Takahasi and Masatake Mori in 1974. It uses the change of variables
to transform an integral on (−1, 1) to an integral on the entire real line. For a given step size h, the integral is approximated by the sum
with the abscissas
and the weights
![w_k = \frac{\pi/2 \,\cosh kh}{\cosh^2[\pi/2 \,\sinh kh]}.](../../../math/6/e/a/6ea8d1e12c236838042236f8eed551c9.png)
Like Gaussian quadrature, tanh-sinh quadrature is well suited for arbitrary-precision integration, where an accuracy of hundreds or thousands of digits is desired. The convergence is quadratic for sufficiently well-behaved integrands: doubling the number of evaluation points roughly doubles the number of correct digits.
Tanh-sinh quadrature is less efficient than Gaussian quadrature for smooth integrands, but unlike Gaussian quadrature tends to work equally well with integrands having singularities or infinite derivatives at one or both endpoints of the integration interval. A further advantage is that the abscissas and weights are relatively easy to compute. The cost of calculating abscissa-weight pairs for n-digit accuracy is roughly n2 log2 n compared to n3 log n for Gaussian quadrature.
Upon comparing the scheme to Gaussian quadrature and error function quadrature, Bailey et al. found that the tanh-sinh scheme "appears to be the best for integrands of the type most often encountered in experimental math research".
[edit] References
- Jonathan Borwein, David H. Bailey, and Roland Girgensohn, Experimentation in Mathematics - Computational Paths to Discovery. A K Peters, 2003. ISBN 1-56881-136-5.
- David H. Bailey, Karthik Jeyabalan, and Xiaoye S. Li, "A comparison of three high-precision quadrature schemes". Experimental Mathematics, 14.3 (2005).
