Problem 2. Consider a double integral ∬f(x,y)dxdy on a rectangular domain Ω = [a, b] × [c, d]. Based on its conversion to an iterated integral, develop a Gaussian quadrature for the double integral. You may use the same 1-dimensional Gaussian quadrature in the x- and y-directions.