rMDPHstar function

Simulate MDPH*