We present a dynamical model of heavy quark evolution in the quark–gluon plasma (QGP) based on the Fokker–Planck equation. We then apply this model to the case of central ultra-relativistic nucleus–nucleus collisions performed at RHIC and estimate the component of $J/\Psi$ production (integrated and differential) stemming from $c–\bar{\rm c}$ pairs that are initially uncorrelated.