How to calculate $P(\sum_{i=1}^{m}(A_i+S_i)\le L)$ with $A_i,L\sim\text{exp}(\lambda),S_i\sim\text{exp}(\mu)$ and positive integers $\lambda\neq\mu$?

Question

Recently I was stumped by the calculation of the probability $$\mathbb{P} \big(\sum_{i=1}^{m} (A_i + S_i) \le L < \sum_{i=1}^{m+1} (A_i + S_i) \big)$$ where $A_i \sim \text{exp}(\lambda), S_i \sim \text{exp}(\mu), L \sim \text{exp}(\lambda)$ are mutually independent, $\lambda \neq \mu$ are two positive integers, and $m$ is an integer parameter.

---

My attempt:

I failed at the very start while calculating the density function of $\sum_{i=1}^{m} (A_i + S_i)$: $$\sum_{i=1}^{m} (A_i + S_i) = \sum_{i=1}^{m}{A_i} + \sum_{i=1}^{m}{S_i} = \big( A \sim \Gamma(m, \lambda) \big) + \big( S \sim \Gamma(m, \mu) \big)$$ Considering the two Gamma random variable $A$ and $S$, their convolution is: $$p_{A+S}(a) = p_{A} \ast p_{S} (a) = \int_{0}^{a} f_{A}(a-y) f_{S}(y) dy \\ = \int_{0}^{a} \frac{\lambda e^{-\lambda (a-y)} (\lambda (a-y))^{m-1}}{\Gamma(m)} \frac{\mu e^{-\mu y} (\mu y)^{m-1}}{\Gamma(m)} dy \\ = e^{-\lambda a} \frac{\lambda^m \mu^m}{\Gamma(m) \Gamma(m)} \int_{0}^{a} e^{(\lambda - \mu) y} \large((a-y)y \large)^{m-1} dy$$

---

And my failure:

I was not able to go any further with the integral. And you can image what a messy calculation we will run into when we are going to take the second step involving $\sum_{i=1}^{m} (A_i + S_i) \le L$.

---

My questions:

Generally, my question is how to make the calculation tractable.

Specifically,

1. Has the probability been studied in research papers?
   People interested in its background in "alternating renewal process" are referred to the following part of the post.

2. What is the density function of $\big( A \sim \Gamma(m, \lambda) \big) + \big( S \sim \Gamma(m, \mu) \big)$ given that $m$ is an integer and $\lambda \neq \mu$ are positive integers?
   If there is no nice closed form, can we make some good approximation?
   Related post: at Math.SE (without proper answers yet).

3. For the whole calculation, is it feasible to obtain numerical solution using mathematical softwares?
   Are there any off-the-shelf libraries for this purpose?

---

---

For people who want to know the background of the probability in "alternating renewal process":

Background: Consider an alternating renewal process, in which a system can be in one of two states: on or off. Whenever it is off, it takes a time $A_i \sim \text{exp}(\lambda)$ before turning on. Whenever it becomes on, it remains on for a time $S_i \sim \text{exp}(\mu)$. Initially the system is off. A cycle $c_i$ consists of the $i$th off-state and the $i$th on-state. All $A_i$ and $S_i$ are mutually independent and $\lambda \neq \mu$ are two positive integers.

Problem: Given a time period $L \sim \text{exp}(\lambda)$, what is the probability that exactly $m$ cycles occur in $L$?
Related post: at MO: its origin; many thanks to @user137846.

Answer 1

Call the left resp. right hand sum $R_m$ resp. $R_{m+1}$. As $L$ is $\exp(\lambda)$ and independent of $(R_m,R_{m+1})$ , taking expection with resp. to $L$ first gives $$\mathbb{P}(R_m\leq L < R_{m+1})=\mathbb{E}(e^{-\lambda R_m}-e^{-\lambda R_{m+1}})= \left(\dfrac{1}{2}\dfrac{\mu}{\mu+\lambda}\right)^m\left(1-\dfrac{1}{2}\dfrac{\mu}{\mu+\lambda}\right)$$ where we have used that the individual summands of $R_m$ resp. $R_{m+1}$ have Laplace-transform $$\mathbf{E}(e^{-p(S_i+A_i)})=\dfrac{\lambda}{\lambda+p}\dfrac{\mu}{\mu+p}$$

Answer 2

You are looking at summation of independent gamma random variables $\sum_i X_i$, in the hard case where the $X_i$ have different scale parameters. It was studied in:

P. G. Moschopoulos (1985) The distribution of the sum of independent gamma random variables, Annals of the Institute of Statistical Mathematics, 37, 541-544

A. M. Mathai (1982) Storage capacity of a dam with gamma type inputs, Annals of the Institute of Statistical Mathematics, 34, 591-597