Вычисление интеграла методом трапеций

[Rikus]

Здравствуйте! Есть программа, которая вычисляет интеграл методом трапеций. Вроде все работает, но когда дана функция arcsin начинает выдавать ошибку. Что не так и как исправить?

unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Math, jpeg, ExtCtrls, TeeProcs, TeEngine, Chart;

type
  TForm1 = class(TForm)
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Edit4: TEdit;
    Button1: TButton;
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Label4: TLabel;
    Chart1: TChart;
    Image1: TImage;
    Label5: TLabel;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation

{$R *.dfm}

var
a,b,e,i:real;

function Fx(x:real):real;
begin
  Fx:=arcsin(sqrt(x/1+x)); //функция для вычисления интеграла
end;

procedure Metod_Trapeciy;
var
x1,x2,x3:real;
c:longint;
begin
  i:=0;
  for c:=1 to round(abs(b-a)/e) do
    begin
      x1:=Fx(a+c*e);
      x2:=Fx(a+c*e+e);
      if x2>x1 then x3:=x1 else x3:=x2;
      i:=i+abs(x2-x1)*e+abs(x3)*e;
    end;
end;


procedure TForm1.Button1Click(Sender: TObject);
begin
  a:=StrToFloat(Edit1.Text);
  b:=StrToFloat(Edit2.Text);
  e:=StrToFloat(Edit3.Text);

  Metod_Trapeciy;

  Edit4.Text:=FloatToStr(i);
end;

end.

[Аватар]

А точно: Fx:=arcsin(sqrt(x/1+x)); ?
Может: Fx:=arcsin(sqrt(x/(1+x))); и x, я надеюсь, не отрицательный

[Rikus]

Да. Спасибо вам, Аватар!