function E = eql_curves(fs, n)

% E = eql_curves(fs, n)
% 
% Funkcia počíta krivky rovnakej hlasitosti.
% 
% fs    - vzorkovacia frekvencia
% n     - počet vzorkov na krivku
% E     - výstupný vektor obsahujúci krivky

%% Výpočet frekvencií
f = linspace(1, ceil(fs/2) ,n);
o = 2*pi*f;

%% Výpočet kriviek
E = ((o.^4).*((o.^2) + 56.9e6))./((((o.^2) + 6.3e6).^2).* ...
    ((o.^2) + 379.4e6).*((o.^6) + 9.6e26));
E = E./(max(E));