rphasetype function

Simulate phase-type