mdph function

Constructor function for multivariate discrete phase-type distributions