I know this may be a stupid question, but still hope someone can help me. I was trying to prove the Fourier inversion formula for which $f$ and $\hat{f}=\int_{\mathbb{R}}f(x)e^{-i2\pi xy}dx$ both lie in $L^1$, and then $f=\int_{\mathbb{R}}\hat{f}(y)e^{i2\pi xy}dy$. I assume the knowing of FT of the Gaussian function $g_c(x)=e^{-c\pi x^2}$, which is $\frac{1}{\sqrt{c}}e^{-\frac{1}{c}\pi x^2}$. Then I have $$ \int_{\mathbb{R}}\hat{f}(y)e^{i2\pi xy}dy=\lim_{c\to0}\int_{\mathbb{R}}\hat{f}(y)g_c(y)e^{i2\pi xy}dy\\ $$
Using the definition of $\hat{f}$ and Fubini, this becomes $$ \lim_{c\to0}\int_{\mathbb{R}}\int_{\mathbb{R}}f(t)e^{-i2\pi ty}dtg_c(y)e^{i2\pi xy}dy= \lim_{c\to0}\int_{\mathbb{R}}f(t)\hat{g}_c(t-x)dt $$
Now I came to the problem. Formally one can guess that $\hat{g}_c(t)$ is just the Dirac function as $c\to0$ so that convolution gives $f$ itself, but how can one show that strictly, especially for an $L^1$ function? On the other hand, I also directly calculated $g_c(t-x)$ which gives me $$ \lim_{c\to0}\int_{\mathbb{R}}f(t)\frac{1}{\sqrt{c}}e^{-\frac{1}{c}\pi(t-x)^2}dt. $$Change of variables gives me $$ \lim_{c\to0}\int_{\mathbb{R}}f(x-\sqrt{c}y)e^{-\pi y^2}dy $$
Then why is it okay for me to pass the limit directly to $f$, since $f$ is just absolutely integrable instead of being continuous? I know this seems true formally, but, again, how can one show this strictly? Thank you for help!