If $X$ and $Y$ have joint density function
$$
f(x,y)=\begin{cases}
\frac{1}{2}(x+y)e^{-x-y}&\text{if }x,y>0,\\
0&\text{otherwise,}
\end{cases}
$$
find the density function of $X$ and $Y$. I was considering two things:
1) I could check whether $X$ and $Y$ are independent, and if so, use the convolution formula $f_Z(z)=\int_{-\infty}^\infty f_X(x)f_Y(z-x)dx$,
2) I could differentiate $\mathbb P(X+Z\leq z)=\int_{-\infty}^\infty\int_{-\infty}^{z-x}\frac{1}{2}(x+y)e^{-x-y}dydx$.
I'm not sure if 2) is reasonably doable. Are these my two only options? Or is there a better way?