function h = kostka3(N)
%
h = zeros(1,6);
r = randi(6,1,N);
for i=1:N
  idx = r(i);
  h(idx) = h(idx)+1;
end
