\rho_P(r) = \bigg(\frac{3M}{4\pi a^3}\bigg)\bigg(1+\frac{r^2}{a^2}\bigg)^{-\frac{5}{2}}\,, .
where M is the total mass of the cluster, and a is the Plummer radius, a scale parameter which sets the size of the cluster core.