Интеграл методом прямоугольников

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls;
 
type
  TForm1 = class(TForm)
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Label4: TLabel;
    Edit4: TEdit;
    Button1: TButton;
    Memo1: TMemo;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
  a,b,n:integer;
  I,y:real;
implementation
 
{$R *.dfm}
function summ(y:real):real;
 var
 h,x,Q:real;
 begin
  h:=(b-a)/n;
  x:=a;
  y:=0;
  Q:=0;
  while b>=x do
   begin
   y:=sqrt(1-0.5*sqr(sin(x)));
   Q:=Q+y;
   x:=x+h;
   end;
 end;
procedure TForm1.Button1Click(Sender: TObject);
var
 h,x:real;
begin
Memo1.Clear;
a:=strtoint(Edit1.Text);
b:=StrToInt(Edit2.Text);
n:=StrToInt(Edit3.Text);
  {h:=(b-a)/n;
  x:=a;
  y:=0;
  Q:=0;
  while b>=x do
   begin
   y:=sqrt(1-0.5*sqr(sin(x)));
   Memo1.Lines.Add(FloatToStr(x)+'    '+FloatToStr(y));
   Q:=Q+y;
   x:=x+h;
   end;}
I:=summ(y)*((b-a)/n);
Edit4.Text:=FloatToStr(I);
end;
 
end.

function summ(Q:real):real;
 var
 h,x,y:real;
 begin
  h:=(b-a)/n;
  x:=a;
  y:=0;
  Q:=0;
  while b>=x do
   begin
   y:=sqrt(1-0.5*sqr(sin(x)));
   Q:=Q+y;
   x:=x+h;
   end;
 end;

result := y;