  We compute by Monte Carlo methods the individual distributions of the $k$th
smallest Dirac operator eigenvalues in QCD, and compare them with recent
analytical predictions. We do this for both massless and massive quarks in an
SU(3) gauge theory with staggered fermions. Very precise agreement is found in
all cases. As a simple by-product we also extract the microscopic spectral
density of the Dirac operator in SU(3) gauge theory with dynamical massive
fermions for $N_f=1$ and 2, and obtain high-accuracy agreement with analytical
expressions.
