I don't have a general expression as a function of $n$, but there is one in terms of harmonic numbers $H_k$ as a function of $k\in\mathbb{R}^+$ for given $n$; for example
$$I_{2,k}=\tfrac{1}{2}\Gamma(k+1)\left(\frac{2 H_{k+3}-3}{ \Gamma(k+4)}+\frac{6\left(25-12 H_{k+5}\right)}{\Gamma (k+6)}\right)$$ $$I_{3,k}=3 \Gamma (k+1) \left(\frac{12 H_{k+5}-25}{\Gamma (k+6)}+\frac{30 \left(49-20 H_{k+7}\right)}{\Gamma (k+8)}\right)$$ $$I_{4,k}=\tfrac{3}{4} \Gamma (k+1) \left(\frac{25-12 H_{k+5}}{\Gamma (k+6)}+\frac{180 \left(20 H_{k+7}-49\right)}{\Gamma (k+8)}+\frac{840 \left(761-280 H_{k+9}\right)}{\Gamma (k+10)}\right)$$ $$I_{5,k}=\tfrac{45}{2} \Gamma (k+1) \left(\frac{147-60 H_{k+7}}{\Gamma (k+8)}+\frac{56 \left(280 H_{k+9}-761\right)}{\Gamma (k+10)}-\frac{504 \left(2520 H_{k+11}-7381\right)}{\Gamma (k+12)}\right)$$
Carlo Beenakker
- 188.1k
- 18
- 448
- 651