Delphi метод простых итерация

[world12_tk]

Здраствуйте уважаемые форумчанины.
Мне нужно сделать прогу, которая считает линейные уравнения методом простых итераций... и возникла такая проблема.
На данный момент я пытаюсь привести матрицу 3x3 к рабочиму ввиду.
Ошибка в том, что когда я начинаю делить строки матрицы на диагональный ее элемент, то почему то у меня счетчики i и j идут в обратную сторону от 3 до 1. и из-за этого у меня некоторые элементы не делятся. Ошибка в этой процедуре:

procedure CalculationOfMaximum;
var
  i,j,tempJ:integer;
  temp:extended;
begin
  for j:=1 to count do
      vsego[j]:=vsego[j]/matrix[j,j];

  for i:=1 to count do
      for j:=1 to count do
         matrix[i,j]:=-1*matrix[i,j]/matrix[i,i];
   for j:=1 to count do
      matrix[j,j]:=vsego[j];
   for i:=2 to count do
      for j:=1 to count do begin
       if i=j then begin
        tempJ:=j;
       repeat
         temp:=matrix[i,tempJ-1];
         matrix[i,tempJ-1]:=matrix[i,tempJ];
         matrix[i,tempJ]:=temp;
         dec(tempJ);
       until tempJ<1;
       end;

      end;
end;

Сам проект приложен в архиве....
Помогите мне исправить ее.....
заранее спасибочки