Расчёт системы линейных уравнений методом Гаусса

procedure TForm1.SBSolveClick(Sender: TObject);
var
   a: Matrix;
   b,x: Vector;
   n,i,j,k: Integer;
   h: Double;
begin //Ввод данных и размерности системы
  n := StrToIntDef(Text, SG1.ColCount); //количество столбцов
  for i:=1 to n-1 do //Прямой ход
begin
     j:=i;
     while (a[j,i]=0) and (j<=n) do inc(j);
     if j<>i then
     begin
          for k:=1 to n do
          begin
 
          end;
          h:=b[i];
          b[i]:=b[j];
          b[j]:=h;
     end;
     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;
  for i:=n 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;
  LB1.Items.Clear;
  for i:=1 to n do
                  begin
                  with LB1 do
                             begin
                             Height:=n*28;
                             Items.Append('x['+IntToStr(i)+']='+FloatToStr(x[i])); //Вывод результата
                             If frac(x[i])<>0 then
                                                  with LB1 do
                                                             begin
                                                             Height:=28;
                                                             Items.Clear;
                                                             Items.Add('Система несовместна');
                                                             end;
 
                             end;
                  end;
end;