
function xc=cheb(a,b,n)

for m=1:n
    xc(m)=(a+b)/2-(b-a)/2*cos(pi*(m-1)/(n-1));
end
