No es elemental (en el sentido de no utilizar análisis complejos), pero así es como yo lo haría:
Dejemos que $$f(z) = \frac{1}{z^s(z+1)},$$ donde $z^s$ denota la "rama natural", es decir, elegir $\phi \arg z \in (0,2\pi)$ y poner $(re^{i\phi})^s = r^s e^{is\phi}$ . Toma un "contorno de ojo de cerradura" C:
e integrar $f$ a lo largo de $C$ utilizando el teorema del residuo:
$$\int_C f(z)\,dz = 2\pi i \operatorname{Res}_{z=-1} f(z).$$
Estimación de $\int_\gamma f(z)\,dz$ y $\int_\Gamma f(z)\,dz$ tenemos
$$\left| \int_\gamma f(z)\,dz \right| \le \frac{M}{r^{\mathrm{Re}(s)}} \cdot 2\pi r \to 0 \qquad\text{as }r\to0$$ y $$\left| \int_\Gamma f(z)\,dz \right| \le \frac{M}{R^{1+\mathrm{Re}(s)}} \cdot 2\pi R \to 0 \qquad\text{as }R\to\infty.$$
(Para la primera estimación, queremos $\mathrm{Re}(s) < 1$ y para el segundo $\mathrm{Re}(s) > 0$ .
Para los dos segmentos de línea restantes, en el segmento "superior" obtenemos la integral que buscamos como $r \to 0$ y $R\to\infty$ . En el segmento "inferior", obtenemos (recuerde la elección de la rama) $$-\int_0^\infty \frac{1}{x^s e^{2\pi i s} (1+x)}\,dx.$$
Ponerlo todo junto: $$(1-e^{-2\pi i s}) \int_0^\infty \frac{1}{x^s(1+x)}\,dx = 2\pi i \operatorname{Res}_{z=-1} f(z) = 2\pi i (-1)^{-s} = 2\pi i e^{-\pi i s}$$ así que $$ \int_0^\infty \frac{1}{x^s(1+x)}\,dx = 2\pi i \frac{e^{-\pi i s}}{1-e^{-2\pi i s}} = \frac{\pi}{\sin s\pi}.$$