rMPHstar function

Simulate a MPH* random vector