optim.probs function

performs the M step of the EM algorithm for the probability parameters