numberExperiments = 1000;
bins = zeros(6,1); 

for count = 1:numberExperiments
%   Generate a randomn number k with value 1, 2, 3, 4, 5, or 6, with
%   probability 1/6.  
    r = rand;
    if (r < 1.0/6)
     k = 1;
    elseif (r < 2.0/6)
     k = 2;
    elseif (r < 3.0/6)
     k = 3;
    elseif (r < 4.0/6)
     k = 4;
    elseif (r < 5.0/6)
     k = 5;
    else
     k = 6;
    end
    bins(k) = bins(k) + 1;      % Count times each number generated
end

bar(bins)                       % Plot the bins