Построение графика на компоненте Image

begin
 Image1.Canvas.Pen.Color:=clRed;
 Image1.Canvas.Pen.Width:=3;
 for n:=1 to Form3.Table.RowCount-1  do
  begin
   X := StrToInt(Form3.Table.Cells[0,n]);
   Y:=StrToFloat(Form3.Table.Cells[1,n]);
   Image1.Canvas.LineTo(Round(X), Round(Y));
  end;
  Form3.BitBtn4.Enabled:=False;
end;

n:integer;
begin
 Image1.Canvas.MoveTo(0,Image1.Height-StrToInt(Form3.Table.Cells[0,1]));
  for n:=1 to Form3.Table.RowCount-1 do
    Image1.Canvas.LineTo(Round(StrToInt(Form3.Table.Cells[0,n])), Round(Image1.Height-StrToFloat(Form3.Table.Cells[1,n])));
 
  end;

procedure TForm1.Button1Click(Sender: TObject);
var n,i,x0,y0:integer;
    max,mx,my:double;
begin
//найдем максимальное по модулю Y
with StringGrid1 do
 begin
  n:=rowcount-1;
  max:=0;
  for i:=1 to n do
  if abs(strtofloat(cells[1,i]))>max then max:=abs(strtofloat(cells[1,i]));
 end;
with Image1 do
 begin
  x0:=30;
  y0:=height div 2;
  mx:=(width-x0-10)/n; //масхтаб по Х
  my:=(y0-20)/max;  //масштаб по Y
  with Canvas do
   begin
    moveto(x0,y0);lineto(width-30,y0); //оси
    moveto(x0,10);lineto(x0,height-10);
    for i:=1 to n do //засечки и подписи по оси Х
     begin
      moveto(x0+round(mx*(i-1)),y0-3);
      lineto(x0+round(mx*(i-1)),y0+3);
      textout(x0+round(mx*(i-1))+5,y0+10,inttostr(i));
     end;
    for i:=1 to 3 do //засечки и подписи по оси Y
     begin
      moveto(x0-3,y0-round(my*i/10));
      lineto(x0+3,y0-round(my*i/10));
      textout(x0-20,y0-round(my*i/10),formatfloat('0.0',i/10));
      moveto(x0-3,y0+round(my*i/10));
      lineto(x0+3,y0+round(my*i/10));
      textout(x0-25,y0+round(my*i/10),formatfloat('0.0',-i/10));
     end;
    textout(width-40,y0-20,'X');
    textout(x0+5,5,'Y');
    pen.Color:=clRed;   //график
    moveto(x0,y0-round(my*strtofloat(StringGrid1.Cells[1,1])));
    for i:=2 to n do
    lineto(x0+round(mx*(i-1)),y0-round(my*strtofloat(StringGrid1.Cells[1,i])));
   end;
 end;
end;

 with Canvas do
   begin
    pen.Color:=clBlack;
    fillrect(cliprect);
    .......................

y:= power((a+b*x1),2.5)/(1+power(cos(a*x1),3));