usage: int = quadc('Fun',xlow,xhigh)
or
int = quadc('Fun',xlow,xhigh,tol)
or
int = quadc('Fun',xlow,xhigh,tol,trace,p1,p2,....)
This function works just like QUAD or QUAD8 but uses a Gaussian-Chebyshev
quadrature integration scheme.
The Gauss-Chebyshev Quadrature integrates an integral of the form
xhigh
Int ((1/sqrt(1-x^2)) fun(x)) dx
-xhigh
This routine ignores xlow