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.