Intégrale méthode Romberg

[Le chat]

Bonjour,
Je programme actuellement une table du Khi2. Tout fonctionne parfaitement, sauf quand la borne supérieure de l'intégrale citée ci-dessous est un nombre décimal:
....... integration(0.398942280401*e^(-x/2)*x^(3/2-1),x,0,borne supérieure)

Nota: le 3 de la fraction 3/2 est le degré de liberté saisi lors de l'appel de la fonction.

Si la borne supérieure est un nombre décimal j'ai le message suivant:
"Unable to find numeric integral using Romberg method, returning the last computed line of approximations"

Merci

[parisse]

Cela vient de la singularité de l'intégrand en x=0 qui empêche la convergence de la méthode numérique d'intégration de Romberg. Lorsque les bornes sont exactes, la primitive est calculée et évaluée entre les bornes donc pas de problèmes. Vous pouvez vous y ramener en remplaçant la borne sup par exact(...) et en faisant un evalf(...) sur le résultat.