function [X]=chebyshev(a,b,N)
% return Chebyshev's interpolation points on the interval [a,b]
I = 1:N;
X=(b+a)/2 + (b-a)/2*cos((2*I-1)*pi/(2*N));