integrate (1/2) (1+xy(x^2-y^2)) dy, x,1