x=1:1933556; %obyvatelia
t=rand(1,1933556); t=t.^4; z=t*98800+1200; %mzdy
ht=[1200, 2000, 6000, 12000, 16000, 22000, 30000, 42000, 100000]; %hranice tried prijmov
pt=histc(z,ht); %vypise, kolko ludi patri do kazdej triedy
n=length(pt); pt=[pt(1:n-2),pt(n-1)+pt(n)]
%vypocita vysku stlpca ako podiel pocetnosti triedy jej sirkou 
ns=length(ht); st=ht(2:ns)-ht(1:ns-1) %vypocet sirky tried
vt=pt./st %vypocita spominanu vysku stlpca podielom pocetnosti a sirky 
bar(ht,[vt,0],'histc')
%vypisanie histogramu


   
