tvr_dph function

Performs TVR for discrete phase-type distributions