By Taylor expanding
$$\frac{x}{e^{x}-1} = \sum_{n=0}^\infty \frac{B_n}{n!}x^n$$
in the Zeta function
$$\zeta(s) = \frac{1}{\Gamma(s)} \int_0^\infty x^{s-2} (\frac{x}{e^{x}-1})dx$$
we find
\begin{align} \zeta(s) &= \frac{1}{\Gamma(s)} \int_0^\infty x^{s-2} (\frac{x}{e^{x}-1})dx = \frac{1}{\Gamma(s)} \int_0^\infty x^{s-2} \sum_{n=0}^\infty \frac{B_n}{n!}x^n dx \\ &= \frac{1}{\Gamma(s)} \sum_{n=0}^\infty \frac{B_n}{n!} \int_0^\infty x^{s + n - 2} dx \end{align}
what exactly do we do next to arrive at the formulas for $\zeta(-n)$ and $\zeta(2n)$?