Fastest way to integrate $\int_0^1 x^{2}\sqrt{x+x^2} \hspace{2mm}dx $ This integral looks simple, but it appears that its not so.
All Ideas are welcome, no Idea is bad, it may not work in this problem, but may be useful in some other case some other day ! :)

Let $x=\sinh^2 u$. (This is the same transformation as CountIblis used, but I'll employ it slightly differently.) Observe that $dx=2\cosh u \sinh u \, du$ and $$x+x^2=\sinh^2 u+\sinh^4 u=\sinh^2 u(1+\sinh^2 u)=\sinh^2 \cosh^2 u$$ since $\cosh^2 u-\sinh^2 u=1$.
Therefore \begin{align}
\int_0^1 x^2 \sqrt{x+x^2}\,dx
&=\int_0^{\sinh^{-1}1}\sinh^4 u\cdot \cosh u\sinh u \cdot \cosh u\sinh u\,du\\
&=\int_0^{\sinh^{-1} 1} \cosh^2 u \sinh^6 u\,du
\end{align}
To compute this integral, we recall that $\cosh u = \dfrac{1}{2}(e^u+e^{-u}),$ $\sinh u = \dfrac{1}{2}(e^u-e^{-u})$. One could expand the product in terms of exponentials and integrate term by term. For a bit less tedious route, first let $z=e^{-u}$ so that the integral takes the form
$$\int_{e^{-\sinh^{-1}(1)}}^{1} \frac{1}{2^8}\left(z+\frac{1}{z}\right)^2 \left(z-\frac{1}{z}\right)^6\frac{dz}{z}$$
We then expand this product into a sum of terms which can each be integrated by the power rule. The one obstacle is the weird-looking lower endpoint of $e^{-\sinh^{-1}1}$. However, one can show from $\sinh\sinh^{-1}(1)=1$ that this is equal to $\sqrt{2}-1$. That means that all that remains is the (admittedly rather tedious process) of integrating term-by-term and juggling square-roots.
