Invertible Matrices

Can a matrix have a left inverse but no right inverse? This question sits at the heart of understanding square matrices. In infinite-dimensional settings the answer is yes. A shift operator can have a one-sided inverse without having a two-sided one. But in finite dimensions something remarkable happens. For square matrices, one-sidedness collapses. A single equation BA=IBA = I forces AB=IAB = I as well. This asymmetry between the finite and infinite worlds is not a technicality. It is a load-bearing fact used constantly, often without comment, throughout linear algebra. Invertible matrices matter because they are exactly the matrices representing bijective linear maps. These are the algebraic shadows of coordinate changes, symmetries, and reversible processes. Solving Ax=bAx = b uniquely for every bb depends on this. So does changing basis, defining determinants and eigenvalues, and building the machinery of Rank of a Matrix. All of these lean on a precise understanding of what invertibility means and how to detect it. The set of invertible n×nn\times n matrices also carries its own rich structure. It forms a group, the general linear group GLn(F)GL_n(F), under matrix multiplication. This group is the algebraic home for symmetry itself. The big picture is this: invertibility of a square matrix AA is a single yes/no condition. It is equivalent to a long list of others. These include nonzero determinant, full rank, trivial kernel, columns forming a basis, and more. This article isolates the cleanest version of this equivalence. It proves the deep one-sided fact rigorously, using only Matrix Multiplication as Composition and known facts about linear maps on finite-dimensional spaces. It also shows exactly where the finite-dimensionality hypothesis is used, and what breaks without it. By the end you will have a complete proof of several key facts. For square matrices over a field, injectivity, surjectivity, and bijectivity of the associated linear map are all equivalent. A one-sided inverse is automatically two-sided. The invertible matrices form a group. We also preview how row reduction (Elementary Matrices and Row Operations) gives a computational test for invertibility.

Key takeaways
  • A square matrix AA is invertible if there exists BB with AB=BA=IAB = BA = I; such BB is unique and written A−1A^{-1}.
  • In finite dimensions, a one-sided inverse of a square matrix is automatically a two-sided inverse — this fails in infinite dimensions.
  • Invertibility of A∈Fn×nA \in F^{n\times n} is equivalent to: det⁡A≠0\det A \ne 0, rank⁡(A)=n\operatorname{rank}(A) = n, columns (or rows) form a basis, and the map x↦Axx \mapsto Ax is injective (equivalently surjective).
  • The invertible n×nn \times n matrices over a field FF form a group GLn(F)GL_n(F) under multiplication, generally non-abelian for n≥2n \ge 2.
  • (AB)−1=B−1A−1(AB)^{-1} = B^{-1}A^{-1} whenever AA and BB are invertible — order reverses, exactly as for composing invertible functions.
  • Key idea: invertibility of a square matrix is a rigid, all-or-nothing property fully witnessed by either side of the inverse equation, thanks to finite dimensionality.

Where this lives today

