Article ID Journal Published Year Pages File Type
4645108 Applied Numerical Mathematics 2014 10 Pages PDF
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.
Related Topics
Physical Sciences and Engineering Mathematics Computational Mathematics
Authors
,