k=39.5;
nuly=[exp(-i*pi*2/3) exp(-i*pi*5/6)]; % zadanie nul
pol=[0.45*exp(-i*pi*0.25) 0.65*exp(-i*pi*0.43)]; %zadanie polov

nuly_all=[nuly conj(nuly)]; nuly_poly=real(poly(nuly_all)); % dodanie komplexne zduzenych nul a vypocet koeficientov polynomu-citatela
pol_all=[pol conj(pol)]; pol_poly=real(poly(pol_all));% dodanie komplexne zduzenych polov a vypocet koeficientov polynomu- menovatela

ft=abs(fft(nuly_poly,1000))./(k*abs(fft(pol_poly,1000)));   % furierova tr.-> frekv. char.

plot([0:500]*pi/500,ft(1:501)); grid on% vykreslenie magnitudy do pi!