$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$ $$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$

# 9-6. Solving Systems with Gaussian Elimination

Solving Systems with Gaussian Elimination
In this section, you will:
• Write the augmented matrix of a system of equations.
• Write the system of equations from an augmented matrix.
• Perform row operations on a matrix.
• Solve a system of linear equations using matrices.
<figure class="small" id="Figure_09_06_001" style="color: rgb(0, 0, 0); font-family: 'Times New Roman'; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: 1; word-spacing: 0px; -webkit-text-stroke-width: 0px;"> <figcaption>German mathematician Carl Friedrich Gauss (1777–1855).</figcaption> </figure>

Carl Friedrich Gauss lived during the late 18th century and early 19th century, but he is still considered one of the most prolific mathematicians in history. His contributions to the science of mathematics and physics span fields such as algebra, number theory, analysis, differential geometry, astronomy, and optics, among others. His discoveries regarding matrix theory changed the way mathematicians have worked for the last two centuries.

We first encountered Gaussian elimination in Systems of Linear Equations: Two Variables. In this section, we will revisit this technique for solving systems, this time using matrices.

# Writing the Augmented Matrix of a System of Equations

A matrix can serve as a device for representing and solving a system of equations. To express a system in matrix form, we extract the coefficients of the variables and the constants, and these become the entries of the matrix. We use a vertical line to separate the coefficient entries from the constants, essentially replacing the equal signs. When a system is written in this form, we call it an augmented matrix.

For example, consider the following<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>2</mn><mtext> </mtext><mo>×</mo><mtext> </mtext><mn>2</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]system of equations.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>3</mn><mi>x</mi><mo>+</mo><mn>4</mn><mi>y</mi><mo>=</mo><mn>7</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 4x−2y=5

We can write this system as an augmented matrix:

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 4 4 −2   |    7 5 ]

We can also write a matrix containing just the coefficients. This is called the coefficient matrix.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 4 4 −2 ]

A three-by-three system of equations such as

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>3</mn><mi>x</mi><mo>−</mo><mi>y</mi><mo>−</mo><mi>z</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]         x+y=5      2x−3z=2

has a coefficient matrix

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 −1 −1 1 1 0 2 0 −3 ]

and is represented by the augmented matrix

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 −1 −1 1 1 0 2 0 −3   |    0 5 2 ]

Notice that the matrix is written so that the variables line up in their own columns: x-terms go in the first column, y-terms in the second column, and z-terms in the third column. It is very important that each equation is written in standard form<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>a</mi><mi>x</mi><mo>+</mo><mi>b</mi><mi>y</mi><mo>+</mo><mi>c</mi><mi>z</mi><mo>=</mo><mi>d</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]so that the variables line up. When there is a missing variable term in an equation, the coefficient is 0.

Given a system of equations, write an augmented matrix.

1. Write the coefficients of the x-terms as the numbers down the first column.
2. Write the coefficients of the y-terms as the numbers down the second column.
3. If there are z-terms, write the coefficients as the numbers down the third column.
4. Draw a vertical line and write the constants to the right of the line.
Writing the Augmented Matrix for a System of Equations

Write the augmented matrix for the given system of equations.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>   </mtext><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>−</mo><mi>z</mi><mo>=</mo><mn>3</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]  2x−y+2z=6   x−3y+3z=4

The augmented matrix displays the coefficients of the variables, and an additional column for the constants.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 2 −1 2 −1 2 1 −3 3   |    3 6 4 ]

Write the augmented matrix of the given system of equations.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>4</mn><mi>x</mi><mn>−3</mn><mi>y</mi><mo>=</mo><mn>11</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 3x+2y=4

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 4 −3 3   2 | 11   4 ]

# Writing a System of Equations from an Augmented Matrix

We can use augmented matrices to help us solve systems of equations because they simplify operations when the systems are not encumbered by the variables. However, it is important to understand how to move back and forth between formats in order to make finding solutions smoother and more intuitive. Here, we will use the information in an augmented matrix to write the system of equations in standard form.

Writing a System of Equations from an Augmented Matrix Form

Find the system of equations from the augmented matrix.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 −3 −5 2 −5 −4 −3 5 4   |    −2 5 6 ]

When the columns represent the variables<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>x</mi><mo>,</mo><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mi>y</mi><mo>,</mo><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]and<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>z</mi><mo>,</mo></mrow></annotation-xml></semantics>[/itex]

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 −3 −5 2 −5 −4 −3 5 4   |    −2 5 6 ]→        x−3y−5z=−2     2x−5y−4z=5 −3x+5y+4z=6

Write the system of equations from the augmented matrix.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 −1   1 2 −1   3 0    1   1    |      5   1 −9 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable><mtr><mtd><mi>x</mi><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>−</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>y</mi><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>+</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>z</mi><mo>=</mo><mn>5</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 2x  −  y   +  3z=1             y   +    z=−9

# Performing Row Operations on a Matrix

Now that we can write systems of equations in augmented matrix form, we will examine the various row operations that can be performed on a matrix, such as addition, multiplication by a constant, and interchanging rows.

