sim-ph-method function

Simulation method for phase-type distributions