International Journal of Engineering Mathematics

Volume 2016 (2016), Article ID 9382739, 14 pages

http://dx.doi.org/10.1155/2016/9382739

## On the Extension of Sarrus’ Rule to Matrices: Development of New Method for the Computation of the Determinant of Matrix

Department of Mechanical Engineering, University of Lagos, Lagos, Nigeria

Received 14 June 2016; Revised 8 August 2016; Accepted 30 August 2016

Academic Editor: Giuseppe Carbone

Copyright © 2016 M. G. Sobamowo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The determinant of a matrix is very powerful tool that helps in establishing properties of matrices. Indisputably, its importance in various engineering and applied science problems has made it a mathematical area of increasing significance. From developed and existing methods of finding determinant of a matrix, basketweave method/Sarrus’ rule has been shown to be the simplest, easiest, very fast, accurate, and straightforward method for the computation of the determinant of 3 × 3 matrices. However, its gross limitation is that this method/rule does not work for matrices larger than 3 × 3 and this fact is well established in literatures. Therefore, the state-of-the-art methods for finding the determinants of 4 × 4 matrix and larger matrices are predominantly founded on non-basketweave method/non-Sarrus’ rule. In this work, extension of the simple, easy, accurate, and straightforward approach to the determinant of larger matrices is presented. The paper presents the developments of new method with different schemes based on the basketweave method/Sarrus’ rule for the computation of the determinant of 4 × 4. The potency of the new method is revealed in generalization of the basketweave method/non-Sarrus’ rule for the computation of the determinant of () matrices. The new method is very efficient, very consistence for handy calculations, highly accurate, and fastest compared to other existing methods.

#### 1. Introduction

Over the years, the subject, linear algebra has been shown to be the most fundamental component in mathematics as it presents powerful tools in wide varieties of areas from theoretical science to engineering, including computer science. Its important role and abilities in solving real life problems and in data clarification [1] have led it to be frequently applied in all the branches of science, engineering, social science, and management. During the applications and analysis in such areas of studies, a system of linear equations can be written in matrix form and solving the system of linear equations and the inversion of matrices is necessary which is mainly dependent on determinant (a real number or a function of the elements of an matrix that yields a single number that well determines something about the matrix). Therefore, the importance of finding the determinant in linear algebra cannot be overemphasized as it does not only help in finding solution to systems of linear equations but it also helps determine whether the system has a unique solution and helps establish relationship and properties of matrices. Undoubtedly, the computation of such single number called the determinant is fundamental in linear algebra. It is one of the basic concepts in linear algebra which has major applications in various branches of engineering and applied science problems such as in the solutions systems of linear equations and also in finding the inverse of an invertible matrix. Also, many complicated expressions of electrical and mechanical systems can be conveniently handled by expressing them in “determinant form.” Therefore, it has become a mathematical area of increasing significance as the computation of the determinant of an matrix of numbers or polynomials is a classical problem and challenge for both numerical and symbolic methods. Consequently, various direct and nondirect methods such as butterfly method, Sarrus’ rule, triangle’s rule, Gaussian elimination procedure, permutation expansion or expansion by the elements of whatever row or column, pivotal or Chio’s condensation method, Dodgson’s condensation method, LU decomposition method, QR decomposition method, Cholesky decomposition method, Hajrizaj’s method, and Salihu and Gjonbalaj’s method [1–35] have been proposed for finding the determinant of matrices. In the gamut of the methods or rules for finding the determinant of the matrices, Sarrus’ rule (a method of finding the determinant of 3 × 3 matrices named after a French mathematician, Pierre Frédéric Sarrus (1798–1861)) has been shown to be the simplest, easiest, fastest, and very straightforward method. Although the wide range of applications of the rule for the computation of the determinant of 3 × 3 matrices is well established, it is grossly limited in applications since it cannot be used for finding the determinants of 4 × 4 matrices and larger matrices. Moreover, the combined idea of finding determinant of 2 × 2 matrices using butterfly method which is the conventional idea in all literatures and using Sarrus’ rule for finding the determinant of 3 × 3 matrices is termed basketweave method. However, the basketweave method does not work on matrices larger than 3 × 3 [1]. Therefore, for larger matrices, the computations of determinants are carried out by methods such as row reduction or column reduction, Laplace expansion method, Dodgson’s condensation method, Chio’s condensations, triangle’s rule, Gaussian elimination procedure, LU decomposition, QR decomposition, and Cholesky decomposition. However, these methods are not as simple, easy, fast, and very straightforward as basketweave method/Sarrus’ rule. Additionally, the cost of the computation of the determinant of a matrix of order is about arithmetic operations using Gauss elimination; if the order of the matrix is large enough, then the computation is not feasible. Therefore, Rezaifar and Rezaee [1] developed a recursion technique to evaluate the determinant of a matrix. In their quests for establishing a new scheme for the generalization of Rezaifar and Rezaee’s procedure, Dutta and Pal [36] pointed out the limitation of Rezaifar and Rezaee’s procedure as it fails to evaluate the values of the determinants of matrices in some cases. Therefore, in this paper, a new method using different schemes based on Sarrus’ rule was developed to carry out the computation of the determinant of 4 × 4 matrices. The developed method is shown to be very quick, easy, efficient, very usable, and highly accurate. It creates opportunities to find other new methods based on Sarrus’ rule to compute determinants of higher orders. Also, the new approach has been shown to be applicable to the computation of determinants of larger matrices such as 5 × 5, 6 × 6, and all other () matrices.

