Numerical quadrature over a rectangular domain in two or more dimensions. I. Quadrature over a square, using up to sixteen equally spaced points

2. Integration over a Square. We restrict the domain of integration to be a square; this covers any rectangular domain by change of scale. We take the center of the square as origin of coordinates, and in the first place take the side of the square to be 2h; and the integral as h h