Вычислить значения функции f(x) на заданном интервале, с шагом h

const xn=-2*pi;
      xk=pi;
      h=0.4;
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
function F(x:real):real;
begin
F:=sqr(sin(x))-sqr(cos(x));
end;
//создание таблицы значений
procedure TForm1.Button1Click(Sender: TObject);
var x:real;
    k:integer;
begin
with Stringgrid1 do
 begin
  colcount:=2;
  rowcount:=2;
  cells[0,0]:='X';
  cells[1,0]:= 'F(x)';
 end;
x:=xn;
k:=1;
while x<xk+h/2 do
 begin
  Stringgrid1.Cells[0,k]:=FloatTostrF(x,ffFixed,4,2);
  Stringgrid1.Cells[1,k]:=FloatTostrF(F(x),ffFixed,6,3);
  x:=x+h;
  if x<=xk then
   begin
    k:=k+1;
    Stringgrid1.RowCount:=Stringgrid1.RowCount+1;
   end;
 end;
end;
//построение графика
procedure TForm1.Button2Click(Sender: TObject);
var xc,yc,i:integer;
    x,mx,my:real;
begin
with Image1 do
 begin
  xc:=Width*2 div 3;
  mx:=(Width-40)/(xk-xn);
  yc:=Height div 2;
  my:=yc-20;
  with Canvas do
   begin
    moveto(0,yc);lineto(Width,yc);
    moveto(xc,0);lineto(xc,Height);
    textout(Width-20,yc-15,'X');
    textout(xc+5,20,'Y');
    textout(xc+5,yc+10,'0');
    for i:=1 to 10 do
     begin
      moveto(xc-3,yc+round(i*my/10));lineto(xc+3,yc+round(i*my/10));
      moveto(xc-3,yc-round(i*my/10));lineto(xc+3,yc-round(i*my/10));
      if i mod 2=0 then
       begin
        textout(xc-30,yc+round(i*my/10)-5,format('%.1f',[-i/10]));
        textout(xc-25,yc-round(i*my/10)-5,format('%.1f',[i/10]));
       end;
      moveto(xc+round(i*mx),yc-3);lineto(xc+round(i*mx),yc+3);
      textout(xc+round(i*mx),yc+10,inttostr(i));
      moveto(xc-round(i*mx),yc-3);lineto(xc-round(i*mx),yc+3);
      textout(xc-round(i*mx),yc+10,inttostr(-i));
     end;
    //график
    x:=xn;
    pen.Color:=clGreen;
    moveto(xc+round(x*mx),yc-round(F(x)*my));
    while x<=xk+h/2 do
     begin
      lineto(xc+round(x*mx),yc-round(F(x)*my));
      x:=x+0.1;
     end;
    //точки
    x:=xn;
    brush.Color:=clRed;
    while x<=xk+h/2 do
     begin
      ellipse(xc+round(x*mx)-3,yc-round(F(x)*my)-3,
            xc+round(x*mx)+3,yc-round(F(x)*my)+3);
      x:=x+h
     end;
   end;
  end;
end;
procedure TForm1.BitBtn1Click(Sender: TObject);
begin
close
end;

rowcount:=2;
  cells[0,0]:='X';
  cells[1,0]:= 'F(x)';
 end;
x:=xn;
k:=1;
while x<xk+h/2 do
 begin
  Stringrid1.Cells[0,k]:=FloatTostrF(x,ffFixed,4,2);