% fair die demo

N = 256;
disp(['Running ', num2str(N),' random trials']);
x = ceil(6*rand(N,1));
figure(1)
plot(x,'+');
figure(2)
bar(hist(x,[1:6])./N);