Given that
$$f(x;\theta) = \frac{\pi x}{2 \theta^2} e^{-\frac{\pi x^2}{4 \theta^2}} \text{ for } x>0$$
I wish to find the method of moments estimators for $\theta$ and $\theta^2$.
I have calculated the $k$-th raw moment as $$E(X^k)=\Gamma\left(\frac{k}{2}+1\right) \frac{2^k \theta^k}{\pi^{\frac{k}{2}}}$$
Hence, substituting $k=1$, we obtain that $\hat \theta_{\text{mom}}=\bar X$. Now substituting $k=2$ $$\theta^2=\frac{\pi}{4}E(X^2)=\frac{\pi}{4}(Var(X)+E(X)^2)\Rightarrow\hat\theta^{2}_{\text{mom}}=\frac{\pi}{4}(S_X^2+\bar X^2 )$$
I believe this is how one should find the method of moments estimator for $\theta^2$ however I have come across a solution where they simply say that $$\hat \theta^2_{\text{mom}}=\bar X^2$$
Can anyone shed some light on this?