p=chevypol(n)

function p=chevypol(n)

k=2;

p(1,1)=1;
p(2,1:2)=[1 0];
for k=2:n
   p(k+1,1:k+1)=2*[p(k,1:k) 0]-[0 0 p(k-1,1:k-1)];
end