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