Let \( f \) be a polynomial function such that \( f(x^{2}+1)=x^{4}+5x^{2}+2 \), for all \( x \in \mathbb{R} \). Then \( \displaystyle \int_{0}^{3} f(x) \, dx \) is equal to:
\( \dfrac{27}{2} \)
\( \dfrac{41}{3} \)
\( \dfrac{5}{3} \)
\( \dfrac{33}{2} \)