Math graphic
πŸ“ Concept diagram

08-07 β€” Determinants (Deep)

Phase: 8 β€” Linear Algebra (Rigorous) Subject: 08-07 Prerequisites: 08-06 β€” Orthogonal Projections Next subject: 08-08 β€” Eigenvalues and Eigenvectors


Learning Objectives

By the end of this subject, you will be able to:

  1. Define the determinant of an n Γ— n matrix via permutations and parity (Leibniz formula), and compute small determinants by hand
  2. Compute determinants using cofactor expansion (Laplace expansion) along any row or column
  3. State and apply the key properties of determinants: multilinearity, alternating property, multiplicativity, and effect of row operations
  4. Use Cramer's rule to solve small linear systems and understand when it applies
  5. Interpret the determinant geometrically as signed volume/area in R^n and use it to test invertibility and orientation

Core Content

1. Definition via Permutations (Leibniz Formula)

The determinant of an n Γ— n matrix A = [a_ij] is:

$det(A) = Ξ£_{Οƒ ∈ S_n} sgn(Οƒ) Β· a_{1,Οƒ(1)} a_{2,Οƒ(2)} ... a_{n,Οƒ(n)}

> **CRITICAL β€” Foundational:** The determinant is uniquely defined by: (1) det(I) = 1, (2) swapping rows flips sign, (3) linear in each row. These three imply ALL other properties.
$

where: - S_n is the set of all n! permutations of {1, 2, ..., n} - sgn(Οƒ) = (+1 if Οƒ is even, βˆ’1 if Οƒ is odd) - A permutation is even/odd based on the number of transpositions (swaps) needed to return it to the identity

Examples: - n = 1: det([a]) = a - n = 2: det([a b; c d]) = aΒ·d βˆ’ bΒ·c (Sβ‚‚ = {identity (even), swap (odd)}) - n = 3: det([a₁₁ a₁₂ a₁₃; a₂₁ aβ‚‚β‚‚ a₂₃; a₃₁ a₃₂ a₃₃]) = a₁₁aβ‚‚β‚‚a₃₃ + a₁₂a₂₃a₃₁ + a₁₃a₂₁a₃₂ βˆ’ a₁₁a₂₃a₃₂ βˆ’ a₁₂a₂₁a₃₃ βˆ’ a₁₃aβ‚‚β‚‚a₃₁ (6 terms: 3 even permutations, 3 odd)

