Gráfico de la función G(z) de Barnes G también conocida como doble gamma en el plano complejo desde -2-2i hasta 2+2i con colores creados con la función ComplexPlot3D de Mathematica 13.1 La función G de Barnes a lo largo de una parte del eje real En matemáticas , la función G de Barnes G ( z ) es una función que es una extensión de los superfactoriales a los números complejos . Está relacionada con la función gamma , la función K y la constante de Glaisher-Kinkelin , y recibió su nombre en honor al matemático Ernest William Barnes . [1] Puede escribirse en términos de la función gamma doble .
Formalmente, la función G de Barnes se define en la siguiente forma de producto de Weierstrass :
GRAMO ( 1 + el ) = ( 2 π ) el / 2 exp ( − el + el 2 ( 1 + gamma ) 2 ) ∏ a = 1 ∞ { ( 1 + el a ) a exp ( el 2 2 a − el ) } {\displaystyle G(1+z)=(2\pi )^{z/2}\exp \left(-{\frac {z+z^{2}(1+\gamma )}{2}}\right)\,\prod _{k=1}^{\infty }\left\{\left(1+{\frac {z}{k}}\right)^{k}\exp \left({\frac {z^{2}}{2k}}-z\right)\right\}} donde es la constante de Euler-Mascheroni , exp ( x ) = e x es la función exponencial y Π denota multiplicación ( notación pi mayúscula ). gamma {\estilo de visualización \,\gamma }
La representación integral, que puede deducirse de la relación con la función gamma doble , es
registro GRAMO ( 1 + el ) = el 2 registro ( 2 π ) + ∫ 0 ∞ d a a [ 1 − mi − el a 4 pecado 2 a 2 + el 2 2 mi − a − el a ] {\displaystyle \log G(1+z)={\frac {z}{2}}\log(2\pi )+\int _{0}^{\infty }{\frac {dt}{t}}\left[{\frac {1-e^{-zt}}{4\sinh ^{2}{\frac {t}{2}}}}+{\frac {z^{2}}{2}}e^{-t}-{\frac {z}{t}}\right]} Como función completa , G es de orden dos y de tipo infinito. Esto se puede deducir de la expansión asintótica que se muestra a continuación.
Ecuación funcional y argumentos enteros La función G de Barnes satisface la ecuación funcional
GRAMO ( el + 1 ) = Γ ( el ) GRAMO ( el ) {\displaystyle G(z+1)=\Gamma (z)\,G(z)} con normalización G (1) = 1. Nótese la similitud entre la ecuación funcional de la función G de Barnes y la de la función gamma de Euler :
Γ ( el + 1 ) = el Γ ( el ) . {\displaystyle \Gamma (z+1)=z\,\Gamma (z).} La ecuación funcional implica que G toma los siguientes valores en argumentos enteros :
GRAMO ( norte ) = { 0 si norte = 0 , − 1 , − 2 , … ∏ i = 0 norte − 2 i ! si norte = 1 , 2 , … {\displaystyle G(n)={\begin{cases}0&{\text{si }}n=0,-1,-2,\puntos \\\prod _{i=0}^{n-2}i!&{\text{si }}n=1,2,\puntos \end{cases}}} (en particular, ) y por lo tanto GRAMO ( 0 ) = 0 , GRAMO ( 1 ) = 1 {\displaystyle \,G(0)=0,G(1)=1}
GRAMO ( norte ) = ( Γ ( norte ) ) norte − 1 K ( norte ) {\displaystyle G(n)={\frac {(\Gamma (n))^{n-1}}{K(n)}}} donde denota la función gamma y K denota la función K. La ecuación funcional define de forma única la función G de Barnes si se cumple la condición de convexidad, Γ ( incógnita ) {\displaystyle \,\Gamma (x)}
( ∀ incógnita ≥ 1 ) d 3 d incógnita 3 registro ( GRAMO ( incógnita ) ) ≥ 0 {\displaystyle (\para todo x\geq 1)\,{\frac {\mathrm {d} ^{3}}{\mathrm {d} x^{3}}}\log(G(x))\geq 0} se añade. [2] Además, la función G de Barnes satisface la fórmula de duplicación, [3]
GRAMO ( incógnita ) GRAMO ( incógnita + 1 2 ) 2 GRAMO ( incógnita + 1 ) = mi 1 4 A − 3 2 − 2 incógnita 2 + 3 incógnita − 11 12 π incógnita − 1 2 GRAMO ( 2 incógnita ) {\displaystyle G(x)G(x+{\frac {1}{2}}\right)^{2}G(x+1)=e^{\frac {1}{4}}A^{-3}2^{-2x^{2}+3x-{\frac {11}{12}}}\pi ^{x-{\frac {1}{2}}}G(2x\right)} ,¿Dónde está la constante de Glaisher-Kinkelin ? A {\estilo de visualización A}
Caracterización De manera similar al teorema de Bohr-Mollerup para la función gamma , para una constante , tenemos para [4] do > 0 {\displaystyle c>0} F ( incógnita ) = do GRAMO ( incógnita ) {\displaystyle f(x)=cG(x)}
F ( incógnita + 1 ) = Γ ( incógnita ) F ( incógnita ) {\displaystyle f(x+1)=\Gamma(x)f(x)}
y para incógnita > 0 {\displaystyle x>0}
F ( incógnita + norte ) ∼ Γ ( incógnita ) norte norte ( incógnita 2 ) F ( norte ) {\displaystyle f(x+n)\sim \Gamma (x)^{n}n^{x \choose 2}f(n)}
como . norte → ∞ {\displaystyle n\to \infty}
La ecuación diferencial para la función G, junto con la ecuación funcional para la función gamma , se puede utilizar para obtener la siguiente fórmula de reflexión para la función G de Barnes (probada originalmente por Hermann Kinkelin ):
registro GRAMO ( 1 − el ) = registro GRAMO ( 1 + el ) − el registro 2 π + ∫ 0 el π incógnita cuna π incógnita d incógnita . {\displaystyle \log G(1-z)=\log G(1+z)-z\log 2\pi +\int _{0}^{z}\pi x\cot \pi x\,dx. } La integral logaritmo-tangente del lado derecho se puede evaluar en términos de la función de Clausen (de orden 2), como se muestra a continuación:
2 π registro ( GRAMO ( 1 − el ) GRAMO ( 1 + el ) ) = 2 π el registro ( pecado π el π ) + Cl 2 ( 2 π el ) {\displaystyle 2\pi \log \left({\frac {G(1-z)}{G(1+z)}}\right)=2\pi z\log \left({\frac {\sin \pi z}{\pi }}\right)+\operatorname {Cl} _{2}(2\pi z)} La prueba de este resultado depende de la siguiente evaluación de la integral cotangente: introduciendo la notación para la integral logarítmica-cotangente y utilizando el hecho de que , una integración por partes da Lc ( z ) {\displaystyle \operatorname {Lc} (z)} ( d / d x ) log ( sin π x ) = π cot π x {\displaystyle \,(d/dx)\log(\sin \pi x)=\pi \cot \pi x}
Lc ( z ) = ∫ 0 z π x cot π x d x = z log ( sin π z ) − ∫ 0 z log ( sin π x ) d x = z log ( sin π z ) − ∫ 0 z [ log ( 2 sin π x ) − log 2 ] d x = z log ( 2 sin π z ) − ∫ 0 z log ( 2 sin π x ) d x . {\displaystyle {\begin{aligned}\operatorname {Lc} (z)&=\int _{0}^{z}\pi x\cot \pi x\,dx\\&=z\log(\sin \pi z)-\int _{0}^{z}\log(\sin \pi x)\,dx\\&=z\log(\sin \pi z)-\int _{0}^{z}{\Bigg [}\log(2\sin \pi x)-\log 2{\Bigg ]}\,dx\\&=z\log(2\sin \pi z)-\int _{0}^{z}\log(2\sin \pi x)\,dx.\end{aligned}}} Realizando la sustitución integral se obtiene y = 2 π x ⇒ d x = d y / ( 2 π ) {\displaystyle \,y=2\pi x\Rightarrow dx=dy/(2\pi )}
z log ( 2 sin π z ) − 1 2 π ∫ 0 2 π z log ( 2 sin y 2 ) d y . {\displaystyle z\log(2\sin \pi z)-{\frac {1}{2\pi }}\int _{0}^{2\pi z}\log \left(2\sin {\frac {y}{2}}\right)\,dy.} La función de Clausen – de segundo orden – tiene la representación integral
Cl 2 ( θ ) = − ∫ 0 θ log | 2 sin x 2 | d x . {\displaystyle \operatorname {Cl} _{2}(\theta )=-\int _{0}^{\theta }\log {\Bigg |}2\sin {\frac {x}{2}}{\Bigg |}\,dx.} Sin embargo, dentro del intervalo , el signo de valor absoluto dentro del integrando se puede omitir, ya que dentro del rango la función 'semi-seno' en la integral es estrictamente positiva y estrictamente distinta de cero. Comparando esta definición con el resultado anterior para la integral logaritmo-tangente, se cumple claramente la siguiente relación: 0 < θ < 2 π {\displaystyle \,0<\theta <2\pi }
Lc ( z ) = z log ( 2 sin π z ) + 1 2 π Cl 2 ( 2 π z ) . {\displaystyle \operatorname {Lc} (z)=z\log(2\sin \pi z)+{\frac {1}{2\pi }}\operatorname {Cl} _{2}(2\pi z).} Así, después de un ligero reordenamiento de términos, la prueba está completa:
2 π log ( G ( 1 − z ) G ( 1 + z ) ) = 2 π z log ( sin π z π ) + Cl 2 ( 2 π z ) . ◻ {\displaystyle 2\pi \log \left({\frac {G(1-z)}{G(1+z)}}\right)=2\pi z\log \left({\frac {\sin \pi z}{\pi }}\right)+\operatorname {Cl} _{2}(2\pi z)\,.\,\Box } Usando la relación y dividiendo la fórmula de reflexión por un factor de se obtiene la forma equivalente: G ( 1 + z ) = Γ ( z ) G ( z ) {\displaystyle \,G(1+z)=\Gamma (z)\,G(z)} 2 π {\displaystyle \,2\pi }
log ( G ( 1 − z ) G ( z ) ) = z log ( sin π z π ) + log Γ ( z ) + 1 2 π Cl 2 ( 2 π z ) {\displaystyle \log \left({\frac {G(1-z)}{G(z)}}\right)=z\log \left({\frac {\sin \pi z}{\pi }}\right)+\log \Gamma (z)+{\frac {1}{2\pi }}\operatorname {Cl} _{2}(2\pi z)} Adamchik (2003) ha dado una forma equivalente de la fórmula de reflexión , pero con una prueba diferente. [5]
Reemplazando z por 1 / 2 − z en la fórmula de reflexión anterior da, después de cierta simplificación, la fórmula equivalente que se muestra a continuación (que involucra polinomios de Bernoulli ):
log ( G ( 1 2 + z ) G ( 1 2 − z ) ) = log Γ ( 1 2 − z ) + B 1 ( z ) log 2 π + 1 2 log 2 + π ∫ 0 z B 1 ( x ) tan π x d x {\displaystyle \log \left({\frac {G\left({\frac {1}{2}}+z\right)}{G\left({\frac {1}{2}}-z\right)}}\right)=\log \Gamma \left({\frac {1}{2}}-z\right)+B_{1}(z)\log 2\pi +{\frac {1}{2}}\log 2+\pi \int _{0}^{z}B_{1}(x)\tan \pi x\,dx}
Expansión de la serie de Taylor Por el teorema de Taylor , y considerando las derivadas logarítmicas de la función de Barnes, se puede obtener el siguiente desarrollo en serie:
log G ( 1 + z ) = z 2 log 2 π − ( z + ( 1 + γ ) z 2 2 ) + ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 . {\displaystyle \log G(1+z)={\frac {z}{2}}\log 2\pi -\left({\frac {z+(1+\gamma )z^{2}}{2}}\right)+\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}.} Es válido para . Aquí está la función zeta de Riemann : 0 < z < 1 {\displaystyle \,0<z<1} ζ ( x ) {\displaystyle \,\zeta (x)}
ζ ( s ) = ∑ n = 1 ∞ 1 n s . {\displaystyle \zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}.} Exponenciando ambos lados de la expansión de Taylor obtenemos:
G ( 1 + z ) = exp [ z 2 log 2 π − ( z + ( 1 + γ ) z 2 2 ) + ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 ] = ( 2 π ) z / 2 exp [ − z + ( 1 + γ ) z 2 2 ] exp [ ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 ] . {\displaystyle {\begin{aligned}G(1+z)&=\exp \left[{\frac {z}{2}}\log 2\pi -\left({\frac {z+(1+\gamma )z^{2}}{2}}\right)+\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}\right]\\&=(2\pi )^{z/2}\exp \left[-{\frac {z+(1+\gamma )z^{2}}{2}}\right]\exp \left[\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}\right].\end{aligned}}} Comparando esto con la forma del producto de Weierstrass de la función de Barnes se obtiene la siguiente relación:
exp [ ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 ] = ∏ k = 1 ∞ { ( 1 + z k ) k exp ( z 2 2 k − z ) } {\displaystyle \exp \left[\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}\right]=\prod _{k=1}^{\infty }\left\{\left(1+{\frac {z}{k}}\right)^{k}\exp \left({\frac {z^{2}}{2k}}-z\right)\right\}}
Al igual que la función gamma, la función G también tiene una fórmula de multiplicación: [6]
G ( n z ) = K ( n ) n n 2 z 2 / 2 − n z ( 2 π ) − n 2 − n 2 z ∏ i = 0 n − 1 ∏ j = 0 n − 1 G ( z + i + j n ) {\displaystyle G(nz)=K(n)n^{n^{2}z^{2}/2-nz}(2\pi )^{-{\frac {n^{2}-n}{2}}z}\prod _{i=0}^{n-1}\prod _{j=0}^{n-1}G\left(z+{\frac {i+j}{n}}\right)} donde es una constante dada por: K ( n ) {\displaystyle K(n)}
K ( n ) = e − ( n 2 − 1 ) ζ ′ ( − 1 ) ⋅ n 5 12 ⋅ ( 2 π ) ( n − 1 ) / 2 = ( A e − 1 12 ) n 2 − 1 ⋅ n 5 12 ⋅ ( 2 π ) ( n − 1 ) / 2 . {\displaystyle K(n)=e^{-(n^{2}-1)\zeta ^{\prime }(-1)}\cdot n^{\frac {5}{12}}\cdot (2\pi )^{(n-1)/2}\,=\,(Ae^{-{\frac {1}{12}}})^{n^{2}-1}\cdot n^{\frac {5}{12}}\cdot (2\pi )^{(n-1)/2}.} Aquí está la derivada de la función zeta de Riemann y es la constante de Glaisher-Kinkelin . ζ ′ {\displaystyle \zeta ^{\prime }} A {\displaystyle A}
Valor absoluto Es cierto que , por lo tanto . A partir de esta relación y por la forma del producto de Weierstrass presentada anteriormente se puede demostrar que G ( z ¯ ) = G ( z ) ¯ {\displaystyle G({\overline {z}})={\overline {G(z)}}} | G ( z ) | 2 = G ( z ) G ( z ¯ ) {\displaystyle |G(z)|^{2}=G(z)G({\overline {z}})}
| G ( x + i y ) | = | G ( x ) | exp ( y 2 1 + γ 2 ) 1 + y 2 x 2 ∏ k = 1 ∞ ( 1 + y 2 ( x + k ) 2 ) k + 1 exp ( − y 2 k ) . {\displaystyle |G(x+iy)|=|G(x)|\exp \left(y^{2}{\frac {1+\gamma }{2}}\right){\sqrt {1+{\frac {y^{2}}{x^{2}}}}}{\sqrt {\prod _{k=1}^{\infty }\left(1+{\frac {y^{2}}{(x+k)^{2}}}\right)^{k+1}\exp \left(-{\frac {y^{2}}{k}}\right)}}.} Esta relación es válida para valores arbitrarios de , y . Si , entonces la siguiente fórmula es válida: x ∈ R ∖ { 0 , − 1 , − 2 , … } {\displaystyle x\in \mathbb {R} \setminus \{0,-1,-2,\dots \}} y ∈ R {\displaystyle y\in \mathbb {R} } x = 0 {\displaystyle x=0}
| G ( i y ) | = y exp ( y 2 1 + γ 2 ) ∏ k = 1 ∞ ( 1 + y 2 k 2 ) k + 1 exp ( − y 2 k ) {\displaystyle |G(iy)|=y\exp \left(y^{2}{\frac {1+\gamma }{2}}\right){\sqrt {\prod _{k=1}^{\infty }\left(1+{\frac {y^{2}}{k^{2}}}\right)^{k+1}\exp \left(-{\frac {y^{2}}{k}}\right)}}} para realidad arbitraria .
Expansión asintótica El logaritmo de G ( z + 1) tiene la siguiente expansión asintótica, según lo establecido por Barnes:
log G ( z + 1 ) = z 2 2 log z − 3 z 2 4 + z 2 log 2 π − 1 12 log z + ( 1 12 − log A ) + ∑ k = 1 N B 2 k + 2 4 k ( k + 1 ) z 2 k + O ( 1 z 2 N + 2 ) . {\displaystyle {\begin{aligned}\log G(z+1)={}&{\frac {z^{2}}{2}}\log z-{\frac {3z^{2}}{4}}+{\frac {z}{2}}\log 2\pi -{\frac {1}{12}}\log z\\&{}+\left({\frac {1}{12}}-\log A\right)+\sum _{k=1}^{N}{\frac {B_{2k+2}}{4k\left(k+1\right)z^{2k}}}~+~O\left({\frac {1}{z^{2N+2}}}\right).\end{aligned}}} Aquí están los números de Bernoulli y es la constante de Glaisher-Kinkelin . (Nótese que, de manera algo confusa en la época de Barnes [7], el número de Bernoulli se habría escrito como , pero esta convención ya no es actual). Esta expansión es válida para en cualquier sector que no contenga el eje real negativo con grandes. B k {\displaystyle B_{k}} A {\displaystyle A} B 2 k {\displaystyle B_{2k}} ( − 1 ) k + 1 B k {\displaystyle (-1)^{k+1}B_{k}} z {\displaystyle z} | z | {\displaystyle |z|}
Relación con la integral log-gamma El log-gamma paramétrico se puede evaluar en términos de la función G de Barnes: [5]
∫ 0 z log Γ ( x ) d x = z ( 1 − z ) 2 + z 2 log 2 π + z log Γ ( z ) − log G ( 1 + z ) {\displaystyle \int _{0}^{z}\log \Gamma (x)\,dx={\frac {z(1-z)}{2}}+{\frac {z}{2}}\log 2\pi +z\log \Gamma (z)-\log G(1+z)} Una prueba de la fórmula
La prueba es algo indirecta e implica considerar primero la diferencia logarítmica de la función gamma y la función G de Barnes:
z log Γ ( z ) − log G ( 1 + z ) {\displaystyle z\log \Gamma (z)-\log G(1+z)} dónde
1 Γ ( z ) = z e γ z ∏ k = 1 ∞ { ( 1 + z k ) e − z / k } {\displaystyle {\frac {1}{\Gamma (z)}}=ze^{\gamma z}\prod _{k=1}^{\infty }\left\{\left(1+{\frac {z}{k}}\right)e^{-z/k}\right\}} y es la constante de Euler-Mascheroni . γ {\displaystyle \,\gamma }
Tomando el logaritmo de las formas del producto de Weierstrass de la función G de Barnes y la función gamma se obtiene:
z log Γ ( z ) − log G ( 1 + z ) = − z log ( 1 Γ ( z ) ) − log G ( 1 + z ) = − z [ log z + γ z + ∑ k = 1 ∞ { log ( 1 + z k ) − z k } ] − [ z 2 log 2 π − z 2 − z 2 2 − z 2 γ 2 + ∑ k = 1 ∞ { k log ( 1 + z k ) + z 2 2 k − z } ] {\displaystyle {\begin{aligned}&z\log \Gamma (z)-\log G(1+z)=-z\log \left({\frac {1}{\Gamma (z)}}\right)-\log G(1+z)\\[5pt]={}&{-z}\left[\log z+\gamma z+\sum _{k=1}^{\infty }{\Bigg \{}\log \left(1+{\frac {z}{k}}\right)-{\frac {z}{k}}{\Bigg \}}\right]\\[5pt]&{}-\left[{\frac {z}{2}}\log 2\pi -{\frac {z}{2}}-{\frac {z^{2}}{2}}-{\frac {z^{2}\gamma }{2}}+\sum _{k=1}^{\infty }{\Bigg \{}k\log \left(1+{\frac {z}{k}}\right)+{\frac {z^{2}}{2k}}-z{\Bigg \}}\right]\end{aligned}}} Una pequeña simplificación y reordenamiento de términos da como resultado la expansión de la serie:
∑ k = 1 ∞ { ( k + z ) log ( 1 + z k ) − z 2 2 k − z } = − z log z − z 2 log 2 π + z 2 + z 2 2 − z 2 γ 2 − z log Γ ( z ) + log G ( 1 + z ) {\displaystyle {\begin{aligned}&\sum _{k=1}^{\infty }{\Bigg \{}(k+z)\log \left(1+{\frac {z}{k}}\right)-{\frac {z^{2}}{2k}}-z{\Bigg \}}\\[5pt]={}&{-z}\log z-{\frac {z}{2}}\log 2\pi +{\frac {z}{2}}+{\frac {z^{2}}{2}}-{\frac {z^{2}\gamma }{2}}-z\log \Gamma (z)+\log G(1+z)\end{aligned}}} Finalmente, tome el logaritmo del producto de Weierstrass de la función gamma e integre sobre el intervalo para obtener: [ 0 , z ] {\displaystyle \,[0,\,z]}
∫ 0 z log Γ ( x ) d x = − ∫ 0 z log ( 1 Γ ( x ) ) d x = − ( z log z − z ) − z 2 γ 2 − ∑ k = 1 ∞ { ( k + z ) log ( 1 + z k ) − z 2 2 k − z } {\displaystyle {\begin{aligned}&\int _{0}^{z}\log \Gamma (x)\,dx=-\int _{0}^{z}\log \left({\frac {1}{\Gamma (x)}}\right)\,dx\\[5pt]={}&{-(z\log z-z)}-{\frac {z^{2}\gamma }{2}}-\sum _{k=1}^{\infty }{\Bigg \{}(k+z)\log \left(1+{\frac {z}{k}}\right)-{\frac {z^{2}}{2k}}-z{\Bigg \}}\end{aligned}}} Igualando las dos evaluaciones se completa la prueba:
∫ 0 z log Γ ( x ) d x = z ( 1 − z ) 2 + z 2 log 2 π + z log Γ ( z ) − log G ( 1 + z ) {\displaystyle \int _{0}^{z}\log \Gamma (x)\,dx={\frac {z(1-z)}{2}}+{\frac {z}{2}}\log 2\pi +z\log \Gamma (z)-\log G(1+z)} Y desde entonces, G ( 1 + z ) = Γ ( z ) G ( z ) {\displaystyle \,G(1+z)=\Gamma (z)\,G(z)}
∫ 0 z log Γ ( x ) d x = z ( 1 − z ) 2 + z 2 log 2 π − ( 1 − z ) log Γ ( z ) − log G ( z ) . {\displaystyle \int _{0}^{z}\log \Gamma (x)\,dx={\frac {z(1-z)}{2}}+{\frac {z}{2}}\log 2\pi -(1-z)\log \Gamma (z)-\log G(z)\,.}
Referencias ^ EW Barnes, "La teoría de la función G", Quarterly Journ. Pure and Appl. Math. 31 (1900), 264–314. ^ MF Vignéras, L'équation fonctionelle de la fonction zêta de Selberg du groupe mudulaire SL ( 2 , Z ) {\displaystyle (2,\mathbb {Z} )} , Astérisque 61 , 235-249 (1979). ^ Park, Junesang (1996). "Una fórmula de duplicación para la función gamma doble $Gamma_2$". Boletín de la Sociedad Matemática Coreana . 33 (2): 289–294. ^ Marichal, Jean Luc. Una generalización del teorema de Bohr-Mollerup para funciones convexas de orden superior (PDF) . Springer. pág. 218. ^ ab Adamchik, Viktor S. (2003). "Contribuciones a la teoría de la función de Barnes". arXiv : math/0308086 . ^ I. Vardi, Determinantes de laplacianos y funciones gamma múltiples , SIAM J. Math. Anal. 19 , 493–507 (1988). ^ ET Whittaker y GN Watson , " Un curso de análisis moderno ", CUP.Askey, RA; Roy, R. (2010), "Función G de Barnes", en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), Manual del NIST de funciones matemáticas , Cambridge University Press, ISBN 978-0-521-19225-5 , Sr. 2723248 .