\left[\begin{array}{c}
Y_1 \\
Y_2 \\
\vdots \\
Y_n
\end{array}\right]=\left[\begin{array}{ccccc}
1 & X_{11} & X_{12} & \cdots & X_{1 k} \\
\vdots & X_{21} & X_{22} & \cdots & X_{2 k} \\
\vdots & \vdots & \vdots & & \vdots \\
1 & X_{n 1} & X_{n 2} & \cdots & X_{n k}
\end{array}\right] \cdot\left(\begin{array}{c}
\beta_0 \\
\beta_1 \\
\beta_2 \\
\vdots \\
\beta_k
\end{array}\right)+\left[\begin{array}{c}
\mu_1 \\
\mu_2 \\
\vdots \\
\mu_n
\end{array}\right]