Perché Integrate sbaglia?

Sep 24 2020

Perché Integratesi sbaglia?

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 *)

Risposte

1 Andreas Sep 24 2020 at 01:19

Mathematica non è in grado di determinare il limite f [b] per b-> 2.

Se non esegui l'integrazione in Infinity, ma solo fino a un cutoff c ottieni

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 buon accordo.