En realidad, también tenemos
$$ \int_{0}^{\infty} t e^{-nct} (1-e^{-ct})^{m} \, dt = \frac{1}{c^2}\frac{(n-1)!m!}{(n+m)!}\sum_{k=0}^{m}\frac{1}{n+k}.$$
Esto es fácilmente demostrado por la integración por partes mediante la introducción de
$$ J(n, m) = c^{2} \frac{(n+m)!}{(n-1)!m!} \int_{0}^{\infty} t e^{-cnt} (1 - e^{-ct})^{m} \, dt$$
y darse cuenta de que $J(n, m) = \frac{1}{n} + J(n+1,m-1)$$J(n, 0) = \frac{1}{n}$.
También podemos invocar la función digamma $\psi_0$ a evaluar directamente. Deje $x = e^{-ct}$, de modo que
\begin{align*}
\int_{0}^{\infty} t e^{-nct} (1-e^{-ct})^{m} \, dt
&= - \frac{1}{c^2} \int_{0}^{1} x^{n-1} (1-x)^{m} \log x \, dx \\
&= - \frac{1}{c^2} \frac{\partial}{\partial n} \int_{0}^{1} x^{n-1} (1-x)^{m} \, dx \\
&= - \frac{1}{c^2} \frac{\partial}{\partial n} \frac{\Gamma(n)\Gamma(m+1)}{\Gamma(n+m+1)} \\
&= \frac{1}{c^2} \frac{\Gamma(n)\Gamma(m+1)}{\Gamma(n+m+1)} [ \psi_{0}(n+m+1) - \psi_{0}(n) ] \\
&= \frac{1}{c^2} \frac{(n-1)!m!}{(n+m)!} [ H_{n+m} - H_{n-1} ],
\end{align*}
cual es el resultado deseado.