a=ceil(6.*rand(1,100));
%wynik cytowany: a=[18 20 18 11 18 15];
sum(a)
stem(hist(a,1:6))
