bivdph function

Constructor function for bivariate discrete phase-type distributions