Consider $x:\mathbb{R}_+ \to \mathbb{R}_+$, where $x\not\in\mathcal{L}^1$, i.e., $$\int_0^\infty x(s)ds = \infty.$$ Consider the convolution $$h(t) = \int_0^t e^{-(t-s)}x(s)ds.$$ I am wondering if it is possible that $h\in\mathcal{L}^1$, i.e., $$\int_0^\infty h(s)ds = C < \infty.$$ I guess that yes, but I cannot find an example.
I know that for $x\in\mathcal{L}^1$ it holds $h \in \mathcal{L}^1$ but is the opposite true?
Since the functions under consideration are non-negative, you can do the following computation in $[0,+\infty]$ : \begin{align*} \int_0^\infty h(t) dt &= \int_0^\infty \left(\int_0^t e^{-(t-s)} x(s) d s \right) d t \\ &= \int_0^\infty \left(\int_s^\infty e^{-(t-s)} d t \right) x(s) d s \\ &= \int_0^\infty x(s) d s \end{align*}
So $x$ is integrable if, and only if, $h$ is.