Intégration des fonctions rationnelles

Le calcul d'une intégrale \(\displaystyle{\int_a^bR(t)dt}\)\(R\) est une fonction rationnelle qui n'a aucun pôle dans l'intervalle \([a, b]\) comporte des calculs des types suivants:

  • \(\displaystyle{\int_a^bE(t)dt}\)\(E\) est une fonction polynomiale,

  • \(\displaystyle{\int_a^b\frac{dt}{(t-\alpha)^r},r\in N^*}\),

  • \(\displaystyle{\int_a^b\frac{tdt}{(t^2+pt+q)^s}\textrm{ et }\int_a^b\frac{dt}{(t^2+pt+q)^s},p^2-4q<0}\).

Les deux premiers ne posent aucun problème. Pour résoudre les deux autres, qui présentent au dénominateur un trinôme du second degré sans racine réelle, on commence, suivant une démarche très fréquente dans ce cas à mettre le trinôme sous forme canonique; ainsi

\(\displaystyle{t^2+pt+q=\left(t+\frac{p}{2}\right)^2+q-\frac{p^2}{4}}\)

En posant alors :

\(\displaystyle{q-\frac{p^2}{4}=\lambda^2\textrm{ et }t+\frac{p}{2}=\lambda u}\)

le changement de variable ainsi défini conduit à se ramener à des calculs d'intégrales du type :

\(\displaystyle{\int_a^b\frac{udu}{(u^2+1)^s}\textrm{ et }\int_a^b\frac{du}{(u^2+1)^s},s\in N^*}\)

Or nous avons déjà rencontré des intégrales de ces deux types :

  • le calcul de \(\displaystyle{\int_a^b\frac{udu}{(u^2+1)^s},s\in N^*}\) se fait en remarquant que le numérateur est au facteur \(1/2\) près la dérivée de \((u^2 + 1)\).

  • celui de \(\displaystyle{\int_a^b\frac{du}{(u^2+1)^s},s\in N^*}\) a été fait en application de la méthode d'intégration par parties, par récurrence.

Rassurez vous, on n'a pas l'habitude de demander de calculer : \(\displaystyle{\int_0^1\frac{dt}{(t^2+1)^{25}}}\). Le calcul est très ennuyeux bien avant 25!

Le calcul avec Maple

On utilise la commande \(\textrm{int}\) (regarder l'aide pour voir la différence avec \(\textrm{Int}\)) et on indique qu'on travaille par rapport à la variable \(x\).

Remarquer l'absence de constante d'intégration. On obtient une primitive.

> int(1/(x^2+1)^25,x) ;

\(\displaystyle{\begin{array}{l} \frac{2687300306925}{35184372088832}\frac{x}{(x^2+1)^2}+\frac{8061900920775}{70368744177664}\frac{x}{x^2+1}+\frac{537460061385}{8796093022208}\frac{x}{(x^2+1)^3}\\ +\frac{1612380184155}{30786325577728}\frac{x}{(x^2+1)^4}+\frac{179153353795}{3848290697216}\frac{x}{(x^2+1)^5}+\frac{895766768975}{21165598834688}\frac{x}{(x^2+1)^6}\\ +\frac{206715408225}{5291399708672}\frac{x}{(x^2+1)^7}+\frac{13781027215}{377957122048}\frac{x}{(x^2+1)^8}+\frac{13781027215}{401579442176}\frac{x}{(x^2+1)^9}\\ +\frac{124029244935}{3815004700672}\frac{x}{(x^2+1)^{10}}+\frac{206715408225}{6676258226176}\frac{x}{(x^2+1)^{11}}+\frac{206715408225}{6979724509184}\frac{x}{(x^2+1)^{12}}\\ +\frac{24805848987}{872465563648}\frac{x}{(x^2+1)^{13}}+\frac{35830670759}{130868345472}\frac{x}{(x^2+1)^{14}}+\frac{1235540371}{46739226624}\frac{x}{(x^2+1)^{15}}\\ +\frac{199280705}{7789871104}\frac{x}{(x^2+1)^{16}}+\frac{199280705}{8033304576}\frac{x}{(x^2+1)^{17}}+\frac{39856141}{1653915648}\frac{x}{(x^2+1)^{18}}\\ +\frac{3231579}{137826304}\frac{x}{(x^2+1)^{19}}+\frac{82861}{3627008}\frac{x}{(x^2+1)^{20}}+\frac{10105}{453376}\frac{x}{(x^2+1)^{21}}+\frac{705}{32384}\frac{x}{(x^2+1)^{22}}\\ +\frac{47}{2208}\frac{x}{(x^2+1)^{23}}+\frac{1}{48}\frac{x}{(x^2+1)^{24}}+\frac{8061900920775}{70368744177664}\arctan(x) \end{array}}\)

Si on veut une intégrale définie, on précise les bornes.

Exemple (la première ligne affiche le résultat sans calculer, la deuxième calcule)

> Int(1/(x^2+1)^25,x=0..1);

\(\displaystyle{\int_0^1\frac{1}{(x^2+1)^{25}}dx}\)

> int(1/(x^2+1)^25,x=0..1);

\(\displaystyle{\frac{8061900920775}{281474976710656}\pi+\frac{215888008265}{2399280300032}}\)