Multiplication of matrices — taking the dot product of the $i$th row of the first matrix and the $j$th column of the second to yield the $ij$th entry of the product — is not a very intuitive operation: if you were to ask someone how to mutliply two matrices, he probably would not think of that method. Of course, it turns out to be very useful: matrix multiplication is precisely the operation that represents composition of transformations. But it's not intuitive. So my question is where it came from. Who thought of multiplying matrices in that way, and why? (Was it perhaps multiplication of a matrix and a vector first? If so, who thought of multiplying them in that way, and why?) My question is intact no matter whether matrix multiplication was done this way only after it was used as representation of composition of transformations, or whether, on the contrary, matrix multiplication came first. (Again, I'm not asking about the utility of multiplying matrices as we do: this is clear to me. I'm asking a question about history.)

  • 1
    $\begingroup$ As to the last, matrix multiplication definitely came first (centuries first), and I'm reasonably certain from a compact representation of systems of linear equations. Leibniz already had a determinant formula. As I have no historic sources for first use, this doesn't answer your question though. $\endgroup$ Commented Jan 7, 2013 at 4:19
  • 1
    $\begingroup$ Matrices are linear operators and have meaning only when it acts on vectors. Given matrices $A$ and $B$, what would we want the operator/matrix $BA$ to mean? Ideally, we would want $BA$ to mean the following. For all vectors $x$, we want $(BA)x = B(Ax)$ Once we have this i.e. $(BA)x = B(Ax)$ for all $x$, then we are forced to live with the way we currently multiply matrices. And as to why matrix-vector product is defined in the way it is, the primary reason for introducing matrices was to handle linear transformation in a notationally convenient way. $\endgroup$
    – user17762
    Commented Jan 7, 2013 at 4:21
  • 1
    $\begingroup$ There is another question, of course, which is not so much why matrix multiplication was defined like this, but why it stuck - why this apparently curious definition took off, and didn't die the death of so many putative definitions. And that was because it proved mathematically fruitful. $\endgroup$ Commented Jan 7, 2013 at 8:47
  • 1
    $\begingroup$ Possible duplicate of math.stackexchange.com/questions/192835/…. $\endgroup$
    – lhf
    Commented Jun 2, 2014 at 12:45
  • $\begingroup$ Why is it "not intuitive"? If you ask someone how to multiply two matrices and they think about what that multiplication is supposed to mean, they absolutely will come up with the usual definition. $\endgroup$ Commented Nov 30, 2016 at 15:25

3 Answers 3


Matrix multiplication is a symbolic way of substituting one linear change of variables into another one. If $x' = ax + by$ and $y' = cx+dy$, and $x'' = a'x' + b'y'$ and $y'' = c'x' + d'y'$ then we can plug the first pair of formulas into the second to express $x''$ and $y''$ in terms of $x$ and $y$: $$ x'' = a'x' + b'y' = a'(ax + by) + b'(cx+dy) = (a'a + b'c)x + (a'b + b'd)y $$ and $$ y'' = c'x' + d'y' = c'(ax+by) + d'(cx+dy) = (c'a+d'c)x + (c'b+d'd)y. $$ It can be tedious to keep writing the variables, so we use arrays to track the coefficients, with the formulas for $x'$ and $x''$ on the first row and for $y'$ and $y''$ on the second row. The above two linear substitutions coincide with the matrix product $$ \left( \begin{array}{cc} a'&b'\\c'&d' \end{array} \right) \left( \begin{array}{cc} a&b\\c&d \end{array} \right) = \left( \begin{array}{cc} a'a+b'c&a'b+b'd\\c'a+d'c&c'b+d'd \end{array} \right). $$ So matrix multiplication is just a bookkeeping device for systems of linear substitutions plugged into one another (order matters). The formulas are not intuitive, but it's nothing other than the simple idea of combining two linear changes of variables in succession.

Matrix multiplication was first defined explicitly in print by Cayley in 1858, in order to reflect the effect of composition of linear transformations. See paragraph 3 at http://darkwing.uoregon.edu/~vitulli/441.sp04/LinAlgHistory.html. However, the idea of tracking what happens to coefficients when one linear change of variables is substituted into another (which we view as matrix multiplication) goes back further. For instance, the work of number theorists in the early 19th century on binary quadratic forms $ax^2 + bxy + cy^2$ was full of linear changes of variables plugged into each other (especially linear changes of variable that we would recognize as coming from ${\rm SL}_2({\mathbf Z})$). For more on the background, see the paper by Thomas Hawkins on matrix theory in the 1974 ICM. Google "ICM 1974 Thomas Hawkins" and you'll find his paper among the top 3 hits.

  • 2
    $\begingroup$ Brilliant way to explain this! I've been searching for a good reasoning for why matrix multiplication is defined the way it is and everyone seems miss the point of the question. This answer, thankfully, explains it in the correct way (for me to understand). $\endgroup$
    – Sia
    Commented Dec 29, 2022 at 20:51
  • 1
    $\begingroup$ Interestingly, Wikipedia mentions that matrix multiplication was first described by Binet in 1812 to represent the composition of linear maps that are represented by matrices. en.wikipedia.org/wiki/Matrix_multiplication#cite_note-2 $\endgroup$ Commented May 1, 2023 at 20:18