Invertible matrices are the algebraic backbone of change of basis: passing between coordinate systems in a vector space requires an invertible transition matrix, used throughout diagonalization and canonical form theory. In differential equations and dynamical systems, the invertibility of the flow map (or of a linear system's coefficient matrix) governs whether trajectories can be uniquely reversed in time. In cryptography, invertible matrices over finite fields (e.g., GLn(F2)GL_n(\mathbb{F}_2)) underlie linear codes and classical cipher constructions such as the Hill cipher, where decryption literally requires computing A−1A^{-1}. In computer graphics and robotics, invertible matrices represent reversible geometric transformations (rotations, scalings, shears), and GLn(R)GL_n(\mathbb{R}) together with its subgroups (like O(n)O(n), SLn(R)SL_n(\mathbb{R})) parametrizes the symmetries used in rendering pipelines and kinematics. In statistics, invertibility of the design matrix's Gram matrix XTXX^TX is precisely the condition for a unique least-squares solution in linear regression. In quantum mechanics, invertible (indeed unitary) matrices represent reversible time evolution of finite-dimensional quantum states.

Ask the tutor: Why exactly does the finite-dimensionality of FnF^n break down for the shift operator example, and where precisely in the Rank–Nullity proof does that assumption get used?

DefinitionInvertible matrix

Let FF be a field and let A∈Fn×nA \in F^{n\times n} be a square matrix. Say AA is invertible (or nonsingular) if there exists a matrix B∈Fn×nB \in F^{n\times n} such that

AB=InandBA=In,AB = I_n \quad \text{and} \quad BA = I_n,

where InI_n is the n×nn\times n identity matrix. If such BB exists it is called the inverse of AA and denoted A−1A^{-1}. A matrix that is not invertible is called singular. The set of all invertible n×nn\times n matrices over FF is denoted

GLn(F)={A∈Fn×n:A is invertible},GL_n(F) = \{A \in F^{n\times n} : A \text{ is invertible}\},

called the general linear group of degree nn over FF.

Intuition

Think of a square matrix AA as encoding a linear map TA:Fn→FnT_A : F^n \to F^n, TA(x)=AxT_A(x) = Ax, in a fixed basis (see Matrix Multiplication as Composition). Composition of linear maps corresponds to matrix multiplication, and the identity map corresponds to II. So "AA is invertible" means exactly "TAT_A is an invertible function," i.e., TAT_A has a two-sided inverse function. The content of this article is that for square matrices on a finite-dimensional space, checking only one direction of invertibility — say that TAT_A is injective — is already enough to conclude the whole story: TAT_A is bijective, and the inverse map is also linear, hence represented by an honest inverse matrix A−1A^{-1}. Picture FnF^n as a rigid frame of nn independent directions (a basis). Applying AA moves this frame to the frame given by the columns of AA. If those images stay independent (never collapsing into a lower-dimensional subspace), the map cannot lose information, so it must also be onto: there's no "room" left over in an nn-dimensional target for the image of an nn-dimensional domain to be a proper subspace and still capture everything only if it fails to be injective. This finite counting argument — dimension in equals dimension out — is the engine behind the one-sided-inverse theorem.

Figure
Columns of AA as images of the standard basis; invertibility means these images stay independent and span the plane
e₁ e₂ domain Fⁿ A = [[2,1],[1,1]] Ae₁=(2,1) Ae₂=(1,1) image spans Fⁿ: invertible
TheoremOne-sided inverses suffice for square matrices

Let FF be a field, let n≥1n \ge 1 be an integer, and let A∈Fn×nA \in F^{n\times n}. Suppose there exists B∈Fn×nB \in F^{n\times n} such that

BA=In.BA = I_n.

Then AA is invertible, and in fact AB=InAB = I_n as well, so B=A−1B = A^{-1}. Equivalently (dually): if there exists C∈Fn×nC \in F^{n\times n} with AC=InAC = I_n, then also CA=InCA = I_n and C=A−1C = A^{-1}. Consequently, for a square matrix AA over FF, the following are equivalent: 1. AA is invertible (has a two-sided inverse). 2. There exists BB with BA=InBA = I_n (a left inverse). 3. There exists CC with AC=InAC = I_n (a right inverse). 4. The linear map TA:Fn→FnT_A : F^n \to F^n, TA(x)=AxT_A(x) = Ax, is injective. 5. The linear map TAT_A is surjective. 6. The columns of AA form a basis of FnF^n.

LemmaRank–nullity restated for square matrices

Let VV be a finite-dimensional vector space with dim⁡V=n\dim V = n, and let T:V→VT: V \to V be linear. Then

dim⁡(ker⁡T)+dim⁡(range⁡T)=n.\dim(\ker T) + \dim(\operatorname{range} T) = n.

Consequently TT is injective if and only if TT is surjective.

Proof

This is the Rank–Nullity Theorem applied with domain and codomain both equal to VV; we use it as a known prerequisite result about linear maps on finite-dimensional spaces (see Invertibility of Linear Maps). Injective ⇒\Rightarrow surjective. If TT is injective, ker⁡T={0}\ker T = \{0\}, so dim⁡(ker⁡T)=0\dim(\ker T) = 0. By Rank–Nullity, dim⁡(range⁡T)=n=dim⁡V\dim(\operatorname{range} T) = n = \dim V. Since range⁡T\operatorname{range} T is a subspace of VV with the same finite dimension as VV, we have range⁡T=V\operatorname{range} T = V (a subspace of a finite-dimensional space that has the full dimension must be the whole space, since any basis of the subspace is a linearly independent set of size nn in VV and hence a basis of VV). Thus TT is surjective. Surjective ⇒\Rightarrow injective. If TT is surjective, range⁡T=V\operatorname{range} T = V, so dim⁡(range⁡T)=n\dim(\operatorname{range} T) = n. By Rank–Nullity, dim⁡(ker⁡T)=0\dim(\ker T) = 0, so ker⁡T={0}\ker T = \{0\}, which means TT is injective (linear maps with trivial kernel are injective: if T(u)=T(v)T(u) = T(v) then T(u−v)=0T(u - v) = 0, so u−v∈ker⁡T={0}u - v \in \ker T = \{0\}, so u=vu = v). ■\blacksquare

∎
Proof

Proof of the main theorem. We prove the chain (2)⇒(4)⇒(5)⇒(3)⇒(1)⇒(2)(2) \Rightarrow (4) \Rightarrow (5) \Rightarrow (3) \Rightarrow (1) \Rightarrow (2), and separately establish (6)(6)'s equivalence, which will make all six statements equivalent, and along the way we extract the two one-sided-inverse claims. Throughout, identify AA with the linear map TA:Fn→FnT_A: F^n \to F^n, TA(x)=AxT_A(x) = Ax; matrix multiplication corresponds to composition of these linear maps, and InI_n corresponds to the identity map id⁡\operatorname{id} (this correspondence is Matrix Multiplication as Composition, assumed known).

Step 1: (2)⇒(4)(2) \Rightarrow (4). Suppose BA=InBA = I_n. Let TB,TAT_B, T_A be the associated linear maps; then TB∘TA=id⁡T_B \circ T_A = \operatorname{id}. Suppose TA(x)=TA(y)T_A(x) = T_A(y) for some x,y∈Fnx, y \in F^n. Apply TBT_B to both sides:

TB(TA(x))=TB(TA(y))  ⟹  (TB∘TA)(x)=(TB∘TA)(y)  ⟹  x=y.T_B(T_A(x)) = T_B(T_A(y)) \implies (T_B \circ T_A)(x) = (T_B \circ T_A)(y) \implies x = y.

Hence TAT_A is injective.

Step 2: (4)⇒(5)(4) \Rightarrow (5). This is exactly the Lemma above, applied with V=FnV = F^n (finite-dimensional, dim⁡Fn=n\dim F^n = n) and T=TAT = T_A.

Step 3: (5)⇒(3)(5) \Rightarrow (3). Suppose TAT_A is surjective. We construct a right inverse matrix CC. Let e1,…,ene_1, \dots, e_n be the standard basis of FnF^n. Since TAT_A is surjective, for each j=1,…,nj = 1, \dots, n there exists cj∈Fnc_j \in F^n with Acj=ejA c_j = e_j. Let CC be the matrix whose jj-th column is cjc_j. By the column interpretation of matrix multiplication, the jj-th column of ACAC is Acj=ejA c_j = e_j. Since this holds for every jj, AC=InAC = I_n.

Step 4: (3)⇒(1)(3) \Rightarrow (1). This is the substantive step and the core of the theorem — we must upgrade a right inverse to a full two-sided inverse. Suppose AC=InAC = I_n. First, apply Steps 1–3 to the matrix CC in place of AA: from AC=InAC = I_n we read this as "AA is a left inverse of CC," i.e. CC satisfies condition (2) with the roles of A,BA, B replaced by C,AC, A (that is, AC=InAC = I_n plays the role of "BA=InBA = I_n" with B:=AB := A, A:=CA := C). By Steps 1–2 applied to CC, the map TCT_C is injective, hence (by the Lemma) surjective, hence by Step 3 there exists D∈Fn×nD \in F^{n\times n} with CD=InCD = I_n. Now we show D=AD = A. We have

A=AIn=A(CD)=(AC)D=InD=D.A = A I_n = A(CD) = (AC)D = I_n D = D.

Here we used associativity of matrix multiplication (which reflects associativity of function composition) and the hypothesis AC=InAC = I_n. So D=AD = A, and therefore CD=InCD = I_n becomes CA=InCA = I_n. Combining AC=InAC = I_n (given) with CA=InCA = I_n (just shown), CC is a genuine two-sided inverse of AA. Hence AA is invertible, i.e. condition (1) holds, with A−1=CA^{-1} = C.

Step 5: (1)⇒(2)(1) \Rightarrow (2). Trivial: if AA is invertible with inverse A−1A^{-1}, then B:=A−1B := A^{-1} satisfies BA=InBA = I_n. This closes the cycle (2)⇒(4)⇒(5)⇒(3)⇒(1)⇒(2)(2) \Rightarrow (4) \Rightarrow (5) \Rightarrow (3) \Rightarrow (1) \Rightarrow (2), so all of (1), (2), (3), (4), (5) are equivalent.

Recovering the one-sided statements. Suppose BA=InBA = I_n. By (2)⇒(1)(2)\Rightarrow(1) (the chain just proved), AA is invertible with some two-sided inverse A−1A^{-1}. Then

B=BIn=B(AA−1)=(BA)A−1=InA−1=A−1,B = BI_n = B(AA^{-1}) = (BA)A^{-1} = I_n A^{-1} = A^{-1},

so B=A−1B = A^{-1}, and in particular AB=AA−1=InAB = AA^{-1} = I_n. This proves the first displayed claim of the theorem. The dual claim (right inverse implies two-sided) is exactly Step 4 above, verbatim.

Step 6: (1)⇔(6)(1) \Leftrightarrow (6). Suppose AA is invertible. If ∑jxj(Aej)=0\sum_j x_j (Ae_j) = 0 for scalars xjx_j, i.e. Ax=0A x = 0 where x=(x1,…,xn)x = (x_1,\dots,x_n), apply A−1A^{-1}: x=A−1(Ax)=A−10=0x = A^{-1}(Ax) = A^{-1}0 = 0. So the columns Ae1,…,AenAe_1,\dots,Ae_n are linearly independent; being nn independent vectors in the nn-dimensional space FnF^n, they form a basis (a linearly independent set of size dim⁡V\dim V in a finite-dimensional space VV is automatically a basis, since it must also span: otherwise it could be extended to a larger independent set, contradicting that any independent set has size at most dim⁡V\dim V). Conversely, if the columns of AA form a basis of FnF^n, then for every y∈Fny \in F^n there are unique scalars expressing yy as their combination, i.e. Ax=yAx = y has a solution xx for every yy; this says TAT_A is surjective, and by the already-established equivalence, AA is invertible. All six conditions are therefore equivalent, and the two one-sided-inverse claims are established. ■\blacksquare

∎
CorollaryGLn(F)GL_n(F) is a group

Let FF be a field. The set GLn(F)GL_n(F) of invertible n×nn \times n matrices over FF, equipped with matrix multiplication, is a group: it is closed under multiplication, contains the identity InI_n, multiplication is associative, and every element has an inverse in GLn(F)GL_n(F).

Proof

Associativity holds for all matrices, invertible or not, hence in particular in GLn(F)GL_n(F) (matrix multiplication is always associative, inherited from associativity of function composition, Matrix Multiplication as Composition). Identity. In∈GLn(F)I_n \in GL_n(F) since InIn=InI_n I_n = I_n, so InI_n is its own inverse. Inverses. By definition, if A∈GLn(F)A \in GL_n(F) then A−1A^{-1} exists with AA−1=A−1A=InAA^{-1} = A^{-1}A = I_n; this same equation read with the roles of AA and A−1A^{-1} swapped shows A−1A^{-1} is invertible with inverse AA, so A−1∈GLn(F)A^{-1} \in GL_n(F). Closure. Let A,B∈GLn(F)A, B \in GL_n(F). We claim (AB)−1=B−1A−1(AB)^{-1} = B^{-1}A^{-1}, which in particular shows AB∈GLn(F)AB \in GL_n(F). Compute, using associativity:

(AB)(B−1A−1)=A(BB−1)A−1=AInA−1=AA−1=In,(AB)(B^{-1}A^{-1}) = A(BB^{-1})A^{-1} = A I_n A^{-1} = AA^{-1} = I_n,
(B−1A−1)(AB)=B−1(A−1A)B=B−1InB=B−1B=In.(B^{-1}A^{-1})(AB) = B^{-1}(A^{-1}A)B = B^{-1} I_n B = B^{-1}B = I_n.

So B−1A−1B^{-1}A^{-1} is a two-sided inverse of ABAB; by uniqueness of inverses (if AA is invertible, its inverse is unique: if B1A=AB1=InB_1A = AB_1 = I_n and B2A=AB2=InB_2 A = A B_2 = I_n then B1=B1In=B1(AB2)=(B1A)B2=InB2=B2B_1 = B_1 I_n = B_1(AB_2) = (B_1A)B_2 = I_n B_2 = B_2), (AB)−1=B−1A−1(AB)^{-1} = B^{-1}A^{-1} and AB∈GLn(F)AB \in GL_n(F). All group axioms hold, so GLn(F)GL_n(F) is a group under matrix multiplication. ■\blacksquare

∎

Worked examples

Example 1: Verifying invertibility via a genuine one-sided check. Let

A=(210111012)∈R3×3.A = \begin{pmatrix} 2 & 1 & 0 \\ 1 & 1 & 1 \\ 0 & 1 & 2 \end{pmatrix} \in \mathbb{R}^{3\times 3}.

We test invertibility by checking injectivity of TAT_A directly — i.e., showing Ax=0⇒x=0Ax = 0 \Rightarrow x = 0 — which by the theorem suffices, with no need to separately verify surjectivity or exhibit a two-sided inverse up front. Suppose A(x1,x2,x3)T=0A(x_1,x_2,x_3)^T = 0:

2x1+x2=0,x1+x2+x3=0,x2+2x3=0.2x_1 + x_2 = 0, \qquad x_1 + x_2 + x_3 = 0, \qquad x_2 + 2x_3 = 0.

From the first equation, x2=−2x1x_2 = -2x_1. Substitute into the third: −2x1+2x3=0-2x_1 + 2x_3 = 0, so x3=x1x_3 = x_1. Substitute both into the second: x1−2x1+x1=0x_1 - 2x_1 + x_1 = 0, which gives 0=00 = 0 — no new information, and the system is consistent for any x1x_1. Set x1=tx_1 = t: then x2=−2tx_2 = -2t, x3=tx_3 = t. Substitute back into the original equations to double check: first equation 2t+(−2t)=02t + (-2t) = 0 ✓; third equation (−2t)+2t=0(-2t) + 2t = 0 ✓; second: t+(−2t)+t=0t + (-2t) + t = 0 ✓. So every choice of tt gives a solution (t,−2t,t)(t, -2t, t). Taking t=1t = 1 gives the nonzero vector (1,−2,1)T(1,-2,1)^T with A(1,−2,1)T=0A(1,-2,1)^T = 0. Hence ker⁡TA≠{0}\ker T_A \ne \{0\}, so TAT_A is not injective. By the theorem, AA is not invertible — no right or left inverse can exist, without our needing to search for one. As a cross-check via determinant (previewing Rank of a Matrix and cofactor expansion):

det⁡A=2(1⋅2−1⋅1)−1(1⋅2−1⋅0)+0=2(1)−1(2)+0=0,\det A = 2(1\cdot 2 - 1\cdot 1) - 1(1\cdot 2 - 1\cdot 0) + 0 = 2(1) - 1(2) + 0 = 0,

confirming singularity independently.

Example 2: Exhibiting a one-sided inverse and getting the other side for free. Let

A=(1203),B=(1−2/301/3)∈Q2×2.A = \begin{pmatrix} 1 & 2 \\ 0 & 3 \end{pmatrix}, \qquad B = \begin{pmatrix} 1 & -2/3 \\ 0 & 1/3 \end{pmatrix} \in \mathbb{Q}^{2\times 2}.

Compute BABA directly:

BA=(1−2/301/3)(1203)=(1⋅1+(−2/3)⋅01⋅2+(−2/3)⋅30⋅1+(1/3)⋅00⋅2+(1/3)⋅3)=(1001).BA = \begin{pmatrix} 1 & -2/3 \\ 0 & 1/3 \end{pmatrix}\begin{pmatrix} 1 & 2 \\ 0 & 3 \end{pmatrix} = \begin{pmatrix} 1\cdot 1 + (-2/3)\cdot 0 & 1\cdot 2 + (-2/3)\cdot 3 \\ 0\cdot 1 + (1/3)\cdot 0 & 0\cdot 2 + (1/3)\cdot 3 \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}.

