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

Chapter 1: Linear Systems

153 problems · 16 figures

One.I.1.17 §I.1 source ↗

Use Gauss’s Method to find the unique solution for each system.

(a)

2x+3y=13x−y=−1

(b)

x−z=03x+y=1−x+y+z=4

Answer (2 parts)

(a)

Gauss’s Method

⟶−(1/2)ρ1+ρ22x+3y=13−(5/2)y=−15/2

gives that the solution is y=3 and x=2⁠.

(b)

Gauss’s Method here

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

gives x=−1⁠, y=4⁠, and z=−1⁠.

One.I.1.18 §I.1 source ↗

Each system is in echelon form. For each, say whether the system has a unique solution, no solution, or infinitely many solutions.

(a)

−3x+2y=0−2y=0

(b)

x+y=4y−z=0

(c)

x+y=4y−z=00=0

(d)

x+y=40=4

(e)

3x+6y+z=−0.5−z=2.5

(f)

x−3y=20=0

(g)

2x+2y=4y=10=4

(h)

2x+y=0

(i)

x−y=−10=00=4

(j)

x+y−3z=−1y−z=2z=00=0

Answer (11 parts)

If a system has a contradictory equation then it has no solution. Otherwise, if there are any variables that are not leading a row then it has infinitely many solution. In the final case, where there is no contradictory equation and every variable leads some row, it has a unique solution.

(a)

Unique solution

(b)

Infinitely many solutions

(c)

Infinitely many solutions

(d)

No solution

(e)

Infinitely many solutions

(f)

Infinitely many solutions

(g)

No solution

(h)

Infinitely many solutions

(i)

No solution

(j)

Unique solution

One.I.1.19 §I.1 source ↗

Use Gauss’s Method to solve each system or conclude ‘many solutions’ or ‘no solutions’.

(a)

2x+2y=5x−4y=0

(b)

−x+y=1x+y=2

(c)

x−3y+z=1x+y+2z=14

(d)

−x−y=1−3x−3y=2

(e)

4y+z=202x−2y+z=0x+z=5x+y−z=10

(f)

2x+z+w=5y−w=−13x−z−w=04x+y+2z+w=9

Answer (6 parts)

(a)

Gaussian reduction

⟶−(1/2)ρ1+ρ22x+2y=5−5y=−5/2

shows that y=1/2 and x=2 is the unique solution.

(b)

Gauss’s Method

⟶ρ1+ρ2−x+y=12y=3

gives y=3/2 and x=1/2 as the only solution.

(c)

Row reduction

⟶−ρ1+ρ2x−3y+z=14y+z=13

shows, because the variable z is not a leading variable in any row, that there are many solutions.

(d)

Row reduction

⟶−3ρ1+ρ2−x−y=10=−1

shows that there is no solution.

(e)

Gauss’s Method

⟶ρ1↔ρ4x+y−z=102x−2y+z=0x+z=54y+z=20⟶−ρ1+ρ3−2ρ1+ρ2x+y−z=10−4y+3z=−20−y+2z=−54y+z=20⟶ρ2+ρ4−(1/4)ρ2+ρ3x+y−z=10−4y+3z=−20(5/4)z=04z=0

gives the unique solution (x,y,z)=(5,5,0)⁠.

(f)

Here Gauss’s Method gives

⟶−2ρ1+ρ4−(3/2)ρ1+ρ32x++z+w=5y−w=−1−(5/2)z−(5/2)w=−15/2y−w=−1⟶−ρ2+ρ42x++z+w=5y−w=−1−(5/2)z−(5/2)w=−15/20=0

which shows that there are many solutions.

One.I.1.20 §I.1 source ↗

Solve each system or conclude ‘many solutions’ or ‘no solutions’. Use Gauss’s Method.

(a)

x+y+z=5x−y=0y+2z=7

(b)

3x+z=7x−y+3z=4x+2y−5z=−1

(c)

x+3y+z=0−x−y=2−x+y+2z=8

Answer (3 parts)

(a)

Gauss’s Method

x+y+z=5x−y=0y+2z=7⟶−ρ1+ρ2x+y+z=5−2y−z=−5y+2z=7⟶(1/2)ρ2+ρ3x+y+z=5−2y−1=−5(3/2)z=9/2

followed by back-substitution gives x=1⁠, y=1⁠, and z=3⁠.

(b)

Here Gauss’s Method

3x+z=7x−y+3z=4x+2y−5z=−1⟶−(1/3)ρ1+ρ3−(1/3)ρ1+ρ23x+z=7−y+(8/3)z=5/32y−(16/3)z=−(10/3)⟶2ρ2+ρ33x+z=7−y+(8/3)z=5/30=0

finds that the variable z does not lead a row. There are infinitely many solutions.

(c)

The steps

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

give (x,y,z)=(−1,−1,4)⁠.

One.I.1.21 §I.1 source ↗

We can solve linear systems by methods other than Gauss’s. One often taught in high school is to solve one of the equations for a variable, then substitute the resulting expression into other equations. Then we repeat that step until there is an equation with only one variable. From that we get the first number in the solution and then we get the rest with back-substitution. This method takes longer than Gauss’s Method, since it involves more arithmetic operations, and is also more likely to lead to errors. To illustrate how it can lead to wrong conclusions, we will use the system

x+3y=12x+y=−32x+2y=0

from Example 1.13.

(a)

Solve the first equation for x and substitute that expression into the second equation. Find the resulting y⁠.

(b)

Again solve the first equation for x⁠, but this time substitute that expression into the third equation. Find this y⁠.

What extra step must a user of this method take to avoid erroneously concluding a system has a solution?

Answer (3 parts)

(a)

From x=1−3y we get that 2(1−3y)+y=−3⁠, giving y=1⁠.

(b)

From x=1−3y we get that 2(1−3y)+2y=0⁠, leading to the conclusion that y=1/2⁠.

Users of this method must check any potential solutions by substituting back into all the equations.

One.I.1.22 §I.1 source ↗

For which values of k are there no solutions, many solutions, or a unique solution to this system?

x−y=13x−3y=k

Answer

Do the reduction

⟶−3ρ1+ρ2x−y=10=−3+k

to conclude this system has no solutions if k≠3 and if k=3 then it has infinitely many solutions. It never has a unique solution.

One.I.1.23 §I.1 source ↗

This system is not linear in that it says sin⁡α instead of α

2sin⁡α−cos⁡β+3tan⁡γ=34sin⁡α+2cos⁡β−2tan⁡γ=106sin⁡α−3cos⁡β+tan⁡γ=9

and yet we can apply Gauss’s Method. Do so. Does the system have a solution?

Answer

Let x=sin⁡α⁠, y=cos⁡β⁠, and z=tan⁡γ⁠:

2x−y+3z=34x+2y−2z=106x−3y+z=9⟶−3ρ1+ρ3−2ρ1+ρ22x−y+3z=34y−8z=4−8z=0

gives z=0⁠, y=1⁠, and x=2⁠. Note that no α∈ℝ satisfies that sin⁡α=2⁠.

One.I.1.24 §I.1 source ↗

What conditions must the constants, the b⁠’s, satisfy so that each of these systems has a solution? Hint. Apply Gauss’s Method and see what happens to the right side.

(a)

x−3y=b13x+y=b2x+7y=b32x+4y=b4

(b)

x1+2x2+3x3=b12x1+5x2+3x3=b2x1+8x3=b3

Answer (2 parts)

(a)

Gauss’s Method

⟶−ρ1+ρ3−2ρ1+ρ4−3ρ1+ρ2x−3y=b110y=−3b1+b210y=−b1+b310y=−2b1+b4⟶−ρ2+ρ4−ρ2+ρ3x−3y=b110y=−3b1+b20=2b1−b2+b30=b1−b2+b4

shows that this system is consistent if and only if both b3=−2b1+b2 and b4=−b1+b2⁠.

(b)

Reduction

⟶−ρ1+ρ3−2ρ1+ρ2x1+2x2+3x3=b1x2−3x3=−2b1+b2−2x2+5x3=−b1+b3⟶2ρ2+ρ3x1+2x2+3x3=b1x2−3x3=−2b1+b2−x3=−5b1+2b2+b3

shows that each of b1⁠, b2⁠, and b3 can be any real number — this system always has a unique solution.

One.I.1.25 §I.1 source ↗

True or false: a system with more unknowns than equations has at least one solution. (As always, to say ‘true’ you must prove it, while to say ‘false’ you must produce a counterexample.)

Answer

This system with more unknowns than equations

x+y+z=0x+y+z=1

has no solution.

One.I.1.26 §I.1 source ↗

Must any Chemistry problem like the one that starts this subsection — a balance the reaction problem — have infinitely many solutions?

Answer

Yes. For example, the fact that we can have the same reaction in two different flasks shows that twice any solution is another, different, solution (if a physical reaction occurs then there must be at least one nonzero solution).

One.I.1.27 §I.1 source ↗

Find the coefficients a⁠, b⁠, and c so that the graph of f(x)=ax2+bx+c passes through the points (1,2)⁠, (−1,6)⁠, and (2,3)⁠.

Answer

Because f(1)=2⁠, f(−1)=6⁠, and f(2)=3 we get a linear system.

1a+1b+c=21a−1b+c=64a+2b+c=3

Gauss’s Method

⟶−4ρ1+ρ3−ρ1+ρ2a+b+c=2−2b=4−2b−3c=−5⟶−ρ2+ρ3a+b+c=2−2b=4−3c=−9

shows that the solution is f(x)=1x2−2x+3⁠.

One.I.1.28 §I.1 source ↗

After Theorem 1.5 we note that multiplying a row by 0 is not allowed because that could change a solution set. Give an example of a system with solution set S0 where after multiplying a row by 0 the new system has a solution set S1 and S0 is a proper subset of S1⁠, that is, S0≠S1⁠. Give an example where S0=S1⁠.

Answer

Here S0={(1,1)}

x+y=2x−y=0⟶0ρ2x+y=20=0

while S1 is a proper superset because it contains at least two points: (1,1) and (2,0)⁠. In this example the solution set does not change.

x+y=22x+2y=4⟶0ρ2x+y=20=0

One.I.1.29 §I.1 source ↗

Gauss’s Method works by combining the equations in a system to make new equations.

(a)

Can we derive the equation 3x−2y=5 by a sequence of Gaussian reduction steps from the equations in this system?

x+y=14x−y=6

(b)

Can we derive the equation 5x−3y=2 with a sequence of Gaussian reduction steps from the equations in this system?

2x+2y=53x+y=4

(c)

Can we derive 6x−9y+5z=−2 by a sequence of Gaussian reduction steps from the equations in the system?

2x+y−z=46x−3y+z=5

Answer (3 parts)

(a)

Yes, by inspection the given equation results from −ρ1+ρ2⁠.

(b)

No. The pair (1,1) satisfies the given equation. However, that pair does not satisfy the first equation in the system.

(c)

Yes. To see if the given row is c1ρ1+c2ρ2⁠, solve the system of equations relating the coefficients of x⁠, y⁠, z⁠, and the constants:

2c1+6c2=6c1−3c2=−9−c1+c2=54c1+5c2=−2

and get c1=−3 and c2=2⁠, so the given row is −3ρ1+2ρ2⁠.

One.I.1.30 §I.1 source ↗

Prove that, where a,b,c,d,e are real numbers with a≠0⁠, if this linear equation

ax+by=c

has the same solution set as this one

ax+dy=e

then they are the same equation. What if a=0⁠?

Answer

If a≠0 then the solution set of the first equation is this.

{(x,y)∈ℝ2|x=(c−by)/a=(c/a)−(b/a)⋅y}(∗)

Thus, given y we can compute the associated x⁠. Taking y=0 gives the solution (c/a,0)⁠, and since the second equation ax+dy=e is supposed to have the same solution set, substituting into it gives that a(c/a)+d⋅0=e⁠, so c=e⁠. Taking y=1 in (∗⁠) gives a((c−b)/a)+d⋅1=e⁠, and so b=d⁠. Hence they are the same equation.

When a=0 the equations can be different and still have the same solution set: e.g., 0x+3y=6 and 0x+6y=12⁠.

One.I.1.31 §I.1 source ↗

Show that if ad−bc≠0 then

ax+by=jcx+dy=k

has a unique solution.

Answer

We take three cases: that a≠0⁠, that a=0 and c≠0⁠, and that both a=0 and c=0⁠.

For the first, we assume that a≠0⁠. Then the reduction

⟶−(c/a)ρ1+ρ2ax+by=j(−(cb/a)+d)y=−(cj/a)+k

shows that this system has a unique solution if and only if −(cb/a)+d≠0⁠; remember that a≠0 so that back substitution yields a unique x (observe, by the way, that j and k play no role in the conclusion that there is a unique solution, although if there is a unique solution then they contribute to its value). But −(cb/a)+d=(ad−bc)/a and a fraction is not equal to 0 if and only if its numerator is not equal to 0⁠. Thus, in this first case, there is a unique solution if and only if ad−bc≠0⁠.

In the second case, if a=0 but c≠0⁠, then we swap

cx+dy=kby=j

to conclude that the system has a unique solution if and only if b≠0 (we use the case assumption that c≠0 to get a unique x in back substitution). But — where a=0 and c≠0 — the condition “⁠b≠0⁠” is equivalent to the condition “⁠ad−bc≠0⁠”. That finishes the second case.

Finally, for the third case, if both a and c are 0 then the system

0x+by=j0x+dy=k

might have no solutions (if the second equation is not a multiple of the first) or it might have infinitely many solutions (if the second equation is a multiple of the first then for each y satisfying both equations, any pair (x,y) will do), but it never has a unique solution. Note that a=0 and c=0 gives that ad−bc=0⁠.

One.I.1.32 §I.1 source ↗

In the system

ax+by=cdx+ey=f

each of the equations describes a line in the xy⁠-plane. By geometrical reasoning, show that there are three possibilities: there is a unique solution, there is no solution, and there are infinitely many solutions.

Answer

Recall that if a pair of lines share two distinct points then they are the same line. That’s because two points determine a line, so these two points determine each of the two lines, and so they are the same line.

Thus the lines can share one point (giving a unique solution), share no points (giving no solutions), or share at least two points (which makes them the same line).

One.I.1.33 §I.1 source ↗

Finish the proof of Theorem 1.5.

Answer

For the reduction operation of multiplying ρi by a nonzero real number k⁠, we have that (s1,…,sn) satisfies this system

a1,1x1+a1,2x2+⋯+a1,nxn=d1⋮kai,1x1+kai,2x2+⋯+kai,nxn=kdi⋮am,1x1+am,2x2+⋯+am,nxn=dm

if and only if a1,1s1+ a1,2s2+ ⋯+ a1,nsn= d1 and … kai,1s1+ kai,2s2+ ⋯+ kai,nsn= kdi and … am,1s1+ am,2s2+ ⋯+ am,nsn= dm by the definition of ‘satisfies’. Because k≠0⁠, that’s true if and only if a1,1s1+ a1,2s2+ ⋯+ a1,nsn= d1 and … ai,1s1+ ai,2s2+ ⋯+ ai,nsn= di and … am,1s1+ am,2s2+ ⋯+ am,nsn= dm (this is straightforward canceling on both sides of the i⁠-th equation), which says that (s1,…,sn) solves

a1,1x1+a1,2x2+⋯+a1,nxn=d1⋮ai,1x1+ai,2x2+⋯+ai,nxn=di⋮am,1x1+am,2x2+⋯+am,nxn=dm

as required.

For the combination operation kρi+ρj⁠, the tuple (s1,…,sn) satisfies

a1,1x1+⋯+a1,nxn=d1⋮ai,1x1+⋯+ai,nxn=di⋮(kai,1+aj,1)x1+⋯+(kai,n+aj,n)xn=kdi+dj⋮am,1x1+⋯+am,nxn=dm

if and only if a1,1s1+⋯+a1,nsn=d1 and … ai,1s1+⋯+ai,nsn=di and … (kai,1+aj,1)s1+ ⋯+ (kai,n+aj,n)sn= kdi+ dj and … am,1s1+ am,2s2+ ⋯+ am,nsn= dm again by the definition of ‘satisfies’. Subtract k times the equation i from equation j⁠. (Here is where we need i≠j⁠; if i=j then the two di⁠’s above are not equal.) The previous compound statement holds if and only if a1,1s1+⋯+a1,nsn=d1 and … ai,1s1+⋯+ai,nsn=di and… (kai,1+aj,1)s1+ ⋯+ (kai,n+aj,n)sn− (kai,1s1+⋯+kai,nsn)= kdi+ dj− kdi and… am,1s1+⋯+am,nsn=dm⁠, which after cancellation says that (s1,…,sn) solves

a1,1x1+⋯+a1,nxn=d1⋮ai,1x1+⋯+ai,nxn=di⋮aj,1x1+⋯+aj,nxn=dj⋮am,1x1+⋯+am,nxn=dm

as required.

One.I.1.34 §I.1 source ↗

Is there a two-unknowns linear system whose solution set is all of ℝ2⁠?

Answer

Yes, this one-equation system:

0x+0y=0

is satisfied by every (x,y)∈ℝ2⁠.

One.I.1.35 §I.1 source ↗

Are any of the operations used in Gauss’s Method redundant? That is, can we make any of the operations from a combination of the others?

Answer

Yes. This sequence of operations swaps rows i and j

⟶ρi+ρj⟶−ρj+ρi⟶ρi+ρj⟶−1ρi

so the row-swap operation is redundant in the presence of the other two.

One.I.1.36 §I.1 source ↗

Prove that each operation of Gauss’s Method is reversible. That is, show that if two systems are related by a row operation S1→S2 then there is a row operation to go back S2→S1⁠.

Answer

Reverse a row swap ρi↔ρj by swapping back ρj↔ρi⁠. Reverse the kρi step of multiplying k≠0 on both sides of a row by dividing through (1/k)ρi⁠.

The row combination case is the nontrivial one. The operation kρi+ρj results in this j⁠-th row.

k⋅ai,1+aj,1+⋯+k⋅ai,n+aj,n=k⋅di+dj

The i⁠-th row unchanged because of the i≠j restriction. Because the i⁠-th row is unchanged, the operation −kρi+ρj returns the j⁠-th row to its original state.

(Observe that the i=j condition on the kρi+ρj is needed, or else this could happen

3x+2y=7⟶2ρ1+ρ19x+6y=21⟶−2ρ1+ρ1−9x−6y=−21

and so the result wouldn’t hold.)

One.I.1.37 §I.1 source ↗

[Anton] A box holding pennies, nickels and dimes contains thirteen coins with a total value of 83 cents. How many coins of each type are in the box? (These are US coins; a penny is 1⁠ cent, a nickel is 5⁠ cents, and a dime is 10⁠ cents.)

Answer

