Integrate[ x y Sqrt[x^2 + y^2], {x, 0, 1}, {y, 0, 1}]