Here is an answer directly reflecting the historical perspective from the paper Memoir on the theory of matrices By Authur Cayley, 1857. This paper is available here.

This paper is credited with "containing the first abstract definition of a matrix" and "a matrix algebra defining addition, multiplication, scalar multiplication and inverses" (source).

In this paper a nonstandard notation is used. I will do my best to place it in a more "modern" (but still nonstandard) notation. The bulk of the contents of this post will come from pages 20-21.

To introduce notation, $$ (X,Y,Z)= \left( \begin{array}{ccc} a & b & c \\ a' & b' & c' \\ a'' & b'' & c'' \end{array} \right)(x,y,z)$$

will represent the set of linear functions $(ax + by + cz, a'x + b'y + c'z, a''x + b''y + c''z)$ which are then called $(X,Y,Z)$.

Cayley defines addition and scalar multiplication and then moves to matrix multiplication or "composition". He specifically wants to deal with the issue of:

$$(X,Y,Z)= \left( \begin{array}{ccc} a & b & c \\ a' & b' & c' \\ a'' & b'' & c'' \end{array} \right)(x, y, z) \quad \text{where} \quad (x, y, z)= \left( \begin{array}{ccc} \alpha & \beta & \gamma \\ \alpha' & \beta' & \gamma' \\ \alpha'' & \beta'' & \gamma'' \\ \end{array} \right)(\xi,\eta,\zeta)$$

He now wants to represent $(X,Y,Z)$ in terms of $(\xi,\eta,\zeta)$. He does this by creating another matrix that satisfies the equation:

$$(X,Y,Z)= \left( \begin{array}{ccc} A & B & C \\ A' & B' & C' \\ A'' & B'' & C'' \\ \end{array} \right)(\xi,\eta,\zeta)$$

He continues to write that the value we obtain is:

$$\begin{align}\left( \begin{array}{ccc} A & B & C \\ A' & B' & C' \\ A'' & B'' & C'' \\ \end{array} \right) &= \left( \begin{array}{ccc} a & b & c \\ a' & b' & c' \\ a'' & b'' & c'' \end{array} \right)\left( \begin{array}{ccc} \alpha & \beta & \gamma \\ \alpha' & \beta' & \gamma' \\ \alpha'' & \beta'' & \gamma'' \\ \end{array} \right)\\[.25cm] &= \left( \begin{array}{ccc} a\alpha+b\alpha' + c\alpha'' & a\beta+b\beta' + c\beta'' & a\gamma+b\gamma' + c\gamma'' \\ a'\alpha+b'\alpha' + c'\alpha'' & a'\beta+b'\beta' + c'\beta'' & a'\gamma+b'\gamma' + c'\gamma'' \\ a''\alpha+b''\alpha' + c''\alpha'' & a''\beta+b''\beta' + c''\beta'' & a''\gamma+b''\gamma' + c''\gamma''\end{array} \right)\end{align}$$

This is the standard definition of matrix multiplication. I must believe that matrix multiplication was defined to deal with this specific problem. The paper continues to mention several properties of matrix multiplication such as non-commutativity, composition with unity and zero and exponentiation.

Here is the written rule of composition:

Any line of the compound matrix is obtained by combining the corresponding line of the first component matrix successively with the several columns of the second matrix (p. 21)

  • 2
    $\begingroup$ Should the set of linear functions $(ax + by + cz, a'z + b'y + c'z, a''z + b''y + c''z)$ be $(ax + by + cz, a'x + b'y + c'z, a''x + b''y + c''z)$? $\endgroup$ Commented May 11, 2017 at 18:11
  • $\begingroup$ Brad's is THE answer, I think. The eventual coefficients obtained from a double transformation of (x, y, z) require row elements of the left transform matrices to be multiplied by their corresponding column element of the right transform matrix. My own old idea was that if multiplication of a vector by a matrix is done like A (x, y, z) = (x', y', z'), i.e. using row vectors, then the resulting system will not "read" as clearly as if the (x, y, z) & (x', y', z') were written as a column vectors. And this is only for a 3-D system. Higher order systems harder still. But Brad's idea is dead-on. $\endgroup$
    – Trunk
    Commented Dec 28, 2017 at 15:57

\begin{align} u & = 3x + 7y \\ v & = -2x + 11y \\ \\ \\ \\ p & =13u-20v \\ q & = 2u+6v \end{align} Given $x$ and $y$, how do you find $p$ and $q$? How do you write: \begin{align} p & = \bullet\, x + \bullet\, y \\ q & = \bullet\, x+\bullet\, y\quad\text{?} \end{align} What numbers go where the four $\bullet$'s are?

That is what matrix multiplication is. The rationale is mathematical, not historical.


You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .