If $X_1, X_2$ and $X_3$ are i.i.d Gamma($\theta$, 1) random variables, and $\underline{\mathbf{Y}}$ = ($X_1$ + $X_3$, $X_2$ + $X_3$) then what will be the joint density function of $\underline{\mathbf{Y}}$?
I started by defining three new random variables as $$Z_1 = X_1 + X_3,\ Z_2 = X_2 + X_3,\ Z_3 = X_3$$ and found the joint density function of $(Z_1, Z_2, Z_3)$; i.e. $f(z_1,z_2,z_3)$, where $z_1 \gt 0, z_2 \gt 0$, and $0 \lt z_3 \lt \min\{z_1, z_2\}$. After that, I obtained the joint density of $Z_1$ and $Z_2$ as $f(z_1,z_2)$ by integrating out $z_3$ over its given support. But I am unable to solve the integral, which is $$\int_{0}^{\min\{z_1, z_2\}} e^{z_3} (1 - z_3/z_1)^{\theta - 1} (1 - z_3/z_2)^{\theta - 1} z_3^{\theta - 1} dz_3.$$ Can anyone please help me in simplifying this integral?