I would likely begin by stating:
$$I(n,s-1)=\lim_{t\to\infty}\left(\int_0^t e^{-nx}x^{s-1}\,dx\right)$$
Okay, now, for the integral, I would use IBP, where:
$$u=x^{s-1}\,\therefore\,du=(s-1)x^{s-2}\,dx$$
$$dv=e^{-nx}\,dx\,\therefore\,v=-\frac{1}{n}e^{-nx}$$
And so we now have:
$$I(n,s-1)=\lim_{t\to\infty}\left(-\frac{1}{n}\left[x^{s-1}e^{-nx}\right]_0^t+\frac{s-1}{n}\int_0^t e^{-nx}x^{s-2}\,dx\right)$$
$$I(n,s-1)=\frac{s-1}{n}I(s-2,n)+\frac{1}{n}\lim_{t\to\infty}\left(\frac{x^{t-1}}{e^{nt}}\right)$$
We observe that the limit goes to zero (polynomial over an exponential) and we have:
$$I(n,s-1)=\frac{s-1}{n}I(n,s-2)$$
Now, if we repeat this process another $s-2$ times, we obtain:
$$I(n,s-1)=\frac{(s-1)!}{n^{s-1}}I(n,0)=\frac{(s-1)!}{n^{s-1}}\lim_{t\to\infty}\left(\int_0^t e^{-nx}\,dx\right)=\frac{(s-1)!}{n^s}$$