dphcdf function

Discrete phase-type cdf