Goodman_pi_matrix function

Computes the full matrix of model-based cell probabilities.