Suppose there are two correlated random variables $X_1$ and $X_2$ both are gamma distributed but having different shape and scale parameters with correlation coefficient $\rho$. What will be the distribution of Y? where $Y=(2X_1 X_2)/(X_1+X_2)$.