code: def quad2d(func=lambda x: 1, lower=(-10, -10), upper=(10, 10)):
      def fun(x, y):
          x = np.column_stack((x, y))
          return func(x)
  
      from scipy.integrate import dblquad
  
      return dblquad(
          fun,
          lower[0],
          upper[0],
          lambda y: lower[1],
          lambda y: upper[1],
      )
  
imports: import numpy
  from scipy import special
  from statsmodels.sandbox.distributions.multivariate import mvstdtprob
  from .extras import mvnormcdf
label_desc: Use image node_id 4 for calling a global function with example usage: quad2d(func, lower, upper) and returns: dblquad, func