I am dealing with a mean-reverting Vasicek process defined as:
\begin{equation} S_t = S_0 e^{-at} + b(1-e^{(-at)}) + \sigma e^{(-at)} \int_{0}^{t} e^{(-as)} \ W_t \end{equation}
I want to determine the following covariance:
\begin{equation} Cov[(S_{t+i}),(S_{t})] \end{equation}
Could someone help me with the analytical derivation? Thanks in advance!