EMstep_bivph function

EM for bivariate phase-type distributions using Pade for matrix exponential