Аппроксимация временного ряда [Delphi]

[StudeHt]

Здравствуйте уважаемые эксперты! Вообщем нужно аппроксимировать ряд, котрый задается координатами точек из StringGrid'а. Аппроксимирую методом наименьших квадратов. У меня получился следующий код для процедуры аппроксимации:

procedure TForm2.Button2Click(Sender: TObject);
var i : integer;
    Sum_x : integer;
    Sum_y : integer;
    Sum_x_2 : integer;
    Proizv_x_y : integer;
begin
  //загружаем координаты y StringGrid в отдельный массив
  try
    for i := 1 to Count do
      Pixels_y[i] := StrToInt(StringGrid1.Cells[1,i]);
  except on EConvertError do
    begin
      ShowMessage('Ошибка! Координаты точек заданы не верно либо пусты!!!');
      Exit;
    end;
  end;
  //рисуем ряд черным цветом
  for i := 1 to Count do
    Series1.AddXY(i, Pixels_y[i], '', clBlack);
  //аппроксимируем исходный ряд
  for i := 1 to Count do
    Pixels_x[i] := StrToInt(StringGrid1.Cells[0,i]);
  Sum_x := 0; Sum_y := 0;
  for i := 1 to Count do
    begin
      Sum_x := Sum_x + Pixels_x[i];
      Sum_y := Sum_y + Pixels_y[i];
    end;
  Proizv_x_y := 0;
  for i := 1 to Count do
    Proizv_x_y := Proizv_x_y + Pixels_x[i] * Pixels_y[i];
  Sum_x_2 := 0;
  for i := 1 to Count do
    begin
      Pixels_x[i] := Pixels_x[i] * i;
      Sum_x_2 := Sum_x_2+Pixels_x[i];
    end;
  a_0 := 0; b_0 := 0;
  b_0 := (Sum_x * Sum_y - Proizv_x_y) / (Sum_x * (Count - 1));
  a_0 := (Proizv_x_y - (Sum_x * b_0)) / Sum_x_2;
  Series2.Clear;
  //рисуем аппроксимацию ряда красным цветом
  for i := 1 to Count do
    Series2.AddXY(i, a_0 * i + b_0, '', clRed);
end;

Но аппроксимируется временной ряд неверно.. Подскажите в чем ошибка, заранее огромное спасибо!