Remarks

1. We can predict that given a function with two variables, there should be an optimal choice for picking the uniform regular matrices, a_nk, and b_nk.
2. We have the Richardson Extrapolation for the quadrature in one dimension as we mentioned earlier. We should be able to apply the Richardson Extrapolation for the two dimensional quadrature to speed up the convergence.