Every problem in the chapter beside the book’s answer, a chapter at a time as you scroll.

Chapter 4: Determinants

145 problems · 11 figures

Four.I.1.1 §I.1 source ↗

Evaluate the determinant of each.

(a)

(31−11)

(b)

(201311−101)

(c)

(40100113−1)

Answer (3 parts)

(a)

4

(b)

3

(c)

−12

Four.I.1.2 §I.1 source ↗

Evaluate the determinant of each.

(a)

(20−13)

(b)

(21105−21−34)

(c)

(234567891)

Answer (3 parts)

(a)

6

(b)

21

(c)

27

Four.I.1.3 §I.1 source ↗

Verify that the determinant of an upper-triangular 3×3 matrix is the product down the diagonal.

det⁡((abc0ef00i))=aei

Do lower-triangular matrices work the same way?

Answer

For the first, apply the formula in this section, note that any term with a d⁠, g⁠, or h is zero, and simplify. Lower-triangular matrices work the same way.

Four.I.1.4 §I.1 source ↗

Use the determinant to decide if each is singular or nonsingular.

(a)

(2131)

(b)

(011−1)

(c)

(4221)

Answer (3 parts)

(a)

Nonsingular, the determinant is −1⁠.

(b)

Nonsingular, the determinant is −1⁠.

(c)

Singular, the determinant is 0⁠.

Four.I.1.5 §I.1 source ↗

Singular or nonsingular? Use the determinant to decide.

(a)

(211322014)

(b)

(101211413)

(c)

(2103−20100)

Answer (3 parts)

(a)

Nonsingular, the determinant is 3⁠.

(b)

Singular, the determinant is 0⁠.

(c)

Singular, the determinant is 0⁠.

Four.I.1.6 §I.1 source ↗

Each pair of matrices differ by one row operation. Use this operation to compare det⁡(A) with det⁡(B)⁠.

(a)

A=(1223)⁠, B=(120−1)

(b)

A=(310001012)⁠, B=(310012001)

(c)

A=(1−1322−6104)⁠, B=(1−1311−3104)

Answer (3 parts)

(a)

det⁡(B)=det⁡(A) via −2ρ1+ρ2

(b)

det⁡(B)=−det⁡(A) via ρ2↔ρ3

(c)

det⁡(B)=(1/2)⋅det⁡(A) via (1/2)ρ2

Four.I.1.7 §I.1 source ↗

Find the determinant of this 4×4 matrix by following the plan: perform Gauss’s Method and look for the determinant to remain unchanged on a row combination, to change sign on a row swap, to rescale on the rescaling of a row, and such that the determinant of the echelon form matrix is the product down its diagonal.

(1202241000−133−114)

Answer

Gauss’s Method does this.

(1202241000−133−114)⟶−3ρ1+ρ4−2ρ1+ρ2⟶ρ2↔ρ4⟶ρ3+ρ4(12020−71−200−13000−1)

The echelon form matrix has a product down the diagonal of 1⋅(−7)⋅(−1)⋅(−1)=−7⁠. In the course of Gauss’s Method no rows got rescaled but there was a row swap, so to get the determinant we change the sign, giving +7⁠.

Four.I.1.8 §I.1 source ↗

Show this.

det⁡((111abca2b2c2))=(b−a)(c−a)(c−b)

Answer

Using the formula for the determinant of a 3×3 matrix we expand the left side

1⋅b⋅c2+1⋅c⋅a2+1⋅a⋅b2−b2⋅c⋅1−c2⋅a⋅1−a2⋅b⋅1

and by distributing we expand the right side.

(bc−ba−ac+a2)⋅(c−b)=c2b−b2c−bac+b2a−ac2+acb+a2c−a2b

Now we can just check that the two are equal. (Remark. This is the 3×3 case of Vandermonde’s determinant which arises in applications).

Four.I.1.9 §I.1 source ↗

Which real numbers x make this matrix singular?

(12−x488−x)

Answer

This equation

0=det⁡((12−x488−x))=64−20x+x2=(x−16)(x−4)

has roots x=16 and x=4⁠.

Four.I.1.10 §I.1 source ↗

Do the Gaussian reduction to check the formula for 3×3 matrices stated in the preamble to this section.

(abcdefghi) is nonsingular iff aei+ bfg+ cdh− hfa− idb− gec≠0

Answer

We first reduce the matrix to echelon form. To begin, assume that a≠0 and that ae−bd≠0⁠.

⟶(1/a)ρ1(1b/ac/adefghi)⟶−gρ1+ρ3−dρ1+ρ2(1b/ac/a0(ae−bd)/a(af−cd)/a0(ah−bg)/a(ai−cg)/a)⟶(a/(ae−bd))ρ2(1b/ac/a01(af−cd)/(ae−bd)0(ah−bg)/a(ai−cg)/a)

This step finishes the calculation.

⟶((ah−bg)/a)ρ2+ρ3(1b/ac/a01(af−cd)/(ae−bd)00(aei+bgf+cdh−hfa−idb−gec)/(ae−bd))

Now assuming that a≠0 and ae−bd≠0⁠, the original matrix is nonsingular if and only if the 3,3 entry above is nonzero. That is, under the assumptions, the original matrix is nonsingular if and only if aei+ bgf+ cdh− hfa− idb− gec≠0⁠, as required.

We finish by running down what happens if the assumptions that were taken for convenience in the prior paragraph do not hold. First, if a≠0 but ae−bd=0 then we can swap

(1b/ac/a00(af−cd)/a0(ah−bg)/a(ai−cg)/a)⟶ρ2↔ρ3(1b/ac/a0(ah−bg)/a(ai−cg)/a00(af−cd)/a)

and conclude that the matrix is nonsingular if and only if either ah−bg=0 or af−cd=0⁠. The condition ‘⁠ah−bg=0 or af−cd=0⁠’ is equivalent to the condition ‘⁠(ah−bg)(af−cd)=0⁠’. Multiplying out and using the case assumption that ae−bd=0 to substitute ae for bd gives this.

0=ahaf−ahcd−bgaf+bgcd=ahaf−ahcd−bgaf+aegc=a(haf−hcd−bgf+egc)

Since a≠0⁠, we have that the matrix is nonsingular if and only if haf−hcd−bgf+egc=0⁠. Therefore, in this a≠0 and ae−bd=0 case, the matrix is nonsingular when haf− hcd− bgf+ egc− i(ae−bd)= 0⁠.

The remaining cases are routine. Do the a=0 but d≠0 case and the a=0 and d=0 but g≠0 case by first swapping rows and then going on as above. The a=0⁠, d=0⁠, and g=0 case is easy — that matrix is singular since the columns form a linearly dependent set, and the determinant comes out to be zero.

Four.I.1.11 §I.1 source ↗

Show that the equation of a line in ℝ2 through (x1,y1) and (x2,y2) is given by this determinant.

det⁡((xy1x1y11x2y21))=0x1≠x2

Answer

Figuring the determinant and doing some algebra gives this.

0=y1x+x2y+x1y2−y2x−x1y−x2y1(x2−x1)⋅y=(y2−y1)⋅x+x2y1−x1y2y=y2−y1x2−x1⋅x+x2y1−x1y2x2−x1

Note that this is the equation of a line (in particular, in contains the familiar expression for the slope), and note that (x1,y1) and (x2,y2) satisfy it.

Four.I.1.12 §I.1 source ↗

Many people have learned this mnemonic for the determinant of a 3×3 matrix: copy the first two columns to the right side of the matrix, then take the products down the forward diagonals and add them together, and then take the products on the backward diagonals and subtract them. That is, first write

(h1,1h1,2h1,3h1,1h1,2h2,1h2,2h2,3h2,1h2,2h3,1h3,2h3,3h3,1h3,2)

and then calculate this.

h1,1h2,2h3,3+h1,2h2,3h3,1+h1,3h2,1h3,2−h3,1h2,2h1,3−h3,2h2,3h1,1−h3,3h2,1h1,2

(a)

Check that this agrees with the formula given in the preamble to this section.

(b)

Does it extend to other-sized determinants?

Answer (2 parts)

(a)

The comparison with the formula given in the preamble to this section is easy.

(b)

While it holds for 2×2 matrices

(h1,1h1,2h1,1h2,1h2,2h2,1)=h1,1h2,2+h1,2h2,1−h2,1h1,2−h2,2h1,1=h1,1h2,2−h1,2h2,1

it does not hold for 4×4 matrices. An example is that this matrix is singular because the second and third rows are equal

(100101100110−1001)

but following the scheme of the mnemonic does not give zero.

(100110001100110110011−1001−100)=1+0+0+0−(−1)−0−0−0

Four.I.1.13 §I.1 source ↗

The cross product of the vectors

x→=(x1x2x3)y→=(y1y2y3)

is the vector computed as this determinant.

x→×y→=det⁡((e→1e→2e→3x1x2x3y1y2y3))

Note that the first row’s entries are vectors, the vectors from the standard basis for ℝ3⁠. Show that the cross product of two vectors is perpendicular to each vector.

Answer

The determinant is (x2y3−x3y2)e→1+ (x3y1−x1y3)e→2+ (x1y2−x2y1)e→3⁠. To check perpendicularity, we check that the dot product with the first vector is zero

(x1x2x3)∙(x2y3−x3y2x3y1−x1y3x1y2−x2y1)=x1x2y3−x1x3y2+x2x3y1−x1x2y3+x1x3y2−x2x3y1=0

and the dot product with the second vector is also zero.

(y1y2y3)∙(x2y3−x3y2x3y1−x1y3x1y2−x2y1)=x2y1y3−x3y1y2+x3y1y2−x1y2y3+x1y2y3−x2y1y3=0

Four.I.1.14 §I.1 source ↗

Prove that each statement holds for 2×2 matrices.

(a)

The determinant of a product is the product of the determinants det⁡(ST)= det⁡(S)⋅det⁡(T)⁠.

(b)

If T is invertible then the determinant of the inverse is the inverse of the determinant det⁡(T−1)=(det⁡(T))−1⁠.

Matrices T and T′ are similar if there is a nonsingular matrix P such that T′=PTP−1⁠. (We shall look at this relationship in Chapter Five.) Show that similar 2×2 matrices have the same determinant.

Answer (3 parts)

(a)

Plug and chug: the determinant of the product is this

det⁡((abcd)(wxyz))=det⁡((aw+byax+bzcw+dycx+dz))=acwx+adwz+bcxy+bdyz−acwx−bcwz−adxy−bdyz

while the product of the determinants is this.

det⁡((abcd))⋅det⁡((wxyz))=(ad−bc)⋅(wz−xy)

Verification that they are equal is easy.

(b)

Use the prior item.

That similar matrices have the same determinant is immediate from the above two: det⁡(PTP−1)= det⁡(P)⋅det⁡(T)⋅det⁡(P−1)⁠.

Four.I.1.15 §I.1 source ↗

Prove that the area of this region in the plane

is equal to the value of this determinant.

det⁡((x1x2y1y2))

Compare with this.

det⁡((x2x1y2y1))

Answer

One way is to count these areas

by taking the area of the entire rectangle and subtracting the area of A the upper-left rectangle, B the upper-middle triangle, D the upper-right triangle, C the lower-left triangle, E the lower-middle triangle, and F the lower-right rectangle (x1+x2)(y1+y2)− x2y1− (1/2)x1y1− (1/2)x2y2− (1/2)x2y2− (1/2)x1y1− x2y1⁠. Simplification gives the determinant formula.

This determinant is the negative of the one above; the formula distinguishes whether the second column is counterclockwise from the first.

Four.I.1.16 §I.1 source ↗

Prove that for 2×2 matrices, the determinant of a matrix equals the determinant of its transpose. Does that also hold for 3×3 matrices?

Answer

The computation for 2×2 matrices, using the formula quoted in the preamble, is easy. It does also hold for 3×3 matrices; the computation is routine.

Four.I.1.17 §I.1 source ↗

Is the determinant function linear — is det⁡(x⋅T+y⋅S)= x⋅det⁡(T)+ y⋅det⁡(S)⁠?

Answer

No. We illustrate with the 2×2 determinant. Recall that constants come out one row at a time.

det⁡((2426))=2⋅det⁡((1226))=2⋅2⋅det⁡((1213))

This contradicts linearity (here we didn’t need S⁠, i.e., we can take S to be the matrix of zeros).