Performing row operations on a matrix is the method we use for solving a system of equations. In order to solve the system of equations, we want to convert the matrix to row-echelon form, in which there are ones down the main diagonal from the upper left corner to the lower right corner, and zeros in every position below the main diagonal as shown.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable><mtr><mtd><mrow><mtext>Row-echelon form</mtext></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] [ 1 a b 0 1 d 0 0 1 ]

We use row operations corresponding to equation operations to obtain a new matrix that is row-equivalent in a simpler form. Here are the guidelines to obtaining row-echelon form.

1. In any nonzero row, the first nonzero number is a 1. It is called a leading 1.
2. Any all-zero rows are placed at the bottom on the matrix.
3. Any leading 1 is below and to the right of a previous leading 1.
4. Any column containing a leading 1 has zeros in all other positions in the column.

To solve a system of equations we can perform the following row operations to convert the coefficient matrix to row-echelon form and do back-substitution to find the solution.

1. Interchange rows. (Notation:<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><msub/></mrow></annotation-xml></semantics>[/itex] R i  ↔   R j)
2. Multiply a row by a constant. (Notation:<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>c</mi><msub/></mrow></annotation-xml></semantics>[/itex] R i)
3. Add the product of a row multiplied by a constant to another row. (Notation:<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><msub/></mrow></annotation-xml></semantics>[/itex] R i +c R j )

Each of the row operations corresponds to the operations we have already learned to solve systems of equations in three variables. With these operations, there are some key moves that will quickly achieve the goal of writing a matrix in row-echelon form. To obtain a matrix in row-echelon form for finding solutions, we use Gaussian elimination, a method that uses row operations to obtain a 1 as the first entry so that row 1 can be used to convert the remaining rows.

Gaussian Elimination

The Gaussian elimination method refers to a strategy used to obtain the row-echelon form of a matrix. The goal is to write matrix<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>A</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]with the number 1 as the entry down the main diagonal and have all zeros below.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mi>A</mi><mo>=</mo><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ] → After Gaussian elimination A=[ 1    b 12    b 13 0   1    b 23 0   0   1 ]

The first step of the Gaussian strategy includes obtaining a 1 as the first entry, so that row 1 may be used to alter the rows below.

Given an augmented matrix, perform row operations to achieve row-echelon form.

1. The first equation should have a leading coefficient of 1. Interchange rows or multiply by a constant, if necessary.
2. Use row operations to obtain zeros down the first column below the first entry of 1.
3. Use row operations to obtain a 1 in row 2, column 2.
4. Use row operations to obtain zeros down column 2, below the entry of 1.
5. Use row operations to obtain a 1 in row 3, column 3.
6. Continue this process for all rows until there is a 1 in every entry down the main diagonal and there are only zeros below.
7. If any rows contain all zeros, place them at the bottom.
Solving a<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>2</mn><mo>×</mo><mn>2</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]System by Gaussian Elimination

Solve the given system by Gaussian elimination.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>2</mn><mi>x</mi><mo>+</mo><mn>3</mn><mi>y</mi><mo>=</mo><mn>6</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]     x−y= 1 2

First, we write this as an augmented matrix.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 2 3 1 −1   |    6 1 2 ]

We want a 1 in row 1, column 1. This can be accomplished by interchanging row 1 and row 2.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><msub><mi>R</mi></msub></mrow></annotation-xml></semantics>[/itex] 1 ↔ R 2 →[ 1 −1 2 3 | 1 2 6 ]

We now have a 1 as the first entry in row 1, column 1. Now let’s obtain a 0 in row 2, column 1. This can be accomplished by multiplying row 1 by<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>−2</mn><mo>,</mo></mrow></annotation-xml></semantics>[/itex]and then adding the result to row 2.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mn>−2</mn><msub/></mrow></annotation-xml></semantics>[/itex] R 1 + R 2 = R 2 →[ 1 −1 0 5 | 1 2 5 ]

We only have one more step, to multiply row 2 by<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mfrac/></mrow></annotation-xml></semantics>[/itex] 1 5 .

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mfrac><mn>1</mn></mfrac></mrow></annotation-xml></semantics>[/itex] 5 R 2 = R 2 →[ 1 −1 0 1 | 1 2 1 ]

Use back-substitution. The second row of the matrix represents<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mo>=</mo><mn>1.</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]Back-substitute<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mo>=</mo><mn>1</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]into the first equation.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>x</mi><mo>−</mo><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo><mo>=</mo><mfrac/></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 1 2          x= 3 2

The solution is the point<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 2 ,1 ).

Solve the given system by Gaussian elimination.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>4</mn><mi>x</mi><mo>+</mo><mn>3</mn><mi>y</mi><mo>=</mo><mn>11</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]   x−3y=−1

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 2, 1 )

Using Gaussian Elimination to Solve a System of Equations

Use Gaussian elimination to solve the given<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>2</mn><mtext> </mtext><mo>×</mo><mtext> </mtext><mn>2</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex] system of equations.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>  </mtext><mn>2</mn><mi>x</mi><mo>+</mo><mi>y</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 4x+2y=6

Write the system as an augmented matrix.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 2 1 4 2   |    1 6 ]

