multivariate phase type distribution