>> n=9;				%zadefinujeme si n
>> A=ones(n);			%1tkova matica nxn
>> B=cumsum(tril(A));		%kumulativny sucet prvkov spodneho trojuhlnika matice	
>> C=triu(fliplr(flipud(B)),1);	%otoci maticu b hore a dolava+vyberie jej horny troluholnik o diagonalu vyssie
>> vysledok=B+C                 %spocita matice B a C a vypise vyslednu maticu



>> dataV=[161 164 164 164 167 168 169 171 173 175 175 176 177 178 178 180 180 181 181 182 183 184 185 186 188 190 193 197 198];
>> avg=mean(dataV);                 		%arit.priemer
>> x=round(randn(1,5000)*5+avg);   		%vytvori vektor nahodnych 5000 prvkov podla normalneho(gaussovho) rozdelenia 
>> hist(x)                          	     	%nakrelsi histogram
>> [p k]=hist(x)                    	     	%vypise hodnoty p=pocetnosti suboru v jednotlivych triedach    k=stredy tried
>> tr=[100 160 168 174 178 182 188 198 250]; 	%vektor tried
>> ptt=histc(x,tr)				%pocetnosti podla hranic jednotlivych tried
>> ptt=[ptt(1:end-3),ptt(end-2)+ptt(end)]       %zluci nulove triedy (podla toho co da treba upravit bud 3 alebo 2 toto je pre to ked budu na konci dve nuly)
>> n=length(x);					%dlzka x
>> rpt=ptt/n*100;				%relativne pocetnosti
>> dint=tr(2:end)-tr(1:end-1);			%vektor dlzok jednotlivych intervalov
>> vint=ptt./dint(1:end-1);			%vysky intervalov(zas treba davat pozor kolko nul je na konci a podla toho tu triedu v ktorej nie su prvky vyhodit)