# I miss Mathematica

Why? Because Mathics is not up to helping me determine if indeed
$f(\{A_1, \ldots, A_p\}) = \frac{(n-p)!^2}{n! p!} \left(\frac{1}{p} \right)^{n-p} \frac{|A_1|^2 \cdots |A_p|^2}{|A_1|!\cdots |A_p|!}$
is a pmf over the set of partitions of the set $$\{1, \ldots, n\}$$ into $$p \leq n$$ nonempty sets. In particular, the sets $$A_1, \ldots, A_p$$ in the formula above satisfy $$|A_1| + \ldots +|A_p| = n$$ and $$|A_i| \geq 1$$ for all $$i.$$

I feel like I could empirically test this easily in Mathematica, but OMG trying to do it in Matlab is a real pain, so I gave up. Combinatorics or set manipulation in Matlab in general is an exercise in trying to make a smoothie with a grater: you can do it, eventually, but it’s going to take forever and make a mess.