#### 2. Definition of Determinants

The determinant of matrix square matrix is a real number or a function of the elements of the matrix which well determines something about the matrix. It determines whether the system has a unique solution and whether the matrix is singular or not.

The determinant of an -order matrix will be called sum, which has ! different terms which will be formed of matrix elements.

Let be an matrix:Then determinant of iswhere The determinant of matrix could also be written in Laplace cofactor form as

#### 3. Existing Methods of Computation of Determinants

The easiest way to find the determinant of a matrix is to use a computer program which has been optimized so as to reduce the computational time and cost, but there are several ways to do it by hand [37–43]. Therefore, the computation of determinants of matrices has been carried out by some existing methods in literature such as basketweave method, butterfly method, Sarrus’ method, triangle’s rule, Gaussian elimination procedure, permutation expansion or Laplace expansion by the elements of whatever row or column, row reduction method, column reduction method, pivotal or Chio’s condensation method, Dodgson’s condensation method, LU decomposition method, QR decomposition method, Cholesky decomposition method, Hajrizaj’s method, Salihu and Gjonbalaj’s method, Rezaifar and Rezaee’s method, and Dutta and Pal’s method. The simplest among these methods is the basketweave method which could be stated as the combination of butterfly method for determinant computation of 2 × 2 matrices and Sarrus’ rule for determinant computation of 3 × 3 matrices.

##### 3.1. The Butterfly Method

A 2 × 2 matrix is written asIn order to find the determinant of the 2 × 2 matrix, we carry out the diagonal products. We then subtract the diagonal product as we go right to left from the diagonal product of a square matrix as left to right as follows:

*Example 1. *Evaluate :

##### 3.2. Sarrus’ Method

A 3 × 3 matrix is written asSarrus’ rule which is sometimes also called the basketweave method is an alternative way to evaluate the determinant of a 3 × 3 matrix. It is a method that is only applicable to 3 × 3 matrices. It follows the same process as carried out in the 3 × 3 matrix, except that we need to repeat the first two columns to the right of the original matrix and then do the basketweave method. Therefore, a 3 × 5 array is constructed by writing down the entries of the 3 × 3 matrix and then repeating the first two columns at the back of the third column. We calculate the products along the six diagonal lines shown in the diagram. The determinant is equal to the sum of products along diagonals labeled 1, 2, and 3 minus the sum of the products along the diagonals labeled 4, 5, and 6. An example is shown as follows:

*Example 2. *Evaluate :Multiplication of the numbers on the same line, addition of the ones from down-going lines, and subtraction of the ones from up-going lines are an approach that led to the name “the basketweave method.” Unfortunately, the simple weave method does not work on matrices larger than 3 × 3.

