p=1<x*y<3 && x^2/2<y<2*x^2;

RegionPlot[p,{x,0,3},{y,0,3},PlotPoints->100]

Integrate[x^3*Boole[p],{x,0,Infinity},{y,0,Infinity}]