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