sum_dph function

Computes the initial distribution and sub-intensity of the sum of two discrete phase-type distributed random variables