[step:Record how standard matrices behave under linear operations]
Let $e_1,\ldots,e_m$ denote the standard basis of $\mathbb{R}^m$, and let $\varepsilon_1,\ldots,\varepsilon_n$ denote the standard basis of $\mathbb{R}^n$. For $T\in\mathcal{L}(\mathbb{R}^m,\mathbb{R}^n)$, the entries of $A(T)$ are determined by
\begin{align*}
T(e_j)=\sum_{i=1}^{n} A(T)_{ij}\varepsilon_i
\end{align*}
for each $j\in\{1,\ldots,m\}$.
Let $S,T\in\mathcal{L}(\mathbb{R}^m,\mathbb{R}^n)$, and let $c\in\mathbb{R}$. For each $j\in\{1,\ldots,m\}$,
\begin{align*}
(S+T)(e_j)=S(e_j)+T(e_j)=\sum_{i=1}^{n}\bigl(A(S)_{ij}+A(T)_{ij}\bigr)\varepsilon_i.
\end{align*}
By uniqueness of coordinates in the basis $\varepsilon_1,\ldots,\varepsilon_n$,
\begin{align*}
A(S+T)_{ij}=A(S)_{ij}+A(T)_{ij}
\end{align*}
for every $i\in\{1,\ldots,n\}$ and $j\in\{1,\ldots,m\}$. Similarly,
\begin{align*}
(cT)(e_j)=cT(e_j)=\sum_{i=1}^{n} cA(T)_{ij}\varepsilon_i,
\end{align*}
so
\begin{align*}
A(cT)_{ij}=cA(T)_{ij}
\end{align*}
for every $i$ and $j$.
[/step]