Monte Carlo for Estimating Exponential Convolution

We study the numerical stability problem that may take place when calculating the cumulative distribution function of the Hypoexponential random variable. This computation is extensively used during the execution of Monte Carlo network reliability estimation algorithms. In spite of the fact that analytical formulas are available, they can be unstable in practice. This instability occurs frequently when estimating very small failure probabilities (10−30 − 10−40) that can happen for example while estimating the unreliability of telecommunication systems. In order to address this problem, we propose a simple unbiased estimation algorithm that is capable of handling a large number of variables.