EMstep_RK function

EM step for phase-type using Runge-Kutta