So BA=I2BA = I_2: BB is a left inverse of AA. By the theorem's one-sided-inverse claim, AA is automatically invertible and AB=I2AB = I_2 as well, with B=A−1B = A^{-1}. We verify this directly rather than taking it on faith:

AB=(1203)(1−2/301/3)=(1+0−2/3+2/30+00+1)=(1001)=I2.AB = \begin{pmatrix} 1 & 2 \\ 0 & 3 \end{pmatrix}\begin{pmatrix} 1 & -2/3 \\ 0 & 1/3 \end{pmatrix} = \begin{pmatrix} 1 + 0 & -2/3 + 2/3 \\ 0 + 0 & 0 + 1 \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} = I_2.

Indeed AB=I2AB = I_2 too, exactly as the theorem predicted — checking one side was already conclusive, and the arithmetic above is the concrete instance of the general argument in Step 4 of the proof (A=A(CD)=(AC)D=DA = A(CD) = (AC)D = D type identity, here with C=BC = B).

CounterexampleWhy finite dimensionality is needed

Drop the hypothesis that the space is finite-dimensional and one-sided inverses can genuinely fail to be two-sided. Let V=FNV = F^{\mathbb{N}}, the space of all sequences (a0,a1,a2,… )(a_0, a_1, a_2, \dots) of elements of a field FF (an infinite-dimensional vector space). Define the right shift S:V→VS: V \to V and left shift T:V→VT: V \to V by

