Tenemos para$s$$\Re{e}(s) > 1$ :
$\Gamma(s)\zeta(s)=\int_0^\infty\frac{t^{s-1}}{{\rm e}^t-1}\mathrm dt=\int_0^1\frac{t^{s-1}}{\mathrm e^t-1}\mathrm dt+\int_1^\infty\frac{t^{s-1}}{\mathrm e^t-1} dt$.
La segunda integral es holomorphe en $s$. Tomamos la serie de Taylor en la primera integral.
Tenemos para todos los $t$ $|t| < 2\pi$
$\frac t{{\rm e}^t-1}=\sum_{n=0}^\infty\frac{B_n t^n}{n!}$,
donde el $B_n$ son los números de Bernoulli. La integración término a término
$\zeta(s)=\frac1{(s-1)\Gamma(s)}+\frac1{\Gamma(s)}\sum_{n=1}^\infty\frac{B_n}{n!(n+s-1)}+\frac1{\Gamma(s)}\int_1^\infty\frac{t^{s-1}}{\mathrm e^t-1} \mathrm dt$.
La serie converge y es holomorphic para todos los s, pero $s=-n$, $(n \in \mathbb{N})$
debido a que el radio de convergencia de la serie no es modificado por la división por $n + s – 1$.
Al $s \to k$ $Γ(s)$ tiene una simple polo en $s=–k$ , $ζ(s)$ es la suma de un término que tienden a $0$ $\frac1{\Gamma(s)}~\frac{B_{k+1}}{(k+1)!~(s+k)}~\underset{\overset{s\to-k}{}}{\sim}~\frac{s+k}{\text{Res}(\Gamma,-k)}~\frac{B_{k+1}}{(k+1)!~(s+k)}=(-1)^k~k!~\frac{B_{k+1}}{(k+1)!}.$
Tenemos la fórmula de Euler :
$\zeta(-k)=(-1)^k\frac{B_{k+1}}{k+1}.$