function z = chebyzeros(n)

%z=chebyzeros(n) produces an n-by-1 vector containing the n  zeros of the Chebyshev polynomial of the first kind of order n. 

N = (1:n)';

%The Chebyshev polynomial of the first kind of degree n has n zeros in the interval [-1,1] given by: 
z = cos(pi * (N-0.5)/n);

%sort elements of z
z = sort(z);