Integrate[xy sqrt(x^2+y^2),{x,0,1},{y,0,1}]