Obtain a 1 in row 1, column 1. This can be accomplished by multiplying the first row by<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mfrac/></mrow></annotation-xml></semantics>[/itex] 1 2 .

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mfrac><mn>1</mn></mfrac></mrow></annotation-xml></semantics>[/itex] 2 R 1 = R 1 →[ 1 1 2 4 2   |    1 2 6 ]

Next, we want a 0 in row 2, column 1. Multiply row 1 by<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>−4</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]and add row 1 to row 2.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mn>−4</mn><msub/></mrow></annotation-xml></semantics>[/itex] R 1 + R 2 = R 2 →[ 1 1 2 0 0   |    1 2 4 ]

The second row represents the equation<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>0</mn><mo>=</mo><mn>4.</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]Therefore, the system is inconsistent and has no solution.

Solving a Dependent System

Solve the system of equations.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>3</mn><mi>x</mi><mo>+</mo><mn>4</mn><mi>y</mi><mo>=</mo><mn>12</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 6x+8y=24

Perform row operations on the augmented matrix to try and achieve row-echelon form.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mi>A</mi><mo>=</mo><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 4 6 8 | 12 24 ]
<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>−</mo><mfrac/></mrow></mtd></mtr></mtable></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 1 2 R 2 + R 1 = R 1 →[ 0 0 6 8 |    0 24 ] R 1 ↔ R 2 →[ 6 8 0 0 | 24    0 ]

The matrix ends up with all zeros in the last row:<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>0</mn><mi>y</mi><mo>=</mo><mn>0.</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]Thus, there are an infinite number of solutions and the system is classified as dependent. To find the generic solution, return to one of the original equations and solve for<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mo>.</mo></mrow></annotation-xml></semantics>[/itex]

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>3</mn><mi>x</mi><mo>+</mo><mn>4</mn><mi>y</mi><mo>=</mo><mn>12</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]          4y=12−3x            y=3− 3 4 x

So the solution to this system is<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] x,3− 3 4 x ).

Performing Row Operations on a 3×3 Augmented Matrix to Obtain Row-Echelon Form

Perform row operations on the given matrix to obtain row-echelon form.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 −3 4 2 −5 6 −3 3 4   |    3 6 6 ]

The first row already has a 1 in row 1, column 1. The next step is to multiply row 1 by<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>−2</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]and add it to row 2. Then replace row 2 with the result.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mn>−2</mn><msub/></mrow></annotation-xml></semantics>[/itex] R 1 + R 2 = R 2 →[ 1 −3 4 0 1 −2 −3 3 4 | 3 0 6 ]

Next, obtain a zero in row 3, column 1.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mn>3</mn><msub/></mrow></annotation-xml></semantics>[/itex] R 1 + R 3 = R 3 →[ 1 −3 4 0 1 −2 0 −6 16 | 3 0 15 ]

Next, obtain a zero in row 3, column 2.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mn>6</mn><msub/></mrow></annotation-xml></semantics>[/itex] R 2 + R 3 = R 3 →[ 1 −3 4 0 1 −2 0 0 4 | 3 0 15 ]

The last step is to obtain a 1 in row 3, column 3.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mfrac><mn>1</mn></mfrac></mrow></annotation-xml></semantics>[/itex] 2 R 3 = R 3 →[ 1 −3 4 0 1 −2 0 0 1   |    3 −6 21 2 ]

Write the system of equations in row-echelon form.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>  </mtext><mi>x</mi><mo>−</mo><mn>2</mn><mi>y</mi><mo>+</mo><mn>3</mn><mi>z</mi><mo>=</mo><mn>9</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]      −x+3y=−4 2x−5y+5z=17

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex]    1 − 5 2    5 2 ​  0   1 5  0  0   1   | 17 2 9 2 ]

# Solving a System of Linear Equations Using Matrices

We have seen how to write a system of equations with an augmented matrix, and then how to use row operations and back-substitution to obtain row-echelon form. Now, we will take row-echelon form a step farther to solve a 3 by 3 system of linear equations. The general idea is to eliminate all but one variable using row operations and then back-substitute to solve for the other variables.

Solving a System of Linear Equations Using Matrices

Solve the system of linear equations using matrices.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable><mtr><mtd><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>x</mi><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>−</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>y</mi><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo>+</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>z</mi><mo>=</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mn>8</mn></mrow></mtd></mtr></mtable></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 2x  +   3y   −   z=−2 3x   −   2y  −9z=    9

First, we write the augmented matrix.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 −1 1 2 3 −1 3 −2 −9    |    8 −2 9 ]

Next, we perform row operations to obtain row-echelon form.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="right"><mtr columnalign="right"><mtd columnalign="right"><mrow><mo>−</mo><mn>2</mn><msub/></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] R 1 + R 2 = R 2 →[ 1 −1 1 0 5 −3 3 −2 −9 | 8 −18 9 ] −3 R 1 + R 3 = R 3 →[ 1 −1 1 0 5 −3 0 1 −12 | 8 −18 −15 ]

