Нахождение значения функции F(x) на отрезке [a;b] с шагом h

unit Unit1; 
 
{$mode objfpc}{$H+}
 
interface
 
uses
  Classes, SysUtils, FileUtil, Forms, Controls, Graphics, Dialogs, StdCtrls,
  Grids;
 
type
 
  { TForm1 }
 
  TForm1 = class(TForm)
    Button1: TButton;
    Button2: TButton;
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    StringGrid1: TStringGrid;
    procedure Button1Click(Sender: TObject);
  private
    { private declarations }
  public
    { public declarations }
  end; 
 
var
  Form1: TForm1; 
 
implementation
 
{$R *.lfm}
 
{ TForm1 }
 
procedure TForm1.Button1Click(Sender: TObject);
    var
      f:real;
      r:integer;
begin
  r:=1;
  f:=StrToFloat(Edit1.Caption);
  while f<=StrToFloat(Edit2.Caption) do begin
    StringGrid1.RowCount:=r+1;
    StringGrid1.Cells[0,r]:=floatToStr(f);
    StringGrid1.Cells[1,r]:=floatToStr(2.0/(f-3)-sin(f+1));
 
    inc(r);
    f:=f+StrToFloat(Edit3.Caption);
  end;
 
   if f<>StrToFloat(Edit2.Caption) then begin
   f:=StrToFloat(Edit2.Caption);
   StringGrid1.RowCount:=r+1;
   StringGrid1.Cells[0,r]:=floatToStr(f);
   StringGrid1.Cells[1,r]:=floatToStr(2.0/(f-3)-sin(f+1));
   end;
end;
 
end.