Let p⁠, n⁠, and d be the number of pennies, nickels, and dimes. For variables that are real numbers, this system

p+n+d=13p+5n+10d=83⟶−ρ1+ρ2p+n+d=134n+9d=70

has more than one solution; in fact, it has infinitely many of them. However, it has a limited number of solutions in which p⁠, n⁠, and d are non-negative integers. Running through d=0⁠, …, d=8 shows that (p,n,d)=(3,4,6) is the only solution using natural numbers.

One.I.1.38 §I.1 source ↗

[Con. Prob. 1955] Four positive integers are given. Select any three of the integers, find their arithmetic average, and add this result to the fourth integer. Thus the numbers 29, 23, 21, and 17 are obtained. One of the original integers is:

(a)

19

(b)

21

(c)

23

(d)

29

(e)

17

Answer

Solving the system

(1/3)(a+b+c)+d=29(1/3)(b+c+d)+a=23(1/3)(c+d+a)+b=21(1/3)(d+a+b)+c=17

we obtain a=12⁠, b=9⁠, c=3⁠, d=21⁠. Thus the second item, 21, is the correct answer.

One.I.1.39 §I.1 source ↗

[Am. Math. Mon., Jan. 1935] Laugh at this: AHAHA+TEHE=TEHAW⁠. It resulted from substituting a code letter for each digit of a simple example in addition, and it is required to identify the letters and prove the solution unique.

Answer

This is how the answer was given in the cited source. A comparison of the units and hundreds columns of this addition shows that there must be a carry from the tens column. The tens column then tells us that A<H⁠, so there can be no carry from the units or hundreds columns. The five columns then give the following five equations.

A+E=W2H=A+10H=W+1H+T=E+10A+1=T

The five linear equations in five unknowns, if solved simultaneously, produce the unique solution: A=4⁠, T=5⁠, H=7⁠, W=6 and E=2⁠, so that the original example in addition was 47474+5272=52746⁠.

One.I.1.40 §I.1 source ↗

[Wohascum no. 2] The Wohascum County Board of Commissioners, which has 20 members, recently had to elect a President. There were three candidates (A⁠, B⁠, and C⁠); on each ballot the three candidates were to be listed in order of preference, with no abstentions. It was found that 11 members, a majority, preferred A over B (thus the other 9 preferred B over A⁠). Similarly, it was found that 12 members preferred C over A⁠. Given these results, it was suggested that B should withdraw, to enable a runoff election between A and C⁠. However, B protested, and it was then found that 14 members preferred B over C⁠! The Board has not yet recovered from the resulting confusion. Given that every possible order of A⁠, B⁠, C appeared on at least one ballot, how many members voted for B as their first choice?

Answer

This is how the answer was given in the cited source. Some additional material was added from [joriki]. Eight commissioners voted for B⁠. To see this, we will use the given information to study how many voters chose each order of A⁠, B⁠, C⁠.

The six orders of preference are ABC⁠, ACB⁠, BAC⁠, BCA⁠, CAB⁠, CBA⁠; assume they receive a⁠, b⁠, c⁠, d⁠, e⁠, f votes respectively. We know that

a+b+e=11d+e+f=12a+c+d=14

from the number preferring A over B⁠, the number preferring C over A⁠, and the number preferring B over C⁠. Because 20 votes were cast, a+b+⋯+f=20⁠. Subtracting the sum of the three above equations from twice the prior equation gives b+c+f=3⁠. We’ve specified that each preference order got at least one vote, so that means b=c=f=1⁠.

From the above three equations the complete solution is then a=6⁠, b=1⁠, c=1⁠, d=7⁠, e=4⁠, and f=1⁠, as we can find with Gauss’s Method. The number of commissioners voting for B as their first choice is therefore c+d=1+7=8⁠.

Comments. The answer to this question would have been the same had we known only that at least 14 commissioners preferred B over C⁠.

The seemingly paradoxical nature of the commissioner’s preferences (A is preferred to B⁠, and B is preferred to C⁠, and C is preferred to A⁠), an example of “non-transitive dominance”, is common when individual choices are pooled.

One.I.1.41 §I.1 source ↗

[Am. Math. Mon., Jan. 1963] “This system of n linear equations with n unknowns,” said the Great Mathematician, “has a curious property.”

“Good heavens!” said the Poor Nut, “What is it?”

“Note,” said the Great Mathematician, “that the constants are in arithmetic progression.”

“It’s all so clear when you explain it!” said the Poor Nut. “Do you mean like 6x+9y=12 and 15x+18y=21⁠?”

“Quite so,” said the Great Mathematician, pulling out his bassoon. “Indeed, the system has a unique solution. Can you find it?”

“Good heavens!” cried the Poor Nut, “I am baffled.”

Are you?

Answer

This is how the answer was given in the cited source. We have not used “dependent” yet; it means here that Gauss’s Method shows that there is not a unique solution. If n≥3 the system is dependent and the solution is not unique. Hence n<3⁠. But the term “system” implies n>1⁠. Hence n=2⁠. If the equations are

ax+(a+d)y=a+2d(a+3d)x+(a+4d)y=a+5d

then x=−1⁠, y=2⁠.

One.I.2.15 §I.2 source ↗

Find the indicated entry of the matrix, if it is defined.

A=(1312−14)

(a)

a2,1

(b)

a1,2

(c)

a2,2

(d)

a3,1

Answer (4 parts)

(a)

2

(b)

3

(c)

−1

(d)

Not defined.

One.I.2.16 §I.2 source ↗

Give the size of each matrix.

(a)

(104215)

(b)

(11−113−1)

(c)

(510105)

Answer (3 parts)

(a)

2×3

(b)

3×2

(c)

2×2

One.I.2.17 §I.2 source ↗

Do the indicated vector operation, if it is defined.

(a)

(211)+(304)

(b)

5(4−1)

(c)

(151)−(311)

(d)

7(21)+9(35)

(e)

(12)+(123)

(f)

6(311)−4(203)+2(115)

Answer (6 parts)

(a)

(515)

(b)

(20−5)

(c)

(−240)

(d)

(4152)

(e)

Not defined.

(f)

(1284)

One.I.2.18 §I.2 source ↗

Solve each system using matrix notation. Express the solution using vectors.

(a)

3x+6y=18x+2y=6

(b)

x+y=1x−y=−1

(c)

x1+x3=4x1−x2+2x3=54x1−x2+5x3=17

(d)

2a+b−c=22a+c=3a−b=0

(e)

x+2y−z=32x+y+w=4x−y+z+w=1

(f)

x+z+w=42x+y−w=23x+y+z=7

Answer (6 parts)

(a)

This reduction

(3618126)⟶(−1/3)ρ1+ρ2(3618000)

leaves x leading and y free. Making y the parameter, gives x=6−2y and this solution set.

{(60)+(−21)y|y∈ℝ}

(b)

A reduction

(1111−1−1)⟶−ρ1+ρ2(1110−2−2)

gives the unique solution y=1⁠, x=0⁠. The solution set is this.

{(01)}

(c)

Gauss’s Method

(10141−1254−1517)⟶−4ρ1+ρ3−ρ1+ρ2(10140−1110−111)⟶−ρ2+ρ3(10140−1110000)

leaves x1 and x2 leading with x3 free. The solution set is this.

{(4−10)+(−111)x3|x3∈ℝ}

(d)

This reduction

(21−1220131−100)⟶−(1/2)ρ1+ρ3−ρ1+ρ2(21−120−1210−3/21/2−1)⟶(−3/2)ρ2+ρ3(21−120−12100−5/2−5/2)

shows that the solution set is a singleton set.

{(111)}

(e)

This reduction is straightforward

(12−103210141−1111)⟶−ρ1+ρ3−2ρ1+ρ2(12−1030−321−20−321−2)⟶−ρ2+ρ3(12−1030−321−200000)

and ends with x and y leading while z and w are free. Solving for y gives y=(2+2z+w)/3 and substitution shows that x+2(2+2z+w)/3−z=3 so x=(5/3)−(1/3)z−(2/3)w⁠, making this the solution set.

{(5/32/300)+(−1/32/310)z+(−2/31/301)w|z,w∈ℝ}

(f)

The reduction

(10114210−1231107)⟶−3ρ1+ρ3−2ρ1+ρ2(1011401−2−3−601−2−3−5)⟶−ρ2+ρ3(1011401−2−3−600001)

shows that there is no solution — the solution set is empty.

One.I.2.19 §I.2 source ↗

Solve each system using matrix notation. Give each solution set in vector notation.

(a)

2x+y−z=14x−y=3

(b)

x−z=1y+2z−w=3x+2y+3z−w=7

(c)

x−y+z=0y+w=03x−2y+3z+w=0−y−w=0

(d)

a+2b+3c+d−e=13a−b+c+d+e=3

Answer (4 parts)

(a)

The reduction

(21−114−103)⟶−2ρ1+ρ2(21−110−321)

ends with x and y leading, and with z free. Solving for y gives y=(1−2z)/(−3)⁠, and then substitution 2x+(1−2z)/(−3)−z=1 shows that x=((4/3)+(1/3)z)/2⁠. Hence the solution set is this.

{(2/3−1/30)+(1/62/31)z|z∈ℝ}

(b)

This application of Gauss’s Method

(10−101012−13123−17)⟶−ρ1+ρ3(10−101012−13024−16)⟶−2ρ2+ρ3(10−101012−1300010)

leaves x⁠, y⁠, and w leading. The solution set is here.

{(1300)+(1−210)z|z∈ℝ}

(c)

This row reduction

(1−1100010103−23100−10−10)⟶−3ρ1+ρ3(1−110001010010100−10−10)⟶ρ2+ρ4−ρ2+ρ3(1−1100010100000000000)

ends with z and w free. We have this solution set.

{(0000)+(−1010)z+(−1−101)w|z,w∈ℝ}

(d)

Gauss’s Method done in this way

(1231−113−11113)⟶−3ρ1+ρ2(1231−110−7−8−240)

ends with c⁠, d⁠, and e free. Solving for b shows that b=(8c+2d−4e)/(−7) and then substitution a+ 2(8c+2d−4e)/(−7)+ 3c+ 1d− 1e= 1 shows that a= 1− (5/7)c− (3/7)d− (1/7)e and we have the solution set.

{(10000)+(−5/7−8/7100)c+(−3/7−2/7010)d+(−1/74/7001)e|c,d,e∈ℝ}

One.I.2.20 §I.2 source ↗

Solve each system using matrix notation. Express the solution set using vectors.

(a)

3x+2y+z=1x−y+z=25x+5y+z=0

(b)

x+y−2z=0x−y=−33x−y−2z=−62y−2z=3

(c)

2x−y−z+w=4x+y+z=−1

(d)

x+y−2z=0x−y=−33x−y−2z=0

Answer (4 parts)

(a)

This reduction

(32111−1125510)⟶−(5/3)ρ1+ρ3−(1/3)ρ1+ρ2(32110−5/32/35/305/3−2/3−5/3)⟶ρ2+ρ3(32110−5/32/35/30000)

gives this solution set.

{(xyz)=(1−10)+(−3/52/51)z|z∈ℝ}

(b)

This is the reduction.

(11−201−1033−1−2−602−23)⟶−3ρ1+ρ3−ρ1+ρ2(11−200−22−30−44−602−23)⟶ρ2+ρ4−2ρ2+ρ3(11−200−22−300000000)

The solution set is this.

{(−3/23/20)+(111)z|z∈ℝ}

(c)

Gauss’s Method

(2−1−1141110−1)⟶−(1/2)ρ1+ρ2(2−1−11403/23/2−1/2−3)

gives the solution set.

{(1−200)+(0−110)z(−1/31/301)w|z,w∈ℝ}

(d)

Here is the reduction.

(11−201−10−33−1−20)⟶−3ρ1+ρ3−ρ1+ρ2(11−200−22−30−440)⟶−2ρ2+ρ3(11−200−22−30006)

The solution set is empty {}⁠.

One.I.2.21 §I.2 source ↗

The vector is in the set. What value of the parameters produces that vector?

(a)

(5−5)⁠, {(1−1)k|k∈ℝ}

(b)

(−121)⁠, {(−210)i+(301)j|i,j∈ℝ}

(c)

(0−42)⁠, {(110)m+(201)n|m,n∈ℝ}

Answer (4 parts)

For each problem we get a system of linear equations by looking at the equations of components.

(a)

k=5

(b)

The second components show that i=2⁠, the third components show that j=1⁠.

(c)

m=−4⁠, n=2

One.I.2.22 §I.2 source ↗

Decide if the vector is in the set.

(a)

(3−1)⁠, {(−62)k|k∈ℝ}

(b)

(54)⁠, {(5−4)j|j∈ℝ}

(c)

(21−1)⁠, {(03−7)+(1−13)r|r∈ℝ}

(d)

(101)⁠, {(201)j+(−3−11)k|j,k∈ℝ}

Answer (5 parts)

For each problem we get a system of linear equations by looking at the equations of components.

(a)

Yes; take k=−1/2⁠.

(b)

No; the system with equations 5=5⋅j and 4=−4⋅j has no solution.

(c)

Yes; take r=2⁠.

(d)

No. The second components give k=0⁠. Then the third components give j=1⁠. But the first components don’t check.

One.I.2.23 §I.2 source ↗

[Cleary] A farmer with 1200 acres is considering planting three different crops, corn, soybeans, and oats. The farmer wants to use all 1200 acres. Seed corn costs $⁠20 per acre, while soybean and oat seed cost $⁠50 and $⁠12 per acre respectively. The farmer has $⁠40000 available to buy seed and intends to spend it all.

(a)

Use the information above to formulate two linear equations with three unknowns and solve it.

(b)

Solutions to the system are choices that the farmer can make. Write down two reasonable solutions.

(c)

Suppose that in the fall when the crops mature, the farmer can bring in revenue of $⁠100 per acre for corn, $⁠300 per acre for soybeans and $⁠80 per acre for oats. Which of your two solutions in the prior part would have resulted in a larger revenue?

Answer (3 parts)

(a)

Let c be the number of acres of corn, s be the number of acres of soy, and a be the number of acres of oats.

c+s+a=120020c+50s+12a=40000⟶−20ρ1+ρ2c+s+a=120030s−8a=16000

To describe the solution set we can parametrize using a⁠.

{(csa)=(2000/31600/30)+(−19/154/151)a|a∈ℝ}

(b)

There are many answers possible here. For instance we can take a=0 to get c=2000/3≈666.67 and s=1600/3≈533.33⁠. Another example is to set c=0 and solve to get a=10000/19≈526.32⁠. Then s=12800/19≈673.68⁠.

(c)

Plugging in the prior part’s first example answers into 100c+300s+80a gives 680000/3⁠, approximately 226666.67⁠. The second example answer gives 100⋅0+ 300⋅12800/19s+ 80⋅10000/19= 4640000/19⁠, about 244210.53⁠.

One.I.2.24 §I.2 source ↗

Parametrize the solution set of this one-equation system.

x1+x2+⋯+xn=0

Answer

This system has one equation. The leading variable is x1⁠, the other variables are free.

{(−11⋮0)x2+⋯+(−10⋮1)xn|x2,…,xn∈ℝ}

One.I.2.25 §I.2 source ↗

(a)

Apply Gauss’s Method to the left-hand side to solve

x+2y−w=a2x+z=bx+y+2w=c

for x⁠, y⁠, z⁠, and w⁠, in terms of the constants a⁠, b⁠, and c⁠.

(b)

Use your answer from the prior part to solve this.

x+2y−w=32x+z=1x+y+2w=−2

Answer (2 parts)

(a)

Gauss’s Method here gives

(120−1a2010b1102c)⟶−ρ1+ρ3−2ρ1+ρ2(120−1a0−412−2a+b0−103−a+c)⟶−(1/4)ρ2+ρ3(120−1a0−412−2a+b00−1/45/2−(1/2)a−(1/4)b+c)

leaving w free. Solve: z=2a+b−4c+10w⁠, and − 4y= − 2a+ b− (2a+b−4c+10w)− 2w so y=a−c+3w⁠, and x= a− 2(a−c+3w)+ w= − a+ 2c− 5w. Therefore the solution set is this.

{(−a+2ca−c2a+b−4c0)+(−53101)w|w∈ℝ}

(b)

Plug in with a=3⁠, b=1⁠, and c=−2⁠.

{(−75150)+(−53101)w|w∈ℝ}

One.I.2.26 §I.2 source ↗

Why is the comma needed in the notation ‘⁠ai,j⁠’ for matrix entries?

Answer

Leaving the comma out, say by writing a123⁠, is ambiguous because it could mean a1,23 or a12,3⁠.

One.I.2.27 §I.2 source ↗

Give the 4×4 matrix whose i,j⁠-th entry is

(a)

i+j⁠;

(b)

−1 to the i+j power.

Answer (2 parts)

(a)

(2345345645675678)

(b)

(1−11−1−11−111−11−1−11−11)

One.I.2.28 §I.2 source ↗

For any matrix A⁠, the transpose of A⁠, written A𝖳⁠, is the matrix whose columns are the rows of A⁠. Find the transpose of each of these.

(a)

(123456)

(b)

(2−311)

(c)

(510105)

(d)

(110)

Answer (4 parts)

(a)

(142536)

(b)

(21−31)

(c)

(510105)

(d)

(110)

One.I.2.29 §I.2 source ↗

(a)

Describe all functions f(x)=ax2+bx+c such that f(1)=2 and f(−1)=6⁠.

(b)

Describe all functions f(x)=ax2+bx+c such that f(1)=2⁠.

Answer (2 parts)

(a)

Plugging in x=1 and x=−1 gives

a+b+c=2a−b+c=6⟶−ρ1+ρ2a+b+c=2−2b=4

so the set of functions is {f(x)=(4−c)x2−2x+c|c∈ℝ}⁠.

(b)

Putting in x=1 gives

a+b+c=2

so the set of functions is {f(x)= (2− b− c)x2+ bx+ c|b, c∈ℝ}⁠.

One.I.2.30 §I.2 source ↗

Show that any set of five points from the plane ℝ2 lie on a common conic section, that is, they all satisfy some equation of the form ax2+by2+cxy+dx+ey+f=0 where some of a,…,f are nonzero.

Answer

Plugging in the five pairs of the form (x,y) gives a system with five equations in the six unknowns a⁠, … f⁠. There are more unknowns than equations so unless there is some inconsistency between the equations, at least one variable will be free. Thus if there is no inconsistency then the system has not just one solution, it has infinitely many.

No inconsistency can exist — these five equations must have at least one solution in common — because each of these five is solved by a=0⁠, … f=0⁠.

One.I.2.31 §I.2 source ↗

Make up a four equations/four unknowns system having

(a)

a one-parameter solution set;

(b)

a two-parameter solution set;

(c)

a three-parameter solution set.

