bivph function

Constructor function for bivariate phase-type distributions