I assume you want X-Y and not X=Y.

Let Z=X+Y and W=X-Y, do a 2-2 transform, jacobian and all from calc3.

Then the new joint density should factor, giving you independence.

Note that sums of gamma's with same second beta is a gamma.

But I'm not sure how you're writing your exponential/gamma density.