Supongamos que empezar por evalutating la suma de dos en turno, donde el
el parámetro $k\ge 1$.
Para el primer lugar vamos a utilizar la siguiente representación integral:
$${n+k\elegir k}
= \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{n+k}}{z^{k+1}} \; dz.$$
Buscamos $$\sum_{n\ge 1} {n-1+k\choose k} {n+k\choose k} x^n.$$
El uso de la integral nos encontramos
$$\frac{1}{2\pi i}
\int_{|z|=\epsilon}
\sum_{n\ge 1} {n-1+k\elegir k} x^n
\frac{(1+z)^{n+k}}{z^{k+1}} \; dz
\\ = \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^k}{z^{k+1}}
\sum_{n\ge 1} {n-1+k\elegir k} (1+z)^n x^n \; dz
\\= \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{x(1+z)^{k+1}}{z^{k+1}}
\sum_{n\ge 1} {n-1+k\elegir k} (1+z)^{n-1} x^{n-1} \; dz
\\= \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{x(1+z)^{k+1}}{z^{k+1}}
\frac{1}{(1-x(1+z))^{k+1}} \; dz
\\= \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{x(1+z)^{k+1}}{z^{k+1}}
\frac{1}{(1-x-xz))^{k+1}} \; dz
\\= \frac{1}{(1-x)^{k+1}} \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{x(1+z)^{k+1}}{z^{k+1}}
\frac{1}{(1-xz/(1-x)))^{k+1}} \; dz
\\ = \frac{x}{(1-x)^{k+1}}
\sum_{q=0}^k {k+1\elegir k-q}
{q+k\elegir k} \left( \frac{x}{1-x} \right)^q
\\ = \frac{x}{(1-x)^{k+1}}
\sum_{q=0}^k {k+1\elegir q+1}
{q+k\elegir k} \left( \frac{x}{1-x} \right)^q.$$
La aplicación de la representación integral desde el principio una segunda vez
obtenemos de esta suma
$$\frac{x}{(1-x)^{k+1}} \frac{1}{2\pi i}
\int_{|z|=\epsilon}
\sum_{q=0}^k {k+1\elegir q+1}
\frac{(1+z)^{q+k}}{z^{k+1}}
\left( \frac{x}{1-x} \right)^q\; dz
\\ = \frac{x}{(1-x)^{k+1}} \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^k}{z^{k+1}}
\sum_{q=0}^k {k+1\elegir q+1}
(1+z)^q
\left( \frac{x}{1-x} \right)^q\; dz
\\ = \frac{1}{(1-x)^k} \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{k-1}}{z^{k+1}}
\sum_{q=0}^k {k+1\elegir q+1}
(1+z)^{q+1}
\left( \frac{x}{1-x} \right)^{q+1}\; dz
\\ = \frac{1}{(1-x)^k} \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{k-1}}{z^{k+1}}
\left(-1 + \left(1+(1+z)\frac{x}{1-x}\right)^{k+1}\right) \; dz.$$
Tenemos $k+1-(k-1) = 2$, por lo que el primer componente dentro de los paréntesis
se retira, dejando
$$ \frac{1}{(1-x)^k} \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{k-1}}{z^{k+1}}
\left(1+(1+z)\frac{x}{1-x}\right)^{k+1} \; dz
\\ = \frac{1}{(1-x)^{2k+1}} \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{k-1}}{z^{k+1}}
\left(1-x+x(1+z)\right)^{k+1} \; dz
\\ = \frac{1}{(1-x)^{2k+1}} \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{k-1}}{z^{k+1}}
\left(1+xz\right)^{k+1} \; dz.$$
Necesitamos una mayor simplificación en esto y poner $z=1/w$, al pasar
$$\frac{1}{(1-x)^{2k+1}} \frac{1}{2\pi i}
\int_{|w|=\epsilon} \frac{(1+1/n)^{k-1}}{(1/w)^{k+1}}
\left(1+x/w\ \ derecho)^{k+1} \; \frac{1}{w^2} dw
\\ = \frac{1}{(1-x)^{2k+1}} \frac{1}{2\pi i}
\int_{|w|=\epsilon} w^2 (w+1)^{k-1}
\left(\frac{w+x}{w} \right)^{k+1} \; \frac{1}{w^2} dw
\\ = \frac{1}{(1-x)^{2k+1}} \frac{1}{2\pi i}
\int_{|w|=\epsilon} \frac{(w+1)^{k-1}}{w^{k+1}} (w+x)^{k+1}\; ps.$$
El reson en que esto funciona es porque somos esencialmente la evaluación de los residuos en el infinito y los residuos de suma cero.
Esto concluye la evaluación de la primera suma.
Para el segundo vamos a utilizar la siguiente representación integral:
$${k-1\elegir j-1}
= \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{k-1}}{z^j} \; dz.$$
Buscamos $$\sum_{j\ge 1} {k+1\choose j} {k-1\choose j-1} x^j.$$
El uso de la integral nos encontramos
$$\frac{1}{2\pi i}
\int_{|z|=\epsilon}
\sum_{j\ge 1} {k+1\elegir j} x^j
\frac{(1+z)^{k-1}}{z^j} \; dz
\\= \frac{1}{2\pi i}
\int_{|z|=\epsilon} (1+z)^{k-1}
\sum_{j\ge 1} {k+1\elegir j} \frac{x^j}{z^j} \; dz
\\= \frac{1}{2\pi i}
\int_{|z|=\epsilon} (1+z)^{k-1}
\left(-1 + (1+z/x)^{k+1} \right) \; dz.$$
Toda la componente de a gotas, dejando
$$\frac{1}{2\pi i}
\int_{|z|=\epsilon} (1+z)^{k-1}
(1+z/x)^{k+1} \; dz
\\= \frac{1}{2\pi i}
\int_{|z|=\epsilon} \frac{(1+z)^{k-1}}{z^{k+1}}
(z+x)^{k+1} \; dz.$$
Sin embargo, esto es precisamente la integral que teníamos para la primera suma
sin el factor en el frente, de hecho.
El único infinito suma que aparece aquí es la primera con la convergencia
al $|(1+z)x|<1.$ por lo Tanto la elección de $|x|\lt 1/Q$ $|z|\lt 1/Q$
con $Q\ge 2$ tenemos $|(Q+1)/Q/Q|= |1/Q^2 + 1/Q|<1$ y obtener
la convergencia de la primera LHS integral en una vecindad de cero.
Un seguimiento de cuándo este método apareció en el MSE, y por el que se inicia en este
MSE enlace.