Compute $$\iint_D e^{(x+y)^2} \ dxdy,$$
where $D=\{ y\leq3x, \ x\leq 3y, \ 0 \leq x+y \leq 2 \}.$
The area $D$ is easily drawn:
From this image, it's easy to seet that we can split the integral into two parts:
$$\iint_D e^{(x+y)^2} \ dxdy = \int_0^{1/2}\int_{x/3}^{3x}e^{(x+y)^2} \ dydx+\int_{1/2}^{3/2}\int_{x/3}^{2-x}e^{(x+y)^2} \ dydx = \frac{e^4-1}{4}.$$
This is the correct answer. The problem here is that I actually did not compute those two integrals myself, I entered them in Maple and for each one of them I got the answers in terms of the error function, however they canceled out luckily and I was left with the correct answer.
Can anyone help me evaluating this integral by a clever substitution so that I don't have to deal with non-elementary functions?

HINT...You can use the following change of variables: let $$u=x+y$$ and $$v=\frac yx$$
The Jacobian is $$\frac{u}{(1+v)^2}$$ and the integral works out quite easily to get your answer. I shall leave this to do yourself.
I hope this helps.