Quadrature rules for weakly singular, strongly singular, and hypersingular integrals in boundary integral equation methods

We present n-point Gauss-Gegenbauer quadrature rules for weakly singular, strongly singular, and hypersingular integrals that arise in integral equation formulations of potential problems in domains with edges and corners. The rules are tailored to weight functions with algebraic endpoint singularities related to the geometrical singularities of the domain. Each rule has two different expressions involving Legendre functions and hypergeometric functions, respectively. Numerical examples amply demonstrate the accuracy and stability of the proposed algorithms. Application to the solution of a singular integral equation is exemplified.