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