mph function

Constructor function for multivariate phase-type distributions