Parity (sign) of permutations: An inversion is a pair (i, j) with i < j but Οƒ(i) > Οƒ(j). Οƒ is even if it has an even number of inversions; odd otherwise. Equivalently, sgn(Οƒ) = (βˆ’1)^(#inversions).

Example: Οƒ = (2, 3, 1). Inversions: (2,1), (3,1) β†’ 2 inversions β†’ sgn = +1 (even).

2. Cofactor Expansion (Laplace Expansion)

The (i, j) minor M_ij is the determinant of the (nβˆ’1) Γ— (nβˆ’1) submatrix obtained by deleting row i and column j.

The (i, j) cofactor is C_ij = (βˆ’1)^{i+j} M_ij.

Cofactor expansion along any row i:

$det(A) = Ξ£_{j=1}^n a_ij C_ij
$

Or along any column j:

$det(A) = Ξ£_{i=1}^n a_ij C_ij
$

Example (3 Γ— 3):

$A = [1  2  3]
    [4  5  6]
    [7  8  9]
$

Expand along row 1: det(A) = 1Β·C₁₁ + 2Β·C₁₂ + 3Β·C₁₃ = 1Β·(+1)Β·det[5 6; 8 9] + 2Β·(βˆ’1)Β·det[4 6; 7 9] + 3Β·(+1)Β·det[4 5; 7 8] = 1Β·(45βˆ’48) βˆ’ 2Β·(36βˆ’42) + 3Β·(32βˆ’35) = 1Β·(βˆ’3) βˆ’ 2Β·(βˆ’6) + 3Β·(βˆ’3) = βˆ’3 + 12 βˆ’ 9 = 0.

3. Properties of Determinants

Fundamental properties:

  1. det(I) = 1

  2. Row interchange: Swapping two rows multiplies det by βˆ’1.

  3. Multilinearity (rows): det is linear in each row when others are fixed. det(..., cΒ·row_i, ...) = cΒ·det(..., row_i, ...) det(..., u+v, ...) = det(..., u, ...) + det(..., v, ...)

  4. Alternating: If two rows are equal, det = 0.

  5. Row operation (add multiple of row): Adding a multiple of one row to another leaves det unchanged.

  6. Triangular matrix: det = product of diagonal entries.

  7. Singularity: det(A) = 0 iff A is singular (not invertible). Equivalently, det(A) β‰  0 iff columns/rows are independent.

  8. Multiplicativity: det(AB) = det(A)Β·det(B)

  9. Transpose: det(A^T) = det(A)

  10. Inverse: det(A⁻¹) = 1/det(A) (when A is invertible)

  11. Block triangular: det([A B; 0 C]) = det(A)Β·det(C)

Effect of row operations on determinant: - Row swap: det β†’ βˆ’det - Row scaling by c: det β†’ cΒ·det - Row addition: det unchanged

Example using row operations:

$A = [1  2  3]
    [2  1  3]
    [3  3  6]
$

Rβ‚‚βˆ’2R₁: [1 2 3; 0 βˆ’3 βˆ’3; 3 3 6] Rβ‚ƒβˆ’3R₁: [1 2 3; 0 βˆ’3 βˆ’3; 0 βˆ’3 βˆ’3] Rβ‚ƒβˆ’Rβ‚‚: [1 2 3; 0 βˆ’3 βˆ’3; 0 0 0] Triangular: det = 1Β·(βˆ’3)Β·0 = 0. A is singular. βœ“

4. Cramer's Rule

For an n Γ— n invertible matrix A, the unique solution to Ax = b has components:

$x_i = det(A_i) / det(A)
$

where A_i is A with column i replaced by b.

Why it works: Write x = A⁻¹b. By the adjugate formula A⁻¹ = adj(A)/det(A). Then x_i = (1/det(A)) Σ_j adj(A)_ij b_j = det(A_i)/det(A). (The sum is precisely the cofactor expansion of A with column i replaced by b.)

Example:

$2x + y = 5
x βˆ’ 3y = βˆ’1
$

A = [2 1; 1 βˆ’3], det(A) = βˆ’6 βˆ’ 1 = βˆ’7.

x = det([5 1; βˆ’1 βˆ’3]) / (βˆ’7) = (βˆ’15 βˆ’ (βˆ’1))/(βˆ’7) = (βˆ’14)/(βˆ’7) = 2. y = det([2 5; 1 βˆ’1]) / (βˆ’7) = (βˆ’2 βˆ’ 5)/(βˆ’7) = (βˆ’7)/(βˆ’7) = 1.

Limitations: Cramer's rule requires O(n!) operations. It's computationally impractical for n > 4 but theoretically useful.

Common Pitfall: Cramer's rule is elegant but EXPENSIVE β€” never use it for n > 4. Use Gaussian elimination (O(nΒ³)) for practical computation.

5. Geometric Interpretation

In R^n, the absolute value |det(A)| is the n-dimensional volume of the parallelepiped spanned by the column vectors of A.

Sign: det(A) > 0 means the columns form a right-handed (positively oriented) basis; det(A) < 0 means left-handed.

Area/volume after linear transformation: If T(x) = Ax, then volume(T(S)) = |det(A)|Β·volume(S) for any region S.

Connection to independence: det(A) β‰  0 iff columns are independent iff the parallelepiped has nonzero n-dimensional volume.



Key Terms

Worked Examples

Example 1: Determinant by Cofactor Expansion

Problem: Compute det(A) where A = [2 βˆ’1 0; 3 1 2; 1 0 βˆ’1].

Solution: Expand along row 3 (has a zero for convenience).

$det(A) = 1Β·C₃₁ + 0Β·C₃₂ + (βˆ’1)Β·C₃₃
$

C₃₁ = (+)det[βˆ’1 0; 1 2] = 1Β·(βˆ’1Β·2 βˆ’ 0Β·1) = βˆ’2. (sign: (βˆ’1)^{3+1} = +1) C₃₃ = (+)det[2 βˆ’1; 3 1] = 1Β·(2Β·1 βˆ’ (βˆ’1)Β·3) = 5. (sign: (βˆ’1)^{3+3} = +1)

det(A) = 1Β·(βˆ’2) + (βˆ’1)Β·5 = βˆ’2 βˆ’ 5 = βˆ’7.

Check via row operations: Rβ‚‚ βˆ’ (3/2)R₁: [2 βˆ’1 0; 0 5/2 2; 1 0 βˆ’1] R₃ βˆ’ (1/2)R₁: [2 βˆ’1 0; 0 5/2 2; 0 1/2 βˆ’1] R₃ βˆ’ (1/5)Rβ‚‚: [2 βˆ’1 0; 0 5/2 2; 0 0 βˆ’7/5] det = 2 Β· (5/2) Β· (βˆ’7/5) = βˆ’7. βœ“

Example 2: Determinant Properties in Action

Problem: Without computing, find the relationship between det(A) and det(B): A = [a b c; d e f; g h i] B = [a b c; 2d 2e 2f; g h i]

Solution: Row 2 of B is 2 Γ— row 2 of A. By multilinearity, det(B) = 2Β·det(A).

Another: If C = [a b c; d+a e+b f+c; g h i] (row 2 = original row 2 + row 1 of A), then det(C) = det(A) since adding a multiple of one row to another doesn't change the determinant.

Another: If D = [a b c; g h i; d e f] (rows 2 and 3 swapped), then det(D) = βˆ’det(A).

Example 3: Cramer's Rule

Problem: Solve using Cramer's rule:

$x + 2y + z = 5
2x βˆ’ y + 3z = 7
x + y + z = 4
$

Solution: A = [1 2 1; 2 βˆ’1 3; 1 1 1].

det(A) = 1Β·det[βˆ’1 3; 1 1] βˆ’ 2Β·det[2 3; 1 1] + 1Β·det[2 βˆ’1; 1 1] = 1Β·(βˆ’1βˆ’3) βˆ’ 2Β·(2βˆ’3) + 1Β·(2βˆ’(βˆ’1)) = βˆ’4 βˆ’ 2Β·(βˆ’1) + 3 = βˆ’4 + 2 + 3 = 1.

x = det([5 2 1; 7 βˆ’1 3; 4 1 1])/1 = 5Β·det[βˆ’1 3; 1 1] βˆ’ 2Β·det[7 3; 4 1] + 1Β·det[7 βˆ’1; 4 1] = 5Β·(βˆ’4) βˆ’ 2Β·(7βˆ’12) + 1Β·(7βˆ’(βˆ’4)) = βˆ’20 βˆ’ 2Β·(βˆ’5) + 11 = βˆ’20 + 10 + 11 = 1.

y = det([1 5 1; 2 7 3; 1 4 1])/1 = 1Β·det[7 3; 4 1] βˆ’ 5Β·det[2 3; 1 1] + 1Β·det[2 7; 1 4] = 1Β·(7βˆ’12) βˆ’ 5Β·(2βˆ’3) + 1Β·(8βˆ’7) = βˆ’5 βˆ’ 5Β·(βˆ’1) + 1 = βˆ’5 + 5 + 1 = 1.

z = det([1 2 5; 2 βˆ’1 7; 1 1 4])/1 = 1Β·det[βˆ’1 7; 1 4] βˆ’ 2Β·det[2 7; 1 4] + 5Β·det[2 βˆ’1; 1 1] = 1Β·(βˆ’4βˆ’7) βˆ’ 2Β·(8βˆ’7) + 5Β·(2+1) = βˆ’11 βˆ’ 2 + 15 = 2.

Solution: (x, y, z) = (1, 1, 2). Check: 1+2+2=5, 2βˆ’1+6=7, 1+1+2=4. βœ“

Example 4: Geometric Interpretation

Problem: Find the area of the triangle with vertices (0,0), (3,1), (1,4).

Solution: The triangle has half the area of the parallelogram spanned by v = (3,1) and w = (1,4).

Area = (1/2) |det([3 1; 1 4])| = (1/2) |3Β·4 βˆ’ 1Β·1| = (1/2)Β·11 = 5.5.

The parallelogram has area 11.

Quiz

Q1: How many terms are in the Leibniz expansion of a 4 Γ— 4 determinant?

A) 4 B) 16 C) 24 D) 256

