Numerical method for hypersingular integrals of highly oscillatory functions on the positive semiaxis

This paper deals with a quadrature rule for the numerical evaluation of hypersingular integrals of highly oscillatory functions on the positive semiaxis. The rule is of product type and consists in approximating the density function f by a truncated interpolation process based on the zeros of generalized Laguerre polynomials and an additional point. We prove the stability and the convergence of the rule, giving error estimates for functions belonging to weighted Sobolev spaces equipped with uniform norm. We also show how the proposed rule can be used for the numerical solution of hypersingular integral equations. Numerical tests which confirm the theoretical estimates and comparisons with other existing quadrature rules are presented.