Joint distribution of two functions of two random variables 1) Suppose I have two random variables $A > 0$ and $B > 0$ with joint p.d.f. $f_{A,B}(a,b)$, and two random variables $X = g_1(A,B)$ and $Y = g_2(A,B)$. What is the general procedure for determining the joint p.d.f. $f_{X,Y}(x,y)$ ?
2) More specifically, suppose that:


*

*$A$ and $B$ are i.i.d. with $A \sim \Gamma (k, m)$.

*$X = AB$

*$Y = A + B$


Clearly it this case $X$ takes a K distribution and $Y$ takes a Gamma distribution.  Due to struggling with 1) above, I've unfortunately been unable to determine whether in this case a closed non-integral expression for the joint p.d.f. $f_{X,Y}(x,y|k,m)$ even exists, let alone what it is...
Any assistance / suggestions would be greatly appreciated.
 A: Let me first state the result:
$$
   f_{X,Y|k,m}(x,y) = \frac{2}{\sqrt{y^2-4 x}} \frac{x^{k-1} e^{-\frac{y}{m}}}{(k-1)!^2 m^{2 k} } I(x>0,y>0,y^2 > 4 x)
$$
Many posts are devoted to the subject of determining the probability density function for the function of random variables, i.e. this, that, and other posts referenced therein.
The mapping $(A,B) \mapsto (A B, A+B)$ is unto. Specifically, for 
$$
   x= a b, y = a+b \implies a = \frac{1}{2} \left( y \pm \sqrt{y^2-4x}\right), b =  \frac{1}{2} \left( y \mp \sqrt{y^2-4x}\right) 
$$
The Jacobians are easy to compute:
$$
    J_{\pm} = \left| \begin{array}{cc} \frac{\partial a}{\partial x} & \frac{\partial a}{\partial y} \cr \frac{\partial b}{\partial x} & \frac{\partial b}{\partial y} \end{array} \right| = \left| \begin{array}{cc} \mp \frac{1}{\sqrt{y^2-4x}} & \pm \frac{a}{\sqrt{y^2-4x}} \cr \pm \frac{1}{\sqrt{y^2-4x}} & \mp\frac{b}{\sqrt{y^2-4x}}\end{array} \right| = \frac{b-a}{y^2- 4 x} = \mp \frac{1}{\sqrt{y^2-4x}}
$$
We can now read off $f_{X,Y}(x,y)$ from the product measure for $(A,B)$:
$$ \begin{eqnarray}
   f_A(a) f_B(b) \mathrm{d} a \mathrm{d} b &=& \left( \frac{1}{\Gamma(k)} \frac{a^{k-1}}{m^k} \mathrm{e}^{-\frac{a}{m}} \mathbf{1}_{a>0} \right) \cdot \left( \frac{1}{\Gamma(k)} \frac{b^{k-1}}{m^k} \mathrm{e}^{-\frac{b}{m}} \mathbf{1}_{b>0} \right) \mathrm{d} a \mathrm{d} b \\
   &=& \frac{(a b)^{k-1}}{\Gamma^2(k) m^{2k}} \mathrm{e}^{-\frac{a+b}{m}} \mathbf{1}_{a>0}  \mathbf{1}_{b>0} \mathrm{d} a \mathrm{d} b \\ &=&
   \frac{x^{k-1}}{\Gamma^2(k) m^{2k}} \mathrm{e}^{-\frac{y}{m}} I(x>0,y>0,y^2>4x) \left(| J_+ | + | J_- | \right) \mathrm{d} x \mathrm{d} y 
\end{eqnarray}
$$