Four.I.1.18 §I.1 source ↗

Show that if A is 3×3 then det⁡(c⋅A)=c3⋅det⁡(A) for any scalar c⁠.

Answer

Bring out the c⁠’s one row at a time.

Four.I.1.19 §I.1 source ↗

Which real numbers θ make

(cos⁡θ−sin⁡θsin⁡θcos⁡θ)

singular? Explain geometrically.

Answer

There are no real numbers θ that make the matrix singular because the determinant of the matrix cos2⁡θ+sin2⁡θ is never 0⁠, it equals 1 for all θ⁠. Geometrically, with respect to the standard basis, this matrix represents a rotation of the plane through an angle of θ⁠. Each such map is one-to-one — for one thing, it is invertible.

Four.I.1.20 §I.1 source ↗

[Am. Math. Mon., Apr. 1955] If a third order determinant has elements 1⁠, 2⁠, …, 9⁠, what is the maximum value it may have?

Answer

This is how the answer was given in the cited source. Let P be the sum of the three positive terms of the determinant and −N the sum of the three negative terms. The maximum value of P is

9⋅8⋅7+6⋅5⋅4+3⋅2⋅1=630.

The minimum value of N consistent with P is

9⋅6⋅1+8⋅5⋅2+7⋅4⋅3=218.

Any change in P would result in lowering that sum by more than 4⁠. Therefore 412 the maximum value for the determinant and one form for the determinant is

|942386517|.

Four.I.2.9 §I.2 source ↗

Find each determinant by performing one row operation.

(a)

|1−2122−41000−100005|

(b)

|11−200403−6|

Answer (2 parts)

(a)

Do 2ρ1+ρ2 to get echelon form, and then multiply down the diagonal. The determinant is 0⁠.

(b)

Swapping the second and third rows brings the system to echelon form (and changes the sign of the determinant). Multiplying down the diagonal gives 12⁠, so the determinant of the given matrix is −12⁠.

Four.I.2.10 §I.2 source ↗

Use Gauss’s Method to find each determinant.

(a)

|312310014|

(b)

|10012110−10101110|

Answer (2 parts)

(a)

|312310014|= |31200−2014|= − |31201400−2|= 6

(b)

|10012110−10101110|= |1001011−20011011−1|= |1001011−200110001|= 1

Four.I.2.11 §I.2 source ↗

Use Gauss’s Method to find each.

(a)

|2−1−1−1|

(b)

|110302522|

Answer (2 parts)

(a)

|2−1−1−1|=|2−10−3/2|=−3⁠;

(b)

|110302522|= |1100−320−32|= |1100−32000|= 0

Four.I.2.12 §I.2 source ↗

For which values of k does this system have a unique solution?

x+z−w=2y−2z=3x+kz=4z−w=2

Answer

When is the determinant not zero?

|101−101−2010k0001−1|=|101−101−2000k−11001−1|

Obviously, k=1 gives nonsingularity and hence a nonzero determinant. If k≠1 then we get echelon form with a (−1/k−1)ρ3+ρ4 combination.

=|101−101−2000k−11000−1−(1/k−1)|

Multiplying down the diagonal gives (k−1)(−1−(1/k−1))= − (k−1)− 1= − k⁠. Thus the matrix has a nonzero determinant, and so the system has a unique solution, if and only if k≠0⁠.

Four.I.2.13 §I.2 source ↗

Express each of these in terms of |H|⁠.

(a)

|h3,1h3,2h3,3h2,1h2,2h2,3h1,1h1,2h1,3|

(b)

|−h1,1−h1,2−h1,3−2h2,1−2h2,2−2h2,3−3h3,1−3h3,2−3h3,3|

(c)

|h1,1+h3,1h1,2+h3,2h1,3+h3,3h2,1h2,2h2,35h3,15h3,25h3,3|

Answer (3 parts)

(a)

Condition (2) of the definition of determinants applies via the swap ρ1↔ρ3⁠.

|h3,1h3,2h3,3h2,1h2,2h2,3h1,1h1,2h1,3|=−|h1,1h1,2h1,3h2,1h2,2h2,3h3,1h3,2h3,3|

(b)

Condition (3) applies.

|−h1,1−h1,2−h1,3−2h2,1−2h2,2−2h2,3−3h3,1−3h3,2−3h3,3|=(−1)⋅(−2)⋅(−3)⋅|h1,1h1,2h1,3h2,1h2,2h2,3h3,1h3,2h3,3|=(−6)⋅|h1,1h1,2h1,3h2,1h2,2h2,3h3,1h3,2h3,3|

(c)

|h1,1+h3,1h1,2+h3,2h1,3+h3,3h2,1h2,2h2,35h3,15h3,25h3,3|=5⋅|h1,1+h3,1h1,2+h3,2h1,3+h3,3h2,1h2,2h2,3h3,1h3,2h3,3|=5⋅|h1,1h1,2h1,3h2,1h2,2h2,3h3,1h3,2h3,3|

Four.I.2.14 §I.2 source ↗

Find the determinant of a diagonal matrix.

Answer

A diagonal matrix is in echelon form, so the determinant is the product down the diagonal.

Four.I.2.15 §I.2 source ↗

Describe the solution set of a homogeneous linear system if the determinant of the matrix of coefficients is nonzero.

Answer

It is the trivial subspace.

Four.I.2.16 §I.2 source ↗

Show that this determinant is zero.

|y+zx+zx+yxyz111|

Answer

Adding the second row to the first gives a matrix whose first row is x+y+z times its third row.

Four.I.2.17 §I.2 source ↗

(a)

Find the 1×1⁠, 2×2⁠, and 3×3 matrices with i,j entry given by (−1)i+j⁠.

(b)

Find the determinant of the square matrix with i,j entry (−1)i+j⁠.

Answer (2 parts)

(a)

(1)⁠, (1−1−11)⁠, (1−11−11−11−11)

(b)

The determinant in the 1×1 case is 1⁠. In every other case the second row is the negative of the first, and so matrix is singular and the determinant is zero.

Four.I.2.18 §I.2 source ↗

(a)

Find the 1×1⁠, 2×2⁠, and 3×3 matrices with i,j entry given by i+j⁠.

(b)

Find the determinant of the square matrix with i,j entry i+j⁠.

Answer (2 parts)

(a)

(2)⁠, (2334)⁠, (234345456)

(b)

The 1×1 and 2×2 cases yield these.

|2|=2|2334|=−1

And n×n matrices with n≥3 are singular, e.g.,

|234345456|=0

because twice the second row minus the first row equals the third row. Checking this is routine.

Four.I.2.19 §I.2 source ↗

Show that determinant functions are not linear by giving a case where |A+B|≠|A|+|B|⁠.

Answer

This one

A=B=(1234)

is easy to check.

|A+B|=|2468|=−8|A|+|B|=−2−2=−4

By the way, this also gives an example where scalar multiplication is not preserved |2⋅A|≠2⋅|A|⁠.

Four.I.2.20 §I.2 source ↗

The second condition in the definition, that row swaps change the sign of a determinant, is somewhat annoying. It means we have to keep track of the number of swaps, to compute how the sign alternates. Can we get rid of it? Can we replace it with the condition that row swaps leave the determinant unchanged? (If so then we would need new 1×1⁠, 2×2⁠, and 3×3 formulas, but that would be a minor matter.)

Answer

No, we cannot replace it. Remark 2.2 shows that the four conditions after the replacement would conflict — no function satisfies all four.

Four.I.2.21 §I.2 source ↗

Prove that the determinant of any triangular matrix, upper or lower, is the product down its diagonal.

Answer

A upper-triangular matrix is in echelon form.

A lower-triangular matrix is either singular or nonsingular. If it is singular then it has a zero on its diagonal and so its determinant (namely, zero) is indeed the product down its diagonal. If it is nonsingular then it has no zeroes on its diagonal, and we can reduce it by Gauss’s Method to echelon form without changing the diagonal.

Four.I.2.22 §I.2 source ↗

Refer to the definition of elementary matrices in the Mechanics of Matrix Multiplication subsection.

(a)

What is the determinant of each kind of elementary matrix?

(b)

Prove that if E is any elementary matrix then |ES|=|E||S| for any appropriately sized S⁠.

(c)

(This question doesn’t involve determinants.) Prove that if T is singular then a product TS is also singular.

(d)

Show that |TS|=|T||S|⁠.

(e)

Show that if T is nonsingular then |T−1|=|T|−1⁠.

Answer (5 parts)

(a)

The properties in the definition of determinant show that |Mi(k)|=k⁠, |Pi,j|=−1⁠, and |Ci,j(k)|=1⁠.

(b)

The three cases are easy to check by recalling the action of left multiplication by each type of matrix.

(c)

If TS is invertible (TS)M=I then the associative property of matrix multiplication T(SM)=I shows that T is invertible. So if T is not invertible then neither is TS⁠.

(d)

If T is singular then apply the prior answer: |TS|=0 and |T|⋅|S|=0⋅|S|=0⁠. If T is not singular then we can write it as a product of elementary matrices |TS|= |Er⋯E1S|= |Er|⋯|E1|⋅|S|= |Er⋯E1||S|= |T||S|⁠.

(e)

1=|I|=|T⋅T−1|=|T||T−1|

Four.I.2.23 §I.2 source ↗

Prove that the determinant of a product is the product of the determinants |TS|=|T||S| in this way. Fix the n×n matrix S and consider the function d:ℳn×n→ℝ given by T↦|TS|/|S|⁠.

(a)

Check that d satisfies condition (1) in the definition of a determinant function.

(b)

Check condition (2).

(c)

Check condition (3).

(d)

Check condition (4).

(e)

Conclude the determinant of a product is the product of the determinants.

Answer (5 parts)

(a)

We must show that if

T⟶kρi+ρjTˆ

then d(T)= |TS|/|S|= |TˆS|/|S|= d(Tˆ)⁠. We will be done if we show that combining rows first and then multiplying to get TˆS gives the same result as multiplying first to get TS and then combining (because the determinant |TS| is unaffected by the combination so we’ll then have |TˆS|=|TS|⁠, and hence d(Tˆ)=d(T)⁠). That argument runs: after adding k times row i of TS to row j of TS⁠, the j,p entry is (kti,1+tj,1)s1,p+ ⋯+ (kti,r+tj,r)sr,p⁠, which is the j,p entry of TˆS⁠.

(b)

We need only show that swapping T⟶ρi↔ρjTˆ and then multiplying to get TˆS gives the same result as multiplying T by S and then swapping (because, as the determinant |TS| changes sign on the row swap, we’ll then have |TˆS|=−|TS|⁠, and so d(Tˆ)=−d(T)⁠). That argument runs just like the prior one.

(c)

Not surprisingly by now, we need only show that multiplying a row by a scalar T⟶kρiTˆ and then computing TˆS gives the same result as first computing TS and then multiplying the row by k (as the determinant |TS| is rescaled by k the multiplication, we’ll have |TˆS|=k|TS|⁠, so d(Tˆ)=kd(T)⁠). The argument runs just as above.

(d)

Clear.

(e)

Because we’ve shown that d(T) is a determinant and that determinant functions (if they exist) are unique, we have that so |T|=d(T)=|TS|/|S|⁠.

Four.I.2.24 §I.2 source ↗

A submatrix of a given matrix A is one that we get by deleting some of the rows and columns of A⁠. Thus, the first matrix here is a submatrix of the second.

(3125)(34109−22−15)

Prove that for any square matrix, the rank of the matrix is r if and only if r is the largest integer such that there is an r×r submatrix with a nonzero determinant.

Answer

We will first argue that a rank r matrix has a r×r submatrix with nonzero determinant. A rank r matrix has a linearly independent set of r rows. A matrix made from those rows will have row rank r and thus has column rank r⁠. Conclusion: from those r rows we can extract a linearly independent set of r columns, and so the original matrix has a r×r submatrix of rank r⁠.

We finish by showing that if r is the largest such integer then the rank of the matrix is r⁠. We need only show, by the maximality of r⁠, that if a matrix has a k×k submatrix of nonzero determinant then the rank of the matrix is at least k⁠. Consider such a k×k submatrix. Its rows are parts of the rows of the original matrix, clearly the set of whole rows is linearly independent. Thus the row rank of the original matrix is at least k⁠, and the row rank of a matrix equals its rank.

