Cómo obtener la expansión de Laurent de la función gamma alrededor de $z=0$ ?

Question

Quiero probar, la expansión de Laurent de la función gamma. $$\begin{align} \Gamma(z) = \frac1z-\gamma+\frac12\left(\gamma^2+\frac {\pi^2}6\right)z-\frac16\left(\gamma^3+\frac {\gamma\pi^2}2+2 \zeta(3)\right)z^2+O(z^3). \end{align}$$

En primer lugar, mi suposición de obtener la expansión anterior, es a partir de las definiciones de la función gamma $$\begin{align} \Gamma(z) &= \int_0^{\infty} dt e^{-t} t^{z-1} \\ & = \int_1^\infty dt e^{-t}t^{z-1} + \int_0^1dt e^{-t} t^{z-1} \\ & = \int_1^\infty dte^{-t}t^{z-1} + \int_0^1 dt t^{z-1} \sum_{n=0}^{\infty}\frac{(-1)^n}{n!}t^n \\ & = \int_1^\infty dt e^{-t}t^{z-1} + \sum_{n=0}^{\infty} \frac{(-1)^n}{n!}\frac{1}{z+n} \end{align}$$ Esto sólo da la función gamma en función de $\frac{1}{z}$ ...

¿O debería empezar con $$\begin{align} \Gamma(z) = \lim_{n \rightarrow \infty} \frac{n! n^z}{z(z+1) \cdots(z+n)} \end{align}$$

Answer 1

Dejemos que $\Gamma(z)$ ser representado por la integral

$$\Gamma(z)=\int_0^\infty x^{z-1}e^{-x}\,dx\tag1$$

para $\text{Re}(z)>0$ . Integrar por partes la integral en $(1)$ con $u=e^{-x}$ y $v=\frac{x^z}{z}$ revela

$$\Gamma(z)=\frac1z\int_0^\infty x^ze^{-x}\,dx\tag2$$

A continuación, ampliamos $x^z$ en una serie de potencias de $z$ para obtener

$$\begin{align} \Gamma(z)&=\frac1z\sum_{n=0}^\infty \frac{z^n}{n!}\int_0^\infty e^{-x}\log^n(x)\,dx\\\\ &=\frac1z+\underbrace{\int_0^\infty \log(x)e^{-x}\,dx}_{=-\gamma}+\frac12 z\underbrace{\int_0^\infty \log^2(x)e^{-x}\,dx}_{\gamma^2+\frac{\pi^2}6}+\frac16z^2 \underbrace{\int_0^\infty \log^3(x)e^{-x}\,dx}_{-\gamma^3-\gamma\pi^2/2-2\zeta(3)}+O(z^3) \end{align}$$

como se iba a demostrar.

---

NOTA:

Los coeficientes $\int_0^\infty \log^n(x)e^{-x}\,dx$ para $n=2,3$ se puede encontrar utilizando la relación $\Gamma'(x)=\Gamma(x)\psi(x)$ entre la derivada logarítmica de la función Gamma y la función Digamma, junto con los valores de $\psi'(1)=\zeta(2)$ y $\psi''(1)=-2\zeta(3)$ .

Answer 2

Una idea para que la desarrolles:

La fórmula de Weierstrass nos dice que

$$\frac1{\Gamma(z)}=ze^{\gamma z}\prod_{n=1}^\infty\left(1+\frac zn\right)e^{-z/n}$$

Ahora toma logaritmos en ambos lados para obtener una relación más o menos conocida:

$$-\log\Gamma(z)=\log z+\gamma z+\sum_{n=1}^\infty\left[\log\left(1+\frac zn\right)-\frac zn\right]$$

Ahora diferencia lo anterior para obtener la derivada logarítmica de la función Gamma:

$$\frac{\Gamma'(z)}{\Gamma(z)}=-\frac1z-\gamma-\sum_{n=1}^\infty\frac1n\left[\frac n{z+n}-1\right]=-\frac1z-\gamma+\sum_{n=1}^\infty\frac z{n(z+n)}$$

y etc. Puede intentar integrar el $\;-\dfrac1z\;$ término en la serie, también.