[motivation]
### The search for well-behaved transformations
Consider a function $f \colon \mathbb{R}^2 \to \mathbb{R}^2$ that transforms the plane. Among all possible such [functions](/page/Function) --- reflections, translations, dilations, polynomials, wild discontinuous maps --- which ones can we hope to understand completely? The answer begins with a structural observation: the functions that respect the algebraic structure of $\mathbb{R}^2$ are precisely the ones amenable to systematic analysis.
What does "respect the structure" mean? A vector space has two operations: we can add vectors and scale them. A function $\alpha \colon \mathbb{R}^2 \to \mathbb{R}^2$ that preserves both operations satisfies
\begin{align*}
\alpha(u + v) &= \alpha(u) + \alpha(v), \\
\alpha(\lambda u) &= \lambda \, \alpha(u),
\end{align*}
for all vectors $u, v \in \mathbb{R}^2$ and all scalars $\lambda \in \mathbb{R}$. These two conditions are the definition of a linear map.
### What fails without linearity
To appreciate why these conditions are so powerful, consider what goes wrong without them. Take the squaring function $g \colon \mathbb{R} \to \mathbb{R}$ defined by $g(x) = x^2$. This map is decidedly nonlinear: $g(2 + 3) = 25$, while $g(2) + g(3) = 13$. As a consequence, $g$ has no well-defined kernel in the algebraic sense. The preimage $g^{-1}(\{0\}) = \{0\}$ is a single point, but there is no analogue of the rank-nullity theorem: knowing $\dim(\ker g)$ tells us nothing about the dimension of the image.
More critically, a nonlinear function cannot be represented by a matrix. The map $h \colon \mathbb{R}^2 \to \mathbb{R}^2$ given by $h(x, y) = (x^2, xy)$ sends lines through the origin to curves, not lines. There is no hope of finding a $2 \times 2$ matrix $A$ such that $h(v) = Av$ for all $v$. Without matrix representation, we lose the entire computational apparatus of determinants, eigenvalues, and canonical forms.
### The matrix connection
For a linear map $\alpha \colon U \to V$ between finite-dimensional vector spaces, choosing a basis for $U$ and a basis for $V$ produces a matrix $A$ such that $\alpha(u)$ is computed by matrix multiplication. Different choices of bases yield different matrices, but all such matrices are related by a change-of-basis transformation. This means that intrinsic properties of $\alpha$ --- its rank, nullity, determinant (when $U = V$), and eigenvalues --- are invariants that do not depend on the coordinate system. The single definition of linearity thus unifies three perspectives: the geometric (transformations of space), the computational (matrix arithmetic), and the abstract (structure-preserving homomorphisms between algebraic objects).
[/motivation]