The use of Laplace cofactor expansion along either the row or column is a common method for the computation of the determinant of 3 × 3, 4 × 4, and 5 × 5 matrices. The evaluation of the determinant of an matrix using the definition involves the summation of ! terms, with each term being a product of factors. As increases, this computation becomes too cumbersome. This drawback is not only peculiar to Laplace cofactor expansion method as other common methods developed in literatures also required additional computational cost and time for the computation of determinant. Therefore, in recent times, different techniques have been devised in literatures. However, these techniques are not as simple, easy, fast, and very straightforward as the basketweave method/Sarrus’ rule. Additionally, many of them come with relatively high computational cost and time.

#### 4. The Development of the New Methods for the Computation of Determinants

Consider a 4 × 4 matrix whose determinant is required, given as follows:Following the definition given in Section 2, the conventional method of finding the determinant by Laplace cofactor expansion method is carried out as follows.

Expanding along the first row, we haveAgain, expanding each of the 3 × 3 matrices along the first row, we haveNow, we haveSo, it is shown that 4! different terms will be needed to compute the determinant of the forth-order matrices.

In order to generate these 4! different terms (24 terms), we have the following 3 different 4 × 4 matrices as follows: In the arrangements, , , , and represent the first, the second, the third, and the fourth columns, respectively, as given in the original 4 × 4 matrix. We could see in (15) that the first arrangement of matrix remains the same as given in the original matrix . To get the second arrangement of another 4 × 4 matrix, remove and transfer the second column in the first arrangement to the last column of the given 4 × 4 matrix . To get the third arrangement of another new 4 × 4 matrix, remove and transfer the second column in the second arrangement to be the last column of the second 4 × 4 matrix. This forms the third 4 × 4 matrix. After the third step, we need not go further to perform the procedure of removing and transferring the second column in the first arrangement to the last column of the given 4 × 4 matrix because if we do we will end up repeating the first step or getting the first original 4 × 4 matrix in this procedure. In fact, this approach helps us know when to stop the procedures. That is why the last arrangement was cancelled.

Following the procedure, we have 4 × 4 matrices.

The first 4 × 4 matrix isThe second 4 × 4 matrix isThe third 4 × 4 matrix isFrom the above, 10 new schemes based on Sarrus’ rule were developed for the computation of the determinant of the 4 × 4 matrix.

In the new method/scheme, the next step to find after the arrangements is as follows:(1)In the first submatrix , rewrite the 1st, 2nd, and 3rd columns on the right-hand side of matrix (as columns 5, 6, and 7). To the resulting 4 × 7 augmented matrix, assign “+” to the leading element in the odd numbered columns and assign “−” sign to the leading element in the even numbered columns. This gives This is the first part of the solution of the computation of determinant of the given 4 × 4 matrix.(2)In the second submatrix , rewrite the 1st, 2nd, and 3rd columns on the right-hand side of matrix ( as columns 5, 6, and 7). As in the first step, to the augmented matrix, assign “+” to the leading element in the odd numbered columns and assign “−” sign to the leading element in the even numbered columns and then apply Sarrus’ rule. This is the second part of the computation of determinant of the given 4 × 4 matrix.(3)In the third submatrix , rewrite the newest 1st, 2nd, and 3rd columns on the right-hand side of matrix (as columns 5, 6, and 7). And again, to the augmented matrix, assign “+” to the leading element in the odd numbered columns and assign “−” sign to the leading element in the even numbered columns and then apply Sarrus’ rule. This is the third part of the computation of determinant of the given 4 × 4 matrix.(4)For each of augmented matrices , , and , apply Sarrus’ rule by adding the products along the four full diagonals that extend from upper left to lower right and subtract the products along the four full diagonals that extend from the lower left to the upper right. After finding the determinant of the augmented matrices , , and , the addition of the results after applying Sarrus’ rule on the augmented matrices , , and is the determinant of . This is shown as follows:Therefore, we have () which is equivalent in all entireties to (14) when Laplace cofactor expansion method is used:Alternatively, the new scheme could be carried out in another way. In the alternative way, the algorithm still remains the same but the difference is in the manner where the submatrices , , and are constructed. In this scheme, we rewrite the 1st, 2nd, and 3rd columns on the left-hand side of matrix (as columns 0, −1, and −2) to form the required 4 × 7 augmented matrix. Therefore, we have the submatrices , , and given as

As before,