Integral definida en delphi

[orodriguezca]

Redondeado lo dicho por ecfisa: los números tipo real, double, extended, etc son imprecisos y para compararlo lo mejor es acudir a las aproximaciones. En el caso de tu algoritmo yo haría lo siguiente, aunque no es exactamente lo más recomendado:

Código Delphi [-]
...until abs(x1-b) <= 0,0001;

Espero te sirva.

[[nlsgarcia]]

guiller130292,

Empezado por guiller130292

...tengo un problema en Delphi que al hacer una función...para calcular la integral definida de y=(x^2)*cos(2x)+3 en un intervalo (a,b) el programa se queda pegado...

Te comento:

1- El problema del código del Msg #1 el cual hace que entre en un loop infinito, esta en la condición de la estructura de control repeat-until :

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

Como fue señalado en el Msg #6, Msg #10 y Msg #11 por variación de decimales en la comparación de números de punto flotante.

2- El código funciona (No entra en un loop infinito) con la siguiente modificación:

Código Delphi [-]
  repeat
  ...
  until CompareValue(x1, b, 0.0000001) = 0;

Revisa esta información:

1- CompareValue

2- Trunc

3- Floating point numbers

Espero sea útil

Nelson.