optim.const.ordi function

performs the M step for the measurement distribution parameters in multinomial case, with an ordinal constraint on the parameters