We propose to use the Balescu-Lenard-Vlasov (BLV) equation to describe relativistic heavy ion collisions. We use an inter-quark Richardson's potential consistent with the indications of Lattice QCD calculations. The color degrees of freedom are explicitly taken into account. We explicitly demonstrate that the Vlasov approach alone is insufficient in the hadronization region. In order to overcome this problem we prepare the initial condition for many events using molecular dynamics with frictional cooling and a Thomas-Fermi approximation to the Fermi motion. These events are averaged and propagated in time using the Vlasov approach. The fluctuations are used to evaluate the collision term and in turn the number of q-$\bar q$ created.