from
ChebyshevTools
by Michael Watson Tool box for solving ODE/PDEs using spectral Chebyshev differentiation matrices.
x=cheb_grid(N);
function x=cheb_grid(N);
% this function returns a chebyshev grid on the interval [-1,1]
% with entry x(1)=1 and x(end)=-1
x=cos(pi.*linspace(0,1,N))';