Исправьте пожалуйста некоторые недочеты в программе по решению СЛАР методом Гаусса

program project2;
var
     a:array of array of real;
       b,x:array of real;
         h:real;
           i,j,k,n:integer;
begin
 writeln('Vvedute koluchestvo uravnenuy: ');
 readln(n);
setlength(a,n+1);
setlength(b,n+1);
setlength(x,n+1);
  for i:=1 to n do
begin
setlength(a[i],n+1);
writeln('uravnenue nomer ',i);
writeln('Vvedute koeficientu');
  for j:=1 to n do
begin 
write('a[',i,',',j,']=');read(a[i,j]);
   end;
writeln('Vvedite pravuiu chast uravnenuya');
write('b[',i,']=');read(b[i]); 
   end;
writeln('Matruca');
  for i:=1 to n do
begin
write(i);
  for j:=1 to n do
write(a[i,j]);
writeln(b[i]) 
   end;
  for i:=1 to n-1 do
  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;
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;
writeln('Kornu uravnenuya');
  for i:=1 to n do writeln('x(',i,')=',x[i]);
end.

type
  Tmass=array of Real;
  Tmatrix=array of Tmass;
  
procedure Per(k,n:integer;var a:Tmatrix;var b:Tmass);
//перестановка данной строки и строки с наибольшим главным элементом
var j,i:integer;z:Real;
begin
   z:=a[k,k];i:=k;
   for j:=k+1 to n-1 do
   begin
     if abs(a[j,k])>z then
     begin
     z:=a[j,k];i:=j;
     end;
   end;
     if i>k then
     for j:=k to n-1 do
     begin
     z:=a[i,j];a[i,j]:=a[k,j];a[k,j]:=z;//перестановка коэффициентов
     end;
     z:=b[i];b[i]:=b[k];b[k]:=z;//перестановка свободных членов
end;
 //прямой ход Гаусса
Procedure PrChod(n:integer;var a:Tmatrix;var b:Tmass);
var k,j,i:integer;r:real;
begin
  for k:=0 to n-1 do
    begin
     if a[k,k]=0 then Per(k,n,a,b);//перестановка строк
     for j:=k+1 to n-1 do //пересчет коэффициентов
        begin
          r:=a[j,k]/a[k,k];
          for i:=k to n-1 do
            begin
              a[j,i]:=a[j,i]-r*a[k,i];
            end;
          b[j]:=b[j]-r*b[k];
        end;
     end;
end;
  //обратный ход Гаусса
Procedure ObrChod(n:integer;a:Tmatrix;b:Tmass;var x:Tmass);
var k,j:integer;r,g:real;
begin
  for k:=n-1 downto 0 do
     begin
       r:=0;
       for j:=k+1 to n-1 do
          begin
            g:=a[k,j]*x[j];
            r:=r+g;
          end;
       x[k]:=(b[k]-r)/a[k,k];
     end;
end;