What is the distribution (probability density function) of the difference of two Gammas with the same scale? The parameters of each Gamma would be positive integers.
That is, I'd like to know the density function ($f$) of $Z=X-Y$, where $X\sim {\rm Gamma}(a_1,b_1)$ and $Y\sim {\rm Gamma}(a_2, b_1)$. Assume $c$, $d$ are the parameters of $f(Z)$, particularly, I want to know the relationships between $c$, $d$, and $a_1$, $a_2$, $b_1$.