This Markov Chain Monte Carlo simulation converges at large enough times to the uniform distribution over partitions of the number $N$. The limit shape of the distributions for large $N$ is $e^{-x/a}+e^{-y/a}=1$, where $a=\sqrt{6N}/\pi$.

N = redraws =

MC steps per redraw =

show limit shape show means

This work is licensed under a Creative Commons Attribution 3.0 United States License