In this paper, we concerned about the Cauchy principal value in oscillatory kind function which , specifically have the form like
where is smooth and , generally is much large. Li  uses Levin method to solve the integral with Cauchy principal value. In that case, the principal value point is and , but the function . And the Levin method uses Chebyshev-Gauss-Labatt nodes to make a matrix equation in order to get the approximate solution. Then Wang  has mentioned a technique of put the integral to complex plane so that we can use a loop integral. In this paper, the author chooses several paths to make it available tot Gauss-Laguerre quadrature. It is a great idea for us to solve Cauchy principal value problem and we are using this idea to solve problem in this paper. Furthermore, Daan  puts forward the form of integral form. But in his paper, the integral without Cauchy principal value. Importantly, we take his idea of choose path, which make it easy to have a good path to integral. Alfrddo  then put the question to a more general situation. With the function not smooth, maybe has stationary points and so on, they use Gauss-Fraud quadrature to solve the path near the stationary points.
In this paper, we firstly transformed the equation into the form which we known before, then we take the path to get the solution with Gauss-Laguerre quadrature. Finally two numerical examples are given to prove the method. One is for real part, and the other is for imaginary part.
2. Integral with Gauss-Laguerre Quadrature
As we can see, the variable x has a different order. First we need to transform x to let them have the same order. The specific process is
In last equation, we use the variable substation of to transform it. The equation can be calculated in some methods like Li  has put forward.
In such complex plane, the integral can be divided into several pieces which can be seen in above figure. Therefore, we only need to calculate the integral of path so that we can have the integral of path .
The path and have similar situation. have mentioned a method about select a path to infinity as following
where is the path start in x and end in infinity, is the function , and i is the imaginary unit, which is , p is the variable. Now we prove the start point and the end point. Because is smooth, so
When , , when , . In this case, because of the known form of integral, the path can be written as following 
We need to notice that, p is the variable of , and x is the start point of the path.
where in the last equation.
We use polar coordinates to solve . Let , where r is the radio and θ is the variable.
When , the equation above can be transformed into
As we discussed above, the integral () now can be written as following
Fortunately, the infinite integral can solve with Gauss-Laguerre quadrature with n-points. The represents Gauss-Laguerre points and the represents Gauss-Laguerre weights.
Remark 1 is the quadrature rule that we mentioned before with n points and weights. The approximation error behaves like
3. Numerical Example
In order to prove this method is effective, we will use two examples to show the approximate solution and exact solution. All calculation is computed in Matlab2019a, and exact solution is given by Mathematica 12.0 with 64-digit.
Example1. The integral is analytic expression. Using above method, we can use to approach it which have the form like following
Table 1 demonstrates the errors in real part. We can draw a conclusion that, with ω and n increasing, the accuracy of method is increasing. We take as example, when the results approach the max accuracy. When , this process moves up when .
Table 1. The real part error of .
Table 2. The imaginary part error of .
Example 2. The integral of only have imaginary part, so in this example, we can see the performance in imaginary part. This time we take , and we still take to compare the efficiency about the method.
As we can see from Table 2, with n-points Gauss-Laguerre quadrature, the convergence is fast and only small n can obtain satisfying solution.
In this paper, we discuss the oscillatory kind function with Cauchy principal value which has a different form. The integral can be transformed into some known situation so that we can handle it easier. The method that puts the integral in complex plane in this paper comes from Wang  and the numerical solution is shown in Table 1 and Table 2. Within the two examples, we compare the real part and imaginary part between numerical solution and exact solution. And the absolute error represents the method for this integral is available.
The author would like to acknowledge to my tutor for teaching and useful discussion, and anonymous referees for suggestions.
 Li, J.B., Wang, X.S. and Wang, T. (2010) Evaluation of Cauchy Principal Value Integrals of Oscillatory Kind. Applied Mathematics and Computation, 217, 2390-2396. https://doi.org/10.1016/j.amc.2010.07.039
 Wang, H.Y. and Xiang, S.H. (2009) On the Evaluation of Cauchy Principal Value Integrals of Oscillatory Functions. Journal of Computational and Applied Mathematics, 234, 95-100. https://doi.org/10.1016/j.cam.2009.12.007
 Huybrechs, D. and Vandewalle, S. (2006) On the Evaluation of Highly Oscillatory Integrals by Analytic Continuation. SIAM Journal on Numerical Analysis, 44, 1026-1048. https://doi.org/10.1137/050636814