High-order quadratures for integral operators with singular kernels

A numerical integration method that has rapid convergence for integrands with known singularities is presented. Based on endpoint corrections to the trapezoidal rule, the quadratures are suited for the discretization of a variety of integral equations encountered in mathematical physics. The quadratures are based on a technique introduced by Rokhlin (1990). The present modification controls the growth of the quadrature weights and permits higher-order rules in practice. Several numerical examples are included.