Economical symmetrical quadrature rules for complete polynomials over a square domain

It is of interest in numerical analysis to develop symmetrical quadrature rules for integration of complete polynomial functions over a square domain with minimum computational effort. Gaussian product quadrature rules integrate such functions with maximum effort. Symmetrical quadrature rules are developed and presented for integration of complete polynomial functions up to 21st order with minimum computational effort.