Analysis of errors in some recent numerical quadrature formulas for periodic singular and hypersingular integrals via regularization

Abstract

Recently, we derived some new numerical quadrature formulas of trapezoidal rule type for the singular integrals I(1)[u]=∫ab(cotπ(x−t)T)u(x)dx and I(2)[u]=∫ab(csc2π(x−t)T)u(x)dx, with b−a=T and u(x) a T-periodic continuous function on R. These integrals are not defined in the regular sense, but are defined in the sense of Cauchy Principal Value and Hadamard Finite Part, respectively. With h=(b−a)/n, n=1,2,…, the numerical quadrature formulas Qn(1)[u] for I(1)[u] and Qn(2)[u] for I(2)[u] areQn(1)[u]=h∑j=1nf(t+jh−h/2),f(x)=(cotπ(x−t)T)u(x), andQn(2)[u]=h∑j=1nf(t+jh−h/2)−T2u(t)h−1,f(x)=(csc2π(x−t)T)u(x). We provided a complete analysis of the errors in these formulas under the assumption that u∈C∞(R) and is T-periodic. We actually showed that,I(1)[u]−Qn(1)[u]=O(n−μ)andI(2)[u]−Qn(2)[u]=O(n−μ)as n→∞,∀μ>0. In this note, we analyze the errors in these formulas under the weaker assumption that u∈Cs(R) for some finite integer s. By first regularizing these integrals, we prove that, if u(s+1) is piecewise continuous, thenI(1)[u]−Qn(1)[u]=o(n−s−1/2)as n→∞, if s⩾1,andI(2)[u]−Qn(2)[u]=o(n−s+1/2)as n→∞, if s⩾2. We also extend these results by imposing different smoothness conditions on u(s+1). Finally, we append suitable numerical examples.