>> n=6;                                	%zadame n                  
>> A=ones(2*n);                        	%vytvori jednotkovu maticu rozmerov 2n x 2n 
>> neparne=1:2:2*n;                    	%vytvori vektor neparnych cisiel 1:2n
>> parne=2:2:2*n;                      	%vytvori vektor parnych cisiel 2:2n
>> A(:,parne)=0;                       	%vynuluje neparne stlpce
>> A(parne,:)=0;                       	%vynuluje parne riadky
>> Q=[1.5:0.5:(n+1);ones(n-1,2*n)];    	%pomocna matica vrchny riadok od 1.5 po n+1 s krokom 0.5 zvysok jednotky
>> P=cumsum(Q);				%spravi kumulativny sucet prvkov matice Q
>> P(:,neparne)=0;			%vynuluje neparne stlpce matice P
>> A(2:2:end,:)=P(1:end,:)              %na parne riadky matice A vlozi riadky matice P a vypise vyslednu maticu A



>> dataV=[161 164 164 164 167 168 169 171 173 175 175 176 177 178 178 180 180 181 181 182 183 184 184 185 186 188 190 193 197 198];  %definicia datoveho suboru
>> avg=mean(dataH); 				%priemer dataH
>> x=round(randn(1,5000)*5+avg)		  	%vytvori vektor nahodnych 5000 prvkov podla normalneho(gaussovho) rozdelenia
>> poc=9;			       	        %definicia poctu tried
>> sirka=(max(x)-min(x))/poc;                   %sirka jednej triedy
>> hr=[min(x):sirka:max(x)];                    %zadefinuje hranice tried
>> pt=histc(x,hr);				%pocetnosti v triedach hr
>> pt=[pt(1:end-2),pt(end-1)+pt(end)];		%pripocita poslednu triedu k predposlednej
>> cl=[100 160 168  174 178 182 188 198 250];	%definicia tried
>> ptt=histc(x,cl)				%pocetnosti v triedach cl
>> ptt=[ptt(2:end-3),ptt(end-2)+ptt(end)]   	%odstrani nulove triedy(zas pozor na cisla podla toho co da povodne ptt)
>> n=length(x);					%dlzka x
>> rpt=ptt/n*100;				%relativne pocetnosti
>> dint=cl(2:end)-cl(1:end-1);			%vektor dlzok jednotlivych intervalov
>> vint=ptt./dint(2:end-1);			%vysky intervalov(zas treba davat pozor kolko nul je na konci a podla toho tu triedu v ktorej nie su prvky vyhodit)