I need to integrate the product of two polynomial functions defined on an arbitrary (convex) planar quadrilateral defined by 4 points in $\mathbb{R}^3$.
I was trying to firstly rotate the system of reference so to make the problem 2D and then transform the quadrilateral into the unit square and perform the integration there with the correct change of coordinates. This is the reason I opened this question.
However some users suggested different strategies, so I move the question to a higher level.
Is there any theorem or formula to perform such integration?
Thank you very much