arxiv: v1 [math.ca] 20 Aug 2008

Transcription

1 A CLASS OF LOGARITHMIC INTEGRALS arxiv:8875v [mathca] Aug 8 LUIS A MEDINA AND VICTOR H MOLL Abstract We present a systematic study of integrals of the form Z I Q = Qxloglog x dx, where Q is a rational function Introduction The classical table of integrals by I S Gradshteyn and I M Ryzhik [8] contains very few evaluations of the form I Q := where Q is a rational function The example Qxloglog/xdx, x + loglog/xdx = π log written in its trigonometric version 3 x + loglog/xdx = π/ π/4 πγ 3 4 Γ 4 loglogtanxdx, is the subject of Vardi s remarkable paper [3] This example appears as 497 in [8] and also in the equivalent form logxdx πγ coshx = πlog Γ, 4 as 437 in the same table We present here a systematic study of the logarithmic integrals Throughout the paper we indicate whether Mathematica 6 is capable of evaluating the integrals considered For example, a direct symbolic evaluation gives as 5 x + loglog/xdx = π 4 log 4π 3 Γ 4 4 The reader should be aware that the question of whether a definite integral is computable by a symbolic language depends on the form in which the integrand is expressed For instance, Mathematica 6 is unable to evaluate the trigonometric version of given as the right-hand side of 3, Date: October 3, 8 99 Mathematics Subject Classification Primary 33E, M35 Key words and phrases Integrals, L-functions, zeta functions, polylogarithm, Lerch phi function

2 LUIS A MEDINA AND VICTOR H MOLL The idea exploited here, introduced by I Vardi in [3], is to associate to each function Q a gamma factor 6 Γ Q s := Qx log x s dx, so that, the integral is given by 7 I Q = Γ Q An explicit evaluation of I Q is achieved in the case where Qx is analytic at x = Starting with the expansion 8 Qx = a n x n, we associate an L-series 9 L Q s := The integral is now evaluated as n= n= a n n+ s I Q = γl Q +L Q, where γ is the Euler-Mascheroni constant The interesting story of this fundamental constant can be found in [] The identity is essentially Vardi s method for the evaluation of Naturally, to obtain an explicit evaluation of I Q, one needs to express L Q and L Q in terms of special functions We employ here the Riemann zeta function ζs = n= ns, s >, and its alternating form n ζ a s = n s = s ζs, s > The relations n= 3 ζn = n+ π n B n, n N := N {}, n! ζ n = n+ B n, n N, n ζ n = nn!ζn+ π n, n N, ζ = log π, where B n are the Bernoulli numbers, will be used to simplify the integrals discussed below A second function that is used in the evaluations described here is the polylogarithm function defined by 7 PolyLog[c, x] := n= x n n c

3 LOGARITHMIC INTEGRALS 3 and its derivative 8 PolyLog, [c,x] := d dc PolyLog[c,x] = n= logn n c xn A complete description of I Q is determined here in the case where Qx is a rational function The method of partial fractions shows that it suffices to consider three types of integrals: the first type is 9 P j := that gives the polynomial part of Q, the second type is R m,j a := that treats the real poles of Q, and the third type is C m,j a,b := x j loglog/xdx, x j loglog/x x+a m+ dx, x j loglog/x x dx, +ax+b m+ with a 4b < This last case deals with the non-real poles of Q Integrals of first type These are simple They are evaluated in as P j := x j loglog/xdx = γ +logj + j + Integrals of second type The special case R m, is evaluated first We introduce the polynomial m T m x := j A m+,j+ x j, with A m,j the Eulerian numbers given in 46 Then the integral 3 E m := is evaluated, for m >, as j= T m x loglog/x x+ m+ 4 E m = m ζ m+ m γ m + m log B m m The sequence E m is then used to produce a recurrence for R m, The initial condition 5 R, = log dx

4 4 LUIS A MEDINA AND VICTOR H MOLL is given in Example 4 Then, the values of R m, are obtained from 6 R m, = E m m b m b k m b m R m k,, where k= m 7 b k m = k j= j A m,j+ k This is described in Corollary 47 The evaluation of R m, a, for a, appears in Proposition 5 and Corollary 53: the value of R, a is given by 8 R, a := and, for m >, we have loglog/xdx x+a γ R m, a = a m +am γ a m+ m! a m m! = γlog+/a PolyLog, [, /a], m j= S m,jt j /a +/a j m S m,jpolylog, [ j, /a] j= Here S m,j are the signless Stirling numbers of the first kind defined by the expansion m 9 t m = S m,jt j, j= where t m = tt+t+ t+m is the Pochhammer symbol The function PolyLog, [c,x] is defined in 8 For j >, the value of R m,j a is now obtained from the recursion in Theorem 6, written here as r 3 R m, a = α j,r ar m r+j,j a, where j= 3 α j,r a := j r j a r This can be used for increasing values of the free parameter r, to obtain analytic expressions for R m,j a For instance, r = gives 3 R m, a = α, ar m, a+α, ar m, a, that determines R m, a in terms of R m, a and R m, a, that were previously computed The value r = gives 33 R m, a = α, ar m, a+α, ar m, a+α, ar m, a, that determines R m, a in terms of previously computed integrals This procedure determines all the integrals R m,j a

