 In linear algebra, the '''determinant''' is a value associated with a square matrix. The determinant provides important information when the matrix is that of the coefficients of a system of linear equations, or when it corresponds to a linear transformation of a vector space: in the former case the system has a unique solution if and only if the determinant is nonzero, in the latter case that same condition means that the transformation has an inverse operation. An intuitive interpretation can be given to the value of the determinant of a square matrix with real entries: the absolute value of the determinant gives the scale factor by which area or volume is multiplied under the associated linear transformation, while its sign indicates whether the transformation preserves orientation. Thus a 2&amp;nbsp;&amp;times;&amp;nbsp;2 matrix with determinant −2, when applied to a region of the plane with finite area, will transform that region into one with twice the area, while reversing its orientation.   Determinants occur throughout mathematics. They appear in calculus as the Jacobian determinant in the substitution rule for integrals of functions of several variables. They are used to define the characteristic polynomial of a matrix that is an essential tool in eigenvalue problems in linear algebra. In some cases they are used just as a compact notation for expressions that would otherwise be unwieldy to write down.   The determinant of a matrix '''A''' is denoted det('''A'''), or without parentheses: det&amp;nbsp;'''A'''. An alternative notation, used for brevity, especially in the case where the matrix entries are written out in full, is to denote the determinant of a matrix by surrounding the matrix entries by vertical bars instead of the usual brackets or parentheses. For instance :The determinant of the matrix  \begin{bmatrix}a&amp;b&amp;c\\d&amp;e&amp;f\\g&amp;h&amp;i\end{bmatrix} is written \begin{vmatrix} a &amp; b &amp; c\\d &amp; e &amp; f\\g &amp; h &amp; i \end{vmatrix} and has the value aei-afh-bdi+bfg+cdh-ceg\,.  Although most often used for matrices whose entries are real or complex numbers, the definition of the determinant only involves addition, subtraction and multiplication, and so it can be defined for square matrices with entries are taken from any commutative ring. Thus for instance the determinant of a matrix with integer coefficients will be an integer, and the matrix has an inverse with integer coefficients if and only if this determinant is 1 or −1 (these being the only invertible elements of the integers). For square matrices with entries in a non-commutative ring, for instance the quaternions, there is no unique definition for the determinant, and no definition that has all the usual properties of determinants over commutative rings.  == Definition ==  The determinant of a square matrix ''A'', one with the same number of rows and columns, is a value that can be obtained by multiplying certain sets of entries of ''A'', and adding and subtracting such products, according to a given rule: it is a polynomial expression of the matrix entries. This expression grows rapidly with the size of the matrix (an ''n''-by-''n'' matrix contributes ''n''! terms), so it will first be given explicitly for the case of 2-by-2 matrices and 3-by-3 matrices, followed by the rule for arbitrary size matrices, which subsumes these two cases.  Assume  is a square matrix with  rows and  columns, so that it can be written as : A = \begin{bmatrix} a_{1,1} &amp; a_{1,2} &amp; \dots &amp; a_{1,n} \\ a_{2,1} &amp; a_{2,2} &amp; \dots &amp; a_{2,n} \\ \vdots &amp; \vdots &amp; \vdots &amp; \vdots \\ a_{n,1} &amp; a_{n,2} &amp; \dots &amp; a_{n,n} \end{bmatrix}.\, The entries can be numbers or expressions (as happens when the determinant is used to define a characteristic polynomial); the definition of the determinant depends only on the fact that they can be added and multiplied together in a commutative manner.  The determinant of  is denoted as , or it can be denoted directly in terms of the matrix entries by writing enclosing bars instead of brackets: :\begin{vmatrix}  a_{1,1} &amp; a_{1,2} &amp; \dots &amp; a_{1,n} \\ a_{2,1} &amp; a_{2,2} &amp; \dots &amp; a_{2,n} \\ \vdots &amp; \vdots &amp; \vdots &amp; \vdots \\ a_{n,1} &amp; a_{n,2} &amp; \dots &amp; a_{n,n} \end{vmatrix}.\,  ===2-by-2 matrices===  The determinant of a 2&amp;times;2 matrix  is defined by :\begin{vmatrix} a &amp; b\\c &amp; d \end{vmatrix}=ad - bc.\   If the matrix entries are real numbers, the matrix  can be used to represent two linear mappings:  one that maps the standard basis vectors to the rows of , and one that maps them to the columns of .  In either case, the images of the basis vectors form a parallelogram that represents the image of the unit square under the mapping.  The parallelogram defined by the rows of the above matrix is the one with vertices at (0,0), (''a'',''b''), (''a'' + ''c'', ''b'' + ''d''), and (''c'',''d''), as shown in the accompanying diagram.  The absolute value of ad - bc is the area of the parallelogram, and thus represents the scale factor by which areas are transformed by .  (The parallelogram formed by the columns of  is in general a different parallelogram, but since the determinant is symmetric with respect to rows and columns, the area will be the same.)  The absolute value of the determinant together with the sign becomes the ''oriented area'' of the parallelogram.  The oriented area is the same as the usual area, except that it is negative when the angle from the first to the second vector defining the parallelogram turns in a clockwise direction (which is opposite to the direction one would get for the identity matrix).  Thus the determinant gives the scaling factor and the orientation induced by the mapping represented by .  When the determinant is equal to one, the linear mapping defined by the matrix represents is equi-areal and orientation-preserving.  ===3-by-3 matrices===  is the absolute value of the determinant of the matrix formed by the rows r1, r2, and r3.]]  The determinant of a 3&amp;times;3 matrix :\begin{vmatrix}a&amp;b&amp;c\\d&amp;e&amp;f\\g&amp;h&amp;i\end{vmatrix} = aei + bfg + cdh - afh - bdi - ceg.\    The rule of Sarrus is a mnemonic for this formula: the sum of the products of three diagonal north-west to south-east lines of matrix elements, minus the sum of the products of three diagonal south-west to north-east lines of elements when the copies of the first two columns of the matrix are written beside it as in the illustration at the right.  This formula does not carry over into higher dimensions.  For example, the determinant of :A = \begin{bmatrix}-2&amp;2&amp;3\\ -1&amp; 1&amp; 3\\ 2 &amp;0 &amp;-1\end{bmatrix} is calculated using this rule: {| border=&quot;0&quot; |- |\det(A)\, |=\, |((-2) \cdot 1 \cdot (-1)) + (2 \cdot 3 \cdot 2 ) + (3 \cdot (-1) \cdot 0) |- | | |-\,(2 \cdot 1 \cdot 3) - (0 \cdot 3 \cdot (-2) ) - ((-1) \cdot (-1) \cdot 2) |- | |=\, | 2 + 12 + 0 - 6 - 0 - 2 = 6\, |}  ===''n''-by-''n'' matrices=== The determinant of a matrix of arbitrary size can be defined by the Leibniz formula or the Laplace formula.  The Leibniz formula for the determinant of an ''n''-by-''n'' matrix '''A''' is  :\det(A) = \sum_{\sigma \in S_n} \sgn(\sigma) \prod_{i=1}^n A_{i,\sigma_i}.\   ﻿Here the sum is computed over all permutations σ of the set  A permutation is a function that reorders this set of integers. The position of the element ''i'' after the reordering σ is denoted σ''i''. For example, for ''n'' = 3, the original sequence 1, 2, 3 might be reordered to ''S'' = [2, 3, 1], with S1 = 2, S2 = 3, S3 = 1.  The set of all such permutations (also known as the symmetric group on ''n'' elements) is denoted ''S''''n''. For each permutation σ, sgn(σ) denotes the signature of σ; it is +1 for even σ and &amp;minus;1 for odd σ. Evenness or oddness can be defined as follows: the permutation is even (odd) if the new sequence can be obtained by an even number (odd, respectively) of switches of numbers. For example, starting from [1, 2, 3] and switching the positions of 2 and 3 yields [1, 3, 2], switching once more yields [3, 1, 2], and finally, after a total of three (an odd number) switches, [3, 2, 1] results. Therefore [3, 2, 1] is an odd permutation. Similarly, the permutation [2, 3, 1] is even ([1, 2, 3] → [2, 1, 3] → [2, 3, 1], with an even number of switches). It is explained in the article on parity of a permutation why a permutation cannot be simultaneously even and odd.  In any of the n! summands, the term  :\prod_{i=1}^n A_{i, \sigma_i}\   is notation for the product of the entries at positions (''i'', σ''i''), where ''i'' ranges from 1 to ''n'':  :A_{1, \sigma_1} \cdot A_{2, \sigma_2} \cdots  A_{n, \sigma_n}.\   For example, the determinant of a 3 by 3 matrix ''A'' (''n'' = 3) is  :\begin{align}  \sum_{\sigma \in S_n} \sgn(\sigma) \prod_{i=1}^n A_{i,\sigma_i}  &amp;=\sgn([1,2,3]) \prod_{i=1}^n A_{i,[1,2,3]_i} + \sgn([1,3,2]) \prod_{i=1}^n A_{i,[1,3,2]_i} + \sgn([2,1,3]) \prod_{i=1}^n A_{i,[2,1,3]_i} \\ &amp;+ \sgn([2,3,1]) \prod_{i=1}^n A_{i,[2,3,1]_i} + \sgn([3,1,2]) \prod_{i=1}^n A_{i,[3,1,2]_i} + \sgn([3,2,1]) \prod_{i=1}^n A_{i,[3,2,1]_i}  \\  &amp;=\prod_{i=1}^n A_{i,[1,2,3]_i} - \prod_{i=1}^n A_{i,[1,3,2]_i} - \prod_{i=1}^n A_{i,[2,1,3]_i} + \prod_{i=1}^n A_{i,[2,3,1]_i} + \prod_{i=1}^n A_{i,[3,1,2]_i} - \prod_{i=1}^n A_{i,[3,2,1]_i}  \\  &amp;=A_{1,1}A_{2,2}A_{3,3}-A_{1,1}A_{2,3}A_{3,2}-A_{1,2}A_{2,1}A_{3,3}+A_{1,2}A_{2,3}A_{3,1}+A_{1,3}A_{2,1}A_{3,2}-A_{1,3}A_{2,2}A_{3,1}  \end{align}  This agrees with the rule of Sarrus given in the previous section.  The formal extension to arbitrary dimensions was made by Tullio Levi-Civita, see (Levi-Civita symbol) using a pseudo-tensor symbol. ==== Levi-Civita symbol ====  The determinant for an ''n''-by-''n'' matrix can be expressed in terms of the totally antisymmetric Levi-Civita symbol as follows: :: \det A = \sum_{i_1,i_2,\ldots,i_n=1}^n \varepsilon_{i_1\cdots i_n}  a_{1,i_1} \cdots a_{n,i_n}  ==Properties characterizing the determinant== The determinant has the following properties:  If ''A'' is a triangular matrix, i.e. ''a''''i'',''j'' = 0 whenever ''i'' &gt; ''j'' or, alternatively, whenever ''i'' \det(A) =  a_{1,1} a_{2,2} \cdots a_{n,n} = \prod_{i=1}^n a_{i,i}\,, the product of the diagonal entries of ''A''. For example, the determinant of the identity matrix  ::\textstyle\mathrm{I}_n = \begin{bmatrix}1&amp;0&amp;\ldots&amp;0\\  0 &amp; 1 &amp; \ldots &amp; 0\\  \vdots &amp; \vdots &amp; \ddots &amp; \vdots\\  0 &amp;0&amp;\ldots&amp;1\end{bmatrix}  is one.  If ''B'' results from ''A'' by interchanging two rows or two columns, then det(''B'') = &amp;minus;det(''A''). The determinant is called ''alternating'' (as a function of the rows or columns of the matrix).  If ''B'' results from ''A'' by multiplying one row or column with a number ''c'', then det(''B'') = ''c'' · det(''A''). As a consequence, multiplying the whole matrix by ''c'' yields :\det(c A) = c^n\det(A).\ \,  If ''B'' results from ''A'' by adding a multiple of one row to another row, or a multiple of one column to another column, then \det(B) = \det(A). \,   These properties can be shown by inspecting the definition via the Leibniz formula. For example, the first property is because, for triangular matrices, the product \prod_{i=1}^n a_{i,\sigma(i)} is zero for any permutation σ different from the identity permutation (the one not changing the order of the numbers 1, 2, ..., ''n''), since then at least one ''a''''i'',&amp;sigma;(''i'') is zero.  These four properties can be used to compute determinants of any matrix, using Gaussian elimination. This is an algorithm that transforms any given matrix to a triangular matrix, only by using the operations in the last three items. Since the effect of these operations on the determinant can be traced, the determinant of the original matrix is known, once Gaussian elimination is performed. For example, the determinant of A = \begin{bmatrix}-2&amp;2&amp;-3\\ -1&amp; 1&amp; 3\\ 2 &amp;0 &amp;-1\end{bmatrix}  can be computed using the following matrices:  B = \begin{bmatrix}-2&amp;2&amp;-3\\ 0 &amp; 0 &amp; 4.5\\ 2 &amp;0 &amp;-1\end{bmatrix},  C = \begin{bmatrix}-2&amp;2&amp;-3\\ 0 &amp; 0 &amp; 4.5\\ 0 &amp; 2 &amp;-4\end{bmatrix},  D = \begin{bmatrix}-2&amp;2&amp;-3\\ 0 &amp; 2 &amp;-4\\ 0 &amp; 0 &amp; 4.5 \end{bmatrix}   Here, ''B'' is obtained from ''A'' by adding &amp;minus;1/2 &amp;times; the first row to the second, so that det(''A'') = det(''B''). ''C'' is obtained from ''B'' by adding the first to the third row, so that det(''C'') = det(''B''). Finally, ''D'' is obtained from ''C'' by exchanging the second and third row, so that det(''D'') = &amp;minus;det(''C''). The determinant of the (upper) triangular matrix ''D'' is the product of its entries on the main diagonal: (&amp;minus;2) · 2 · 4.5 = &amp;minus;18. Therefore det(''A'') = +18.  ==Further properties== In addition to the above-mentioned properties characterizing the determinant, there are a number of further basic properties. For example, a matrix and its transpose have the same determinant:  :\det(A^\mathrm{T}) = \det (A).\   These properties are chiefly important from a theoretical point of view. For instance, the relation of the determinant and eigenvalues is not typically used to numerically compute either one, especially for large matrices, because of efficiency and numerical stability considerations.  In this section all matrices are assumed to be ''n''-by-''n'' matrices.  ===Multiplicativity and matrix groups=== The determinant of a matrix product of square matrices equals the product of their determinants:  :\det(AB) = \det (A) \det (B).\   Thus the determinant is a ''multiplicative map''. This property is a consequence of the characterization given above of the determinant as the unique ''n''-linear alternating function of the columns with value&amp;nbsp;1 on the identity matrix, since the function  that maps  can easily be seen to be ''n''-linear and alternating in the columns of ''M'', and takes the value  at the identity. The formula can be generalized  to (square) products of rectangular matrices, giving the Cauchy-Binet formula, which also provides an independent proof of the multiplicative property.  The determinant det(''A'') of a matrix ''A'' is non-zero if and only if ''A'' is invertible or, yet another equivalent statement, if its rank equals the size of the matrix. If so, the determinant of the inverse matrix is given by :\det (A^{-1}) = \frac 1 {\det (A)}.  In particular, products and inverses of matrices with determinant one still have this property. Thus, the set of such matrices (of fixed size ''n'') form a group known as the special linear group. More generally, the word &quot;special&quot; indicates the subgroup of another matrix group of matrices of determinant one. Examples include the special orthogonal group (which if ''n'' is 2 or 3 consists of all rotation matrices), and the special unitary group.  ===Relation to eigenvalues and trace===   Determinants can be used to find the eigenvalues of the matrix ''A'': they are the solutions of the characteristic equation :\det(A - xI) = 0, \,  where ''I'' is the identity matrix of the same dimension as ''A''. Conversely, det(''A'') is the product of the eigenvalues of ''A'', counted with their algebraic multiplicities. The product of all non-zero eigenvalues is referred to as pseudo-determinant.  An Hermitian matrix is positive definite if all its eigenvalues are positive. Sylvester's criterion asserts that this is equivalent to the determinants of the submatrices :A_k := \begin{bmatrix} a_{1,1} &amp; a_{1,2} &amp; \dots &amp; a_{1,k} \\ a_{2,1} &amp; a_{2,2} &amp; \dots &amp; a_{2,k} \\ \vdots &amp;  \vdots &amp; \ddots &amp; \vdots \\ a_{k,1} &amp; a_{k,2} &amp; \dots &amp; a_{k,k} \end{bmatrix}  being positive, for all ''k'' between 1 and ''n''.  The trace tr(''A'') is by definition the sum of the diagonal entries of ''A'' and also equals the sum of the eigenvalues. Thus, for complex matrices ''A'', :\det(\exp(A)) = \exp(\mathrm{tr}(A))\,  or, for real matrices ''A'', :\mathrm{tr}(A) = \log(\det(\exp(A))). \, Here exp(''A'') denotes the matrix exponential of ''A'', because every eigenvalue &amp;lambda; of ''A'' corresponds to the eigenvalue exp(&amp;lambda;) of exp(''A''). In particular, given any logarithm of ''A'', that is, any matrix ''L'' satisfying :\exp(L) = A\, the determinant of ''A'' is given by :\det(A) = \exp(\mathrm{tr}(L)). \, For example, for ''n'' = 2 and ''n'' = 3, respectively, :\det(A) = (\mathrm{tr}(A)^2 - \mathrm{tr}(A^2))/2, \,  :\det(A) = (\mathrm{tr}(A)^3 - 3 \mathrm{tr}(A) \mathrm{tr}(A^2) + 2 \mathrm{tr}(A^3))/6. \, These formulae are closely related to Newton's identities.  A generalization of the above identities can be obtained from the following Taylor series expansion of the determinant:  : \begin{align}  \det(\mathsf{I} + \mathsf{A}) = \sum_{k=0}^{\infty} \frac{1}{k!} \left( - \sum_{j=1}^{\infty} \frac{(-1)^j}{j}\mathrm{tr}(\mathsf{A}^j) \right) ^k\, ,  \end{align}   where '''I''' is the identity matrix.  ===Laplace's formula and the adjugate matrix===  Laplace's formula expresses the determinant of a matrix in terms of its minors. The minor ''M''''i'',''j'' is defined to be the the determinant of the (''n''&amp;minus;1)&amp;times;(''n''&amp;minus;1)-matrix that results from ''A'' by removing the ''i''-th row and the ''j''-th column. The expression (&amp;minus;1)''i''+''j''''M''''i'',''j'' is known as cofactor. The determinant of ''A'' is given by  :\det(A) = \sum_{j=1}^n (-1)^{i+j} a_{i,j} M_{i,j} = \sum_{i=1}^n (-1)^{i+j} a_{i,j} M_{i,j}  Calculating det(''A'') by means of that formula is referred to as expanding the determinant along a row or column. For the example 3-by-3 matrix A = \begin{bmatrix}-2&amp;2&amp;-3\\ -1&amp; 1&amp; 3\\ 2 &amp;0 &amp;-1\end{bmatrix} , Laplace expansion along the second column (''j'' = 2, the sum runs over ''i'') yields: {| border=&quot;0&quot; |- |\det(A)\, |=\, |(-1)^{1+2}\cdot 2 \cdot \det \begin{bmatrix}-1&amp;3\\ 2 &amp;-1\end{bmatrix} + (-1)^{2+2}\cdot 1 \cdot \det \begin{bmatrix}-2&amp;-3\\ 2&amp;-1\end{bmatrix} + (-1)^{3+2}\cdot 0 \cdot \det \begin{bmatrix}-2&amp;-3\\ -1&amp;3\end{bmatrix}  |- | |=\, |(-2)\cdot((-1)\cdot(-1)-2\cdot3)+1\cdot((-2)\cdot(-1)-2\cdot(-3)) |- | |=\, |(-2)\cdot(-5)+8 = 18.\, |- | |} However, Laplace expansion is efficient for small matrices only.  The adjugate matrix adj(''A'') is the transpose of the matrix consisting of the cofactors, i.e., :(\operatorname{adj}(A))_{i,j} = (-1)^{i+j} M_{j,i}.\,   ===Cramer's rule=== For a matrix equation : Ax = b\,  the solution is given by Cramer's rule: : x_i = \frac{\det(A_i)}{\det(A)} \qquad i = 1, \ldots, n \,  where ''A''''i'' is the matrix formed by replacing the ''i''-th column of ''A'' by the column vector ''b''. This fact is implied by the following identity  :A\, \mathrm{adj}(A) = \mathrm{adj}(A)\, A = \det(A)\, I_n \qquad,   It has recently been shown that Cramer's rule can be implemented in O(''n''3) time, which is comparable to more common methods of solving systems of linear equations, such as LU, QR, or singular value decomposition.  ===Sylvester's determinant theorem===  Sylvester's determinant theorem states that for ''A'', an ''m''-by-''n'' matrix, and ''B'', an ''n''-by-''m'' matrix,  :\det(I_\mathit{m} + AB) = \det (I_\mathit{n} + BA),  where I_m and I_n are the ''m''-by-''m'' and ''n''-by-''n'' identity matrices, respectively.  For the case of column vector ''c'' and row vector ''r'', each with ''m'' components, the formula allows the quick calculation of the determinant of a matrix that differs from the identity matrix by a matrix of rank 1:  :\det(I_\mathit{m} + cr) = 1 + rc.  More generally, for any invertible ''m''-by-''m'' matrix ''X'',Proofs can be found in   :\det(X + AB) = \det(X) \det(I_\mathit{n} + BX^{-1}A), and :\det(X + cr) = \det(X) (1 + rX^{-1}c).  ===Block matrices===  Suppose A, B, C, and D are ''n''&amp;times;''n''-, ''n''&amp;times;''m''-, ''m''&amp;times;''n''-, and ''m''&amp;times;''m''-matrices, respectively. Then  :\det\begin{pmatrix}A&amp; 0\\ C&amp; D\end{pmatrix} = \det\begin{pmatrix}A&amp; B\\ 0&amp; D\end{pmatrix} = \det(A) \det(D) .  This can be seen from the Leibniz formula or by induction on n. When A is invertible, employing the following identity  :\begin{pmatrix}A&amp; B\\ C&amp; D\end{pmatrix} = \begin{pmatrix}A&amp; 0\\ C&amp; I\end{pmatrix} \begin{pmatrix}I&amp; A^{-1} B\\ 0&amp; D - C A^{-1} B\end{pmatrix}  leads to  :\det\begin{pmatrix}A&amp; B\\ C&amp; D\end{pmatrix} = \det(A) \det(D - C A^{-1} B) .  When D is invertible, a similar identity with \det(D) factored out can be derived analogously,These identities were taken  that is,  :\det\begin{pmatrix}A&amp; B\\ C&amp; D\end{pmatrix} = \det(D) \det(A - B D^{-1} C) .  When the blocks are square matrices of the same order further formulas hold. For example, if ''C'' and ''D'' commute (i.e., ''CD'' = ''DC''), then the following formula comparable to the determinant of a 2-by-2 matrix holds:Proofs are given at  :\det\begin{pmatrix}A&amp; B\\ C&amp; D\end{pmatrix} = \det(AD - BC)..  ===Derivative===  By definition, e.g., using the Leibniz formula, the determinant of real (or analogously for complex) square matrices is a polynomial function from '''R'''''n''&amp;times;''n'' to '''R'''. As such it is everywhere differentiable. Its derivative can be expressed using Jacobi's formula:  :\frac{\mathrm{d} \det(A)}{\mathrm{d} \alpha} = \operatorname{tr}\left(\operatorname{adj}(A) \frac{\mathrm{d} A}{\mathrm{d} \alpha}\right).  where adj(A) denotes the adjugate of A. In particular, if A is invertible, we have  :\frac{\mathrm{d} \det(A)}{\mathrm{d} \alpha} = \det(A) \operatorname{tr}\left(A^{-1} \frac{\mathrm{d} A}{\mathrm{d} \alpha}\right).  Expressed in terms of the entries of ''A'', these are  :  \frac{\partial \det(A)}{\partial A_{ij}}= \operatorname{adj}(A)_{ji}= \det(A)(A^{-1})_{ji}.  Yet another equivalent formulation is  :\det(A + \epsilon X) - \det(A) = \operatorname{tr}(\operatorname{adj}(A) X) \epsilon + O(\epsilon^2) = \det(A) \operatorname{tr}(A^{-1} X) \epsilon + O(\epsilon^2),  using big O notation. The special case where A = I, the identity matrix, yields  :\det(I + \epsilon X) = 1 + \operatorname{tr}(X) \epsilon + O(\epsilon^2).  This identity is used in describing the tangent space of certain matrix Lie groups.  If the matrix A is written as A = \begin{bmatrix}\mathbf{a} &amp; \mathbf{b} &amp; \mathbf{c}\end{bmatrix} where '''a''', '''b''', '''c''' are vectors, then the gradient over one of the three vectors may be written as the cross product of the other two:  :  \begin{align}  \nabla_\mathbf{a}\det(A) &amp;= \mathbf{b} \times \mathbf{c} \\  \nabla_\mathbf{b}\det(A) &amp;= \mathbf{c} \times \mathbf{a} \\  \nabla_\mathbf{c}\det(A) &amp;= \mathbf{a} \times \mathbf{b}.  \end{align}   ==Abstract algebraic aspects== ===Determinant of an endomorphism=== The above identities concerning the determinant of a products and inverses of matrices imply that similar matrices have the same determinant: two matrices ''A'' and ''B'' are similar, if there exists an invertible matrix ''X'' such that ''A'' = ''X''&amp;minus;1''BX''. Indeed, repeatedly applying the above identities yields  :\det(A) = \det(X)^{-1} \det(BX) = \det(X)^{-1} \det(B)\det(X) = \det(B) \det(X)^{-1} \det(X) = \det(B).\   The determinant is therefore also called a similarity invariant. The determinant of a linear transformation :T : V \rightarrow V\, for some finite dimensional vector space ''V'' is defined to be the determinant of the matrix describing it, with respect to an arbitrary choice of basis in ''V''. By the similarity invariant, this determinant is independent of the choice of the basis for ''V'' and therefore only depends on the endomorphism ''T''.  === Exterior algebra === The determinant can also be characterized as the unique function :D: M_n(K) \to K\,  from the set of all ''n''-by-''n'' matrices with entries in a field ''K'' to this field satisfying the following three properties: first, ''D'' is an ''n''-linear function: considering all but one column of ''A'' fixed, the determinant is linear in the remaining column, that is :D (v_1, \dots, v_{i-1}, a v_i + b w, v_{i+1}, \dots, v_n) = a D (v_1, \dots, v_{i-1}, v_i, v_{i+1}, \dots, v_n) + b D (v_1, \dots, v_{i-1}, w, v_{i+1}, \dots, v_n)\, for any column vectors ''v''1, ..., ''v''''n'', and ''w'' and any scalars (elements of ''K'') a and b. Second, ''D'' is an alternating function: for any matrix ''A'' with two identical columns  0}}. Finally, ''D''(''I''''n'') = 1. Here ''I''''n'' is the identity matrix.  This fact also implies that any every other ''n''-linear alternating function  satisfies :F(M)=F(I)D(M).\  The last part in fact follows from the preceding statement: one easily sees that if ''F'' is nonzero it satisfies , and function that associates  to  satisfies all conditions of the theorem. The importance of stating this part is mainly that it remains validRoger Godement, ''Cours d'Algèbre'', seconde édition, Hermann (1966), §23, Théorème 5, p.&amp;nbsp;303 if ''K'' is any commutative ring rather than a field, in which case the given argument does not apply.  The determinant of a linear transformation ''A'' : ''V'' &amp;rarr; ''V'' of an ''n''-dimensional vector space ''V'' can be formulated in a coordinate-free manner by considering the ''n''-th exterior power Λ''n''''V'' of ''V''. ''A'' induces a linear map :\Lambda^n A: \Lambda^n V \rightarrow \Lambda^n V, v_1 \wedge v_2 \wedge \dots v_n \mapsto A v_1 \wedge A v_2 \wedge \dots A v_n. As Λ''n''''V'' is one-dimensional, the map &amp;Lambda;''n''A is given by multiplying with some scalar. This scalar coincides with the determinant of ''A'', that is to say :(\Lambda^n A)(v_1 \wedge \dots v_n) = \det(A) \cdot v_1 \wedge \dots v_n. This definition agrees with the more concrete coordinate-dependent definition. This follows from the above characterization of the determinant given above. For example, switching two columns changes the parity of the determinant; likewise, permuting the vectors in the exterior product ''v''1 ∧ ''v''2  ∧ ...  ∧ ''v''''n'' to ''v''2 ∧ ''v''1  ∧ ''v''3 ∧ ...  ∧ ''v''''n'', say, also alters the parity.  For this reason, the highest non-zero exterior power &amp;Lambda;''n''(''V'') is sometimes also called the determinant of ''V'' and similarly for more involved objects such as vector bundles or chain complexes of vector spaces. Minors of a matrix can also be cast in this setting, by considering lower alternating forms Λ''k''''V'' with ''k'' r \cdot s = s \cdot r is supposed to hold for all elements ''r'' and ''s'' of the ring. For example, the integers form a commutative ring.  Many of the above statements and notions carry over mutatis mutandis to determinants of these more general matrices: the determinant is multiplicative in this more general situation, and Cramer's rule also holds. A square matrix over a commutative ring ''R'' is invertible if and only if its determinant is a unit in ''R'', that is, an element having a (multiplicative) inverse. (If ''R'' is a field, this latter condition is equivalent to the determinant being nonzero, thus giving back the above characterization.) For example, a matrix ''A'' with entries in '''Z''', the integers, is invertible (in the sense that the inverse matrix has again integer entries) if the determinant is +1 or &amp;minus;1. Such a matrix is called unimodular.  The determinant defines a mapping between :\mathrm{GL}_n(R) \rightarrow R^\times, \, the group of invertible ''n''×''n'' matrices with entries in ''R'' and the multiplicative group of units in . Since it respects the multiplication in both groups, this map is a group homomorphism. Secondly, given a ring homomorphism , there is a map  given by replacing all entries in  by their images under . The determinant respects these maps, i.e., given a matrix  (''a''''i'',''j'')}}  with entries in , the identity :f(\det((a_{i,j}))) = \det ((f(a_{i,j})))\, holds. For example, the determinant of the complex conjugate of a complex matrix (which is also the determinant of its conjugate transpose) is the complex conjugate of its determinant, and for integer matrices: the reduction modulo&amp;nbsp; of the determinant of such a matrix is equal to the determinant of the matrix reduced modulo&amp;nbsp; (the latter determinant being computed using modular arithmetic). In the more high-brow parlance of category theory, the determinant is a natural transformation between the two functors  and . Adding yet another layer of abstraction, this is captured by saying that the determinant is a morphism of algebraic groups, from the general linear group to the multiplicative group, :\det: \mathrm{GL}_n \rightarrow \mathbb G_m.\,  ==Generalizations and related notions== === Infinite matrices === For matrices with an infinite number of rows and columns, the above definitions of the determinant do not carry over directly. For example, in Leibniz' formula, an infinite sum (all of whose terms are infinite products) would have to be calculated. Functional analysis provides different extensions of the determinant for such infinite-dimensional situations, which however only work for particular kinds of operators.  The Fredholm determinant defines the determinant for operators known as trace class operators by an appropriate generalization of the formula :\det(I+A) = \exp(\mathrm{tr}(\log(I+A))). \,  Another infinite-dimensional notion of determinant is the functional determinant.  === Notions of determinant over non-commutative rings ===  For square matrices with entries in a non-commutative ring, there are various difficulties in defining determinants in a manner analogous to that for commutative rings. A meaning can be given to the Leibniz formula provided the order for the product is specified, and similarly for other ways to define the determinant, but non-commutativity then leads to the loss of many fundamental properties of the determinant, for instance the multiplicative property or the fact that the determinant is unchanged under transposition of the matrix. Over non-commutative rings, there is no reasonable notion of a multilinear form (if a bilinear form exists with a regular element of ''R'' as value on some pair of arguments, it can be used to show that all elements of ''R'' commute). Nevertheless various notions of non-commutative determinant have been formulated, which preserve some of the properties of determinants, notably quasideterminants and the Dieudonné determinant.  === Further variants === Determinants of matrices in superrings (that is, '''Z'''/2-graded rings) are known as Berezinians or superdeterminants.[  The permanent of a matrix is defined as the determinant, except that the factors sgn(&amp;sigma;) occurring in Leibniz' rule are omitted. The immanant generalizes both by introducing a character of the symmetric group S''n'' in Leibniz' rule.  ==Calculation== Naive methods of implementing an algorithm to compute the determinant include using Leibniz' formula or Laplace's formula. Both these approaches are extremely inefficient for large matrices, though, since the number of required operations grows very quickly: it is of order ''n''! (''n'' factorial) for an ''n''&amp;times;''n'' matrix ''M''. For example, Leibniz' formula requires to calculate ''n''! products. Therefore, more involved techniques have been developed for calculating determinants.  ===Decomposition methods=== Given a matrix ''A'', some methods compute its determinant by writing ''A'' as a product of matrices whose determinants can be more easily computed. Such techniques are referred to as decomposition methods. Examples include the LU decomposition, Cholesky decomposition or the QR decomposition. These methods are of order O(''n''3), which is a significant improvement over O(''n''!)  The LU decomposition expresses ''A'' in terms of a lower triangular matrix ''L'', an upper triangular matrix ''U'' and a permutation matrix ''P'': : A = PLU \, The determinants of ''L'' and ''U'' can be quickly calculated, since they are the products of the respective diagonal entries.  The determinant of ''P'' is just the sign \varepsilon of the corresponding permutation. The determinant of ''A'' is then  : \det(A) = \varepsilon \det(L)\cdot\det(U), \,   Moreover, the decomposition can be chosen such that ''L'' is a unitriangular matrix and therefore has determinant&amp;nbsp;1, in which case the formula further simplifies to  : \det(A) = \varepsilon\det(U).  ===Further methods=== If the determinant of ''A'' and the inverse of ''A'' have already been computed, the matrix determinant lemma allows to quickly calculate the determinant of , where ''u'' and ''v'' are column vectors.  Since the definition of the determinant does not need divisions, a question arises: do fast algorithms exist that do not need divisions? This is especially interesting for matrices over rings. Indeed algorithms with run-time proportional to ''n''4 exist. An algorithm of Mahajan and Vinay, and Berkowitz is based on closed ordered walks (short ''clow''). It computes more products than the determinant definition requires, but some of these products cancel and the sum of these products can be computed more efficiently. The final algorithm looks very much like an iterated product of triangular matrices.  If two matrices of order ''n'' can be multiplied in time ''M''(''n''), where ''M''(''n'')≥''n''''a'' for some ''a''&gt;2, then the determinant can be computed in time O(''M''(''n'')).J.R. Bunch and J.E. Hopcroft, Triangular factorization and inversion by fast matrix multiplication, ''Mathematics of Computation'', 28 (1974) 231–236. This means, for example, that an O(''n''2.376) algorithm exists based on the Coppersmith–Winograd algorithm.  Algorithms can also be assessed according to their bit complexity, i.e., how many bits of accuracy are needed to store intermediate values occuring in the computation. For example, the Gaussian elimination (or LU decomposition) methods is of order O(''n''3), but the bit length of intermediate values can become exponentially long. The Bareiss Algorithm, on the other hand, is an exact-division method based on Sylvester's identity is also of order ''n''3, but the bit complexity roughly the bit size of the original entries in the matrix times ''n''.  ==History== Historically, determinants were considered without reference to matrices: originally, a determinant was defined as a property of a system of linear equations. The determinant &quot;determines&quot; whether the system has a unique solution (which occurs precisely if the determinant is non-zero). In this sense, determinants were first used in the Chinese mathematics textbook ''The Nine Chapters on the Mathematical Art'' (九章算術, Chinese scholars, around the 3rd century BC). In Europe, two-by-two determinants were considered by Cardano at the end of the 16th century and larger ones by Leibniz.Eves, H: &quot;An Introduction to the History of Mathematics&quot;, pages 405, 493&amp;ndash;494, Saunders College Publishing, 1990.A Brief History of Linear Algebra and Matrix Theory :  F. [ ''A History of Mathematics'' p. 80]  In Europe, Cramer (1750) added to the theory, treating the subject in relation to sets of equations. The recurrent law was first announced by Bézout (1764).  It was Vandermonde (1771) who first recognized determinants as independent functions.Campbell, H: &quot;Linear Algebra With Applications&quot;, pages 111-112. Appleton Century Crofts, 1971 Laplace (1772) Expansion of determinants in terms of minors: Laplace, Pierre-Simon (de) &quot;Researches sur le calcul intégral et sur le systéme du monde,&quot; ''Histoire de l'Académie Royale des Sciences'' (Paris), seconde partie, pages 267-376 (1772).Muir, Sir Thomas, ''The Theory of Determinants in the historical Order of Development'' [London, England: Macmillan and Co., Ltd., 1906]. gave the general method of expanding a determinant in terms of its complementary minors: Vandermonde had already given a special case. Immediately following, Lagrange (1773) treated determinants of the second and third order. Lagrange was the first to apply determinants to questions of elimination theory; he proved many special cases of general identities.  Gauss (1801) made the next advance. Like Lagrange, he made much use of determinants in the theory of numbers. He introduced the word '''''determinants''''' (Laplace had used ''resultant''), though not in the present signification, but rather as applied to the discriminant of a quantic. Gauss also arrived at the notion of reciprocal (inverse) determinants, and came very near the multiplication theorem.  The next contributor of importance is Binet (1811, 1812), who formally stated the theorem relating to the product of two matrices of ''m'' columns and ''n'' rows, which for the special case of ''m'' = ''n'' reduces to the multiplication theorem. On the same day (November 30, 1812) that Binet presented his paper to the Academy, Cauchy also presented one on the subject. (See Cauchy-Binet formula.) In this he used the word '''''determinant''''' in its present sense,The first use of the word &quot;determinant&quot; in the modern sense appeared in: Cauchy, Augustin-Louis “Memoire sur les fonctions qui ne peuvent obtenir que deux valeurs égales et des signes contraires par suite des transpositions operées entre les variables qu'elles renferment,&quot; which was first read at the Institute de France in Paris on November 30, 1812, and which was subsequently published in the ''Journal de l'Ecole Polytechnique'', Cahier 17, Tome 10, pages 29-112 (1815).Origins of mathematical terms:  summarized and simplified what was then known on the subject, improved the notation, and gave the multiplication theorem with a proof more satisfactory than Binet's.History of matrices and determinants:  With him begins the theory in its generality.  The next important figure was Jacobi (from 1827). He early used the functional determinant which Sylvester later called the Jacobian, and in his memoirs in ''Crelle'' for 1841 he specially treats this subject, as well as the class of alternating functions which Sylvester has called ''alternants''. About the time of Jacobi's last memoirs, Sylvester (1839) and Cayley began their work.The first use of vertical lines to denote a determinant appeared in: Cayley, Arthur &quot;On a theorem in the geometry of position,&quot; ''Cambridge Mathematical Journal'', vol. 2, pages 267-271 (1841).History of matrix notation:   The study of special forms of determinants has been the natural result of the completion of the general theory. Axisymmetric determinants have been studied by Lebesgue, Hesse, and Sylvester; persymmetric determinants by Sylvester and Hankel; circulants by Catalan, Spottiswoode, Glaisher, and Scott; skew determinants and Pfaffians, in connection with the theory of orthogonal transformation, by Cayley; continuants by Sylvester; Wronskians (so called by Muir) by Christoffel and Frobenius; compound determinants by Sylvester, Reiss, and Picquet; Jacobians and Hessians by Sylvester; and symmetric gauche determinants by Trudi. Of the text-books on the subject Spottiswoode's was the first. In America, Hanus (1886), Weld (1893), and Muir/Metzler (1933) published treatises.  Axler in 1995 attacked determinant's place in Linear Algebra.  He saw it as something to be derived from the core principles of Linear Algebra, not to be used to derive the core principles.Down with Determinants:   ==Applications==  ===Linear independence=== As mentioned above, the determinant of a matrix (with real or complex entries, say) is zero if and only if the column vectors of the matrix are linearly dependent. Thus, determinants can be used to characterize linearly dependent vectors. For example, given two vectors ''v''1, ''v''2 in '''R'''3, a third vector ''v''3 lies in the plane spanned by the former two vectors exactly if the determinant of the 3-by-3 matrix consisting of the three vectors is zero. The same idea is also used in the theory of differential equations: given ''n'' functions ''f''1(''x''), ..., ''f''''n''(''x'') (supposed to be ''n''&amp;minus;1 times differentiable), the Wronskian is defined to be : W(f_1, \ldots, f_n) (x)= \begin{vmatrix}  f_1(x) &amp; f_2(x) &amp; \cdots &amp; f_n(x) \\ f_1'(x) &amp; f_2'(x) &amp; \cdots &amp; f_n' (x)\\ \vdots &amp; \vdots &amp; \ddots &amp; \vdots \\ f_1^{(n-1)}(x)&amp; f_2^{(n-1)}(x) &amp; \cdots &amp; f_n^{(n-1)}(x) \end{vmatrix}.  It is zero (for some ''x'') if and only if the given functions and all their derivatives up to order ''n''&amp;minus;1 are linearly independent.  ===Orientation of a basis===  The determinant can be thought of as assigning a number to every sequence of ''n'' in '''R'''''n'', by using the square matrix whose columns are the given vectors. For instance, an orthogonal matrix with entries in '''R'''''n'' represents an orthonormal basis in Euclidean space. The determinant of such a matrix determines whether the orientation of the basis is consistent with or opposite to the orientation of the standard basis. Namely, if the determinant is +1, the basis has the same orientation. If it is −1, the basis has the opposite orientation.  More generally, if the determinant of ''A'' is positive, ''A'' represents an orientation-preserving linear transformation (if ''A'' is an orthogonal 2×2 or 3×3 matrix, this is a rotation), while if it is negative, ''A'' switches the orientation of the basis.  ===Volume and Jacobian determinant=== As pointed out above, the absolute value of the determinant of real vectors is equal to the volume of the parallelepiped spanned by those vectors. As a consequence, if  is the linear map represented by the matrix ''A'', and ''S'' is any measurable subset of '''R'''''n'', then the volume of ''f''(''S'') is given by |det(''A'')| times the volume of ''S''. More generally, if the linear map  is represented by the ''m''-by-''n'' matrix ''A'', then the ''m''-dimensional volume of ''f''(''S'') is given by: :\operatorname {volume} (f(S)) = \sqrt{\det(A^\mathrm{T} A)} \times \operatorname{volume}(S).  By calculating the volume of the tetrahedron bounded by four points, they can be used to identify skew lines. The volume of any tetrahedron, given its vertices '''a''', '''b''', '''c''', and '''d''', is (1/6)·|det('''a'''&amp;nbsp;&amp;minus;&amp;nbsp;'''b''',&amp;nbsp;'''b'''&amp;nbsp;&amp;minus;&amp;nbsp;'''c''', '''c'''&amp;nbsp;&amp;minus;&amp;nbsp;'''d''')|, or any other combination of pairs of vertices that would form a spanning tree over the vertices.  For a general differentiable function, much of the above carries over by considering the Jacobian matrix of ''f''. For :f: \mathbf R^n \rightarrow \mathbf R^n, the Jacobian is the ''n''-by-''n'' matrix whose entries are given by :D(f) = \left (\frac {\partial f_i}{\partial x_j} \right )_{1 \leq i, j \leq n}. \, Its determinant, the Jacobian determinant appears in the higher-dimensional version of integration by substitution: for suitable functions ''f'' and an open subset ''U'' of '''R''''''n'' (the domain of ''f''), the integral over ''f''(''U'') of some other function  is given by : \int_{f(U)} \varphi(\mathbf{v})\, d \mathbf{v} = \int_U \phi(f(\mathbf{u})) \left|\det(\operatorname{D}f)(\mathbf{u})\right| \,d \mathbf{u}. The Jacobian also occurs in the inverse function theorem.  