OU-process in VAR form

For equidistant observation times $t_i=i\,\delta$ the OU-process takes the VAR(1) form

$$\xi_{i+1} -\mu_\infty = e^{\delta\,A}(\xi_i-\mu_\infty) + V_\delta\,\epsilon_i$$ (73)

where the $\epsilon_i$ are iid $\mathcal N(0,1)$.