5 LOGARITHMIC INTEGRALS 5 Integrals of third type These are integrals where the corresponding quadratic factor has non-real zeros The expression 34 x +ax+b = x cx c = x rxcosθ +r is used to define 35 D m,j r,θ = x j loglog/x x rxcosθ +r dx m+ Naturally C m,j a,b = D m,j r,θ, we are simply emphasizing the polar representation of the poles The value D, r,θ is computed first Theorem 7 treats the case r =, with the value D,,θ = π [ ] Γ θ/π θ/πlogπ+log sinθ Γθ/π The case r is given in Theorem 7 as D, r,θ = γ rsinθ tan sinθ r cosθ + PolyLog, [,e iθ /r] PolyLog, [,e iθ /r] risinθ The next step is to compute D, r,θ This is described at the end of Section 7 The result is expressed in terms of the Lerch zeta function z k 36 Φz, s, a := k +a s, as k= D, r,θ = γ r log rcosθ + r sinθ γcotθ tan r cosθ ] e + [Φ,, iθ e Φ risinθ r,,,, iθ r,, The reader will find in [] information about this function The values of D m,j r,θ for m, j >, are determined by the recurrences 37 D m,j r,θ = rmsinθ and 38 D m,j r,θ = θ D m,j r,θ, mcosθ r D m,j r,θ+rmd m,j r,θ ThesefollowdirectlyfromthedefinitionofD m,j r,θ DetailsaregiveninSection7 Comment Integration by parts, shows that the integrals I Q in include those of the form 39 J Q = Qx dx logx

6 6 LUIS A MEDINA AND VICTOR H MOLL This class was originally studied by V Adamchik [] They were considered by Baxter, Tempereley and Ashley [3] in their work on the so-called Potts model for the triangular lattice In that model, the generating function has the form 4 P 3 t := 3 sinhπ tx sinh tx 3 x sinhπx coshtx V Adamchik determined analytic expressions for 39 in the case where the denominator of Q is a cyclotomic polynomial The expressions involve derivatives of the Hurwitz zeta function 4 ζz, q := n= n+q z Zhang Nan-Hue and K Williams [4], [5] used the Hurwitz zeta function to evaluate definite integrals similar to the class considered here Examples of integrals that involve ζz,q in the integrand are given in [5], [6] and [7] The evaluation of the integral I Q := The main tool Qxloglog/xdx, for a general function Qx, is achieved by introducing the function Γ Q s := The next result is elementary Lemma The integral I Q is given by 3 I Q = Γ Q Qxlog/x s dx Example The simplest case is Qx Here we obtain 4 where 5 Γs = loglog/xdx = Γ, log x s dx = t s e t dt, is the classical gamma function The reader will find in [4] the identity 6 Γ = γ, where γ is the Euler-Mascheroni constant This example appears as 49 in [8] Example Consider now the case Qx = x a, for a R Observe that 7 Γ x as = dx e a+t t s dt = Γs a+ s, for a > and s > Differentiate with respect to s at s = to produce 8 x a loglog/xdx = γ +loga+ a+

7 LOGARITHMIC INTEGRALS 7 Differentiating 7 n times with respect to s, yields 9 x a log n xloglog/xdx = n+ n! +a n+ log+a+γ H n, where H n is the n-th harmonic number Mathematica 6 is unable to evaluate 7 if both a and n are entered as parameters The same holds for 9 Note The expression 8, with a = m N, provides the evaluation of the integral P m in 9: P m := This appears as 4358 in [8] x m loglog/xdx = γ +logm+ m+ 3 The case where Q is analytic at x = In this section we consider the evaluation of the integral 3 I Q := where Q is admits an expansion 3 Qx = Qxloglog/xdx, a n x n n= The expression for I Q is expressed in terms of the associated L-function defined by 33 L Q s := n= a n n+ s The idea for the next lemma comes from [3] Lemma 3 The function Γ Q satisfies Γ Q s = ΓsL Q s Proof The linearity of Γ Q s in the Q-argument shows that 34 Γ Q s = a n Γ x ns n= The result now follows from the value of Γ x ns in 7 Theorem 3 Assume Q is given by 3 Then 35 I Q := Qxloglog/xdx = γl Q +L Q Proof Differentiate the expression for Γ Q in the previous lemma and use the result of Lemma The theorem reduces the evaluation of I Q to the evaluation of L Q and L Q The first series of examples come from prescribing the coefficients a n of Qx so that the L Q function is relatively simple

8 8 LUIS A MEDINA AND VICTOR H MOLL Example 3 Choose a n = /n+ Then x n 36 Qx = n+ = log x x Then 37 L Q s = n= n= = ζs+, n+ s+ where ζs is the classical Riemann zeta function Theorem 3 gives 38 log x x loglog/xdx = γπ 6 ζ Mathematica 6 is unable to check this directly, but the change of variables x = e t transforms 38 to 39 logtlog e t dt = γπ 6 ζ This is computable by Mathematica 6 The constant 3 ζ = n= logn n can be expressed in terms of the Glaisher constant 3 loga := ζ, by 3 ζ = π γ +logπ loga 6 This gives 33 log x x as an alternative form for 38 loglog/xdx = π 6 loga logπ Example 3 We now consider the alternating version of Example 3 and choose a n = n /n+ In this case n x n 34 Qx = = log+x, n+ x n= and the corresponding L-series is n 35 L Q s = n+ s+ = s ζs+ n= Theorem 3 and the evaluations 36 L Q = π and L Q = π log + ζ,

9 LOGARITHMIC INTEGRALS 9 now yield 37 log+x x This can also be expressed as 38 log+x x loglog/xdx = π log γ+ ζ loglog/xdx = π log4π loga As in the previous example, Mathematica 6 is unable to produce this evaluation, but it succeeds with the alternate version 39 logtlog+e t dt = π log γ+ ζ Example 33 Adding the results of the first two examples yields log x 3 x Their difference produces loglog/xdx = π log+γ ζ 3 logtlogtanhtdt = γπ ζ + π log This cannot be evaluated symbolically Example 34 This example generalizes Example 3 The integrand involves the polylogarithm function defined in 7 The choice a n = /n+ c produces the function 3 Qx = n= x n n+ c = x PolyLog[c,x], and the corresponding L-function is 33 L Q s = = ζs+c n+ s+c Then 35 gives 34 PolyLog[c, x] x n= loglog/xdx = γζc++ζ c+ Example 35 Choosing now a n = n /n+ c gives n x n 35 Qx = n+ c = x PolyLog[c, x] and 36 L Q s = We conclude that PolyLog[c, x] x n= n= n n+ s+c = s c ζs+c loglog/xdx = γ c c log ζc+ c ζ c+

10 LUIS A MEDINA AND VICTOR H MOLL In the special casewhere c is a negativeinteger, the function Q reduces to a rational function Details are provided in section 4 Example 36 Interesting integrands can be produced when a n that are periodic sequences For example, the choice 37 a n = πn n cos 3 gives the function 38 Qx = n= πn n cos 3 and the corresponding L-function 39 L Q s = Then 35 gives 33 n= log+x+x x cos πn 3 n s+ x n = log+x+x, x = 3 s ζs+ loglog/xdx = γπ π log3+ 3 ζ Example 37 This example presents a second periodic sequence The choice 33 a n = πn n cos 5 gives the function 33 Qx = n= n= πn n cos 5 x n = log ϕx+x, x where ϕ = 5 / The corresponding L-function is cos πn 5 L Q s = n s+ = PolyLog[s+,e πi/5 ]+PolyLog[s+,e πi/5 ] The values L Q = π /5 and 333 L Q = and Theorem 3 give the identity 334 log ϕx+x x n= logn n cos πn loglog/xdx = γπ n=, logn n cos πn Example 38 TheresultofTheorem3reducestheevaluationofacertainclassof integrals to the evaluation of the corresponding L-functions Many natural choices ofthefunctionqleadtoseriesthattheauthorsareunabletoevaluate Forexample, Qx = e x produces the identity 335 e x loglog/xdx = γe n= logn, n! 5

11 LOGARITHMIC INTEGRALS and we have been unable to procude an analytic expression for the series above The same is true for the series in Evaluation of integrals with real poles An expression for R m, We now turn to the evaluation of the integrals 4 R m,j a = x j loglog/x x+a m+ dx The method of partial fractions can then be used to produce explicit formulas for integrals of the type 4 I Q = Qx loglog/xdx where Q is a rational function with only real poles In this section we introduce a special family of polynomials T m x and produce an explicit analytic expression for 43 E m := T m x loglog/x x+ m+ These are then employed to evaluate R m, dx Definition 4 The Eulerian polynomials A m are defined by the generating function 44 x xexp[t x] = m= A m x tm m! Note 4 The Eulerian polynomials appear in many combinatorial problems The coefficients A m,j in 45 A m x = m A m,j x j j= are the Eulerian numbers They count the number of permutations of {,,,n} which show exactly j increases between adjacent elements, the first element always being counted as a jump The numbers A m,j have an explicit formula 46 A m,j = and a recurrence relation j m+ k j k m, k k= 47 A m,j = ja m,j +m j +A m,j that follows from 48 A m+ x = x x d dx A mx+m+xa m x,

12 LUIS A MEDINA AND VICTOR H MOLL with A x = The recurrence 48 follows directly from 44 and it immediatly implies that A m x is a polynomial of degree m The first few are A x =, A x = x, A x = x +x, A 3 x = x 3 +4x +x, A 4 x = x 4 +x 3 +x +x More information about these polynomials can be found in [9] We now present the relation between Eulerian polynomials and the polylogarithm function Lemma 43 Let m N The polynomial A m x satisfies 49 PolyLog[ m, x] = A m x x+ m+ Proof The identity 4 PolyLog[ m, x] = n+ n+ m x n+, n= shows that 4 PolyLog[ m, x] = ϑ m, +x where ϑ = x d dx The claim now follows by using 48 and an elementary induction We now employ Eulerian polynomials to evaluate an auxiliary family of integrals Proposition 44 Let m N Define 4 T m x := A m+ x x m = j A m+,j+ x j, j= and 43 E m := T m x loglog/x x+ m+ Then 44 E m = m ζ m+γ m + m logζ m, for m dx Proof This is a special case of Example 34 Note 45 The relation 4 gives 45 E m = m ζ m+ m+ γ m + m log B m, for m m

13 LOGARITHMIC INTEGRALS 3 For example, 46 E = γ + logπ log, E = 4 γ 4 log +3logA, 3 where A is the Glaisher constant defined in 3 The expression for E m in 45 is now used to provide a recurrence for the integrals R m,, where R m,j a is defined in We begin with an example that will provide an initial condition for the recurrence in Corollary 47 Example 4 The integral R, is given by 47 R, = loglog/xdx +x = log The choice a n = n in Theorem 3 produces Qx = /+x and L Q s = s ζs Passing to the limit as s and using R, = γl Q +L Q gives the result Theorem 46 The integrals E m in 45 satisfy 48 E m = where m k= b k mr m k,, m 49 b k m = k j=k and A m,j are the Eulerian numbers given in 46 Proof In the expression 4 E m = use 4 to obtain j A m,j+ k T m x loglog/x x+ m+ m 4 E m = j A m,j+ j= dx x j loglog/x x+ m+ dx Now write x = x+, expand the resulting binomial and reverse the order of summation to obtain the result Corollary 47 The integrals R m, satisfy the recurrence 4 43 R, = E m E m R m, = b m b k m b m R m k, k= Proof First observe that b m Indeed, m 44 b m = A m,j = A m j=

14 4 LUIS A MEDINA AND VICTOR H MOLL Using the recurrrence 48 we conclude that A m+ = m+a m Therefore b m = A m = m! Example 4 The previous result provides the values R, = γ +logπ log, R, = 3 9γ log+36loga+6logπ 4 R 3, = 3 7γ 8log+36logA+4logπ+7ζ3/π 4 These integrals are computable using Mathematica 6 5 An expression for R m, a In this section we present an analytic expression for 5 R m, a := loglog/x dx x+a m+ The result is given in terms of the polylogarithm function defined in 7 and the derivative PolyLog, [c,x] defined in 8 The evaluation employs the expression for R m, given the previous section Proposition 5 The integral R, a is given by 5 loglog/xdx R, a := = γlog+/a PolyLog, [, /a] x+a Proof The expansion 53 Qx = x+a = a n= n a n x n, n= produces the L-function n 54 L Q s = a n+ n+ s = PolyLog[s, /a] Theorem 3 gives the result The evaluation of R m, a for m employs the signless Stirling numbers of the first kind S m,j defined by the expansion m 55 t m = S m,jt j, where t m = tt+t+ t+m is the Pochhammer symbol j= Theorem 5 Let m N and a > The L-function associated to Qx = /x+a m+ is 56 L Q s = m a m S m,jpolylog[s j, /a] m! j=

15 LOGARITHMIC INTEGRALS 5 Proof The identity 57 is used in the expansion 58 to produce 59 β = k β k k k! x+a β = a β +x/a β = a β x+a β = a β k= k= β a k x k k /a k β k k!k + s Now choose β = m+ and use the elementary identity m+ k 5 = k + m, k! m! to write the L-function corresponding to Qx = /x+a m+ as 5 L Q s = a m+ /a k k + m m! k + s Finally use the expression 55 to write m 5 L Q s = a m+ /a k S m,j m! k + s j, j= k= k= and identify the series as a polylogarithm to produce the result Corollary 53 Let m N Then the integral 53 R m, a := is given by γ R m, a = a m +am γ a m+ m! a m m! m j= loglog/x dx x+a m+ S m,jt j /a +/a j m S m,jpolylog, [ j, /a], j= where T j is the polynomial defined in 4 Proof The result follows from Theorem 3 and the expression for polylogarithms given in 43 Example 5 The choice a = and m = gives 54 R, := loglog/xdx x+ = γ 6 PolyLog, [, ] Mathematica 6 is unable to compute these integral The expansion of the polylogarithm function gives the identity loglog/xdx 55 x+ = γ 6 + n logn n+ n=

16 6 LUIS A MEDINA AND VICTOR H MOLL 6 An algorithm for the evaluation of R m,j a In this section we present an analytic expression for 6 R m,j a := x j loglog/x x+a m+ dx These results provide the evaluation of integrals of the type 6 I Q = Qxloglog/xdx where Q has real poles The case of non-real poles is discussed in the next section Theorem 6 The integrals R m,j a satisfy the recurrence r 63 R m, a = α j,r ar m r+j,j a, for any r m Here j= 64 α j,r a := j r j a r Note 6 The recurrence 63 can now be used for increasing values of the free parameter r, to obtain analytic expressions for R m,j a For instance, r = gives 65 R m, a = α, ar m, a+α, ar m, a, that determines R m, a in terms of R m, a and R m, a, that were previously computed The proof of Theorem 6 employs a recurrence for the functions α j,r a that is established first Lemma 63 Let k,r N and α j,r a as in 64 Then 66 α,r a x+a k + α,rax x+a k+ + + α r,rax r x+a k+r = x+a k+r Proof Expand the identity = x+ar a r x r x+a Proof of Theorem 6 Multiply the relation in Lemma 63 by loglog/x and integrate over [,] Example 6 We now use the method described above to check that 67 R, := The recursion 63 gives x loglog/x x+ dx = log +γ logπ+log 68 R, = α, R, +α, R, Using the values α, = and α, = and the integrals 69 R, = log

17 LOGARITHMIC INTEGRALS 7 given in Example 4 and 6 R, = γ +logπ log, computed in Example 4, we obtain the result Example 6 The computation of R, can be obtained from the recurrence 6 R, = α, R, +α, R,, and the previously computed values 6 R, = γ 6 PolyLog, [, ], and 63 R, = γlog 3 PolyLog, [, ] It follows that R, := x loglog/x x+ dx = γ 3 γlog 3 +PolyLog, [, ] PolyLog, [, ] = γ 3 γlog 3 n= n logn n /n Example 63 We now illustrate the recurrence 63 to obtain the value 64 R 3, 5 = We first let m = 3 in 33 to obtain x loglog/x x R 3, 5 = α, 5R, 5+α, 5R, 5+α, 5R 3, 5 The integrals with second index are given in 53 by dx 66 R, 5 = γ 3 5 PolyLog, [, 5 ] and R 3, 5 = 9γ 8 75 PolyLog, [, 5 ] 5 PolyLog, [, 5 ] 375 PolyLog, [, 5 ] The next step is to put m = in 3 to obtain 67 R, 5 = α, 5R, 5+α, 5R, 5 The values R, 5 = γ 8 5 PolyLog, [, 5 ] 5 PolyLog, [, 5 ] and R, 5 is given in 66 These come from Corollary53 Equation 67 now gives R, 5 = γ 36 + PolyLog, [, 5 ] PolyLog, [, 5 ]

18 8 LUIS A MEDINA AND VICTOR H MOLL Finally we obtain R 3, 5 = γ 34 3 PolyLog, [, 5 ] + PolyLog, [, 5 ] 5 PolyLog, [, 5 ] from 65 Note 64 The integrals R m,j a are computable by Mathematica 6 for a =, but not for a 7 Evaluation of integrals with non-real poles The integrals C m,j a,b = D m,j r,θ We consider now the evaluation of integrals 7 C m,j a,b := x j x loglog/xdx, +ax+b m+ where a 4b <, so that the quadratic factor has non-real zeros This is written as 7 x +ax+b = x cx c = x rxcosθ +r and we write 73 D m,j r,θ = x j x rxcosθ +r loglog/xdx m+ Naturally C m,j a,b = D m,j r,θ, we are simply emphasizing the polar representation of the poles Plan of evaluation: the computation of D m,j r,θ can be reduced to the range m and j m + by dividing x j by x rxcosθ + m+, in case j m + The fact is that the recurrences 37 and 38 determine all the integrals D m,j r,θ from D, r,θ and D, r,θ This is illustrated with the four integrals D,j r,θ : j 3 Begin with 37 with m = j = This gives 74 D, r,θ = rsinθ θ D,r,θ and then 38 with m = j = gives 75 D, r,θ = cosθ r D,r,θ+rD, r,θ, and this determines D, r,θ Now use m =, j = in 37 to obtain 76 D, r,θ = rsinθ θ D,r,θ Finally, 38 with m = and j = 3 yields 77 D,3 r,θ = cosθ r D,r,θ+rD, r,θ,

19 LOGARITHMIC INTEGRALS 9 Dividing x by x rxcosθ + r expresses D, r,θ as a linear combination of D, r,θ and D, r,θ This process determines D,3 r,θ We compute first the integral D, r,θ This task is is divided into two cases according to whether r = or not Theorem 7 gives the result for the first case and Theorem 7 describes the case r The evaluation of the integrals D m,j r,θ are then obtained by using the recurrences 37 and 38 Calculation of D,,θ This is stated in the next theorem Theorem 7 Assume < θ < π Then D,,θ := = Proof Consider the function 78 Qx = with the classical expansion 79 Qx = loglog/x x xcosθ + dx [ π θ/πlogπ+log sinθ sinθ x xcosθ + sink +θx k k= The corresponding L-function is given by 7 L Q s = and its value at s = is given by 7 L Q = k= k= sink +θ k + s, sink +θ k + while the closed form of the derivative at s = is 7 Γ θ/π = π θ, Γθ/π ] L Q = sink +θ logk + k + k= = π Γθ/π θ log +γ +logπ Γ θ/π π This identity can be found in [], page 5, #3 The result now follows from Theorem 3 Calculation of D, r,θ in the case r This is stated in the theorem below

20 LUIS A MEDINA AND VICTOR H MOLL Theorem 7 Assume < θ < π and r Then D, r,θ := loglog/x x rxcosθ +r dx sinθ r cosθ = γ rsinθ tan + risinθ Proof Consider the function 73 Qx = PolyLog, [,e iθ /r] PolyLog, [,e iθ /r] r sinθ x rxcosθ +r = and the corresponding L-function sink +θ L Q s := r k k + s The identity k= k= sink +θ r k x k e ik+θ e ik+θ = ir k k + s k= = r PolyLog[s,e it /r] PolyLog[s,e it /r] i 74 PolyLog[,a] = log a, yields the value 75 L Q = r [ log e iθ /r+log e iθ /r ] i The value of L Q can be written as sinθ 76 L Q = rtan r cosθ This follows by checking that both expressions for L Q match at θ = and their derivatives with respect to θ match We now present several special cases of these evaluations Many of them appear in the table of integrals [8] Example 7 The value D,,θ in Theorem 7 appears as 4357 in [8] Example 7 Replacing θ by θ +π, we obtain the evaluation loglog/x x +xcosθ + dx = π [ θlogπ +log sinθ π This appears as 43 in [8] Γ/+θ/π Γ/ θ/π Example 73 The value θ = π/ provides D,, π loglog/x = +x dx = π 4 logπ + π Γ3/4 log Γ/4 This is the example discussed by Vardi in [3] ]

21 LOGARITHMIC INTEGRALS Example 74 The angle θ = π/3 yields the value D,, π loglog/x πlogπ = dx = 3 x+x 3 + π log Γ5/6 3 3 Γ/6 This appears as 4356 in [8] where the answer is written in the equivalent form [ ] loglog/x π 5 dx = x+x 3 6 logπ logγ/6 The equivalent form logxdx e x +e x = π [ ] logπ logγ/6 appears as 433 in [8] Example 75 The angle θ = π/3 provides an evaluation of 4355 in [8]: D,, π loglog/x πlogπ = dx = 3 +x+x 3 + π log Γ/3 3 3 Γ/3 The equivalent form 77 logxdx e x +e x + = π Γ/3 π log 3 Γ/3 appears incorrectly as 433 in [8] The correct result is obtained by replacing π by 3 π Example 76 The limit of D,,θ as θ π gives the evaluation of D,, π loglog/x = +x dx = logπ log γ This is 4353 of [8] Example 77 It is easy to choose an angle and produce an integral that cannot be evaluated by Mathematica 6 For example, θ = 3π/4 gives D,, 3π loglog/x = 4 + x+x dx = π logπ Γ5/8 +log 4 Γ3/8 Note 73 The evaluation of L Q yields the identity 78 dx x rxcosθ +r = rsinθ tan sinθ r cosθ Note 74 The polylogarithm terms appearing in the expression for D, r,θ can be written as in terms of the sum logk + 79 Ur,θ := sink +θ r k k + k= The authors are unable to express the function Ur,θ in terms of special functions with real arguments We now proceed to a systematic determination of the integrals D m,j r,θ for m, j > For that, we use the recurrences 37 and 38

22 LUIS A MEDINA AND VICTOR H MOLL Proposition 75 Assume < θ < π Then D,,θ := x loglog/x x xcosθ+ dx = logπ 4sin θ +π θcotθ+ 8sin θ ψθ/π+ψ θ/π + π Γ θ/π 4 csc θcotθlog Γθ/π Proof The result follows directly from 37 and Theorem 7 Note 76 For the case r the value of D, r,θ can be obtained by differentiating the expression for D, r,θ in Theorem 7 Particular cases of this result are stated next Example 78 The angle θ = π/ produces D,, π := x loglog/x x + dx = logπ ψ ψ 4 Here ψx = Γ x/γx is the polygamma function Using the values ψ = γ π 3 4 3log and ψ = γ + π 4 3log that appear in [8] as and respectively, we obtain 7 The equivalent version 7 x loglog/x x + dx = logπ log γ 4 appears as 4373 in [8] logxdx cosh x = logπ log γ, Example 79 The angle θ = π/3 gives the evaluation D,, π x loglog/x := 3 x x+ dx = logπ + πlogπ π Γ5/ log Γ/6 + 6 ψ ψ 6 Using elementary properties of the ψ function and the values 7 ψ = γ π log3 log, and 73 ψ 5 = γ + π log3 log,

23 LOGARITHMIC INTEGRALS 3 that appear in [], page, we obtain x loglog/x x x+ dx = γ 3 log + 7πlog 3 9 log3 3 + logπ 3 + 8πlogπ 9 3 4π 3 3 logγ πlog Example 7 The recurrences 37 and 38 yield x loglog/x x 3x+ dx = Γ/ πlogπ+5πlog Γ/ This can be written as x loglog/x 3γ x dx = 3 3x ψ/+ψ/+ 8π ψ / ψ / 3 3log πlog 9 9 3log3+ 4 log 3 5πlog π ψ /+ 8π ψ / 55 logπ + πlogπ πlogγ/ 6 Example 7 The values r = and θ = π/3 yields the evaluation loglog/xdx x x+4 = γπ 6 3 i PolyLog,[ ], +i 3 4 PolyLog,[ ], i Mathematica 6 is unable to evaluate this integral Calculation of D, r,θ The integral 74 D, r,θ = corresponds to 75 Qx = xloglog/xdx x rxcosθ+r x x rxcosθ+r To evaluate the integral D, r,θ we employ the expansion 76 x x rxcosθ +r = sinkθ r k+ sinθ xk k= We conclude that the L-function associated to this Q is 77 L Q s = sinkθ sinθ r k+ k + s Therefore 78 L Q = sinθ k= k= sink +θ r k+ k +

24 4 LUIS A MEDINA AND VICTOR H MOLL To evaluate this sum, integrate 76 from to to produce 79 L Q = xdx x rxcosθ +r Observe that xdx x rxcosθ+r = x rcosθdx x rxcosθ +r + dx + rcosθ x rxcosθ +r Both integrals are elementary, the latter is given in 78 Therefore, 73 L Q = r log rcosθ + sinθ r +cotθ tan r cosθ The L-series 77 can be expressed in terms of the Lerch Φ-function 73 Φz, s, a := Indeed, L Q s = = = isinθ irsinθ risinθ k= z k k +a s e ikθ e ikθ r k+ k + s k= e iθ /r k k + s e iθ /r k k + s k= [ ] e iθ e iθ Φ r,s, Φ r,s, The next statement gives the value of D, r,θ Theorem 77 The integral 73 D, r,θ = xloglog/xdx x xrcosθ+r is given by D, r,θ = γ r log rcosθ + r sinθ γcotθ tan r cosθ ] e + [Φ,, iθ e Φ risinθ r,,,, iθ r,, Note 78 This evaluation completes the algorithm to evaluate all the integrals D m,j r,θ

25 LOGARITHMIC INTEGRALS 5 8 Adamchik s integrals V Adamchik presented in [] a series of beautiful evaluations of integrals of the form, where the denominator has the form + x n m for n, m N The results are expressed in terms of the Hurwitz zeta function and its derivatives For example, in Proposition 3, it is shown that x p +x γ +logn loglog/xdx = n n ζ, + n followed by Proposition 4 that states that x p x p ψ n p n γ +logn loglog/xdx = xn n + n ζ, p n ψ n+p n ζ, n+p n p ψ ψ n, p+ n ζ, p+ n The expressions become more complicated as the exponent of the denominator increases For instance, Proposition 5 gives x p n p p n+p +x n loglog/xdx = n logn+γ ψ ψ n n Γ p n γ +logn log n Γ n+p n and in Proposition 6 we find + n p n ζ, x p 3n p +x n loglog/xdx = 3 n log p n Γ p n ζ, n+p n, Γ n+p 5n plogn+γ 8n n + n pn plogn+γ p n+p 4n ψ ψ n n + ζ p, ζ, n+p n n n + n pn p 4n 3 ζ p, ζ, n+p n n We now describe some examples on how to use Theorem 3 to obtain some of the specific examples in [] Example 8 From Proposition 3 it follows that 8 x n logn loglog/xdx = log +xn n This appears as 7 in [] To check this evaluation, observe that 8 Qx = xn +x n = k x kn, k=

26 6 LUIS A MEDINA AND VICTOR H MOLL so the corresponding L-function is k 83 L Q s = k s n s k= A direct calculation shows that L Q = log n and 84 L Q = γlog n Then 8 follows from 35 log n = s ζs n s log logn n Example 8 Formula 8 in [] is also obtained from Proposition 3 and it states that x n 85 +x n loglog/xdx = log +log logn γ n To prove this, consider 86 Qx = xn +x n = k x +kn, whose L-function is 87 L Q s = k= k= k +k s n s = n s s ζs Replacing the values L Q = log and L n Q = log γlog logn+loglogn, n in 35 we obtain the result Example 83 The identity 88 x +x 4 loglog/xdx = π 4 log πγ 3 4 Γ 4 appears as formula 3 in [] To establish it, consider the function 89 Qx = x +x 4 = k x 4k+, with L-function 8 L Q s = k= k= k 4k + s = 3s ζs, 4 ζs, 3 4 The result follows from Theorem 3 by using the values 8 L Q = 8 ψ 3 4 ψ 4, and 8 L Q = 3log 8 The special values ψ 3 4 ψ 4, + 8 ζ, 4 ζ, ψ 4 = γ π 3log and ψ3 4 = γ + π 3log,

27 LOGARITHMIC INTEGRALS 7 and the value 84 ζ, 4 ζ, 3 4 = π γ +log+3logπ 4logΓ 4, are used to simplified the result This last expression comes from ζ, p q ζ, p πp q = πcot logπq +γ q q π log Γ j q sin πjp q This identity, established in [], follows directly from the classical Rademacher formula q 85 ζ z, p q = Γ zπq z πz sin + jpz q ζ z, j q An alternative evaluation of this integral comes from the partial fraction decomposition 86 x x 4 + a x + = a j= j= x ax+ ax +ax+ We assume a < and write a = cosθ Then 86 yields 87 x loglog/x x 4 + 4cos θx + dx = 4cosθ D,,θ D,,π+θ Using the result of Theorem 7 we obtain 88 x loglog/xdx x 4 + 4cos θx + = 4π 3 log sinθ The special case θ = π/4 produces 88 π 4sinθ θ π logπ log [ θ Γ Γ π + θ ] π Example 84 The case θ = π/ in the previous example reduces to Example 78 Example 85 The angle θ = π/3 in 88 yields 89 x loglog/x x 4 +x + dx = π 6log 3log3+8logπ logγ 3 3 Example 86 The angle θ = π/8 in 88 yields 8 x loglog/x x 4 x + dx = π 8 7logπ 4logsin π 8 8logΓ 8 9 A hyperbolic example The method introduced here can be used to provide an analytic expression for the family 9 LC n := logtdt cosh n+ t

28 8 LUIS A MEDINA AND VICTOR H MOLL The table of integrals [8] contains 9 LC = as formula 437 and 93 LC = as 4373 logtdt cosht = π log+3logπ 4logΓ 4, logtdt cosh = γ +logπ log, t The change of variables x = e t shows that 94 LC n = n+ x n loglog/x x dx, + n+ that identifies this integral as 95 LC n = n+ D n,n, π The recurrence 37, for r = and j = m, become 96 D m,m,θ = msinθ θ D m,m,θ, and the initial condition D,,θ = π sinθ [ θ/πlogπ+log ] Γ θ/π, Γθ/π provides a systematic procedure to compute LC n For instance, it follows that and LC = 3 D,, π = π Catalan+ π 4 LC 3 = 4 D 3,3, π = γ 3 4log 3 log+3logπ 4logΓ 4, + logπ 3 The Catalan constant appearing above is defined by 97 Catalan = n= n n ζ Small sample of a new type of evaluations We have introduced here a systematic method to deal with integrals of the form I Q = Qxloglog/xdx

29 LOGARITHMIC INTEGRALS 9 Extensions of this technique provides examples such as π +logxlogx+ loglog/xdx = +γ π + log ζ, +logx tan x loglog/xdx = γ π 48 π 4 + log tanh x x Details will presented elsewhere loglog/xdx = γπ 4 + π log ζ Conclusions We have developed an algorithm to evaluate integrals of the form I Q = Qx loglog/xdx + ζ 8, In the case where Qx is analytic at x =, with power series expansion Qx = a n x n we associate to Q its L-function 3 L Q s = Then n= n= a n n+ s 4 I Q = γl Q +L Q In the case Qx is a rational function, we provide explicit expressions for I Q in terms of special values of the logarithm, the Riemann zeta function, the polylogarithm PolyLog[c, x], its first derivative with respect to c and the Lerch Φ-function Acknowledgments The second author was partially funded by NSF-DMS 7567 The first author was supported by VIGRE grant NSF-DMS as a graduate student References [] V Adamchik A class of logarithmic integrals Proceedings ISSAC, 8, 997 [] A Apelblat Tables of Integrals and Series Verlag Harry Deutsch, Thun; Frankfurt am Main, 996 [3] R J Baxter, H N V Temperley, and S E Ashley Triangular potts model at its transition temperature, and related models Proc Royal Soc London, A 358: , 978 [4] G Boros and V Moll Irresistible Integrals Cambridge University Press, New York, st edition, 4 [5] O Espinosa and V Moll On some definite integrals involving the Hurwitz zeta function Part The Ramanujan Journal, 6:59 88, [6] O Espinosa and V Moll On some definite integrals involving the Hurwitz zeta function Part The Ramanujan Journal, 6: ,

30 3 LUIS A MEDINA AND VICTOR H MOLL [7] OEspinosa and VMollThe evaluation of Tornheimdouble sumspartjournal of Number Theory, 6: 9, 6 [8] I S Gradshteyn and I M Ryzhik Table of Integrals, Series, and Products Edited by A Jeffrey and D Zwillinger Academic Press, New York, 7th edition, 7 [9] R Graham, D Knuth, and O Patashnik Concrete Mathematics Addison Wesley, Boston, nd edition, 994 [] J Havil Gamma: exploring Euler s constant Princeton University Press, st edition, 3 [] A Laurincikas and R Garunkstis The Lerch zeta function Springer-Verlag, st edition, 3 [] H M Srivastava and J Choi Series associated with the zeta and related functions Kluwer Academic Publishers, st edition, [3] I Vardi Integrals, an Introduction to Analytic Number Theory Amer Math Monthly, 95:38 35, 988 [4] Z Yue and KS Williams Applications of the Hurwitz zeta function to the evaluation of certain integrals Canad Math Bull, 36: , 993 [5] Z Yue and KS Williams Values of the Riemann zeta function and integrals involving logsinh θ and logsin θ Pac Jour Math, 68:7 89, 995 Department of Mathematics, Tulane University, New Orleans, LA 78 address: lmedina@mathtulaneedu Department of Mathematics, Tulane University, New Orleans, LA 78 address: vhm@mathtulaneedu