Four.I.2.25 §I.2 source ↗

Prove that a matrix with rational entries has a rational determinant.

Answer

A matrix with only rational entries reduces with Gauss’s Method to an echelon form matrix using only rational arithmetic. Thus the entries on the diagonal must be rationals, and so the product down the diagonal is rational.

Four.I.2.26 §I.2 source ↗

[Am. Math. Mon., Feb. 1953] Find the element of likeness in (a) simplifying a fraction, (b) powdering the nose, (c) building new steps on the church, (d) keeping emeritus professors on campus, (e) putting B⁠, C⁠, D in the determinant

|1aa2a3a31aa2Ba31aCDa31|.

Answer

This is how the answer was given in the cited source. The value (1−a4)3 of the determinant is independent of the values B⁠, C⁠, D⁠. Hence operation (e) does not change the value of the determinant but merely changes its appearance. Thus the element of likeness in (a), (b), (c), (d), and (e) is only that the appearance of the principle entity is changed. The same element appears in (f) changing the name-label of a rose, (g) writing a decimal integer in the scale of 12⁠, (h) gilding the lily, (i) whitewashing a politician, and (j) granting an honorary degree.

Four.I.3.18 §I.3 source ↗

For this matrix, find the term associated with each 3⁠-permutation.

M=(123456789)

That is, fill in the rest of this table.

permutation ϕiϕ1ϕ2ϕ3ϕ4ϕ5ϕ6
term m1,ϕi(1)m2,ϕi(2)m3,ϕi(3)1⋅5⋅9

Answer

Call the matrix M⁠.

permutationϕ1ϕ2ϕ3ϕ4ϕ5ϕ6
term1⋅5⋅91⋅6⋅82⋅4⋅92⋅6⋅73⋅4⋅83⋅5⋅7

Four.I.3.19 §I.3 source ↗

For each 3⁠-permutation ϕ find |Pϕ|⁠.

Answer

We can swap each Pϕ to the identity matrix. If the number of swaps is even then it’s determinant is +1⁠, while if the number of swaps is odd then the determinant is −1

permutationϕ1ϕ2ϕ3ϕ4ϕ5ϕ6
number of swaps011221
|Pϕi|+1−1−1+1+1−1

Remark. In that table we simply swapped until we found a number that brought us back to ‘⁠1,2,3⁠’. But is a different number of swaps possible? If one person found 2 swaps and another found 4 that would be OK, since both give a determinant of +1⁠. But if we can swap in two different ways and one of them is even and one is odd then that would be a problem. The next section shows that a mix of even and odd is not possible.

Four.I.3.20 §I.3 source ↗

This determinant is 7 by the 2×2 formula. Compute it with the permutation expansion.

|2315|

Answer

|2315|=2⋅5⋅|Pϕ1|+1⋅3⋅|Pϕ2|=2⋅5⋅|1001|+1⋅3⋅|0110|=2⋅5⋅1+1⋅3⋅(−1)=7

Four.I.3.21 §I.3 source ↗

This determinant is 0 because the first two rows add to the third. Compute the determinant using the permutation expansion.

|−101314215|

Answer

|−101314215|=[t](−1)(1)(5)|Pϕ1|+(−1)(4)(1)|Pϕ2|+(0)(3)(5)|Pϕ3|+(0)(4)(2)|Pϕ4|+(1)(3)(1)|Pϕ5|+(1)(1)(2)|Pϕ6|=[t](−1)(1)(5)⋅1+(−1)(4)(1)⋅(−1)+(0)(3)(5)⋅(−1)+(0)(4)(2)⋅1+(1)(3)(1)⋅1+(1)(1)(2)⋅(−1)=−5+4+0+0+3−2=0

Four.I.3.22 §I.3 source ↗

Compute the determinant by using the permutation expansion.

(a)

|123456789|

(b)

|2213−10−205|

Answer (2 parts)

(a)

This matrix is singular.

|123456789|=[t](1)(5)(9)|Pϕ1|+(1)(6)(8)|Pϕ2|+(2)(4)(9)|Pϕ3|+(2)(6)(7)|Pϕ4|+(3)(4)(8)|Pϕ5|+(7)(5)(3)|Pϕ6|=0

(b)

This matrix is nonsingular.

|2213−10−205|=[t](2)(−1)(5)|Pϕ1|+(2)(0)(0)|Pϕ2|+(2)(3)(5)|Pϕ3|+(2)(0)(−2)|Pϕ4|+(1)(3)(0)|Pϕ5|+(−2)(−1)(1)|Pϕ6|=−42

Four.I.3.23 §I.3 source ↗

Compute these both with Gauss’s Method and the permutation expansion formula.

(a)

|2131|

(b)

|014023151|

Answer (2 parts)

(a)

Gauss’s Method gives this

|2131|=|210−1/2|=−1

and permutation expansion gives this.

|2131|=|2001|+|0130|=(2)(1)|1001|+(1)(3)|0110|=−1

(b)

Gauss’s Method gives this

|014023151|=−|151023014|=−|151023005/2|=−5

and the permutation expansion gives this.

|014023151|=[t](0)(2)(1)|Pϕ1|+(0)(3)(5)|Pϕ2|+(1)(0)(1)|Pϕ3|+(1)(3)(1)|Pϕ4|+(4)(0)(5)|Pϕ5|+(4)(2)(1)|Pϕ6|=−5

Four.I.3.24 §I.3 source ↗

Use the permutation expansion formula to derive the formula for 3×3 determinants.

Answer

Following Example 3.6 gives this.

|t1,1t1,2t1,3t2,1t2,2t2,3t3,1t3,2t3,3|=[t]t1,1t2,2t3,3|Pϕ1|+t1,1t2,3t3,2|Pϕ2|+t1,2t2,1t3,3|Pϕ3|+t1,2t2,3t3,1|Pϕ4|+t1,3t2,1t3,2|Pϕ5|+t1,3t2,2t3,1|Pϕ6|=[t]t1,1t2,2t3,3(+1)+t1,1t2,3t3,2(−1)+t1,2t2,1t3,3(−1)+t1,2t2,3t3,1(+1)+t1,3t2,1t3,2(+1)+t1,3t2,2t3,1(−1)

Four.I.3.25 §I.3 source ↗

List all of the 4⁠-permutations.

Answer

This is all of the permutations where ϕ(1)=1

ϕ1=⟨1,2,3,4⟩ϕ2=⟨1,2,4,3⟩ϕ3=⟨1,3,2,4⟩ϕ4=⟨1,3,4,2⟩ϕ5=⟨1,4,2,3⟩ϕ6=⟨1,4,3,2⟩

the ones where ϕ(1)=1

ϕ7=⟨2,1,3,4⟩ϕ8=⟨2,1,4,3⟩ϕ9=⟨2,3,1,4⟩ϕ10=⟨2,3,4,1⟩ϕ11=⟨2,4,1,3⟩ϕ12=⟨2,4,3,1⟩

the ones where ϕ(1)=3

ϕ13=⟨3,1,2,4⟩ϕ14=⟨3,1,4,2⟩ϕ15=⟨3,2,1,4⟩ϕ16=⟨3,2,4,1⟩ϕ17=⟨3,4,1,2⟩ϕ18=⟨3,4,2,1⟩

and the ones where ϕ(1)=4⁠.

ϕ19=⟨4,1,2,3⟩ϕ20=⟨4,1,3,2⟩ϕ21=⟨4,2,1,3⟩ϕ22=⟨4,2,3,1⟩ϕ23=⟨4,3,1,2⟩ϕ24=⟨4,3,2,1⟩

Four.I.3.26 §I.3 source ↗

A permutation, regarded as a function from the set {1,..,n} to itself, is one-to-one and onto. Therefore, each permutation has an inverse.

(a)

Find the inverse of each 2⁠-permutation.

(b)

Find the inverse of each 3⁠-permutation.

Answer (3 parts)

Each of these is easy to check.

(a)
permutationϕ1ϕ2
inverseϕ1ϕ2
(b)
permutationϕ1ϕ2ϕ3ϕ4ϕ5ϕ6
inverseϕ1ϕ2ϕ3ϕ5ϕ4ϕ6

Four.I.3.27 §I.3 source ↗

Prove that f is multilinear if and only if for all v→,w→∈V and k1,k2∈ℝ⁠, this holds.

f(ρ→1,…,k1v→1+k2v→2,…,ρ→n)=k1f(ρ→1,…,v→1,…,ρ→n)+k2f(ρ→1,…,v→2,…,ρ→n)

Answer

For the ‘if’ half, the first condition of Definition 3.2 follows from taking k1=k2=1 and the second condition follows from taking k2=0⁠.

The ‘only if’ half also routine. From f(ρ→1, …, k1v→1+ k2v→2, …, ρ→n) the first condition of Definition 3.2 gives = f(ρ→1,…,k1v→1,…,ρ→n)+ f(ρ→1,…,k2v→2,…,ρ→n) and the second condition, applied twice, gives the result.

Four.I.3.28 §I.3 source ↗

How would determinants change if we changed property (4) of the definition to read that |I|=2⁠?

Answer

They would all double.

Four.I.3.29 §I.3 source ↗

Verify the second and third statements in Corollary 3.16.

Answer

For the second statement, given a matrix, transpose it, swap rows, and transpose back. The result is swapped columns, and the determinant changes by a factor of −1⁠. The third statement is similar: given a matrix, transpose it, apply multilinearity to what are now rows, and then transpose back the resulting matrices.

Four.I.3.30 §I.3 source ↗

Show that if an n×n matrix has a nonzero determinant then we can express any column vector v→∈ℝn as a linear combination of the columns of the matrix.

Answer

An n×n matrix with a nonzero determinant has rank n so its columns form a basis for ℝn⁠.

Four.I.3.31 §I.3 source ↗

[Strang 80] True or false: a matrix whose entries are only zeros or ones has a determinant equal to zero, one, or negative one.

Answer

False.

|011101110|=2

Four.I.3.32 §I.3 source ↗

(a)

Show that there are 120 terms in the permutation expansion formula of a 5×5 matrix.

(b)

How many are sure to be zero if the 1,2 entry is zero?

Answer (2 parts)

(a)

For the column index of the entry in the first row there are five choices. Then, for the column index of the entry in the second row there are four choices. Continuing, we get 5⋅4⋅3⋅2⋅1=120⁠. (See also the next question.)

(b)

Once we choose the second column in the first row, we can choose the other entries in 4⋅3⋅2⋅1=24 ways.

Four.I.3.33 §I.3 source ↗

How many n⁠-permutations are there?

Answer

n⋅(n−1)⋯2⋅1=n!

Four.I.3.34 §I.3 source ↗

Show that the inverse of a permutation matrix is its transpose.

Answer

[Schmidt] We will show that PP𝖳=I⁠; the P𝖳P=I argument is similar. The i,j entry of PP𝖳 is the sum of terms of the form pi,kqk,j where the entries of P𝖳 are denoted with q⁠’s, that is, qk,j=pj,k⁠. Thus the i,j entry of PP𝖳 is the sum ∑k=1npi,kpj,k⁠. But pi,k is usually 0⁠, and so Pi,kPj,k is usually 0⁠. The only time Pi,k is nonzero is when it is 1⁠, but then there are no other i′≠i such that Pi′,k is nonzero (i is the only row with a 1 in column k⁠). In other words,

