Further extension of a class of periodizing variable transformations for numerical integration

摘要

Class Sm variable transformations with integer m, for accurate numerical computation of finite-range integrals via the trapezoidal rule, were introduced and studied by the author. A representative of this class is the sinm-transformation. In a recent work of the author, this class was extended to arbitrary noninteger values of m, and it was shown that exceptionally high accuracies are achieved by the trapezoidal rule in different circumstances with suitable values of m. In another recent work by Monegato and Scuderi, the sinm-transformation was generalized by introducing two integers p and q, instead of the single integer m; we denote this generalization as the sinp,q-transformation here. When p=q=m, the sinp,q-transformation becomes the sinm-transformation. Unlike the sinm-transformation which is symmetric, the sinp,q-transformation is not symmetric when p≠q, and this offers an advantage when the behavior of the integrand at one endpoint is quite different from that at the other endpoint. In view of the developments above, in the present work, we generalize the class Sm by introducing a new class of nonsymmetric variable transformations, which we denote as Sp,q, where p and q can assume arbitrary noninteger values, such that the sinp,q-transformation is a representative of this class and Sm⊂Sm,m. We provide a detailed analysis of the trapezoidal rule approximation following a variable transformation from the class Sp,q, and show that, with suitable and not necessarily integer p and q, it achieves an unusually high accuracy when the integrand has algebraic endpoint singularities. We also illustrate our results with numerical examples via the sinp,q-transformation. Finally, we discuss the computation of surface integrals in R3 containing point singularities with the help of class Sp,q transformations.