AJCMAmerican Journal of Computational Mathematics2161-1203Scientific Research Publishing10.4236/ajcm.2016.62011AJCM-67263ArticlesPhysics&Mathematics The Rectangle Rule for Computing Cauchy Principal Value Integral on Circle JinLi1BenxueGong2WeiLiu3School of Science, Shandong Jianzhu University, Jinan, ChinaSchool of Science, Shandong University of Technology, Zibo, ChinaSchool of Mathematics and Statistics Science, Ludong University, Yantai, China270420160602981078 March 2016accepted 10 June 13 June 2016© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

The classical composite rectangle (constant) rule for the computation of Cauchy principle value integral with the singular kernel is discussed. We show that the superconvergence rate of the composite midpoint rule occurs at certain local coodinate of each subinterval and obtain the corresponding superconvergence error estimate. Then collation methods are presented to solve certain kind of Hilbert singular integral equation. At last, some numerical examples are provided to validate the theoretical analysis.

Cauchy Principal Value Integral Extrapolation Method Composite Rectangle Rule Superconvergence Error Expansion
1. Introduction

Consider the Cauchy principle integral

where denotes a Cauchy principle value integral and s is the singular point.

There are several different definitions which can be proved equally, such as the definition of subtraction of the singularity, regularity definition, direct definition and so on. In this paper we adopt the following one

Cauchy principal value integrals have recently attracted a lot of attention  -  . The main reason for this interest is probably due to the fact that integral equations with Cauchy principal value integrals have shown to be an adequate tool for the modeling of many physical situations, such as acoustics, fluid mechanics, elasticity, fracture mechanics and electromagnetic scattering problems and so on. It is the aim of this paper to investigate the superconvergence phenomenon of rectangle rule for it and, in particular, to derive error estimates.

The superconvergence of composite Newton-Cotes rules for Hadamard finite-part integrals was studied in  -  , where the superconvergence rate and the superconvergence point were presented, respectively. Lyness  derived the Euler-Maclaurin formula for Cauchy principal value integrals. Elliott and Venturino  employed sigmoidal transformations to obtain better approximation to Cauchy principal value integrals. In the reference Avram Sidi   and  presented high-accuracy numerical quadrature methods for integrals of singular periodic functions. The classical Euler-Maclaurin summation formula  expressed the difference between a definite integral over and its approximation using the trapezoidal rule with step length as an asymptotic expansion in powers of h together with a remainder term.

The extrapolation method for the computation of Hadamard finite-part integrals on the interval and in a circle is studied in  and  which focus on the asymptotic expansion of error function. Based on the asymptotic expansion of the error functional, algorithm with theoretical analysis of the generalized extrapolation is given.

In this paper, the density function f(x) is replaced by the approximation function fC(x) while the singular kernel

is computed analysis in each subinterval, where fC(x) is the midpoint rectangle rule. This methods

may be considered as the semi-discrete methods and the order of singularity kernel can be reduced somehow. This idea was firstly presented by Linz  in the paper to calculated the hypersingular integral on interval. He used the trapezoidal rule and Simpson rule to approximate the density function f(x) and the convergence rate was when the singular point was always located at the middle of certain subinterval. This paper focuses on the superconvergence of mid-rectangle rule for Cauchy principle integrals. We prove both theoreti- cally and numerically that the composite mid-rectangle rule reaches the superconvergence rate when the local coordinate of the singular point s is. Then a collation methods is presented to solve certain kind of Hilbert singular integral equation.

The rest of this paper is organized as follows. In Sect. 2, after introducing some basic formulas of the rectangle rule, we present the main resluts. In Sect. 3, we perform the proof. Finally, several numerical examples are provided to validate our analysis.

2. Main Result

Let be a uniform partition of the interval with mesh size. Define by the piecewise constant interpolant for

and a linear transformation

from the reference element to the subinterval. Replacing in (2) with gives the composite rectangle rule:

where denotes the error functional and is the Cote coefficients given by

We also define

Theorem 1: Assume. For the rectangle rule defined as (5). Assume that, there exist a positive constant C, independent of h and s, such that

where

Proof: Let, then we have As

For the first part of (10), we have

For the second part of (10),

Combining (11) and (13) together, the proof is completed.

Setting

Lemma 1: Assume with. Let be defined by (14), then there holds that

Proof: For, by the definition of cauchy principal value integral, we have

For, taking integration by parts on the correspondent Riemann integral, we have

Now, by using the well-known identity

and

The proof is completed.

By the identity in 

then we get

and

For, by the definition of cauchy principal value integral, we have

Let be the function of the second kind associated with the Legendre polynomial, defined by (cf.  )

We also define

Then, by the definition of W,

it follows that

Theorem 2: Assume. For the rectangle rule defined as (5). Assume that, there exist a positive constant C, independent of h and s, such that

