Convolution of two functions from set $C$ We have space $L^1(\mathbb{R}^k) $ and a set $C =$ {$f \in L^1 (\mathbb{R}^k) : f(x) = e^{-a\cdot ||x||^2}, a>0 $}. Prove that if $f, g \in C$, then $f * g \in C$ (convolution of functions).
So the beginning of my solutions is: 
$$(f * g )(x) = \int\limits_{\mathbb{R}^k} f (x-y) g(y) dl_k(y) = \int\limits_{\mathbb{R}^k} e ^{-a||x-y||^2} e^ {-b ||y||^2} dl_k(y)
$$
$$=\int\limits_{\mathbb{R}^k} e^{-a ((y_1 -x_1)^2 + ... + (y_k - x_k)^2) - b(y_1^2 +... y_k^2)} dl_k(y)
$$ 
and...? I'm stuck here. Could somebody help?
 A: Edit: This is false. But almost true. Up to multiplication by a constant. 
Steps: Do the case $k=1$ first. Develop the quadratic in $y$. Complete the square. Change the variable appropriately to use $\int_\mathbb{R}e^{-t^2}dt=\sqrt{\pi}$. Then apply Fubini in the general $k$ case.
Details: For $k=1$, $f(x)=e^{-ax^2}$, and $g(x)=e^{-bx^2}$, we have
$$
(f*g)(x)=\int_{\mathbb{R}}e^{-a(x-y)^2}e^{-by^2}dy=\int_{\mathbb{R}}\exp{\left(-(a+b)\left(y^2-\frac{2ax}{a+b}y+\frac{ax^2}{a+b}\right)\right)}dy
$$
$$
=\int_{\mathbb{R}}\exp{\left(-(a+b)\left(\left(y-\frac{ax}{a+b}\right)^2+\frac{abx^2}{(a+b)^2}\right)\right)}dy
$$
$$
=e^{-\frac{abx^2}{a+b}}\int_{\mathbb{R}}\exp{\left(-(a+b)\left(y-\frac{ax}{a+b}\right)^2\right)}dy
$$
$$
=e^{-\frac{abx^2}{a+b}}\int_{\mathbb{R}}e^{-u^2}\frac{du}{\sqrt{a+b}}=\frac{\sqrt{\pi}}{\sqrt{a+b}}e^{-\frac{abx^2}{a+b}}.
$$
Back to the general $k$ case, and applying Fubini, we obtain
$$
(f*g)(x)=\prod_{j=1}^k\int_{\mathbb{R}}e^{-a(x_j-y_j)^2}e^{-by_j^2}dy=\prod_{j=1}^k\frac{\sqrt{\pi}}{\sqrt{a+b}}e^{-\frac{abx^2}{a+b}}
$$
$$
=\left(\frac{\pi}{a+b} \right)^\frac{k}{2}e^{-\frac{a}{a+b}\|x\|^2}.
$$
