Программа, выводящая обратную матрицу с вещественными коэффициентами

uses crt;
const
  //
  nmax = 10;
// перечисление переменных
//
var
  // A - матрица
  A: array[,] of real;
  // i, j, k - цикловые переменные
  // n - введенная пользователем размерность
  i, j, k, n, l, m, b, y, x: integer;
  z, buff: real;
 
begin
  WriteLn('Программа обращения матрицы методом Жордана - Гаусса');
  WriteLn;
  Write('Введите размерность матрицы А ( < 11). n = ');
  // чтение размерности
  ReadLn(n);
  if (n < 2) or (n > 10) then
    WriteLn('Нарушено ограничение на размерность матрицы.')
  else
  begin
    // Создаем расширенную матрицы, дополняя справа единичной матрицей
    A := new real[n, 2 * n];
    WriteLn('Ввод значений элементов матрицы A.');
    for i := 0 to n - 1 do 
    begin
      
      A[i, n + i] := 1;
      for j:= 0 to n - 1 do 
        Begin
        Write('Значение [',i,',',j,'] элемента: ');
        Read(A[i, j]); ReadLn;
        end;
    end;
 
    WriteLn;
    WriteLn('Расширенная матрица');
    for i := 0 to n - 1 do
    begin
      WriteLn; 
      for j := 0 to 2 * n - 1 do
        Write(A[i, j]:10:5);
    end;
    WriteLn;
 
    for i := 0 to n - 1 do
    begin
      for j := 0 to n - 1 do 
      begin
        if j = i then 
          if A[i, j] = 0 then
          begin
            repeat
              y := random(n);
            until (y <> i) and (A[y, j] <> 0);
            for x := 0 to 2 * n - 1 do
            begin
              buff := A[x, j];
              A[x, j] := A[b, j];
              A[b, j] := buff;
            end; 
          
          end; 
      end;
    end;
    
    for i := 0 to n - 1 do
    begin
      for j := 0 to n - 1 do 
      begin
        if j = i then 
          z := A[i, i]
        else 
          z := -A[j, i] / A[i, i];
        if z <> 0 then
        begin
          for k := i to 2 * n - 1 do
          begin
            if j = i then
              A[j, k] := A[j, k] / z   
            else
              A[j, k] := A[j, k] + A[i, k] * z;
          end;
        end;
      end;
    end;
    Writeln;
    WriteLn('Обратная матрица');
    for i := 0 to n - 1 do
    begin
      WriteLn; 
      for j := n to 2 * n - 1 do 
        Write(A[i, j]:10:5);
    end;
    Writeln;
  end;
end.