I don't know how to obtain the autocovariance function of the following process, having a multiplication makes it difficult for me.
$X_t = Z_t + \theta Z_tZ_{t-1}$
with $Z_i \sim N(0, \sigma^2)$ (white noise)
pd: the answer is:
$\gamma_X(k) = \sigma^2(1+\theta\sigma^2), \hspace{1cm} k = 0$