Equations for VAR model with GARCH errors How do we actually write the VAR-GARCH model? I don't really understand it.
Here is the VAR($m$) process:
$$ x_t = \phi_0 + \sum_{i=1}^m \phi_i x_{t-i} + \epsilon_t; $$
and here is the GARCH(p,q):
\begin{eqnarray*} 
Var[\boldsymbol{x_t}|\boldsymbol{I}_{t-1}] &=& Var[\boldsymbol{\epsilon_t}|I_t] \\
\sigma_t^2 &=& \alpha_0 + \sum_{l=1}^{p} \alpha_l \sigma^2_{t-l} + \sum_{m=1}^{q} \beta_m x^2_{t-m} \\
\end{eqnarray*}
As I understand, the variance is modeled as GARCH and we use VAR to model the conditional mean of the variables. How do we actually write the complete equation of VAR($m$)-GARCH($p,q$)?
 A: You are pretty close to your goal already, but you should replace $x$s with $\epsilon$s in the conditional variance equation. Note also that VAR is a multivariate model, i.e. the dependent variable is a vector rather than a scalar. Then the GARCH model needs to be multivariate as well because it needs to characterize the conditional covariance matrix of a vector of errors. 
So for VAR($m$)-GARCH($p,q$) you have
$$
\begin{aligned}
x_t &= \phi_0 + \sum_{i=1}^m \phi_i x_{t-i} + \epsilon_t \\
\epsilon_t &= \Sigma^{1/2}_t u_t \\
\Sigma_t &= \alpha_0 + \sum_{l=1}^{p} \alpha_l \Sigma_{t-l} + \sum_{m=1}^{q} \beta_m \epsilon_{t-m}\epsilon'_{t-m} \\
u_t &\sim i.i.d(0,I)
\end{aligned}
$$
where $x_t$ and $\epsilon_t$ are vectors; $\phi_0$ and $\alpha_0$ are parameter vectors; $\Sigma_t$ is the conditional variance matrix of $\epsilon_t$; $I$ is an identity matrix; and $\phi_1$ through $\phi_m$, $\alpha_1$ through $\alpha_p$, and $\beta_1$ through $\beta_q$ are parameter matrices. Some restrictions in the equation of $\Sigma_t$ are needed to ensure $\Sigma_t$ is a valid covariance matrix (positive definite). Also, the equation for $\Sigma_t$ could be replaced by an equation for vech-torized $\Sigma_t$ (stacking the columns below and including the diagonal) to avoid having a matrix as a dependent variable; that would change the representation but not the basic idea.
I was trying to stick to your notation as closely as possible; you would encounter somewhat different notation in time series textbooks (e.g. using capital letters for matrices, vech-torizing $\Sigma_t$, etc.).
