dmodifiedCount_bi function

Compute count probabilities based on modified renewal process (bi)