where

is defined as (9).

It is known that the global convergence rate of the composite rectangle rule is lower than Riemann integral.

3. Proof of the Theorem

In this section, we study the superconvergence of the composite rectangle rule for Cauchy principle integrals.

Preliminaries

In the following analysis, C will denote a generic constant that is independent of h and s and it may have different values in different places.

Lemma 2: Under the same assumptions of Theorem 2, it holds that

where.

Proof: Performing Taylor expansion of at the point x, we have

and

Combining (29) and (30) together we get the results.

Proof of Theorem 2: we have

For, we have

Putting (31) and (32) together yields

Here

with the linear transformation from to the identity interval. As for the last part of

which can be considered as the error estimate of left rectangle rule for the definite integral. Obviously,by the Theorem, it can be expanded by the Euler-Maclaurin expansions and we have

It is easy to see that there are not relation with the singular point which can be written as

The proof is complete.

We actually obtain the error expansion of the rectangle rule and moreover, get the explicit expression of the first order term. So it is easy for us to get the superconvergence point with, which means that is the superconvergence point in subinterval not near the end of the interval.

Based on the theorem 1, we present the modify rectangle rule

4. Numerical Example

In this section, computational results are reported.

Example 1: We consider the Hilbert singular integral with. with is the superconvergence point.

From Table 1 and Table 2, we know that the superconvergence point is with the coordinate location of singular point equal zero, while for the local coordinate of singular point do not equal zero,it is not convergence in general which coincides with our analysis.

For the modify classical rectangle rule, from Table 3 and Table 4, for the non-superconvergence point and the supersonvergence point, we all get the supercocergence phenomenon.

In this section, we consider the integral equation

with the compatibility condition

As in  , under the condition of (38), there exists a unique solution for the integral Equation (37). In order to get a unique solution, we adopt the following condition

By choosing the middle points, we get the composite rectangle rule to approximate the Hilbert singular integral in (37), then the following linear system is obtained

An error estimate of the rectangle rule<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1100516x102.png" xlink:type="simple"/></inline-formula>
8−1.7764e−01506.6963e−0011.7335e+000
16−1.7764e−01502.2690e+0004.1887e+000
32002.9882e+0005.2932e+000
64−5.3291e−0152.6645e−0153.3190e+0005.8039e+000
128−4.4409e−014−4.3521e−0143.4762e+0006.0477e+000
2561.3323e−014−9.7700e−0153.5526e+0006.1665e+000
An error estimate of the rectangle rule<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1100516x107.png" xlink:type="simple"/></inline-formula>
8−1.1102e−015−6.2172e−0155.0863e+0008.7150e+000
16−3.1086e−015−7.1054e−0154.6011e+0007.8365e+000
32−8.8818e−015−4.3521e−0144.1701e+0007.1371e+000
64−2.6645e−015−1.7764e−0143.9119e+0006.7284e+000
128−2.6645e−015−3.8192e−0143.7729e+0006.5102e+000
2562.9310e−0141.0658e−0133.7010e+0006.3978e+000
An error estimate of the modify rectangle rule<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1100516x112.png" xlink:type="simple"/></inline-formula>
8−1.7764e−0150−1.7764e−0153.5527e−015
16−1.7764e−0150−5.3291e−0157.1054e−015
3200−3.5527e−0152.3093e−014
64−5.3291e−0152.6645e−0151.0658e−0143.0198e−014
128−4.4409e−014−4.3521e−014−4.7962e−014−1.3323e−013
2561.3323e−014−9.7700e−015−6.0396e−0142.8422e−014
An error estimate of the modify rectangle rule<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1100516x117.png" xlink:type="simple"/></inline-formula>
8−1.1102e−015−6.2172e−015−8.8818e−0160
16−3.1086e−015−7.1054e−015−2.8866e−015−3.5527e−015
32−8.8818e−015−4.3521e−014−9.5479e−015−1.5099e−014
64−2.6645e−015−1.7764e−014−7.3275e−015−3.6637e−015
128−2.6645e−015−3.8192e−0144.4409e−015−2.6645e−014
2562.9310e−0141.0658e−0133.1974e−0141.2212e−014

and written as the matrix expression as

where

here denotes the numerical solution of f at. By directly calculation, we get that is not only a symmetric Toeplitz matrix but also a circulant matrix. As for any,

from (43), we know that is singular matrix, then we cannot use system (40) or (41) to solve the integral Equation (37).

In order to get a well-conditioned definite system, we introduce a regularizing factor in (40), which leads to linear system

where defined by

Then the matrix form of system (44) can be presented as

where

and.

Example 2: Now we consider an example of solving Hilbert integral equation by collocation scheme. Let, the exact solution is.

We examine the maximal nodal error, defined by