Answer (3 parts)

(a)

Here is one — the fourth equation is redundant but still OK.

x+y−z+w=0y−z=02z+2w=0z+w=0

(b)

Here is one.

x+y−z+w=0w=0w=0w=0

(c)

This is one.

x+y−z+w=0x+y−z+w=0x+y−z+w=0x+y−z+w=0

One.I.2.32 §I.2 source ↗

[Shepelev] This puzzle is from a Russian web-site and there are many solutions to it, but mine uses linear algebra and is very naive. There’s a planet inhabited by arbuzoids (watermeloners, to translate from Russian). Those creatures are found in three colors: red, green and blue. There are 13⁠ red arbuzoids, 15⁠ blue ones, and 17⁠ green. When two differently colored arbuzoids meet, they both change to the third color.

The question is, can it ever happen that all of them assume the same color?

Answer

This is how the answer was given in the cited source. My solution was to define the numbers of arbuzoids as 3⁠-dimensional vectors, and express all possible elementary transitions as such vectors, too.

R:13
G:15
B:17

Operations: (−1−12)⁠, (−12−1)⁠, and (2−1−1)

Now, it is enough to check whether the solution to one of the following systems of linear equations exists:

(131517)+x(−1−12)+y(−12−1)+z(2−1−1)=(0045)(or (0450) or (4500))

Solving

(−1−12−13−12−1−152−1−128)⟶2ρ1+ρ3−ρ1+ρ2⟶ρ2+ρ3(−1−12−1303−3−20000)

gives y+2/3=z so if the number of transformations z is an integer then y is not. The other two systems give similar conclusions so there is no solution.

One.I.2.33 §I.2 source ↗

[USSR Olympiad no. 174]

(a)

Solve the system of equations.

ax+y=a2x+ay=1

For what values of a does the system fail to have solutions, and for what values of a are there infinitely many solutions?

(b)

Answer the above question for the system.

ax+y=a3x+ay=1

Answer (3 parts)

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

(a)

Formal solution of the system yields

x=a3−1a2−1y=−a2+aa2−1.

If a+1≠0 and a−1≠0⁠, then the system has the single solution

x=a2+a+1a+1y=−aa+1.

If a=−1⁠, or if a=+1⁠, then the formulas are meaningless; in the first instance we arrive at the system

