Метод Гаусса, исправить код

var
  i,j,k,n:integer;
  A,A2:array[1..3,1..3] of real;
  x,f,f2,sum:array[1..3] of real;
  Aii,c: real;
begin
  write('Введие порядок матрицы n=');
  readln(n);
  A[1,1]:=1; A[1,2]:=2; A[1,3]:=1; f[1]:=-1;
  A[2,1]:=4; A[2,2]:=-1; A[2,3]:=3; f[2]:=2;
  A[3,1]:=-1; A[3,2]:=4; A[3,3]:=-1; f[3]:=-5;
  writeln(A);
  
  for i:=1 to n do
  begin   
   k:=0;
   for j:=1 to n do
   begin
     if A[i,j]=0 then k:=k+1;
   end;
   if k=n then exit;// Решения нет
  end;
  for j:=1 to n do
  begin   
    k:=0;
    for i:=1 to n do
    begin
      if A[i,j]=0 then k:=k+1;
    end;
     if k=n then exit;// Решения нет
  end;
  //Прямой ход.
  for i:=1 to n do 
  begin
    Aii:=A[i,i];
    for k:=i to n do A[i,k]:=A[i,k]/Aii;
    f[i]:= f[i]/Aii;
    if i<n then 
    begin
      for j:=i+1 to n do 
      begin
        c:=A[j,i];
        for k:=i to n do 
        begin
          A[j,k]:=A[j,k]-c*A[i,k];
        end;
      end;
    end;
  end;
     //Обратный ход
  for i:=n downto 1 do
  begin
    x[i]:=f[i];
    for j:=n downto i+1 do
    begin
      x[i]:=x[i]-A[i,j]*x[j];
    end;
  end;
  writeln(x);
  writeln(A);
  writeln(f);
     //Проверка решения
  for i:=1 to n do
  begin
    sum[i]:=0;
    for j:=1 to n do
    begin
      sum[i]:=sum[i]+A[i,j]*x[j];
    end;
    sum[i]:= sum[i]-f[i];
  end;
  writeln(sum);
     //if f<>x then writeln ('Работает');
end.