Il est clair que :
$$(X_{1}+X_{2})(\Omega)=\llbracket0,n_{1}+n_{2}\rrbracket$$Posons $q=1-p$. Soit donc $k\in\llbracket0,n_{1}+n_{2}\rrbracket$. D'après le produit de convolution :
$$\begin{array}{rcl} \mathbb{P}(X_{1}+X_{2}=k) & = & \ds\sum_{\substack{i\in\llbracket0,n_{1}\rrbracket\\k-i\in\llbracket0,n_{2}\rrbracket}}^{n_{1}}{\mathbb{P}(X_{1}=i)\mathbb{P}(X_{2}=k-i)} \\ & = & \ds\sum_{i=\max\{0,j-n_{2}\}}^{\min\{j,n_{1}\}}{\left[\binom{n_{1}}{i}p^{i}q^{n_{1}-i}\binom{n_{2}}{k-i}p^{k-i}q^{n_{2}-k+i}\right]} \\ & = & \ds p^{k}q^{n_{1}+n_{2}-k}\sum_{i=\max\{0,j-n_{2}\}}^{\min\{j,n_{1}\}}{\binom{n_{1}}{i}\binom{n_{2}}{k-i}} \\ & = & \ds p^{k}q^{n_{1}+n_{2}-k}\binom{n_{1}+n_{2}}{k} \end{array}$$selon la formule de Vandermonde.