Closed expressions for the integral
∫∞0xm−1logk(ax)(x2u+1)(x3u+1)dx
are given where the variables a, k, m and u are general complex numbers. Some of these closed expressions are given in [
Citation: Robert Reynolds, Allan Stauffer. Note on an integral by Anatolii Prudnikov[J]. AIMS Mathematics, 2021, 6(3): 2680-2689. doi: 10.3934/math.2021162
[1] | Robert Reynolds, Allan Stauffer . Extended Prudnikov sum. AIMS Mathematics, 2022, 7(10): 18576-18586. doi: 10.3934/math.20221021 |
[2] | 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 |
[3] | Robert Reynolds, Allan Stauffer . Note on an integral by Fritz Oberhettinger. AIMS Mathematics, 2021, 6(1): 564-568. doi: 10.3934/math.2021034 |
[4] | Robert Reynolds, Allan Stauffer . Derivation of some integrals in Gradshteyn and Ryzhik. AIMS Mathematics, 2021, 6(2): 1816-1821. doi: 10.3934/math.2021109 |
[5] | Maliha Rashid, Amna Kalsoom, Maria Sager, Mustafa Inc, Dumitru Baleanu, Ali S. Alshomrani . Mellin transform for fractional integrals with general analytic kernel. AIMS Mathematics, 2022, 7(5): 9443-9462. doi: 10.3934/math.2022524 |
[6] | 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 |
[7] | Shahid Mubeen, Rana Safdar Ali, Iqra Nayab, Gauhar Rahman, Thabet Abdeljawad, Kottakkaran Sooppy Nisar . Integral transforms of an extended generalized multi-index Bessel function. AIMS Mathematics, 2020, 5(6): 7531-7547. doi: 10.3934/math.2020482 |
[8] | D. L. Suthar, A. M. Khan, A. Alaria, S. D. Purohit, J. Singh . Extended Bessel-Maitland function and its properties pertaining to integral transforms and fractional calculus. AIMS Mathematics, 2020, 5(2): 1400-1410. doi: 10.3934/math.2020096 |
[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] | Waleed Mohamed Abd-Elhameed, Ahad M. Al-Sady, Omar Mazen Alqubori, Ahmed Gamal Atta . Numerical treatment of the fractional Rayleigh-Stokes problem using some orthogonal combinations of Chebyshev polynomials. AIMS Mathematics, 2024, 9(9): 25457-25481. doi: 10.3934/math.20241243 |
Closed expressions for the integral
∫∞0xm−1logk(ax)(x2u+1)(x3u+1)dx
are given where the variables a, k, m and u are general complex numbers. Some of these closed expressions are given in [
In 1986 Prudnikov et al [4] published a famous book containing a significant number of integrals. In this manuscript we focus on the integral
∫∞0xm−1logk(ax)(x2u+1)(x3u+1)dx |
which has a closed form in in terms of the Lerch function. We appeal to our method to derive a closed form solution which is continued analytically along with some special cases. A consequence of using our method with this integral is the derivation of the Mellin transform of the product of the logarithmic function and a rational function which is not found in current literature. This work also looks to expand the Table 2.5.1 in [1] where similar and relevant formula are listed. In our case the constants in the formulas are general complex numbers subject to the restrictions given below. The derivations follow the method used by us in [5], [6], [7], [8], [9] and [10]. The generalized Cauchy's integral formula is given by
ykk!=12πi∫Cewywk+1dy. | (1.1) |
This method involves using a form of Eq (1.1) then multiplies both sides by a function, then takes 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 [9]. The variable of integration in the contour integral is α=m+w. The cut and contour are in the first quadrant of the complex α-plane with 0<ℜ(α)<5u. The cut approaches the origin from the interior of the first quadrant and the contour goes round the origin with zero radius and is on opposite sides of the cut. Using a generalization of Cauchy's integral formula, we first replace y by log(ax) followed by multiplying both sides by xm−1(x2u+1)(x3u+1), then taking the definite integral with respect x∈[0,1] to get
∫∞0xm−1logk(ax)(x2u+1)(x3u+1)k!dx=12πi∫∞0∫Caww−k−1xα−1(x2u+1)(x3u+1)dαdx=12πi∫C∫∞0aww−k−1xα−1(x2u+1)(x3u+1)dxdα=12πi∫Cπaww−k−16u(4sin(π(u−4α)6u)+3√2sin(πα2u+π4)+1)csc(παu)dα | (2.1) |
from Eq (2.2.7.7) in [4], where 0<ℜ(α)<5u. We are able to switch the order of integration over α and x using Fubini's theorem since the integrand is of bounded measure over the space C×[0,∞).
Using Eq (1.1) and replacing y by log(a)+iπ(2y+1)u, then multiplying both sides by −2iπ6uexp((2y+1)iπmu), we get
−2iπk+1(iu)ke(2y+1)iπmu(−iulog(a)π+2y+1)k6uk!=−12πi∫C2iπw−k−1exp(wlog(a)+iπα(2y+1)u)dα6u | (3.1) |
We then take the infinite sum over y∈[0,∞) to get
−i(2π)k+1(iu)keiπmuΦ(e2imπu,−k,π−iulog(a)2π)6uk!=−12πi∞∑y=0∫C2iπw−k−1exp(wlog(a)+iπα(2y+1)u)dα6u=−12πi∫C∞∑y=02iπw−k−1exp(wlog(a)+iπα(2y+1)u)dα6u=12πi∫Cπaww−k−1csc(παu)dα6u | (3.2) |
from (1.232.3) in [3], where csch(ix)=−icsc(x) from (4.5.10) in [2] and ℑ(α)>0.
Next using Eq (1.1), we obtain the left-hand side replacing y by (y+iπ/(2u)) and multiplying by iexp(iπ(2m+u)/(4u)) to get the first equation and by replacing y by (y−iπ/(2u)) and multiplying by iexp(−iπ(2m+u)/(4u)) to get the second equation and subtracting to get
−e−iπ(2m+u)4u(y−iπ2u)k−eiπ(2m+u)4u(y+iπ2u)kk!=12πi∫C2iw−k−1ewysin(π(u+2α)4u)dα | (3.3) |
Next using Eq (3.3), we replace y by log(a)+iπ(2y+1)u then multiply both sides by −121/2πuexp(iπm(2y+1)/u) take the infinite sum over y∈[0,∞) to get
∞∑y=0πexp(iπm(2y+1)u)√2uk!(e−iπ(2m+u)4u(log(a)+iπ(2y+1)u−iπ2u)k−eiπ(2m+u)4u(log(a)+iπ(2y+1)u+iπ2u)k)=−12πi∞∑y=0∫Ci√2πw−k−1usin(π(u+2α)4u)awexp(iπw(2y+1)u+2iπmyu+iπmu)dα=−12πi∫C∞∑y=0i√2πw−k−1usin(π(u+2α)4u)awexp(iπw(2y+1)u+2iπmyu+iπmu)dα | (3.4) |
Next we simplify to get the Lerch function contour integral representation given by
2k−12πk+1(iu)kuk!(eiπ(2m−u)4uΦ(e2imπu,−k,π−2iulog(a)4π)−e−iπ(6m+u)4uΦ(e2imπu,−k,34−iulog(a)2π))=12πi∫Cπaww−k−1√2usin(π(u+2α)4u)csc(παu)dα | (3.5) |
from (1.232.3) in [3], where csch(ix)=−icsc(x) from (4.5.10) in [2] and ℑ(α)>0.
Here we use Eq (1.1) and replacing y by y+2iπ3u and multiply by exp(−iπ(u−4m)6u) for the first equation and then by y−2iπ3u and multiply by exp(iπ(u−4m)6u) for the second equation. Then we subtract the two to get
eiπ(u−4m)6u(y−2iπ3u)k−e−iπ(u−4m)6u(y+2iπ3u)kk!=12πi∫C2iw−k−1ewysin(π(u−4α)6u)dα | (3.6) |
Next we replace y by log(a)+iπ(2y+1)u then multiply both sides by −2iπexp(iπm(2y+1)u) to get
−2iπexp(2iπmyu+iπmu)k!(πk(iu)keiπ(u−4m)6u(−iulog(a)π+2y+13)k−πke−iπ(u−4m)6u(iu)k(−iulog(a)π+2y+53)k)=12πi∫C4πw−k−1sin(π(u−4α)6u)exp(w(log(a)+iπ(2y+1)u)+2iπmyu+iπmu)dw | (3.7) |
Next we take the infinite sum over y∈[0,∞) and multiply by 1/6u to get
−(2π)k+1(iu)k3uk!(−e−iπ6e−iπm3uΦ(e2imπu,−k,56−iulog(a)2π)−e−iπ6eiπm3uΦ(e2imπu,−k,π−3iulog(a)6π))=112uπi∞∑y=0∫C4πw−k−1sin(π(u−4α)6u)awexp(iπw(2y+1)u+iπm(2y+1)u)dα=112uπi∫C∞∑y=04πw−k−1sin(π(u−4α)6u)awexp(iπw(2y+1)u+iπm(2y+1)u)dα=12π∫C2πaww−k−13usin(π(u−4α)6u)csc(παu)dα | (3.8) |
from (1.232.3) in [3], where csch(ix)=−icsc(x) from (4.5.10) in [3] and ℑ(α)>0.
Since the right-hand side of Eq (2.1) is equal to the sum of Eqs (3.2), (3.5) and (3.8), we may equate the left hand sides to get
∫∞0xm−1logk(ax)(x2u+1)(x3u+1)dx=2k−12πk+1(iu)keiπ(2m−u)4uu(Φ(e2imπu,−k,π−2iulog(a)4π)−ieiπmuΦ(e2imπu,−k,34−iulog(a)2π))+e−iπ/6(2π)k+1(iu)keiπm3u3u(e4iπm3uΦ(e2imπu,−k,56−iulog(a)2π)−eπi/3Φ(e2imπu,−k,π−3iulog(a)6π))−132kπk+1(iu)k+1eiπmuΦ(e2imπu,−k,π−iulog(a)2π) | (4.1) |
The Lerch function has a series representation given by
Φ(z,s,v)=∞∑n=0(v+n)−szn | (4.2) |
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 | (4.3) |
where Re(v)>0, or |z|≤1,z≠1,Re(s)>0, or z=1,Re(s)>1.
Using Eq (4.1) and setting k=0 and simplifying, we get
∫∞0xm−1(x2u+1)(x3u+1)dx=πcsc(πmu)6u(−2√3sin(2πm3u)+3√2sin(π(2m+u)4u)+2cos(2πm3u)+1) | (5.1) |
Using Eq (4.1), then taking the first partial derivative with respect to m, then set m=1, k=1 and a=1/a on the right-hand side followed by using L'Hopital's rule on the right-hand side as u→1, we get
∫∞0log(x)log(xa)(x2+1)(x3+1)dx=1432π2(27π−37log(1a)) | (6.1) |
Using Eq (4.1), then taking the first partial derivative with respect to m, then set m=2, k=1 and a=1/a followed by using L'Hopital's rule on the right-hand side as u→1, we get
∫∞0xlog(x)log(xa)(x2+1)(x3+1)dx=π2((160√3−243)π−45log(1a))3888 | (7.1) |
Using Eq (4.1) setting k=a=1 and m=1/2, we get
∫∞0log(x)√x(x2u+1)(x3u+1)dx=−π2csc2(π2u)72u2(20cos(π6u)+27cos(π4u)+12cos(π2u)+9cos(3π4u)+4(−9sin3(π4u)+5√3sin(π6u)−√3sin(5π6u)+cos(5π6u))) | (8.1) |
Using L'Hopital's rule as u→1/2 to the right-hand side of equation (8.1), we get
∫∞0log(x)x3+x2+x3/2+√xdx=−37π2108 | (8.2) |
We get
∫∞0log(x)√x(x2+1)(x3+1)dx=−2π23√3 | (8.3) |
Using Eq (4.1), we first set a=1 and replace u by 2m and m by m/2, then take the first partial derivative with respect to k, then set k=0 simplify to get
∫∞0xm2−1log(log(x))(x2m+1)(x3m+1)dx=−π6m(3√2log(−35imΓ(−78)Γ(−58)12πΓ(−38)Γ(−18))+log(900iπΓ(−34)2Γ(−512)2Γ(−112)25929mΓ(−1112)2Γ(−712)2Γ(−14)2)−4i√3coth−1(√3)) | (9.1) |
Next we set m=1 in Eq (9.1) simplify to get
∫∞0log(log(x))√x(x2+1)(x3+1)dx=−16π(−4i√3coth−1(√3)+3√2log(−35iΓ(−78)Γ(−58)12πΓ(−38)Γ(−18))+log(900iπΓ(−34)2Γ(−512)2Γ(−112)25929Γ(−1112)2Γ(−712)2Γ(−14)2)) | (9.2) |
Note there is a singularity at x=1 which is removable using the principal value of the integral.
In this section, we will list some interesting examples by evaluating Eq (4.1) for various values of the parameters k, a, m and u in terms of fundamental constants and trigonometric functions.
Using Eq (4.1) and setting k=a=−1, u=1 and m=1/2 and rationalizing the real and imaginary parts, we get
∫∞0log(x)√x(x2+1)(x3+1)(log2(x)+π2)dx=√32−π3 | (10.1) |
and
∫∞01√x(x2+1)(x3+1)(log2(x)+π2)dx=−1+2√2+log(6−4√2)2π | (10.2) |
Using Eq (4.1) and setting k=−1, a=i, u=1 and m=1/2 and rationalizing the real and imaginary parts, we get
∫∞0log(x)√x(x2+1)(x3+1)(4log2(x)+π2)dx=18(4√3−√2π+√2log(5−2√6)) | (10.3) |
and
∫∞01√x(x2+1)(x3+1)(4log2(x)+π2)dx=−48+9√2π+2√2log(8(99+70√2))48π | (10.4) |
Using Eq (4.1) and setting k=−1, a=−1, and u=1/2 followed by taking the first partial derivative with respect to m, then apply L'Hopital's rule as m→1/2 simplify the real and imaginary parts, we get
∫∞0log(x)(x3/2+x3+x2+√x)(log2(x)+π2)dx=4−3π4−log(6) | (10.5) |
and
∫∞0log2(x)(x3/2+x3+x2+√x)(log2(x)+π2)dx=π2−12πlog(2) | (10.6) |
Using Eq (4.1) and setting k=−1, a=−1, and u=1 and m=1/4 for the first equation, then replace m=1/2 for the second equation, followed by subtracting the two and rationalizing the real and imaginary parts, we get
∫∞0(√x−4√x)log(x)(x2+1)(x3+1)(log2(x)+π2)dx=148√2(2(4√2−5)π+4(−30−6√3+6√6+12√2+√2+log(8)+cosh−1(1351)−12tanh−1(√1−1√2))) | (10.7) |
and
∫∞0√x−4√x(x2+1)(x3+1)(log2(x)+π2)dx=5√2π−248π(log(4096)+√2log((√3−2)4(1+√2+√2)68(2+√3)2(√2+√2−1)6)+6(−2+√6+4√2sin(π8)−2√13sin(12tan−1(512)))) | (10.8) |
Using Eq (4.1) and setting k=−1, a=1, and u=1/2 simplifying, we get
∫∞0xm−1(x+1)(x3/2+1)log(x)dx=13(−3itan−1(eiπm)−2tanh−1(e2iπm3)+2tanh−1(3√−1e2iπm3)+4tanh−1((−1)2/3e2iπm3)−3tanh−1(eiπm)−tanh−1(e2iπm)) | (10.9) |
Then we form a second equation by replacing m by n then taking the difference. Then we set n=1/4 followed by using L'Hopital's rule as m→1/2 simplify to get
∫∞04√x−1x3/4(x+1)(x3/2+1)log(x)dx=12log(29(7+4√3)) | (10.10) |
In this article, we derived some interesting definite integrals in [4]. We found that we are able to achieve a wider range of computation using one formula as opposed to previous works. 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 research is supported by NSERC Canada under grant 504070.
All authors declare no conflicts of interest in this paper.
[1] | Yu. A. Brychkov, O. I. Marichev, N. V. Savischenko, Handbook of Mellin Transforms, CRC Press, Taylor & Francis Group, Boca Raton, Fl, 2019. |
[2] | M. Abramowitz, I. A. Stegun, (Eds), Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing, New York, Dover, 1982. |
[3] | I. S. Gradshteyn, I. M. Ryzhik, Tables of Integrals, Series and Products, 6 Ed, Academic Press, USA, 2000. |
[4] | A. P. Prudnikov, Yu. A. Brychkov, O. I. Marichev, Integrals and Series, More Special Functions, USSR Academy of Sciences, Vol. 1, Moscow, 1990. |
[5] |
R. Reynolds, A. Stauffer, Definite Integral of Arctangent and Polylogarithmic Functions Expressed as a Series, Mathematics, 7 (2019), 1099. doi: 10.3390/math7111099
![]() |
[6] |
R. Reynolds, A. Stauffer, A Definite Integral Involving the Logarithmic Function in Terms of the Lerch Function, Mathematics, 7 (2019), 1148. doi: 10.3390/math7121148
![]() |
[7] |
R. Reynolds, A. Stauffer, Derivation of Logarithmic and Logarithmic Hyperbolic Tangent Integrals Expressed in Terms of Special Functions, Mathematics, 8 (2020), 687. doi: 10.3390/math8050687
![]() |
[8] |
R. Reynolds, A. Stauffer, Definite integrals involving product of logarithmic functions and logarithm of square root functions expressed in terms of special functions, AIMS Mathematics, 5 (2020), 5724–5733. doi: 10.3934/math.2020367
![]() |
[9] |
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
![]() |
[10] |
R. Reynolds, A. Stauffer, Integrals in Gradshteyn and Ryzhik: hyperbolic and algebraic functions, International Mathematical Forum, 15 (2020), 255–263. doi: 10.12988/imf.2020.91279
![]() |
1. | Robert Reynolds, Allan Stauffer, Barbara Martinucci, The Mellin Transform of Logarithmic and Rational Quotient Function in terms of the Lerch Function, 2021, 2021, 2314-4785, 1, 10.1155/2021/8814756 |