X = zeros(30,1); for i=1:10000 D = sum( ceil( 6*rand(1,5) ) ); X(D) = X(D) + 1; end bar(X) xlabel('Die Total'); ylabel('Frequency'); title('10000 die rolls')