Уточнение элементов обратной матрицы.

uses crt;
const  eps=0.0001;dd=51;maxiter=500;
label inp1,inp2;
type matr=array[1..10,1..10] of real;
var a,b,a0:matr; eps1:double;
    i,j,imx,n,np:byte;
    s0,s1:real;
 
Procedure Init (var m,m1:matr);
 var i,j:byte;
 begin
  for i:=1 to n do
  begin
  for j:=1 to n do
  begin
  m1[i,j]:=0;
   m[i,j]:=0 + random(55);
   end;
   readln;
 m1[i,i]:=1;
 end; end;
 
procedure PrintMatr(m,m1:matr);
var i,j:byte;
  begin
  for i:=1 to n do
    begin
    if (i=1) then write(np:2,':')
             else write('   ');
    for j:=1 to n do
      write(m[i,j]:6:1);
    for j:=1 to n do
      write(m1[i,j]:6:1);
    writeln;
    end;
  inc(np);
  end;
 
procedure MultString(var a,b:matr;i1:integer;r:real);
var j:byte;
  begin
  for j:=1 to n do
    begin
    a[i1,j]:=a[i1,j]*r;
    b[i1,j]:=b[i1,j]*r;
    end;
  end;
 
procedure AddStrings(var a,b:matr;i1,i2:integer;r:real);
{Процедура прибавляет к i1 строке матриц А i2-ю умноженную на r}
var j:byte;
  begin
  for j:=1 to n do
    begin
    a[i1,j]:=a[i1,j]+r*a[i2,j];
    b[i1,j]:=b[i1,j]+r*b[i2,j];
    end;
  end;
 
procedure MultMatr(a,b:matr;var c:matr);
{Перемножение матриц}
var i,j,k:byte;
    s:real;
  begin
  for i:=1 to n do
  for j:=1 to n do
    begin
    s:=0;
    for k:=1 to n do
    s:=s+a[i,k]*b[k,j];
    c[i,j]:=s;
    end;
  end;
 
function sign(r:real):shortint;
  begin
  if (r>=1) then sign:=1 else sign:=-1;
  end;
 
{========================================================================}
{========================================================================}
{Начало основной программы}
var k:byte; c:matr;
 
begin
Clrscr;
Randomize;
textcolor(15);
 
inp1: writeln('Введите размерн.массива от 1 до 5 ');readln(n);
if (n<1) or (n>5) then
   begin
   writeln('Я же сказал от 1 до 5'); goto inp1
   end
 else goto inp2;
 
inp2: writeln ('Введите 1,2,3');
  writeln ('1: заполнение массива с клавиатуры');
  writeln ('2: заполнение массива генератором');
  writeln ('3: выход');
  readln(i);
case i of
 1:for i:=1 to n do
  begin
  for j:=1 to n do
    begin
    b[i,j]:=0;
    write('Введите ',j,'-ый элемент ', i ,'-строки матрицы ');
       readln(a[i,j]);
       end; readln;
  b[i,i]:=1;
  end;
  2:init(a,b);
  3:exit;
  end;
 
  if n=1 then
  begin
  writeln('Обратной матрицы не существует ');
  goto inp1;
  end;
 
  writeln(' Элементы основной матрицы ');
  for i:=1 to n do
    begin
    for j:=1 to n do
      write( ' ', a[i,j]:6:0,' ');
    writeln;
    end;
 
for i:=1 to n do
for j:=1 to n do
  a0[i,j]:=a[i,j];
writeln('Инициализация матриц -основная и единичная :'); np:=1;
PrintMatr(a,b);
 
for i:=1 to n do
  begin
    for j:=i+1 to n do
    AddStrings(a,b,i,j,sign(a[i,i])*sign(a[j,i]));
     if (abs(a[i,i])>eps) then
       begin
       MultString(a,b,i,1/a[i,i]);
       for j:=i+1 to n do
       AddStrings(a,b,j,i,-a[j,i]);
       end
  else
     begin
     writeln('Обратной матрицы не существует.');  goto inp1;
     end
  end;
if (a[n,n]>eps) then
  begin
  for i:=n downto 1 do
    for j:=1 to i-1 do
      begin
      AddStrings(a,b,j,i,-a[j,i]);
      end;
  end;
 
MultMatr(a0,b,a);
 
writeln('начальная матрица, обратная к ней матрица:');
PrintMatr(a0,b);
writeln('_____________________________________________________');
writeln('Уточнение элементов обратной матрицы');
 
 
Multmatr(a0,b,c);
for i:=1 to n do
 begin
for j:=1 to n do
  write(' ',c[i,j]:6:1,' ');
  writeln;
  end;
 
 readln;
 readkey;
 
end.