This might be a very simple question, but I couldn't find an answer yet.
Consider a multidimensional stochastic SDE of the form
$S_t = x + b(S_t) dt + \sigma(S_t) dW_t $
where $b: \mathbb{R}^n \to \mathbb{R}^n $, $\sigma: \mathbb{R}^{n \times n} \to \mathbb{R}^{n \times n}$ are measurable functions, $x \in \mathbb{R}^n$, $W_t $ is a vector of independent Wiener processes and it is assumed that a strong solution for $S_t$ exists. Is it then possible to calculate the covariance matrix of $S_{\tilde{t}}$ for some time $\tilde{t}$ by
$Cov(S_{\tilde{t}}) = \int_{0}^{\tilde{t}} \frac{d Cov(S_t)}{dt}dt = \mathbb{E} [\int_{0}^{\tilde{t}} \sigma(S_t) \sigma(S_t)^T dt] \quad .$
Heuristically the change of the covariance matrix should depend on
$\sigma(S_t) dW_t \sim N(0, \sigma(S_t) \sigma(S_t)^T dt) \quad .$
So my first question would be wether the formula for $Cov(S_{\tilde{t}})$ is correct and if that's the case, wether someone knows a more formal argument for that.