It is well-known applying the monotone convergence theorem that for all $x>1$ we have the following functional equation : $\zeta(x)\Gamma(x)= \displaystyle \int \limits_{[0,+\infty]} \dfrac{t^{x-1}}{\exp(t)-1}\mathrm{d}\lambda(t)$ or $\zeta(x)=\dfrac{1}{\Gamma(x)}\displaystyle \int \limits_{[0,+\infty]} \dfrac{t^{x-1}}{\exp(t)-1}\mathrm{d}\lambda(t) \ $ (where $\lambda$ is the Lebesgue's measure).
Recall that : $\Gamma(x) = \displaystyle \int \limits_{[0,+\infty]}\exp(-t)t^{x-1}\mathrm{d}\lambda(t)$ and $\zeta(x)= \displaystyle \sum\limits_{n=1}^{+\infty} \dfrac{1}{n^x}$ and are well-defined for all $x>1$.
Is it possible to find a link between this functional equation and Bernoulli's numbers ?
I was thinking about the following result (which is not trivial) : $\dfrac{t}{\exp(t)-1} = \displaystyle \sum\limits_{n=0}^{+\infty}B_n \dfrac{t^n}{n!} \ $ where the sequence $(B_n)_{n\ge 0}$ represents the Bernoulli's numbers.
Then we will obtain for all $x > 1$ : $\zeta(x)\Gamma(x)= \displaystyle \int \limits_{[0,1]}t^{x-2}\left( \displaystyle \sum\limits_{n=0}^{+\infty}B_n \dfrac{t^n}{n!} \right) \mathrm{d}\lambda(t) + \displaystyle \int \limits_{[1,+\infty]} \dfrac{t^{x-1}}{\exp(t)-1} \mathrm{d}\lambda(t) \\ = \displaystyle \int \limits_{[0,1]}\left( \displaystyle \sum\limits_{n=0}^{+\infty}B_n \dfrac{t^{n+x-2}}{n!} \right) \mathrm{d}\lambda(t) + \displaystyle \int \limits_{[1,+\infty]} \dfrac{t^{x-1}}{\exp(t)-1} \mathrm{d}\lambda(t) \\ = \displaystyle\sum \limits_{n=0}^{+\infty} \dfrac{B_n}{n!(n+x-1)} + \displaystyle \int \limits_{[1,+\infty]} \dfrac{t^{x-1}}{\exp(t)-1} \mathrm{d}\lambda(t)$
The last equality is obtained using the fact that the series converges uniformly on every compact set included in $[-2\pi, 2\pi]$.
I doubt that we could go further ?
Thnaks in advance !