Compact Numerical Quadrature Formulas for Hypersingular Integrals and Integral Equations

In the first part of this work, we derive compact numerical quadrature formulas for finite-range integrals , where f(x)=g(x)|x-t| (beta) , beta being real. Depending on the value of beta, these integrals are defined either in the regular sense or in the sense of Hadamard finite part. Assuming that gaC (a)[a,b], or gaC (a)(a,b) but can have arbitrary algebraic singularities at x=a and/or x=b, and letting h=(b-a)/n, n an integer, we derive asymptotic expansions for , where x (j) =a+jh and ta{x (1),aEuro broken vertical bar,x (n-1)}. These asymptotic expansions are based on some recent generalizations of the Euler-Maclaurin expansion due to the author (A. Sidi, Euler-Maclaurin expansions for integrals with arbitrary algebraic endpoint singularities, in Math. Comput., 2012), and are used to construct our quadrature formulas, whose accuracies are then increased at will by applying to them the Richardson extrapolation process. We pay particular attention to the case in which beta=-2 and f(x) is T-periodic with T=b-a and , which arises in the context of periodic hypersingular integral equations. For this case, we propose the remarkably simple and compact quadrature formula , and show that as h -> 0 is not an element of mu > 0, and that it is exact for a class of singular integrals involving trigonometric polynomials of degree at most n. We show how can be used for solving hypersingular integral equations in an efficient manner. In the second part of this work, we derive the Euler-Maclaurin expansion for integrals , where f(x)=g(x)(x-t) (beta) , with g(x) as before and beta=-1,-3,-5,aEuro broken vertical bar, from which suitable quadrature formulas can be obtained. We revisit the case of beta=-1, for which the known quadrature formula satisfies as h -> 0 is not an element of mu > 0, when f(x) is T-periodic with T=b-a and . We show that this formula too is exact for a class of singular integrals involving trigonometric polynomials of degree at most n-1. We provide numerical examples involving periodic integrands that confirm the theoretical results.