I am trying to figure out how to prove that this function is Lipschitz. I tried the following (for $h>0$):
$$ \left|e^{-(x+h)^2}-e^{-x^2}\right|=\left|\frac{e^{-(x+h)^2}-e^{-x^2}}{h}\right|\cdot|h| $$
By playing with Desmos, I noticed that $\displaystyle\left|\frac{e^{-(x+h)^2}-e^{-x^2}}{h}\right|$ increases as $h$ decreases. Thus, we get an upper bound if we let $h$ approach $0.$ So we get $\left|-2xe^{-x^2}\right|$ and, by elementary calculus, we get a uniform maximum of $\sqrt{\frac{2}{e}},$ which we may choose as our Lipschitz constant. What I wonder is if there is a better (more analytic) argument.