∑k=1npi,kpj,k={1i=j0otherwise

and this is exactly the formula for the entries of the identity matrix.

Four.I.3.35 §I.3 source ↗

A matrix A is skew-symmetric if A𝖳=−A⁠, as in this matrix.

A=(03−30)

Show that n×n skew-symmetric matrices with nonzero determinants exist only for even n⁠.

Answer

In |A|=|A𝖳|=|−A|=(−1)n|A| the exponent n must be even.

Four.I.3.36 §I.3 source ↗

What is the smallest number of zeros, and the placement of those zeros, needed to ensure that a 4×4 matrix has a determinant of zero?

Answer

Showing that no placement of three zeros suffices is routine. Four zeroes does suffice; put them all in the same row or column.

Four.I.3.37 §I.3 source ↗

If we have n data points (x1,y1), (x2,y2), …, (xn,yn) and want to find a polynomial p(x)= an−1xn−1+ an−2xn−2+ ⋯+ a1x+ a0 passing through those points then we can plug in the points to get an n⁠ equation/⁠n⁠ unknown linear system. The matrix of coefficients for that system is the Vandermonde matrix. Prove that the determinant of the transpose of that matrix of coefficients

|11…1x1x2…xnx12x22…xn2⋮x1n−1x2n−1…xnn−1|

equals the product, over all indices i,j∈{1,…,n} with i<j⁠, of terms of the form xj−xi⁠. (This shows that the determinant is zero, and the linear system has no solution, if and only if the xi⁠’s in the data are not distinct.)

Answer

The n=3 case shows what to do. The row combination operations of −x1ρ2+ρ3 and −x1ρ1+ρ2 give this.

|111x1x2x3x12x22x32|=|111x1x2x30(−x1+x2)x2(−x1+x3)x3|=|1110−x1+x2−x1+x30(−x1+x2)x2(−x1+x3)x3|

Then the row combination operation of x2ρ2+ρ3 gives the desired result.

=|1110−x1+x2−x1+x300(−x1+x3)(−x2+x3)|=(x2−x1)(x3−x1)(x3−x2)

Four.I.3.38 §I.3 source ↗

We can divide a matrix into blocks, as here,

(12034000−2)

which shows four blocks, the square 2×2 and 1×1 ones in the upper left and lower right, and the zero blocks in the upper right and lower left. Show that if a matrix is such that we can partition it as

T=(JZ2Z1K)

where J and K are square, and Z1 and Z2 are all zeroes, then |T|=|J|⋅|K|⁠.

Answer

Let T be n×n⁠, let J be p×p⁠, and let K be q×q⁠. Apply the permutation expansion formula

|T|=∑permutations ϕt1,ϕ(1)t2,ϕ(2)…tn,ϕ(n)|Pϕ|

Because the upper right of T is all zeroes, if a ϕ has at least one of p+1,…,n among its first p column numbers ϕ(1),…,ϕ(p) then the term arising from ϕ is 0 (e.g., if ϕ(1)=n then t1,ϕ(1)t2,ϕ(2)…tn,ϕ(n) is 0⁠). So the above formula reduces to a sum over all permutations with two halves: first rearrange 1,…,p and after that comes a permutation of p+1,…,p+q⁠. To see this gives |J|⋅|K|⁠, distribute.

[∑perms ϕ1of 1,…,pt1,ϕ1(1)⋯tp,ϕ1(p)|Pϕ1|]⋅[∑perms ϕ2of p+1,…,p+qtp+1,ϕ2(p+1)⋯tp+q,ϕ2(p+q)|Pϕ2|]

Four.I.3.39 §I.3 source ↗

Prove that for any n×n matrix T there are at most n distinct reals r such that the matrix T−rI has determinant zero (we shall use this result in Chapter Five).

Answer

The n=3 case shows what happens.

|T−rI|=|t1,1−xt1,2t1,3t2,1t2,2−xt2,3t3,1t3,2t3,3−x|

Each term in the permutation expansion has three factors drawn from entries in the matrix (e.g., (t1,1− x)(t2,2− x)(t3,3− x) and (t1,1−x)(t2,3)(t3,2)⁠), and so the determinant is expressible as a polynomial in x of degree 3⁠. Such a polynomial has at most 3 roots.

In general, the permutation expansion shows that the determinant is a sum of terms, each with n factors, giving a polynomial of degree n⁠. A polynomial of degree n has at most n roots.

Four.I.3.40 §I.3 source ↗

[Math. Mag., Jan. 1963, Q307] The nine positive digits can be arranged into 3×3 arrays in 9! ways. Find the sum of the determinants of these arrays.

Answer

This is how the answer was given in the cited source. When two rows of a determinant are interchanged, the sign of the determinant is changed. When the rows of a three-by-three determinant are permuted, 3 positive and 3 negative determinants equal in absolute value are obtained. Hence the 9! determinants fall into 9!/6 groups, each of which sums to zero.

Four.I.3.41 §I.3 source ↗

[Math. Mag., Jan. 1963, Q237] Show that

|x−2x−3x−4x+1x−1x−3x−4x−7x−10|=0.

Answer

This is how the answer was given in the cited source. When the elements of any column are subtracted from the elements of each of the other two, the elements in two of the columns of the derived determinant are proportional, so the determinant vanishes. That is,

|21x−442x−363x−10|=|1x−3−12x−1−23x−7−3|=|x−2−1−2x+1−2−4x−4−3−6|=0.

Four.I.3.42 §I.3 source ↗

[Am. Math. Mon., Jan. 1949] Let S be the sum of the integer elements of a magic square of order three and let D be the value of the square considered as a determinant. Show that D/S is an integer.

Answer

This is how the answer was given in the cited source. Let

abcdefghi

have magic sum N=S/3⁠. Then

N=(a+e+i)+(d+e+f)+(g+e+c)−(a+d+g)−(c+f+i)=3e

and S=9e⁠. Hence, adding rows and columns,

D=|abcdefghi|=|abcdef3e3e3e|=|ab3ede3e3e3e9e|=|abedee111|S.

Four.I.3.43 §I.3 source ↗

[Am. Math. Mon., Jun. 1931] Show that the determinant of the n2 elements in the upper left corner of the Pascal triangle

1111..123..13..1....

has the value unity.

Answer

This is how the answer was given in the cited source. Denote by Dn the determinant in question and by ai,j the element in the i⁠-th row and j⁠-th column. Then from the law of formation of the elements we have

ai,j=ai,j−1+ai−1,j,a1,j=ai,1=1.

Subtract each row of Dn from the row following it, beginning the process with the last pair of rows. After the n−1 subtractions the above equality shows that the element ai,j is replaced by the element ai,j−1⁠, and all the elements in the first column, except a1,1=1⁠, become zeroes. Now subtract each column from the one following it, beginning with the last pair. After this process the element ai,j−1 is replaced by ai−1,j−1⁠, as shown in the above relation. The result of the two operations is to replace ai,j by ai−1,j−1⁠, and to reduce each element in the first row and in the first column to zero. Hence Dn=Dn+i and consequently

Dn=Dn−1=Dn−2=⋯=D2=1.

Four.I.4.12 §I.4 source ↗

Give the permutation expansion of a general 2×2 matrix and its transpose.

Answer

This is the permutation expansion of the determinant of a 2×2 matrix

|abcd|=ad⋅|1001|+bc⋅|0110|

and the permutation expansion of the determinant of its transpose.

|acbd|=ad⋅|1001|+cb⋅|0110|

As with the 3×3 expansions described in the subsection, the permutation matrices from corresponding terms are transposes (although this is disguised by the fact that each is self-transpose).

Four.I.4.13 §I.4 source ↗

This problem appears also in the prior subsection.

(a)

Find the inverse of each 2⁠-permutation.

(b)

Find the inverse of each 3⁠-permutation.

Answer (3 parts)

Each of these is easy to check.

(a)
permutationϕ1ϕ2
inverseϕ1ϕ2
(b)
permutationϕ1ϕ2ϕ3ϕ4ϕ5ϕ6
inverseϕ1ϕ2ϕ3ϕ5ϕ4ϕ6

Four.I.4.14 §I.4 source ↗

(a)

Find the signum of each 2⁠-permutation.

(b)

Find the signum of each 3⁠-permutation.

Answer (2 parts)

(a)

sgn⁡(ϕ1)=+1⁠, sgn⁡(ϕ2)=−1

(b)

sgn⁡(ϕ1)=+1⁠, sgn⁡(ϕ2)=−1⁠, sgn⁡(ϕ3)=−1⁠, sgn⁡(ϕ4)=+1⁠, sgn⁡(ϕ5)=+1⁠, sgn⁡(ϕ6)=−1

Four.I.4.15 §I.4 source ↗

Find the only nonzero term in the permutation expansion of this matrix.

|0100101001010010|

Compute that determinant by finding the signum of the associated permutation.

Answer

To get a nonzero term in the permutation expansion we must use the 1,2 entry and the 4,3 entry. Having fixed on those two we must also use the 2,1 entry and the 3,4 entry. The signum of ⟨2,1,4,3⟩ is +1 because from

(0100100000010010)

the two row swaps ρ1↔ρ2 and ρ3↔ρ4 will produce the identity matrix.

Four.I.4.16 §I.4 source ↗

[Strang 80] What is the signum of the n⁠-permutation ϕ=⟨n,n−1,…,2,1⟩⁠?

Answer

The pattern is this.

i123456…
sgn⁡(ϕi)+1−1−1+1+1−1…

So to find the signum of ϕn!⁠, we subtract one n!−1 and look at the remainder on division by four. If the remainder is 1 or 2 then the signum is −1⁠, otherwise it is +1⁠. For n>4⁠, the number n! is divisible by four, so n!−1 leaves a remainder of −1 on division by four (more properly said, a remainder or 3⁠), and so the signum is +1⁠. The n=1 case has a signum of +1⁠, the n=2 case has a signum of −1 and the n=3 case has a signum of −1⁠.

Four.I.4.17 §I.4 source ↗

Prove these.

(a)

Every permutation has an inverse.

(b)

sgn⁡(ϕ−1)=sgn⁡(ϕ)

(c)

Every permutation is the inverse of another.

Answer (3 parts)

(a)

We can view permutations as maps ϕ:{1,…,n}→{1,…,n} that are one-to-one and onto. Any one-one and onto map has an inverse.

(b)

If it always takes an odd number of swaps to get from Pϕ to the identity, then it always takes an odd number of swaps to get from the identity to Pϕ (any swap is reversible).

(c)

This is the first question again.

Four.I.4.18 §I.4 source ↗

Prove that the matrix of the permutation inverse is the transpose of the matrix of the permutation Pϕ−1=Pϕ𝖳⁠, for any permutation ϕ⁠.

Answer

If ϕ(i)=j then ϕ−1(j)=i⁠. The result now follows on the observation that Pϕ has a 1 in entry i,j if and only if ϕ(i)=j⁠, and Pϕ−1 has a 1 in entry j,i if and only if ϕ−1(j)=i⁠,

Four.I.4.19 §I.4 source ↗

Show that a permutation matrix with m inversions can be row swapped to the identity in m steps. Contrast this with Corollary 4.5.

Answer

This does not say that m is the least number of swaps to produce an identity, nor does it say that m is the most. It instead says that there is a way to swap to the identity in exactly m steps.

Let ιj be the first row that is inverted with respect to a prior row and let ιk be the first row giving that inversion. We have this interval of rows.

(⋮ιkιr1⋮ιrsιj⋮)j<k<r1<⋯<rs

Swap.

(⋮ιjιr1⋮ιrsιk⋮)

The second matrix has one fewer inversion because there is one fewer inversion in the interval (s vs. s+1⁠) and inversions involving rows outside the interval are not affected.

Proceed in this way, at each step reducing the number of inversions by one with each row swap. When no inversions remain the result is the identity.

The contrast with Corollary 4.5 is that the statement of this exercise is a ‘there exists’ statement: there exists a way to swap to the identity in exactly m steps. But the corollary is a ‘for all’ statement: for all ways to swap to the identity, the parity (evenness or oddness) is the same.

Four.I.4.20 §I.4 source ↗

For any permutation ϕ let g(ϕ) be the integer defined in this way.

g(ϕ)=∏i<j[ϕ(j)−ϕ(i)]

(This is the product, over all indices i and j with i<j⁠, of terms of the given form.)

(a)

Compute the value of g on all 2⁠-permutations.

(b)

Compute the value of g on all 3⁠-permutations.

(c)

Prove that g(ϕ) is not 0⁠.

(d)

Prove this.

sgn⁡(ϕ)=g(ϕ)|g(ϕ)|

Many authors give this formula as the definition of the signum function.

Answer (4 parts)

(a)

First, g(ϕ1) is the product of the single factor 2−1 and so g(ϕ1)=1⁠. Second, g(ϕ2) is the product of the single factor 1−2 and so g(ϕ2)=−1⁠.

(b)
permutation ϕϕ1ϕ2ϕ3ϕ4ϕ5ϕ6
g(ϕ)2−2−222−2
(c)

It is a product of nonzero terms.

(d)

Note that ϕ(j)−ϕ(i) is negative if and only if ιϕ(j) and ιϕ(i) are in an inversion of their usual order.

Four.II.1.8 §II.1 source ↗

Is

(412)

inside of the box formed by these three?

(331)(261)(105)

Answer

Solving

c1(331)+c2(261)+c3(105)=(412)

gives the unique solution c3=11/57⁠, c2=−40/57 and c1=99/57⁠. Because c1>1⁠, the vector is not in the box.

Four.II.1.9 §II.1 source ↗

Find the volume of the region defined by the vectors.

(a)

⟨(13),(−14)⟩

(b)

⟨(210),(3−24),(8−38)⟩

(c)

⟨(1201), (2222), (−1305), (0107)⟩

Answer (4 parts)

For each, find the determinant and take the absolute value.

(a)

7

(b)

0

(c)

−58

Four.II.1.10 §II.1 source ↗

In this picture the rectangle on the left is defined by the points (3,1) and (1,2)⁠. Apply the matrix to get the rectangle on the right, defined by (7,1) and (4,2)⁠. Why doesn’t this contradict Theorem 1.5?

⟶(2101)
area is 2determinant is 2area is 3

Answer

We have drawn that picture to mislead. The picture on the left is not the box formed by two vectors. If we slide it to the origin then it becomes the box formed by this sequence.

⟨(01),(20)⟩

Then the image under the action of the matrix is the box formed by this sequence.

⟨(11),(40)⟩

which has an area of 4⁠.

Four.II.1.11 §II.1 source ↗

Find the volume of this region.

Answer

Move the parallelepiped to start at the origin, so that it becomes the box formed by

⟨(30),(21)⟩

and now the absolute value of this determinant is easily computed as 3⁠.

|3201|=3

Four.II.1.12 §II.1 source ↗

Suppose that |A|=3⁠. By what factor do these change volumes?

(a)

A

(b)

A2

(c)

A−2

Answer (3 parts)

(a)

3

(b)

9

(c)

1/9

Four.II.1.13 §II.1 source ↗

Consider the linear transformation t of ℝ3 represented with respect to the standard bases by this matrix.

(10−1311−103)

(a)

Compute the determinant of the matrix. Does the transformation preserve orientation or reverse it?

(b)

Find the size of the box defined by these vectors. What is its orientation?

(1−12)(20−1)(110)

(c)

Find the images under t of the vectors in the prior item and find the size of the box that they define. What is the orientation?

Answer (3 parts)

(a)

Gauss’s Method

⟶ρ1+ρ3−3ρ1+ρ2(10−1014002)

gives the determinant as +2⁠. The sign is positive so the transformation preserves orientation.

(b)

The size of the box is the value of this determinant.

|121−1012−10|=+6

The orientation is positive.

(c)

Since this transformation is represented by the given matrix with respect to the standard bases, and with respect to the standard basis the vectors represent themselves, to find the image of the vectors under the transformation just multiply them, from the left, by the matrix.

(1−12)↦(−145)(20−1)↦(35−5)(110)↦(14−1)

Then compute the size of the resulting box.

|−1314545−5−1|=+12

The starting box is positively oriented, the transformation preserves orientations (since the determinant of the matrix is positive), and the ending box is also positively oriented.

Four.II.1.14 §II.1 source ↗

By what factor does each transformation change the size of boxes?

(a)

(xy)↦(2x3y)

(b)

(xy)↦(3x−y−2x+y)

(c)

(xyz)↦(x−yx+y+zy−2z)

Answer (4 parts)

Express each transformation with respect to the standard bases and find the determinant.

(a)

6

(b)

−1

(c)

−5

Four.II.1.15 §II.1 source ↗

What is the area of the image of the rectangle [2..4]×[2..5] under the action of this matrix?

(234−1)

Answer

The starting area is 6 and the matrix changes sizes by −14⁠. Thus the area of the image is 84⁠.

Four.II.1.16 §II.1 source ↗

If t:ℝ3→ℝ3 changes volumes by a factor of 7 and s:ℝ3→ℝ3 changes volumes by a factor of 3/2 then by what factor will their composition changes volumes?

Answer

By a factor of 21/2⁠.

Four.II.1.17 §II.1 source ↗

In what way does the definition of a box differ from the definition of a span?

Answer

For a box we take a sequence of vectors (as described in the remark, the order of the vectors matters), while for a span we take a set of vectors. Also, for a box subset of ℝn there must be n vectors; of course for a span there can be any number of vectors. Finally, for a box the coefficients t1⁠, …, tn are in the interval [0..1]⁠, while for a span the coefficients are free to range over all of ℝ⁠.

Four.II.1.18 §II.1 source ↗

Does |TS|=|ST|⁠? |T(SP)|=|(TS)P|⁠?

Answer

Yes to both. For instance, the first is |TS|= |T|⋅|S|= |S|⋅|T|= |ST|⁠.

Four.II.1.19 §II.1 source ↗

Show that there are no 2×2 matrices A and B satisfying these.

AB=(1−120)BA=(2111)

Answer

Because |AB|=|A|⋅|B|=|BA| and these two matrices have different determinants.

Four.II.1.20 §II.1 source ↗

(a)

Suppose that |A|=3 and that |B|=2⁠. Find |A2⋅B𝖳⋅B−2⋅A𝖳|⁠.

(b)

Assume that |A|=0⁠. Prove that |6A3+5A2+2A|=0⁠.

Answer (2 parts)

(a)

If it is defined then it is (32)⋅(2)⋅(2−2)⋅(3)⁠.

(b)

|6A3+ 5A2+ 2A|= |A|⋅|6A2+ 5A+ 2I|⁠.

Four.II.1.21 §II.1 source ↗

Let T be the matrix representing (with respect to the standard bases) the map that rotates plane vectors counterclockwise through θ radians. By what factor does T change sizes?

Answer

|cos⁡θ−sin⁡θsin⁡θcos⁡θ|= 1

Four.II.1.22 §II.1 source ↗

Must a transformation t:ℝ2→ℝ2 that preserves areas also preserve lengths?

Answer

No, for instance the determinant of

T=(2001/2)

is 1 so it preserves areas, but the vector Te→1 has length 2⁠.

Four.II.1.23 §II.1 source ↗

What is the volume of a parallelepiped in ℝ3 bounded by a linearly dependent set?

Answer

It is zero.

Four.II.1.24 §II.1 source ↗

Find the area of the triangle in ℝ3 with endpoints (1,2,1)⁠, (3,−1,4)⁠, and (2,2,2)⁠. (This asks for area, not volume. The triangle defines a plane; what is the area of the triangle in that plane?)

Answer

Two of the three sides of the triangle are formed by these vectors.

(222)−(121)=(101)(3−14)−(121)=(2−33)

One way to find the area of this triangle is to find half the volume of the parallelogram formed by those two vectors and a length-one vector orthogonal to those two. To find the family of vectors orthogonal to those two we can start with the two relations

(101)∙(xyz)=(000)(2−33)∙(xyz)=(000)

and solve the system

x+z=02x−3y+3z=0⟶−2ρ1+ρ2x+z=0−3y+z=0

to get this solution set.

{(−11/31)z|z∈ℝ}

Here is a length one solution.

119/9⋅(−11/31)=(−3/191/193/19)

Thus the area of the triangle is half of the absolute value of this determinant.

|12−3/190−31/19133/19|=−19/19

Half of the absolute value is 19/2⁠.

Four.II.1.25 §II.1 source ↗

An alternate proof of Theorem 1.5 uses the definition of determinant functions.

(a)

Note that the vectors forming S make a linearly dependent set if and only if |S|=0⁠, and check that the result holds in this case.

(b)

For the |S|≠0 case, to show that |TS|/|S|=|T| for all transformations, consider the function d:ℳn×n→ℝ given by T↦|TS|/|S|⁠. Show that d has the first property of a determinant.

(c)

Show that d has the remaining three properties of a determinant function.

(d)

Conclude that |TS|=|T|⋅|S|⁠.

Answer (4 parts)

(a)

Because the image of a linearly dependent set is linearly dependent, if the vectors forming S make a linearly dependent set, so that |S|=0⁠, then the vectors forming t(S) make a linearly dependent set, so that |TS|=0⁠, and in this case the equation holds.

(b)

We must check that if T⟶kρi+ρjTˆ then d(T)= |TS|/|S|= |TˆS|/|S|= d(Tˆ)⁠. We can do this by checking that combining rows first and then multiplying to get TˆS gives the same result as multiplying first to get TS and then combining (because the determinant |TS| is unaffected by the combining rows so we’ll then have that |TˆS|=|TS| and hence that d(Tˆ)=d(T)⁠). This check runs: after adding k times row i of TS to row j of TS⁠, the j,p entry is (kti,1+tj,1)s1,p+ ⋯+ (kti,r+tj,r)sr,p⁠, which is the j,p entry of TˆS⁠.

(c)

For the second property, we need only check that swapping T⟶ρi↔ρjTˆ and then multiplying to get TˆS gives the same result as multiplying T by S first and then swapping (because, as the determinant |TS| changes sign on the row swap, we’ll then have |TˆS|=−|TS|⁠, and so d(Tˆ)=−d(T)⁠). This check runs just like the one for the first property.

For the third property, we need only show that performing T⟶kρiTˆ and then computing TˆS gives the same result as first computing TS and then performing the scalar multiplication (as the determinant |TS| is rescaled by k⁠, we’ll have |TˆS|=k|TS| and so d(Tˆ)=kd(T)⁠). Here too, the argument runs just as above.

The fourth property, that if T is I then the result is 1⁠, is obvious.

(d)

Determinant functions are unique, so |TS|/|S|=d(T)=|T|⁠, and so |TS|=|T||S|⁠.

Four.II.1.26 §II.1 source ↗

Give a non-identity matrix with the property that A𝖳=A−1⁠. Show that if A𝖳=A−1 then |A|=±1⁠. Does the converse hold?

Answer

Any permutation matrix has the property that the transpose of the matrix is its inverse.

For the implication, we know that |A𝖳|=|A|⁠. Then 1= |A⋅A−1|= |A⋅A𝖳|= |A|⋅|A𝖳|= |A|2⁠.

The converse does not hold; here is an example.

(3121)

Four.II.1.27 §II.1 source ↗

The algebraic property of determinants that factoring a scalar out of a single row will multiply the determinant by that scalar shows that where H is 3×3⁠, the determinant of cH is c3 times the determinant of H⁠. Explain this geometrically, that is, using Theorem 1.5. (The observation that increasing the linear size of a three-dimensional object by a factor of c will increase its volume by a factor of c3 while only increasing its surface area by an amount proportional to a factor of c2 is the Square-cube law [Wikipedia, Square-cube Law].)

Answer

Where the sides of the box are c times longer, the box has c3 times as many cubic units of volume.

Four.II.1.28 §II.1 source ↗

We say that matrices H and G are similar if there is a nonsingular matrix P such that H=P−1GP (we will study this relation in Chapter Five). Show that similar matrices have the same determinant.

Answer

If H=P−1GP then |H|= |P−1||G||P|= |P−1||P||G|= |P−1P||G|= |G|⁠.

Four.II.1.29 §II.1 source ↗

We usually represent vectors in ℝ2 with respect to the standard basis so vectors in the first quadrant have both coordinates positive.

Repℰ2⁡(v→)=(+3+2)

Moving counterclockwise around the origin, we cycle through four regions:

⋯⟶(++)⟶(−+)⟶(−−)⟶(+−)⟶⋯.

Using this basis

B=⟨(01),(−10)⟩

gives the same counterclockwise cycle. We say these two bases have the same orientation.

(a)

Why do they give the same cycle?

(b)

What other configurations of unit vectors on the axes give the same cycle?

(c)

Find the determinants of the matrices formed from those (ordered) bases.

(d)

What other counterclockwise cycles are possible, and what are the associated determinants?

(e)

What happens in ℝ1⁠?

(f)

What happens in ℝ3⁠?

A fascinating general-audience discussion of orientations is in [Gardner].

Answer (6 parts)

(a)

The new basis is the old basis rotated by π/4⁠.

(b)

⟨(−10),(0−1)⟩⁠, ⟨(0−1),(10)⟩

(c)

In each case the determinant is +1 (we say that these bases have positive orientation).

(d)

Because only one sign can change at a time, the only other cycle possible is

⋯⟶(++)⟶(+−)⟶(−−)⟶(−+)⟶⋯.

Here each associated determinant is −1 (we say that such bases have a negative orientation).

(e)

There is one positively oriented basis ⟨(1)⟩ and one negatively oriented basis ⟨(−1)⟩⁠.

(f)

There are 48 bases (6 half-axis choices are possible for the first unit vector, 4 for the second, and 2 for the last). Half are positively oriented like the standard basis on the left below, and half are negatively oriented like the one on the right

In ℝ3 positive orientation is sometimes called ‘right hand orientation’ because if a person places their right hand with their fingers curling from e→1 to e→2 then the thumb will point with e→3⁠.

Four.II.1.30 §II.1 source ↗

This question uses material from the optional Determinant Functions Exist subsection. Prove Theorem 1.5 by using the permutation expansion formula for the determinant.

Answer

We will compare det⁡(s→1,…,s→n) with det⁡(t(s→1),…,t(s→n)) to show that the second differs from the first by a factor of |T|⁠. We represent the s→⁠’s with respect to the standard bases

Repℰn⁡(s→i)=(s1,is2,i⋮sn,i)

and then we represent the map application with matrix-vector multiplication

Repℰn⁡(t(s→i))=(t1,1t1,2…t1,nt2,1t2,2…t2,n⋮tn,1tn,2…tn,n)(s1,js2,j⋮sn,j)=s1,j(t1,1t2,1⋮tn,1)+s2,j(t1,2t2,2⋮tn,2)+⋯+sn,j(t1,nt2,n⋮tn,n)=s1,jt→1+s2,jt→2+⋯+sn,jt→n

where t→i is column i of T⁠. Then det⁡(t(s→1),…,t(s→n)) equals det⁡(s1,1t→1+ s2,1t→2+ …+ sn,1t→n, …, s1,nt→1+ s2,nt→2+ …+ sn,nt→n)⁠.

As in the derivation of the permutation expansion formula, we apply multilinearity, first splitting along the sum in the first argument

det⁡(s1,1t→1,…,s1,nt→1+s2,nt→2+⋯+sn,nt→n)+⋯+det⁡(sn,1t→n,…,s1,nt→1+s2,nt→2+⋯+sn,nt→n)

and then splitting each of those n summands along the sums in the second arguments, etc. We end with, as in the derivation of the permutation expansion, nn summand determinants, each of the form det⁡(si1,1t→i1, si2,2t→i2, …, sin,nt→in)⁠. Factor out each of the si,j⁠’s = si1,1si2,2…sin,n⋅det⁡(t→i1, t→i2, …, t→in)⁠.

As in the permutation expansion derivation, whenever two of the indices in i1⁠, …, in are equal then the determinant has two equal arguments, and evaluates to 0⁠. So we need only consider the cases where i1⁠, …, in form a permutation of the numbers 1⁠, …, n⁠. We thus have

det⁡(t(s→1),…,t(s→n))=∑permutations ϕsϕ(1),1…sϕ(n),ndet⁡(t→ϕ(1),…,t→ϕ(n)).

Swap the columns in det⁡(t→ϕ(1),…,t→ϕ(n)) to get the matrix T back, which changes the sign by a factor of sgn⁡ϕ⁠, and then factor out the determinant of T⁠.

=∑ϕsϕ(1),1…sϕ(n),ndet⁡(t→1,…,t→n)⋅sgn⁡ϕ=det⁡(T)∑ϕsϕ(1),1…sϕ(n),n⋅sgn⁡ϕ.

As in the proof that the determinant of a matrix equals the determinant of its transpose, we commute the s⁠’s to list them by ascending row number instead of by ascending column number (and we substitute sgn⁡(ϕ−1) for sgn⁡(ϕ)⁠).

=det⁡(T)∑ϕs1,ϕ−1(1)…sn,ϕ−1(n)⋅sgn⁡ϕ−1=det⁡(T)det⁡(s→1,s→2,…,s→n)

Four.II.1.31 §II.1 source ↗

(a)

Show that this gives the equation of a line in ℝ2 through (x2,y2) and (x3,y3)⁠.

|xx2x3yy2y3111|=0

(b)

[Petersen] Prove that the area of a triangle with vertices (x1,y1)⁠, (x2,y2)⁠, and (x3,y3) is

12|x1x2x3y1y2y3111|.

(c)

[Math. Mag., Jan. 1973] Prove that the area of a triangle with vertices at (x1,y1)⁠, (x2,y2)⁠, and (x3,y3) whose coordinates are integers has an area of N or N/2 for some positive integer N⁠.

Answer (3 parts)

(a)

An algebraic check is easy.

0=xy2+x2y3+x3y−x3y2−xy3−x2y=x⋅(y2−y3)+y⋅(x3−x2)+x2y3−x3y2

simplifies to the familiar form

y=x⋅(x3−x2)/(y3−y2)+(x2y3−x3y2)/(y3−y2)

(the y3−y2=0 case is easily handled).

For geometric insight, this picture shows that the box formed by the three vectors. Note that all three vectors end in the z=1 plane. Below the two vectors on the right is the line through (x2,y2) and (x3,y3)⁠.

The box will have a nonzero volume unless the triangle formed by the ends of the three is degenerate. That only happens (assuming that (x2,y3)≠(x3,y3)⁠) if (x,y) lies on the line through the other two.

(b)

This is how the answer was given in the cited source. We find the altitude through (x1,y1) of a triangle with vertices (x1,y1) (x2,y2) and (x3,y3) in the usual way from the normal form of the above:

1(x2−x3)2+(y2−y3)2|x1x2x3y1y2y3111|.

Another step shows the area of the triangle to be

12|x1x2x3y1y2y3111|.

This exposition reveals the modus operandi more clearly than the usual proof of showing a collection of terms to be identical with the determinant.

(c)

This is how the answer was given in the cited source. Let

D=|x1x2x3y1y2y3111|

then the area of the triangle is (1/2)|D|⁠. Now if the coordinates are all integers, then D is an integer.

Four.III.1.10 §III.1 source ↗

Find the cofactor.

T=(102−11302−1)

(a)

T2,3

(b)

T3,2

(c)

T1,3

Answer (3 parts)

(a)

(−1)2+3|1002|=−2

(b)

(−1)3+2|12−13|=−5

(c)

(−1)4|−1102|=−2

Four.III.1.11 §III.1 source ↗

Find the adjoint to this matrix.

T=(102−11302−1)

Answer

adj⁡(T)=(|132−1|−|022−1||0213|−|−130−1||120−1|−|12−13||−1102|−|1002||10−11|)=(−74−2−1−1−5−2−21)

Four.III.1.12 §III.1 source ↗

This determinant is 0⁠. Compute that by expanding on the first row.

|123456789|

Answer

|123456789|=1⋅(+1)⋅|5689|+2⋅(−1)⋅|4679|+3⋅(+1)⋅|4578|

The formula for 2×2 matrices gives 1⋅(−3)−2⋅(−6)+3⋅(−3)=0⁠.

Four.III.1.13 §III.1 source ↗

Find the determinant by expanding

|301122−130|

(a)

on the first row

(b)

on the second row

(c)

on the third column.

Answer (3 parts)

(a)

3⋅(+1)|2230|+ 0⋅(−1)|12−10|+ 1⋅(+1)|12−13|= − 13

(b)

1⋅(−1)|0130|+ 2⋅(+1)|31−10|+ 2⋅(−1)|30−13|= − 13

(c)

1⋅(+1)|12−13|+ 2⋅(−1)|30−13|+ 0⋅(+1)|3012|= − 13

Four.III.1.14 §III.1 source ↗

Find the adjoint of the matrix in Example 1.5.

Answer

This is adj⁡(T)⁠.

(T1,1T2,1T3,1T1,2T2,2T3,2T1,3T2,3T3,3)=(|5689|−|2389|+|2356|−|4679|+|1379|−|1346|+|4578|−|1278|+|1245|)=(−36−36−126−36−3)

Four.III.1.15 §III.1 source ↗

Find the matrix adjoint to each.

(a)

(214−102101)

(b)

(3−124)

(c)

(1150)

(d)

(143−103189)

Answer (4 parts)

(a)

This is the adjoint.

(T1,1T2,1T3,1T1,2T2,2T3,2T1,3T2,3T3,3)=(|0201|−|1401||1402|−|−1211||2411|−|24−12||−1010|−|2110||21−10|)=(0−123−2−8011)

(b)

The minors are 1×1⁠. (T1,1T2,1T1,2T2,2)= (|4|−|−1|−|2||3|)= (41−23)

(c)

(0−1−51)

(d)

The minors are 2×2⁠.

(T1,1T2,1T3,1T1,2T2,2T3,2T1,3T2,3T3,3)=(|0389|−|4389||4303|−|−1319||1319|−|13−13||−1018|−|1418||14−10|)=(−24−1212126−6−8−44)

Four.III.1.16 §III.1 source ↗

Find the inverse of each matrix in the prior question with Theorem 1.8.

Answer (4 parts)

(a)

(1/3)⋅(0−123−2−8011)= (0−1/32/31−2/3−8/301/31/3)

(b)

(1/14)⋅(41−23)= (2/71/14−1/73/14)

(c)

(1/−5)⋅(0−1−51)= (01/51−1/5)

(d)

The matrix has a zero determinant, and so has no inverse.

Four.III.1.17 §III.1 source ↗

Find the matrix adjoint to this one.

(2100121001210012)

Answer

(T1,1T2,1T3,1T4,1T1,2T2,2T3,2T4,2T1,3T2,3T3,3T4,3T1,4T2,4T3,4T4,4)= (4−32−1−36−422−46−3−12−34)

Four.III.1.18 §III.1 source ↗

Expand across the first row to derive the formula for the determinant of a 2×2 matrix.

Answer

The determinant

|abcd|

expanded on the first row gives a⋅(+1)|d|+ b⋅(−1)|c|= ad− bc (note the two 1×1 minors).

Four.III.1.19 §III.1 source ↗

Expand across the first row to derive the formula for the determinant of a 3×3 matrix.

Answer

The determinant of

(abcdefghi)

is this.

a⋅|efhi|−b⋅|dfgi|+c⋅|degh|=a(ei−fh)−b(di−fg)+c(dh−eg)

Four.III.1.20 §III.1 source ↗

(a)

Give a formula for the adjoint of a 2×2 matrix.

(b)

Use it to derive the formula for the inverse.

Answer (2 parts)

(a)

(T1,1T2,1T1,2T2,2)= (|t2,2|−|t1,2|−|t2,1||t1,1|)= (t2,2−t1,2−t2,1t1,1)

(b)

(1/t1,1t2,2− t1,2t2,1)⋅(t2,2−t1,2−t2,1t1,1)

Four.III.1.21 §III.1 source ↗

Can we compute a determinant by expanding down the diagonal?

Answer

No. Here is a determinant whose value

|100010001|=1

doesn’t equal the result of expanding down the diagonal.

1⋅(+1)|1001|+1⋅(+1)|1001|+1⋅(+1)|1001|=3

Four.III.1.22 §III.1 source ↗

Give a formula for the adjoint of a diagonal matrix.

Answer

Consider this diagonal matrix.

D=(d100…0d2000d3⋱dn)

If i≠j then the i,j⁠ minor is an (n−1)×(n−1) matrix with only n−2 nonzero entries, because we have deleted both di and dj⁠. Thus, at least one row or column of the minor is all zeroes, and so the cofactor Di,j is zero. If i=j then the minor is the diagonal matrix with entries d1⁠, …, di−1⁠, di+1⁠, …, dn⁠. Its determinant is obviously (−1)i+j=(−1)2i=1 times the product of those.

adj⁡(D)=(d2⋯dn000d1d3⋯dn0⋱d1⋯dn−1)

By the way, Theorem 1.8 provides a slicker way to derive this conclusion.

Four.III.1.23 §III.1 source ↗

Prove that the transpose of the adjoint is the adjoint of the transpose.

Answer

Just note that if S=T𝖳 then the cofactor Sj,i equals the cofactor Ti,j because (−1)j+i=(−1)i+j and because the minors are the transposes of each other (and the determinant of a transpose equals the determinant of the matrix).

Four.III.1.24 §III.1 source ↗

Prove or disprove: adj⁡(adj⁡(T))=T⁠.

Answer

It is false; here is an example.

T=(123456789)adj⁡(T)=(−36−36−126−36−3)adj⁡(adj⁡(T))=(000000000)

Four.III.1.25 §III.1 source ↗

A square matrix is upper triangular if each i,j entry is zero in the part above the diagonal, that is, when i>j⁠.

(a)

Must the adjoint of an upper triangular matrix be upper triangular? Lower triangular?

(b)

Prove that the inverse of a upper triangular matrix is upper triangular, if an inverse exists.

Answer (2 parts)

(a)

An example

M=(123045006)

suggests the right answer.

adj⁡(M)=(M1,1M2,1M3,1M1,2M2,2M3,2M1,3M2,3M3,3)=(|4506|−|2306||2345|−|0506||1306|−|1305||0400|−|1200||1204|)=(24−12−206−5004)

The result is indeed upper triangular.

This check is detailed but not hard. The entries in the upper triangle of the adjoint are Ma,b where a>b⁠. We need to verify that the cofactor Ma,b is zero if a>b⁠. With a>b⁠, row a and column b of M⁠,

(m1,1…m1,b…m2,1…m2,b⋮⋮ma,1…ma,b…ma,n⋮⋮mn,b)

when deleted, leave an upper triangular minor, because entry i,j of the minor is either entry i,j of M (this happens if a>i and b>j⁠; in this case i<j implies that the entry is zero) or it is entry i,j+1 of M (this happens if i<a and j>b⁠; in this case, i<j implies that i<j+1⁠, which implies that the entry is zero), or it is entry i+1,j+1 of M (this last case happens when i>a and j>b⁠; obviously here i<j implies that i+1<j+1 and so the entry is zero). Thus the determinant of the minor is the product down the diagonal. Observe that the a−1,a entry of M is the a−1,a−1⁠ entry of the minor (it doesn’t get deleted because the relation a>b is strict). But this entry is zero because M is upper triangular and a−1<a⁠. Therefore the cofactor is zero, and the adjoint is upper triangular. (The lower triangular case is similar.)

(b)

This is immediate from the prior part, by Theorem 1.8.

Four.III.1.26 §III.1 source ↗

This question requires material from the optional Determinants Exist subsection. Prove Theorem 1.4 by using the permutation expansion.

Answer

We will show that each determinant can be expanded along row i⁠. The argument for column j is similar.

Each term in the permutation expansion contains one and only one entry from each row. As in Example 1.1, factor out each row i entry to get |T|= ti,1⋅Tˆi,1+ ⋯+ ti,n⋅Tˆi,n⁠, where each Tˆi,j is a sum of terms not containing any elements of row i⁠. We will show that Tˆi,j is the i,j cofactor.

Consider the i,j=n,n case first:

tn,n⋅Tˆn,n=tn,n⋅∑ϕt1,ϕ(1)t2,ϕ(2)…tn−1,ϕ(n−1)sgn⁡(ϕ)

where the sum is over all n⁠-permutations ϕ such that ϕ(n)=n⁠. To show that Tˆi,j is the minor Ti,j⁠, we need only show that if ϕ is an n⁠-permutation such that ϕ(n)=n and σ is an n−1⁠-permutation with σ(1)=ϕ(1)⁠, …, σ(n−1)=ϕ(n−1) then sgn⁡(σ)=sgn⁡(ϕ)⁠. But that’s true because ϕ and σ have the same number of inversions.

Back to the general i,j case. Swap adjacent rows until the i⁠-th is last and swap adjacent columns until the j⁠-th is last. Observe that the determinant of the i,j⁠-th minor is not affected by these adjacent swaps because inversions are preserved (since the minor has the i⁠-th row and j⁠-th column omitted). On the other hand, the sign of |T| and Tˆi,j changes n−i plus n−j times. Thus Tˆi,j= (−1)n−i+n−j|Ti,j|= (−1)i+j|Ti,j|⁠.

Four.III.1.27 §III.1 source ↗

Prove that the determinant of a matrix equals the determinant of its transpose using Laplace’s expansion and induction on the size of the matrix.

Answer

This is obvious for the 1×1 base case.

For the inductive case, assume that the determinant of a matrix equals the determinant of its transpose for all 1×1⁠, …, (n−1)×(n−1) matrices. Expanding on row i gives |T|=ti,1Ti,1+…+ti,nTi,n and expanding on column i gives |T𝖳|= t1,i(T𝖳)1,i+ ⋯+ tn,i(T𝖳)n,i Since (−1)i+j=(−1)j+i the signs are the same in the two summations. Since the j,i minor of T𝖳 is the transpose of the i,j minor of T⁠, the inductive hypothesis gives |(T𝖳)i,j|=|Ti,j|⁠.

Four.III.1.28 §III.1 source ↗

Show that

Fn=|1−11−11−1…110101…011010…001101…......…|

where Fn is the n⁠-th term of 1,1,2,3,5,…,x,y,x+y,…⁠, the Fibonacci sequence, and the determinant is of order n−1⁠. [Am. Math. Mon., Jun. 1949]

Answer

This is how the answer was given in the cited source. Denoting the above determinant by Dn⁠, it is seen that D2=1⁠, D3=2⁠. It remains to show that Dn=Dn−1+Dn−2,n≥4⁠. In Dn subtract the (n−3)⁠-th column from the (n−1)⁠-th, the (n−4)⁠-th from the (n−2)⁠-th, …, the first from the third, obtaining

Fn=|1−10000…11−1000…011−100…0011−10…......…|.

By expanding this determinant with reference to the first row, there results the desired relation.

1 Topic: Cramer’s Rule source ↗

Use Cramer’s Rule to solve each for each of the variables.

(a)

x−y=4−x+2y=−7

(b)

−2x+y=−2x−2y=−2

Answer (2 parts)

(a)

Solve for the variables separately.

x=|4−1−72||1−1−12|=11=1y=|14−1−7||1−1−12|=−31=−3

(b)

x=2⁠, y=2

2 Topic: Cramer’s Rule source ↗

Use Cramer’s Rule to solve this system for z⁠.

2x+y+z=13x+z=4x−y−z=2

Answer

z=1

3 Topic: Cramer’s Rule source ↗

Prove Cramer’s Rule.

Answer

Determinants are unchanged by combinations, including column combinations, so det⁡(Bi)= det⁡(a→1, …, x1a→1+ ⋯+ xia→i+ ⋯+ xna→n, …, a→n)⁠. Use the operation of taking −x1 times the first column and adding it to the i⁠-th column, etc., to see this is equal to det⁡(a→1,…,xia→i,…,a→n)⁠. In turn, that is equal to xi⋅det⁡(a→1, …, a→i, …, a→n)= xi⋅det⁡(A)⁠, as required.

4 Topic: Cramer’s Rule source ↗

Here is an alternative proof of Cramer’s Rule that doesn’t overtly contain any geometry. Write Xi for the identity matrix with column i replaced by the vector x→ of unknowns x1⁠, …, xn⁠.

(a)

Observe that AXi=Bi⁠.

(b)

Take the determinant of both sides.

Answer (2 parts)

(a)

Here is the case of a 2×2 system with i=2⁠.

a1,1x1+a1,2x2=b1a2,1x1+a2,2x2=b2⟺(a1,1a1,2a2,1a2,2)(1x10x2)=(a1,1b1a2,1b2)

(b)

The determinant function is multiplicative det⁡(Bi)= det⁡(AXi)= det⁡(A)⋅det⁡(Xi)⁠. The Laplace expansion shows that det⁡(Xi)=xi⁠, and solving for xi gives Cramer’s Rule.

5 Topic: Cramer’s Rule source ↗

Suppose that a linear system has as many equations as unknowns, that all of its coefficients and constants are integers, and that its matrix of coefficients has determinant 1⁠. Prove that the entries in the solution are all integers. (Remark. This is often used to invent linear systems for exercises.)

Answer

Because the determinant of A is nonzero, Cramer’s Rule applies and shows that xi=|Bi|/1⁠. Since Bi is a matrix of integers, its determinant is an integer.

6 Topic: Cramer’s Rule source ↗

Use Cramer’s Rule to give a formula for the solution of a two equations/two unknowns linear system.

Answer

The solution of

ax+by=ecx+dy=f

is

x=ed−fbad−bcy=af−ecad−bc

provided of course that the denominators are not zero.

7 Topic: Cramer’s Rule source ↗

Can Cramer’s Rule tell the difference between a system with no solutions and one with infinitely many?

Answer

Of course, singular systems have |A| equal to zero, but we can characterize the infinitely many solutions case is by the fact that all of the |Bi| are zero as well.

8 Topic: Cramer’s Rule source ↗

The first picture in this Topic (the one that doesn’t use determinants) shows a unique solution case. Produce a similar picture for the case of infinitely many solutions, and the case of no solutions.

Answer

We can consider the two nonsingular cases together with this system

x1+2x2=6x1+2x2=c

where c=6 of course yields infinitely many solutions, and any other value for c yields no solutions. The corresponding vector equation

x1⋅(11)+x2⋅(22)=(6c)

gives a picture of two overlapping vectors. Both lie on the line y=x⁠. In the c=6 case the vector on the right side also lies on the line y=x but in any other case it does not.

1 Topic: Speed of Calculating Determinants source ↗

To get an idea of what happens for typical matrices we can use the ability of computer systems to generate random numbers (of course, these are only pseudo-random in that they come from an algorithm but they pass a number of reasonable statistical tests for randomness).

(a)

Fill a 5×5 array with random numbers say, in the range [0..1)⁠). See if it is singular. Repeat that experiment a few times. Are singular matrices frequent or rare in this sense?

