usage: 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)
