Given $X \sim \exp(\lambda)$ and $Y\sim\exp(\mu)$ with $\lambda \neq \mu$ I want to derive the cumulative distribution function of $Z := X + Y$.
My text does not use the convolution of probability density functions (in fact computing the CDF is a means to get to the PDF of $Z$).
It says, that for the CDF $F_Z$ the following holds: $$ F_Z(a) = P(Z < a) \stackrel{(*)}{=} \int_0^aP(Y < a - x)f_X(x)~dx, $$ where $f_X(x) = \lambda e^{-\lambda x}$ is the probability density function of $X$.
I have difficulty with $(*)$. I found a related question in the discrete case which I can understand but here it is in particular unclear where this product of a probability and the PDF of $X$ comes from.
Any insight and/or intermediate steps would be appreciated.