Does a distribution on $\mathbb{R}$ exists such that for a random variable $X$ on $\mathbb{R}$ with said distribution we have $$ \mathbb{E}[X^{2l}] = l! \ \ \text{ and } \ \ \mathbb{E}[X^{2l+1}] = 0 $$ for all $l \in \mathbb{N}_0$?
The equality $\mathbb{E}[X^{2l+1}] = 0$ clearly holds for every symmetric distribution (where the moments exist) but I don't know of a distribution with such even moments.
Any help is much apprechiated.
From comments:
You could look for a $Y$ with $\mathbb E\left[Y^n\right]=n!$ and then let $X=\pm\sqrt{Y}$ with equal probability. A little bit of extra work would let you specify the distribution more simply.
Such a $Y$ is well known: it is the exponential distribution with parameter $1$.
This leads to $X$ having density $$f(x)=|x|\, e^{-x^2}$$ over the reals, as illustrated below. This is clearly symmetric about $0$, so the odd moments are $0$. A relatively simple integration shows the even moments are the desired factorials.