I've got a question about the distribution of the random variable $Y = X_{1}X_{2}$, where $X_i \sim \operatorname{N}(0,1)$, and $X_1$, $X_2$ independent.
Specifically, I need to find the MGF of $Y$.
Here's what I've tried:
The $r^\text{th}$ moment of $Y$ is given by:
$$\begin{align} \operatorname{E} (Y^r) & = \operatorname{E} (X_1^r X_2^r) \\ & = \operatorname{E} (X_1^r) \operatorname{E} (X_2^r) & \text{because } X_1 \text{ and } X_2 \text{ are independent} \\ \end{align}$$
Since $X_i \sim \operatorname{N}(0,1)$, I know that the $r^\text{th}$ central moment is given by:
$$ \operatorname{E}(X_i^r) = \begin{cases} 0 & r \text{ odd} \\ (2r-1)! & r \text{ even} \\ \end{cases}$$
Because $\operatorname{E}(X_i^r)$ does not depend on $i$, $\operatorname{E} (X_1^r) = \operatorname{E} (X_2^r)$ and so
$$\operatorname{E} (X_1^r) \operatorname{E} (X_2^r) = ((2r-1)!)^2 \quad \text{for } r \text{ even}$$
But I can't really derive a formula for the $r^{\text{th}}$ integral, especially if I don't know how the formula for $\operatorname{E}(X_i^r)$ looks before plugging in $t=0$.
Are there other ways to derive the MGF of $Y$?
EDIT:
I also need to show that $Y$ can be expressed as a difference of two independent, Gamma-distributed variables.