Research article

Note on an integral by Anatolii Prudnikov

  • Received: 17 November 2020 Accepted: 30 December 2020 Published: 04 January 2021
  • MSC : 01A55, 11M06, 11M35, 30-02, 30D10, 30D30, 30E20

  • Closed expressions for the integral

    0xm1logk(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 [4]. Some special cases of the integral are derived and discussed.

    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

    Related Papers:

    [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

    0xm1logk(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 [4]. Some special cases of the integral are derived and discussed.



    In 1986 Prudnikov et al [4] published a famous book containing a significant number of integrals. In this manuscript we focus on the integral

    0xm1logk(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πiCewywk+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 xm1(x2u+1)(x3u+1), then taking the definite integral with respect x[0,1] to get

    0xm1logk(ax)(x2u+1)(x3u+1)k!dx=12πi0Cawwk1xα1(x2u+1)(x3u+1)dαdx=12πiC0awwk1xα1(x2u+1)(x3u+1)dxdα=12πiCπawwk16u(4sin(π(u4α)6u)+32sin(πα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πiC2iπwk1exp(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πiy=0C2iπwk1exp(wlog(a)+iπα(2y+1)u)dα6u=12πiCy=02iπwk1exp(wlog(a)+iπα(2y+1)u)dα6u=12πiCπawwk1csc(πα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 (yiπ/(2u)) and multiplying by iexp(iπ(2m+u)/(4u)) to get the second equation and subtracting to get

    eiπ(2m+u)4u(yiπ2u)keiπ(2m+u)4u(y+iπ2u)kk!=12πiC2iwk1ewysin(π(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!(eiπ(2m+u)4u(log(a)+iπ(2y+1)uiπ2u)keiπ(2m+u)4u(log(a)+iπ(2y+1)u+iπ2u)k)=12πiy=0Ci2πwk1usin(π(u+2α)4u)awexp(iπw(2y+1)u+2iπmyu+iπmu)dα=12πiCy=0i2πwk1usin(π(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

    2k12πk+1(iu)kuk!(eiπ(2mu)4uΦ(e2imπu,k,π2iulog(a)4π)eiπ(6m+u)4uΦ(e2imπu,k,34iulog(a)2π))=12πiCπawwk12usin(π(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π(u4m)6u) for the first equation and then by y2iπ3u and multiply by exp(iπ(u4m)6u) for the second equation. Then we subtract the two to get

    eiπ(u4m)6u(y2iπ3u)keiπ(u4m)6u(y+2iπ3u)kk!=12πiC2iwk1ewysin(π(u4α)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π(u4m)6u(iulog(a)π+2y+13)kπkeiπ(u4m)6u(iu)k(iulog(a)π+2y+53)k)=12πiC4πwk1sin(π(u4α)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!(eiπ6eiπm3uΦ(e2imπu,k,56iulog(a)2π)eiπ6eiπm3uΦ(e2imπu,k,π3iulog(a)6π))=112uπiy=0C4πwk1sin(π(u4α)6u)awexp(iπw(2y+1)u+iπm(2y+1)u)dα=112uπiCy=04πwk1sin(π(u4α)6u)awexp(iπw(2y+1)u+iπm(2y+1)u)dα=12πC2πawwk13usin(π(u4α)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

    0xm1logk(ax)(x2u+1)(x3u+1)dx=2k12πk+1(iu)keiπ(2mu)4uu(Φ(e2imπu,k,π2iulog(a)4π)ieiπmuΦ(e2imπu,k,34iulog(a)2π))+eiπ/6(2π)k+1(iu)keiπm3u3u(e4iπm3uΦ(e2imπu,k,56iulog(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,v0,1,.. and is continued analytically by its integral representation given by

    Φ(z,s,v)=1Γ(s)0ts1evt1zetdt=1Γ(s)0ts1e(v1)tetzdt (4.3)

    where Re(v)>0, or |z|1,z1,Re(s)>0, or z=1,Re(s)>1.

    Using Eq (4.1) and setting k=0 and simplifying, we get

    0xm1(x2u+1)(x3u+1)dx=πcsc(πmu)6u(23sin(2πm3u)+32sin(π(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 u1, 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 u1, we get

    0xlog(x)log(xa)(x2+1)(x3+1)dx=π2((1603243)π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)+53sin(π6u)3sin(5π6u)+cos(5π6u))) (8.1)

    Using L'Hopital's rule as u1/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π233 (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

    0xm21log(log(x))(x2m+1)(x3m+1)dx=π6m(32log(35imΓ(78)Γ(58)12πΓ(38)Γ(18))+log(900iπΓ(34)2Γ(512)2Γ(112)25929mΓ(1112)2Γ(712)2Γ(14)2)4i3coth1(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π(4i3coth1(3)+32log(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

    01x(x2+1)(x3+1)(log2(x)+π2)dx=1+22+log(642)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(432π+2log(526)) (10.3)

    and

    01x(x2+1)(x3+1)(4log2(x)+π2)dx=48+92π+22log(8(99+702))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 m1/2 simplify the real and imaginary parts, we get

    0log(x)(x3/2+x3+x2+x)(log2(x)+π2)dx=43π4log(6) (10.5)

    and

    0log2(x)(x3/2+x3+x2+x)(log2(x)+π2)dx=π212π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(x4x)log(x)(x2+1)(x3+1)(log2(x)+π2)dx=1482(2(425)π+4(3063+66+122+2+log(8)+cosh1(1351)12tanh1(112))) (10.7)

    and

    0x4x(x2+1)(x3+1)(log2(x)+π2)dx=52π248π(log(4096)+2log((32)4(1+2+2)68(2+3)2(2+21)6)+6(2+6+42sin(π8)213sin(12tan1(512)))) (10.8)

    Using Eq (4.1) and setting k=1, a=1, and u=1/2 simplifying, we get

    0xm1(x+1)(x3/2+1)log(x)dx=13(3itan1(eiπm)2tanh1(e2iπm3)+2tanh1(31e2iπm3)+4tanh1((1)2/3e2iπm3)3tanh1(eiπm)tanh1(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 m1/2 simplify to get

    04x1x3/4(x+1)(x3/2+1)log(x)dx=12log(29(7+43)) (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
  • This article has been cited by:

    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
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1451) PDF downloads(32) Cited by(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog