Lista 1 - Exercício 3 - item b

set.seed(666)

noise <- rnorm(1000, 0, 1)

t = 1
Z <- noise[t]

t = 2
Z[t] <- .3*Z[t-1] + noise[t]

for (t in 3:length(noise)) {
  Z[t] <- .3*Z[t-1] - .585*Z[t-2] + noise[t]
}