function L6Bthreedice_comparison(n)
if nargin<1, n=100; end
  
subplot(1,2,1)
pexact = L6Bthreedice_exact;
bar(pexact)
title('Exact distribution')

subplot(1,2,2)
psamp = L6Bthreedice_sampling(n);
bar(psamp)
title(sprintf('Empirical from n=%d', n))