Determinants are defined here.
Suppose $A$ is a matrix. The $(i, j)$-entry of $A$ is denoted as $[A]_{i,j}$.
Suppose that $a_1$, $a_2$, $\dots$, $a_n$ are $m \times 1$ matrices. Then $[a_1, a_2, \dots, a_n]$ is the $m \times n$ matrix whose $(i, j)$-entry is equal to $[a_j]_{i,1}$. The piece of notation allows one to display a matrix using its columns.
How does one prove the multilinear property?
Suppose that $j$ is a positive integer less than or equal to $n$. Suppose that $a_1$, $\dots$, $a_{j-1}$, $a_{j+1}$, $\dots$, $a_n$ are $n \times 1$ matrices. Suppose that $x$, $y$ are $n \times 1$ matrices. Suppose that $s$, $t$ are numbers. Then $$ \begin{aligned} & \det {[a_1, \dots, a_{j-1}, \overset{\text{column}\,j}{sx + ty}, a_{j+1}, \dots, a_n]} \\ = {} & s \det {[a_1, \dots, a_{j-1}, x, a_{j+1}, \dots, a_n]} + t \det {[a_1, \dots, a_{j-1}, y, a_{j+1}, \dots, a_n]}. \end{aligned} $$
The proof is indirect, which is based on no more than the definition of the determinant of a square matrix.
Let the proposition $P(n)$ be as follows:
It will be proved that $P(n)$ is true for $n = 1$, $2$, $3$, $\dots$ by mathematical induction.
$P(1)$ is true, which is left as an exercise.
Suppose that $P(n-1)$ is true. It will be proved that $P(n)$ is true under the hypothesis.
Define $$ \begin{aligned} & A = [a_1, \dots, a_{j-1}, x, a_{j+1}, \dots, a_n], \\ & B = [a_1, \dots, a_{j-1}, y, a_{j+1}, \dots, a_n], \\ & C = [a_1, \dots, a_{j-1}, sx + ty, a_{j+1}, \dots, a_n]. \end{aligned} $$ If $k \neq j$, then $[A]_{i,k} = [B]_{i,k} = [C]_{i,k}$. Hence $A(i|j) = B(i|j) = C(i|j)$. Note that $[C]_{i,j} = s[A]_{i,j} + t[B]_{i,j}$.
Suppose first that $j = 1$. Then $$ \begin{aligned} & \det {(C)} \\ = {} & \sum_{i = 1}^{n} { (-1)^{i+1} [C]_{i,1} \det {(C(i|1))} } \\ = {} & \sum_{i = 1}^{n} { (-1)^{i+1} (s[A]_{i,1} + t[B]_{i,1}) \det {(C(i|1))} } \\ = {} & s\sum_{i = 1}^{n} { (-1)^{i+1} [A]_{i,1} \det {(C(i|1))} } + t\sum_{i = 1}^{n} { (-1)^{i+1} [B]_{i,1} \det {(C(i|1))} } \\ = {} & s\sum_{i = 1}^{n} { (-1)^{i+1} [A]_{i,1} \det {(A(i|1))} } + t\sum_{i = 1}^{n} { (-1)^{i+1} [B]_{i,1} \det {(B(i|1))} } \\ = {} & s \det{(A)} + t \det{(B)}. \end{aligned} $$
Suppose then that $j > 1$. Note that for $\ell \neq j - 1$, column $\ell$ of $A(i|1)$, that of $B(i|1)$ and that of $C(i|1)$ are equal. Suppose that column $j-1$ of $A(i|1)$, that of $B(i|1)$ and that of $C(i|1)$ are $p$, $q$ and $r$, respectively. Then $r = sp + tq$. By the hypothesis, one has $$ \det {(C(i|1))} = s\det {(A(i|1))} + t\det {(B(i|1))}. $$ Hence $$ \begin{align*} & \det {(C)} \\ = {} & \sum_{i = 1}^{n} {(-1)^{i+1} [C]_{i,1} \det {(C(i|1))}} \\ = {} & \sum_{i = 1}^{n} {(-1)^{i+1} [C]_{i,1} (s\det {(A(i|1))} + t\det {(B(i|1))})} \\ = {} & s \sum_{i = 1}^{n} {(-1)^{i+1} [C]_{i,1} \det {(A(i|1))}} + t \sum_{i = 1}^{n} {(-1)^{i+1} [C]_{i,1} \det {(B(i|1))}} \\ = {} & s \sum_{i = 1}^{n} {(-1)^{i+1} [A]_{i,1} \det {(A(i|1))}} + t \sum_{i = 1}^{n} {(-1)^{i+1} [B]_{i,1} \det {(B(i|1))}} \\ = {} & s \det {(A)} + t \det {(B)}. \end{align*} $$
Hence by mathematical induction, $P(n)$ is true for $n = 1$, $2$, $3$, $\dots$.