Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)

 
 

 

Euler-Maclaurin expansions for integrals with arbitrary algebraic endpoint singularities


Author: Avram Sidi
Journal: Math. Comp. 81 (2012), 2159-2173
MSC (2010): Primary 30E15, 40A25, 41A60, 65B15, 65D30
DOI: https://doi.org/10.1090/S0025-5718-2012-02597-X
Published electronically: April 10, 2012
MathSciNet review: 2945150
Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: In this paper, we provide the Euler-Maclaurin expansions for (offset) trapezoidal rule approximations of the divergent finite-range integrals $ \int ^b_af(x)\,dx$, where $ f\in C^{\infty }(a,b)$ but can have arbitrary algebraic singularities at one or both endpoints. We assume that $ f(x)$ has asymptotic expansions of the general forms

  $\displaystyle f(x)\sim K\,(x-a)^{-1}+\sum ^{\infty }_{s=0}c_s(x-a)^{\gamma _s}$$\displaystyle \quad \text {as}\ x\to a+,$    
  $\displaystyle f(x)\sim L\,(b-x)^{-1}+\sum ^{\infty }_{s=0}d_s(b-x)^{\delta _s}$$\displaystyle \quad \text {as}\ x\to b-,$    

where $ K,L$, and $ c_s, d_s$, $ s=0,1,\ldots ,$ are some constants, $ \vert K\vert+\vert L\vert\neq 0,$ and $ \gamma _s$ and $ \delta _s$ are distinct, arbitrary and, in general, complex, and different from $ -1$, and satisfy

