  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{c}$ pairs that are initially uncorrelated.
