Is there a general formula for $$ I_n = \int_0^1 (x-1)(x-2)\dots (x-n) dx $$ in terms of familiar sequences of number (Harmonic numbers, Bell numbers, Bernoulli numbers, etc) ?
Define $$ P_n(x) = \int_0^x (u-1)(u-2)\dots (u-n) du $$ Have the polynomials $P_n(x)$ appeared in the literature? Is there a simple recursive formula for $P_n(x)$? What are some interesting properties of $P_n(x)$?
The explicit formula for $P_n(x)$ can be obtained by expanding $(u-1)(u-2)\dots(u-n)$ and integrating terms by terms $$ P_n(x) = \frac{x^{n+1}}{n+1} + \sum_{k=1}^n \frac{x^{n-k+1}}{n-k+1} \sum_{1 \leq j_1 < \dots < j_k \leq n} (-1)^k j_1 j_2 \dots j_k $$ Also, it is easy to see that $$ P_{n+1}'(x) = (x-n-1)P'_n(x) $$