prob_categories function

Probability for certain degree of saturation