Изменить заполнение матрицы (Pascal)

[peggass2]

Есть программа которая решает СЛАУ методом Гаусса.
Заполнение матриц идет следующим образом:

То есть отдельно записывается каждый элемент матрицы с новой строки.
Надо что-бы заполняли по строкам пример вот:

Вот процедура заполнения матрицы:

procedure zapolnenie;
  var i,j:integer;
  begin
    write(' Какого порядка СЛАУ: ');
    readln(n);  
    for i:=1 to n do 
      begin
        writeln('Уравнение:',i);
        for j:=1 to n do
          begin
            write('Введите A(',i,',',j,'):');
            readln(a[i,j]);
          end;
        write('Свободный член: B(',i,'):');
        readln(b[i]);
      end;
  end;

Полный код программы (если понадобится):

const max=20;

var i,j,k,n,vubor:integer;
    a:array[1..max,1..max] of real;
    b:array[1..max] of real;
    x:array[1..max] of real;
    h:real;
    
procedure zapolnenie;
  var i,j:integer;
  begin
    write(' Какого порядка СЛАУ: ');
    readln(n);  
    for i:=1 to n do 
      begin
        writeln('Уравнение:',i);
        for j:=1 to n do
          begin
            write('Введите A(',i,',',j,'):');
            readln(a[i,j]);
          end;
        write('Свободный член: B(',i,'):');
        readln(b[i]);
      end;
  end;
  
procedure treugolniy_vid;
  var i,j,k:integer;
  begin
    for i:=1 to n do
      begin
        for j:=i+1 to n do
          begin
            a[j,i]:=-a[j,i]/a[i,i];
            for k:=i+1 to n do
            a[j,k]:=a[j,k]+a[j,i]*a[i,k];
            b[j]:=b[j]+a[j,i]*b[i];
          end;
      end;
  end;
  
procedure obratniy_hod;
  var i,j:integer;
  begin
    x[n]:=b[n]/a[n,n];
    for i:=n-1 downto 1 do
      begin
        h:=b[i];
        for j:=I+1 to n do
        h:=h-x[j]*a[i,j];
        x[i]:=h/a[i,i];
      end;
  end;
  
procedure vuvod;
  var i:integer;
  begin
    writeln('Решение системы:');
    for i:=1 to n do
    writeln('X(',i,')= ',x[i]:3:2);
    readln;
  end;
  
begin
   
  repeat
    writeln('-------------------------------------------------------------------------------');
    write('Использовать систему по умолчанию (1) или ввести вручную (2), для выхода (3):');
    readln(vubor);
    case vubor of
      1: begin
         
         end;
      2: begin 
          zapolnenie;
          treugolniy_vid;
          obratniy_hod;
          vuvod;
        end;
    end;
    until vubor=3;
end.