The easiest way to obtain a 1 in row 2 of column 1 is to interchange<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><msub/></mrow></annotation-xml></semantics>[/itex] R 2  and<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><msub/></mrow></annotation-xml></semantics>[/itex] R 3 .

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext>Interchange</mtext><mtext> </mtext><msub/></mrow></annotation-xml></semantics>[/itex] R 2  and  R 3 →[ 1 −1 1 8 0 1 −12 −15 0 5 −3 −18 ]

Then

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mrow/></mtd></mtr><mtr><mtd><mtable columnalign="right"><mtr columnalign="right"><mtd columnalign="right"><mrow><mn>−5</mn><msub/></mrow></mtd></mtr></mtable></mtd></mtr></mtable></annotation-xml></semantics>[/itex] R 2 + R 3 = R 3 →[ 1 −1 1 0 1 −12 0 0 57 | 8 −15 57 ] − 1 57 R 3 = R 3 →[ 1 −1 1 0 1 −12 0 0 1 | 8 −15 1 ]

The last matrix represents the equivalent system.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mi>x</mi><mo>−</mo><mi>y</mi><mo>+</mo><mi>z</mi><mo>=</mo><mn>8</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]    y−12z=−15              z=1

Using back-substitution, we obtain the solution as<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 4,−3,1 ).

Solving a Dependent System of Linear Equations Using Matrices

Solve the following system of linear equations using matrices.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="right"><mtr columnalign="right"><mtd columnalign="right"><mrow><mo>−</mo><mi>x</mi><mn>−2</mn><mi>y</mi><mo>+</mo><mi>z</mi><mo>=</mo><mn>−1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]  2x+3y=2     y−2z=0

Write the augmented matrix.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] −1 −2 1 2 3 0 0 1 −2   |    −1 2 0 ]

First, multiply row 1 by<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>−1</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]to get a 1 in row 1, column 1. Then, perform row operations to obtain row-echelon form.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mo>−</mo><msub/></mrow></annotation-xml></semantics>[/itex] R 1 →[ 1 2 −1 1 2 3 0 2 0 1 −2 0 ]
<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><msub><mi>R</mi></msub></mrow></annotation-xml></semantics>[/itex] 2 ↔ R 3 →[ 1 2 −1 0 1 −2 2 3 0   | 1 0 2 ]
<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mn>−2</mn><msub/></mrow></annotation-xml></semantics>[/itex] R 1 + R 3 = R 3 →[ 1 2 −1 0 1 −2 0 −1 2 | 1 0 0 ]
<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><msub><mi>R</mi></msub></mrow></annotation-xml></semantics>[/itex] 2 + R 3 = R 3 →[ 1 2 −1 0 1 −2 0 0 0 | 2 1 0 ]

The last matrix represents the following system.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>−</mo><mi>z</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]         y−2z=0                0=0

We see by the identity<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>0</mn><mo>=</mo><mn>0</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]that this is a dependent system with an infinite number of solutions. We then find the generic solution. By solving the second equation for<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]and substituting it into the first equation we can solve for<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>z</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]in terms of<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>x</mi><mo>.</mo></mrow></annotation-xml></semantics>[/itex]

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>     </mtext><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>−</mo><mi>z</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]                    y=2z x+2(2z)−z=1            x+3z=1                     z= 1−x 3

Now we substitute the expression for<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>z</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]into the second equation to solve for<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]in terms of<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>x</mi><mo>.</mo></mrow></annotation-xml></semantics>[/itex]

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>          </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mi>y</mi><mo>−</mo><mn>2</mn><mi>z</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]                         z= 1−x 3  y−2( 1−x 3 )=0                        y= 2−2x 3

The generic solution is<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] x, 2−2x 3 , 1−x 3 ).

Solve the system using matrices.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable><mtr><mtd><mi>x</mi><mo>+</mo><mn>4</mn><mi>y</mi><mo>−</mo><mi>z</mi><mo>=</mo><mn>4</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 2x+5y+8z=15 x+3y−3z=1

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1,  1,  1 )

Can any system of linear equations be solved by Gaussian elimination?

Yes, a system of linear equations of any size can be solved by Gaussian elimination.

Given a system of equations, solve with matrices using a calculator.

1. Save the augmented matrix as a matrix variable<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mo stretchy="false">[</mo><mi>A</mi><mo stretchy="false">]</mo><mo>,</mo><mtext> </mtext><mo stretchy="false">[</mo><mi>B</mi><mo stretchy="false">]</mo><mo>,</mo><mtext> </mtext><mo stretchy="false">[</mo><mi>C</mi><mo stretchy="false">]</mo><mtext>,</mtext><mo> </mo><mo>…</mo><mo>.</mo></mrow></annotation-xml></semantics>[/itex]
2. Use the ref( function in the calculator, calling up each matrix variable as needed.
Solving Systems of Equations with Matrices Using a Calculator

Solve the system of equations.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="right"><mtr columnalign="right"><mtd columnalign="right"><mrow><mo> </mo><mn>5</mn><mi>x</mi><mo>+</mo><mn>3</mn><mi>y</mi><mo>+</mo><mn>9</mn><mi>z</mi><mo>=</mo><mn>−1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] −2x+3y−z=−2 −x−4y+5z=1

