A partir de la definición de $P_k(n)$ tenemos
$$
P_k(n) = \sum_{i=1}^n i^k \Rightarrow
P_0(n) = \sum_{i=1}^n 1 = n.
$$
Ahora lo usamos para completar la inducción sobre $k$ a prueba la declaración de $S(k)$
$$
P_k(x) = k \int\limits_0^x P_{k-1}(t) \, dt + C_k \, x \quad (*)
$$
Caso Base
Para $k=1$ tenemos $S(1)$:
$$
1 \int\limits_0^x P_0(t) \, dt + C_1 \, x
=
\int\limits_0^x t \, dt + C_1 \, x
=
\frac{1}{2} x^2 + C_1 \, x
$$
que corresponde a la conocida Gauss suma fórmula $P_1(n)$ si $C_1 = 1/2$.
Paso inductivo
Suponiendo que la ecuación de $(*)$ es cierto para $\{ 1, \ldots, k \}$ realizamos la recursividad, comenzando con el $S(k)$ y, a continuación, aplicar el $S(k-1), S(k-2), \ldots$ hasta que llegamos a la parte inferior con $S(1)$ y consigue $P_0 = \mbox{id}$. Que integrando luego es lo suficientemente simple para realizar el la $k$ a afirmar integraciones:
\begin{align}
P_k(x)
&=
k \int\limits_0^x P_{k-1}(t_k) \, dt_k + C_k \, x \\
&=
k \int\limits_0^x \left( (k-1) \int\limits_0^{t_k}
P_{k-2}(t_{k-1}) \, dt_{k-1} + C_{k-1} \, t_k \right) \, dt_k + C_k \, x \\
&=
k (k-1) \int\limits_0^x \int\limits_0^{t_k}
P_{k-2}(t_{k-1}) \, dt_{k-1} \, dt_k + \frac{k}{2} C_{k-1} \, x^2 + C_k \, x \\
&=
k! \int\limits_0^x \cdots \int\limits_0^{t_2}
P_0(t_1) \, dt_1 \cdots \,dt_k
+ \sum_{j=1}^k C_{k-j+1} \frac{k!}{(k-j+1)!j!} x^j \\
&=
k! \int\limits_0^x \cdots \int\limits_0^{t_2}
t_1 \, dt_1 \cdots \,dt_k
+ \sum_{j=1}^k \binom{k}{j} \frac{C_{k-j+1}}{k-j+1} x^j \\
&=
k! \int\limits_0^x \cdots \int\limits_0^{t_3}
\frac{1}{2}t_2^2 \, dt_2 \cdots \,dt_k
+ \sum_{j=1}^k \binom{k}{j} \frac{C_{k-j+1}}{k-j+1} x^j \\
&=
k! \int\limits_0^x
\frac{1}{k!}t_k^k \,dt_k
+ \sum_{j=1}^k \binom{k}{j} \frac{C_{k-j+1}}{k-j+1} x^j \\
&=
\frac{1}{k+1}x^{k+1}
+ \sum_{j=1}^k \binom{k}{j} \frac{C_{k-j+1}}{k-j+1} x^j \\
\end{align}
A continuación, tratamos de llegar a $S(k+1)$:
\begin{align}
(k+1) \int\limits_0^x P_k(t) \, dt + C_{k+1} \, x
&=
(k+1) \int\limits_0^x \left(
\frac{1}{k+1}t^{k+1} +
\sum_{j=1}^k \binom{k}{j} \frac{C_{k-j+1}}{k-j+1} t^j \right)
\, dt \\
& + C_{k+1} \, x \\
&=
(k+1) \left(
\frac{1}{(k+1)(k+2)}x^{k+2}
+ \right. \\
& \left. \sum_{j=1}^k \binom{k}{j} \frac{C_{k-j+1}}{(k-j+1)(j+1)} x^{j+1}
\right) + C_{k+1} \, x \\
&=
\frac{1}{k+2}x^{k+2}
+ \sum_{j=1}^k \binom{k}{j} \frac{(k+1)C_{k-j+1}}{(k-j+1)(j+1)} x^{j+1}
+ C_{k+1} \, x \\
&=
\frac{1}{k+2}x^{k+2}
+ \sum_{j=1}^k \binom{k+1}{j+1} \frac{C_{k-j+1}}{k-j+1} x^{j+1}
+ C_{k+1} \, x \\
&=
\frac{1}{k+2}x^{k+2}
+ \sum_{j=2}^{k+1} \binom{k+1}{j} \frac{C_{k+1-j+1}}{k+1-j+1} x^j
+ C_{k+1} \, x \\
&=
\frac{1}{k+2}x^{k+2}
+ \sum_{j=1}^{k+1} \binom{k+1}{j} \frac{C_{k+1-j+1}}{k+1-j+1} x^j \\
&=
P_{k+1}(x)
\end{align}
Por lo tanto $S(k+1)$ sigue.
Por el principio de inducción $(*)$ tiene para todos los $k \in \mathbb{N} \setminus \{ 0\}$.