2017-02-22 2 views
1

J'intègre la fonction f (x) = a (n-1)/(1 + ax)^n de x = 0 à x = infini, où 1 < n < 2. Pour grand a, le gradient près de x = 0 devient grand et gsl_integration commence à perdre de sa précision (même en utilisant l'intégration adaptative). Est-il possible de maintenir la même précision qu'une augmentation? Est-ce que diviser l'intégrale jusqu'à deux parties avec un près de x = 0 aide?Exactitude de l'intégration GSL

Répondre

0

L'intégrale de ces limites doit être indépendante de la valeur de a.

Si vous faites la substitution u = ax, vous trouverez que

enter image description here

Il n'y aura donc pas besoin de se soucier de la perte de précision à très grande échelle a.