About the Project
5 Gamma FunctionProperties

§5.17 Barnes’ G-Function (Double Gamma Function)

5.17.1
G\left(z+1\right)=\Gamma\left(z\right)G\left(z\right),
G\left(1\right)=1,

In this equation (and in (5.17.5) below), the \operatorname{Ln}’s have their principal values on the positive real axis and are continued via continuity, as in §4.2(i).

When z\to\infty in |\operatorname{ph}z|\leq\pi-\delta\;(<\pi),

5.17.5 \operatorname{Ln}G\left(z+1\right)\sim\tfrac{1}{4}z^{2}+z\operatorname{Ln}%
\Gamma\left(z+1\right)-\left(\tfrac{1}{2}z(z+1)+\tfrac{1}{12}\right)\ln z-\ln A%
+\sum_{k=1}^{\infty}\frac{B_{2k+2}}{2k(2k+1)(2k+2)z^{2k}}.

For error bounds and an exponentially-improved extension, see Nemes (2014a). Here B_{2k+2} is the Bernoulli number (§24.2(i)), and A is Glaisher’s constant, given by

5.17.6 A=e^{C}=1.28242\;71291\;00622\;63687\;\ldots,

where

and \zeta' is the derivative of the zeta function (Chapter 25).

For Glaisher’s constant see also Greene and Knuth (1982, p. 100) and §2.10(i).