{−x+y=1x−y=1

which is a contradictory system. In the second instance we have

{x+y=1x+y=1

which has an infinite number of solutions (for example, for x arbitrary, y=1−x⁠).

(b)

Solution of the system yields

x=a4−1a2−1y=−a3+aa2−1.

Here, is a2−1≠0⁠, the system has the single solution x=a2+1⁠, y=−a⁠. For a=−1 and a=1⁠, we obtain the systems

{−x+y=−1x−y=1{x+y=1x+y=1

both of which have an infinite number of solutions.

One.I.2.34 §I.2 source ↗

[Math. Mag., Sept. 1952] In air a gold-surfaced sphere weighs 7588 grams. It is known that it may contain one or more of the metals aluminum, copper, silver, or lead. When weighed successively under standard conditions in water, benzene, alcohol, and glycerin its respective weights are 6588⁠, 6688⁠, 6778⁠, and 6328 grams. How much, if any, of the forenamed metals does it contain if the specific gravities of the designated substances are taken to be as follows?

Aluminum2.7Alcohol0.81
Copper8.9Benzene0.90
Gold19.3Glycerin1.26
Lead11.3Water1.00
Silver10.8

Answer

This is how the answer was given in the cited source. Let u⁠, v⁠, x⁠, y⁠, z be the volumes in cm3 of Al, Cu, Pb, Ag, and Au, respectively, contained in the sphere, which we assume to be not hollow. Since the loss of weight in water (specific gravity 1.00⁠) is 1000 grams, the volume of the sphere is 1000 cm3. Then the data, some of which is superfluous, though consistent, leads to only 2 independent equations, one relating volumes and the other, weights.

u+v+x+y+z=10002.7u+8.9v+11.3x+10.5y+19.3z=7558

Clearly the sphere must contain some aluminum to bring its mean specific gravity below the specific gravities of all the other metals. There is no unique result to this part of the problem, for the amounts of three metals may be chosen arbitrarily, provided that the choices will not result in negative amounts of any metal.

If the ball contains only aluminum and gold, there are 294.5 cm3 of gold and 705.5 cm3 of aluminum. Another possibility is 124.7 cm3 each of Cu, Au, Pb, and Ag and 501.2 cm3 of Al.

One.I.3.14 §I.3 source ↗

Solve this system. Then solve the associated homogeneous system.

x+y−2z=0x−y=−33x−y−2z=−62y−2z=3

Answer

This reduction solves the system.

(11−201−1033−1−2−602−23)⟶−3ρ1+ρ3−ρ1+ρ2(11−200−22−30−44−602−23)⟶ρ2+ρ4−2ρ2+ρ3(11−200−22−300000000)

The solution set is this.

{(−3/23/20)+(111)z|z∈ℝ}

Similarly we can reduce the associated homogeneous system

(11−201−1003−1−2002−20)⟶−3ρ1+ρ3−ρ1+ρ2(11−200−2200−44002−20)⟶ρ2+ρ4−2ρ2+ρ3(11−200−22000000000)

to get its solution set.

{(111)z|z∈ℝ}

One.I.3.15 §I.3 source ↗

Solve each system. Express the solution set using vectors. Identify a particular solution and the solution set of the homogeneous system.

(a)

3x+6y=18x+2y=6

(b)

x+y=1x−y=−1

(c)

x1+x3=4x1−x2+2x3=54x1−x2+5x3=17

(d)

2a+b−c=22a+c=3a−b=0

(e)

x+2y−z=32x+y+w=4x−y+z+w=1

(f)

x+z+w=42x+y−w=23x+y+z=7

Answer (6 parts)

(a)

Here is the Gauss-Jordan reduction.

(3618126)⟶(−1/3)ρ1+ρ2(3618000)⟶(1/3)ρ1(126000)

This is the solution set.

S={(60)+(−21)y|y∈ℝ}

Here is a particular solution, and the solution set to the associated homogeneous system.

(60)and{(−21)y|y∈ℝ}

(As with the solution set description, there is more than one correct answer for a particular solution. For instance, in the description given above of S⁠, setting the parameter y to be 2 gives a vector with first component 2 and second component also 2⁠.)

(b)

The reduction

(1111−1−1)⟶−ρ1+ρ2(1110−2−2)⟶(−1/2)ρ2(111011)⟶−ρ2+ρ1(100011)

gives that this is the solution set.

{(01)}

Here is a particular solution, and the solution set for the associated homogeneous system.

(01){(00)}

(c)

The solution set is infinite.

{(4−10)+(−111)x3|x3∈ℝ}

This is a particular solution and the solution set for the associated homogeneous system.

(4−10){(−111)x3|x3∈ℝ}

(d)

The solution set is a singleton.

{(111)}

A particular solution and the solution set for the associated homogeneous system are here.

(111){(000)}

(e)

The solution set is infinite.

{(5/32/300)+(−1/32/310)z+(−2/31/301)w|z,w∈ℝ}

A particular solution and the solution set for the associated homogeneous system are here.

(5/32/300){(−1/32/310)z+(−2/31/301)w|z,w∈ℝ}

(f)

This system’s solution set is empty. Thus, there is no particular solution. The solution set of the associated homogeneous system is this.

{(−1210)z+(−1301)w|z,w∈ℝ}

One.I.3.16 §I.3 source ↗

Solve each system, giving the solution set in vector notation. Identify a particular solution and the solution of the homogeneous system.

(a)

2x+y−z=14x−y=3

(b)

x−z=1y+2z−w=3x+2y+3z−w=7

(c)

x−y+z=0y+w=03x−2y+3z+w=0−y−w=0

(d)

a+2b+3c+d−e=13a−b+c+d+e=3

Answer (5 parts)

The answers from the prior subsection show the row operations. Each answer here just lists the solution set, the particular solution, and the homogeneous solution.

(a)

The solution set is this.

{(2/3−1/30)+(1/62/31)z|z∈ℝ}

A particular solution and the solution set for the associated homogeneous system are here.

(2/3−1/30){(1/62/31)z|z∈ℝ}

(b)

The solution set is infinite.

{(1300)+(1−210)z|z∈ℝ}

Here are a particular solution and the solution set for the associated homogeneous system.

(1300){(1−210)z|z∈ℝ}

(c)

This is the solution set.

{(0000)+(−1010)z+(−1−101)w|z,w∈ℝ}

Here is a particular solution and the solution set for the associated homogeneous system.

(0000){(−1010)z+(−1−101)w|z,w∈ℝ}

(d)

The solution set is this.

{(10000)+(−5/7−8/7100)c+(−3/7−2/7010)d+(−1/74/7001)e|c,d,e∈ℝ}

And, this is a particular solution and the solution set for the associated homogeneous system.

(10000){(−5/7−8/7100)c+(−3/7−2/7010)d+(−1/74/7001)e|c,d,e∈ℝ}

One.I.3.17 §I.3 source ↗

For the system

2x−y−w=3y+z+2w=2x−2y−z=−1

which of these can be used as the particular solution part of some general solution?

(a)

(0−350)

(b)

(2110)

(c)

(−1−48−1)

Answer (4 parts)

Just plug them in and see if they satisfy all three equations.

(a)

No.

(b)

Yes.

(c)

Yes.

One.I.3.18 §I.3 source ↗

Lemma 3.7 says that we can use any particular solution for p→⁠. Find, if possible, a general solution to this system

x−y+w=42x+3y−z=0y+z+w=4

that uses the given vector as its particular solution.

(a)

(0004)

(b)

(−51−710)

(c)

(2−111)

Answer (4 parts)

Gauss’s Method on the associated homogeneous system

(1−101023−10001110)⟶−2ρ1+ρ2(1−101005−1−2001110)⟶−(1/5)ρ2+ρ3(1−101005−1−20006/57/50)

gives this is the solution to the homogeneous problem.

{(−5/61/6−7/61)w|w∈ℝ}

(a)

That vector is indeed a particular solution, so the required general solution is this.

{(0004)+(−5/61/6−7/61)w|w∈ℝ}

(b)

That vector is a particular solution so the required general solution is this.

{(−51−710)+(−5/61/6−7/61)w|w∈ℝ}

(c)

That vector is not a solution of the system since it does not satisfy the third equation. No such general solution exists.

One.I.3.19 §I.3 source ↗

One is nonsingular while the other is singular. Which is which?

(a)

(134−12)

(b)

(13412)

Answer

The first is nonsingular while the second is singular. Just do Gauss’s Method and see if the echelon form result has non-⁠0 numbers in each entry on the diagonal.

One.I.3.20 §I.3 source ↗

Singular or nonsingular?

(a)

(1213)

(b)

(12−3−6)

(c)

(121131)

(d)

(121113347)

(e)

(221105−114)

Answer (5 parts)

(a)

Nonsingular:

⟶−ρ1+ρ2(1201)

ends with each row containing a leading entry.

(b)

Singular:

⟶3ρ1+ρ2(1200)

ends with row 2 without a leading entry.

(c)

Neither. A matrix must be square for either word to apply.

(d)

Singular.

(e)

Nonsingular.

One.I.3.21 §I.3 source ↗

Is the given vector in the set generated by the given set?

(a)

(23), {(14),(15)}

(b)

(−101), {(210),(101)}

(c)

(130), {(104), (215), (330), (421)}

(d)

(1011), {(2101),(3002)}

Answer (5 parts)

In each case we must decide if the vector is a linear combination of the vectors in the set.

(a)

Yes. Solve

c1(14)+c2(15)=(23)

with

(112453)⟶−4ρ1+ρ2(11201−5)

to conclude that there are c1 and c2 giving the combination.

(b)

No. The reduction

(21−1100011)⟶−(1/2)ρ1+ρ2(21−10−1/21/2011)⟶2ρ2+ρ3(21−10−1/21/2002)

shows that

c1(210)+c2(101)=(−101)

has no solution.

(c)

Yes. The reduction

(123410132345010)⟶−4ρ1+ρ3(12341013230−3−12−15−4)⟶3ρ2+ρ3(123410132300−3−95)

shows that there are infinitely many ways

{(c1c2c3c4)=(−108−5/30)+(−97−31)c4|c4∈ℝ}

to write a combination.

(130)=c1(104)+c2(215)+c3(330)+c4(421)

(d)

No. Look at the third components.

One.I.3.22 §I.3 source ↗

Prove that any linear system with a nonsingular matrix of coefficients has a solution, and that the solution is unique.

Answer

Because the matrix of coefficients is nonsingular, Gauss’s Method ends with an echelon form where each variable leads an equation. Back substitution gives a unique solution.

(Another way to see that the solution is unique is to note that with a nonsingular matrix of coefficients the associated homogeneous system has a unique solution, by definition. Since the general solution is the sum of a particular solution with each homogeneous solution, the general solution has at most one element.)

One.I.3.23 §I.3 source ↗

In the proof of Lemma 3.6, what happens if there are no non-⁠0=0 equations?

Answer

In this case the solution set is all of ℝn and we can express it in the required form.

{c1(10⋮0)+c2(01⋮0)+⋯+cn(00⋮1)|c1,…,cn∈ℝ}

One.I.3.24 §I.3 source ↗

Prove that if s→ and t→ satisfy a homogeneous system then so do these vectors.

(a)

s→+t→

(b)

3s→

(c)

ks→+mt→ for k,m∈ℝ

What’s wrong with this argument: “These three show that if a homogeneous system has one solution then it has many solutions — any multiple of a solution is another solution, and any sum of solutions is a solution also — so there are no homogeneous systems with exactly one solution.”?

Answer (5 parts)

Assume s→,t→∈ℝn and write them as here.

s→=(s1⋮sn)t→=(t1⋮tn)

Also let ai,1x1+⋯+ai,nxn=0 be the i⁠-th equation in the homogeneous system.

(a)

The check is easy.

ai,1(s1+t1)+⋯+ai,n(sn+tn)=(ai,1s1+⋯+ai,nsn)+(ai,1t1+⋯+ai,ntn)=0+0

(b)

This is similar to the prior one.

ai,1(3s1)+⋯+ai,n(3sn)=3(ai,1s1+⋯+ai,nsn)=3⋅0=0

(c)

This one is not much harder.

ai,1(ks1+mt1)+⋯+ai,n(ksn+mtn)=k(ai,1s1+⋯+ai,nsn)+m(ai,1t1+⋯+ai,ntn)=k⋅0+m⋅0

What is wrong with that argument is that any linear combination involving only the zero vector yields the zero vector.

One.I.3.25 §I.3 source ↗

Prove that if a system with only rational coefficients and constants has a solution then it has at least one all-rational solution. Must it have infinitely many?

Answer

First the proof.

If the system has all rational coefficients then Gauss’s Method with back substitution will use only rationals (e.g., −(m/n)ρi+ρj⁠). Thus we can express the solution set using only rational numbers as the components of each vector. Therefore, the particular solution is all rational.

As to the question about infinitely many, there are infinitely many rational vector solutions if and only if the associated homogeneous system has infinitely many real vector solutions. That’s because setting any parameters to be rationals will produce an all-rational solution.

One.I.4.3 §I.4 source ↗

Decide if the vector is a member of the set.

(a)

(23)⁠, {(12)k|k∈ℝ}

(b)

(−33)⁠, {(1−1)k|k∈ℝ}

(c)

(−334)⁠, {(1−12)k|k∈ℝ}

(d)

(−334)⁠, {(1−12)k+(002)m|k,m∈ℝ}

(e)

(1414)⁠, {(225)k+(−102)m|k,m∈ℝ}

(f)

(146)⁠, {(225)k+(−102)m|k,m∈ℝ}

Answer (6 parts)

(a)

No.

(b)

Yes.

(c)

No.

(d)

Yes.

(e)

Yes; use Gauss’s Method to get k=4 and m=−3⁠.

(f)

No; use Gauss’s Method to conclude that there is no solution.

One.I.4.4 §I.4 source ↗

Produce two descriptions of this set that are different than this one.

{(2−5)k|k∈ℝ}

Answer

One easy thing to do is to double and triple the vector:

{(4−10)k|k∈ℝ}and{(6−55)k|k∈ℝ}.

One.I.4.5 §I.4 source ↗

Show that the three descriptions given at the start of this subsection all describe the same set.

Answer

Instead of showing all three equalities, we can show that the first equals the second, and that the second equals the third. Both equalities are easy, using the methods of this subsection.

One.I.4.6 §I.4 source ↗

Show that these sets are equal

{(1411)+(−1010)z|z∈ℝ}and{(0421)+(−1010)k|k∈ℝ},

and that both describe the solution set of this system.

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

Answer

That system reduces like this:

⟶−ρ1+ρ2x−y+z+w=−1y−w=3y+w=5⟶−ρ2+ρ3x−y+z+w=−1y−w=32w=2

showing that w=1⁠, y=4 and x=2−z⁠.

One.I.4.7 §I.4 source ↗

Decide if the sets are equal.

(a)

{(12)+(03)t|t∈ℝ} and {(18)+(0−1)s|s∈ℝ}

(b)

{(131)t+(215)s|t,s∈ℝ} and {(477)m+ (−4−2−10)n|m, n∈ℝ}

(c)

{(12)t|t∈ℝ} and {(24)m+(48)n|m,n∈ℝ}

(d)

{(102)s+(−110)t|s,t∈ℝ} and {(−111)m+(013)n|m,n∈ℝ}

(e)

{(131)t+(246)s|t,s∈ℝ} and {(377)t+(131)s|t,s∈ℝ}

Answer (6 parts)

For each item, we call the first set S1 and the other S2⁠.

(a)

They are equal.

To see that S1⊆S2⁠, we must show that any element of the first set is in the second, that is, for any vector of the form

v→=(12)+(03)t

there is an appropriate s such that

v→=(18)+(0−1)s.

Restated, given t we must find s so that this holds.

1+0s=1+0t8−1s=2+3t

That system reduces to

1=1s=6−3t

That is,

(12)+(03)t=(18)+(0−1)(6−3t)

and so we can state any vector in the form for S1 can also in the form needed for inclusion in S2⁠.

For S2⊆S1⁠, we look for t so that these equations hold.

1+0t=1+0s2+3t=8−1s

Rewrite that as

1=1t=2−(1/3)s

and so

(18)+(0−1)s=(12)+(03)(2−(1/3)s).

(b)

These two are equal.

To show that S1⊆S2⁠, we check that for any t,s we can find an appropriate m,n so that these hold.

4m−4n=1t+2s7m−2n=3t+1s7m−10n=1t+5s

Use Gauss’s Method

(4−41t+2s7−23t+1s7−101t+5s)⟶(−7/4)ρ1+ρ3(−7/4)ρ1+ρ2(4−41t+2s05(5/4)t−(10/4)s0−3−(3/4)t+(6/4)s)⟶(3/5)ρ2+ρ3(4−41t+2s05(5/4)t−(10/4)s000)

to conclude that

(131)t+(215)s=(477)((1/2)t)+(−4−2−10)((1/4)t−(1/2)s)

and so S1⊆S2⁠.

For S2⊆S1⁠, solve

1t+2s=4m−4n3t+1s=7m−2n1t+5s=7m−10n

with Gaussian reduction

(124m−4n317m−2n157m−10n)⟶−ρ1+ρ3−3ρ1+ρ2(124m−4n0−5−5m+10n033m−6n)⟶(3/5)ρ2+ρ3(124m−4n0−5−5m+10n000)

to get

(477)m+(−4−2−10)n=(131)(2m)+(215)(m−2n)

and so we can express any member of S2 in the form needed for S1⁠.

(c)

These sets are equal.

To prove that S1⊆S2⁠, we must be able to solve

2m+4n=1t4m+8n=2t

for m and n in terms of t⁠. Apply Gaussian reduction

(241t482t)⟶−2ρ1+ρ2(241t000)

to conclude that any pair m,n where 2m+4n=t will do. For instance,

(12)t=(24)((1/2)t)+(48)(0)

or

(12)t=(24)((−3/2)t)+(48)(t).

Thus S1⊆S2⁠.

For S2⊆S1⁠, we solve

1t=2m+4n2t=4m+8n

with Gauss’s Method

(12m+4n24m+8n)⟶−2ρ1+ρ2(12m+4n00)

to deduce that any vector in S2 is also in S1⁠.

(24)m+(48)n=(12)(2m+4n)∈S1.

(d)

Neither set is a subset of the other.

For S1⊆S2 to hold we must be able to solve

−1m+0n=1s−1t1m+1n=0s+1t1m+3n=2s+0t

for m and n in terms of t and s⁠. Gauss’s Method

(−101s−1t110s+1t132s+0t)⟶ρ1+ρ3ρ1+ρ2(−101s−1t011s+0t033s−1t)⟶−3ρ2+ρ3(−101s−1t011s+0t030s−1t)

shows that we can only find an appropriate pair m,n when t=0⁠. That is,

(−110)

has no expression of the form

(−111)m+(013)n.

Having shown that S1 is not a subset of S2⁠, we know S1≠S2 so, strictly speaking, we need not go further. But we shall also show that S2 is not a subset of S1⁠.

For S2⊆S1 to hold, we must be able to solve

1s−1t=−1m+0n0s+1t=1m+1n2s+0t=1m+3n

for s and t⁠. Apply row reduction

(1−1−1m+0n011m+1n201m+3n)⟶−2ρ1+ρ3(1−1−1m+0n011m+1n023m+3n)⟶−2ρ2+ρ3(1−1−1m+0n011m+1n001m+1n)

to deduce that the only vectors from S2 that are also in S1 are of the form

(−111)m+(013)(−m).

For instance,

(−111)

is in S2 but not in S1⁠.

(e)

These sets are equal.

First we change the parameters:

S2={(377)m+(131)n|m,n∈ℝ}.

Now, to show that S1⊆S2⁠, we solve

3m+1n=1t+2s7m+3n=3t+4s7m+1n=1t+6s

with Gauss’s Method

(311t+2s733t+4s711t+6s)⟶(−7/3)ρ1+ρ3(−7/3)ρ1+ρ2(311t+2s02/3(2/3)t−(2/3)s0−4/3(−4/3)t+(4/3)s)⟶2ρ2+ρ3(311t+2s02/3(2/3)t−(2/3)s000)

to get that

(131)t+(246)s=(377)(s)+(131)(t−s)

and so S1⊆S2⁠.

The proof that S2⊆S1 involves solving

1t+2s=3m+1n3t+4s=7m+3n1t+6s=7m+1n

with Gaussian reduction

(123m+1n347m+3n167m+1n)⟶−ρ1+ρ3−3ρ1+ρ2(123m+1n0−2−2m044m)⟶2ρ2+ρ3(123m+1n0−2−2m000)

to conclude

(377)m+(131)n=(131)(m+n)+(246)(m)

and so any vector in S2 is also in S1⁠.

One.II.1.1 §II.1 source ↗

Find the canonical name for each vector.

(a)

the vector from (2,1) to (4,2) in ℝ2

(b)

the vector from (3,3) to (2,5) in ℝ2

(c)

the vector from (1,0,6) to (5,0,3) in ℝ3

(d)

the vector from (6,8,8) to (6,8,8) in ℝ3

Answer (4 parts)

(a)

(21)

(b)

(−12)

(c)

(40−3)

(d)

(000)

One.II.1.2 §II.1 source ↗

Decide if the two vectors are equal.

(a)

the vector from (5,3) to (6,2) and the vector from (1,−2) to (1,1)

(b)

the vector from (2,1,1) to (3,0,4) and the vector from (5,1,4) to (6,0,7)

Answer (2 parts)

(a)

No, their canonical positions are different.

(1−1)(03)

(b)

Yes, their canonical positions are the same.

(1−13)

One.II.1.3 §II.1 source ↗

Does (1,0,2,1) lie on the line through (−2,1,1,0) and (5,10,−1,4)⁠?

Answer

That line is this set.

{(−2110)+(79−24)t|t∈ℝ}

Note that this system

−2+7t=11+9t=01−2t=20+4t=1

has no solution. Thus the given point is not in the line.

One.II.1.4 §II.1 source ↗

(a)

Describe the plane through (1,1,5,−1)⁠, (2,2,2,0)⁠, and (3,1,0,4)⁠.

(b)

Is the origin in that plane?

Answer (2 parts)

(a)

Note that

(2220)−(115−1)=(11−31)(3104)−(115−1)=(20−55)

and so the plane is this set.

{(115−1)+(11−31)t+(20−55)s|t,s∈ℝ}

(b)

No; this system

1+1t+2s=01+1t=05−3t−5s=0−1+1t+5s=0

has no solution.

One.II.1.5 §II.1 source ↗

Give a vector description of each.

(a)

the plane subset of ℝ3 with equation x−2y+z=4

(b)

the plane in ℝ3 with equation 2x+y+4z=−1

(c)

the hyperplane subset of ℝ4 with equation x+y+z+w=10

Answer (3 parts)

(a)

Think of x−2y+z=4 as a one-equation linear system and parametrize with the variables y and z to get x=4+2y−z⁠. That gives this vector description of the plane.

{(xyz)=(400)+(210)⋅y+(−101)⋅z|y,z∈ℝ}

(b)

Parametrizing gives x=−(1/2)−(1/2)y−2z⁠, so this is the vector description.

(xyz)=(−1/200)+(−1/210)⋅y+(−201)⋅z

(c)

Here x=10−y−z−w and so we get a vector description with three parameters.

(xyzw)=(10000)+(−1100)⋅y+(−1010)⋅z+(−1001)⋅w

One.II.1.6 §II.1 source ↗

Describe the plane that contains this point and line.

(203){(−10−4)+(112)t|t∈ℝ}

Answer

The vector

(203)

is not in the line. Because

(203)−(−10−4)=(307)

we can describe that plane in this way.

{(−10−4)+m(112)+n(307)|m,n∈ℝ}

One.II.1.7 §II.1 source ↗

Intersect these planes.

{(111)t+(013)s|t,s∈ℝ}{(110)+(030)k+(204)m|k,m∈ℝ}

Answer

The points of coincidence are solutions of this system.

t=1+2mt+s=1+3kt+3s=4m

Move the unknowns to the left and do Gauss’s Method.

(100−2111−301130−40)⟶−ρ1+ρ3−ρ1+ρ2(100−2101−320030−2−1)⟶−3ρ2+ρ3(100−2101−320009−8−1)

The result is k=−(1/9)+(8/9)m⁠, so s=−(1/3)+(2/3)m and t=1+2m⁠. The intersection is this.

{(110)+(030)(−19+89m)+(204)m|m∈ℝ}={(12/30)+(28/34)m|m∈ℝ}

One.II.1.8 §II.1 source ↗

Intersect each pair, if possible.

(a)

{(112)+t(011)|t∈ℝ}⁠,  {(13−2)+s(012)|s∈ℝ}

(b)

{(201)+t(11−1)|t∈ℝ}⁠,  {s(012)+w(041)|s,w∈ℝ}

Answer (2 parts)

(a)

The system

1=11+t=3+s2+t=−2+2s

gives s=6 and t=8⁠, so this is the solution set.

{(1910)}

(b)

This system

2+t=0t=s+4w1−t=2s+w

gives t=−2⁠, w=−1⁠, and s=2 so their intersection is this point.

(0−23)

One.II.1.9 §II.1 source ↗

How should we define ℝ0⁠?

Answer

We shall later define it to be a set with one element — an “origin”.

One.II.1.10 §II.1 source ↗

[Math. Mag., Jan. 1957] A person traveling eastward at a rate of 3 miles per hour finds that the wind appears to blow directly from the north. On doubling his speed it appears to come from the north east. What was the wind’s velocity?

Answer

This is how the answer was given in the cited source. The vector triangle is as follows, so w→=32 from the north west.

One.II.1.11 §II.1 source ↗

Euclid describes a plane as “a surface which lies evenly with the straight lines on itself”. Commentators such as Heron have interpreted this to mean, “(A plane surface is) such that, if a straight line pass through two points on it, the line coincides wholly with it at every spot, all ways”. (Translations from [Heath], pp. 171-172.) Do planes, as described in this section, have that property? Does this description adequately define planes?

Answer

Euclid no doubt is picturing a plane inside of ℝ3. Observe, however, that both ℝ1 and ℝ2 also satisfy that definition.

One.II.2.11 §II.2 source ↗

Find the length of each vector.

(a)

(31)

(b)

(−12)

(c)

(411)

(d)

(000)

(e)

(1−110)

Answer (5 parts)

(a)

32+12=10

(b)

5

(c)

18

(d)

0

(e)

3

One.II.2.12 §II.2 source ↗

Find the angle between each two, if it is defined.

(a)

(12),(14)

(b)

(120),(041)

(c)

(12),(14−1)

Answer (3 parts)

(a)

arccos⁡(9/85)≈ 0.22 radians

(b)

arccos⁡(8/85)≈ 0.52 radians

(c)

Not defined.

One.II.2.13 §II.2 source ↗

[Ohanian] During maneuvers preceding the Battle of Jutland, the British battle cruiser Lion moved as follows (in nautical miles): 1.2 miles north, 6.1 miles 38 degrees east of south, 4.0 miles at 89 degrees east of north, and 6.5 miles at 31 degrees east of north. Find the distance between starting and ending positions. (Ignore the earth’s curvature.)

Answer

We express each displacement as a vector, rounded to one decimal place because that’s the accuracy of the problem’s statement, and add to find the total displacement (ignoring the curvature of the earth).

(0.01.2)+(3.8−4.8)+(4.00.1)+(3.35.6)=(11.12.1)

The distance is 11.12+2.12≈11.3⁠.

One.II.2.14 §II.2 source ↗

Find k so that these two vectors are perpendicular.

(k1)(43)

Answer

Solve (k)(4)+(1)(3)=0 to get k=−3/4⁠.

One.II.2.15 §II.2 source ↗

Describe the set of vectors in ℝ3 orthogonal to the one with entries 1⁠, 3⁠, and −1⁠.

Answer

We could describe the set

{(xyz)|1x+3y−1z=0}

with parameters in this way.

{(−310)y+(101)z|y,z∈ℝ}

One.II.2.16 §II.2 source ↗

(a)

Find the angle between the diagonal of the unit square in ℝ2 and any one of the axes.

(b)

Find the angle between the diagonal of the unit cube in ℝ3 and one of the axes.

(c)

Find the angle between the diagonal of the unit cube in ℝn and one of the axes.

(d)

What is the limit, as n goes to ∞⁠, of the angle between the diagonal of the unit cube in ℝn and any one of the axes?

Answer (4 parts)

(a)

We can use the x⁠-axis.

arccos⁡((1)(1)+(0)(1)12)≈0.79 radians

(b)

Again, use the x⁠-axis.

arccos⁡((1)(1)+(0)(1)+(0)(1)13)≈0.96 radians

(c)

The x⁠-axis worked before and it will work again.

arccos⁡((1)(1)+⋯+(0)(1)1n)=arccos⁡(1n)

(d)

Using the formula from the prior item, limn→∞⁡arccos⁡(1/n)= π/2 radians⁠.

One.II.2.17 §II.2 source ↗

Is there any vector in ℝn that is perpendicular to itself?

Answer

Definition 2.7 stipulates that the angle between two vectors is a right angle if either is the zero vector. (It stipulates that separately since we cannot rely on the formula because the fraction’s denominator is 0⁠.) Thus, by definition, the zero vector is perpendicular to itself.

It is the only such member of ℝn since u1u1+⋯+unun is 0 if and only if each ui is 0⁠.

One.II.2.18 §II.2 source ↗

Describe the algebraic properties of dot product.

(a)

Is it right-distributive over addition: (u→+v→)∙w→=u→∙w→+v→∙w→⁠?

(b)

Is it left-distributive (over addition)?

(c)

Does it commute?

(d)

Associate?

(e)

How does it interact with scalar multiplication?

As always, you must back any assertion with a suitable argument.

Answer (6 parts)

In each item below, assume that the vectors u→,v→,w→∈ℝn have components u1,…,un,v1,…,wn⁠.

(a)

Dot product is right-distributive.

(u→+v→)∙w→=[(u1⋮un)+(v1⋮vn)]∙(w1⋮wn)=(u1+v1⋮un+vn)∙(w1⋮wn)=(u1+v1)w1+⋯+(un+vn)wn=(u1w1+⋯+unwn)+(v1w1+⋯+vnwn)=u→∙w→+v→∙w→

(b)

Dot product is also left distributive: w→∙(u→+v→)=w→∙u→+w→∙v→⁠. The proof is just like the prior one.

(c)

Dot product commutes.

(u1⋮un)∙(v1⋮vn)=u1v1+⋯+unvn=v1u1+⋯+vnun=(v1⋮vn)∙(u1⋮un)

(d)

Because u→∙v→ is a scalar, not a vector, the expression (u→∙v→)∙w→ makes no sense; the dot product of a scalar and a vector is not defined.

(e)

This is a vague question so it has many answers. Some are (1) k(u→∙v→)=(ku→)∙v→ and k(u→∙v→)=u→∙(kv→)⁠, (2) k(u→∙v→)≠(ku→)∙(kv→) (in general; an example is easy to produce), and (3) |kv→|=|k||v→| (the connection between length and dot product is that the square of the length is the dot product of a vector with itself).

One.II.2.19 §II.2 source ↗

Verify the equality condition in Corollary 2.6, the Cauchy-Schwarz Inequality.

(a)

Show that if u→ is a negative scalar multiple of v→ then u→∙v→ and v→∙u→ are less than or equal to zero.

(b)

Show that |u→∙v→|=|u→||v→| if and only if one vector is a scalar multiple of the other.

Answer (2 parts)

(a)

Verifying that (kx→)∙y→= k(x→∙y→)= x→∙(ky→) for k∈ℝ and x→,y→∈ℝn is easy. Now, for k∈ℝ and v→,w→∈ℝn⁠, if u→=kv→ then u→∙v→=(kv→)∙v→=k(v→∙v→)⁠, which is k times a nonnegative real.

The v→=ku→ half is similar (actually, taking the k in this paragraph to be the reciprocal of the k above gives that we need only worry about the k=0 case).

(b)

We first consider the u→∙v→≥0 case. From the Triangle Inequality we know that u→∙v→=|u→||v→| if and only if one vector is a nonnegative scalar multiple of the other. But that’s all we need because the first part of this exercise shows that, in a context where the dot product of the two vectors is positive, the two statements ‘one vector is a scalar multiple of the other’ and ‘one vector is a nonnegative scalar multiple of the other’, are equivalent.

We finish by considering the u→∙v→<0 case. Because 0< |u→∙v→|= − (u→∙v→)= (−u→)∙v→ and |u→||v→|=|−u→||v→|⁠, we have that 0<(−u→)∙v→=|−u→||v→|⁠. Now the prior paragraph applies to give that one of the two vectors −u→ and v→ is a scalar multiple of the other. But that’s equivalent to the assertion that one of the two vectors u→ and v→ is a scalar multiple of the other, as desired.

One.II.2.20 §II.2 source ↗

Suppose that u→∙v→=u→∙w→ and u→≠0→⁠. Must v→=w→⁠?

Answer

No. These give an example.

u→=(10)v→=(10)w→=(11)

One.II.2.21 §II.2 source ↗

Does any vector have length zero except a zero vector? (If “yes”, produce an example. If “no”, prove it.)

Answer

We prove that a vector has length zero if and only if all its components are zero.

Let u→∈ℝn have components u1,…,un⁠. Recall that the square of any real number is greater than or equal to zero, with equality only when that real is zero. Thus |u→|2=u12+⋯+un2 is a sum of numbers greater than or equal to zero, and so is itself greater than or equal to zero, with equality if and only if each ui is zero. Hence |u→|=0 if and only if all the components of u→ are zero.

One.II.2.22 §II.2 source ↗

Find the midpoint of the line segment connecting (x1,y1) with (x2,y2) in ℝ2⁠. Generalize to ℝn⁠.

Answer

We can easily check that

(x1+x22,y1+y22)

is on the line connecting the two, and is equidistant from both. The generalization is obvious.

One.II.2.23 §II.2 source ↗

Show that if v→≠0→ then v→/|v→| has length one. What if v→=0→⁠?

Answer

Assume that v→∈ℝn has components v1,…,vn⁠. If v→≠0→ then we have this.

(v1v12+⋯+vn2)2+⋯+(vnv12+⋯+vn2)2=(v12v12+⋯+vn2)+⋯+(vn2v12+⋯+vn2)=1

If v→=0→ then v→/|v→| is not defined.

One.II.2.24 §II.2 source ↗

Show that if r≥0 then rv→ is r times as long as v→⁠. What if r<0⁠?

Answer

For the first question, assume that v→∈ℝn and r≥0⁠, take the root, and factor.

|rv→|=(rv1)2+⋯+(rvn)2=r2(v12+⋯+vn2=r|v→|

For the second question, the result is r times as long, but it points in the opposite direction in that rv→+(−r)v→=0→⁠.

One.II.2.25 §II.2 source ↗

A vector v→∈ℝn of length one is a unit vector. Show that the dot product of two unit vectors has absolute value less than or equal to one. Can ‘less than’ happen? Can ‘equal to’?

Answer

Assume that u→,v→∈ℝn both have length 1⁠. Apply Cauchy-Schwarz: |u→∙v→|≤|u→||v→|=1⁠.

To see that ‘less than’ can happen, in ℝ2 take

u→=(10)v→=(01)

and note that u→∙v→=0⁠. For ‘equal to’, note that u→∙u→=1⁠.

One.II.2.26 §II.2 source ↗

When a plane does not pass through the origin, performing operations on vectors whose bodies lie in it is more complicated than when the plane does pass through the origin. Consider the picture in this subsection of the plane

{(200)+(−0.510)y+(−0.501)z|y,z∈ℝ}

and the three vectors with endpoints (2,0,0)⁠, (1.5,1,0)⁠, and (1.5,0,1)⁠.

(a)

Redraw the picture, including the vector starting at (2,0,0) whose body is in the plane, and that is twice as long as the vector shown in the plane whose endpoint is (1.5,1,0)⁠. The endpoint of this vector is not (3,2,0)⁠; what is it?

(b)

Redraw the picture, including the parallelogram in the plane that shows the sum of the vectors ending at (1.5,0,1) and (1.5,1,0)⁠. The endpoint of the sum, on the diagonal, is not (3,1,1)⁠; what is it?

Answer (2 parts)

(a)

The vector shown

is not the result of doubling

(200)+(−0.510)⋅1=(1.510)

instead it is the result of doubling the parameter.

(200)+(−0.510)⋅2=(120)

This compares the lengths.

12+22+02=(2⋅0.5)2+(2⋅1)2+(2⋅0)2=22⋅(0.52+12+02)=2⋅0.52+12+02

(b)

The vector

is not the result of adding

((200)+(−0.510)⋅1)+((200)+(−0.501)⋅1)

instead it is

(200)+(−0.510)⋅1+(−0.501)⋅1=(111)

which adds the parameters.

One.II.2.27 §II.2 source ↗

Show that the line segments (a1,a2)(b1,b2) and (c1,c2)(d1,d2) have the same lengths and slopes if b1−a1=d1−c1 and b2−a2=d2−c2⁠. Is that only if?

Answer

The “if” half is straightforward. If b1−a1=d1−c1 and b2−a2=d2−c2 then

(b1−a1)2+(b2−a2)2=(d1−c1)2+(d2−c2)2

so they have the same lengths, and the slopes are just as easy:

b2−a2b1−a1=d2−c2d1−a1

(if the denominators are 0 they both have undefined slopes).

For “only if”, assume that the two segments have the same length and slope (the case of undefined slopes is easy; we will do the case where both segments have a slope m⁠). Also assume, without loss of generality, that a1<b1 and that c1<d1⁠. The first segment is (a1,a2)(b1,b2)= {(x, y)|y= mx+ n1, x∈[a1..b1]} (for some intercept n1⁠) and the second segment is (c1,c2)(d1,d2)= {(x, y)|y= mx+ n2, x∈[c1..d1]} (for some n2⁠). Then the lengths of those segments are

(b1−a1)2+((mb1+n1)−(ma1+n1))2=(1+m2)(b1−a1)2

and, similarly, (1+m2)(d1−c1)2⁠. Therefore, |b1−a1|=|d1−c1|⁠. Thus, as we assumed that a1<b1 and c1<d1⁠, we have that b1−a1=d1−c1⁠.

The other equality is similar.

One.II.2.28 §II.2 source ↗

Is |u→1+ ⋯+ u→n|≤|u→1|+ ⋯+ |u→n|⁠? If it is true then it would generalize the Triangle Inequality.

Answer

Yes; we can prove this by induction.

Assume that the vectors are in some ℝk⁠. Clearly the statement applies to one vector. The Triangle Inequality is this statement applied to two vectors. For an inductive step assume the statement is true for n or fewer vectors. Then this

|u→1+⋯+u→n+u→n+1|≤|u→1+⋯+u→n|+|u→n+1|

follows by the Triangle Inequality for two vectors. Now the inductive hypothesis, applied to the first summand on the right, gives that as less than or equal to |u→1|+⋯+|u→n|+|u→n+1|⁠.

One.II.2.29 §II.2 source ↗

What is the ratio between the sides in the Cauchy-Schwarz inequality?

Answer

By definition

u→∙v→|u→||v→|=cos⁡θ

where θ is the angle between the vectors. Thus the ratio is |cos⁡θ|⁠.

One.II.2.30 §II.2 source ↗

Describe the angle between two vectors in ℝ1⁠.

Answer

We can find the angle between (a) and (b) (for a,b≠0⁠) with

arccos⁡(aba2b2).

If a or b is zero then the angle is π/2 radians. Otherwise, if a and b are of opposite signs then the angle is π radians, else the angle is zero radians.

One.II.2.31 §II.2 source ↗

Give a simple necessary and sufficient condition to determine whether the angle between two vectors is acute, right, or obtuse.

Answer

The angle between u→ and v→ is acute if u→∙v→>0⁠, is right if u→∙v→=0⁠, and is obtuse if u→∙v→<0⁠. That’s because, in the formula for the angle, the denominator is never negative.

One.II.2.32 §II.2 source ↗

Generalize to ℝn the converse of the Pythagorean Theorem, that if u→ and v→ are perpendicular then |u→+v→|2=|u→|2+|v→|2.

Answer

Suppose that u→,v→∈ℝn⁠. If u→ and v→ are perpendicular then

|u→+v→|2=(u→+v→)∙(u→+v→)=u→∙u→+2u→∙v→+v→∙v→=u→∙u→+v→∙v→=|u→|2+|v→|2

(the third equality holds because u→∙v→=0⁠).

One.II.2.33 §II.2 source ↗

Show that |u→|=|v→| if and only if u→+v→ and u→−v→ are perpendicular. Give an example in ℝ2.

Answer

Where u→,v→∈ℝn⁠, the vectors u→+v→ and u→−v→ are perpendicular if and only if 0= (u→+v→)∙(u→−v→)= u→∙u→− v→∙v→⁠, which shows that those two are perpendicular if and only if u→∙u→=v→∙v→⁠. That holds if and only if |u→|=|v→|⁠.

One.II.2.34 §II.2 source ↗

Show that if a vector is perpendicular to each of two others then it is perpendicular to each vector in the plane they generate. (Remark. They could generate a degenerate plane — a line or a point — but the statement remains true.)

Answer

Suppose u→∈ℝn is perpendicular to both v→∈ℝn and w→∈ℝn⁠. Then, for any k,m∈ℝ we have this.

u→∙(kv→+mw→)=k(u→∙v→)+m(u→∙w→)=k(0)+m(0)=0

One.II.2.35 §II.2 source ↗

Prove that, where u→,v→∈ℝn are nonzero vectors, the vector

u→|u→|+v→|v→|

bisects the angle between them. Illustrate in ℝ2.

Answer

We will show something more general: if |z→1|=|z→2| for z→1,z→2∈ℝn⁠, then z→1+z→2 bisects the angle between z→1 and z→2

(we ignore the case where z→1 and z→2 are the zero vector).

The z→1+z→2=0→ case is easy. For the rest, by the definition of angle, we will be finished if we show this.

z→1∙(z→1+z→2)|z→1||z→1+z→2|=z→2∙(z→1+z→2)|z→2||z→1+z→2|

But distributing inside each expression gives

z→1∙z→1+z→1∙z→2|z→1||z→1+z→2|z→2∙z→1+z→2∙z→2|z→2||z→1+z→2|

and z→1∙z→1= |z→1|2= |z→2|2= z→2∙z→2⁠, so the two are equal.

One.II.2.36 §II.2 source ↗

Verify that the definition of angle is dimensionally correct: (1) if k>0 then the cosine of the angle between ku→ and v→ equals the cosine of the angle between u→ and v→⁠, and (2) if k<0 then the cosine of the angle between ku→ and v→ is the negative of the cosine of the angle between u→ and v→⁠.

Answer

We can show the two statements together. Let u→,v→∈ℝn⁠, write

u→=(u1⋮un)v→=(v1⋮vn)

and calculate.

cos⁡θ=ku1v1+⋯+kunvn(ku1)2+⋯+(kun)2b12+⋯+bn2=k|k|u→⋅v→|u→||v→|=±u→∙v→|u→||v→|

One.II.2.37 §II.2 source ↗

Show that the inner product operation is linear: for u→,v→,w→∈ℝn and k,m∈ℝ⁠, u→∙(kv→+mw→)= k(u→∙v→)+ m(u→∙w→)⁠.

Answer

Let

u→=(u1⋮un),v→=(v1⋮vn)w→=(w1⋮wn)

and then

u→∙(kv→+mw→)=(u1⋮un)∙((kv1⋮kvn)+(mw1⋮mwn))=(u1⋮un)∙(kv1+mw1⋮kvn+mwn)=u1(kv1+mw1)+⋯+un(kvn+mwn)=ku1v1+mu1w1+⋯+kunvn+munwn=(ku1v1+⋯+kunvn)+(mu1w1+⋯+munwn)=k(u→∙v→)+m(u→∙w→)

as required.

One.II.2.38 §II.2 source ↗

[Cleary] Astrologers claim that they can recognize trends in personality and fortune that depend on an individual’s birthday by incorporating where the stars were 20⁠ centuries ago. Suppose that instead of star-gazers coming up with stuff, math teachers who like linear algebra (we’ll call them vectologers) had come up with a similar system as follows: Consider your birthday as a row vector (monthday)⁠. For instance, I was born on July 12 so my vector would be (712)⁠. Vectologers assert that how well individuals get along with each other depends on the angle between vectors. The smaller the absolute value of that angle, the more harmonious the relationship.

(a)

Find the angle between your vector and mine, in radians.

(b)

Would you get along better with me, or with a professor born on September 19⁠?

(c)

For maximum harmony in a relationship, when should the other person be born?

(d)

Is there a person with whom you have a “worst case” relationship, i.e., your vector and theirs are orthogonal? If so, what are the birthdate(s) for such people? If not, explain why not.

Answer (4 parts)

(a)

My birthday of October 12 gives this.

θ=arccos⁡((712)∙(1012)|(712)|⋅|(1012)|)=arccos⁡(214244193)≈0.17 rad

(b)

Applying the same equation to (919) gives about 0.09⁠ radians.

(c)

The angle will measure 0⁠ radians if the other person is born on the same day. It will also measure 0 if one birthday is a scalar multiple of the other. For instance, a person born on Mar 6 would be harmonious with a person born on Feb 4⁠.

Given a birthday, we can get Sage to plot the angle for other dates. This example shows the relationship of all dates with July 12.

  sage: plot3d(lambda x, y: math.acos((x*7+y*12)/(math.sqrt(7**2+12**2)*math.sqrt(x**2+y**2))),
                                                                 (1,12),(1,31))

(d)

We want to maximize this.

θ=arccos⁡((712)∙(md)|(712)|⋅|(md)|)

Of course, we cannot take m or d negative and so we cannot get a vector orthogonal to the given one. This Python script finds the largest angle by brute force.

  import math
  days={1:31,  # Jan
        2:29, 3:31, 4:30, 5:31, 6:30, 7:31, 8:31, 9:30, 10:31, 11:30, 12:31}
  BDAY=(7,12)
  max_res=0
  max_res_date=(-1,-1)
  for month in range(1,13):
      for day in range(1,days[month]+1):
          num=BDAY[0]*month+BDAY[1]*day
          denom=math.sqrt(BDAY[0]**2+BDAY[1]**2)*math.sqrt(month**2+day**2)
          if denom>0:
              res=math.acos(min(num*1.0/denom,1))
              print "day:",str(month),str(day)," angle:",str(res)
              if res>max_res:
                  max_res=res
                  max_res_date=(month,day)
  print "For ",str(BDAY),"worst case",str(max_res),"rads, date",str(max_res_date)
  print "  That is ",180*max_res/math.pi,"degrees"

The result is

  For  (7, 12) worst case 0.95958064648 rads, date (12, 1)
    That is  54.9799211457 degrees

A more conceptual approach is to consider the relation of all points (month,day) to the point (7,12)⁠. The picture below makes clear that the answer is either Dec 1 or Jan 31⁠, depending on which is further from the birthdate. The dashed line bisects the angle between the line from the origin to Dec 1⁠, and the line from the origin to Jan 31⁠. Birthdays above the line are furthest from Dec 1 and birthdays below the line are furthest from Jan 31⁠.

One.II.2.39 §II.2 source ↗

[Am. Math. Mon., Feb. 1933] A ship is sailing with speed and direction v→1⁠; the wind blows apparently (judging by the vane on the mast) in the direction of a vector a→⁠; on changing the direction and speed of the ship from v→1 to v→2 the apparent wind is in the direction of a vector b→⁠.

Find the vector velocity of the wind.

Answer

This is how the answer was given in the cited source. The actual velocity v→ of the wind is the sum of the ship’s velocity and the apparent velocity of the wind. Without loss of generality we may assume a→ and b→ to be unit vectors, and may write

v→=v→1+sa→=v→2+tb→

where s and t are undetermined scalars. Take the dot product first by a→ and then by b→ to obtain

s−ta→∙b→=a→∙(v→2−v→1)sa→∙b→−t=b→∙(v→2−v→1)

Multiply the second by a→∙b→⁠, subtract the result from the first, and find

s=[a→−(a→∙b→)b→]∙(v→2−v→1)1−(a→∙b→)2.

Substituting in the original displayed equation, we get

v→=v→1+[a→−(a→∙b→)b→]∙(v→2−v→1)a→1−(a→∙b→)2.

One.II.2.40 §II.2 source ↗

Verify the Cauchy-Schwarz inequality by first proving Lagrange’s identity:

(∑1≤j≤najbj)2=(∑1≤j≤naj2)(∑1≤j≤nbj2)−∑1≤k<j≤n(akbj−ajbk)2

and then noting that the final term is positive. This result is an improvement over Cauchy-Schwarz because it gives a formula for the difference between the two sides. Interpret that difference in ℝ2⁠.

Answer

We use induction on n⁠.

In the n=1 base case the identity reduces to

(a1b1)2=(a12)(b12)−0

and clearly holds.

For the inductive step assume that the formula holds for the 0⁠, …, n cases. We will show that it then holds in the n+1 case. Start with the right-hand side

(∑1≤j≤n+1aj2)(∑1≤j≤n+1bj2)−∑1≤k<j≤n+1(akbj−ajbk)2=[(∑1≤j≤naj2)+an+12][(∑1≤j≤nbj2)+bn+12]−[∑1≤k<j≤n(akbj−ajbk)2+∑1≤k≤n(akbn+1−an+1bk)2]=(∑1≤j≤naj2)(∑1≤j≤nbj2)+∑1≤j≤nbj2an+12+∑1≤j≤naj2bn+12+an+12bn+12−[∑1≤k<j≤n(akbj−ajbk)2+∑1≤k≤n(akbn+1−an+1bk)2]=(∑1≤j≤naj2)(∑1≤j≤nbj2)−∑1≤k<j≤n(akbj−ajbk)2+∑1≤j≤nbj2an+12+∑1≤j≤naj2bn+12+an+12bn+12−∑1≤k≤n(akbn+1−an+1bk)2

and apply the inductive hypothesis.

=(∑1≤j≤najbj)2+∑1≤j≤nbj2an+12+∑1≤j≤naj2bn+12+an+12bn+12−[∑1≤k≤nak2bn+12−2∑1≤k≤nakbn+1an+1bk+∑1≤k≤nan+12bk2]=(∑1≤j≤najbj)2+2(∑1≤k≤nakbn+1an+1bk)+an+12bn+12=[(∑1≤j≤najbj)+an+1bn+1]2

to derive the left-hand side.

One.III.1.8 §III.1 source ↗

Use Gauss-Jordan reduction to solve each system.

(a)

x+y=2x−y=0

(b)

x−z=42x+2y=1

(c)

3x−2y=16x+y=1/2

(d)

2x−y=−1x+3y−z=5y+2z=5

Answer (5 parts)

These answers show only the Gauss-Jordan reduction. With it, describing the solution set is easy.

(a)

The solution set contains only a single element.

(1121−10)⟶−ρ1+ρ2(1120−2−2)⟶−(1/2)ρ2(112011)⟶−ρ2+ρ1(101011)

(b)

The solution set has one parameter.

(10−142201)⟶−2ρ1+ρ2(10−14022−7)⟶(1/2)ρ2(10−14011−7/2)

(c)

There is a unique solution.

(3−21611/2)⟶−2ρ1+ρ2(3−2105−3/2)⟶(1/5)ρ2(1/3)ρ1(1−2/31/301−3/10)⟶(2/3)ρ2+ρ1(102/1501−3/10)

(d)

A row swap in the second step makes the arithmetic easier.

(2−10−113−150125)⟶−(1/2)ρ1+ρ2(2−10−107/2−111/20125)⟶ρ2↔ρ3(2−10−1012507/2−111/2)⟶−(7/2)ρ2+ρ3(2−10−1012500−8−12)⟶−(1/8)ρ2(1/2)ρ1(1−1/20−1/201250013/2)⟶−2ρ3+ρ2(1−1/20−1/201020013/2)⟶(1/2)ρ2+ρ1(1001/201020013/2)

One.III.1.9 §III.1 source ↗

Do Gauss-Jordan reduction.

(a)

x+y−z=32x−y−z=13x+y+2z=0

(b)

x+y+2z=02x−y+z=14x+y+5z=1

Answer (2 parts)

(a)

(11−132−1−113120)⟶−3ρ1+ρ3−2ρ1+ρ2(11−130−31−50−25−9)⟶−(2/3)ρ2+ρ3(11−130−31−50013/3−17/3)⟶(3/13)ρ3−(1/3)ρ2(11−1301−1/35/3001−17/13)⟶(1/3)ρ3+ρ2ρ3+ρ1(11022/1301016/13001−17/13)⟶−ρ2+ρ1(1006/1301016/13001−17/13)

(b)

(11202−1114151)⟶−4ρ1+ρ3−2ρ1+ρ2(11200−3−310−3−31)⟶−ρ2+ρ3(11200−3−310000)⟶−(1/3)ρ2(1120011−1/30000)⟶−ρ2+ρ1(1011/3011−1/30000)

One.III.1.10 §III.1 source ↗

Find the reduced echelon form of each matrix.

(a)

(2113)

(b)

(131204−1−3−3)

(c)

(103121421534812)

(d)

(013200561515)

Answer (5 parts)

Use Gauss-Jordan reduction.

(a)

The reduced echelon form is all zeroes except for a diagonal of ones.

⟶−(1/2)ρ1+ρ2(2105/2)⟶(2/5)ρ2(1/2)ρ1(11/201)⟶−(1/2)ρ2+ρ1(1001)

(b)

As in the prior problem, the reduced echelon form is all zeroes but for a diagonal of ones.

⟶ρ1+ρ3−2ρ1+ρ2(1310−6200−2)⟶−(1/2)ρ3−(1/6)ρ2(13101−1/3001)⟶−ρ3+ρ1(1/3)ρ3+ρ2(130010001)⟶−3ρ2+ρ1(100010001)

(c)

There are more columns than rows so we must get more than just a diagonal of ones.

⟶−3ρ1+ρ3−ρ1+ρ2(1031204−10304−1−2−4)⟶−ρ2+ρ3(1031204−103000−2−7)⟶−(1/2)ρ3(1/4)ρ2(1031201−1/403/400017/2)⟶−ρ3+ρ1(1030−3/201−1/403/400017/2)

(d)

As in the prior item, this is not a square matrix.

⟶ρ1↔ρ3(151500560132)⟶ρ2↔ρ3(151501320056)⟶(1/5)ρ3(151501320016/5)⟶−ρ3+ρ1−3ρ3+ρ2(15019/5010−8/50016/5)⟶−5ρ2+ρ1(10059/5010−8/50016/5)

One.III.1.11 §III.1 source ↗

Get the reduced echelon form of each.

(a)

(0212−11−2−10)

(b)

(131262−100)

Answer (2 parts)

(a)

Swap first.

⟶ρ1↔ρ2⟶ρ1+ρ3⟶ρ2+ρ3⟶(1/2)ρ2(1/2)ρ3(1/2)ρ1⟶(−1/2)ρ3+ρ2(−1/2)ρ3+ρ1⟶(1/2)ρ2+ρ1(100010001)

(b)

Here the swap is in the middle.

⟶ρ1+ρ3−2ρ1+ρ2⟶ρ2↔ρ3⟶(1/3)ρ2⟶−3ρ2+ρ1(100011/3000)

One.III.1.12 §III.1 source ↗

Find each solution set by using Gauss-Jordan reduction and then reading off the parametrization.

(a)

2x+y−z=14x−y=3

(b)

x−z=1y+2z−w=3x+2y+3z−w=7

(c)

x−y+z=0y+w=03x−2y+3z+w=0−y−w=0

(d)

a+2b+3c+d−e=13a−b+c+d+e=3

Answer (5 parts)

For the Gauss’s halves, see the answers to Chapter One’s section I.2 question Exercise 2.19.

(a)

The “Jordan” half goes this way.

⟶−(1/3)ρ2(1/2)ρ1(11/2−1/21/201−2/3−1/3)⟶−(1/2)ρ2+ρ1(10−1/62/301−2/3−1/3)

The solution set is this

{(2/3−1/30)+(1/62/31)z|z∈ℝ}

(b)

The second half is

⟶ρ3+ρ2(10−1010120300010)

so the solution is this.

{(1300)+(1−210)z|z∈ℝ}

(c)

This Jordan half

⟶ρ2+ρ1(10110010100000000000)

gives

{(0000)+(−1010)z+(−1−101)w|z,w∈ℝ}

(of course, we could omit the zero vector from the description).

(d)

The “Jordan” half

⟶−(1/7)ρ2(1231−11018/72/7−4/70)⟶−2ρ2+ρ1(105/73/71/71018/72/7−4/70)

ends with this solution set.

{(10000)+(−5/7−8/7100)c+(−3/7−2/7010)d+(−1/74/7001)e|c,d,e∈ℝ}

One.III.1.13 §III.1 source ↗

Give two distinct echelon form versions of this matrix.

(211364121515)

Answer

Routine Gauss’s Method gives one:

⟶−(1/2)ρ1+ρ3−3ρ1+ρ2(211301−2−709/21/27/2)⟶−(9/2)ρ2+ρ3(211301−2−70019/235)

and any cosmetic change, such as multiplying the bottom row by 2⁠,

(211301−2−7001970)

gives another.

One.III.1.14 §III.1 source ↗

List the reduced echelon forms possible for each size.

(a)

2×2

(b)

2×3

(c)

3×2

(d)

3×3

Answer (5 parts)

In the cases listed below, we take a,b∈ℝ⁠. Thus, some canonical forms listed below actually include infinitely many cases. In particular, they includes the cases a=0 and b=0⁠.

(a)

(0000)⁠, (1a00)⁠, (0100)⁠, (1001)

(b)

(000000)⁠, (1ab000)⁠, (01a000)⁠, (001000)⁠, (10a01b)⁠, (1a0001)⁠, (010001)

(c)

(000000)⁠, (1a0000)⁠, (010000)⁠, (100100)

(d)

(000000000)⁠, (1ab000000)⁠, (01a000000)⁠, (010001000)⁠, (001000000)⁠, (10a01b000)⁠, (1a0001000)⁠, (100010001)

One.III.1.15 §III.1 source ↗

What results from applying Gauss-Jordan reduction to a nonsingular matrix?

Answer

A nonsingular homogeneous linear system has a unique solution. So a nonsingular matrix must reduce to a (square) matrix that is all 0⁠’s except for 1⁠’s down the upper-left to lower-right diagonal, such as these.

(1001)(100010001)

One.III.1.16 §III.1 source ↗

Decide whether each relation is an equivalence on the set of 2×2 matrices.

(a)

two matrices are related if they have the same entry in the first row and first column

(b)

two matrices are related if they have the same entry in the first row and first column, or the same entry in the second row and second column

Answer (2 parts)

(a)

This is an equivalence.

We can write M1∼M2 if they are related. The ∼⁠ relation is reflexive because any matrix has the same 1,1 entry as itself. The relation is symmetric because if M1 has the same 1,1 entry as M2 then clearly also M2 has the same 1,1 entry as M1⁠. Finally, the relation is transitive because if M1∼M2 so they have the same 1,1 entry as each other, and M2∼M3 so they have the same as each other, then all three have the same 1,1 entry, and M1∼M3⁠.

(b)

This is not an equivalence because it is not transitive. The first and second matrix below are related by their 1,1 entries, and the second and third are related by their 2,2⁠ entries. But the first and third are not related.

(1000)(100−1)(000−1)

One.III.1.17 §III.1 source ↗

[Cleary] Consider the following relationship on the set of 2×2 matrices: we say that A is sum-what like B if the sum of all of the entries in A is the same as the sum of all the entries in B⁠. For instance, the zero matrix would be sum-what like the matrix whose first row had two sevens, and whose second row had two negative sevens. Prove or disprove that this is an equivalence relation on the set of 2×2 matrices.

Answer

It is an equivalence relation. To prove that we must check that the relation is reflexive, symmetric, and transitive.

Assume that all matrices are 2×2⁠. For reflexive, we note that a matrix has the same sum of entries as itself. For symmetric, we assume A has the same sum of entries as B and obviously then B has the same sum of entries as A⁠. Transitivity is no harder — if A has the same sum of entries as B and B has the same sum of entries as C then A has the same as C⁠.

One.III.1.18 §III.1 source ↗

The proof of Lemma 1.5 contains a reference to the i≠j condition on the row combination operation.

(a)

Write down a 2×2 matrix with nonzero entries, and show that the −1⋅ρ1+ρ1 operation is not reversed by 1⋅ρ1+ρ1⁠.

(b)

Expand the proof of that lemma to make explicit exactly where it uses the i≠j condition on combining.

Answer (2 parts)

(a)

For instance,

(1234)⟶−ρ1+ρ1(0034)⟶ρ1+ρ1(0034)

leaves the matrix changed.

(b)

This operation

(⋮ai,1⋯ai,n⋮aj,1⋯aj,n⋮)⟶kρi+ρj(⋮ai,1⋯ai,n⋮kai,1+aj,1⋯kai,n+aj,n⋮)

leaves the i⁠-th row unchanged because of the i≠j restriction. Because the i⁠-th row is unchanged, this operation

⟶−kρi+ρj(⋮ai,1⋯ai,n⋮−kai,1+kai,1+aj,1⋯−kai,n+kai,n+aj,n⋮)

returns the j⁠-th row to its original state.

One.III.1.19 §III.1 source ↗

[Cleary] Consider the set of students in a class. Which of the following relationships are equivalence relations? Explain each answer in at least a sentence.

(a)

Two students x,y are related if x has taken at least as many math classes as y⁠.

(b)

Students x,y are related if they have names that start with the same letter.

Answer (3 parts)

To be an equivalence, each relation must be reflexive, symmetric, and transitive.

(a)

This relation is not symmetric because if x has taken 4⁠ classes and y has taken 3 then x is related to y but y is not related to x⁠.

(b)

This is reflexive because x⁠’s name starts with the same letter as does x⁠’s. It is symmetric because if x⁠’s name starts with the same letter as y⁠’s then y⁠’s starts with the same letter as does x⁠’s. And it is transitive because if x⁠’s name starts with the same letter as does y⁠’s and y⁠’s name starts with the same letter as does z⁠’s then x⁠’s starts with the same letter as does z⁠’s. So it is an equivalence.

One.III.1.20 §III.1 source ↗

Show that each of these is an equivalence on the set of 2×2 matrices. Describe the equivalence classes.

(a)

Two matrices are related if they have the same product down the diagonal, that is, if the product of the entries in the upper left and lower right are equal.

(b)

Two matrices are related if they both have at least one entry that is a 1⁠, or if neither does.

Answer (3 parts)

For each we must check the three conditions of reflexivity, symmetry, and transitivity.

(a)

Any matrix clearly has the same product down the diagonal as itself, so the relation is reflexive. The relation is symmetric because if A has the same product down its diagonal as does B⁠, if a1,1⋅a2,2=b1,1⋅b2,2⁠, then B has the same product as does A⁠.

Transitivity is similar: suppose that A⁠’s product is r and that it equals B⁠’s product. Suppose also that B⁠’s product equals C⁠’s. Then all three have a product of r⁠, and A⁠’s equals C⁠’s.

There is an equivalence class for each real number, namely the class contains all 2×2 matrices whose product down the diagonal is that real.

(b)

For reflexivity, if the matrix A has a 1 entry then it is related to itself while if it does not then it is also related to itself. Symmetry also has two cases: suppose that A and B are related. If A has a 1 entry then so does B⁠, and thus B is related to A⁠. If A has no 1 then neither does B⁠, and again B is related to A⁠.

For transitivity, suppose that A is related to B and B to C⁠. If A has a 1 entry then so does B⁠, and because B is related to C⁠, therefore so does C⁠, and hence A is related to C⁠. Likewise, if A has no 1 then neither does B⁠, and consequently neither does C⁠, giving the conclusion that A is related to C⁠.

There are exactly two equivalence classes, one containing any 2×2 matrix that has at least one entry that is a 1⁠, and the other containing all the matrices that have no 1⁠’s.

One.III.1.21 §III.1 source ↗

Show that each is not an equivalence on the set of 2×2 matrices.

(a)

Two matrices A,B are related if a1,1=−b1,1⁠.

(b)

Two matrices are related if the sum of their entries are within 5⁠, that is, A is related to B if |(a1,1+⋯+a2,2)− (b1,1+⋯+b2,2)|< 5⁠.

Answer (2 parts)

(a)

This relation is not reflexive. For instance, any matrix with an upper-left entry of 1 is not related to itself.

(b)

This relation is not transitive. For these three, A is related to B⁠, and B is related to C⁠, but A is not related to C⁠.

A=(0000),B=(4000),C=(8000),

One.III.2.11 §III.2 source ↗

Decide if the matrices are row equivalent.

(a)

(1248),(0112)

(b)

(1023−115−15), (1020210204)

(c)

(21−111043−1),(1020210)

(d)

(111−122),(03−1225)

(e)

(111003),(0121−11)

Answer (6 parts)

Bring each to reduced echelon form and compare.

(a)

The first gives

⟶−4ρ1+ρ2(1200)

while the second gives

⟶ρ1↔ρ2(1201)⟶−2ρ2+ρ1(1001)

The two reduced echelon form matrices are not identical, and so the original matrices are not row equivalent.

(b)

The first is this.

⟶−5ρ1+ρ3−3ρ1+ρ2(1020−1−50−1−5)⟶−ρ2+ρ3(1020−1−5000)⟶−ρ2(102015000)

The second is this.

⟶−2ρ1+ρ3(1020210000)⟶(1/2)ρ2(102015000)

These two are row equivalent.

(c)

These two are not row equivalent because they have different sizes.

(d)

The first,

⟶ρ1+ρ2(111033)⟶(1/3)ρ2(111011)⟶−ρ2+ρ1(100011)

and the second.

⟶ρ1↔ρ2(22503−1)⟶(1/3)ρ2(1/2)ρ1(115/201−1/3)⟶−ρ2+ρ1(1017/601−1/3)

These are not row equivalent.

(e)

Here the first is

⟶(1/3)ρ2(111001)⟶−ρ2+ρ1(110001)

while this is the second.

⟶ρ1↔ρ2(1−11012)⟶ρ2+ρ1(103012)

These are not row equivalent.

One.III.2.12 §III.2 source ↗

Which of these matrices are row equivalent to each other?

(a)

(1324)

(b)

(15210)

(c)

(1−130)

(d)

(26410)

(e)

(01−10)

(f)

(3322)

Answer (7 parts)

Perform Gauss-Jordan reduction on each. Two matrices are row-equivalent if and only if they have the same reduced echelon form. Here is the reduced form for each.

(a)

(1001)

(b)

(1500)

(c)

(1001)

(d)

(1001)

(e)

(1001)

(f)

(1100)

One.III.2.13 §III.2 source ↗

Produce three other matrices row equivalent to the given one.

(a)

(134−1)

(b)

(012111234)

Answer (3 parts)

For each you can just perform some row operations on the starting matrix.

(a)

Multiplying the first row by 3 gives this.

(394−1)

(There is no sense to this particular choice of row operation; it is just the first thing that came to mind.) Two other row operations are a row swap ρ1↔ρ2 and adding ρ1+ρ2⁠.

(4−113)(1352)

(b)

Doing the same three arbitrary row operations gives these three.

(036111234)(111012234)(012123234)

One.III.2.14 §III.2 source ↗

Perform Gauss’s Method on this matrix. Express each row of the final matrix as a linear combination of the rows of the starting matrix.

(1213−10040)

Answer

The Gaussian reduction is routine.

(1213−10040)⟶−3ρ1+ρ2(1210−7−3040)⟶(4/7)ρ2+ρ3(1210−7−300−12/7)

Denoting those matrices A⁠, D⁠, and B respectively, we have this.

(α1α2α3)⟶−3ρ1+ρ2(δ1=α1δ2=−3α1+α2δ3=α3)⟶(4/7)ρ2+ρ3(β1=α1β2=−3α1+α2β3=−(12/7)α1+(4/7)α2+α3)

One.III.2.15 §III.2 source ↗

Describe the matrices in each of the classes represented in Example 2.10.

Answer

First, the only matrix row equivalent to the matrix of all 0⁠’s is itself (since row operations have no effect).

Second, the matrices that reduce to

(1a00)

have the form

(bbacca)

(where a,b,c∈ℝ⁠, and b and c are not both zero).

Next, the matrices that reduce to

(0100)

have the form

(0a0b)

(where a,b∈ℝ⁠, and not both are zero).

Finally, the matrices that reduce to

(1001)

are the nonsingular matrices. That’s because a linear system for which this is the matrix of coefficients will have a unique solution, and that is the definition of nonsingular. (Another way to say the same thing is to say that they fall into none of the above classes.)

One.III.2.16 §III.2 source ↗

Describe all matrices in the row equivalence class of these.

(a)

(1000)

(b)

(1224)

(c)

(1113)

Answer (3 parts)

(a)

They have the form

(a0b0)

where at least one of a,b∈ℝ is nonzero.

(b)

They have this form

(a2ab2b)

where at least one of a,b∈ℝ is nonzero.

(c)

The given matrix is nonsingular. So the row equivalence class consists of all nonsinglar 2×2 matrices. (To give a formula, they have the form

(abcd)

for a,b,c,d∈ℝ⁠) where ad−bc≠0⁠. We will see in Chapter Four that this formula determines when a 2×2 matrix is nonsingular.)

One.III.2.17 §III.2 source ↗

How many row equivalence classes are there?

Answer

Infinitely many. For instance, in

(1k00)

each k∈ℝ gives a different class.

One.III.2.18 §III.2 source ↗

Can row equivalence classes contain different-sized matrices?

Answer

No. Row operations do not change the size of a matrix.

One.III.2.19 §III.2 source ↗

How big are the row equivalence classes?

(a)

Show that for any matrix of all zeros, the class is finite.

(b)

Do any other classes contain only finitely many members?

Answer (2 parts)

(a)

A row operation on a matrix of zeros has no effect. Thus each such matrix is alone in its row equivalence class.

(b)

No. Any nonzero entry can be rescaled.

One.III.2.20 §III.2 source ↗

Give two reduced echelon form matrices that have their leading entries in the same column but that are not row equivalent.

Answer

Here are two.

(110001)and(100001)

One.III.2.21 §III.2 source ↗

Show that any two n×n nonsingular matrices are row equivalent. Are any two singular matrices row equivalent?

Answer

Any two n×n nonsingular matrices have the same reduced echelon form, namely the matrix with all 0⁠’s except for 1⁠’s down the diagonal.

(100010⋱001)

Two same-sized singular matrices need not be row equivalent. For example, these two 2×2 singular matrices are not row equivalent.

(1100)and(1000)

One.III.2.22 §III.2 source ↗

Describe all of the row equivalence classes containing these.

(a)

2×2⁠ matrices

(b)

2×3⁠ matrices

(c)

3×2⁠ matrices

(d)

3×3⁠ matrices

Answer

Since there is one and only one reduced echelon form matrix in each class, we can just list the possible reduced echelon form matrices.

For that list, see the answer for Exercise 1.14.

One.III.2.23 §III.2 source ↗

(a)

Show that a vector β→0 is a linear combination of members of the set {β→1,…,β→n} if and only if there is a linear relationship 0→=c0β→0+⋯+cnβ→n where c0 is not zero. (Hint. Watch out for the β→0=0→ case.)

(b)

Use that to simplify the proof of Lemma 2.5.

Answer (2 parts)

(a)

If there is a linear relationship where c0 is not zero then we can subtract c0β→0 from both sides and divide by −c0 to get β→0 as a linear combination of the others. (Remark: if there are no other vectors in the set — if the relationship is, say, 0→=3⋅0→ — then the statement is still true because the zero vector is by definition the sum of the empty set of vectors.)

Conversely, if β→0 is a combination of the others β→0=c1β→1+⋯+cnβ→n then subtracting β→0 from both sides gives a relationship where at least one of the coefficients is nonzero; namely, the −1 in front of β→0⁠.

(b)

The first row is not a linear combination of the others for the reason given in the proof: in the equation of components from the column containing the leading entry of the first row, the only nonzero entry is the leading entry from the first row, so its coefficient must be zero. Thus, from the prior part of this exercise, the first row is in no linear relationship with the other rows.

Thus, when considering whether the second row can be in a linear relationship with the other rows, we can leave the first row out. But now the argument just applied to the first row will apply to the second row. (That is, we are arguing here by induction.)

One.III.2.24 §III.2 source ↗

[Trono] Three truck drivers went into a roadside cafe. One truck driver purchased four sandwiches, a cup of coffee, and ten doughnuts for $⁠8.45⁠. Another driver purchased three sandwiches, a cup of coffee, and seven doughnuts for $⁠6.30⁠. What did the third truck driver pay for a sandwich, a cup of coffee, and a doughnut?

Answer

We know that 4s+c+10d=8.45 and that 3s+c+7d=6.30⁠, and we’d like to know what s+c+d is. Fortunately, s+c+d is a linear combination of 4s+c+10d and 3s+c+7d⁠. Calling the unknown price p⁠, we have this reduction.

(41108.453176.30111p)⟶−(1/4)ρ1+ρ3−(3/4)ρ1+ρ2(41108.4501/4−1/2−0.037503/4−3/2p−2.1125)⟶−3ρ2+ρ3(41108.4501/4−1/2−0.0375000p−2.00)

The price paid is $⁠2.00⁠.

One.III.2.25 §III.2 source ↗

The Linear Combination Lemma says which equations can be gotten from Gaussian reduction of a given linear system.

1. Produce an equation not implied by this system.
3x+4y=82x+y=3

2. Can any equation be derived from an inconsistent system?

Answer

1. An easy answer is this.
0=3
For a less wiseguy-ish answer, solve the system:
(348213)⟶−(2/3)ρ1+ρ2(3480−5/3−7/3)
gives y=7/5 and x=4/5⁠. Now any equation not satisfied by (7/5,4/5) will do, e.g., 5x+5y=10⁠.

2. Every equation can be derived from an inconsistent system. For instance, here is how to derive 3x+2y=4 from 0=5⁠. First,
0=5⟶(3/5)ρ10=3⟶xρ10=3x
(validity of the x=0 case is separate but clear). Similarly, 0=2y⁠. Ditto for 0=4⁠. But now, 0+0=0 gives 3x+2y=4⁠.

One.III.2.26 §III.2 source ↗

[Hoffman & Kunze] Extend the definition of row equivalence to linear systems. Under your definition, do equivalent systems have the same solution set?

Answer

Define linear systems to be equivalent if their augmented matrices are row equivalent. The proof that equivalent systems have the same solution set is easy.

One.III.2.27 §III.2 source ↗

In this matrix

(123303145)

the first and second columns add to the third.

(a)

Show that remains true under any row operation.

(b)

Make a conjecture.

(c)

Prove that it holds.

Answer (3 parts)

(a)

The three possible row swaps are easy, as are the three possible rescalings. One of the six possible row combinations is kρ1+ρ2⁠:

(123k⋅1+3k⋅2+0k⋅3+3145)

and again the first and second columns add to the third. The other five combinations are similar.

(b)

The obvious conjecture is that row operations do not change linear relationships among columns.

(c)

A case-by-case proof follows the sketch given in the first item.

1 Topic: Computer Algebra Systems source ↗

Use the computer to solve the two problems that opened this chapter.

(a)

This is the Statics problem.

40h+15c=10025c=50+50h

(b)

This is the Chemistry problem.

7h=7j8h+1i=5j+2k1i=3j3i=6j+1k

Answer (2 parts)

(a)

Sage does this.

sage: var('h,c')
(h, c)
sage: statics = [40*h + 15*c == 100,
....:            25*c == 50 + 50*h]
sage: solve(statics, h,c)
[[h == 1, c == 4]]

Other Computer Algebra Systems have similar commands. These Maple commands

> A:=array( [[40,15],
             [-50,25]] );
> u:=array([100,50]);
> linsolve(A,u);

get the answer [1,4].

(b)

Here there is a free variable. Sage gives this.

sage: var('h,i,j,k')
(h, i, j, k)
sage: chemistry = [7*h == 7*j,
....:              8*h + 1*i == 5*j + 2*k,
....:              1*i == 3*j,
....:              3*i == 6*j + 1*k]
sage: solve(chemistry, h,i,j,k)
[[h == 1/3*r1, i == r1, j == 1/3*r1, k == r1]]

Similarly, this Maple session

> A:=array( [[7,0,-7,0],
             [8,1,-5,2],
             [0,1,-3,0],
             [0,3,-6,-1]] );
> u:=array([0,0,0,0]);
> linsolve(A,u);

prompts the same reply (but with parameter t1⁠).

2 Topic: Computer Algebra Systems source ↗

Use the computer to solve these systems from the first subsection, or conclude ‘many solutions’ or ‘no solutions’.

(a)

2x+2y=5x−4y=0

(b)

−x+y=1x+y=2

(c)

x−3y+z=1x+y+2z=14

(d)

−x−y=1−3x−3y=2

(e)

4y+z=202x−2y+z=0x+z=5x+y−z=10

(f)

2x+z+w=5y−w=−13x−z−w=04x+y+2z+w=9

Answer (6 parts)

(a)

The answer is x=2 and y=1/2⁠. A Sage session does this.

sage: var('x,y')
(x, y)
sage: system = [2*x + 2*y == 5,
....:              x - 4*y == 0]
sage: solve(system, x,y)
[[x == 2, y == (1/2)]]

A Maple session

> A:=array( [[2,2],
             [1,-4]] );
> u:=array([5,0]);
> linsolve(A,u);

gets the same answer, of course.

(b)

The answer is x=1/2 and y=3/2⁠.

sage: var('x,y')
(x, y)
sage: system = [-1*x + y == 1,
....:              x + y == 2]
sage: solve(system, x,y)
[[x == (1/2), y == (3/2)]]
(c)

This system has infinitely many solutions. In the first subsection, with z as a parameter, we got x=(43−7z)/4 and y=(13−z)/4⁠. Sage gets the same.

sage: var('x,y,z')
(x, y, z)
sage: system = [x - 3*y + z == 1,
....:           x + y + 2*z == 14]
sage: solve(system, x,y)
[[x == -7/4*z + 43/4, y == -1/4*z + 13/4]]

Maple responds with (−12+7t1,t1,13−4t1)⁠, preferring y as a parameter.

(d)

There is no solution to this system. Sage gives an empty list.

sage: var('x,y')
(x, y)
sage: system = [-1*x - y == 1,
....:           -3*x - 3*y == 2]
sage: solve(system, x,y)
[]

Similarly, When the array A and vector u are given to Maple and it is asked to linsolve(A,u), it returns no result at all; that is, it responds with no solutions.

(e)

Sage finds

sage: var('x,y,z')
(x, y, z)
sage: system = [       4*y + z == 20,
....:            2*x - 2*y + z == 0,
....:             x +        z == 5,
....:             x +    y - z == 10]
sage: solve(system, x,y,z)
[[x == 5, y == 5, z == 0]]

that the solutions is (x,y,z)=(5,5,0)⁠.

(f)

There are infinitely many solutions. Sage does this.

sage: var('x,y,z,w')
(x, y, z, w)
sage: system = [ 2*x +        z + w == 5,
....:                   y -       w == -1,
....:            3*x -        z - w == 0,
....:            4*x +  y + 2*z + w == 9]
sage: solve(system, x,y,z,w)
[[x == 1, y == r2 - 1, z == -r2 + 3, w == r2]]

Maple gives (1,−1+t1,3−t1,t1)⁠.

3 Topic: Computer Algebra Systems source ↗

Use the computer to solve these systems from the second subsection.

(a)

3x+6y=18x+2y=6

(b)

x+y=1x−y=−1

(c)

x1+x3=4x1−x2+2x3=54x1−x2+5x3=17

(d)

2a+b−c=22a+c=3a−b=0

(e)

x+2y−z=32x+y+w=4x−y+z+w=1

(f)

x+z+w=42x+y−w=23x+y+z=7

Answer (6 parts)

(a)

This system has infinitely many solutions. In the second subsection we gave the solution set as

{(60)+(−21)y|y∈ℝ}

and Sage finds

sage: var('x,y')
(x, y)
sage: system = [ 3*x + 6*y == 18,
....:              x + 2*y == 6]
sage: solve(system, x,y)
[[x == -2*r3 + 6, y == r3]]

while Maple responds with (6−2t1,t1)⁠.

(b)

The solution set has only one member.

{(01)}

Sage gives this.

sage: var('x,y')
(x, y)
sage: system = [ x + y == 1,
....:            x - y == -1]
sage: solve(system, x,y)
[[x == 0, y == 1]]
(c)

This system’s solution set is infinite

{(4−10)+(−111)x3|x3∈ℝ}

Sage gives this

sage: var('x1,x2,x3')
(x1, x2, x3)
sage: system = [   x1 +        x3 == 4,
....:              x1 - x2 + 2*x3 == 5,
....:            4*x1 - x2 + 5*x3 == 17]
sage: solve(system, x1,x2,x3)
[[x1 == -r4 + 4, x2 == r4 - 1, x3 == r4]]

and Maple gives (t1,−t1+3,−t1+4)⁠.

(d)

There is a unique solution

{(111)}

and Sage finds it.

sage: var('a,b,c')
(a, b, c)
sage: system = [ 2*a +  b - c == 2,
....:            2*a +      c == 3,
....:              a - b      == 0]
sage: solve(system, a,b,c)
[[a == 1, b == 1, c == 1]]
(e)

This system has infinitely many solutions; in the second subsection we described the solution set with two parameters.

{(5/32/300)+(−1/32/310)z+(−2/31/301)w|z,w∈ℝ}

Sage does the same.

sage: var('x,y,z,w')
(x, y, z, w)
sage: system = [   x + 2*y - z     == 3,
....:            2*x +   y     + w == 4,
....:              x -   y + z + w == 1]
sage: solve(system, x,y,z,w)
[[x == r6, y == -r5 - 2*r6 + 4, z == -2*r5 - 3*r6 + 5, w == r5]]

as does Maple (3− 2t1+ t2, t1, t2, − 2+ 3t1− 2t2)⁠.

(f)

The solution set is empty.

sage: var('x,y,z,w')
(x, y, z, w)
sage: system = [   x +      z + w == 4,
....:            2*x +  y     - w == 2,
....:            3*x +  y + z     == 7]
sage: solve(system, x,y,z,w)
[]

4 Topic: Computer Algebra Systems source ↗

What does the computer give for the solution of the general 2×2 system?

ax+cy=pbx+dy=q

Answer

Sage does this.

sage: var('x,y,a,b,c,d,p,q')
(x, y, a, b, c, d, p, q)
sage: system = [ a*x +  c*y == p,
....:            b*x +  d*y == q]
sage: solve(system, x,y)
[[x == -(d*p - c*q)/(b*c - a*d), y == (b*p - a*q)/(b*c - a*d)]]

In response to this prompting

> A:=array( [[a,c],
             [b,d]] );
> u:=array([p,q]);
> linsolve(A,u);

Maple gave this reply.

[−−dp+qc−bc+ad,−bp+aq−bc+ad]

1 Topic: Input-Output Analysis source ↗

With the steel-auto system given above, estimate next year’s total productions in these cases.

(a)

Next year’s external demands are up 200 from this year for steel and are unchanged for autos.

(b)

Next year’s external demands are up 100 for steel and are up 200 for autos.

(c)

Next year’s external demands are up 200 for steel and are up 200 for autos.

Answer (4 parts)

These answers are rounded to the number of decimal places given in the data, including when the data gives zero places.

(a)

Sage gets s≈25952 and a≈30312⁠.

sage: var('s,a')
(s, a)
sage: system = [(20053/25448)*s -  (2664/30346)*a == 17789,
....:             -(48/25448)*s + (21316/30346)*a == 21243]
sage: solve(system, s,a)
[[s == (2772443022712/106830469), a == (6476439541923/213660938)]]
sage: n(2772443022712/106830469), n(6476439541923/213660938)
(25951.8005365305, 30311.7621898814)
(b)

Sage gets s≈25857 and a≈30596⁠.

sage: system = [(20053/25448)*s -  (2664/30346)*a == 17689,
....:             -(48/25448)*s + (21316/30346)*a == 21443]
sage: solve(system, s,a)
[[s == (2762271457112/106830469), a == (6537219545323/213660938)]]
sage: n(2762271457112/106830469), n(6537219545323/213660938)
(25856.5883213711, 30596.2316112410)
(c)

Sage gets s≈25984 and a≈30597⁠.

sage: system = [(20053/25448)*s -  (2664/30346)*a == 17789,
....:             -(48/25448)*s + (21316/30346)*a == 21443]
sage: solve(system, s,a)
[[s == (2775832696312/106830469), a == (6537292375723/213660938)]]
sage: n(2775832696312/106830469), n(6537292375723/213660938)
(25983.5300012771, 30596.5724802865)

2 Topic: Input-Output Analysis source ↗

For the steel-auto system, with the external demand estimates of 17,589 and 21,243 discussed above, what will be the value of steel used by steel, the value of steel used by auto, etc.?

Answer

This Sage session

sage: var('a,s')
(a, s)
sage: eqns=[(20053/25448)*s - (2664/30346)*a == 17589,
....:       (-48/25448)*s + (21316/30346)*a == 21243] 
sage: solve(eqns, s, a)
[[s == (2745320544312/106830469), a == (6476293881123/213660938)]]
sage: s_by_s = (5395/25448)*2745320544312/106830469
sage: s_by_s
582010544505/106830469
sage: n(s_by_s)
5447.98267716114
sage: s_by_a = (2664/30346)*6476293881123/213660938  
sage: n(s_by_a)
2660.93449955743
sage: a_by_s = (48/25448)*2745320544312/106830469
sage: n(a_by_s)
48.4713936058822
sage: a_by_a = (9030/30346)*6476293881123/213660938
sage: n(a_by_a)
9019.60905818451

gives this table.

used by
steel
used by
auto
used by
others
total
value of
steel
544826611758925698
value of
auto
4890202124330311

For comparison here is the original table.

used by
steel
used by
auto
used by
others
total
value of
steel
539526641738925448
value of
auto
4890302126830346

3 Topic: Input-Output Analysis source ↗

In the steel-auto system, the ratio for the use of steel by the auto industry is 2,664/30,346⁠, about 0.0878⁠. Imagine that a new process for making autos reduces this ratio to .0500⁠.

(a)

How will the predictions for next year’s total productions change compared to the first example discussed above, that is, where next year’s external demands are 17,589 for steel and 21,243 for autos?

(b)

Predict next year’s totals if in addition, the external demand for autos rises to be 21,500 because the new process makes cars cheaper.

Answer (2 parts)

(a)

Sage gives s=24,244⁠, a=30,307⁠.

 sage: var('s,a')
(s, a)
sage: system = [(20053/25448)*s -  0.0500*a == 17589,
....:             -(48/25448)*s + (21316/30346)*a == 21243]
sage: solve(system, s,a)
[[s == (12951731858499/534221147), a == (32381469405615/1068442294)]]
sage: n(12951731858499/534221147), n(32381469405615/1068442294)
(24244.1392131918, 30307.1767071166)
(b)

Sage gives s=24,267⁠, a=30,673⁠.

sage: system = [(20053/25448)*s -  0.0500*a == 17589,
....:             -(48/25448)*s + (21316/30346)*a == 21500]
sage: solve(system, s,a)
[[s == (12964136043940/534221147), a == (16386224431390/534221147)]]
sage: n(12964136043940/534221147), n(16386224431390/534221147)
(24267.3584090448, 30673.1107957993)

4 Topic: Input-Output Analysis source ↗

This table gives the numbers for the auto-steel system from a different year, 1947 (see [Leontief 1951]). The units here are billions of 1947 dollars.

used by
steel
used by
auto
used by
others
total
value of
steel
6.901.2818.69
value of
autos
04.4014.27
(a)

Solve for total output if next year’s external demands are: steel’s demand up 10% and auto’s demand up 15%.

(b)

How do the ratios compare to those given above in the discussion for the 1958 economy?

(c)

Illustrate the effect of how the ratios change over time by using the 1947 ratios from the first item with the 1958 external demands. Compare the predictions for total output. (Inflation was 30% between the two years, but note that inflation does not affect the ‘⁠used by/total⁠’ ratios because the factor of 1.30 is in both the numerator and the denominator.)

Answer (3 parts)

(a)

These are the equations, where the units are 1958 billions.

(11.79/18.69)s−(1.28/14.27)a=17.589−(0/18.69)s+(9.87/14.27)a=21.243

Sage gives this result (including an inital check).

sage: var('a,s')
(a, s)
sage: eqns=[(11.79/18.69)*s - (1.28/14.27)*a == 10.51,
....:       (-0)*s         + (9.87/14.27)*a == 9.87]
sage: solve(eqns,s,a)
[[s == (1869/100), a == (1427/100)]]
sage: n(1869/100)
18.6900000000000
sage: n(1427/100)
14.2700000000000
sage: eqns=[(11.79/18.69)*s - (1.28/14.27)*a == 11.561,
....:       (-0)*s         + (9.87/14.27)*a == 11.3505]
sage: solve(eqns,s,a)
[[s == (8119559/393000), a == (32821/2000)]]
sage: n(8119559/393000)
20.6604554707379
sage: n(32821/2000)
16.4105000000000
(b)

These are the 1947 ratios, to three decimal places.

1947by steelby autos
use of steel(6.90/18.69)=0.369(1.28/14.27)=0.090
use of autos(0/18.69)=0.00(4.40/14.27)=0.308

These are the 1958 ratios.

1958by steelby autos
use of steel(5395/25448)=0.212(2664/30346)=0.088
use of autos(48/25448)=0.002(9030/30346)=0.298
(c)

The external demands in 1958 are 17.589 and 21.243⁠, in billions of 1958 dollars. In 1948 dollars those are (17.589/1.3)=13.53 and (21.243/1.3)=16.34⁠. Sage gives these numbers.

sage: eqns=[(11.79/18.69)*s - (1.28/14.27)*a == 13.53,
....:       (-0)*s         + (9.87/14.27)*a == 16.34]
sage: solve(eqns,s,a)
[[s == (137466107/5541300), a == (1165859/49350)]]
sage: n(137466107/5541300)
24.8075554472777
sage: n(1165859/49350)
23.6242958459980

Converting back to 1958 dollars by multiplying by 1.3 gives 32.25 and 30.71⁠.

5 Topic: Input-Output Analysis source ↗

Predict next year’s total productions of each of the three sectors of the hypothetical economy shown below

used by
farm
used by
rail
used by
shipping
used by
others
total
value of
farm
2550100800
value of
rail
255050300
value of
shipping
15100500

if next year’s external demands are as stated.

(a)

625 for farm, 200 for rail, 475 for shipping

(b)

650 for farm, 150 for rail, 450 for shipping

Answer (2 parts)

(a)

These commands in Sage

var('f, r, s')
farm_by_farm = 25
farm_by_rail = 50
farm_by_ship = 100
farm_total = 800
# farm_by_others = farm_total -farm_by_farm - farm_by_rail - farm_by_ship
rail_by_farm = 25
rail_by_rail = 50
rail_by_ship = 50
rail_total = 300
# rail_by_others = rail_total -rail_by_farm - rail_by_rail - rail_by_ship
ship_by_farm = 15
ship_by_rail = 10
ship_by_ship = 0
ship_total = 500
# ship_by_others = ship_total - ship_by_farm - ship_by_rail - ship_by_ship

farm_by_others = 625
rail_by_others = 200
ship_by_others = 475
sys = [(farm_by_farm/farm_total)*f + (farm_by_rail/rail_total)*r + (farm_by_ship/ship_total)*s + farm_by_others == f,
(rail_by_farm/farm_total)*f + (rail_by_rail/rail_total)*r + (rail_by_ship/ship_total)*s + rail_by_others == r,
(ship_by_farm/farm_total)*f + (ship_by_rail/rail_total)*r + (ship_by_ship/ship_total)*s + ship_by_others == s ]
sols = solve(sys, f, r, s, solution_dict=True)
farm = sols[0][f]
rail = sols[0][r]
ship = sols[0][s]

give the conclusion that f=805⁠, r=330⁠, and s=501⁠.

(b)

Changing the lines in the prior item to these.

farm_by_others = 650
rail_by_others = 150
ship_by_others = 450

gives the solution f=815⁠, r=267⁠, and s=474⁠.

6 Topic: Input-Output Analysis source ↗

This table gives the interrelationships among three segments of an economy (see [Clark & Coupe]).

used by
food
used by
wholesale
used by
retail
used by
others
total
value of
food
02318467911869
value of
wholesale
393108922459122242
value of
retail
35375116041

We will do an Input-Output analysis on this system.

(a)

Fill in the numbers for this year’s external demands.

(b)

Solve the system where we get next year’s external demands by taking this year’s external demands and inflating them 10%. Do all three sector totals increase by 10%? Do they even increase at the same rate?

(c)

Solve the system where we get next year’s external demands by taking this year’s external demands and reducing them 7%. (The study source for these numbers concluded that because of the closing of a local military facility, overall personal income in the area would fall 7%, so this might be a first guess at what would happen.)

Answer (4 parts)

This is the basic setup in Sage.

var('f, w, r')
food_by_food = 0
food_by_whole = 2318
food_by_retail = 4679
food_total = 11869
food_by_others = food_total - food_by_food - food_by_whole - food_by_retail
whole_by_food = 393
whole_by_whole = 1089
whole_by_retail = 22459
whole_total = 122242
whole_by_others = whole_total - whole_by_food - whole_by_whole - whole_by_retail
retail_by_food = 3
retail_by_whole = 53
retail_by_retail = 75
retail_total = 116041
retail_by_others = retail_total - retail_by_food - retail_by_whole - retail_by_retail

sys = [(food_by_food/food_total)*f + (food_by_whole/whole_total)*w + (food_by_retail/retail_total)*r + food_by_others == f,
(whole_by_food/food_total)*f + (whole_by_whole/whole_total)*w + (whole_by_retail/retail_total)*r + whole_by_others == w,
(retail_by_food/food_total)*f + (retail_by_whole/whole_total)*w + (retail_by_retail/retail_total)*r + retail_by_others == r ]
sols = solve(sys, f, w, r, solution_dict=True)
food = sols[0][f]
whole = sols[0][w]
retail = sols[0][r]
(a)

With that, Sage gives that external food is 4872⁠, external wholesale is 98301⁠, and external retail is 115910⁠.

(b)

With ten percent inflation in the external demand numbers, the solutions are f=13056⁠, w=134466⁠, and r=127645⁠.

(c)

The seven percent deflation in external demand gives f=11038⁠, w=113685⁠, and r=107918⁠.

1 Topic: Accuracy of Computations source ↗

Using two decimal place numbers, add 253 and 2/3⁠.

Answer

Scientific notation is convenient for expressing the two-place restriction: .25×102+.67×100=.25×102. Adding the 2/3 has no effect.

2 Topic: Accuracy of Computations source ↗

Solve this system (from [Rice]) by rounding at each step to four digits.

0.0003x+1.556y=1.5590.3454x−2.346y=1.108

Answer

Reduction uses the step −(0.3454/0.0003)ρ1+ρ2⁠. First, − 0.3454/0.0003= − 1151.33333333333 and rounding to four digits gives −1151=−1.115×103. We take that the leading entry in the second row is now zero (that’s the point of the row reduction step). For the second entry, 1.556⋅(−1151)=−1790.956⁠, which rounds to four digits as −1791=−1.791×103. Adding −1791−2.346 rounds to −1793=−1.793×103. As to the third entry, 1.559⋅(−1151) rounds to −1794⁠, and adding rounds to −1793=−1.793×103.

(0.00031.5561.5590−1793−1793)

The conclusion is that y is −1793/−1793⁠, which rounds to 1.000=1.000×100. Substituting gives 1.556⋅1.000⁠, which rounds to 1.556⁠, and then 1.559−1.556 rounds to 0.003⁠. Solving 0.0003x=0.003 gives 0.003/0.0003⁠, which is y=10=1.000×101. Checking that by substitution into the original equations shows that it is correct.

3 Topic: Accuracy of Computations source ↗

Assume that the machine uses four digit numbers.

(a)

Show that it will compute (x+y)+z as unequal to x+(y+z) when x=−1.000×105⁠, and y=1.000×105⁠, and z=1.000×100⁠.

(b)

Compare the computer’s versions of (1/3)x+y=0 and (2/3)x+2y=0⁠. Is twice the first equation the same as the second?

Answer (2 parts)

(a)

For the first, x+y is zero and so (x+y)+z is 1⁠. For the other, (y+z) is 10000+1 rounded to four digits, 1.000×105. Then x+(y+z) is zero.

(b)

The first equation is 0.3333⋅x+ 1.0000⋅y= 0.0000 while the second is 0.6667⋅x+ 2.0000⋅y= 0.0000⁠. Twice the first is 0.6666⋅x+ 2.0000⋅y= 0.0000⁠, so they aren’t equal.

4 Topic: Accuracy of Computations source ↗

Ill-conditioning is not only dependent on the matrix of coefficients. This example [Hamming] shows that it can arise from an interaction between the left and right sides of the system. Let ε be a small real.

3x+2y+z=62x+2εy+2εz=2+4εx+2εy−εz=1+ε

(a)

Solve the system by hand. Notice that the ε⁠’s divide out only because there is an exact cancellation of the integer parts on the right side as well as on the left.

(b)

Solve the system by hand, rounding to two decimal places, and with ε=0.001⁠.

Answer (2 parts)

(a)

Here is the calculation.

⟶−(1/3)ρ1+ρ3−(2/3)ρ1+ρ2(32160−(4/3)+2ε−(2/3)+2ε−2+4ε0−(2/3)+2ε−(1/3)−ε−1+ε)⟶−(1/2)ρ2+ρ3(32160−(4/3)+2ε−(2/3)+2ε−2+4ε0ε−2ε−ε)

The third equation is εy−2εz=−ε⁠. Cancel the ε⁠’s and solve to get y=−1+2z⁠. Substitute into the second equation.

((−4/3)+2ε)(−1+2z)+((−2/3)+2ε)z=−2+4ε(4/3)−2ε−(10/3)z+6εz=−2+4ε((−10/3)+6ε)z=(−10/3)+6ε

Thus z=1 and so y=1⁠. Then the first equation says that x=1⁠.

(b)

As above, scientific notation is convenient to express the restriction on the number of digits. With two digits

(3.0×1002.0×1001.0×1006.0×1001.0×1002.0×10−42.0×10−4.20×1013.0×1002.0×10−4−1.0×10−41.0×100)⟶−(1/3)ρ1+ρ3−(2/3)ρ1+ρ2(3.0×1002.0×1001.0×1006.0×1000−1.3×100−6.7×10−1−2.0×1000−6.7×10−1−3.3×10−1−1.0×100)⟶−(.67/1.3)ρ2+ρ3(3.0×1002.0×1001.0×1006.0×1000−1.3×100−6.7×10−1−2.0×100001.5×10−33.1×10−3)

we get z=2.1⁠, y=2.6⁠, and x=−0.43⁠. These differ significantly from the true values found in the first part.

1 Topic: Analyzing Networks source ↗

Calculate the amperages in each part of each network.

(a)

This is a simple network.

(b)

Compare this one with the parallel case discussed above.

(c)

This is a reasonably complicated network.

Answer (3 parts)

(a)

The total resistance is 7⁠ ohms. With a 9⁠ volt potential, the flow will be 9/7⁠ amperes. The voltage drops will then be: 27/7⁠ volts across the 3⁠ ohm resistor, and 18/7⁠ volts across each of the two 2⁠ ohm resistors.

(b)

The 2⁠ ohm resistor on the left has a voltage drop of 9⁠ volts (and hence the flow through it is 9/2⁠ amperes), and the remaining portion on the right also has a voltage drop of 9⁠ volts, and so we can analyze it as in the prior item. We can also use linear systems.

The Current Law applied to junction in the the upper middle gives i0=i1+i2⁠. Applied to the lower middle junction it gives i1=i2+i3⁠. The Voltage Law applied to the leftmost circuit is that the 9⁠ volts up equals the voltage drop across the 2⁠ ohm resistor, 9=2i1⁠. The rightmost circuit has a total resistance of 3+2+2=7⁠ ohm, so the Voltage Law gives 9=7i2⁠. Here is the linear system.

i0−i1−i2=0i1+i2−i3=02i1=97i2=9

Gauss’s Method

⟶−2ρ2+ρ3(1−1−100011−1000−22900709)⟶(7/2)ρ3+ρ4⟶(1/7)ρ4(−1/2)ρ3⟶ρ4+ρ3ρ4+ρ2⟶−ρ3+ρ2ρ3+ρ1⟶ρ2+ρ1(100081/1401009/200109/7000181/14)

gives the unique solution i0=81/14⁠, i1=9/2⁠, i2=9/7⁠, and i3=81/14⁠.

(c)

Using these variables,

below is one linear system sufficient to yield the solution. The first four equations are from the Current Law. The last three are from applying the Voltage Law to the circuit involving i0⁠-⁠i1⁠-⁠i6⁠, the circuit involving i0⁠-⁠i2⁠-⁠i4⁠-⁠i5⁠-⁠i6⁠, and the circuit with i0⁠-⁠i2⁠-⁠i3⁠-⁠i5⁠-⁠i6⁠.

i0−i1−i2=0i2−i3−i4=0i3+i4−i5=0i1+i5−i6=03i1=93i2+2i4+2i5=93i2+9i3+2i5=9

Gauss’s Method gives i0=318/73≈4.36⁠, i1=3⁠, i2=99/73≈1.36⁠, i3=18/73≈0.25⁠, i4=81/73≈1.11⁠, i5=99/73≈1.36⁠, and i6=318/73≈4.36⁠.

2 Topic: Analyzing Networks source ↗

In the first network that we analyzed, with the three resistors in series, we just added to get that they acted together like a single resistor of 10⁠ ohms. We can do a similar thing for parallel circuits. In the second circuit analyzed,

we found that the current through the battery is 25/6⁠ amperes. Thus, the parallel portion is equivalent to a single resistor of 20/(25/6)=4.8⁠ ohms.

(a)

What is the equivalent resistance if we change the 12⁠ ohm resistor to 5⁠ ohms?

(b)

What is the equivalent resistance if the two are each 8⁠ ohms?

(c)

Find the formula for the equivalent resistance if the two resistors in parallel are r1⁠ ohms and r2⁠ ohms.

Answer (3 parts)

(a)

Using the variables from the earlier analysis,

i0−i1−i2=0−i0+i1+i2=05i1=208i2=20−5i1+8i2=0

Gauss’s Method gives that the current flowing in each branch is i0=13/2=6.5⁠, i1=4⁠, and i2=5/2=2.5⁠, all in amperes. Thus the parallel portion is acting as a single resistor of 20/(13/2)≈3.08⁠ ohms.

(b)

A similar analysis

i0−i1−i2=0−i0+i1+i2=08i1=208i2=20−8i1+8i2=0

gives that io=5⁠, i1=5/2=2.5⁠, and i2=5/2=2.5⁠ amperes. The equivalent resistance is 20/5=4⁠ ohms.

(c)

Another analysis like the prior ones

i0−i1−i2=0−i0+i1+i2=0r1i1=20r2i2=20−r1i1+r2i2=0

gives i2=20/r2⁠, i1=20/r1⁠, and i0=20(r1+r2)/(r1r2)⁠, all in amperes. The parallel portion is acting as a single resistor of size 20/i0=r1r2/(r1+r2)⁠ ohms. (This equation is often stated as: the equivalent resistance r satisfies 1/r=(1/r1)+(1/r2)⁠.)

3 Topic: Analyzing Networks source ↗

A Wheatstone bridge is used to measure resistance.

Show that in this circuit if the current flowing through rg is zero then r4=r2r3/r1⁠. (To operate the device, put the unknown resistance at r4⁠. At rg is a meter that shows the current. We vary the three resistances r1⁠, r2⁠, and r3 — typically they each have a calibrated knob — until the current in the middle reads 0⁠. Then the equation gives the value of r4⁠.)

Answer

Kirchoff’s Current Law, applied to the node where r1⁠, r2⁠, and rg come together, and also applied to the node where r3⁠, r4⁠, and rg come together gives these.

i1−i2−ig=0i3−i4+ig=0

Kirchoff’s Voltage law, applied to the loop in the top right, and to the loop in the bottom right, gives these.

i3r3−igrg−i1r1=0i4r4−i2r2+igrg=0

Assuming that ig is zero gives that i1=i2⁠, that i3=i4⁠, that i1r1=i3r3⁠, and that i2r2=i4r4⁠. Then rearranging the last equality

r4=i2r2i4⋅i3r3i1r1

and cancelling the i⁠’s gives the desired conclusion.

4 Topic: Analyzing Networks source ↗

Consider this traffic circle.

This is the traffic volume, in units of cars per ten minutes.

NorthPierMain
into
out of
100
75
150
150
25
50

We can set up equations to model how the traffic flows.

(a)

Adapt Kirchhoff’s Current Law to this circumstance. Is it a reasonable modeling assumption?

(b)

Label the three between-road arcs in the circle with a variable: let i1 be the number of cars going from North Avenue to Main, let i2 be the number of cars between Main and Pier, and let i3 be the number between Pier and North. Using the adapted law, for each of the three in-out intersections state an equation describing the traffic flow at that node.

(c)

Solve that system.

(d)

Interpret your solution.

(e)

Restate the Voltage Law for this circumstance. How reasonable is it?

Answer (5 parts)

(a)

An adaptation is: in any intersection the flow in equals the flow out. It does seem reasonable in this case, unless cars are stuck at an intersection for a long time.

(b)

Because 50 cars leave via Main while 25⁠ cars enter, i1−25=i2⁠. Similarly Pier’s in/out balance means that i2=i3 and North gives i3+25=i1⁠. We have this system.

i1−i2=25i2−i3=0−i1+i3=−25

(c)

The row operations ρ1+ρ2 and ρ2+ρ3 lead to the conclusion that there are infinitely many solutions. With i3 as the parameter, this is the solution set.

{(25+i3i3i3)|i3∈ℝ}

Since the problem is stated in number of cars we might restrict i3 to be a natural number.

(d)

If we picture an initially-empty circle with the given input/output behavior, we can superimpose i3⁠-many cars circling endlessly to get a new solution.

(e)

A suitable restatement might be: the number of cars entering the circle must equal the number of cars leaving. The reasonableness of this one is not as clear. Over the five minute time period we could find that a half dozen more cars entered than left, although the problem statement’s into/out table does satisfy this property. In any event, it is of no help in getting a unique solution since for that we would need to know the number of cars circling endlessly.

5 Topic: Analyzing Networks source ↗

This is a network of streets.

We can observe the hourly flow of cars into this network’s entrances, and out of its exits.

east Winooskiwest WinooskiWillowJayShelburne
into
out of
80
30
50
5
65
70
–
55
40
75

(Note that to reach Jay a car must enter the network via some other road first, which is why there is no ‘into Jay’ entry in the table. Note also that over a long period of time, the total in must approximately equal the total out, which is why both rows add to 235⁠ cars.) Once inside the network, the traffic may flow in different ways, perhaps filling Willow and leaving Jay mostly empty, or perhaps flowing in some other way. Kirchhoff’s Laws give the limits on that freedom.

(a)

Determine the restrictions on the flow inside this network of streets by setting up a variable for each block, establishing the equations, and solving them. Notice that some streets are one-way only. (Hint: this will not yield a unique solution, since traffic can flow through this network in various ways; you should get at least one free variable.)

(b)

Suppose that someone proposes construction for Winooski Avenue East between Willow and Jay, and traffic on that block will be reduced. What is the least amount of traffic flow that can we can allow on that block without disrupting the hourly flow into and out of the network?

Answer (2 parts)

(a)

Here is a variable for each unknown block; each known block has the flow shown.

We apply Kirchhoff’s principle that the flow into the intersection of Willow and Shelburne must equal the flow out to get i1+25=i2+125⁠. Doing the intersections from right to left and top to bottom gives these equations.

i1−i2=10−i1+i3=15i2+i4−i5=5−i3−i4+i6=−50i5−i7=−10−i6+i7=30

The row operation ρ1+ρ2 followed by ρ2+ρ3 then ρ3+ρ4 and ρ4+ρ5 and finally ρ5+ρ6 result in this system.

i1−i2=10−i2+i3=25i3+i4−i5=30−i5+i6=−20i6−i7=−300=0

Since the free variables are i4 and i7 we take them as parameters.

i6=i7−30i5=i6+20=(i7−30)+20=i7−10i3=−i4+i5+30=−i4+(i7−10)+30=−i4+i7+20i2=i3−25=(−i4+i7+20)−25=−i4+i7−5i1=i2+10=(−i4+i7−5)+10=−i4+i7+5

Obviously i4 and i7 have to be positive, and in fact the first equation shows that i7 must be at least 30⁠. If we start with i7⁠, then the i2⁠ equation shows that 0≤i4≤i7−5⁠.

(b)

We cannot take i7 to be zero or else i6 will be negative (this would mean cars going the wrong way on the one-way street Jay). We can, however, take i7 to be as small as 30⁠, and then there are many suitable i4⁠’s. For instance, the solution

(i1,i2,i3,i4,i5,i6,i7)=(35,25,50,0,20,0,30)

results from choosing i4=0⁠.

Loading chapter 2…

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.