I bumped into the following problem and can't see how I am supposed to calculate divide the intergral into correct areas of definition.
The problem is as follows: $X$ and $Y$ are random variables with joint density as
$$f_{X,Y} = \begin{cases}(2/5)(2x+3y) & 0 < x, y < 1 \\ 0 & \text{o.w.} \end{cases},$$
I want to calculate the distribution of $2X+3Y$.
So I did the following:
Set $U = X$ and $V = 2X + 3Y$ from the latter I derived $Y = (V - 2U)/3$. Here I see that if $X$ and $Y$ are defined in $[0,1]$ then $V$ is defined in $[0,5]$.
Now: $$\begin{align*} f_v(v) &= \int f_{U,V}(u,v) du \\ &= \int f_{X,Y}\left(u, \frac{v - 2u}{3}\right) |J| du \\ &= \frac{2}{15} \int 2u + 3\cdot\frac{v - 2u}{3} du \\ &=\frac{2 v}{15} \int du \\ &= \frac{2 v}{15} \left[ \, u \, \big|_0^v \, \right] \\ &= \frac{2v^2}{ 15} \end{align*}$$
But the answer is something else.

The method you want to use for the sum of random variables is called convolution. It goes like this:
You want
$$f_{2X + 3Y}(z) = P(2X + 3Y = z) = P\left(X = x, Y = \dfrac{z - 2x}{3}\right) = f_{X,Y}\left(x, \dfrac{z - 2x}{3}\right)$$
for a given $z$. You are right that $0 < 2X + 3Y < 5$. So, using the joint distribution of $X$ and $Y$ to get
$$ \begin{align*} f_{2X + 3Y}(z) &= \int_{0}^{5} f_{X,Y}\left(x, \dfrac{z - 2x}{3}\right) \, dx \\ \\ &= \frac{2}{5} \int_{0}^{5} 2x + 3\left(\frac{z-2x}{3}\right) \, dx \\ \\ &= \frac{2}{5} \int_{0}^{5} 3z \, dx \\ \\ &= 6z. \end{align*}$$