% Lecture 2A, example 1: Estimating P(X<=2) in a six-sided die
%
for n=[100 10000 1e6 1e8]
    fprintf('n=%-8g ...  ', n);
    tic;
    X = randi(6, 1, n);
    I = (X <= 2);
    elapsed = toc;
    relfreq = sum(I)/n;
    fprintf('relfreq=%10.8f    time=%.3f s\n',...
        relfreq, elapsed);
end