Write the augmented matrix for the system of equations.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 5 3 9 −2 3 −1 −1 −4 5   |    5 −2 −1 ]

On the matrix page of the calculator, enter the augmented matrix above as the matrix variable<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] A ].

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mo stretchy="false">[</mo><mi>A</mi><mo stretchy="false">]</mo><mo>=</mo><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 5 3 9 −1 −2 3 −1 −2 −1 −4 5 1 ]

Use the ref( function in the calculator, calling up the matrix variable<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] A ].

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext>ref</mtext><mo stretchy="false">(</mo><mo stretchy="false">[</mo><mi>A</mi><mo stretchy="false">]</mo><mo stretchy="false">)</mo></mrow></annotation-xml></semantics>[/itex]

Evaluate.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mrow><mo>[</mo></mrow></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 1    3 5    9 5 1 5 0   1    13 21 − 4 7 0   0   1 − 24 187 ]→ x+ 3 5 y+ 9 5 z=− 1 5        y+ 13 21 z=− 4 7                   z=− 24 187

Using back-substitution, the solution is<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 61 187 ,− 92 187 ,− 24 187 ).

Applying 2 × 2 Matrices to Finance

Carolyn invests a total of $12,000 in two municipal bonds, one paying 10.5% interest and the other paying 12% interest. The annual interest earned on the two investments last year was$1,335. How much was invested at each rate?

We have a system of two equations in two variables. Let<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>x</mi><mo>=</mo><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]the amount invested at 10.5% interest, and<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mo>=</mo><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]the amount invested at 12% interest.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>                </mtext><mi>x</mi><mo>+</mo><mi>y</mi><mo>=</mo><mn>12,000</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 0.105x+0.12y=1,335

As a matrix, we have

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 1 0.105 0.12   |    12,000 1,335 ]

Multiply row 1 by<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>−0.105</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]and add the result to row 2.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 1 0 0.015   |    12,000 75 ]

Then,

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>0.015</mn><mi>y</mi><mo>=</mo><mn>75</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]          y=5,000

So<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mn>12,000</mn><mn>−5,000</mn><mo>=</mo><mn>7,000.</mn></mrow></annotation-xml></semantics>[/itex]

Thus, $5,000 was invested at 12% interest and$7,000 at 10.5% interest.

Applying 3 × 3 Matrices to Finance

Ava invests a total of $10,000 in three accounts, one paying 5% interest, another paying 8% interest, and the third paying 9% interest. The annual interest earned on the three investments last year was$770. The amount invested at 9% was twice the amount invested at 5%. How much was invested at each rate?

We have a system of three equations in three variables. Let<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>x</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]be the amount invested at 5% interest, let<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]be the amount invested at 8% interest, and let<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>z</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]be the amount invested at 9% interest. Thus,

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>                     </mtext><mi>x</mi><mo>+</mo><mi>y</mi><mo>+</mo><mi>z</mi><mo>=</mo><mn>10</mn><mo>,</mo><mn>000</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 0.05x+0.08y+0.09z=770                          2x−z=0

As a matrix, we have

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 1 1 0.05 0.08 0.09 2 0 −1   |    10,000 770 0 ]

Now, we perform Gaussian elimination to achieve row-echelon form.

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>−0.05</mn><msub/></mrow></mtd></mtr></mtable></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] R 1 + R 2 = R 2 →[ 1 1 1 0 0.03 0.04 2 0 −1 | 10,000 270 0 ]       −2 R 1 + R 3 = R 3 →[ 1 1 1 0 0.03 0.04 0 −2 −3 |10,000 270 −20,000 ]                  1 0.03 R 2 = R 2 →[ 0 1 1 0 1 4 3 0 −2 −3 | 10,000 9,000 −20,000 ]             2 R 2 + R 3 = R 3→[ 1 1 1 0 1 4 3 0 0 − 1 3 | 10,000 9,000 −2,000 ]

The third row tells us<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mo>−</mo><mfrac/></mrow></annotation-xml></semantics>[/itex] 1 3 z=−2,000; thus<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>z</mi><mo>=</mo><mn>6,000.</mn></mrow></annotation-xml></semantics>[/itex]

The second row tells us<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mo>+</mo><mfrac/></mrow></annotation-xml></semantics>[/itex] 4 3 z=9,000. Substituting<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>z</mi><mo>=</mo><mn>6,000</mn><mo>,</mo></mrow></annotation-xml></semantics>[/itex]we get

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="right"><mtr columnalign="right"><mtd columnalign="right"><mrow><mi>y</mi><mo>+</mo><mfrac/></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 4 3 (6,000)=9,000 y+8,000=9,000 y=1,000

The first row tells us<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>x</mi><mo>+</mo><mi>y</mi><mo>+</mo><mi>z</mi><mo>=</mo><mn>10</mn><mo>,</mo><mn>000.</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]Substituting<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>y</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>000</mn><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]and<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>z</mi><mo>=</mo><mn>6</mn><mo>,</mo><mn>000</mn><mo>,</mo></mrow></annotation-xml></semantics>[/itex]we get

