Integrate[Integrate[x^2 Cos[x^2 - x y], {x, y, 1}], {y, 0, 1}]