function [x,LP]=KLP(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 legendrovy 
%polynomy, vyjádřené vektorem koeficientů
LP=zeros(m+1,m+1);
%tvorba matice do jejíž řádků se budou ukládat koeficienty polynomů 
LP(1,m+1)=1; %polynom nultého stupně
LP(2,m)=1; %polynom prvního stupně
i=1;
	while i<=m-1 %cyklus pro výpočet polynomů vyšších řádů
		r=[(2*i+1)/(i+1) 0];
		p1=conv(r,LP(i+1,:));
		p1(1)=[];
		p0=i/(i+1)*LP(i,:);
		LP(i+2,:)= p1 - p0;
		i=i+1; 
	end
x=roots(LP(m+1,:)); %výpočet kořenů posledního polynomu
end