Решение СЛАУ (Системы Линейных Алгебраических Уравнений) методом Гаусса

const n=3;n1=4;
type stroka=array[1..n1] of real;
     matr=array[1..n] of stroka;
     stolb=array[1..n] of real;
 
var  a:matr;
     x:stolb;
     maxabs,v:real;
     i,j,k,l:integer;
 
begin
 
{“Прямой ход”}
{Формирование матрицы А 3 х 4, в 4-ом столбце находятся свободные члены системы}
for i:=1 to n do                                                     
    for j:=1 to n1 do
     begin 
        write('a[',i,'][',j,']=');
        read(a[i][j]);
      end;
for i:=1 to n do
begin 
    maxabs:=abs(a[i][i]);
    k:=i;
    for l:=i+1 to n do
        if abs(a[l][i])>maxabs then 
            begin
              maxabs:=abs(a[l][i]);
              k:=l;
            end;
if k<>i then  
  for j:=i to n1 do 
    begin 
    v:=a[i][j];
    a[i][j]:=a[k][j];
    a[k][j]:=v;  {Выбор рабочей строки матрицы}
    end;
 
{Этап исключения хi из l-го уравнения  с помощью  i-го}
v:=a[i][i];
for j:=i to n1 do a[i][j]:=a[i][j]/v;
    for l:=i+1 to n do begin 
        v:=a[l][i];
        for j:=i+1 to n1 do 
          a[l][j]:=a[l][j]-a[i][j]*v;
    end;
 
{«Обратный ход»}
x[n]:=a[n][n1];
for i:=n-1 downto 1 do 
    begin 
    x[i]:=a[i][n1];
    end;
for j:=i+1 to n do  
    x[i]:=x[i]-a[i][j]*x[j];
end;
 
{Вывод полученного решения системы уравнений}
for i:=1 to n do 
    writeln('x',i,'=',x[i]:5:2);
end.

for j:=n-1 downto 1 do