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