Estoy tratando de calcular la covarianza de un proceso de Ornstein-Uhlenbeck unidimensional 1$dx_t=-\theta x_t dt+ \sigma dW_t$,$\theta>0$ y estoy en la etapa,$$\text{Cov }(x_s,x_t)=\sigma^2 e^{-\theta(t+s)} \mathbb{E}\left[ \int_0^s e^{\theta u}dW_u \int_0^t e^{\theta v} dW_v\right].$ $
¿Es posible evaluar las integrales estocásticas explícitamente y, si no, cómo se hace para simplificar esto? En Wikipedia, dicen que esto es igual a$\frac{\sigma^2}{2\theta} e^{-\theta(t+s)} (e^{2\theta s\wedge t } -1)$ pero no puedo ver cómo llegan a esta conclusión.
Gracias.