Реализация метода Гаусса для СЛАУ комплексных переменных (Delphi)!!!

[Vitek13]

Здравствуйте!!! Не могли бы вы , уважаемые форумчане ,подсказать в чем заключается моя ошибка. Учитывая то, что процедуры работы с компл.числами верны!

procedure TForm3.Button1Click(Sender: TObject);
  var a,a1: Matrix;
      b,b1,x,x1: Vector;
      h,h1,Q,W,e: Complex;
      i,i1,j,j1,k,k1,n,n1:integer;
begin
  //Ввод данных
  //Размерность системы
  n := StrToIntDef(Text, StringGrid1.ColCount);
  n1 := StrToIntDef(Text, StringGrid3.ColCount);
  //Коэффициенты
  for j := 0 to n - 1 do
    for i := 0 to n - 1 do begin
      a[i + 1, j + 1].re:= StrToFloatDef(StringGrid1.Cells[j, i], 0);
      a[i + 1, j + 1].im:= StrToFloatDef(StringGrid3.Cells[j, i], 0);
      end;
  //Правая часть уравнения
  for I := 0 to n - 1 do begin
  b[i + 1].re := StrToFloatDef(StringGrid2.Cells[0, i], 0);
  b[i + 1].im := StrToFloatDef(StringGrid4.Cells[0, i], 0);
  end;
  //Прямой ход - исключения переменных 
  for i:=1 to n-1 do
    for j:=i+1 to n do
    begin
      del(a[j,i],a[i,i],a[j,i]);
      a[j,i].re:=-a[j,i].re;
      a[j,i].im:=-a[j,i].im;
      end;
      for k:=i+1 to n do  begin
         Multicomplex(a[j,i],a[i,k],Q);
         AddComplex(a[j,k],Q,a[j,k]);
         Multicomplex(a[j,i],b[i],W);
         AddComplex(b[j],W,b[j]);
    end;
    del(b[n],a[n,n],x[n]);
    //Обратный ход - нахождения корней
    for i:=n-1 downto 1 do
    begin
      h:=b[i];
      for j:=i+1 to n do begin
        Multicomplex(x[j],a[i,j],E);
        Sub(h,E,h);
       del(h,a[i,i],x[i]);
    end;         end;
    for i:=1 to n do ListBox1.Items.Append('x(' + IntToStr(i) + ')=' + floattostr(x[i].re));
    for i:=1 to n do ListBox2.Items.Append('x(' + IntToStr(i) + ')=' + floattostr(x[i].im));