(b)

Time your computer algebra system at finding the determinant of ten 10×10 arrays of random numbers. Find the average time per array. Repeat the prior item for 20×20 arrays, 30×30 arrays, … 100×100 arrays, and compare to the numbers given above. (Notice that, when an array is singular, we can sometimes decide that quickly, for instance if the first row equals the second. In the light of your answer to the first part, do you expect that singular systems play a large role in your average?)

(c)

Graph the input size versus the average time.

Answer

Your timing will depend in part on the computer algebra system that you use, and in part on the power of the computer on which you do the calculation. But you should get a curve that is similar to the one shown.

2 Topic: Speed of Calculating Determinants source ↗

Compute the determinant of each of these by hand using the two methods discussed above.

(a)

|215−3|

(b)

|311−105−12−2|

(c)

|210013200−1−2100−21|

Count the number of multiplications and divisions used in each case, for each of the methods.

Answer (4 parts)

The number of operations depends on exactly how we do the operations.

(a)

The determinant is −11⁠. To row reduce takes a single row combination with two multiplications (−5/2 times 2 plus 5⁠, and −5/2 times 1 plus −3⁠) and the product down the diagonal takes one more multiplication. The permutation expansion takes two multiplications (2 times −3 and 5 times 1⁠).

(b)

The determinant is −39⁠. Counting the operations is routine.

