Метод обратной матрицы

program obr_matrica;
var a1,a,c:array[1..3,1..6] of real;
       b,x:array[1..3]  of real;
     j,i,n,k:integer;s,s1:real;
begin
  writeln('Введите n');
  read(n);
  writeln('‚Введите матрицу A1');
  for i:=1 to n do
    for j:=1 to n do read(a1[i,j]);
   writeln('‚ Введите матрицу B');
  for i:=1 to n do read (b[i]);
 
  for i:=1 to n do
                  begin
                   for j:=1 to n do a[i,j]:=a1[i,j];
                   for k:=n+1 to 2*n do if k=i+n then a[i,k]:=1
                                                 else a[i,k]:=0;
                  end;
for i:=1 to n do begin
                 s:=a[i,i];
                 for j:=1 to 2*n do a[i,j]:=a[i,j]/s;
                 for k:=1 to n do
                     if k<>i then
                             begin
                             s1:=a[k,i];
                             for j:=1 to 2*n do a[k,j]:=a[k,j]-s1*a[i,j];
                             end;
                 end;
for i:=1 to n do
for j:=1 to n do c[i,j]:=a[i,n+j];
writeln('*********** Обратная матрица С  ‘**********');
for i:=1 to n do
                 begin
                   for j:=1 to n do write(c[i,j]:10:6,'  ');
                   writeln;
                 end;
for i:=1 to n do
                begin
                 x[i]:=0;
                 for j:=1 to n do x[i]:=x[i]+c[i,j]*b[j];
                end;
writeln('*********** Решение *************');
for i:=1 to n do write (x[i]:4:2,'  ');
end.