integrate (x^3+1)*y, x from 0 to 1, y from 0 to x