Citation: Robert Reynolds, Allan Stauffer. Derivation of some integrals in Gradshteyn and Ryzhik[J]. AIMS Mathematics, 2021, 6(2): 1816-1821. doi: 10.3934/math.2021109
[1] | Robert Reynolds, Allan Stauffer . Definite integrals involving product of logarithmic functions and logarithm of square root functions expressed in terms of special functions. AIMS Mathematics, 2020, 5(6): 5724-5733. doi: 10.3934/math.2020367 |
[2] | Robert Reynolds, Allan Stauffer . Definite integral of the logarithm hyperbolic secant function in terms of the Hurwitz zeta function. AIMS Mathematics, 2021, 6(2): 1324-1331. doi: 10.3934/math.2021082 |
[3] | Ye-Cong Han, Chuan-Yu Cai, Ti-Ren Huang . Monotonicity, convexity properties and inequalities involving Gaussian hypergeometric functions with applications. AIMS Mathematics, 2022, 7(4): 4974-4991. doi: 10.3934/math.2022277 |
[4] | A. Belafhal, N. Nossir, L. Dalil-Essakali, T. Usman . Integral transforms involving the product of Humbert and Bessel functions and its application. AIMS Mathematics, 2020, 5(2): 1260-1274. doi: 10.3934/math.2020086 |
[5] | Robert Reynolds, Allan Stauffer . Extended Prudnikov sum. AIMS Mathematics, 2022, 7(10): 18576-18586. doi: 10.3934/math.20221021 |
[6] | Marta Na Chen, Wenchang Chu . Yabu's formulae for hypergeometric $ _3F_2 $-series through Whipple's quadratic transformations. AIMS Mathematics, 2024, 9(8): 21799-21815. doi: 10.3934/math.20241060 |
[7] | Robert Reynolds, Allan Stauffer . A quintuple integral involving the product of Hermite polynomial $ H_{n}(\beta x) $ and parabolic cylinder function $ D_{v}(\alpha t) $: derivation and evaluation. AIMS Mathematics, 2022, 7(5): 7464-7470. doi: 10.3934/math.2022418 |
[8] | Abdelhamid Zaidi, Saleh Almuthaybiri . Explicit evaluations of subfamilies of the hypergeometric function $ _3F_2(1) $ along with specific fractional integrals. AIMS Mathematics, 2025, 10(3): 5731-5761. doi: 10.3934/math.2025264 |
[9] | Robert Reynolds . A short note on a extended finite secant series. AIMS Mathematics, 2023, 8(11): 26882-26895. doi: 10.3934/math.20231376 |
[10] | Mohamed Akel, Muajebah Hidan, Salah Boulaaras, Mohamed Abdalla . On the solutions of certain fractional kinetic matrix equations involving Hadamard fractional integrals. AIMS Mathematics, 2022, 7(8): 15520-15531. doi: 10.3934/math.2022850 |
We will derive integrals as indicated in the abstract in terms of the Lerch function. The derivations follow the method used by us in [6]. A generalization of Cauchy's integral formula is given by
ykΓ(k+1)=12πi∫Cewywk+1dw | (1.1) |
where C will be defined below. This method involves using a form of Eq (1.1) then multiplying both sides by a function, then take a definite integral of both sides. This yields a definite integral in terms of a contour integral. Then we multiply both sides of Eq (1.1) by another function and take the infinite sum of both sides such that the contour integral of both equations are the same.
We use the method in [6]. Using a generalization of Cauchy's integral formula we multiply both sides by in the first term eiay and replacing y by iy+log(z) and in the second by e−iay and replacing y by −iy+log(z) then add these two equations, followed by multiplying both sides by 12cosh(by) to get
e−iay(−iy+log(z))k+eiay(iy+log(z))k2Γ(k+1)cosh(by)=12πi∫Czww−k−1sech(by)cos((a+w)y)dw | (2.1) |
where the logarithmic function is defined in Eq (4.2) in [2]. The variables a, b, k and z are general complex numbers. Define →a=(ℜ(a),ℑ(a)), →b=(ℜ(a),ℑ(b))) and →w=(ℜ(w),ℑ(w)). Then the integral over the contour C has w replaced by w+a. We can define the cut and contour C as lying on the side of →w+→a opposite to →b, coming from infinity initially parallel to →w+→a, around the origin with zero radius and back to infinity parallel to →w+→a. The contour lies on opposite sides of the cut. We then take the definite integral over y∈[0,∞) of both sides to get
1Γ(k+1)∫∞0e−iay(−iy+log(z))k+eiay(iy+log(z))k2cosh(by)dy=14πi∫∞0∫Czww−k−1sech(by)cos((a+w)y)dwdy=14πi∫C(∫∞0sech(by)cos((a+w)y)dy)zwdwwk+1=14ib∫Czww−k−1sech(π(a+w)2b)dw | (2.2) |
from Eq (1.7.7.1) in [1]. In a similar manner we can derive the second integral formula by subtracting the two equations and multiplying by 12isinh(by) we get
eiay(iy+log(z))k−e−iay(−iy+log(z))kΓ(k+1)2isinh(by)=12πi∫Czww−k−1csch(by)sin((a+w)y)dw | (2.3) |
We then take the definite integral over y∈[0,∞) of both sides to get
1k!∫∞0eiay(iy+log(z))k−e−iay(−iy+log(z))k2isinh(by)dy=12πi∫∞0∫Czww−k−1csch(by)sin((a+w)y)dwdy=12πi∫C(∫∞0csch(by)sin((a+w)y)dy)zwdwwk+1=14bi∫Czww−k−1tanh(π(a+w)2b)dw | (2.4) |
from Eq (2.7.7.6) in [1].
Again, using the method in [6], replacing y with π(2p+1)/(2b)+log(z) and multiplying both sides by (−1)p(πb)eaπ(2p+1)/(2b) to yield
(−1)p(πb)ea(π(2p+1)/(2b))(π(2p+1)/(2b)+log(z))kk!=(−1)p2πi(πb)∫Cew(π(2p+1)2b)+aπ(2p+1)2bwk+1zwdw | (3.1) |
followed by taking the infinite sum of both sides of Eq (3.1) with respect to p over [0,∞) to get
b−k−1eaπ2bπk+1k!Φ(−eaπb,−k,12+blog(z)π)=12ib∞∑p=0∫Cew(π(2p+1)2b)+aπ(2p+1)2b(−1)−pwk+1zwdw=12ib∫C∞∑p=0ew(π(2p+1)2b)+aπ(2p+1)2b(−1)−pwk+1zwdw=14ib∫Czww−k−1sech(π(a+w)2b)dw | (3.2) |
from (1.232.1) in [5] where ℜ((w+a)/b)<0 for the convergence of the sum and if the ℜ(k)<0 then the argument of the sum over p cannot be zero for some value of p. We use (9.550) in [5] where Φ(z,s,v) is the Lerch function which is a generalization of the Hurwitz Zeta and polylogarithm functions.
The Lerch function has a series representation given by
Φ(z,s,v)=∞∑n=0(v+n)−szn | (3.3) |
where |z|<1,v≠0,−1,.. and is continued analytically by its integral representation given by
Φ(z,s,v)=1Γ(s)∫∞0ts−1e−vt1−ze−tdt=1Γ(s)∫∞0ts−1e−(v−1)tet−zdt | (3.4) |
where ℜ(v)>0, or |z|≤1,z≠1,ℜ(s)>0, or z=1,ℜ(s)>1. Similarly, using the method in [6], replacing y with 2π(p+1)/(2b)+log(z) and multiplying both sides by (−1)p(πb)e2aπ(p+1)/(2b) to yield
(−1)p(πb)e2aπ(p+1)2b(2π(p+1)2b+log(z))kk!=(−1)p2πi(πb)∫Cew(2π(p+1)/(2b))+2aπ(p+1)/(2b)wk+1zwdw | (3.5) |
followed by taking the infinite sum of both sides of Eq (3.5) with respect to p over [0,∞) to get
eaπbπk+1b−k−1k!Φ(−eaπb,−k,1+blog(z)π)=12bi∞∑p=0∫Cew(2π(p+1)/(2b))+2aπ(p+1)/(2b)(−1)−pwk+1zwdw=12bi∫C∞∑p=0ew(2π(p+1)/(2b))+2aπ(p+1)/(2b)(−1)−pwk+1zwdw=14bi∫C(zww−k−1+zww−k−1tanh(π(a+w)2b))dw | (3.6) |
from Eq (1.232.1) in [5] where ℜ((a+w)/b)<0.
Since the right hand-side of Eq (2.2) is equal to the right-hand side of (3.2), we can equate the left hand-sides and simplify to get
∫∞0e−iay(−iy+log(z))k+eiay(iy+log(z))kcosh(by)dy=2b−k−1eaπ2bπk+1Φ(−eaπb,−k,12+blog(z)π) | (4.1) |
where a, b, k and z are general complex numbers.
To obtain the first contour integral in the last line of Eq (3.6) we use the Cauchy formula by replacing y by log(z) and multiplying both sides by π2b and simplifying we get
πlogk(z)2bk!=14bi∫Czww−1−kdw | (4.2) |
We can write down an equivalent formula for the corresponding Lerch function for the second integral using Eqs (2.4), (3.6), and (4.2)
∫∞0eiay(iy+log(z))k−e−iay(−iy+log(z))ksinh(by)dy=2ieaπbπk+1b−k−1Φ(−eaπb,−k,1+blog(z)π)−πilogk(z)b | (4.3) |
a, b, k, and z are general complex numbers.
.
In this section we will derive definite integrals in terms of the Lerch function which simplify to the Hypergeometric function by Eq (1.11.10) in [3].
Φ(z,1,v)=∞∑n=0znn+v=v−12F1(1,v,1+v;z) | (5.1) |
For an extensive list of special cases of the hypergeometric function see [7]. These can be derived using the contiguous relations in section 15.2 in [2].
Here we form an equation by taking the difference between (4.3) and that equation where a is replaced by −a. Then set k=−1 and z=eβ to get
∫∞0sin(ay)sinh(by)dyβ2+y2=12β(eaπbΦ(−eaπb,1,1+bβπ)−e−aπbΦ(−e−aπb,1,1+bβπ)) | (5.2) |
Here we form an equation by taking the sum of (4.3) and that equation where a is replaced by −a. Then set k=−1 and z=eβ to get
∫∞0ycos(ay)sinh(by)dyβ2+y2=12(πbβ−e−aπbΦ(−e−aπb,1,1+bβπ)−eaπbΦ(−eaπb,1,1+bβπ)) | (5.3) |
We then set k=−1, z=eβ by forming an equation using (4.1) by taking the difference between (4.1) and that equation where a is replaced by −a to get
∫∞0cos(ay)cosh(by)dyβ2+y2=12a(eβπ2bΦ(−eβπb,1,12+bβπ)+e−aπ2bΦ(−e−βπb,1,12+bβπ)) | (5.4) |
Here we form an equation using (4.1) by taking the sum between (4.1) and that equation where a is replaced by −a.
∫∞0ysin(ay)cosh(by)dyβ2+y2=12(eβπ2bΦ(−eβπb,1,12+bβπ)−e−aπ2bΦ(−e−βπb,1,12+bβπ)) | (5.5) |
f(x) | ∫∞0f(x)dx |
sin(ay)sinh(πy)dyy2+β2 | 12β(eaΦ(−ea,1,1+β)−e−aΦ(−e−a,1,1+β)) |
sin(ay)sinh(πy)dyy2+1 | 12(eaΦ(−ea,1,2)−e−aΦ(−e−a,1,2)) |
sin(ay)sinh(π2y)dyy2+1 | 12(e2aΦ(−e2a,1,32)−e−2aΦ(−e−2a,1,32)) |
sin(ay)sinh(π4y)dyy2+1 | 12(e4aΦ(−e4a,1,54)−e−4aΦ(−e−4a,1,54)) |
ysin(ay)cosh(π4y)dy1+y2 | 12(e4aΦ(−e4a,1,34)−e−4aΦ(−e−4a,1,34)) |
ycos(ay)sinh(πy)dyy2+1 | 12(−e−aΦ(−e−a,1,2)−eaΦ(−ea,1,2)+1) |
ycos(ay)sinh(π2y)dyy2+1 | 12(−e−2aΦ(−e−2a,1,32)−e2aΦ(−e2a,1,32)+2) |
cos(ay)cosh(πy)dyβ2+y2 | 12β(ea/2Φ(−ea,1,β+12)+e−a/2Φ(−e−a,1,β+12)) |
cos(ay)cosh(πy)dy1+y2 | 12(e−a/2Φ(−e−a,1,32)+ea/2Φ(−ea,1,32)) |
cos(ay)cosh(π2y)dy1+y2 | 12(e−aΦ(−e−2a,1,1)+eaΦ(−e2a,1,1)) |
cos(ay)cosh(π4y)dy1+y2 | 12(e−2aΦ(−e−4a,1,34)+e2aΦ(−e4a,1,34)) |
Table 4.113 in [5] is correct with the exception of entry (4.113.1).
Note that the second parameter of the Lerch function in every entry in the table is unity so we can write it in terms of the hypergeometric function using (5.1). Many of these can be evaluated in terms of simpler functions from the extensive tables in [7].
In comparing our results with Table 4.113 our formulae have a wider range of the parameters than are listed in the Gradshteyn and Ryzhik book [5] due to the use of the Lerch function in the derivation of these integrals. We also provided correct formula for an integral supplied by Erdélyi et al [4]. We will be looking at other integrals using this contour integral method for future work. The results presented were numerically verified for both real and imaginary values of the parameters in the integrals using Mathematica by Wolfram.
This paper is fully supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) Grant No. 504070.
The authors declare no conflict of interest in this paper.
[1] | F. Oberhettinger, Tables of Fourier Transforms and Fourier Transforms of Distributions, 1st ed., Springer-Verlag, Berlin Heidelberg, 1990. |
[2] | M. Abramowitz, I. A. Stegun (Eds), Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing, New York, Dover, 1982. |
[3] | A. Erdéyli, W. Magnus, F. Oberhettinger, F. G. Tricomi, Higher Transcendental Functions, Vol. I, McGraw-Hill Book Company, Inc., New York-Toronto-London, 1953. |
[4] | A. Erdélyi, W. Magnus, F. Oberhettinger, F. G. Tricomi, Tables of Integral Transforms, Vol. I, McGraw-Hill Book Company, Inc., New York-Toronto-London, 1954. |
[5] | I. S. Gradshteyn, I. M. Ryzhik, Tables of Integrals, Series and Products, 6 Ed, Academic Press, USA, 2000. |
[6] |
R. Reynolds, A. Stauffer, A Method for Evaluating Definite Integrals in Terms of Special Functions with Examples, International Mathematical Forum, 15, 2020, 235–244. doi: 10.12988/imf.2020.91272
![]() |
[7] | Wolfram Research. The Wolfram Functions Site, Gauss Hypergeometric function 2F1, 1998–2020 Wolfram Research, Inc. Available from: https://functions.wolfram.com/HypergeometricFunctions/Hypergeometric2F1/03/. |