quad2dc

error:  int = quad2dc('Fun',xlow,xhigh,ylow,yhigh)
or
        int = quad2dc('Fun',xlow,xhigh,ylow,yhigh,tol)

This function is similar to QUAD or QUAD8 for 2-dimensional integration,
but it uses a Gaussian-Chebyshev quadrature integration scheme.  
     int     -- value of the integral
       Fun     -- Fun(x,y) (function to be integrated)
       xlow    -- lower x limit of integration  (should be -xhigh)
       xhigh   -- upper x limit of integration
       ylow    -- lower y limit of integration  (should be -yhigh)
       yhigh   -- upper y limit of integration
       tol     -- tolerance parameter (optional)
  The Gauss-Chebyshev Quadrature integrates an integral of the form
     yhigh                xhigh
  Int ((1/sqrt(1-y^2)) Int ((1/sqrt(1-x^2)) fun(x,y)) dx dy
    -yhigh               -xlow