A = zeros(100,1);
for n_tries = 1:100,
    numberExperiments = 1000;
    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
        A(n_tries) = A(n_tries) + k;        % Add each value of k
    end
    A(n_tries) = A(n_tries)/numberExperiments; % Divide by 1000 to get the average of k 
end

hist(A)                                     % Plot histogram of average values