how to find the distribution density of the variable X + Y? X and Y are independent and have the same exponential distribution wih density: $$\lambda e^{-\lambda x} 1_{(0,\infty)} (x) $$ The density of the X + Y distribution can be calculated by convolution of the functions: $$f*f(t)=\int_{-\infty}^{\infty} f(x) f(t-x) \, dx $$ but I still have no idea how to solve it, during the lecture we only gave an example for a uniform distribution on (0,1)

\begin{align}
& \int_{-\infty}^\infty \lambda e^{-\lambda x} 1_{(0,\infty)} (x) \cdot \lambda e^{-\lambda (t-x)} 1_{(0,\infty)} (t-x) \, dx \\[12pt]
= {} & \int_0^\infty \lambda e^{-\lambda x}(x) \cdot \lambda e^{-\lambda(t-x)} 1_{(0,\infty)}(t-x)\,dx \\
& \text{because } 1_{(0,\infty)}(x)=0 \text{ when } x<0 \\[12pt]
= {} & \int_0^t \lambda e^{-\lambda x}(x) \cdot \lambda e^{-\lambda(t-x)} \,dx \\
& \text{because } 1_{(0,\infty)}(t-x)=0 \text{ when } x>t \\[12pt]
= {} & \int_0^t \lambda^2 e^{-\lambda t} \,dx
\end{align}
and then observe that $\lambda^2 e^{-\lambda t}$ does not change as $x$ goes from $0$ to $t,$ i.e. you're integrating a constant function.