(c)

The determinant is 4⁠.

3 Topic: Speed of Calculating Determinants source ↗

The use by the timing routine of do_matrix has a bug. That routine does two things, generate a random matrix and then do gauss_method on it, and the timing number returned is for the combination. Produce code that times only the gauss_method routine.

Answer

You should build the matrix once and then time it being reduced many times. Note that gauss_method, as written, changes the input matrix, so be careful not to loop using the output of that routine or all the loops but the first will be on an echelon form matrix.

4 Topic: Speed of Calculating Determinants source ↗

What 10×10 array can you invent that takes your computer the longest time to reduce? The shortest?

Answer

The identity matrix typically does not take long to reduce. For matrices that are slow, in Python you can try ones with numbers that are quite large. (Some computer algebra system also have the ability to generate special matrices; try some of those.)

5 Topic: Speed of Calculating Determinants source ↗

Some computer language specifications requires that arrays be stored “by column,” that is, the entire first column is stored contiguously, then the second column, etc. Does the code fragment given take advantage of this, or can it be rewritten to make it faster, by taking advantage of the fact that computer fetches are faster from contiguous locations?

Answer

No, the above code handles the numbers by row.

1 Topic: Chiò’s Method source ↗

Use Chiò’s Method to find each determinant.

(a)

|123456789|

