Langevin Simulation Including Dynamical Quark Loops

Abstract
A Langevin method is proposed for simulation of full QCD including dynamical quark loops. It is shown that the method works well once a proper discretization of the fictitious time is made. A realistic simulation with this method is perfectly feasible on vector computers presently available.