<math display="block" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>x</mi><mo>+</mo><mn>1</mn><mo>,</mo><mn>000</mn><mo>+</mo><mn>6</mn><mo>,</mo><mn>000</mn><mo>=</mo><mn>10,000</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]                               x=3,000

The answer is $3,000 invested at 5% interest,$1,000 invested at 8%, and $6,000 invested at 9% interest. A small shoe company took out a loan of$1,500,000 to expand their inventory. Part of the money was borrowed at 7%, part was borrowed at 8%, and part was borrowed at 10%. The amount borrowed at 10% was four times the amount borrowed at 7%, and the annual interest on all three loans was $130,500. Use matrices to find the amount borrowed at each rate.$150,000 at 7%, $750,000 at 8%,$600,000 at 10%

Access these online resources for additional instruction and practice with solving systems of linear equations using Gaussian elimination.

# Key Concepts

• An augmented matrix is one that contains the coefficients and constants of a system of equations. See [link].
• A matrix augmented with the constant column can be represented as the original system of equations. See [link].
• Row operations include multiplying a row by a constant, adding one row to another row, and interchanging rows.
• We can use Gaussian elimination to solve a system of equations. See [link], [link], and [link].
• Row operations are performed on matrices to obtain row-echelon form. See [link].
• To solve a system of equations, write it in augmented matrix form. Perform row operations to obtain row-echelon form. Back-substitute to find the solutions. See [link] and [link].
• A calculator can be used to solve systems of equations using matrices. See [link].
• Many real-world problems can be solved using augmented matrices. See [link] and [link].

# Section Exercises

## Verbal

Can any system of linear equations be written as an augmented matrix? Explain why or why not. Explain how to write that augmented matrix.

Yes. For each row, the coefficients of the variables are written across the corresponding row, and a vertical bar is placed; then the constants are placed to the right of the vertical bar.

Can any matrix be written as a system of linear equations? Explain why or why not. Explain how to write that system of equations.

Is there only one correct method of using row operations on a matrix? Try to explain two different row operations possible to solve the augmented matrix<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 9 3 1 −2   |    0 6 ].

No, there are numerous correct methods of using row operations on a matrix. Two possible ways are the following: (1) Interchange rows 1 and 2. Then<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><msub/></mrow></annotation-xml></semantics>[/itex] R 2 = R 2 −9 R 1 . (2)<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><msub/></mrow></annotation-xml></semantics>[/itex] R 2 = R 1 −9 R 2 . Then divide row 1 by 9.

Can a matrix whose entry is 0 on the diagonal be solved? Explain why or why not. What would you do to remedy the situation?

Can a matrix that has 0 entries for an entire row have one solution? Explain why or why not.

No. A matrix with 0 entries for an entire row would have either zero or infinitely many solutions.

## Algebraic

For the following exercises, write the augmented matrix for the linear system.

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>8</mn><mi>x</mi><mn>−37</mn><mi>y</mi><mo>=</mo><mn>8</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 2x+12y=3

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mtext>   </mtext><mn>16</mn><mi>y</mi><mo>=</mo><mn>4</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 9x−y=2

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 0 16 9 −1 | 4 2 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mn>3</mn><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>+</mo><mn>10</mn><mi>z</mi><mo>=</mo><mn>3</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] −6x+2y+5z=13              4x+z=18

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>  </mtext><mi>x</mi><mo>+</mo><mn>5</mn><mi>y</mi><mo>+</mo><mn>8</mn><mi>z</mi><mo>=</mo><mn>19</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]      12x+3y=4 3x+4y+9z=−7

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 5 8 12 3 0 3 4 9 | 16 4 −7 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>6</mn><mi>x</mi><mo>+</mo><mn>12</mn><mi>y</mi><mo>+</mo><mn>16</mn><mi>z</mi><mo>=</mo><mn>4</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]   19x−5y+3z=−9              x+2y=−8

For the following exercises, write the linear system from the augmented matrix.

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] −2 5 6 −18   |    5 26 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>−2</mn><mi>x</mi><mo>+</mo><mn>5</mn><mi>y</mi><mo>=</mo><mn>5</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 6x−18y=26

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 4 10 17   |    10 439 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3 2 0 −1 −9 4 8 5 7   |    3 −1 8 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>3</mn><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>=</mo><mn>13</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] −x−9y+4z=53 8x+5y+7z=80

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 8 29 1 −1 7 5 0 0 3   |    43 38 10 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 4 5 −2 0 1 58 8 7 −3   |    12 2 −5 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>4</mn><mi>x</mi><mo>+</mo><mn>5</mn><mi>y</mi><mn>−2</mn><mi>z</mi><mo>=</mo><mn>12</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]         y+58z=2 8x+7y−3z=−5

For the following exercises, solve the system by Gaussian elimination.

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 0 0 0   |    3 0 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 0 1 0   |    1 2 ]

