2D integration over non-rectangular domain using cubature

It’s quite natural, once you get used to it, to think of as “\int dx” as the operator “integrate the following function with respect to x”