tvr_ph function

Performs TVR for phase-type distributions