No solutions

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 2 4 5   |    3 6 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] −1 2 4 −5   |    −3 6 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mo stretchy="false">(</mo><mn>−1</mn><mo>,</mo><mn>−2</mn><mo stretchy="false">)</mo></mrow></annotation-xml></semantics>[/itex]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] −2 0 0 2   |    1 −1 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>  </mtext><mn>2</mn><mi>x</mi><mo>−</mo><mn>3</mn><mi>y</mi><mo>=</mo><mo>−</mo><mn>9</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 5x+4y=58

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 6,7 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>6</mn><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>=</mo><mn>−4</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 3x+4y=−17

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>2</mn><mi>x</mi><mo>+</mo><mn>3</mn><mi>y</mi><mo>=</mo><mn>12</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]   4x+y=14

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3,2 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>−4</mn><mi>x</mi><mn>−3</mn><mi>y</mi><mo>=</mo><mn>−2</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]    3x−5y=−13

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>−5</mn><mi>x</mi><mo>+</mo><mn>8</mn><mi>y</mi><mo>=</mo><mn>3</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]  10x+6y=5

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 5 , 1 2 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mtext>  </mtext><mn>3</mn><mi>x</mi><mo>+</mo><mn>4</mn><mi>y</mi><mo>=</mo><mn>12</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] −6x−8y=−24

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>−60</mn><mi>x</mi><mo>+</mo><mn>45</mn><mi>y</mi><mo>=</mo><mn>12</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]   20x−15y=−4

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] x, 4 15 (5x+1) )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mn>11</mn><mi>x</mi><mo>+</mo><mn>10</mn><mi>y</mi><mo>=</mo><mn>43</mn></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 15x+20y=65

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mn>2</mn><mi>x</mi><mo>−</mo><mi>y</mi><mo>=</mo><mn>2</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 3x+2y=17

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 3,4 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mtable columnalign="left"><mtr><mtd><mrow/></mtd></mtr><mtr><mtd><mn>−1.06</mn><mi>x</mi><mn>−2.25</mn><mi>y</mi><mo>=</mo><mn>5.51</mn></mtd></mtr></mtable></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] −5.03x−1.08y=5.40

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mfrac><mn>3</mn></mfrac></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 4 x− 3 5 y=4 1 4 x+ 2 3 y=1

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 196 39 ,− 5 13 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mfrac><mn>1</mn></mfrac></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 4 x− 2 3 y=−1 1 2 x+ 1 3 y=3

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 0 0 0 1 1 0 0 1   |    31 45 87 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 31,−42,87 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 0 1 1 1 0 0 1 1   |    50 20 −90 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1 2 3 0 5 6 0 0 8   |    4 7 9 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 21 40 , 1 20 , 9 8 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>[</mo></mrow></mrow></annotation-xml></semantics>[/itex] −0.1 0.3 −0.1 −0.4 0.2 0.1 0.6 0.1 0.7   |    0.2 0.8 −0.8 ]

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mn>−2</mn><mi>x</mi><mo>+</mo><mn>3</mn><mi>y</mi><mo>−</mo><mn>2</mn><mi>z</mi><mo>=</mo><mn>3</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]       4x+2y−z=9      4x−8y+2z=−6

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 18 13 , 15 13 ,− 15 13 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>      </mtext><mi>x</mi><mo>+</mo><mi>y</mi><mo>−</mo><mn>4</mn><mi>z</mi><mo>=</mo><mn>−4</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]   5x−3y−2z=0   2x+6y+7z=30

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>      </mtext><mn>2</mn><mi>x</mi><mo>+</mo><mn>3</mn><mi>y</mi><mo>+</mo><mn>2</mn><mi>z</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]   −4x−6y−4z=−2  10x+15y+10z=5

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] x,y, 1 2 (1−2x−3y) )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>    </mtext><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>−</mo><mi>z</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] −x−2y+2z=−2 3x+6y−3z=5

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext> </mtext><mtext>    </mtext><mi>x</mi><mo>+</mo><mn>2</mn><mi>y</mi><mo>−</mo><mi>z</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] −x−2y+2z=−2  3x+6y−3z=3

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] x,− x 2 ,−1 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>​</mtext><mtext>   </mtext><mi>x</mi><mo>+</mo><mi>y</mi><mo>=</mo><mn>2</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]    x+z=1 −y−z=−3

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>x</mi><mo>+</mo><mi>y</mi><mo>+</mo><mi>z</mi><mo>=</mo><mn>100</mn></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex]     x+2z=125 −y+2z=25

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 125,−25,0 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mtable columnalign="left"><mtr><mtd><mfrac><mn>1</mn></mfrac></mtd></mtr></mtable></annotation-xml></semantics>[/itex] 4 x− 2 3 z=− 1 2 1 5 x+ 1 3 y= 4 7 1 5 y− 1 3 z= 2 9

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>−</mo><mfrac/></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 1 2 x+ 1 2 y+ 1 7 z=− 53 14     1 2 x− 1 2 y+ 1 4 z=3      1 4 x+ 1 5 y+ 1 3 z= 23 15

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 8,1,−2 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>−</mo><mfrac/></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 1 2 x− 1 3 y+ 1 4 z=− 29 6     1 5 x+ 1 6 y− 1 7 z= 431 210 − 1 8 x+ 1 9 y+ 1 10 z=− 49 45

## Extensions

