¿Por qué Integrate se equivoca?

Sep 24 2020

¿Por qué Integratese equivoca?

b = 2;

NIntegrate[E^(-x^2) HypergeometricU[-b/2, 1/2, x^2], {x, -Infinity, Infinity},
           AccuracyGoal -> 4]
(* -7.67473*10^-9 *)

Clear[b];

Integrate[E^(-x^2) HypergeometricU[-b/2, 1/2, x^2], {x, -Infinity, Infinity}]
(* ConditionalExpression[-(((-2 + b) Sqrt[π] Hypergeometric2F1[1, (1 - b)/2, 1/2, 1])/
                           Gamma[1 - b/2]), Re[b] > 0] *)

f[b_] := ConditionalExpression[-(((-2 + b) Sqrt[π]
                                  Hypergeometric2F1[1, (1 - b)/2, 1/2, 1])/
                                 Gamma[1 - b/2]), Re[b] > 0]

f[2]
(* Indeterminate *)

Respuestas

1 Andreas Sep 24 2020 at 01:19

Mathematica no puede determinar el límite f [b] para b-> 2.

Si no se integra a Infinity, pero solo hasta un límite c, obtiene

c = 10; Table[{NIntegrate[Exp[-x^2] HypergeometricU[-b/2, 1/2, x^2], {x, -c, c}],c MeijerG[{{1/2}, {(1 - b)/2}}, {{0, 1/2}, {-(1/2)}}, c^2]}, {b, 0,
3, 0.2}] // N

con buen acuerdo.