A numerical study of Langevin QCD with dynamical quarks