direktes Verfahren zur Berechnung der Inversen A⁻¹ einer nichtsingulären Matrix A ∈ ℝ^n×n.

Dazu löst man die Gleichungssysteme Ax_k = e_k mit der rechten Seite \({e}_{1}^{T}\space =\space (1,\space 0, \ldots,\space 0)\), \({e}_{2}^{T}\space =\space (0,1,\space 0, \ldots,\space 0),\space \ldots \), \({e}_{n}^{T}\space =\space (0,\space \ldots,\space 0,\space 1)\space \) mit dem Gaußschen Algorithmus und setzt die Inverse A⁻¹ aus den Lösungsvektoren x_k, k = 1,…, n als Spaltenvektoren zusammen.

Typischerweise führt man das Lösen der n Gleichungssysteme in einem Schritt durch, indem man statt der Matrizen (A e_k) hier die Matrix (A I) betrachtet, wobei I = (e₁, e₂,…,e_n) die Einheitsmatrix sei. Nun geht man zunächst wie beim Gauß-Verfahren vor und transformiert A in eine obere Dreiecksmatrix R. Dabei wendet man alle Transformationen auch auf I an. Man erhält (R X).

Anschließend eliminiert man nun noch die oberhalb der Diagonale stehenden Elemente von R durch Addition geeigneter Vielfacher der darunterstehenden Zeilen (spaltenweise von rechts nach links). Dadurch erreicht man Diagonalgestalt \begin{eqnarray}D=\left(\begin{array}{llll}{d}_{11} & & & \\ & {d}_{22} & & \\ & & \ddots & \\ & & & {d}_{nn}\end{array}\right).\end{eqnarray}

Wieder wendet man alle Transformationen auch auf X an und erhält so (D A⁻¹). Multipliziert man nun mit der Inversen von D (d. h. dividiert man die j-te Zeile durch d_jj), so erhält man (I A⁻¹).

Die Berechnung der Inversen einer Matrix ist in der Praxis eine eher seltene Aufgabe. Die Lösung des linearen Gleichungssystems Ax = b ist zwar durch x = A⁻¹b gegeben, wird aber besser mittels des Gauß-Verfahrens berechnet, da dies weniger Rechenzeit beansprucht und numerisch günstiger ist.