Behavior of vortices is core to the understanding of phenomenons in nature, such as Jupiter's red spot. The simplification to point vortices provides us with an easier way to model behaviors of large numbers of vortices. Although we know from statistical mechanics that there exists a time invariant measure, it is unclear whether the measure derived from the Mean Field Equation (MFE) reflects the true behavior in nature. We explore the point vortex model on the 2 dimensional disk and test for the convergence to the Mean Field distribution using Kolmogorov-Smirinov test (K-S test). Due to the pairwise interaction of the vortices, we use the Fast Multipole Method (FMM) to decrease the calculations required. We hope that the observed behavior on the disk imply a general result on all simply connected domains via conformal mappings.