 Solution

Convolve three rate-$\lambda$ exponential distributions. The first two have convolution density ${\lambda }^{2}x{e}^{-\lambda x}{1}_{x>0}$; convolving once more givesThus $X$ is an Erlang distribution of shape three, equal in law to the sum of three independent exponentials. Addition of independent means and variances yieldsFor $t\ge \lambda$ the moment-generating function is infinite.

Realize each holding time as a consecutive block of three exponential holding times from one Poisson process $N_t$. The renewal count is then $X_t=\lfloor N_t/3\rfloor$ pathwise under this coupling. If $R_t$ is the remainder of $N_t$ modulo three, $N_t=3X_t+R_t$ and $\mathbb{E}{R}_t=p_1(t)+2p_2(t)$. Taking expectations givesThis proves the formula without assuming that the renewal count itself is Poisson.