Cramer's rule

Cramer's rule is a theorem in linear algebra, which gives the solution of a system of linear equations in terms of determinants.

Computationally, it is generally inefficient and thus not used in practical applications which may involve many equations. However, it is of theoretical importance in that it gives an explicit expression for the solution of the system.

It is named after Gabriel Cramer (1704 - 1752).

The system of equations is represented in matrix multiplication form as:

<math>Ax = c<math>

where the square matrix <math>A<math> is invertible and the vector <math>x<math> is the column vector of the variables: <math>( x_i )<math>.

The theorem then states that:

<math>x_i = { \det(A_i) \over \det(A)}<math>

where <math>A_i<math> is the matrix formed by replacing the ith column of <math>A<math> by the column vector <math>c<math>.

Example

A good way to use Cramer's rule on a 2×2 matrix is to use this formula:

Given

<math>ax+by = e<math> and

<math>cx+dy = f<math>,

<math>x = \frac { \begin{vmatrix} e & b \\ f & d \end{vmatrix} } { \begin{vmatrix} a & b \\ c & d \end{vmatrix} } = { ed - bf \over ad - bc}<math>

and

<math>y = \frac { \begin{vmatrix} a & e \\ b & f \end{vmatrix} } { \begin{vmatrix} a & b \\ c & d \end{vmatrix} } = { af - ec \over ad - bc}<math>

Applications to differential geometry

Cramer's rule is also extremely useful for solving problems in differential geometry. Consider the two equations <math>F(x, y, u, v) = 0\,<math> and <math>G(x, y, u, v) = 0\,<math>. When u and v are independent variables, we can define <math>x = X(u, v)\,<math> and <math>y = Y(u, v)\,<math>.

Finding an equation for <math>\partial x/\partial u<math> is a trivial application of Cramer's rule.

First, calculate the first derivatives of F, G, x and y.

<math>dF = \frac{\partial F}{\partial x} dx + \frac{\partial F}{\partial y} dy +\frac{\partial F}{\partial u} du +\frac{\partial F}{\partial v} dv = 0<math>

<math>dG = \frac{\partial G}{\partial x} dx + \frac{\partial G}{\partial y} dy +\frac{\partial G}{\partial u} du +\frac{\partial G}{\partial v} dv = 0<math>

<math>dx = \frac{\partial X}{\partial u} du + \frac{\partial X}{\partial v} dv<math>

<math>dy = \frac{\partial Y}{\partial u} du + \frac{\partial Y}{\partial v} dv<math>

Substituting dx, dy into dF and dG, we have:

<math>dF = \left(\frac{\partial F}{\partial x} \frac{\partial x}{\partial u} +\frac{\partial F}{\partial y} \frac{\partial y}{\partial u} +\frac{\partial F}{\partial u} \right) du + \left(\frac{\partial F}{\partial x} \frac{\partial x}{\partial v} +\frac{\partial F}{\partial y} \frac{\partial y}{\partial v} +\frac{\partial F}{\partial v} \right) dv = 0<math>

<math>dG = \left(\frac{\partial G}{\partial x} \frac{\partial x}{\partial u} +\frac{\partial G}{\partial y} \frac{\partial y}{\partial u} +\frac{\partial G}{\partial u} \right) du + \left(\frac{\partial G}{\partial x} \frac{\partial x}{\partial v} +\frac{\partial G}{\partial y} \frac{\partial y}{\partial v} +\frac{\partial G}{\partial v} \right) dv = 0<math>

Since u, v are both independent, the coefficients of du, dv must be zero. So we can write out equations for the coefficients:

<math>\frac{\partial F}{\partial x} \frac{\partial x}{\partial u} +\frac{\partial F}{\partial y} \frac{\partial y}{\partial u} = -\frac{\partial F}{\partial u}<math>

<math>\frac{\partial G}{\partial x} \frac{\partial x}{\partial u} +\frac{\partial G}{\partial y} \frac{\partial y}{\partial u} = -\frac{\partial G}{\partial u}<math>

<math>\frac{\partial F}{\partial x} \frac{\partial x}{\partial v} +\frac{\partial F}{\partial y} \frac{\partial y}{\partial v} = -\frac{\partial F}{\partial v}<math>

<math>\frac{\partial G}{\partial x} \frac{\partial x}{\partial v} +\frac{\partial G}{\partial y} \frac{\partial y}{\partial v} = -\frac{\partial G}{\partial v}<math>

Now, by Cramer's rule, we see that:

<math>

\frac{\partial x}{\partial u} = \frac{\begin{vmatrix} -\frac{\partial F}{\partial u} & \frac{\partial F}{\partial y} \\ -\frac{\partial G}{\partial u} & \frac{\partial G}{\partial y}\end{vmatrix}}{\begin{vmatrix}\frac{\partial F}{\partial x} & \frac{\partial F}{\partial y} \\ \frac{\partial G}{\partial x} & \frac{\partial G}{\partial y}\end{vmatrix}} <math>

This is now a formula in terms of two Jacobians:

<math>\frac{\partial x}{\partial u} = - \frac{\left(\frac{\partial\left(F, G\right)}{\partial\left(y, u\right)}\right)}{\left(\frac{\partial\left(F, G\right)}{\partial\left(x, y\right)}\right)}<math>

Similar formulae can be derived for <math>\frac{\partial x}{\partial v}<math>, <math>\frac{\partial y}{\partial u}<math>, <math>\frac{\partial y}{\partial v}<math>.de:Cramersche Regel fr:Règle de Cramer ko:크래머공식 pl:Wzory Cramera sv:Cramers regel