Abstract
This article presents the result of Langevin simulation studies of finite-temperature behavior of QCD for a various number of flavor species. Most of the simulations employ an 83×4 lattice. A full description is made of the data and the identification problem of a first-order phase transition. The systematic bias problem is also investigated.