Correct: C)


Q2: The cofactor C_ij of matrix entry a_ij includes the sign factor:

A) (βˆ’1)^i B) (βˆ’1)^j C) (βˆ’1)^{i+j} D) (βˆ’1)^{iΒ·j}

Correct: C)


Q3: If you swap two rows of a matrix, what happens to its determinant?

A) It stays the same B) It is multiplied by βˆ’1 C) It becomes 0 D) It doubles

Correct: B)


Q4: For an n Γ— n matrix A, det(cA) equals:

A) cΒ·det(A) B) c^nΒ·det(A) C) nΒ·det(A) D) det(A)/c

Correct: B)


Q5: Cramer's rule for solving Ax = b gives the i-th component as:

A) det(A_i)/det(A) where A_i is A with column i replaced by b B) det(A)/det(A_i) C) det(A)Β·det(A_i) D) det(A_i) βˆ’ det(A)

Correct: A)


Q6: Geometrically, |det(A)| for a 2 Γ— 2 matrix A represents:

A) The perimeter of the parallelogram spanned by its columns B) The area of the parallelogram spanned by its columns C) The angle between its columns D) The sum of the lengths of its columns

Correct: B)


Practice Problems

(Answers are below. Try each problem before checking.)

Problem 1: State the sign of each permutation of {1, 2, 3}: (a) (1, 2, 3), (b) (2, 3, 1), (c) (3, 1, 2), (d) (3, 2, 1).