(b)

|2140014011110211|

Answer (2 parts)

(a)

Chiò’s matrix is

C=(−3−6−6−12)

and its determinant is 0

(b)

Start with

C3=(2801−22422)

and then the next step

C2=(−124−284)

with determinant det⁡(C2)=64⁠. The determinant of the original matrix is thus 64/(22⋅21)=8

2 Topic: Chiò’s Method source ↗

What if a1,1 is zero?

Answer

The same construction as was used for the 3×3⁠ case above shows that in place of a1,1 we can select any nonzero entry ai,j⁠. Entry cp,q of Chiò’s matrix is the value of this determinant

|a1,1a1,q+1ap+1,1ap+1,q+1|

where p+1≠i and q+1≠j⁠.

3 Topic: Chiò’s Method source ↗

The Rule of Sarrus is a mnemonic that many people learn for the 3×3 determinant formula. To the right of the matrix, copy the first two columns.

abcabdefdeghigh

Then the determinant is the sum of the three upper-left to lower-right diagonals minus the three lower-left to upper-right diagonals aei+bfg+cdh−gec−hfa−idb⁠. Count the operations involved in Sarrus’s formula and in Chiò’s.

Answer

Sarrus’s formula uses 12⁠ multiplications and 5 additions (including the subtractions in with the additions). Chiò’s formula uses two multiplications and an addition (which is actually a subtraction) for each of the four 2×2⁠ determinants, and another two multiplications and an addition for the 2×2⁠ Chió’s determinant, as well as a final division by a1,1⁠. That’s eleven multiplication/divisions and five addition/subtractions. So Chiò is the winner.

