The question:
$$\begin{equation*} f(t)=e^t+e^t\int_0^te^{-\tau}f(\tau)\mathrm{d}\tau \end{equation*}.$$
I believe we need to take the Laplace transform of all terms. I am getting stuck with this part:
$$e^t\int_0^te^{-\tau}f(\tau)\mathrm{d}\tau.$$
I believe we could utilize the convolution theorem
$$\begin{equation*} (f\ast g)(t)=\int_0^tf(\tau)g(t-\tau)\mathrm{d}\tau \end{equation*}$$
but the $e^t$ in front of the integral is causing problems for me.
Any help is appreciated, thanks.
Hint
Suppose that $f$ is a solution of $$f(t)=e^t+e^t\int_0^te^{-\tau}f(\tau)\,\mathrm d \tau.$$
In one hand, $f(0)=1$. Moreover, if you derivate, you get \begin{align*} f'(t)&=e^t+\underbrace{e^t\int_0^t e^{-\tau}f(\tau)\,\mathrm d \tau}_{=f(t)-e^t}+f(t)=2f(t). \end{align*}
I let you conclude.