hh = 6;
h = 25; %pocet hladin, timto cislem vynasobim nejvyssi hodnotu a pak budu zaokrouhlovat, tim se vytvori hladiny blizke hodnoty se zaokrouhli na jednu hladinu
y = dlmread('HXR_'); % 6659
z = zeros(h*max(y) + 1, 1); 
u = h*y; 
v = round(u);
for k = 2:39999
    if (u(k) > u(k-1)) && (u(k+1) < u(k))
        if v(k)> 2
           z(v(k)) = z(v(k)) + 1;
        end
    end
end

% w = linspace(1,hh*h*round(max(y)), hh*h*round(max(y)));
figure('visible','off');
% figure(2);
hold on;
%plot(z);
m = smooth(z,10,'lowess'); 
ww=6*length(m);
w = linspace(1,6* length(m), length(m));
plot(w, m,'r');
print('-dpng','-r64','spectrum.png'); 
dlmwrite('smooth.txt',m,'delimiter','\t');
hold off;
exit;

% % tohle dělá úplně to samé!!
% figure(3);
% hist(y,200)
% axis([0,10, 0,  200])