Объясните алгоритм нахождения обратной матрицы delphi

[darya_lf]

Понимаю, что тема была создана 1000 раз, но облазила пол форума и не нашла того, что ищу. Вернее, сам алгоритм нашла, спасибо этому форуму огромное. Но вот не могу понять его суть. Например, как в методе с минорами: находишь минор и делишь на определитель, или как в методе Гаусса: элементарными операциями приводим (A|E) к (E|A^(-1)). Объясните пожалуйста, логику действий этого алгоритма.

for k:=0 to mo-1 do
          begin
            for i:=0 to mo-1 do
              for j:=0 to mo-1 do
              begin
                if (i=k) and (j=k) then
                  O[i,j]:=1/OIsh[i,j];
                if (i=k) and (j<>k) then
                  O[i,j]:=-OIsh[i,j]/OIsh[k,k];
                if (i<>k) and (j=k) then
                  O[i,j]:=OIsh[i,k]/OIsh[k,k];
                if (i<>k) and (j<>k) then
                  O[i,j]:=OIsh[i,j]-OIsh[k,j]*OIsh[i,k]/OIsh[k,k];
              end;
            for i:=0 to mo-1 do
              for j:=0 to mo-1 do
                OIsh[i,j]:=O[i,j];
          end;