compute_choice_probabilities function

Compute probit choice probabilities