Function that computes the area's integral for polygonal domains

[carlos-adir]

Create a function that allows computing the integral:

$$I_{a,b} = \int_{D} x^a \cdot y^b \ dx \ dy$$

For polygonal domains