Wolfram Library Archive

Courseware Demos MathSource Technical Notes
All Collections Articles Books Conference Proceedings

Extended Gaussian quadratures for functions with an end-point singularity of logarithmic type

K. Pachucki
M. Puchalski
Organization: Adam Mickiewicz University
V.A. Yerokhin
Organization: St. Petersburg State Polytechnical University
Journal / Anthology

Computer Physics Communications
Year: 2014
Volume: 185
Page range: 29132919

The extended Gaussian quadrature rules are shown to be an efficient tool for numerical integration of wide class of functions with singularities of logarithmic type. The quadratures are exact for the functions pol1n−1(x) + ln x pol2n−1(x), where pol1n−1(x) and pol2n−1(x) are two arbitrary polynomials of degree n − 1 and n is the order of the quadrature formula. We present an implementation of numerical algorithm that calculates the nodes and the weights of the quadrature formulas, provide a Fortran code for numerical integration, and test the performance of different kinds of Gaussian quadratures for functions with logarithmic singularities.

*Science > Physics