How do you evaluate int [x(x^2+1)]dx for [0, 1]?