Апроксимация, не могу понять ошибку

[hunter03]

Здравствуйте! У меня следующая задача:
у меня в компоненте Image загружена картинка с кривой которую нужно апроксимировать Ax+B.
Я считываю в массивы координаты Х и координаты У моей линии
Апроксимация делается по условию таким образом: я ввожу коэффициент апроксимации(например 4), и после апроксимации моя линия представляется в виде 4 линий.
Получается для апроксимации мои массивы нужно разбить на 4 части и в каждой из частей нужно подсчитать Эти формулы:http://piccy.info/view3/2217087/bcfd...7d46cd0421adc/
Пытаюсь сделать так:

procedure TForm2.Button1Click(Sender: TObject);
var n,i,n1,j:Integer;  h,h1,h2:Real;
  a,b:array[1..20] of integer;
  x,y:array of Integer;
  z,z1:array[1..100] of real;
  s1,s2,s3,s4:Integer;
begin

n:=StrToInt(Form2.Edit1.Text) ;


SetLength(x,Form1.Table1.RecordCount+1);
SetLength(y,Form1.Table1.RecordCount+1);

n1:=1;{считываю массив А и В}
  while n1<=Form1.Table1.RecordCount do begin
  x[n1]:=Form1.Table1.fieldbyname('X').AsInteger;
  y[n1]:=Form1.Table1.fieldbyname('Y').AsInteger;
    Form1.ListBox1.Items.Add(IntToStr(x[n1]));
  n1:=n1+1;  Form1.Table1.Next;
  end;

  h:=Form1.Table1.RecordCount;
  
h2:=h/n;
h2:=Round(h2);

  h1:=h2;
{разделяю массив на n частей}
  z[1]:=1;
  //ShowMessage(floatToStr(Z[1]));
  for i:=2 to n do begin
  z[i]:=h1;
  h1:=h1+h2;
  z1[i]:=z[i]-z[i-1];
  //ShowMessage(floatToStr(z[i]));
  end;
  z[n]:=h;
 // ShowMessage(floatToStr(Z[n]));




  for i:=1 to n-1 do begin
  s1:=0;
  s2:=0;
  s3 := 0;
  s4 := 0;
  for j:=StrToInt(FloatToStr(z[i])) to StrToInt(FloatToStr(z[i+1])) do begin
  s1 := s1 + x[j] * y[j];
  s2 := s2 + x[j];
  s3 := s3 + y[j];
  s4 := s4 + x[j] * x[j];
    end;
 A[i] := Trunc((z1[i] * s1 - s2 * s3) / (z1[i] * s4 - s2 * s2));
  b[i] := Trunc((s3 - n * s2) / z1[i]);

  {if a[i]<0 then a[i]:=a[i]*(-1);
  if b[i]<0 then b[i]:=b[i]*(-1);

     ShowMessage('a='+floatToStr(a[i]));
      ShowMessage('b='+floatToStr(b[i]));  }
       end   ;

         with Form1.image2.canvas do begin
      for i:=1 to n-1 do  begin
lineto(Round(a[i]),Round(b[i]));
  end;
  end;
  end;

Но не получается никак...может посоветуете как сделать или тут можно поправить код? Если нужно могу поргу залить