sum_ph function

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