function int = gauss_cheb(f,n,a,b);
x = cos(((2*(1:n) - 1)/(2*n))*pi);
w = pi/n;
int = sum(w.*f(x));