Let $Z\sim N(0,1)$ and $Y=a+bZ+cZ^2$. I want to compute the variance of $Y$. This is what I did: $$\operatorname{Var}(Y)=0+b^2\operatorname{Var}(Z)+c^2\operatorname{Var}(Z^2)=b^2+c^2\operatorname{Var}(Z^2)$$ To get $\operatorname{Var}(Z^2)$, I tried to use the definition $\operatorname{Var}(Z^2)=\mathbb{E}[Z^4]-\mathbb{E}[Z^2]^2$ But im having with this part. If this was a odd for example $\mathbb{E}[Z^3]$ you can say that because of the symmetry of the normal distribution $\mathbb{E}[Z^3]=0$, but in this is pair.
Thank you
$$\mathbb{E}[Z^4]=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}z^4e^{\large-\frac{z^2}{2}}dz=\frac{2}{\sqrt{2\pi}}\int_{0}^{\infty}z^4e^{\large-\frac{z^2}{2}}dz$$ set $\frac{z^2}{2}=u$, we have $z\,dz=du$ and $z^3=2\sqrt{2}u\sqrt{u}$, thus $$\mathbb{E}[Z^4]=\frac{4}{\sqrt{\pi}}\int_{0}^{\infty}u\sqrt{u}\,e^{-u}du=\frac{4}{\sqrt{\pi}}\Gamma\left(\frac{5}{2}\right)=\frac{4}{\sqrt{\pi}}\times \frac32 \times \frac 12\Gamma\left(\frac{1}{2}\right)=3$$ please check this link
Note $$\Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}$$ Generally $$\mathbb{E}[Z^{2k}]=\frac{(2k)!}{2^k k!}$$