Integral definida en delphi

[guiller130292]

Buenos dias amigos, tengo un problema en delphi que al hacer una función (function) para calcular la integral definida de la siguiente función (y=(x^2)*cos(2x)+3) en un intervalo (a,b) el programa se queda pegado, anexo el código para que lo analicen y me puedan ayudar (separé lo más que pude para ver si el error estaba por ahí pero no logré nada).

Código Delphi [-]
unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Math;

type
  TForm1 = class(TForm)
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Button1: TButton;
    Label1: TLabel;
    Label2: TLabel;
  procedure Button1Click(Sender: TObject);

  private
  { Private declarations }
  public
  { Public declarations }
end;

var
  Form1: TForm1;

implementation

{$R *.dfm}
function integral(a,b:extended):extended;
var
  c,d,e,f,g,h,i,j,k,l,m,n,o,x1,x2,y1,y2,at:extended;
begin
  c:=(b-a)/10000;
  x1:=a;
  x2:=x1+c;
  e:=0;
  repeat
    f:=Power(x1,2);
    g:=cos(2*x1);
    h:=f*g;
    y1:=h+3;
    i:=Power(x2,2);
    j:=cos(2*x2);
    k:=i*j;
    y2:=k+3;
    l:=x2-x1;
    m:=y2-y1;
    n:=l*m;
    at:=n/2;
    if y1 < y2 then
      d:=((l)*y1)+at
    else
      d:=((l)*y2)+at;
    e:=e+d;
    x1:=x2;
    x2:=x2+c;
  until x1=b;
  result:=e;
end;

procedure TForm1.Button1Click(Sender: TObject);
var
  a,b,c:extended;
begin
  a:=Strtofloat(Edit1.Text);
  b:=strtofloat(Edit2.Text);
  if b < a then
    Showmessage('El límite inferior debe ser menor al límite superior')
  else
  begin
    c:=Integral(a,b);
    Edit3.Text:=Floattostr(c);
  end;
end;

end.

[[nlsgarcia]]

guiller130292,

¡Bienvenido al Club Delphi!

Te sugiero revisar la Guía de estilo de los foros

¡Gracias por tu cooperación!

Saludos,

Nelson.

[guiller130292]

Amigo Nelson, realmente no sé con qué motivo me has recomendado la guía de estilos de foro, no sé si es para orientarme en las publicaciones o si cometí algún error, de ser así por favor indícamelo más explícitamente que no lo consigo, gracias.

[Casimiro Notevi]

¿Has leído la guía de estilo?
¿En qué quieres que te ayudemos?

[guiller130292]

Mi duda con el código es que realmente no sé dónde está el error. Elaboré una función que según mi criterio debería calcular la integral definida de la función antes especificada en un intervalo introducido por el usuario. El programa compila, pero al pulsar el boton que debe calcular la integral, se queda pegado y no consigo el error en el código y me parece que la lógica está correcta, si por favor alguien puede revisarlo lo agradezco mucho. Gracias por sus respuestas.

[Casimiro Notevi]

¿Pegado?
¿Eso es bloqueado?, seguramente entra en un bucle infinito porque nunca se cumple que:

Código Delphi [-]
...until x1=b;