piFuns function

Compute multinomial cell probabilities