Discrete phase-type distribution