\left(\begin{matrix} \epsilon_i \\ u_i \end{matrix}\right) &\sim& \left[\left(\begin{matrix} 0 \\ 0 \end{matrix}\right), \left(\begin{matrix} \sigma^2_u & \\ \rho \sigma_u \sigma_\epsilon & \sigma^2_\epsilon\\ \end{matrix}\right)\right]