Obtaining cubatures for rectangles and other planar regions by using orthogonal polynomials