Problem 2: Compute det([2 1 3; 4 βˆ’1 2; 1 0 5]) using cofactor expansion.

Problem 3: Use determinant properties to compute det([a b c; d e f; a+g b+h c+i]) in terms of other determinants.

Problem 4: Solve using Cramer's rule: 3x + y = 5, 2x βˆ’ y = 0.

Problem 5: Find the volume of the parallelepiped spanned by (1, 0, 1), (2, 1, 0), and (0, 1, 3).

Problem 6: Prove that det(A⁻¹) = 1/det(A) for invertible A.

Problem 7: Show that det(cA) = c^n det(A) for an n Γ— n matrix A.

Answers (click to expand) **Problem 1:** (a) Identity: 0 inversions β†’ sgn = +1 (even) (b) (2,3,1): inversions (2,1), (3,1) β†’ 2 β†’ sgn = +1 (even) (c) (3,1,2): inversions (3,1), (3,2) β†’ 2 β†’ sgn = +1 (even) (d) (3,2,1): inversions (3,2), (3,1), (2,1) β†’ 3 β†’ sgn = βˆ’1 (odd) Note: (b) and (c) are even β€” they are 3-cycles. All 3-cycles have sgn = +1. **Problem 2:** Expand along row 3 (has a zero): det = 1Β·C₃₁ + 0Β·C₃₂ + 5Β·C₃₃ C₃₁ = (+)det[1 3; βˆ’1 2] = 1Β·(2 βˆ’ (βˆ’3)) = 5. C₃₃ = (+)det[2 1; 4 βˆ’1] = 1Β·(βˆ’2 βˆ’ 4) = βˆ’6. det = 1Β·5 + 5Β·(βˆ’6) = 5 βˆ’ 30 = βˆ’25. **Problem 3:** Row 3 = (a+g, b+h, c+i) = (a,b,c) + (g,h,i) = row 1 + (g,h,i). By multilinearity: det([a b c; d e f; a+g b+h c+i]) = det([a b c; d e f; a b c]) + det([a b c; d e f; g h i]). The first term has rows 1 and 3 equal β†’ det = 0. So the result = det([a b c; d e f; g h i]). **Problem 4:** A = [3 1; 2 βˆ’1], det(A) = βˆ’3 βˆ’ 2 = βˆ’5. A₁ = [5 1; 0 βˆ’1], det = βˆ’5. Aβ‚‚ = [3 5; 2 0], det = βˆ’10. x = βˆ’5/(βˆ’5) = 1, y = βˆ’10/(βˆ’5) = 2. Solution: (1, 2). Check: 3+2=5, 2βˆ’2=0. βœ“ **Problem 5:** A = [1 2 0; 0 1 1; 1 0 3]. det(A) = 1Β·det[1 1; 0 3] βˆ’ 2Β·det[0 1; 1 3] + 0Β·... = 1Β·3 βˆ’ 2Β·(βˆ’1) = 3 + 2 = 5. Volume = |5| = 5. **Problem 6:** det(A A⁻¹) = det(I) = 1. By multiplicativity: det(A)Β·det(A⁻¹) = 1. Since det(A) β‰  0 (A is invertible), det(A⁻¹) = 1/det(A). **Problem 7:** cA has each row multiplied by c. By multilinearity, each row contributes a factor c, and there are n rows. So det(cA) = c^n det(A). Alternative via definition: each term in the Leibniz sum has n factors, each multiplied by c, giving c^n.

