Abstract
We propose a general method for computing indefinite integrals of the form
where g is a smooth function, and k is a function that contains a singular factor or is rapidly oscillatory. The only assumption on k is that it satisfies a linear differential equation with polynomial coefficients. The approximate value of the integral is given in terms of Chebyshev coefficients of functions that form a solution of a certain system of differential equations. As an illustration, we present effective algorithms for computing indefinite integrals of the functions g(t)|t–d|α e iωt, g(t)log|t–d| e iωt, g(t) t α J ν (ct).
Similar content being viewed by others
References
Abramowitz, M., Stegun, I.A.: Handbook of mathematical functions. National Bureau of Standards, Washington DC (1967)
Bakhvalov, N.S., Vasil’eva, L.G.: Evaluation of the integrals of oscillating functions by interpolation at nodes of Gaussian quadrature. Z. Vycisl. Mat. Mat. Fiz. 8, 175–181 (1968)
Branders, M., Piessens, R.: Modified Clenshaw–Curtis method for the computation of Bessel function integrals. BIT 23, 370–381 (1983)
Branders, M., Piessens, R.: A survey of numerical methods for the computation of Bessel function integrals. Special functions: theory and computation, vol. Fasc. spec., Rend. Sem. Mat. Univ. Polit. Torino, pp. 249–265 (1985)
Branders, M., Piessens, R.: On the computation of Fourier transforms of singular functions. J. Comp. Appl. Math. 43, 159–169 (1992)
Chung, K.C., Evans, G.A., Webster, J.R.: A method to generate generalized quadrature rules for oscillatory integrals. Appl. Numer. Math. 34, 85–93 (2000)
Filon, L.N.G.: On a quadrature formula for trigonometric integrals. Proc. Roy. Soc. Edinburgh 49, 38–47 (1929)
Flinn, E.A.: A modification of Filon’s method of numerical integration. J. Assoc. Comp. Mach. 7, 181–184 (1960)
Gautschi, W.: Computational aspects of three-term recurrence relations. SIAM Rev. 9, 24–82 (1967)
Gentleman, W.M.: Implementing Clenshaw–Curtis quadrature, I Methodology and experience. Comm. ACM 15, 337–342 (1972)
Gentleman, W.M.: Implementing Clenshaw–Curtis quadrature, II Computing the cosine transformation. Comm. ACM 15, 343–346 (1972)
Hasegawa, T., Torii, T.: Indefinite integration of oscillatory functions by the Chebyshev series expansion. J. Comp. Appl. Math. 17, 21–29 (1987)
Hasegawa, T., Torii, T.: Application of a modified FFT to product type integration. J. Comp. Appl. Math. 38, 157–168 (1991)
Hasegawa, T., Torii, T.: An algorithm for nondominant solutions of linear second-order inhomogeneous difference equations. Math. Comp. 64, 1199–1204 (1995)
Hasegawa, T., Torii, T., Sugiura, H.: An algorithm based on the FFT for a generalized Chebyshev interpolation. Math. Comp. 54, 195–210 (1990)
Iserles, A., Nørsset, S.: On quadrature methods for highly oscillatory integrals and their implementation. BIT 44, 755–772 (2004)
Keller, P.: Obliczanie całek nieoznaczonych funkcji osobliwych lub oscylujących. PhD thesis, Institute of Computer Science, Wrocław (2003)
Keller, P.: Indefinite integration of oscillatory functions. Applic. Math. 25, 301–311 (1998)
Levin, D.: Procedures for computing one- and two-dimensional integrals of functions with rapid irregular oscillations. Math. Comp. 38, 531–538 (1982)
Levin, D.: Fast integration of rapidly oscillatory functions. J. Comp. Appl. Math. 67, 95–101 (1996)
Lewanowicz, S.: Evaluation of Bessel function integrals with algebraic singularities. J. Comp. Appl. Math. 37, 101–112 (1991)
Lewanowicz, S.: A new approach to the problem of constructing recurrence relations for the Jacobi coefficients. Applic. Math. 21, 303–326 (1991)
Lewanowicz, S.: Quick construction of recurrence relations for the Jacobi coefficients. J. Comp. Appl. Math. 43, 355–372 (1992)
Lozier, D.W.: Numerical solution of linear difference equations. Tech. report, Math. Analysis Division, National Bureau of Standards, Washington, DC (1980)
Luke, Y.L.: On the computation of oscillatory integrals. Proc. Cambridge Philos. Soc. 50, 269–277 (1954)
Luke, Y.L.: The special functions and their approximations, vol. II. Academic Press, New York (1969)
Olver, S.: Moment-free numerical integration of highly oscillatory functions. IMA J. Num. Anal. 26, 213–227 (2006)
Piessens, R., Branders, M.: Computation of Fourier transform integrals using Chebyshev series expansions. Computing 32, 177–186 (1984)
Piessens, R., Poleunis, F.: A numerical method for the integration of oscillatory functions. BIT 11, 317–327 (1971)
Wimp, J.: Computation with recurrence relations. Applicable Mathematics Series, Pitman, London (1984)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Keller, P. A method for indefinite integration of oscillatory and singular functions. Numer Algor 46, 219–251 (2007). https://doi.org/10.1007/s11075-007-9134-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-007-9134-y