I know that $\zeta(2n)$ can be found using the contour integral $$\int_{\Gamma_N} \frac{cot(\pi z)}{z^{2n}}dz$$ and computing its residues as $N\to\infty$.
My question is could other values of the zeta function be derived through an integral such as $$\int_{\gamma_n} \frac{\Gamma(z)}{z^s}dz$$ since $\Gamma(z)$ has poles at every negative integer. I do have alot of doubts about it, since I am not sure that the contour will approach zero as $n\to\infty$, and I find it hard to believe that this was not already considered by someone else.
However, Wolfram Residue Calculator shows interesting values for the residues at zero for some s.
First few values for $s=1\to s=4$ are $$Res_{s=1}=-\gamma$$ $$Res_{s=2}=\frac{1}{12}(6\gamma^2+\pi^2)$$ $$Res_{s=3}=\frac{1}{12}(-2\gamma^3-\gamma\pi^2+2\psi^2(1))$$ $$Res_{s=4}=\frac{1}{480}(20\gamma^4+20\gamma^2\pi^2+3\pi^4-80\gamma\psi^2(1))$$ If these values do not have anything to do with $\zeta(s)$, can anyone at least provide some meaning to why these values appear and if there is any pattern/formula for deriving them?