EMstep_dph function

EM for discrete phase-type