We denote $S _ { r } = \sum _ { n = 1 } ^ { + \infty } \frac { H _ { n } } { ( n + 1 ) ^ { r } }$ for $r \geqslant 2$, and $B ( x ) = \int _ { 0 } ^ { 1 } ( \ln ( 1 - t ) ) ^ { 2 } t ^ { x - 1 } \mathrm {~d} t$. We have shown that $S _ { r } = \frac { ( - 1 ) ^ { r } } { 2 ( r - 2 ) ! } \int _ { 0 } ^ { 1 } \frac { ( \ln t ) ^ { r - 2 } ( \ln ( 1 - t ) ) ^ { 2 } } { t } \mathrm {~d} t$.
Deduce that for every integer $r \geqslant 2 , S _ { r } = \frac { ( - 1 ) ^ { r } } { 2 ( r - 2 ) ! } \lim _ { x \rightarrow 0 ^ { + } } B ^ { ( r - 2 ) } ( x )$.