For the following exercises, use Gaussian elimination to solve the system.

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mfrac><mrow><mi>x</mi><mn>−1</mn></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 7 + y−2 8 + z−3 4 =0                      x+y+z=6       x+2 3 +2y+ z−3 3 =5

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] 1,2,3 )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mfrac><mrow><mi>x</mi><mn>−1</mn></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 4 − y+1 4 +3z=−1    x+5 2 + y+7 4 −z=4          x+y− z−2 2 =1

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>      </mtext><mfrac/></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] x−3 4 − y−1 3 +2z=−1 x+5 2 + y+5 2 + z+5 2 =8                      x+y+z=1

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mrow><mo>(</mo></mrow></mrow></annotation-xml></semantics>[/itex] x, 31 28 − 3x 4 , 1 28 (−7x−3) )

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mfrac><mrow><mi>x</mi><mn>−3</mn></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] 10 + y+3 2 −2z=3   x+5 4 − y−1 8 +z= 3 2 x−1 4 + y+4 2 +3z= 3 2

<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>     </mtext><mfrac/></mrow></mtd></mtr></mtable></mrow></annotation-xml></semantics>[/itex] x−3 4 − y−1 3 +2z=−1 x+5 2 + y+5 2 + z+5 2 =7                      x+y+z=1

No solutions exist.

## Real-World Applications

For the following exercises, set up the augmented matrix that describes the situation, and solve for the desired solution.

Every day, a cupcake store sells 5,000 cupcakes in chocolate and vanilla flavors. If the chocolate flavor is 3 times as popular as the vanilla flavor, how many of each cupcake sell per day?

At a competing cupcake store, $4,520 worth of cupcakes are sold daily. The chocolate cupcakes cost$2.25 and the red velvet cupcakes cost $1.75. If the total number of cupcakes sold per day is 2,200, how many of each flavor are sold each day? 860 red velvet, 1,340 chocolate You invested$10,000 into two accounts: one that has simple 3% interest, the other with 2.5% interest. If your total interest payment after one year was $283.50, how much was in each account after the year passed? You invested$2,300 into account 1, and $2,700 into account 2. If the total amount of interest after one year is$254, and account 2 has 1.5 times the interest rate of account 1, what are the interest rates? Assume simple interest rates.

4% for account 1, 6% for account 2

Bikes’R’Us manufactures bikes, which sell for $250. It costs the manufacturer$180 per bike, plus a startup fee of $3,500. After how many bikes sold will the manufacturer break even? A major appliance store is considering purchasing vacuums from a small manufacturer. The store would be able to purchase the vacuums for$86 each, with a delivery fee of $9,200, regardless of how many vacuums are sold. If the store needs to start seeing a profit after 230 units are sold, how much should they charge for the vacuums?$126

The three most popular ice cream flavors are chocolate, strawberry, and vanilla, comprising 83% of the flavors sold at an ice cream shop. If vanilla sells 1% more than twice strawberry, and chocolate sells 11% more than vanilla, how much of the total ice cream consumption are the vanilla, chocolate, and strawberry flavors?

At an ice cream shop, three flavors are increasing in demand. Last year, banana, pumpkin, and rocky road ice cream made up 12% of total ice cream sales. This year, the same three ice creams made up 16.9% of ice cream sales. The rocky road sales doubled, the banana sales increased by 50%, and the pumpkin sales increased by 20%. If the rocky road ice cream had one less percent of sales than the banana ice cream, find out the percentage of ice cream sales each individual ice cream made last year.

Banana was 3%, pumpkin was 7%, and rocky road was 2%

A bag of mixed nuts contains cashews, pistachios, and almonds. There are 1,000 total nuts in the bag, and there are 100 less almonds than pistachios. The cashews weigh 3 g, pistachios weigh 4 g, and almonds weigh 5 g. If the bag weighs 3.7 kg, find out how many of each type of nut is in the bag.

A bag of mixed nuts contains cashews, pistachios, and almonds. Originally there were 900 nuts in the bag. 30% of the almonds, 20% of the cashews, and 10% of the pistachios were eaten, and now there are 770 nuts left in the bag. Originally, there were 100 more cashews than almonds. Figure out how many of each type of nut was in the bag to begin with.

100 almonds, 200 cashews, 600 pistachios

## Glossary

augmented matrix
a coefficient matrix adjoined with the constant column separated by a vertical line within the matrix brackets
coefficient matrix
a matrix that contains only the coefficients from a system of equations
Gaussian elimination
using elementary row operations to obtain a matrix in row-echelon form
main diagonal
entries from the upper left corner diagonally to the lower right corner of a square matrix
row-echelon form
after performing row operations, the matrix form that contains ones down the main diagonal and zeros at every space below the diagonal
row-equivalent
two matrices<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>A</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]and<math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><annotation-xml encoding="MathML-Content"><mrow><mtext> </mtext><mi>B</mi><mtext> </mtext></mrow></annotation-xml></semantics>[/itex]are row-equivalent if one can be obtained from the other by performing basic row operations
row operations
adding one row to another row, multiplying a row by a constant, interchanging rows, and so on, with the goal of achieving row-echelon form