In[]:=
Integrate[xy*sqrt(x^2+y^2+z^2),Element[{x,y,z},ImplicitRegion[z>=sqrt(x^2+y^2)&&x^2+y^2+z^2<=1,{x,y,z}]]]
​
​
​
​
​
​
​
​
​