Solution (source code)

= Solution

At $q=2$, the preceding <tensor-product quadrature> gives
$$
\boxed{Q(g)=\frac{1+1}{36}=\frac1{18}=I(g).}
$$
Both factors of the <separable integrand> are <quadratic polynomials>, so exactness follows directly from the construction of <Simpson's rule>.