4 Topic: Chiò’s Method source ↗

Prove Chiò’s formula.

Answer

Consider an n×n matrix.

A=(a1,1a1,2⋯a1,n−1a1,na2,1a2,2⋯a2,n−1a2,n⋮an−1,1an−1,2⋯an−1,n−1an−1,nan,1an,2⋯an,n−1an,n)

Rescale every row but the first by a1,1⁠.

⟶a1,1ρ3⋮a1,1ρna1,1ρ2(a1,1a1,2⋯a1,n−1a1,na2,1a1,1a2,2a1,1⋯a2,n−1a1,1a2,na1,1⋮an−1,1a1,1an−1,2a1,1⋯an−1,n−1a1,1an−1,na1,1an,1a1,1an,2a1,1⋯an,n−1a1,1an,na1,1)

That rescales the determinant by a factor of a1,1n−1⁠.

Next perform the row operation −ai,1ρ1+ρi on each row i>1⁠. These row operations don’t change the determinant.

⟶−a3,1ρ1+ρ3⋮−an,1ρ1+ρn−a2,1ρ1+ρ2

The result is a matrix B whose first row is unchanged, whose first column is all zeros (except for the 1,1 entry of a1,1⁠), and whose remaining entries are these.

(a1,2⋯a1,n−1a1,na2,2a1,1−a2,1a1,2⋯a2,n−1an,n−a2,n−1a1,n−1a2,nan,n−a2,na1,n⋮an,2a1,1−an,1a1,2⋯an,n−1a1,1−an,1a1,n−1an,na1,1−an,1a1,n)

The determinant of this n×n matrix B is a1,1n−1⁠ times the determinant of A⁠.

Denote by C the 1,1 minor of the matrix B⁠, that is, the submatrix consisting of the last n−1 rows and columns. The Laplace expansion down the first column of B gives that its determinant is det⁡(B)= (−1)1+1a1,1det⁡(C)⁠.

If a1,1≠0 then setting equal the two expressions for {(}B) and canceling gives det⁡(A)=det⁡(C)/an,nn−2⁠.

1 Topic: Projective Geometry source ↗

What is the equation of this point?

(100)

Answer

From the dot product

0=(100)∙(L1L2L3)=L1

we get that the equation is L1=0⁠.

2 Topic: Projective Geometry source ↗

(a)

Find the line incident on these points in the projective plane.

(123),(456)

(b)

Find the point incident on both of these projective lines.

(123),(456)

Answer (2 parts)

(a)

This determinant

0=|14x25y36z|=−3x+6y−3z

shows that the line is L=(−36−3)⁠.

(b)

(−36−3)

3 Topic: Projective Geometry source ↗

Find the formula for the line incident on two projective points. Find the formula for the point incident on two projective lines.

Answer

The line incident on

u=(u1u2u3)v=(v1v2v3)

comes from this determinant equation.

0=|u1v1xu2v2yu3v3z|=(u2v3−u3v2)⋅x+(u3v1−u1v3)⋅y+(u1v2−u2v1)⋅z

The equation for the point incident on two lines is the same.

4 Topic: Projective Geometry source ↗

Prove that the definition of incidence is independent of the choice of the representatives of p and L⁠. That is, if p1⁠, p2⁠, p3⁠, and q1⁠, q2⁠, q3 are two triples of homogeneous coordinates for p⁠, and L1⁠, L2⁠, L3⁠, and M1⁠, M2⁠, M3 are two triples of homogeneous coordinates for L⁠, prove that p1L1+p2L2+p3L3=0 if and only if q1M1+q2M2+q3M3=0⁠.

Answer

If p1⁠, p2⁠, p3⁠, and q1⁠, q2⁠, q3 are two triples of homogeneous coordinates for p then the two column vectors are in proportion, that is, lie on the same line through the origin. Similarly, the two row vectors are in proportion.

k⋅(p1p2p3)=(q1q2q3)m⋅(L1L2L3)=(M1M2M3)

Then multiplying gives the answer (km)⋅(p1L1+p2L2+p3L3)= q1M1+ q2M2+ q3M3= 0⁠.

5 Topic: Projective Geometry source ↗

Give a drawing to show that central projection does not preserve circles, that a circle may project to an ellipse. Can a (non-circular) ellipse project to a circle?

Answer

The picture of the solar eclipse — unless the image plane is exactly perpendicular to the line from the sun through the pinhole — shows the circle of the sun projecting to an image that is an ellipse. (Another example is that in many pictures in this Topic, we’ve shown the circle that is the sphere’s equator as an ellipse, that is, a viewer of the drawing sees a circle as an ellipse.)

The solar eclipse picture also shows the converse. If we picture the projection as going from left to right through the pinhole then the ellipse I projects through P to a circle S⁠.

6 Topic: Projective Geometry source ↗

Give the formula for the correspondence between the non-equatorial part of the antipodal modal of the projective plane, and the plane z=1⁠.

Answer

A spot on the unit sphere

(p1p2p3)

is non-equatorial if and only if p3≠0⁠. In that case it corresponds to this point on the z=1 plane

(p1/p3p2/p31)

since that is intersection of the line containing the vector and the plane.

7 Topic: Projective Geometry source ↗

(Pappus’s Theorem) Assume that T0⁠, U0⁠, and V0 are collinear and that T1⁠, U1⁠, and V1 are collinear. Consider these three points: (i) the intersection V2 of the lines T0U1 and T1U0⁠, (ii) the intersection U2 of the lines T0V1 and T1V0⁠, and (iii) the intersection T2 of U0V1 and U1V0⁠.

(a)

Draw a (Euclidean) picture.

(b)

Apply the lemma used in Desargue’s Theorem to get simple homogeneous coordinate vectors for the T⁠’s and V0⁠.

(c)

Find the resulting homogeneous coordinate vectors for U⁠’s (these must each involve a parameter as, e.g., U0 could be anywhere on the T0V0 line).

(d)

Find the resulting homogeneous coordinate vectors for V1⁠. (Hint: it involves two parameters.)

(e)

Find the resulting homogeneous coordinate vectors for V2⁠. (It also involves two parameters.)

(f)

Show that the product of the three parameters is 1⁠.

(g)

Verify that V2 is on the T2U2 line.

Answer (7 parts)

(a)

Other pictures are possible, but this is one.

The intersections T0U1∩T1U0=V2⁠, T0V1∩T1V0=U2⁠, and U0V1∩U1V0=T2 are labeled so that on each line is a T⁠, a U⁠, and a V⁠.

(b)

The lemma used in Desargue’s Theorem gives a basis B with respect to which the points have these homogeneous coordinate vectors.

RepB⁡(t→0)=(100)RepB⁡(t→1)=(010)RepB⁡(t→2)=(001)RepB⁡(v→0)=(111)

(c)

First, any U0 on T0V0

RepB⁡(u→0)=a(100)+b(111)=(a+bbb)

has homogeneous coordinate vectors of this form

(u011)

(u0 is a parameter; it depends on where on the T0V0 line the point U0 is, but any point on that line has a homogeneous coordinate vector of this form for some u0∈ℝ⁠). Similarly, U2 is on T1V0

RepB⁡(u→2)=c(010)+d(111)=(dc+dd)

and so has this homogeneous coordinate vector.

(1u21)

Also similarly, U1 is incident on T2V0

RepB⁡(u→1)=e(001)+f(111)=(ffe+f)

and has this homogeneous coordinate vector.

(11u1)

(d)

Because V1 is T0U2∩U0T2 we have this.

g(100)+h(1u21)=i(u011)+j(001)⟹g+h=iu0hu2=ih=i+j

Substituting hu2 for i in the first equation

(hu0u2hu2h)

shows that V1 has this two-parameter homogeneous coordinate vector.

(u0u2u21)

(e)

Since V2 is the intersection T0U1∩T1U0

k(100)+l(11u1)=m(010)+n(u011)⟹k+l=nu0l=m+nlu1=n

and substituting lu1 for n in the first equation

(lu0u1llu1)

gives that V2 has this two-parameter homogeneous coordinate vector.

(u0u11u1)

(f)

Because V1 is on the T1U1 line its homogeneous coordinate vector has the form

p(010)+q(11u1)=(qp+qqu1)(∗)

but a previous part of this question established that V1⁠’s homogeneous coordinate vectors have the form

(u0u2u21)

and so this a homogeneous coordinate vector for V1⁠.

(u0u1u2u1u2u1)(∗∗)

By (∗⁠) and (∗∗⁠), there is a relationship among the three parameters: u0u1u2=1⁠.

(g)

The homogeneous coordinate vector of V2 can be written in this way.

(u0u1u2u2u1u2)=(1u2u1u2)

Now, the T2U2 line consists of the points whose homogeneous coordinates have this form.

r(001)+s(1u21)=(ssu2r+s)

Taking s=1 and r=u1u2−1 shows that the homogeneous coordinate vectors of V2 have this form.

1 Topic: Computer Graphics source ↗

Calculate the product in (∗⁠).

Answer

Straightforward calculation gives this for the product.

(100012001)(0−10−100001)(10001−2001)=(0−12−102001)

2 Topic: Computer Graphics source ↗

Find the matrix that reflects about the line y=2x⁠.

Answer

Working in ℝ2⁠, let the matrix be M⁠. We get these two.

(12)↦(12)(−21)↦(2−1)

(For the second one, the starting vector is on the line through the origin that is perpendicular to y=2x⁠.) We have this.

M(1−221)=(122−1)

Solving gives

M=(122−1)(1−221)−1=(−3/54/54/53/5)

and so for homogeneous coordinates the matrix is this.

(−3/54/504/53/50001)

3 Topic: Computer Graphics source ↗

Find the matrix that reflects about the line y=2x−4⁠.

Answer

Move all points over by 4⁠, reflect about the line y=2x using the prior exercise, and them move them back.

(10001−4001)(−3/54/504/53/50001)(100014001)=(−3/54/516/54/53/5−8/5001)

4 Topic: Computer Graphics source ↗

Rotation and translation are rigid operations. What is the matrix for a rotation followed by a translation?

Answer

(10tx01ty001)(cos⁡θ−sin⁡θ0sin⁡θcos⁡θ0001)=(cos⁡θ−sin⁡θtxsin⁡θcos⁡θty001)

5 Topic: Computer Graphics source ↗

The homogeneous coordinates extend to manipulations of three dimensional space in the obvious way: every coordinate is a set of four-tall nonzero vectors that are related by being scalar multiples of each other. Give the matrix to do rotation about the z⁠ axis, and the matrix for rotation about the y⁠ axis.

Answer

Here are the two 4×4 matrices.

(cos⁡θ−sin⁡θ00sin⁡θcos⁡θ0000000001)(cos⁡θ0−sin⁡θ00000sin⁡θ0cos⁡θ00001)

Loading chapter 5…

Exercises and answers are adapted from Linear Algebra by Jim Hefferon, which is free to read and download, and are used under the CC BY-SA 2.5 licence. They are taken from the book’s LaTeX source and set for the web, so they may differ slightly from the printed page. This adaptation is released under the same licence and is neither affiliated with nor endorsed by Jim Hefferon.