tic
rep = 1e5;
nsuccess = 0;

for i=1:rep
    emp = hist(dice(12), 1:6);
    if all(emp == 2)
        nsuccess = nsuccess + 1;
    end
end
toc

fprintf('%10d trials\n', rep);
fprintf('%10d successes (all bars height 2)\n', nsuccess);
fprintf('%10.5f success rate\n', nsuccess/rep);