arXiv:2102.04891v5 [math.FA] 13 Aug 2021 µ comprehensive truly his In Introduction 1 lands; nteGmafnto,wt plctost h eafunctio Beta the to applications w with from function, pro Weierstrass’s Gamma integral after the Beta perhaps on and Euler’s hand of other the properties on Today, many p the 123-343 knowledge were the time pp. reviewed has numbers Binet page contributions, correct beautiful the but article, ordinary ie a eie w nerlrpeettoso i funct his of representations ( Re integral two derived has Binet ( ∗ 1 z aut fEetia niern,MteaisadComput and Mathematics Engineering, Electrical of Faculty h eeec nWitkradWto 2,font,p 248] p. footnote, [28, Watson and Whittaker in reference The nrlto oteGmafnto ( Γ function Gamma the to relation in ) z ) email xedBntsmto ofcoilsre fLpaetransforms. Laplace of series factorial to method Binet’s extend xml ht ihasm ubro em vlae,teBntgen Binet the evaluated, of terms value of optimized number same an a with that, example ucinadapytoesre otedgmaadplgmafunct polygamma and Stirling’s f digamma with the series series to factorial factorial series generalized corresponding those the apply and compute function We presented. are o h ie function Binet the for fe re eiwo h ie function Binet the of review brief a After > ie’ atra eisadetnin oLpaetransfo Laplace to extensions and series factorial Binet’s eivsiaeagnrlzto fBntsfcoilsre ntepara the in series factorial Binet’s of generalization a investigate We 0, [email protected]. : µ 1 α mmi”[,Scin3 .23 n13,JcusBntdefines Binet Jacques 1839, in 223] p. 3, Section [2, “memoir” ( µ z a ettebs osbeacrc fSiln’ xaso.Fnly w Finally, expansion. Stirling’s of accuracy possible best the beat can ( o ( Γ log = ) z o ( Γ log = ) ef nvriyo Technology of University Delft µ µ ( ( z z z = ) itVnMieghem Van Piet z ) 2 = ) ) 1Ags 2021 August 11 − − z m µ X asymptotic as ) ∞ =1 ( Abstract z z Z z ,svrlpoete fteBntpolynomials Binet the of properties several ), 0 − Q − ∞ bu h eafnto eoe13.Tefclpita that at point focal The 1839. before Beta-function the about 1 m k 1 2 1 2 =0 utadHne’ otu nerl h hoyconcentrate theory the integral, contour Hankel’s and duct arctan − e n. 1 itt oksz raie pr rmadn i own his adding from Apart treatise. book-size a to oint b 2 ( log m log πt z rSine . o 01 60G ef,TeNether- The Delft, GA 2600 5031, Box P.O Science, er ihpoete fteGmafnto eederived. were function Gamma the of properties hich xaso n eosrt yanumerical a by demonstrate and expansion ( + α z − z ) α + ihpg ubr rm13t 4 ugssan suggests 143 to 123 from numbers page with + 1 o 2,p 4-5] 1,p 1-1] for 211-217], p. [18, 248-251], p. [28, ion z t + z ∗ z dt − k − ) 1 2 1 2 o (2 log rtedrvtvso h Binet the of derivatives the or o (2 log rlzdfcoilsre with series factorial eralized π os ecmaeBinet’s compare We ions. meter ) π (1) ) α rms i function his b m ( α e ) (2) s and 1 e zt 1+ e t 2 µ (z)= ∞ − − dt (3) 2 t 1 e t − t Z0 − − There exist more representations2 of µ (z), but here we merely concentrate on elegant, converging factorial series for Re (z) > 0 due to Binet3 in [2, p. 234],
∞ β µ (z)= m (6) (z + 1) (z + 2) (z + m) m=1 X · · · where the coefficients are 1 1 1 β = u u (u + 1) (u + m 1) du (7) m m − 2 · · · − Z0 1 59 29 533 1577 280361 69311 Explicitly, β1 = β2 = 12 , β3 = 360 , β4 = 60 , β5 = 280 , β6 = 168 , β7 = 5040 , β8 = 180 , and βm is positive and rapidly increasing in m. Nearly at the end of his memoir and somewhat hidden, Binet [2, p. 342] has given a second factorial series (31), that is rederived slightly differently in Theorem 1 in Section 2.3 and generalized to Laplace transforms in Section 6.2. Recently, Nemes [13, Theorem 2.1] has generalized Binet’s expansion (6), for 0 a 1 and ≤ ≤ Re (z) > 0,
1 1 ∞ c (a) log Γ (z + a)= z + a log z z + log (2π)+ m (8) − 2 − 2 (z + 1) (z + 2) (z + m) mX=1 · · · where
m 1 1 1 1 1 a − cm (a)= u + a u (u + 1) (u + m 1) du (u a +1+ j) du (9) m 0 − 2 · · · − − m 0 − Z Z jY=0
Clearly, if a = 0, we retrieve Binet’s first factorial expansion (6) and coefficients βm in (7). Nemes’ expansion (8), expressed in terms of Binet’s function µ (z) with the definition (1),
1 z a ∞ cm (a) µ (z)= z log − + a + m − 2 z k=1(z a + k) mX=1 − bears resemblance to Q ∞ bm (α) µ (z)= m 1 (10) m=1 k=0− (z + α + k) X 2Blagouchine [3] lists 7 different formulae for log Γ (z).Q 3In Binet’s notation [2, p. 234], I (1) I (2) I (3) 2µ (z)= + + + · · · (4) z + 1 2(z +1)(z + 2) 3(z +1)(z +2)(z + 3) where the integral form of the coefficients (derived at [2, p. 238]) is
1 I (m)= x (x +1)(x + 2) · · · (x + m − 1) (2x − 1) dx (5) Z0
1 1 59 227 Binet [2, p. 234] lists the first few coefficients, I (1) = 6 , I (2) = 3 , I (3) = 60 and I (4) = 60 .
2 in our main Theorem 3, where the coefficients bm (α) are called Binet polynomials, to honour Jacques Binet. Much earlier, Hermite [11] has deduced the corresponding generalized Stirling asymptotic expansion
K k 1 1 1 ( 1) − B (a) log Γ (z + a) z + a log z z + log (2π)+ − k+1 (11) ≃ − 2 − 2 k (k + 1) zk Xk=1 in terms of the Bernoulli polynomials Bn (a) that reduces for a = 0 to Stirling’s original asymptotic series (23). Starting from a complex integral (20) for Binet’s function µ (z), Stirling’s asymptotic series (23) is derived in Section 2.1.1, where also the meaning of the upper bound K is explained. Section 2.1.2 presents the convergent companion (26) of Stirling’s divergent series (if K in (23)). →∞ The main motivation that led to this article is twofold. Originally, I was confused about Binet’s achievements: he derived two different factorial expansions (6) and (31) for the same function µ (z). While I thought initially that one of them must be wrong, I discovered later with (10) that infinitely many different factorial series exist. The second motivation was my unbelief that Stirling’s asymptotic but divergent expansion seems unbeatable in performance for some optimal, finite K. Our main result is the demonstration that there exist infinitely many different factorial expansions in a complex parameter α for Binet’s function µ (z) in Theorem 3 (Section 3) and for Laplace transforms in Theorem 4 (Section 6.3). For α = 1 in (10), we recover Binet’s first factorial series (6) with bm (1) = βm in (7), while α = 0 in (10) corresponds to Binet’s second factorial series (31) with bm (0) = cm in (32). Interestingly, the α = 1 factorial series has all positive coefficients and any truncation thus lower bounds Binet’s function µ (z), while the α = 0 factorial series is shown in
Theorem 2 to possess coefficients cm < 0 for m> 2. Thus, a truncation of (31) upper bounds µ (z). For 1 α around 10 , which is close to the largest zero of the Binet polynomial bm (α), numerical computations exhibit the fastest convergence. With the same number K of terms evaluated, the variant α = 0 is more accurate than the variant α = 1. Perhaps, the slower convergence of Binet’s expansion (6) has caused that his factorial series is not listed in handbooks of functions, like Abramowitz & Stegun [1] nor in its successor by Olver et al. [15]. We will first discuss the main properties of Binet’s function µ (z) in Section 2 and the deductions from the complex integral (20) in Section 2.1, before we review parts of Binet’s great treatise. In subsection 2.2, we sketch Binet’s route towards his first factorial series (6), that is covered in the literature (see e.g. [28, p. 253], [16, p. 30]). Binet’s second factorial series, that I have not found in later works, is derived in more detail in subsection 2.3, also because I believe that, being the case α = 0 in (10), it is slightly more important than his first series (6) corresponding to α = 1. Moreover, Binet’s method towards his second factorial series, which is a recipe in five steps, enables a far reaching generalization to Laplace transforms as explained in Section 6. Factorial expansions for the derivatives dnµ(z) dzn are derived in Section 4 and applied to the digamma and polygamma function. In particular for Γ′(z) the digamma function ψ (z)= Γ(z) , we thus add a convergent series (60) to its asymptotic counterpart (61). Section 5 discusses and compares, with a same number of terms, the accuracy of Stirling’s asymptotic expansion and the best possible that can be attained by the generalized Binet factorial expansion. The commonly accepted belief about the superiority of Stirling’s asymptotic expansion is demonstrated and plotted in Fig. 1. However, at the zeros of the Binet polynomial bm (α), the accuracy of the generalized Binet factorial expansion (10) improves considerably as drawn in Fig. 2.
3 By a numerical example, we show that Stirling’s series accuracy is not always better than a factorial series (with a same number of terms)! In other words, the generalized Binet expansion (10) can be optimized with respect to the “free” parameter α to achieve, at least, a comparable accuracy with a same computational effort. Perhaps, this observation deserves to list the generalized Binet factorial expansion (10) in handbooks of functions. Computations are deferred to the appendices in order to enhance the readability and focus on the essential parts.
2 Binet’s function µ (z)
The definition (1) of Binet’s function µ (z) directly shows, for a positive integer z = n, that 1 1 µ (n)= n + log (n 1)! n log (n) log (2π) − − − 2 − 2 The sequence µ (1) 0.0811, µ (2) 0.0413, µ (3) 0.0277, µ (4) 0.0208, µ (5) 0.0166,. . . , ≃ ≃ ≃ ≃ ≃ µ (10) 0.0083 demonstrates the slow decay roughly as µ (n) 1 . The precise decay is given in ≃ ≈ 12n (65) below. Maximum at real, positive values of z. Binet’s integral (3) can be rewritten as e zt 1 1 1 µ (z)= ∞ − + dt t et 1 − t 2 Z0 − 1 1 1 1 where 0 et 1 t + 2 2 for real, non-negative t. Since the integrand is positive for real z = x> 0, ≤ − − ≤ we observe that µ (x) > 0 for real, positive x. In addition, for a complex number z = x + iy, the above integral shows that µ (z) is analytic for Re (z) > 0 and that e xt 1 1 1 µ (z) ∞ − + dt = µ (x) | |≤ t et 1 − t 2 Z0 − In other words, the maximum absolute value of Binet’s function µ (z) for Re(z) > 0 is attained at the positive real axis. Moreover, µ (x) strictly decreases with x = Re(z). Another rather straightforward bound is
1 1 1 1 ∞ xt 1 1 1 1 1 µ (z) max + e− dt = max + | |≤ 0 t t et 1 − t 2 x 0 t t et 1 − t 2 ≤ − Z0 ≤ − Since the maximum occurs at t = 0, the generating function of the Bernoulli numbers Bn
1 1 1 ∞ t2n 1 = + B − for t 2π (12) et 1 t − 2 2n (2n)! | |≤ n=1 − X
1 1 1 1 B2 1 illustrates that max0 t t et 1 t + 2 = 2 = 12 . Hence, for x = Re(z) > 0, we find [28, p. 249] ≤ − − 1 µ (z) µ (x) (13) | |≤ ≤ 12x
Difference µ (z + 1) µ (z). Binet’s definition (1) − 1 1 µ (z) = log Γ (z) z log z + z log (2π) − − 2 − 2 4 illustrates that the complex conjugate µ∗ (z) = µ (z∗), by the reflection principle [22, p. 155]. The functional equation of the Gamma function Γ (z +1) = zΓ (z) leads to 1 1 µ (z + 1) = log Γ (z) + log z z + log (z +1)+ z + 1 log (2π) − 2 − 2 The forward difference ∆µ(z)= µ (z + 1) µ (z) equals − 1 z µ (z + 1) µ (z)= z + log + 1 (14) − 2 z + 1 which is valid for any complex z, except at the negative real axis that is a branch cut for µ (z). In particular, around z = 0, the forward difference (14) shows, with µ (1) = 1 1 log (2π), that − 2 1 1 µ (z) log z log (2π) (15) ∼−2 − 2 illustrating that Binet’s function µ (z) possesses a logarithmic singularity at z = 0. Erdelyi et al. [6, p. 24] deduce from Gauss’s multiplication formula that z+1 1 log Γ (t) dt = z log z z + log (2π) − 2 Zz which we rewrite with the definition (1) as z+1 1 z 1 µ (t) dt = z (z + 1) log + z + (16) 2 z + 1 2 Zz Differentiation of (16) with respect to z again leads to the forward difference (14). For z = 0 in (16), 1 1 we find 0 µ (t) dt = 4 . ReflectionR formula of Binet’s function µ (z). We replace z 1 z in Binet’s definition (1) → − 1 1 µ (1 z) = log Γ (1 z) z log (1 z) z + 1 log (2π) − − − 2 − − − − 2 which added to µ (z) in (1), yields 1 z µ (z)+ µ (1 z) = log (Γ (z)Γ(1 z)) z log + 1 log (2π) − − − − 2 1 z − − After invoking the reflection formula of the Gamma function Γ (z) Γ (1 z) = π , we find the − sin πz corresponding “reflection” formula for Binet’s function 1 z µ (1 z) = 1 µ (z) log (2 sin πz) z log (17) − − − − − 2 1 z − which is valid for any complex number z = x + iy, with the exception of the negative real axis and odd integers x = 2n + 1 with n Z at any y. Since µ (z) is analytic for Re (z) > 0, the Binet ∈ reflection formula (17) illustrates that Binet’s function µ (z) has only logarithmic singularities at negative integers z = k (with k N) including z = 0, as shown in (15). − ∈ Duplication and multiplication formula for µ (z). Combining the duplication formula Γ (2z) = 22z−1 1 √π Γ (z)Γ z + 2 in [1, 6.1.18] for the Gamma function Γ (z) and the definition (1) of Binet’s function µ (z ) leads to 1 1 1 µ (2z)= µ z + + µ (z)+ z log 1+ (18) 2 2z − 2 5 n 1 1 (1 n) 1 − 2 nz 2 k which is generalized by Gauss’s multiplication formula Γ (nz) = (2π) − n − Γ z + n in [1, k=0 6.1.20] as Y
n 1 n 1 − k − k 1 k 1 µ (nz)= µ (z)+ µ z + + z + log z + (n 1) + z log z (19) n n − 2 n − − 2 Xk=1 Xk=1 After letting z = 1 in (19), we find for n 2 n ≥ n 2 n 2 1 − k + 1 − k + 1 1 k + 1 1 1 µ = µ log + (n 1) log n n − n − n − 2 n − 2 − n Xk=1 Xk=1 For example, for n = 2, we obtain µ 1 = 1 (1 log 2) 0.1534 but, for n = 3, µ 1 = 1 µ 2 2 2 − ≃ 3 − 3 − 1 log (2) 1 log (3), which is also immediate from Binet’s reflection formula (17). Gauss’s multiplication 6 − 2 formula already contains the information embedded in the reflection formula Γ (z) Γ (1 z)= π . − sin πz The duplication formula (18) yields a closed expression for z = n + 1 with integer n 1, 2 ≥ 2n 1 1 − 1 µ n + = log k n (log (4n + 2) 1) + (log (2) + 1) 2 − − 2 kX=n 2.1 Complex integral for Binet’s function µ (z)
In Appendix A, we deduce the complex integral (in two ways)
1 c+i π ζ (s) µ (z)= ∞ zsds with 1 1 1 sin πx ∞ log(k) sin (2πxk) log Γ (x) = (γ + log (2π)) x log + 2 − − 2 π πk Xk=2 We evaluate the integral (20) along the line s = c + it, where 1