$\displaystyle \Re \gamma _0\leq \Re \gamma _1\leq \cdots , \ \ \lim _{s\to \inf... ...\leq \Re \delta _1\leq \cdots , \ \ \lim _{s\to \infty }\Re \delta _s=+\infty .$

Hence the integral $ \int ^b_af(x)\,dx$ exists in the sense of Hadamard finite part. The results we obtain in this work extend some of the results in [A. Sidi, Numer. Math. 98 (2004), pp. 371-387] that pertain to the cases in which $ K=L=0.$ They are expressed in very simple terms based only on the asymptotic expansions of $ f(x)$ as $ x\to a+$ and $ x\to b-$. With $ h=(b-a)/n$, where $ n$ is a positive integer, one of these results reads

$\displaystyle h\sum ^{n-1}_{i=1}f(a+ih)\sim I[f]$ $\displaystyle +K\,(C -\log h) + \sum ^{\infty }_{\substack {s=0\\ \gamma _s\not \in \{2,4,\ldots \}}}c_s \zeta (-\gamma _s)\,h^{\gamma _s+1}$    
  $\displaystyle +L\,(C -\log h) +\sum ^{\infty }_{\substack {s=0\\ \delta _s\not \in \{2,4,\ldots \}}}d_s\zeta (-\delta _s) h^{\delta _s+1}$$\displaystyle \quad \text {as $h\to 0$},$    

where $ I[f]$ is the Hadamard finite part of $ \int ^b_af(x)\,dx$, $ C$ is Euler's constant and $ \zeta (z)$ is the Riemann Zeta function. We illustrate the results with an example.


References [Enhancements On Off] (What's this?)

  • 1. G.E. Andrews, R. Askey, and R. Roy, Special Functions, Cambridge University Press, Cambridge, 1999. MR 1688958 (2000g:33001)
  • 2. J.S. Brauchart, D.P. Hardin, and E.B. Saff, The Riesz energy of the $ N$th roots of unity: an asymptotic expansion for large $ N$, Bull. London Math. Soc. 41 (2009), 621-633. MR 2521357 (2010g:31001)
  • 3. P.J. Davis and P. Rabinowitz, Methods of Numerical Integration, second ed., Academic Press, New York, 1984. MR 760629 (86d:65004)
  • 4. G. Evans, Practical Numerical Integration, Wiley, New York, 1993. MR 1245265 (94k:65001)
  • 5. P.K. Kythe and M.R. Schäferkotter, Handbook of Computational Methods for Integration, Chapman & Hall/CRC Press, New York, 2005. MR 2129972
  • 6. Y.L. Luke, The Special Functions and Their Approximations, vol. I, Academic Press, New York, 1969. MR 0241700 (39:3039)
  • 7. J.N. Lyness, Finite-part integrals and the Euler-Maclaurin expansion, Approximation and Computation (Boston-Basel-Berlin) (R.V.M. Zahar, ed.), ISNM, no. 119, Birkhäuser Verlag, 1994, pp. 397-407. MR 1333632 (96a:65029)
  • 8. J.N. Lyness and B.W. Ninham, Numerical quadrature and asymptotic expansions, Math. Comp. 21 (1967), 162-178. MR 0225488 (37:1081)
  • 9. G. Monegato and J.N. Lyness, The Euler-Maclaurin expansion and finite-part integrals, Numer. Math. 81 (1998), 273-291. MR 1657780 (99m:65049)
  • 10. I. Navot, An extension of the Euler-Maclaurin summation formula to functions with a branch singularity, J. Math. and Phys. 40 (1961), 271-276. MR 0140876 (25:4290)
  • 11. -, A further extension of the Euler-Maclaurin summation formula, J. Math. and Phys. 41 (1962), 155-163.
  • 12. B.W. Ninham, Generalised functions and divergent integrals, Numer. Math. 8 (1966), 444-457. MR 0196935 (33:5119)
  • 13. F.W.J. Olver, D.W. Lozier, R.F. Boisvert, and C.W. Clark, editors, NIST Handbook of Mathematical Functions, Cambridge University Press, Cambridge, 2010. MR 2723248
  • 14. A. Sidi, Practical Extrapolation Methods: Theory and Applications, Cambridge Monographs on Applied and Computational Mathematics, no. 10, Cambridge University Press, Cambridge, 2003. MR 1994507 (2004e:65005)
  • 15. -, Euler-Maclaurin expansions for integrals with endpoint singularities: a new perspective, Numer. Math. 98 (2004), 371-387. MR 2092747 (2005g:65012)
  • 16. -, Extension of a class of periodizing variable transformations for numerical integration, Math. Comp. 75 (2006), 327-343. MR 2176402 (2006g:41066)
  • 17. -, A novel class of symmetric and nonsymmetric periodizing variable transformations for numerical integration, J. Sci. Comput. 31 (2007), 391-417. MR 2320555 (2008f:65046)
  • 18. -, Further extension of a class of periodizing variable transformations for numerical integration, J. Comp. Appl. Math. 221 (2008), 132-149. MR 2458757 (2010d:65057)
  • 19. J.F. Steffensen, Interpolation, second ed., Dover, New York, 2006.

Similar Articles

Retrieve articles in Mathematics of Computation with MSC (2010): 30E15, 40A25, 41A60, 65B15, 65D30

Retrieve articles in all journals with MSC (2010): 30E15, 40A25, 41A60, 65B15, 65D30


Additional Information

Avram Sidi
Affiliation: Computer Science Department, Technion - Israel Institute of Technology, Haifa 32000, Israel
Email: asidi@cs.technion.ac.il

DOI: https://doi.org/10.1090/S0025-5718-2012-02597-X
Keywords: Euler–Maclaurin expansions, trapezoidal rule, endpoint singularities, Hadamard finite part, asymptotic expansions.
Received by editor(s): November 16, 2010
Received by editor(s) in revised form: April 17, 2011
Published electronically: April 10, 2012
Additional Notes: This research was supported in part by the United States–Israel Binational Science Foundation grant no. 2008399.
Article copyright: © Copyright 2012 American Mathematical Society
The copyright for this article reverts to public domain 28 years after publication.

American Mathematical Society