To solve a system of linear algebraic equations (SLAE)
where x and b are the unknown and known vectors of k size, respectively, is a nonsingular matrix, , the entries of matrices and vectors are considered to be real, a simple iterative method reduces (1) to the form
with the matrix C and vector d. The indexed form of (2) is given by
where . Equation (3) is called a linear difference equation (LDE) or an iterative equation.
The mathematical community considers that this method is approximate and it does not give an exact solution of (1) in a finite number of steps:
where I is the identity matrix. This opinion is unreasonable and it is assumed as a postulate. In numerous literature on mathematics, for example in     , there are particular statements given by mathematicians about the fact that a simple iterative method is approximate. The objective of the work is to disprove the established opinion.
2. The Inconsistency of the Postulate for a Homogeneous System
In control theory, the case d = 0 describes a linear discrete system with respect to deviation from a given state. The reduction of the matrix C to the Frobenius form containing a row of coefficients of the characteristic polynomial and the additive feedback in the form of a row vector with the opposite sign make a row of the matrix equal to zero wherein the matrix becomes nilpotent. We denote such a matrix by N and recall its properties
where σ is the spectrum.
The resulting system goes into an equilibrium state and the derivation is eliminated in a finite number of steps. Such a process may be naturally called finite. This property distinguishes it from the process with a theoretically infinite number of steps in a stable system with the spectrum of the matrix inside the unit circle. This means that a simple iterative method for (1) with respect to the error
with С = N, is exact.
3. The Inconsistency of the Postulate for a Non-Homogeneous System
We pose the question about the nature of the process generated by the non-homogeneous Equation (3)
with the nilpotent matrix.
First, we single out such a process in solving the system (1) in which the matrix is reduced to the triangular form. This approach is used in a number of direct methods, for example, in the well-known Gauss elimination method. The triangular matrix is split into the identity matrix and strictly triangular matrix that is nilpotent. Therefore, the back substitution is a formally iterative process with a nilpotent matrix of a special form.
We write the solution of (6) as
where is the common solution of the homogeneous equation, x0 is an initial vector, is the particular solution of the non-homogeneous equation. The vector xn(h) becomes zero at the step n = k by virtue of (5). This proves the necessity of a nilpotent matrix for a finite process.
An unexpected fact is that this condition is also sufficient. Indeed, in the expression
defines the inverse matrix in (4), with C = N, and establishes finally
This result indicates that an exact solution of the LDE (6) is obtained after performing the iteration with the number equal to the dimension, which completely refutes the postulate that a simple iterative method is approximate.
Next, we needed to answer the question: “Is it possible to transform (1) to (6) without zeroing the elements in advance?” An affirmative answer was received in 2001 for a system with k = 2. Let us show this transformation, for which we represent the system (1) in the form (2)
where h1 and h2 are coefficients. We transfer the term h1x1 to the left hand-side, substitute x2 from the second equation, and divide both sides by a factor of x1. As a result, (10) takes the form (3):
where , , , , , .
The elements c11 and c12 contain the coefficients h1 and h2, by which the spectrum can be changed. To make a zero spectrum, the coefficients of the characteristic polynomial
must satisfy the conditions g1 = 0, g2 = 0, i.e.
Substituting the elements of the matrix C from (11) into (12) gives a linear system
The solution (13) substituted into (11) gives a nilpotent matrix and the right-hand side
, , (14)
where , , with which the equation
generates an exact answer (4) for an arbitrary x0 after the second iteration.
The same difference equation
where , , ,
can be obtained from a system of the form (10) with the differences in the second row.
When x0 = 0, the first iteration of (6) yields d, hence the choice of x0 = d gives an exact answer after the first iteration of (15) and (16).
The generalization of the transformation to SLAE of an arbitrary size made it possible to create a theory of finite iterative processes with the derivation of linear difference equations, with different nilpotent matrices, that generate an exact solution after k iterations for any x0, and the choice of x0 = d reduces the number of iterations by one. The theory has been developed on the basis of a direct method for controlling the spectrum of a matrix   .