1. Why Linear Algebra Is So Important in Modern Engineering Universities?
Linear algebra originated in 18th Century. Until 1950’s, it was only studied by mathematics students and some graduate students of other areas. But after 1960, it was included in the undergraduate curriculum of the non Mathematics Departments of American Universities. In China, the engineering departments offered no linear algebra until 1980s. Why is linear algebra so popular in recent decades? The most fundamental reason is that it is the basis of many algorithms for the machine implementation of numerous technical problems. In Prof. Leontiff’s work leading to the 1973 economics Nobel Prize, he used computer to solve the 42-th order of linear equation group. In early 80’s, linear algebra software package LINPACK was developed. So users can solve the problem of large complex linear algebra without mastering the mathematical details of matrix theory. LINPACK is also the background support of MATLAB. The achievements of linear algebra in the last half of the 20th century also focused on solving the problems of large-scale matrix using computer. It is no longer a pure theory of mathematics but becomes a calculation tool for common engineering students.
Modern society has put forward higher requirements for engineering mathematics. In 1989, Dr. Tsien Hsueshen wrote: “In the future, solving a problem can all rely on the computer. People do not need to calculate one point after another. But until today, the mathematics of engineering is built on the calculation by people themselves, ... The mathematics course of science and engineering must be reformed. It should not focus on solving problems by students themselves, but let them learn to command machine in solving problems and to understand the answers given by the computer...” . Linear algebra is the basis of numerical calculation, which not only takes the lead in the use of computer, but also promotes the computer applications based on other mathematic courses and engineering problems.
The importance of linear algebra is also reflected in the increasing number of new fields using multivariable matrix models. Leontiff’s story was a good example. This kind of “demand pull” for non-mathematics major is in the use of computer to solve high-order complex matrix model, not in the manual deduction of mathematical formula. Its contents are very different to the classical linear algebra. At the beginning, no one had a clear idea, the curriculum of mathematics department was copied to teach the non-mathematics students. But the task for non mathematical students is to solve real world problems. What they face is high order, numerical complex matrices models. Their different requirement can be summed up in Table 1.
After implementation for decades, some famous professors in the USA became seriously concerned. They set up a Linear Algebra Curriculum Study Group (LACSG) in 1990. Five recommendations of linear algebra as the first course were put forward :
1) It is to meet the needs of non-mathematics students;
2) It should be oriented to matrix operation instead of vector space;
3) It should be adapted to the level and needs of students (freshman or sophmore);
4) It should use the latest computing technology;
Since 2009, a project: “Using MATLAB and Modeling to Reform Engineering Linear Algebra Course” was set up by the Ministry of Higher Education in China. It was proposed and hosted by Xidian University. 19 universities participated in and implemented for two years. It involved 200 professors and 45,000 students. Now it is a consensus among these universities to require their engineering freshmen to learn MATLAB and the use of computers. This laid the first stone of computer usage in the undergraduate curriculum.
Under the constraint of credit hour to the course of linear algebra, there is no sufficient time to cover adequately both the classical contents of linear algebra and its MATLAB programming. Thus it is of practical importance to determine the minimum part of linear algebra essential to engineering implementation.
2. How to Meet the Needs of Non Mathematics Majors?
To realize this, one must go deep into the practice of non-mathematics majors. We found that matrix equation is the core. The author had programmed more than 200 application in a dozen or more subsequent courses, this book took 46 cases among them. These including the following: (1) Matrix calculation and inversion in 12 cases (highest order 21); (2) Solving linear equations 16 cases (highest order 7); (3) Vector spatial relationship and motion transform 8 cases, (involving spatial geometry and motion, coplanarity problems); (4) linear transformation and characteristic roots 10 cases (containing complex eigenvalue problems). Other than some higher-order problems, the 46 cases in this book cover practically most application types in undergraduate courses.
Table 1. Different features of linear algebra for math and non-math major.
Method of reverse thinking is used to extract from these examples the minimum required linear algebra theory. All theories needed by these cases were kept and explain in detail. Any contents not directly used are removed. Some complicated theory is simplified if possible. After all, engineering talent is the use of mathematics, not engaged in mathematics itself. This book takes full advantage of MATLAB resources to link up theory with the subroutines, so that students will understand what and why they are doing.
The main points are as follows:
(1) Gaussian elimination method for solving linear equations is the core theory of the book. Determinant is simplified and defined as the product of the pivot elements in the row echelon form . Rank, vector correlation, ... are defined accordingly. Under this definition, subroutine for row echelon form is often used. It can replace almost all determinant calculations.
(2) In order to integrate with computer in large scale numerical computations, students must have the concepts of computational complexity and accuracy which are never seen in the previous textbooks. The computational feasibility knock out two classical definitions of determinant in engineering. Error dependent rank is the key for solving coplanarity of vector group, it cannot be discovered by manual computation with integer elements. These examples show that our new concerns prevail the old ones.
(3) At the freshman level, we emphasized the relationship between linear algebra and 3-D solid analytic geometry. It is not necessary to teach N-D vector space before they clearly grasp 3-D space. Besides we cannot find any application of N-D in the undergraduate engineering courses.
Thus the contents and organization of this book are very different from the classical linear algebra textbooks.
3. The Reform of Determinant Definition
For the computation of determinant, the definition and method of its realization must be paid attention to. In numerical calculation, the computational complexity is usually measured by the number of multiplications. Table 2 gives a comparison of three methods with different matrix order:
Comparing method 1 and 3, we see that only when n = 2, 3, method 1 is better. When n ≥ 4, method 1 and 2 are unable to calculate manually. When n = 25, the multiplication number of these two definitions are too large that it would take hundreds of years to calculate a single determinant using a top computer today. This phenomenon is called “Dimensionality Curse”. Clearly definitions 1 and 2 are useless to high order systems and can be ignored in textbook. For definition 3, not only the amount of calculation is greatly reduced, the formula is also very easy to program. Concepts like number of inversion, cofactors and its expansion, adjoint matrix, ... that discussed in the classical theory can all be omitted. Avoiding these "stumbling block" can greatly reduce the difficulty for the non-mathematics students and save many hours of instruction.
We found that in all the application cases, none of them really need determinant (except area and volume calculation which is not linear algebra problem). Many theorems on determinants and its nature are useless. This is because Gauss elimination method has already used nonzero criteria for pivot elements, which equivalents to determinant criteria. This means that determinant has already been used in Gaussian elimination, its further calculation is redundant. Pivot elements are the natural extension of Gauss elimination, no new concept or formula is needed. By using this definition, it is also easy to prove a variety of useful properties of determinant.
Someone takes the classical definition of determinant as the basis of linear algebra. This is not true. After proved that it has the fatal defect for computing high order matrix, we should not let the engineering students waste time on it. The era of the classical determinant is over! One of the supporting evidence is the software tools. In MATLAB, you cannot find terms or subroutines based upon these classical concepts. This means that real world do not need these concepts. The second evidence is some well known textbooks  . In , Lay made a brilliant comment: “In Cauchy’s day, when life was simple and matrices were small, determinant played a major role in mathematics. Today determinant is of little numerical value in large scale matrix computations.”
Table 2. Multiplications needed for different determinant definitions.
4. Give a Thorough Exposition to 3-D Vector Space
5. Weaken the Underdetermined, Strengthen the Overdetermined
We found that underdetermined applications are rarely seen, because they are mostly due to defect in the problem proposition. Engineer should either reject it or introduce supplement conditions to make it determined. We cannot find any usage of underdetermined solution set in real world as taught by classical linear algebra. For example, balancing chemical equations remains underdetermined unless you add a condition that the solution must be minimum positive integer numbers. It would be funny that you give the solution set as a super line in a virtual space instead of coefficients of each substances.
The over determined problems come from the inevitable interference and measurement errors in practice. They occur in many engineering problems. Finding the least mean square solution using geometric method involves the concept of vector space. But most of the textbook in China do not contain it. So many engineering problems are unsolvable to them. Over emphasizing underdetermined system but ignoring over determined problem is the major weakness. This caused linear algebra go far away from engineering. We should take a balanced way.
6. Reduce the Order of Characteristic Equation but Allow Complex Roots
In the classical textbooks, the characteristic equations are discussed in N-th order, but no book tells how to find the characteristic roots of even 3rd order by hand calculation. This is apparently inconsistent. Moreover, we cannot find their physical meaning and geometric model in undergraduate courses. So we restrict characteristic equations and quadratic forms to 2nd order, but included complex eigenvalues which can help to understand vibration problems in engineering.
7. Brief Introduction to Some Examples That Demand for Course Reform 
The following examples can prove the helplessness and necessity of reforming for classical theory.
Example 1. Circuit as Figure 1, find Ia, Ib, Ic.
The matrix model is:
In real world circuits, the matrix order could be dozens. In AC circuit, the variables I and coefficients R are complex numbers. It is impossible to solve the problem without computer software, while determinant calculation is not needed. We can get the solution easily without knowing anything about its determinant.
Example 2. Signal Flow Graph as Figure 2, find transfer function w4 = x4/u.
The equations and matrix model are as follows:
Figure 1. A resistance circuit.
Figure 2. A signal flow graph.
The result given by MATLAB execution is:
This method can also be used in digital filter (Figure 3) and control system (Figure 4) analysis to find out the overall system function. No matter how high its order is, using matrix operation supported software tools can give the solution in seconds.
The answer of Figure 3 is:
The answer of Figure 4 is :
This system has a 2 × 1 input vector and a 5 × 1 output vector, we obtain 5 × 2 transfer functions in one execution. This is a very important usage of linear algebra in subsequent courses. Using hand calculation could never get this privilege. These examples come from 3 subsequent important courses, so we can see the importance and power of linear algebra after reform to engineers.
Example 3. Evaluate coplanarity of 5 spatial points listed in Table 3.
This is a problem that can be encountered in modern mechanical workshop with coordinate measuring machine. The steps for solving this problem are as follows:
Figure 3. A digital filter.
Figure 4. A control system.
Table 3. Coordinates measured for the 5 points.
(1) Given vector group V = [va, vb, vc, vd, ve], its columns are the coordinates of these points. Let U = V − [ve, ve, ve, ve, ve] be the difference vector group, which means the coordinate origin is moved to point e.
(2) if the points are in a plane, the rank of U must be 2 or less. Because of the limited data precision and computing error, precise plane can never realized. You will always find rank(U) = 3.
(3) The question is that in what extent these points are on the same plane. What is the biggest error? This relates to singular value decomposition which is too difficult. But it can be solved by software tools for non mathematic students. In MATLAB, function RANK (A, TOL) contains an tolerant argument TOL, its default value is 10-16. So very small distortion of the plane will cause rank(U) = 3. Using bigger TOL, we obtain rank (U, 0.01) = 3 and rank(U,0.1) = 2. This means that these 5 points are about on a plane within an error of ±0.1.
(4) To find the equation of the fitted plane and the distance of each point to it, the easiest way is to use qr decomposition. Inputting [Q,R] = qr(U), we get:
The third column of Q is the direction cosine of the normal line of the plane. Numbers on the third row of R are the distances of each points to the fitted plane. Its equation can be obtained using point-slope form of plane:
In solving this kind of problem, we can see that the 3-D concept and computer software are both extremely important. Hand calculation is useless and should be replaced.
Using these reform measures, the amount of theory is greatly reduced, useful theory has increased, the contents of theorem are significantly reduced; but the ability of solving practical problems has been considerably boosted. This book has set up a minimum theoretical limit to meet the needs of the engineering application objectively. Each content has its application background. Removing any part inevitable result in failure in solving some practical problems. The implementation of this book, plus MATLAB, need about 30 class hours, suitable for low class syllabus and continuing education. This book is also beneficial to the high hours curriculum. Because most of the existing textbooks are belong to the classical mathematical system. They cannot solve most of the application problems in this book. We can suggest this book as a reference in solving actual problems. All the textbooks can use these problems as the minimum engineering requirements, then use the spare time in their syllabus for higher requirements.
Postgraduate enrollment examination is an important factor affecting the reform of this course. Many exam problems were coming from classical linear algebra. Obviously the examination should be reformed to reduce useless theory and include some computer applications. Before the examination reform, we can leave these classical content into an follow-up elective course for senior students who want to take the examination.
Appendix: Define Determinant by Product of Diagonal Pivots  
While solving matrix equation AX = b using Gauss elimination, the augmented matrix C = [A, b] is converted to a diagonal matrix P and a column vector d as follows：
If all the pivots p11, p22, ..., pnn on the diagonal do not equal to zero, the solution of the equation is existed.
This condition can be expressed in another way: Only if, the solution is existed. The product
of diagonal pivots D is defined as determinant of matrix A. Here we discard its definition of sign to save paper. Actually we only want to know whether it equals to zero or not. Its exact value and sign are all unimportant. In most real world problems, people only want the solution, none requires determinant. So this definition is enough for most non math majors. In some cases like example-3, when some deeper concept is needed, classical theory won’t help. The only way is to search subroutines in MATLAB.