rdphasetype function

Simulate discrete phase-type