`int xsqrt(x+x^2) dx`

Explanation

`1/3 (x^2+x)^(3//2)-1/8(2x+1)sqrt(x^2 +x) -1/16 log(x+1/2+sqrt(x^2 +x))+c`