Summary

  1. The determinant is defined by the Leibniz formula summing over all n! permutations with signs determined by parity; it's the unique multilinear alternating form with det(I) = 1
  2. Cofactor expansion provides a recursive computation method: pick a row or column, compute minors and cofactors with signs (βˆ’1)^{i+j}
  3. Determinants have powerful algebraic properties: det(AB) = det(A)det(B), row swap negates det, scaling a row scales det, and two equal rows give det = 0
  4. Cramer's rule gives explicit formulas for solving Ax = b when A is invertible: x_i = det(A_i)/det(A); elegant in theory but computationally expensive
  5. Geometrically, |det(A)| is the n-volume of the parallelepiped spanned by its columns; det(A) = 0 ⇔ columns are linearly dependent ⇔ zero n-volume

Pitfalls

  1. Assuming det(A + B) = det(A) + det(B). The determinant is multiplicative, NOT additive. In general det(A + B) β‰  det(A) + det(B). A simple counterexample: A = I, B = βˆ’I, where det(I) + det(βˆ’I) = 1 + 1 = 2 but det(I + (βˆ’I)) = det(0) = 0.

  2. Getting cofactor signs wrong. The cofactor sign is (βˆ’1)^{i+j}, creating a checkerboard pattern. A common error is using (βˆ’1)^i or (βˆ’1)^j alone. For a 3Γ—3 expansion, misplacing a single sign flips the sign of the entire determinant.

  3. Forgetting that det(cA) = c^n det(A), not cΒ·det(A). Each of the n rows is scaled by c, so the determinant picks up a factor of c^n. For a 3Γ—3 matrix, det(2A) = 8Β·det(A), not 2Β·det(A).

  4. Using Cramer's rule numerically for n > 3. Cramer's rule requires computing n+1 determinants of size n, each O(n!). For n = 10, this is ~4Γ—10^8 operations vs. ~10^3 for Gaussian elimination. Cramer's rule is a theoretical tool; use elimination for computation.

  5. Confusing the effect of row operations on the determinant. Row swaps multiply det by βˆ’1. Row scaling multiplies det by the scale factor. Row addition (adding a multiple of one row to another) leaves det UNCHANGED β€” but only if you add to a DIFFERENT row. Adding a multiple of a row to itself doubles it, which scales det.



Next Steps

Move on to 08-08 β€” Eigenvalues and Eigenvectors to discover the characteristic equation, how to find eigenvalues and eigenvectors, the distinction between algebraic and geometric multiplicity, and the conditions for diagonalization.