Numerical simulations of quantum chromodynamics