Comenzamos con la representación integral de la función Gamma dada por
$$\Gamma(z)=\int_0^\infty x^{z-1}e^{-x}\,dx\tag1$$
para $z>0$.
En la siguiente sección, mostramos que $\Gamma(z)$ expresada por $(1)$ puede representarse por el límite
$$\Gamma(z)= \lim_{n\to\infty}\frac{n^z\,n!}{z(z+1)(z+2)\cdots (z+n)}$$
Definición de límite de Gamma
Sea $G_n(z)$ la secuencia de funciones dada por
$$G_n(z)=\int_0^n x^{z-1}\left(1-\frac{x}{n}\right)^n\,dx$$
He demostrado en ESTA RESPUESTA, usando solo la Desigualdad de Bernoulli, que la secuencia $\left(1-\frac{x}{n}\right)^n$ aumenta monótonamente para $x\le n$. Por lo tanto, $\left|x^{z-1} \left(1-\frac{x}{n}\right)^n\right|\le x^{z-1}e^{-x}$ para $x\le n$. El Teorema de Convergencia Dominada garantiza que podemos escribir
$$\begin{align} \lim_{n\to \infty} G_n(z)=&\lim_{n\to \infty}\int_0^n x^{z-1}\left(1-\frac{x}{n}\right)^n\,dx\\\\ &=\lim_{n\to \infty}\int_0^\infty \xi_{[0,n]}\,s^{x-1}\left(1-\frac{s}{n}\right)^n\,ds\\\\ &=\int_0^\infty \lim_{n\to \infty} \left(\xi_{[0,n]}\,\left(1-\frac{x}{n}\right)^n\right)\,x^{z-1}\,\,dx\\\\ &=\int_0^\infty x^{z-1}e^{-x}\,dx\\\\ &=\Gamma(z) \end{align}$$
PRUEBA ALTERNATIVA: Definición de límite de Gamma
Si uno no está familiarizado con el Teorema de Convergencia Dominada, entonces simplemente podemos mostrar que
$$\lim_{n\to \infty}\int_0^n x^{z-1}e^{-x}\left(1-e^x\left(1-\frac{x}{n}\right)^n\right)=0$$
Para hacer esto, recurrimos nuevamente al análisis en ESTA RESPUESTA. Procediendo, tenemos
$$\begin{align} 1-e^x\left(1-\frac{x}{n}\right)^n &\le 1-\left(1+\frac{x}{n}\right)^n\left(1-\frac{x}{n}\right)^n\\\\ &=1-\left(1-\frac{x^2}{n^2}\right)^n\\\\ &\le 1-\left(1-\frac{x^2}{n}\right)\\\\ &=\frac{x^2}{n} \end{align}$$
donde se usó la Desigualdad de Bernoulli para llegar a la última desigualdad. De manera similar, vemos que
$$\begin{align} 1-e^x\left(1-\frac{x}{n}\right)^n &\ge 1-e^xe^{-x}\\\\ &=0 \end{align}$$
Por lo tanto, aplicando el teorema del sandwich se obtiene el deseado límite
$$\lim_{n\to \infty}\int_0^n x^{z-1}e^{-x}\left(1-e^x\left(1-\frac{x}{n}\right)^n\right)=0$$
lo que implica que $\lim_{n\to \infty}G_n(z)=\Gamma(z)$.
Integrando por partes repetidamente la representación integral de $G_n(z)$ revela
$$G_n(z)=\frac{n^z\,n!}{z(z+1)(z+2)\cdots (z+n)}$$
de modo que
$$\bbox[5px,border:2px solid #C0A000]{\Gamma(z)=\lim_{n\to \infty}\frac{n^z\,n!}{z(z+1)(z+2)\cdots (z+n)}}\tag2$$
Ahora, usamos $(2)$ para encontrar una representación de límite de la derivada de $\Gamma(z)$. Para facilitar el análisis, usamos $(2)$ para encontrar el logaritmo de $\Gamma(z)$. Procediendo tenemos
$$\log\left(\Gamma(z)\right)=\lim_{n\to \infty}\left(z\log(n)+\log(n!)-\sum_{k=0}^n \log(z+k)\right)\tag3$$
Diferenciando $(3)$ revela
$$\begin{align} \frac{\Gamma'(z)}{\Gamma(z)}&=\lim_{n\to\infty}\left(\log(n)-\sum_{k=0}^n \frac1{z+k}\right)\\\\ &=\lim_{n\to\infty}\left(\log(n)-\sum_{k=0}^n\frac1{k+1}-\sum_{k=0}^n \left(\frac1{z+k}-\frac1{k+1}\right)\right)\\\\ &=-\gamma-\sum_{k=0}^\infty \left(\frac1{z+k}-\frac1{k+1}\right)\tag4 \end{align}$$
Colocando $z=1/2$ en $(4)$ y usando $\Gamma(1/2)=\sqrt \pi$ obtenemos
$$\begin{align} \Gamma'(1/2)&=\sqrt{\pi}\left(-\gamma-\sum_{k=0}^\infty \left(\frac{1}{k+1/2}-\frac1{k+1}\right)\right)\\\\ &=\sqrt{\pi}\left(-\gamma-2\sum_{k=0}^\infty \left(\frac{1}{2k+1}-\frac1{2k+2}\right)\right)\\\\ &=\sqrt{\pi}\left(-\gamma-2\sum_{k=1}^\infty \left(\frac{1}{2k-1}-\frac1{2k}\right)\right)\\\\ &=\sqrt{\pi}\left(-\gamma-2\sum_{k=1}^\infty \frac{(-1)^{k-1}}{k}\right)\\\\ &=-\sqrt\pi\left(\gamma+\log(4)\right) \end{align}$$
¡como se quería mostrar!