- 3
*x*_{1}+ 2*x*_{2}-*x*_{3}= 1 - 2
*x*_{1}- 2*x*_{2}+ 4*x*_{3}= -2 - -
*x*_{1}+ 1/2*x*_{2}-*x*_{3}= 0

Systems of linear equations belong to the oldest problems in mathematics and they have many applications, such as in digital signal processing, estimation, forecasting and generally in linear programming and in the approximation of non-linear problems in numerical analysis. An efficient way to solve systems of linear equations is given by the Gauss-Jordan elimination or by the Cholesky decomposition.

In general, a system with *m* linear equations and *n* unknowns can be written as

*a*_{11}*x*_{1}+*a*_{12}*x*_{2}+ … +*a*_{1n}*x*_{n}=*b*_{1}*a*_{21}*x*_{1}+*a*_{22}*x*_{2}+ … +*a*_{2n}*x*_{n}=*b*_{2}- :
- :
*a*_{m1}*x*_{1}+*a*_{m2}*x*_{2}+ … +*a*_{mn}*x*_{n}=*b*_{m}

*A***x**=**b**,

If the field is infinite (as in the case of the real or complex numbers), then only the following three cases are possible for any given system of linear equations:

- the system has no solution
- the system has a single solution
- the system has infinitely many solutions

*A***x**=**0**

Especially in view of the abiove applications, several more efficient alternatives to Gauss-Jordan elimination have been developed for a wide diversity of special cases. Many of these improved algorithms are of complexity O(*n*²). Some of the most common special cases are:

- For problems of the form
*A***x**=**b**, where*A*is a symmetric Toeplitz matrix, we can use Levinson recursion or one of its derivatives. One special commonly used Levinson-like derivative is Schur recursion, which is used in many digital signal processing applications. - For problems of the form
*A***x**=*b*, where*A*is a singular matrix or nearly singular, the matrix A is decomposed into the product of three matrices in a process called singular-value decomposition. The left and right hand matrices are left and right hand singular vectors. The middle matrix is a diagonal matrix and contains the singular values. The matrix can then be inverted simply by reversing the order of the three components, transposing the singular vector matrices, and taking the reciprocal of the diagonal elements of the middle matrix. If any of the singular values is too close to zero and therefore close to being singular, they are set to zero.