Solution (source code)

= Solution
{c}

Independence and convolution give
$$
\mathbb P(X+Y=k)=e^{-(\mu+\lambda)}\sum_{j=0}^k\frac{\mu^j\lambda^{k-j}}{j!(k-j)!}
=e^{-(\mu+\lambda)}\frac{(\mu+\lambda)^k}{k!},
$$
by the <binomial theorem>. Thus
$$
\boxed{X+Y\sim\operatorname{Poisson}(\mu+\lambda).}
$$