S(a0,a1,a2,… )=(0,a0,a1,… ),T(a0,a1,a2,… )=(a1,a2,a3,… ).S(a_0, a_1, a_2, \dots) = (0, a_0, a_1, \dots), \qquad T(a_0, a_1, a_2, \dots) = (a_1, a_2, a_3, \dots).

Both are linear. Compute the composition T∘ST \circ S:

T(S(a0,a1,… ))=T(0,a0,a1,… )=(a0,a1,… ),T(S(a_0,a_1,\dots)) = T(0,a_0,a_1,\dots) = (a_0, a_1, \dots),

so T∘S=id⁡VT \circ S = \operatorname{id}_V: here TT is a left inverse of SS (equivalently SS is a right inverse of TT). But the other composition:

S(T(a0,a1,a2,… ))=S(a1,a2,a3,… )=(0,a1,a2,… )≠(a0,a1,a2,… )S(T(a_0,a_1,a_2,\dots)) = S(a_1,a_2,a_3,\dots) = (0, a_1, a_2, \dots) \neq (a_0,a_1,a_2,\dots)

in general — for instance with a0=1a_0 = 1 and all other entries 00, T∘ST\circ S fixes the sequence but S∘TS \circ T sends it to (0,0,0,… )≠(1,0,0,… )(0,0,0,\dots) \ne (1,0,0,\dots). So S∘T≠id⁡VS \circ T \ne \operatorname{id}_V. Thus SS has a left inverse (TT) but SS itself has no right inverse making it two-sided invertible: SS is injective but not surjective (nothing maps to (1,0,0,… )(1,0,0,\dots) under SS since every output of SS has first coordinate 00), and TT is surjective but not injective (its kernel contains all sequences of the form (a0,0,0,… )(a_0,0,0,\dots)). This is exactly the failure mode the finite-dimensional theorem rules out: injectivity and surjectivity split apart once VV is infinite-dimensional, because the Rank–Nullity dimension count that forced them together no longer applies. The hypothesis "AA is a square matrix acting on the finite-dimensional space FnF^n" is therefore not a technicality — it is precisely what makes the one-sided-inverse theorem true.

