Обращение матрицы

PROGRAM OBR_MATRIX_EE_OPREDELITEL;
Uses Crt;
Const n=4; {Размер матрицы}
Type Matrix=array[1..n] of array[1..n] of extended;
         RashMatrix=array[1..n] of array[1..2*n] of extended;
{чтобы активизировать тип extended: options->compiler->8087/80287(поставить X)}
VAR i,j,k,kk:word; { вспомогательные величины}
    A,E,OB:Matrix;   {величины для обозначения данной матрицы, еденичной и обратной соответственно}
    AE:RashMatrix;   { величина для обозначения расширенной матрицы }
    Del,Opredelitel: extended; {вспомогательная величина,величина для обозначения определителя }
 
PROCEDURE Print_matrix( var A: Matrix ); {Вывод матрицы на экран}
    begin
     For i:=1 to n do
               begin
                  For j:=1 to n do
                   If A[i,j]>=0 then
                   write(' ',(Trunc(A[i,j])+round(frac(A[i,j])*10000)*0.0001):4:4, '   ' )
                                else
                   write((Trunc(A[i,j])+round(frac(A[i,j])*10000)*0.0001):4:4, '   ' );
                   writeln;
               end;
   end;
 
 
PROCEDURE Opredelitel_Obratn_Matrix(OB: Matrix;Opredelitel:extended); {Определитель обратной матрицы}
Begin
         Opredelitel:=OB[1,1];
        k:=1;
  While k<=n-1 do
     Begin
           Del:=OB[k,k];
           For j:=k to n do OB[k,j]:=OB[k,j]/Del;
           For i:=k to n do Begin
                         For j:=1 to n do OB[i,j]:=OB[i,j]-OB[k,j]*OB[i,k];
                              End;
           Opredelitel:=Opredelitel*OB[k+1,k+1];
           k:=k+1;
     End;
     write('определитель обратной матрицы равен  ', (Trunc(Opredelitel)+round(frac(Opredelitel)*10000)*0.0001):4:4)
End;
PROCEDURE Obratn_matrix(A,E,OB: Matrix;AE:RashMatrix);
{Обратная матрица, вывод её на экран, вывод на экран её определителя}
Begin
  {Составляем еденичную матрицу}
      For i:=1 to n do
                      For j:=1 to n do
                                      If  j=i then E[i,j]:=1
                                              else E[i,j]:=0;
  {Составляем расширенную матрицу}
       For i:=1 to n do
                    For j:=1 to 2*n do
                                      If  j<=n then AE[i,j]:=A[i,j]
                                                   else AE[i,j]:=E[i,j-n];
 {Преобразования расширенной матрицы}
        k:=1;
   While k<=n do
      Begin
            Del:=AE[k,k];
          For j:=k to 2*n do AE[k,j]:= AE[k,j]/Del;
          For i:=k+1 to n do
                                  For j:=1 to 2*n do AE[i,j]:=AE[i,j]-AE[k,j]*A[i,k];
          k:=k+1;
      End;
      
 
 
 
kk:=n;
   While kk>0 do
      Begin
          For i:=kk downto 1 do
          For j:=1 to 2*n do
          AE[i,j]:=AE[i,j]-AE[kk,j]*AE[i,kk];
          kk:=kk-1;
      end;
  For i:=1 to n do
  For j:=1 to n do
     OB[i,j]:=AE[i,j+n];
     writeln;
 Print_Matrix(OB);{Для вывода обратной матрицы}
 Opredelitel_Obratn_Matrix(OB,Opredelitel)
End;
{тело программы}
BEGIN
ClrScr;
 For i:=1 to n do
   For j:=1 to n do
                   begin
   writeln ('введите элемент',i,'-ой строки ', j,'-го столбца' );
   readln (A[i,j]);
                   end;
writeln (' данная матрица ');
print_matrix(A);
writeln;
write('обратная матрица');
Obratn_matrix(A,E,OB,AE);
END.