quad2dcgen

error:  int = quad2dcgen('Fun','funxlow','funxhigh',ylow,yhigh)
or
        int = quad2dcgen('Fun','funxlow','funxhigh',ylow,yhigh,tol)

This function is similar to QUAD or QUAD8 for 2-dimensional integration
over a general 2-dimensional region, but it uses a Gauss-Chebyshev 
quadrature integration scheme.  
The integral is like:
               yhigh                   funxhigh(y)
      int = Int  (1/sqrt(1-y.^2))   Int  (1/sqrt(1-x.^2))  Fun(x,y)  dx  dy
               ylow                    funxlow(y)

     int     -- value of the integral
       Fun     -- Fun(x,y) (function to be integrated)
       funxlow -- funxlow(y)
       funxhigh-- funxhigh(y)
       ylow    -- lower y limit of integration
       yhigh   -- upper y limit of integration
       tol     -- tolerance parameter (optional)