Scope, edge cases & common traps

SituationWhat happensWhy
AA square, FnF^n finite-dim.One-sided inverse ⇒\Rightarrow two-sidedRank–Nullity forces injective ⇔\Leftrightarrow surjective
AA non-square (m×nm \times n, m≠nm \ne n)One-sided inverses can exist without the other sideDomain/codomain dimensions differ; no contradiction from Rank–Nullity
Infinite-dimensional VVOne-sided inverses need not be two-sidedShift operators; see counterexample
FF not a field but a general ringTheorem can failRank–Nullity and basis arguments use field structure (division by nonzero scalars)
det⁡A=0\det A = 0AA singular, no inverse of either kindEquivalent condition (6): columns dependent

The mathematician behind it

Origins · 1 minArthur Cayley and how this was discovered
AC
Arthur Cayley
1821–1895 · Victorian algebra, matrix theory

Arthur Cayley, a British mathematician and one of the most prolific algebraists of the nineteenth century, is credited with founding matrix theory as an independent algebraic subject. Working initially as a lawyer to support himself before securing an academic post at Cambridge, Cayley treated matrices as objects to be added and multiplied in their own right, rather than mere arrays of coefficients for a linear system. His 1858 memoir on the theory of matrices introduced matrix multiplication, the identity matrix, and — crucially for this article — the notion of a matrix inverse, alongside an early version of what is now called the Cayley–Hamilton theorem.

