function [x,CP]=KCP(m)
%Vstup
%m - počet kořenů které potřebujeme vypočítat,
%neboli stupeň polynomu
%Výstupy
%x - sloupcový vektor kořenů m-tého polynomu
%LP - čtvercová matice o velikosti m+1, řádky matice jsou čebyševovy 
%polynomy vyjádřené vektorem koeficientů
CP=zeros(m+1,m+1);
%tvorba matice do jejíž řádků se budou ukládat koeficienty polynomů 
%vyjádřené vektorem
CP(1,m+1)=1; %polynom nultého stupně
CP(2,m)=1; %polynom druhého stupně
i=1;
	while (i<=m-1) %cyklus pro výpočet polynomů vyšších řádů
		r=[2 0];
		p0=conv(r,CP(i+1,:));
		p0(1)=[];
		CP(i+2,:)=p0 - CP(i,:);
		i=i+1; 
	end
x=roots(CP(m+1,:)); %výpočet kořenů posledního polynomu
end