In the Zeta-urn model, given $M$ urns and $N$ balls, the probability distribution at equilibrium is: \begin{equation} P(n_i)=Z^{-1}(N,M,\beta)e^{-\beta\sum_{i}E_i}, \end{equation} with $E_i=\ln(n_i+1)$, $\beta$ a real parameter and $Z$ a normalization constant.

Now, does anyone know a simple way to implement a simulation of this system in order to obtain the distribution of balls at equilibrium, given $\beta$ and $\rho$?

The question is about the dynamic of the system, i.e. how to define a rule to fill the urns in order to reach an equilibrium state that satisfy the probability distribution.


closed as off-topic by Qmechanic Mar 15 '17 at 20:01

  • This question does not appear to be about physics within the scope defined in the help center.
If this question can be reworded to fit the rules in the help center, please edit the question.

  • $\begingroup$ I'm voting to close this question as off-topic e.g. because it seems to be about statistics & software implementation rather than physics. Would Computational Science be a better home for this question? $\endgroup$ – Qmechanic Mar 15 '17 at 20:01
  • $\begingroup$ @Szabolcs: You are encouraged to improve the post via edits if possible. It still seems to ask for software implementation, though. $\endgroup$ – Qmechanic Mar 15 '17 at 21:12
  • $\begingroup$ @Qmechanic Yes, sorry, I deleted my comment before you responded because I realized it was pointless. I am not arguing for reopening. $\endgroup$ – Szabolcs Mar 15 '17 at 21:18