History
- 1858 — Cayley's "A Memoir on the Theory of Matrices" formalizes matrix algebra, multiplication, and inverses as objects of independent study.
- 1878 — Georg Frobenius develops rank and canonical forms, sharpening the equivalence between invertibility and non-vanishing determinant.
- 1888 — Giuseppe Peano gives the first fully axiomatic definition of a vector space, providing the setting in which "finite-dimensional" becomes precise.
- 1920s — the abstract group-theoretic viewpoint (Emmy Noether's structural algebra) recasts GLn(F)GL_n(F) as a fundamental example of a group acting on a vector space.
- 1940s–1950s — numerical linear algebra (Turing, von Neumann) turns invertibility testing and matrix inversion into central problems of scientific computing.
Does a square matrix need both AB=IAB=I and BA=IBA=I checked to confirm invertibility?

No — this is the central point of this article. For a square matrix AA over a field acting on a finite-dimensional space, if you find any BB with AB=IAB = I or any BB with BA=IBA = I, that single equation already forces AA to be invertible with A−1=BA^{-1} = B, and the other equation follows automatically. The proof runs through Rank–Nullity: one-sided inverse gives injectivity or surjectivity of TAT_A, and on a finite-dimensional space these two properties are equivalent for maps from a space to itself, so bijectivity — and hence full invertibility — comes for free.

Why does this one-sided shortcut fail for non-square matrices?

Because Rank–Nullity in that setting reads dim⁡(ker⁡T)+dim⁡(range⁡T)=n\dim(\ker T) + \dim(\operatorname{range} T) = n where nn is the dimension of the domain, while the codomain may have a different dimension mm. If m≠nm \ne n, injectivity and surjectivity are no longer forced to coincide — an m×nm\times n matrix with m>nm > n can be injective without being surjective (its image is a proper subspace of the larger codomain), and one with m<nm < n can be surjective without injective. One-sided inverses for rectangular matrices exist in different senses (left vs. right pseudo-inverses) but are genuinely different objects.

Is GLn(F)GL_n(F) abelian?

Only when n=1n = 1 (or in degenerate cases like FF having very few elements combined with small nn, but generically no). For n≥2n \ge 2, matrix multiplication is non-commutative: for instance swapping rows via one elementary matrix then scaling via another gives a different result depending on order. So GLn(F)GL_n(F) is a genuinely non-abelian group for n≥2n\ge2, which is part of why it plays such a rich role as a source of examples in group theory and representation theory.

How does the determinant relate to this theorem?

The determinant gives a single-number test: AA is invertible if and only if det⁡A≠0\det A \ne 0. This is consistent with, and can be derived alongside, the equivalences in the theorem — e.g., columns of AA form a basis if and only if they are linearly independent if and only if det⁡A≠0\det A \ne 0. The determinant is often the fastest hand-computation route to check invertibility for small matrices, while the injective/surjective formulation is the conceptual reason the determinant test works at all.

Does the theorem hold over rings that are not fields, like Z\mathbb{Z}?

Not in general, and the argument used here relies on field-specific facts, particularly the basis-extension and Rank–Nullity arguments valid for vector spaces over a field. Over Z\mathbb{Z}, for instance, invertibility as a matrix over Z\mathbb{Z} additionally requires det⁡A=±1\det A = \pm 1 (not merely nonzero), because you need 1/det⁡A1/\det A to lie back in Z\mathbb{Z}. The general theory of modules over rings replaces "basis" with more delicate notions (free modules, torsion), and one-sided invertibility results become subtler.

If AA and BB are both invertible n×nn\times n matrices, is A+BA+B necessarily invertible?

No. Invertibility is not preserved under addition. For example, with A=I2A = I_2 and B=−I2B = -I_2, both are invertible, but A+B=0A + B = 0, the zero matrix, which is singular (it sends every vector to 00, so it is not injective for n≥1n \ge 1). Invertibility is a multiplicative property, tied to the group structure of GLn(F)GL_n(F), and there is no reason to expect it to survive under the unrelated operation of addition.

What is the fastest way in practice to test if a matrix is invertible?

Row reduce AA to row echelon form using elementary row operations (see Elementary Matrices and Row Operations); AA is invertible exactly when this process produces nn nonzero pivot rows, equivalently when the reduced row echelon form is InI_n. This is computationally the standard method, and it simultaneously computes A−1A^{-1} if you augment with the identity matrix and row-reduce the pair [A∣I][A \mid I] to [I∣A−1][I \mid A^{-1}].

Before this

Where this leads