Example 5.  Balance the propane-oxygen equation by solve the homogeneous linear system

        [Graphics:Images/HomogenLinSystemMod_gr_372.gif]

Solution 5.

Enter the equations into Mathematica.  

[Graphics:../Images/HomogenLinSystemMod_gr_373.gif]

Identify the matrix of coefficients A and column vector B for the matrix problem AX = B.  

[Graphics:../Images/HomogenLinSystemMod_gr_374.gif]


[Graphics:../Images/HomogenLinSystemMod_gr_375.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_376.gif]

 

[Graphics:../Images/HomogenLinSystemMod_gr_377.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_378.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_379.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_380.gif]


Form the augmented matrix  M = [A, B]  and perform Gauss-Jordan elimination with row interchanges.

[Graphics:../Images/HomogenLinSystemMod_gr_381.gif]


[Graphics:../Images/HomogenLinSystemMod_gr_382.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_383.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_384.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_385.gif]

Find the reduced row echelon form of the augmented matrix  M = [A, B].  

[Graphics:../Images/HomogenLinSystemMod_gr_386.gif]


[Graphics:../Images/HomogenLinSystemMod_gr_387.gif]
[Graphics:../Images/HomogenLinSystemMod_gr_388.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_389.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_390.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_391.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_392.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_393.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_394.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_395.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_396.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_397.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_398.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_399.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_400.gif]

This linear system is equivalent to:

    [Graphics:../Images/HomogenLinSystemMod_gr_401.gif]

There is one free variable which we choose to be  [Graphics:../Images/HomogenLinSystemMod_gr_402.gif].  It is used in computing  [Graphics:../Images/HomogenLinSystemMod_gr_403.gif].  

Solve the previous equations for  [Graphics:../Images/HomogenLinSystemMod_gr_404.gif].  

    [Graphics:../Images/HomogenLinSystemMod_gr_405.gif]  

Make the substitution  [Graphics:../Images/HomogenLinSystemMod_gr_406.gif].

    [Graphics:../Images/HomogenLinSystemMod_gr_407.gif]  

The solution vector  [Graphics:../Images/HomogenLinSystemMod_gr_408.gif]  is

[Graphics:../Images/HomogenLinSystemMod_gr_409.gif]

We are done.

Aside.  We can verify that this is the solution by direct multiplication A X.  This is just for fun !

[Graphics:../Images/HomogenLinSystemMod_gr_410.gif]


[Graphics:../Images/HomogenLinSystemMod_gr_411.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_412.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_413.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_414.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_415.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_416.gif]

Aside.  We can let Mathematica find the reduced row echelon matrix.  This is just for fun !

[Graphics:../Images/HomogenLinSystemMod_gr_417.gif]


[Graphics:../Images/HomogenLinSystemMod_gr_418.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_419.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_420.gif]

Notice.  Since the last row is entirely zero, the system has reduced to three equations and four unknowns.  

We can add the equation  [Graphics:../Images/HomogenLinSystemMod_gr_421.gif]  to those in the reduced row echelon form and then row reduce one more time to get the solution.

[Graphics:../Images/HomogenLinSystemMod_gr_422.gif]


[Graphics:../Images/HomogenLinSystemMod_gr_423.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_424.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_425.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_426.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_427.gif]

The 4×4 identity matrix appears in the left 3 columns of  M, and the given linear system is equivalent to:  

    [Graphics:../Images/HomogenLinSystemMod_gr_428.gif]

The solution vector is the fourth column of  M.

[Graphics:../Images/HomogenLinSystemMod_gr_429.gif]


[Graphics:../Images/HomogenLinSystemMod_gr_430.gif]


We can verify that this is the solution by direct multiplication A X.  This is just for fun !

[Graphics:../Images/HomogenLinSystemMod_gr_431.gif]



[Graphics:../Images/HomogenLinSystemMod_gr_432.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_433.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_434.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_435.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_436.gif]

[Graphics:../Images/HomogenLinSystemMod_gr_437.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004