% Simulation of Drawing an Ace
%
n = 10000;                % Number of draws.
X = ceil(52*rand(1,n));
aces = (1 <= X & X <= 4);
naces = sum(aces);
fprintf('There were %g aces in %g draws.\n',naces,n)