where denotes the approximation of at. Numerical results presented in Table 5 show that both the maximal nodal errors are as follow.

Acknowledgements

The work of Jin Li was supported by National Natural Science Foundation of China (Grant No. 11471195 and

Errors for the solution of the Hilbert integral equation of first kind
n
322.2204e−16
641.9984e−15
1283.9968e−15
2568.8818e−15
5121.4433e−14

Grant No. 11101247), China Postdoctoral Science Foundation (Grant No. 2013M540541 and 2015T80703). The work of Wei Liu was supported by National Natural Science Foundation of China (Grant No. 11401289).

Cite this paper

Lianju Chen,Xue Su,Yuhong Li,Jin Li,1 1,Benxue Gong,Wei Liu, (2016) The Rectangle Rule for Computing Cauchy Principal Value Integral on Circle. American Journal of Computational Mathematics,06,98-107. doi: 10.4236/ajcm.2016.62011

ReferencesDiethdm K. ,et al. (1995)Asymptotically Sharp Error Bounds for a Quadrature Rule for Cauchy Principal Value Integrals Based on Piecewise Linear Interpolation, Approx Theory of Probability and Its Applications 11, 78-89.Elliott, D. and Venturino, E. (1997) Sigmoidal Transformations and the Euler-Maclaurin Expansion for Evaluating Certain Hadamard Finite-Part Integrals. Numerische Mathematik, 77, 453-465. http://dx.doi.org/10.1007/s002110050295Hasegawa, T. (2004) Uniform Approximations to Finite Hilbert Transform and Its Derivative. Journal of Computational and Applied Mathematics, 163, 127-138. http://dx.doi.org/10.1016/j.cam.2003.08.059Ioakimidis, N.I. (1985) On the Uniform Convergence of Gaussian Quadrature Rules for Cauchy Principal Value Integrals and Their Derivatives. Mathematics of Computation, 44, 191-198. http://dx.doi.org/10.1090/S0025-5718-1985-0771040-8Yu, D.H. (2002) Natural Boundary Integrals Method and Its Applications. Kluwer Academic Publishers, Dordrecht, 50.Yu D.H. ,et al. (1992)The Approximate Computation of Hypersingular Integrals on Interval Numerical Mathematics: A Journal of Chinese Universities (English Series) 1, 114-127.Zhang, X.P., Wu, J.M. and Yu, D.H. (2010) The Superconvergence of Composite Trapezoidal Rule for Hadamard Finite-Part Integral on a Circle and Its Application. International Journal of Computer Mathematics, 87, 855-876. http://dx.doi.org/10.1080/00207160802226517Wu, J.M. and Sun, W.W. (2008) The Superconvergence of Newton-Cotes Rules for the Hadamard Finite-Part Integral on an Interval. Numerische Mathematik, 109, 143-165. http://dx.doi.org/10.1007/s00211-007-0125-7Lyness, J.N. and Ninhan, B.W. (1967) Numerical Quadrature and Asymptotic Expansions. Mathematics of Computation, 21, 162-178. http://dx.doi.org/10.1090/S0025-5718-1967-0225488-XSidi, A. and Israeli, M. (1988) Quadrature Methods for Periodic Singular and Weakly Singular Fredholm Integral Equations. Journal of Scientific Computing, 3, 201-231. http://dx.doi.org/10.1007/BF01061258Sidi, A. (2003) Practical Extrapolation Methods: Theory and Applications. Cambridge University Press, Cambridge. http://dx.doi.org/10.1017/CBO9780511546815Zeng, G., Lei, L. and Huang, J. (2015) A New Construction of Quadrature Formulas for Cauchy Singular Integral. Journal of Computational Analysis and Applications, 17, 426-436.Sidi, A. (2014) Analysis of Errors in Some Recent Numerical Quadrature Formulas for Periodic Singular and Hypersingular Integrals via Regularization. Applied Numerical Mathematics, 81, 30-39. http://dx.doi.org/10.1016/j.apnum.2014.02.011Li, J., Wu, J.M. and Yu, D.H. (2009) Generalized Extrapolation for Computation of Hypersingular Integrals in Boundary Element Methods. CMES: Computer Modeling in Engineering & Sciences, 42, 151-175.Li, J., Zhang, X.P. and Yu, D.H. (2013) Extrapolation Methods to Compute Hypersingular Integral in Boundary Element Methods. Science China Mathematics, 56, 1647-1660. http://dx.doi.org/10.1007/s11425-013-4593-1Linz, P. (1985) On the Approximate Computation of Certain Strongly Singular Integrals. Computing, 35, 345-353. http://dx.doi.org/10.1007/BF02240199Andrews, L.C. (1992) Special Functions of Mathematics for Engineers. 2nd Edition, McGraw-Hill, Inc., Cambridge, 258.