a couple small points: if OP wants to clean up the denominator, setting up some inequalities like the below should sandwich the result.
##\frac {\int_0^t (\arctan x) ^2dx}{t}\frac{t}{t+1}= \frac {\int_0^t (\arctan x) ^2 dx}{(t+1)}= \frac {\int_0^t (\arctan x) ^2 dx}{\sqrt {(t+1)^2}}= \frac {\int_0^t (\arctan x) ^2 dx}{\sqrt {t^2 + 2t +1}} \leq \frac {\int_0^t (\arctan x) ^2 dx}{\sqrt {t^2+1}} \leq \frac {\int_0^t (\arctan x) ^2 dx}{\sqrt {t^2}} = \frac {\int_0^t (\arctan x) ^2 dx}{t}##also the fact that ##\arctan^2 x## has a limit implies that the Cesaro mean evaluates to said limit. We are dealing with integrals -- specifically ##\frac{1}{t} \int_0^t (\arctan x) ^2 dx##
not series here, but it's a closely related point -- and we could work with
##\frac{1}{t} \sum_{x=1}^t (\arctan x) ^2##
if we were so inclined. Normally its the case that people use integrals to estimate and bound series, but amusingly, I think the above integrals can be sandwiched by two cesaro means, to give the result. Given monotonicity / negative convexity for the function which exists for all ##x \geq 1##, it shouldn't be too hard