Sauter au contenu

Section 3.1 Intégrales doubles

Sous-section Tranches verticales

Supposons que nous voulons calculer la masse d’une plaque d’épaisseur constante et très mince qui remplit la région \(\cR\) dans le plan \(Oxy\text{.}\) Supposons de plus que la densité de la plaque, disons en kilogrammes par mètre carré, dépend de la position. Appelons cette densité \(f(x,y)\text{.}\) Pour simplifier, nous supposerons que \(\cR\) est la région entre la courbe inférieure \(y=B(x)\) et la courbe supérieure \(y=T(x)\text{,}\) avec \(x\) allant de \(a\) à \(b\text{,}\) c’est-à-dire :
\begin{gather*} \cR=\big\{\ (x,y)\ \big|\ a\leqslant x\leqslant b,\ B(x)\leqslant y\leqslant T(x)\ \big\}. \end{gather*}
Nous allons exprimer la masse comme une intégrale double. À titre d’échauffement, rappelons la procédure afin de calculer l’aire de la région \(\cR\text{.}\)
  • Choisissons un nombre naturel \(n\) (que nous ferons tendre vers l’infini).
  • Divisons ensuite \(\cR\) en \(n\) tranches verticales minces, chacune de largeur \(\De x=\frac{b-a}{n}\text{.}\) Posons \(x_i = a + i\,\De x\text{,}\) l’abscisse de l’extrémité droite de la tranche \(i\text{.}\)
  • Pour chaque \(i=1,2,\dots,n\text{,}\) la \(i^{\text{ième}}\) tranche s’étend de \(x= x_{i-1}\) à \(x = x_i\text{.}\) Nous approximons son aire par l’aire d’un rectangle. Nous choisissons pour cela un nombre \(x_i^*\) entre \(x_{i-1}\) et \(x_i\text{,}\) et nous nous approchons de la tranche par un rectangle dont l’extrémité supérieure se trouve à \(y=T(x_i^*)\) et dont l’extrémité inférieure se trouve à \(y=B(x_i^*)\text{.}\) Le \(\text{rectangle}_i^*\) est illustré en bleu dans la figure ci-dessous.
  • Ainsi, l’aire de la \(i^{\text{ième}}\) tranche est approximativement \(\big[T(x_i^*)-B(x_i^*)\big]\De x\text{.}\)
  • Nous obtenons de cette façon une approximation de l’aire totale de la région \(\cR\) par une somme de Riemann :
    \begin{align*} \text{Aire} &\approx \sum_{i=1}^n \big[T(x_i^*)-B(x_i^*)\big]\De x . \end{align*}
  • En évaluant la limite lorsque \(n \to \infty\) (de sorte que la largeur des rectangles tendent vers zéro), nous obtenons une intégrale définie, et la valeur exacte de l’aire de la région \(\cR\) est :
    \begin{gather*} \text{Aire} =\lim_{n\rightarrow\infty}\sum_{i=1}^n \big[T(x_i^*)-B(x_i^*)\big]\De x =\int_a^b \big[T(x)-B(x)\big]\dee{x}. \end{gather*}
Nous pouvons étendre cette procédure pour trouver la masse de notre plaque plutôt que son aire. Nous devons simplement remplacer notre approximation \(\big[T(x_i^*)-B(x_i^*)\big]\De x\) de l’aire de la \(i^{\text{ième}}\) tranche par une approximation de sa masse. Pour cela,
  • nous choisissons un nombre naturel \(m\) (que nous ferons tendre vers l’infini), puis
  • nous divisons la \(i^\text{ième}\) tranche en \(m\) petits rectangles, chacun de largeur \(\De x\) et de hauteur \(\De y = \frac{1}{m}\big[T(x_i^*)-B(x_i^*)\big]\text{.}\) Posons \(y_j = B(x_i^*) + j\,\De y\text{,}\) l’ordonnée de l’extrémité supérieure du \(j^\text{ième}\) rectangle.
  • À ce moment, nous pouvons approcher la densité à l’intérieur de chaque rectangle par une constante (en supposant que la fonction de densité est continue, la densité ne varie pas trop brusquement sur un des petits rectangles). Pour chaque \(j=1,2,\dots,m\text{,}\) le \(j^\text{ième}\) rectangle s’étend de \(y=y_{j-1}\) à \(y= y_j\text{.}\) Nous choisissons au hasard un nombre \(y_j^*\) entre \(y_{j-1}\) et \(y_j\text{,}\) puis nous approximons la densité pour ce \(j^\text{ième}\) rectangle de la \(i^\text{ième}\) tranche par la constante \(f\big(x_i^*,y_j^*\big)\text{.}\)
  • Ainsi, la masse du \(j^\text{ième}\) rectangle dans la \(i^\text{ième}\) tranche est approximativement donnée par \(f\big(x_i^*,y_j^*\big)\,\De x\,\De y\text{.}\)
  • De cette façon, la somme de Riemann ci-après fournit une approximation de la masse de la \(i^\text{ième}\) tranche.
    \begin{align*} \text{Masse } i^\text{ième} \text{tranche} &\approx \sum_{j=1}^m f\big(x_i^*,y_j^*\big)\,\De x\,\De y \end{align*}
    Il faut noter que les nombres \(y_j^*\) dépendent de \(i\) et de \(m\text{.}\)
  • En évaluant la limite lorsque \(m \to \infty\) (de sorte que la hauteur des rectangles tendent vers zéro), nous transformons la somme de Riemann en une intégrale définie :
    \begin{align*} \text{Masse } i^\text{ième} \text{ tranche} &\approx \De x \int_{B(x_i^*)}^{T(x_i^*)} f\big(x_i^*,y\big)\,\dee{y} = F(x_i^*)\,\De x, \end{align*}
    \begin{equation*} F(x) = \int_{B(x)}^{T(x)} f\big(x,y\big)\,\dee{y}. \end{equation*}
    Remarquons que, bien que nous ayons débuté avec une fonction de densité \(f(x,y)\) de deux variables, \(x\) et \(y\text{,}\) en évaluant cette intégrale définie par rapport à la variable \(y\text{,}\) nous obtenons comme résultat \(F(x)\text{,}\) une fonction de \(x\) seulement.
  • Finalement, en évaluant la limite lorsque \(n \to \infty\) (de sorte que la largeur des tranches tendent vers zéro), nous obtenons
    \begin{align*} \text{Masse} &=\lim_{n\rightarrow\infty}\sum_{i=1}^n \De x \int_{B(x_i^*)}^{T(x_i^*)} f\big(x_i^*,y\big)\,\dee{y} = \lim_{n\rightarrow\infty}\sum_{i=1}^n F(x_i^*)\,\De x. \end{align*}
    Nous sommes de retour au contexte familier des fonctions d’une variable. La somme de Riemann \(\sum\limits_{i=1}^n F(x_i^*)\,\De x\) est une approximation de l’intégrale \(\int_a^b F(x)\,\dee{x}\text{.}\) Ainsi,
    \begin{align*} \text{Masse} &= \int_a^b F(x)\,\dee{x} =\int_a^b \left[\int_{B(x)}^{T(x)} f\big(x,y\big)\,\dee{y}\right]\dee{x}. \end{align*}
Cette expression est notre première intégrale double. Il existe d’autres notations courantes pour cette intégrale.

Définition 3.1.1.

\begin{align*} \iint_\cR f\big(x,y\big)\,\dee{x}\,\dee{y} &= \int_a^b \left[\int_{B(x)}^{T(x)} f\big(x,y\big)\,\dee{y}\right]\dee{x}\\ &=\int_a^b \int_{B(x)}^{T(x)} f\big(x,y\big)\,\dee{y}\,\dee{x} =\int_a^b \dee{x} \int_{B(x)}^{T(x)} \dee{y}\, f\big(x,y\big) \end{align*}
Les trois dernières intégrales sont dites des intégrales itérées, pour des raisons évidentes.
Notons que :
  • Pour évaluer l’intégrale \(\displaystyle\int_a^b \int_{B(x)}^{T(x)} f\big(x,y\big)\,\dee{y}\,\dee{x}\text{,}\)
    • nous évaluons d’abord l’intégrale à l’intérieur \(\int_{B(x)}^{T(x)} f\big(x,y\big)\,\dee{y}\) en intégrant \(f(x,y)\) par rapport à \(y\) de \(y=B(x)\) à \(y=T(x)\) et en traitant \(x\) comme une constante. Nous intégrons en utilisant les techniques d’intégration à une variable vues en calcul intégral. Le résultat de cette intégrale est une fonction de \(x\) seulement, nommons-la \(F(x)\text{.}\)
    • Nous évaluons ensuite l’intégrale extérieure \(\int_a^b F(x)\,\dee{x}\text{,}\) dont l’intégrande est la réponse à l’intégrale intérieure. Encore une fois, cette intégrale est évaluée en utilisant les techniques d’intégration à une variable, donc en intégrant \(F(x)\) par rapport à \(x\) de \(x=a\) à \(x=b\text{.}\)
    Nous utiliserons cette notation.
  • Pour évaluer l’intégrale \(\displaystyle\int_a^b \dee{x} \int_{B(x)}^{T(x)} \dee{y}\, f\big(x,y\big)\text{,}\)
    • nous évaluons d’abord l’intégrale à l’intérieur \(\int_{B(x)}^{T(x)} \dee{y}\, f\big(x,y\big)\) en intégrant \(f(x,y)\) par rapport à \(y\) de \(B(x)\) à \(T(x)\text{.}\)
    • Puis, nous évaluons l’intégrale extérieure \(\int_a^b \dee{x}\,F(x)\text{,}\) pour laquelle l’intégrande est le résultat de l’intégrale intérieure et dont les bornes d’intégration vont de \(x=a\) à \(x=b\text{.}\)
    Cette façon d’écrire est aussi courante, elle met l’accent sur le point de vue qui considère l’intégration comme un opérateur.
À ce stade-ci, nous pouvons nous demander si nous devons toujours utiliser des tranches verticales ou si nous devons toujours intégrer par rapport à \(y\) en premier. La réponse à ces deux questions est « non ». Cela nous amène dès lors à considérer les tranches horizontales.

Sous-section Tranches horizontales

Nous avons constaté, lors du calcul d’aire de régions dans le plan \(Oxy\text{,}\) qu’il est souvent plus avantageux d’utiliser des tranches horizontales que des tranches verticales. Il en va de même lors de la mise en place d’intégrales doubles. Nous répétons donc la procédure de la section précédente, mais en commençant par des tranches horizontales plutôt que verticales. Cette procédure sera utile avec des régions de la forme
\begin{equation*} \cR = \big\{\ (x,y)\ \big|\ c\leqslant y\leqslant d,\ L(y)\leqslant x\leqslant R(y)\ \big\}. \end{equation*}
Ici, dans la région \(\cR\text{,}\) \(L(y)\) (\(L\) pour « left ») est la plus petite valeur de \(x\) lorsque l’ordonnée est \(y\text{.}\) \(R(y)\) (\(R\) pour « right ») est la plus grande valeur de \(x\) lorsque l’ordonnée est \(y\text{.}\) Supposons que nous souhaitons évaluer la masse d’une plaque qui remplit la région \(\cR\) et que la densité de la plaque est \(f(x,y)\text{.}\) Nous suivons essentiellement la même procédure que celle que nous avons utilisée avec les coupes verticales, mais en interchangeant les rôles de \(x\) et \(y\text{.}\)
  • Choisissons un nombre naturel \(n\) (que nous ferons tendre vers l’infini), puis
  • divisons l’intervalle \(c\leqslant y\leqslant d\) en \(n\) petits sous-intervalles, chacun de largeur \(\De y=\frac{d-c}{n}\text{.}\) Chaque sous-intervalle définit une mince tranche horizontale de la région (voir la figure ci-dessous).
  • Nous bous approchons de la \(i^\text{ième}\) tranche par un mince rectangle horizontal (le rectangle gris foncé dans la figure ci-dessous). Sur cette tranche, l’ordonnée varie dans un très petit intervalle. Nous choisissons un nombre \(y_i^*\) dans cet intervalle. Nous nous approchons donc de la \(i^\text{ième}\) tranche par un rectangle dont l’extrémité gauche est à \(x=L(y_i^*)\) et dont l’extrémité droite est à \(x=R(y_i^*)\text{.}\)
  • Si nous voulions calculer l’aire de la région \(\cR\text{,}\) il faudrait approximer l’aire de la \(i^\text{ième}\) tranche par \(\big[R(x_i^*)-L(x_i^*)\big]\De y\text{,}\) soit l’aire du rectangle de largeur \(\big[R(x_i^*)-L(x_i^*)\big]\) et de hauteur \(\De y\text{.}\)
  • Pour calculer la masse, similairement à ce que nous avons fait auparavant,
    • choisissons un autre nombre naturel \(m\) (que nous ferons tendre vers l’infini), puis
    • divisons la \(i^\text{ième}\) tranche en \(m\) petits rectangles dont chacun a une hauteur de \(\De y\) et une largeur de \(\De x=\frac{1}{m}\big[R(y_i^*)-L(y_i^*)\big]\text{.}\)
    • Pour chaque \(j=1,2,\dots,m\text{,}\) l’abscisse du \(j^\text{ième}\) rectangle varie dans un petit intervalle. Choisissons un nombre \(x_j^*\) dans cet intervalle (voir le petit rectangle noir de la figure ci-dessous).
      Voici une vue agrandie de la \(i^\text{ième}\) tranche.
    • Sur le \(j^\text{ième}\) rectangle de la \(i^\text{ième}\) tranche, nous approximons la densité par \(f\big(x_j^*,y_i^*\big)\text{.}\) Par conséquent, sa masse est approximativement donnée par \(f\big(x_j^*,y_i^*\big)\,\De x\,\De y\text{.}\)
    • Ainsi, l’approximation de la masse de la \(i^\text{ième}\) tranche horizontale au moyen de la somme de Riemann est
      \begin{align*} \text{Masse } i^{\text{ième}} \text{ tranche} &\approx \sum_{j=1}^m f\big(x_j^*,y_i^*\big)\,\De x\,\De y. \end{align*}
    • En évaluant la limite lorsque \(m \to \infty\) (de sorte que la largeur des rectangles tendent vers zéro), nous transformons la somme de Riemann en une intégrale définie :
      \begin{align*} \text{Masse }i^{\text{ième}} \text{ tranche} &\approx \De y \int_{L(y_i^*)}^{R(y_i^*)} f\big(x,y_i^*\big)\,\dee{x} = F(y_i^*)\,\De y, \end{align*}
      \begin{equation*} F(y) = \int_{L(y)}^{R(y)} f\big(x,y\big)\,\dee{x}. \end{equation*}
      Remarquons que l’intégration par rapport à la variable \(x\) élimine celle-ci : \(F(y)\) est une fonction de \(y\) seulement.
  • Finalement, en évaluant la limite lorsque \(n \to \infty\) (de sorte que la largeur de la tranche tendent vers zéro), nous obtenons
    \begin{align*} \text{Masse} &=\lim_{n\rightarrow\infty}\sum_{i=1}^n \De y \int_{L(y_i^*)}^{R(y_i^*)} f\big(x,y_i^*\big)\,\dee{x} = \lim_{n\rightarrow\infty}\sum_{i=1}^n F(y_i^*)\,\De y. \end{align*}
    Or, la somme de Riemann \(\sum\limits_{i=1}^n F(y_i^*)\,\De y\) est une approximation de l’intégrale \(\int_c^d F(y)\,\dee{y}\text{.}\) Ainsi,
    \begin{align*} \text{Masse} &= \int_c^d F(y)\,\dee{y} =\int_c^d \left[\int_{L(y)}^{R(y)} f\big(x,y\big)\,\dee{x}\right]\dee{y}. \end{align*}
Les notations standards de la Définition 3.1.1 s’appliquent aussi à cette intégrale.

Définition 3.1.2.

\begin{align*} \iint_\cR f\big(x,y\big)\,\dee{x}\,\dee{y} &=\int_c^d \left[\int_{L(y)}^{R(y)} f\big(x,y\big)\,\dee{x}\right]\dee{y}\\ &=\int_c^d \int_{L(y)}^{R(y)} f\big(x,y\big)\,\dee{x}\,\dee{y} =\int_c^d \dee{y} \int_{L(y)}^{R(y)} \dee{x}\, f\big(x,y\big) \end{align*}
Notons que :
  • Pour évaluer l’intégrale \(\displaystyle\int_c^d \int_{L(y)}^{R(y)} f\big(x,y\big)\,\dee{x}\,\dee{y}\text{,}\)
    • nous évaluons en premier l’intégrale à l’intérieur \(\int_{L(y)}^{R(y)} f\big(x,y\big)\,\dee{x}\) en intégrant \(f(x,y)\) par rapport à \(x\) de \(x=L(y)\) à \(x=R(y)\text{.}\) Le résultat de cette intégrale est une fonction de \(y\) seulement. Nommons-la \(F(y)\text{.}\)
    • Puis, nous évaluons l’intégrale à l’extérieur \(\int_c^d F(y)\,\dee{y}\text{,}\) dont l’intégrande est le résultat de l’intégrale intérieure.
  • Pour évaluer l’intégrale \(\displaystyle\int_c^d \dee{y} \int_{L(y)}^{R(y)} \dee{x}\, f\big(x,y\big)\text{,}\)
    • nous évaluons en premier l’intégrale \(\int_{L(y)}^{R(y)} \dee{x}\, f\big(x,y\big)\) en utilisant les bornes d’intégration qui se trouvent à la gauche de \(\dee{x}\text{.}\) “ \(\dee{x}\) ” est écrit directement à côté de \(\int_{L(y)}^{R(y)}\) pour mettre en évidence que les bornes d’intégration \(L(y)\) et \(R(y)\) sont pour l’intégrale par rapport à \(x\text{.}\) Cette intégrale définie est habituellement écrite \(\int_{L(y)}^{R(y)} f\big(x,y\big)\ \dee{x}\text{.}\) Le résultat de l’intégrale intérieure est de nouveau une fonction de \(y\) seulement, que nous nommons \(F(y)\text{.}\)
    • Puis, nous évaluons l’intégrale à l’extérieur \(\int_c^d \dee{y}\,F(y)\text{,}\) dont l’intégrande est le résultat de l’intégrale intérieure et dont les bornes d’intégration sont à la gauche de \(\dee{y}\text{.}\)
En résumé, nous avons maintenant deux intégrales qui calculent la masse de la même plaque dans le plan \(Oxy\text{.}\)
Nous avons, implicitement, dans le Théorème 3.1.3, que si
\begin{align*} &\big\{\ (x,y)\ \big|\ a\leqslant x\leqslant b,\ B(x)\leqslant y\leqslant T(x)\ \big\}\\ =& \big\{\ (x,y)\ \big|\ c\leqslant y\leqslant d,\ L(y)\leqslant x\leqslant R(y)\ \big\} \end{align*}
et que si \(f(x,y)\) est continue, alors
\begin{equation*} \int_a^b \int_{B(x)}^{T(x)} f\big(x,y\big)\,\dee{y}\,\dee{x} =\int_c^d \int_{L(y)}^{R(y)} f\big(x,y\big)\,\dee{x}\,\dee{y}. \end{equation*}
Ce résultat est connu sous le nom de théorème de Fubini
 1 
En l’honneur du mathématicien italien Guido Fubini (1879-1943).
. Il sera abordé en détail plus loin.

Définition 3.1.4.

Les intégrales du Théorème 3.1.3 sont notées
\begin{equation*} \iint_\cR f(x,y)\,\dee{x}\dee{y}\text{ ou }\iint_\cR f(x,y)\,\dee{y}\dee{x} \text{ ou } \iint_\cR f(x,y)\,\dee{A}. \end{equation*}
Le symbole \(\dee{A}\) représente l’aire d’une partie infinitésimale de \(\cR\text{.}\)
Voyons un exemple simple. Nous en verrons des plus complexes dans une section ultérieure.

Exemple 3.1.5.

Soit \(\cR\) la région triangulaire située au-dessus de l’axe des \(x\text{,}\) à droite de l’axe des \(y\) et à gauche de la droite \(x+y=1\text{.}\) Trouver la masse de \(\cR\) si la densité est donnée par \(f(x,y)=y\text{.}\)
Solution.
Nous ferons ce calcul de deux façons : avec les tranches verticales, ensuite avec les tranches horizontales. Voici d’abord une esquisse de \(\cR\text{.}\)
Solution avec des tranches verticales.\(\ \ \ \) Nous allons écrire une intégrale double pour calculer la masse en utilisant des rectangles verticaux. Remarquons, dans la figure suivante,
  • que les points le plus à gauche dans \(\cR\) ont comme abscisse \(x=0\) et que le point le plus à droite de \(\cR\) a \(x=1\text{.}\)
  • De plus, si nous choisissons un \(x\) fixe entre \(0\) et \(1\text{,}\) le point \((x,y)\) de \(\cR\) avec la plus petite valeur de \(y\) a \(y=0\text{,}\) et celui avec la plus grande valeur de \(y\) a \(y=1-x\text{.}\)
Ainsi,
\begin{equation*} \cR=\Set{(x,y)}{0=a\leqslant x\leqslant b=1,\ 0= B(x)\leqslant y\leqslant T(x) = 1-x} \end{equation*}
et, en vertu de la partie (a) du Théorème 3.1.3,
\begin{align*} \text{Masse} &= \int_a^b \left[ \int_{B(x)}^{T(x)} f\big(x,y\big)\, \dee{y}\right]\, \dee{x} = \int_0^1 \left[ \int_0^{1-x} y \, \dee{y}\right] \dee{x}. \end{align*}
On a que l’intégrale intérieure est
\begin{gather*} \int_0^{1-x} y \ \dee{y} =\left[ \frac{y^2}{2} \right]_0^{1-x} =\frac{1}{2}(1-x)^2 , \end{gather*}
de sorte que la masse vaut
\begin{align*} \text{Masse} &= \int_0^1 \frac{(1-x)^2}{2}\, \dee{x} =\left[-\frac{(1-x)^3}{6}\right]_0^1 = \frac{1}{6}. \end{align*}
Solution avec des tranches horizontales.\(\ \ \ \) Faisons maintenant le calcul dans l’autre ordre d’intégration, soit en utilisant des rectangles horizontaux. De nouveau, commençons par une esquisse.
On remarque, dans la figure, que
  • les points les plus bas dans \(\cR\) ont comme ordonnée \(y=0\) et que le point le plus haut a \(y=1\text{.}\)
  • De plus, pour chaque valeur de \(y\) fixe entre \(0\) et \(1\text{,}\) le point \((x,y)\) de \(\cR\) avec la valeur la plus petite de \(x\) a \(x=0\text{,}\) et celui avec la valeur la plus grande de \(x\) a \(x=1-y\text{.}\)
Ainsi,
\begin{equation*} \cR=\Set{(x,y)}{0=c\leqslant y\leqslant d=1,\ 0= L(y)\leqslant x\leqslant R(y) = 1-y} \end{equation*}
et, en vertu de la partie (b) du Théorème 3.1.3,
\begin{align*} \text{Masse} &= \int_c^d \left[ \int_{L(y)}^{R(y)}f\big(x,y\big) \ \dee{x}\right] \dee{y} \, = \int_0^1 \left[ \int_0^{1-y} y \, \dee{x} \right] \dee{y}. \end{align*}
On a que l’intégrale intérieure est
\begin{gather*} \int_0^{1-y} y \ \dee{x} =\left[ xy \right]_0^{1-y} =y-y^2 \end{gather*}
puisque l’intégration par rapport à \(y\) considère \(x\) comme une constante. Par conséquent,
\begin{align*} \text{Masse} &= \int_0^1 \big[y-y^2\big]\, \dee{y}\ =\left[\frac{y^2}{2}-\frac{y^3}{3}\right]_0^1 = \frac{1}{2}-\frac{1}{3} =\frac{1}{6}. \end{align*}
Les intégrales doubles vérifient les mêmes propriétés de base que les intégrales simples, telles qu’énumérées dans le théorème suivant.
Dans le cas particulier (mais pas si inhabituel) où \(\cR\) est un rectangle
\begin{equation*} \cR=\Set{(x,y)}{a\leqslant x\leqslant b,\ c\leqslant y\leqslant d} \end{equation*}
et où l’intégrande est un produit \(f(x,y)=g(x)h(y)\text{,}\)
\begin{align*} \iint_\cR f(x,y)\,\dee{A} &=\int_a^b \int_c^d g(x) h(y)\,\dee{y} \dee{x}\\ &=\int_a^b g(x)\, \left[ \int_c^d h(y)\ \dee{y}\right]\dee{x},\\ &\hskip0.25in\text{puisque } g(x) \text{ est constante lorsqu'on intègre par rapport à } y,\\ &=\left[ \int_a^b g(x) \ \dee{x}\right]\ \left[\int_c^d h(y) \ \dee{y}\right],\\ &\hskip0.25in\text{puisque } \int_c^d h(y)\, \dee{y} \text{ est une constante lorsqu'on intègre par rapport à } x. \end{align*}
Cela mérite d’être énoncé dans un théorème.

Exemple 3.1.8.

Soit \(\cR = [0,2] \times [0,1]\text{.}\) Calculer la valeur de \(\displaystyle \iint_{\cR} \frac{yx^3}{2+y^2}\, \dee{A}\text{.}\)
Solution.
Notons que l’intégrande, à savoir la fonction \(f(x,y) = \frac{x^3y}{2+y^2}\text{,}\) peut s’écrire comme produit
\begin{equation*} f(x,y) = g(x)\, h(y) \text{ avec } g(x) = x^3, h(y) = \frac{y}{2+y^2}. \end{equation*}
De plus, le domaine d’intégration est un rectangle. Il vient donc que
\begin{align*} \iint_{\cR} \frac{x^3y}{2+y^2}\, \dee{A} \amp = \left[\int_0^2 x^3 \, \dee{x}\right] \left[\int_0^1 \frac{y}{2+y^2}\, \dee{y}\right]\\ \amp = \left.\frac{1}{4} x^4\right|_{x=0}^{x=2} \left. \times \frac{1}{2} \ln (2+y^2)\right|_{y=0}^{y=1}\\ \amp = \frac{1}{4}(16-0)\cdot \frac{1}{2}(\ln 3 - \ln 2) =2\ln(3/2). \end{align*}
Comme pour les intégrales à une variable, il n’est pas toujours nécessaire de connaître la valeur exacte d’une intégrale double. Nous nous intéressons plutôt à trouver des bornes à sa valeur. Le théorème suivant fournit quelques outils simples pour obtenir de telles bornes. Ces propriétés des intégrales à plusieurs variables sont analogues à celles des intégrales à une variable.

Sous-section Volumes

Voyons maintenant quelques applications des intégrales doubles que nous venons de définir. Le problème du calcul de l’aire de la région sous une courbe et au-dessus de l’axe des \(x\) dans le plan, avec des intégrales simples, se transforme en celui du calcul du volume du solide sous une surface et au-dessus du plan \(Oxy\text{,}\) avec des intégrales doubles. Cette fois, au lieu d’interpréter la fonction \(f(x,y)\) comme une densité, t comme la hauteur d’un solide se trouvant au-dessus d’une région \(\cR\) du plan \(Oxy\text{.}\)
Voyons comment cela s’applique dans la partie (b) du Théorème 3.1.3. Le cas de la partie (a) se traitera de la même façon. Supposons un solide \(\cV\) défini au-dessus de la région suivante,
\begin{equation*} \cR = \big\{\ (x,y)\ \big|\ c\leqslant y\leqslant d,\ L(y)\leqslant x\leqslant R(y)\ \big\}, \end{equation*}
et tel que
\begin{equation*} \cV = \big\{\ (x,y,z)\ \big|\ (x,y)\in\cR,\ 0\leqslant z\leqslant f(x,y)\ \big\}. \end{equation*}
La région \(\cR\) (qui est aussi la région vue au-dessus du solide \(\cV\)) est illustrée dans la figure de gauche, et la portion de \(\cV\) dans le premier octant est illustrée dans la figure de droite.
Afin de calculer le volume de \(\cV\text{,}\) nous allons
  • choisir un nombre naturel \(n\) et découper \(\cR\) en rectangles de largeur \(\De y=\frac{d-c}{n}\text{;}\)
  • diviser le rectangle \(i\) en \(m\) sous-rectangles, chacun de largeur \(\De y\) et de hauteur \(\De x=\frac{1}{m}\cdots\text{;}\)
  • calculer, approximativement, le volume de la portion de \(\cV\) qui se trouve au-dessus de chaque rectangle;
  • calculer la limite lorsque \(m\rightarrow\infty\text{,}\) puis la limite lorsque \(n\rightarrow\infty\text{.}\)
Étant donné que nous avons déjà fait ce type de raisonnement à deux reprises, nous sauterons quelques étapes à notre raisonnement :
  • Il faut découper la région \(\cR\) en rectangles de largeur infinitésimale \(\dee{y}\text{.}\)
  • Il faut diviser chaque rectangle en sous-rectangles « infinitésimaux » de largeur \(\dee{y}\) et de hauteur \(\dee{x}\text{.}\) Voir la figure de gauche ci-dessus.
  • Le volume de la portion de \(\cV\) qui se trouve au-dessus du sous-rectangle centré en \((x,y)\) est approximativement égal à \(f(x,y)\,\dee{x}\,\dee{y}\text{.}\) Voir la figure de droite ci-dessus.
  • Ainsi, le volume de la portion de \(\cV\) qui se trouve au-dessus du rectangle centré en \(y\) est approximativement égal à
     2 
    Pensez à la portion de \(\cV\) qui se trouve au-dessus de ce rectangle comme à une fine tranche de pain. Ainsi, le facteur \(\dee{y}\) dans l’expression \(\dee{y}\int_{L(y)}^{R(y)} f(x,y) \ \dee{x}\) est l’épaisseur de la tranche. Le facteur \(\int_{L(y)}^{R(y)} f(x,y) \ \dee{x}\) est l’aire de la surface transversale \(\Set{(x,z)}{L(y)\leqslant x\leqslant R(y),\ 0\leqslant z\leqslant f(x,y)}\text{,}\) c.-à-d. l’aire de la surface de la tranche de pain.
    \(\dee{y}\int_{L(y)}^{R(y)}f(x,y) \ \dee{x}\text{.}\)
  • Nous arrivons donc à la conclusion suivante :
Similairement :

Sous-section Exemples

Maintenant que nous avons vu beaucoup de formules et de théorie, il est temps de mettre tout cela en pratique. Commençons par un exemple de calcul de masse, puis poursuivons avec un exemple de calcul de volume. Vous remarquerez dans les deux cas que les calculs mathématiques sont vraiment similaires. Seulement l’interprétation du résultat change.

Exemple 3.1.12. Masse.

Soit \(\nu \gt 0\) une constante et \(\cR\) la région au-dessus de la courbe \(x^2=4\nu y\) et à la droite de la courbe \(y^2=\frac{1}{2}\nu x\text{.}\) Trouver la masse de \(\cR\) si sa densité est \(f(x,y)=xy\text{.}\)
Solution.
Pour illustrer les méthodes, nous présentons deux solutions : avec des tranches verticales, puis avec des tranches horizontales. Pour commencer, esquissons la région \(\cR\text{.}\) Notons d’abord que, puisque \(y\ge\frac{x^2}{4\nu}\) et \(x\geqslant \frac{2y^2}{\nu}\text{,}\) les valeurs de \(x\) et de \(y\) sont positives sur \(\cR\text{.}\) Les deux courbes se coupent aux points \((x,y)\) vérifiant les équations
\begin{align*} x = \frac{2y^2}{\nu}\text{ et }y=\frac{x^2}{4\nu} &\implies\quad x = \frac{2y^2}{\nu} = \frac{2}{\nu}\left(\frac{x^2}{4\nu}\right)^2 = \frac{x^4}{8\nu^3}\\ &\iff\quad \left(\frac{x^3}{8\nu^3}-1\right)x=0. \end{align*}
Cette équation admet deux solutions réelles, soit \(x=0\) et \(x=2\nu\text{.}\) Ainsi, la parabole \(y = \frac{x^2}{4\nu}\text{,}\) dont l’ouverture est vers le haut, et la parabole \(x=\frac{2y^2}{\nu}\text{,}\) dont l’ouverture est vers la droite, se coupent en \((0,0)\) et \((2\nu,\nu)\text{.}\)
Solution avec des tranches verticales. Écrivons maintenant une intégrale double qui calcule la masse en utilisant des tranches verticales. Nous allons utiliser le raisonnement abrégé vu à la section précédente sur les volumes. Déduisons de la figure ci-dessus que
\begin{equation*} \cR=\bigg\{(x,y)\ \bigg|\ 0=a\leqslant x\leqslant b=2\nu,\ \frac{x^2}{4\nu}= B(x)\leqslant y\leqslant T(x) = \sqrt{\frac{\nu x}{2}}\ \bigg\}. \end{equation*}
  • Découpons \(\cR\) en rectangles verticaux de largeur infinitésimale \(\dee{x}\text{.}\)
  • Divisons chaque rectangle en sous-rectangles de hauteur infinitésimale \(\dee{y}\) et de largeur \(\dee{x}\text{.}\) Voir la figure ci-dessous.
  • La masse du sous-rectangle centré en \((x,y)\) est approximativement égale à \(f(x,y)\,\dee{x}\,\dee{y}=xy\,\dee{x}\,\dee{y}\text{.}\)
  • Ainsi, la masse du rectangle centré en \(x\) est approximativement égale à \(\ \left[\int_{B(x)}^{T(x)} f(x,y)\, \dee{y}\right]\dee{x} \) (l’intégrale par rapport à \(y\) additionne les masses de tous les sous-rectangles sur un même rectangle vertical).
  • Nous concluons que
    \begin{align*} \text{Masse}(\cR) &=\int_a^b \int_{B(x)}^{T(x)} f(x,y)\, \dee{y}\ \dee{x} =\int_0^{2\nu} \int_{x^2/(4\nu)}^{\sqrt{\nu x/2}} xy \, \dee{y}\ \dee{x}. \end{align*}
    Ici, l’intégrale par rapport à \(x\) additionne les masses de chacun des rectangles verticaux.
    Rappelons que lorsque nous intégrons par rapport à \(y\text{,}\) \(x\) est une constante, de sorte que nous pouvons factoriser la constante \(x\) à l’extérieur de l’intégrale par rapport à \(y\text{.}\)
    \begin{align*} \int_{x^2/(4\nu)}^{\sqrt{\nu x/2}} xy \ \dee{y} &= x\int_{x^2/(4\nu)}^{\sqrt{\nu x/2}} y \ \dee{y}\\ &= x\left[\frac{y^2}{2}\right]_{x^2/(4\nu)}^{\sqrt{\nu x/2}}\\ &=\frac{\nu x^2}{4}-\frac{x^5}{32\nu^2}, \end{align*}
    et alors
    \begin{align*} \text{Masse}(\cR) &=\int_0^{2\nu} \left[\frac{\nu x^2}{4}-\frac{x^5}{32\nu^2}\right] \ \dee{x}\\\ &= \frac{\nu (2\nu)^3}{3\times 4}-\frac{(2\nu)^6}{6\times 32\nu^2} =\frac{\nu^4}{3}. \end{align*}
Avec des tranches horizontales. Déduisons de la figure au début de cet exemple que
\begin{equation*} \cR=\bigg\{(x,y)\ \bigg|\ 0=c\leqslant y\leqslant d=\nu,\ \frac{2 y^2}{\nu} = L(y)\leqslant x\leqslant R(y) = \sqrt{4\nu y} \ \bigg\}. \end{equation*}
  • Découpons \(\cR\) en rectangles horizontaux de hauteur infinitésimale \(\dee{y}\text{.}\)
  • Divisons chaque rectangle en sous-rectangles de hauteur \(\dee{y}\) et de largeur infinitésimale \(\dee{x}\text{.}\) Voir la figure ci-dessous.
  • La masse du sous-rectangle centré en \((x,y)\) est approximativement égale à \(f(x,y)\,\dee{x}\,\dee{y}=xy\,\dee{x}\,\dee{y}\text{.}\)
  • Ainsi, la masse du rectangle horizontal centré en \(y\) est approximativement \(\left[ \int_{L(y)}^{R(y)} f(x,y)\ \dee{x}\ \right] \dee{y} \) (l’intégrale par rapport à \(x\) additionne les masses de chacun des sous-rectangles sur un même rectangle horizontal).
  • Nous concluons que
    \begin{align*} \text{Masse}(\cR) &=\int_c^d \int_{L(y)}^{R(y)}\ f(x,y)\, \dee{x} \ \dee{y} =\int_0^{\nu} \int_{2y^2/\nu}^{\sqrt{4\nu y}} xy\ \dee{x} \ \dee{y}. \end{align*}
    Ici, l’intégrale par rapport à \(y\) additionne les masses de chacun des rectangles horizontaux. Lorsque nous intégrons par rapport à \(x\text{,}\) \(y\) est constante. Ainsi,
    \begin{align*} \text{Masse}(\cR) &=\int_0^{\nu} \ y\left[ \int_{2y^2/\nu}^{\sqrt{4\nu y}} \dee{x}\ x\right]\ \dee{y}\\ &=\int_0^{\nu} y\left[\frac{x^2}{2}\right] _{2y^2/\nu}^{\sqrt{4\nu y}}\ \dee{y}\\\ &=\int_0^{\nu} \left[2\nu y^2-\frac{2y^5}{\nu^2}\right]\ \dee{y}\\\ &= \frac{2\nu (\nu)^3}{3}-\frac{2\nu^6}{6\nu^2} =\frac{\nu^4}{3}. \end{align*}

Exemple 3.1.13. Volume.

Soit \(\cR\) la région du plan \(Oxy\) qui se trouve au-dessus de l’axe des \(x\) et sous la parabole \(y=1-x^2\text{.}\) Calculer le volume du solide au-dessus de la région \(\cR\) et sous la surface \(z=x^2\sqrt{1-y}\text{.}\)
Solution.
Nous allons présenter en détail la solution avec des tranches verticales. Pour les tranches horizontales, nous nous contenterons de poser l’intégrale. De nouveau, commençons par une représentation de \(\cR\text{.}\)
Avec les tranches verticales.\(\ \ \ \) Nous pouvons déduire de la figure ci-dessous
que
  • le point le plus à gauche de \(\cR\) a pour abscisse \(x=-1\) et que celui le plus à droite a \(x=1\text{;}\)
  • pour chaque valeur de \(x\) fixe entre \(-1\) et \(1\text{,}\) le point \((x,y)\) dans \(\cR\) avec la plus petite ordonnée a \(y=0\) et que celui avec la plus grande ordonnée vérifie \(y=1-x^2\text{.}\)
Ainsi,
\begin{equation*} \cR=\Set{(x,y)}{-1=a\leqslant x\leqslant b=1,\ 0= B(x)\leqslant y\leqslant T(x) = 1-x^2} \end{equation*}
et, d’après la Formule 3.1.11,
\begin{align*} \text{Volume} &= \int_a^b \int_{B(x)}^{T(x)} f\big(x,y\big)\ \dee{y}\, \dee{x} = \int_{-1}^1 \int_0^{1-x^2} x^2\sqrt{1-y}\ \dee{y}\ \,\dee{x}\\ &=2\int_0^1 \int_0^{1-x^2} \ x^2\sqrt{1-y}\ \dee{y}\,\dee{x}, \end{align*}
car l’intégrale \(F(x) = \int_0^{1-x^2} x^2\sqrt{1-y}\ \dee{y}\) est une fonction paire de \(x\text{.}\) Or, pour \(x\geqslant 0\text{,}\) l’intégrale intérieure est
\begin{align*} \int_0^{1-x^2} x^2\sqrt{1-y} \ \dee{y} &= x^2 \int_0^{1-x^2} \sqrt{1-y} \ \dee{y} =x^2 \left[ -\frac{2}{3}(1-y)^{3/2} \right]_0^{1-x^2}\\ &=\frac{2}{3}x^2\big(1-x^3\big), \end{align*}
de sorte que
\begin{align*} \text{Volume} &= 2\int_0^1 \dee{x}\ \frac{2}{3}x^2\big(1-x^3\big) =\frac{4}{3}\left[\frac{x^3}{3}-\frac{x^6}{6}\right]_0^1 =\frac{2}{9}. \end{align*}
Avec les tranches horizontales.\(\ \ \ \) Avec cette méthode, nous obtenons la figure ci-dessous.
La région d’intégration est
\begin{equation*} \cR=\Set{(x,y)}{0=c\leqslant y\leqslant d=1,\ -\sqrt{1-y}= L(y)\leqslant x\leqslant R(y) = \sqrt{1-y}}. \end{equation*}
Donc, d’après la Formule 3.1.10,
\begin{align*} \text{Volume} &= \int_c^d \int_{L(y)}^{R(y)}\ f\big(x,y\big)\ \dee{x}\, \dee{y} = \int_0^1 \int_{-\sqrt{1-y}}^{\sqrt{1-y}} x^2\sqrt{1-y} \ \dee{x}\,\dee{y}. \end{align*}
Nous vous laissons le soin de vérifier que le résultat est le même que celui obtenu précédemment.

Remarque 3.1.14.

Il est important de noter que, dans les figures que nous avons faites, nous avons concentré notre attention sur la région d’intégration. L’intégrande \(f(x,y)\) n’a joué aucun rôle dans cette figure, nous l’avons considéré seulement lors de l’évaluation des intégrales.

Exemple 3.1.15. Volume.

Calculer le volume de la partie commune aux deux cylindres \(x^2+y^2=a^2\) et \(x^2+z^2=a^2\text{.}\)
Solution.
La première étape est de bien visualiser le solide spécifié. Remarquons que :
  • La variable \(z\) n’apparaît pas dans l’équation \(x^2+y^2=a^2\text{,}\) de sorte que, pour toute valeur \(z_0\text{,}\) la partie du cylindre \(x^2+y^2=a^2\) dans le plan \(z=z_0\) est le cercle \(x^2+y^2=a^2\text{,}\) \(z=z_0\text{.}\) Ainsi, le cylindre \(x^2+y^2=a^2\) est formé d’une infinité de cercles de rayon \(a\) superposés les uns par-dessus les autres. La partie du cylindre \(x^2+y^2=a^2\) qui se trouve au-dessus du plan \(Oxy\) est illustrée dans la figure de gauche ci-dessous.
  • La variable \(y\) n’apparaît pas dans \(x^2+z^2=a^2\text{.}\) Ainsi, pour toute valeur constante \(y_0\text{,}\) la partie du cylindre \(x^2+z^2=a^2\) dans le plan \(y=y_0\) est le cercle \(x^2+z^2=a^2\text{,}\) \(y=y_0\text{.}\) Ainsi, le cylindre \(x^2+z^2=a^2\) est formé d’une infinité de cercles superposés horizontalement, les uns à côté des autres. La partie du cylindre \(x^2+z^2=a^2\) qui se trouve à la droite du plan \(Oxz\) est illustrée dans la figure de droite ci-dessous.
Nous devons calculer le volume de partie commune à ces deux cylindres.
  • Les équations \(x^2+y^2=a^2\) et \(x^2+z^2=a^2\) restent inchangées si nous substituons \(x\) par \(-x\text{.}\) En conséquence, les deux cylindres, et donc le solide qui nous intéresse, sont symétriques par rapport au plan \(Oyz\text{.}\) En particulier, le volume de la partie du solide dans l’octant \(x\leqslant 0\text{,}\) \(y\geqslant 0\text{,}\) \(z\geqslant 0\) est le même que celui de la partie se trouvant dans le premier octant \(x\geqslant 0\text{,}\) \(y\geqslant 0\text{,}\) \(z\geqslant 0\text{.}\) De la même façon, les équations ne changent pas si \(y\) est remplacée par \(-y\text{,}\) ou si \(z\) est remplacée par \(-z\text{.}\) Notre solide est donc également symétrique par rapport au plan \(Oxz\) et par rapport au plan \(Oxy\text{.}\) Ainsi, le volume du solide est le même dans chacun des huit octants.
  • Nous allons donc calculer le volume du solide dans le premier octant seulement, puis multiplier par huit.
La partie du solide qui se trouve dans le premier octant est illustrée dans la figure de gauche ci-dessous. Un point \((x,y,z)\) se trouve dans le premier cylindre si et seulement si \(x^2+y^2\leqslant a^2\text{.}\)
Ce point se trouve dans le second cylindre si et seulement si \(x^2+z^2\leqslant a^2\text{.}\) Ainsi, la partie du solide qui se trouve dans le premier octant est
\begin{align*} \cV_1&=\Set{(x,y,z)}{x\geqslant 0,\ y\geqslant 0,\ z\geqslant 0,\ x^2+y^2\leqslant a^2,\ x^2+z^2\leqslant a^2}. \end{align*}
Notons que, dans \(\cV_1\text{,}\) \(z^2\leqslant a^2-x^2\text{,}\) de sorte que \(z\le\sqrt{a^2-x^2}\text{.}\) Ainsi,
\begin{align*} \cV_1 &=\Set{(x,y,z)}{x\geqslant 0,\ y\geqslant 0,\ x^2+y^2\leqslant a^2,\ 0\leqslant z\leqslant \sqrt{a^2-x^2}}. \end{align*}
Une vue d’en haut de la partie du solide dans le premier octant est illustrée dans la figure de droite ci-dessus. On voit notamment que \(x\) varie de \(0\) à \(a\) et que, pour chaque \(x\) fixe, \(y\) varie de \(0\) à \(\sqrt{a^2-x^2}\text{.}\) On peut donc redéfinir
\begin{gather*} \cV_1=\Set{(x,y,z)}{(x,y)\in\cR,\ 0\leqslant z\leqslant f(x,y)}, \end{gather*}
\begin{equation*} \cR=\Big\{\ (x,y)\ \Big|\ 0\leqslant x\leqslant a,\ 0\leqslant y\leqslant \sqrt{a^2-x^2}\ \Big\} \quad\text{ et }\quad f(x,y) = \sqrt{a^2-x^2}. \end{equation*}
Remarquons aussi que \(f(x,y)\) ne dépend pas de \(y\text{.}\) Cela facilitera les calculs suivants.
Nous pouvons calculer le volume de \(\cV_1\) avec notre démarche habituelle.
  • On découpe \(\cR\) en rectangles horizontaux, de largeur infinitésimale \(\dee{x}\text{.}\)
  • On divise chaque rectangle en sous-rectangles de hauteur infinitésimale \(\dee{y}\) et de largeur \(\dee{x}\text{.}\) Voir la figure ci-après.
  • Le volume de la partie de \(\cV_1\) au-dessus du sous-rectangle centré en \((x,y)\) est approximativement égal à
    \begin{equation*} f(x,y)\,\dee{x}\,\dee{y}=\sqrt{a^2-x^2}\ \dee{x}\,\dee{y}. \end{equation*}
  • Ainsi, le volume de la partie de \(\cV_1\) au-dessus du rectangle centré en \(x\) est approximativement
    \begin{equation*} \left[\int_0^{\sqrt{a^2-x^2}} \sqrt{a^2-x^2}\ \dee{y}\right] \dee{x}. \end{equation*}
  • On conclut que
    \begin{align*} \text{Volume}(\cV_1) &=\int_0^a \int_0^{\sqrt{a^2-x^2}} \sqrt{a^2-x^2}\ \dee{y}\ \, \dee{x}. \end{align*}
    En gardant à l’esprit que l’intégration par rapport à \(y\) considère \(x\) comme une constante, calculons le volume.
    \begin{align*} \text{Volume}(\cV_1) &=\int_0^a \sqrt{a^2-x^2}\left[ \int_0^{\sqrt{a^2-x^2}} \dee{y} \right] \ \dee{x}\\ &=\int_0^a \big(a^2-x^2\big)\ \dee{x}\\ &=\left[a^2x-\frac{x^3}{3}\right]_0^a\\ &=\frac{2a^3}{3} \end{align*}
    Ainsi, le volume total du solide recherché est
    \begin{equation*} \text{Volume}(\cV) =8\,\text{Volume}(\cV_1) =\frac{16a^3}{3}. \end{equation*}

Exemple 3.1.16. Interprétation géométrique.

Évaluer \(\displaystyle\int_0^2\int_0^a \sqrt{a^2-x^2}\ \dee{x}\,\dee{y}\text{.}\)
Solution.
Cette intégrale représente le volume d’un solide commun, de sorte qu’elle peut être évaluée sans le calcul d’une intégrale. Le domaine d’intégration est
\begin{equation*} \cR=\Set{(x,y)}{0\leqslant y\leqslant 2,\ 0\leqslant x\leqslant a}, \end{equation*}
et l’intégrande est \(f(x,y) = \sqrt{a^2-x^2}\text{.}\) Par conséquent, l’intégrale représente le volume du solide compris entre le plan \(Oxy\) et la surface \(z=\sqrt{a^2-x^2}\text{,}\) avec \((x,y) \in \cR\text{.}\) On peut réécrire l’équation de la surface comme \(x^2+z^2=a^2\text{,}\) qui représente, comme cela a été vu dans l’Exemple 3.1.15, un cylindre de rayon \(a\) dont l’axe est l’axe des \(y\text{.}\) Nous cherchons le volume de la partie de ce cylindre qui se trouve au-dessus de \(\cR\text{,}\) comme illustré dans la figure ci-dessous.
Les sections du solide que coupe tout plan vertical \(y = y_0\) sont toutes des quarts de cercle de rayon \(a\text{,}\) de sorte que leur aire est égale à \(\frac{1}{4}\pi a^2\text{.}\) L’intégrale intérieure \(\int_0^a \sqrt{a^2-x^2}\ \dee{x}\) donne précisément cette aire. Ainsi, lorsque \(y\) varie de \(0\) à \(2\text{,}\)
\begin{equation*} \int_0^2\int_0^a \sqrt{a^2-x^2}\ \dee{x}\,\dee{y} =\frac{1}{4}\pi a^2\times 2 =\frac{\pi a^2}{2}. \end{equation*}

Exemple 3.1.17. Exemple 3.1.16, démarche complexe.

Il est possible, mais fastidieux, d’évaluer l’intégrale \(\int_0^2\int_0^a \sqrt{a^2-x^2}\ \dee{x}\,\dee{y}\) de l’Exemple 3.1.16 avec les techniques d’intégration des fonctions à une variable. Faisons ce calcul à titre de révision.
L’intégrale à l’intérieur est \(\int_0^a \sqrt{a^2-x^2}\ \dee{x}\text{.}\) La stratégie habituelle afin d’éliminer des expressions contenant des racines carrées dans l’intégrande, comme \(\sqrt{a^2-x^2}\text{,}\) est d’effectuer une substitution trigonométrique. Dans ce cas-ci, le changement de variables approprié est de poser
\begin{equation*} x=a\sin\theta,\qquad \dee{x}=a\cos\theta\ \dee{\theta}. \end{equation*}
La borne d’intégration inférieure \(x=0\text{,}\) soit \(a\sin\theta=0\text{,}\) équivaut à \(\theta=0\text{,}\) et la borne supérieure \(x=a\text{,}\) soit \(a\sin\theta=a\text{,}\) équivaut à \(\theta=\frac{\pi}{2}\text{.}\) Ainsi,
\begin{align*} \int_0^a \sqrt{a^2-x^2}\ \dee{x} &=\int_0^{\pi/2} \sqrt{\underbrace{a^2-a^2\sin^2\theta}_{a^2\cos^2\theta}} \ a\cos\theta\ \dee{\theta} =a^2\int_0^{\pi/2} \cos^2\theta\ \dee{\theta}. \end{align*}
La méthode usuelle pour évaluer cette intégrale fait appel à l’identité trigonométrique
\begin{equation*} \cos(2\theta) = 2\cos^2\theta -1,\qquad\text{qu'on réécrit}\qquad \cos^2\theta =\frac{1+\cos(2\theta)}{2}, \end{equation*}
et alors
\begin{align*} \int_0^a \sqrt{a^2-x^2}\ \dee{x} &=a^2\int_0^{\pi/2} \cos^2\theta\ \dee{\theta} =\frac{a^2}{2} \int_0^{\pi/2} \big[1+\cos(2\theta)\big]\ \dee{\theta}\\ &=\frac{a^2}{2}\left[\theta+\frac{\sin(2\theta)}{2}\right]_0^{\pi/2}\\ &=\frac{\pi a^2}{4}. \end{align*}
Remarquons cependant qu’il existe une autre méthode efficace, quoiqu’un peu recherchée, pour évaluer cette intégrale. Regardons les figures ci-dessous.
Nous observons que
\begin{equation*} \int_0^{\pi/2} \cos^2\theta\ \dee{\theta} =\int_0^{\pi/2} \sin^2\theta\ \dee{\theta}. \end{equation*}
Par conséquent,
\begin{align*} \int_0^{\pi/2} \cos^2\theta\ \dee{\theta} &=\int_0^{\pi/2} \sin^2\theta\ \dee{\theta} =\int_0^{\pi/2} \frac{1}{2}\big[\sin^2\theta+\cos^2\theta\big]\,\dee{\theta}\\ &=\frac{1}{2}\int_0^{\pi/2} \dee{\theta} =\frac{\pi}{4}. \end{align*}
Dans tous les cas, nous avons que l’intégrale intérieure est égale à
\begin{equation*} \int_0^a \sqrt{a^2-x^2}\ \dee{x} = a^2\int_0^{\pi/2} \cos^2\theta\ \dee{\theta} = \frac{\pi a^2}{4} \end{equation*}
et que l’intégrale de départ est égale à
\begin{equation*} \int_0^2\int_0^a \sqrt{a^2-x^2}\ \dee{x}\,\dee{y} =\frac{\pi a^2}{4}\int_0^2\dee{y} =\frac{\pi a^2}{2}, \end{equation*}
comme obtenu dans l’Exemple 3.1.16.

Exemple 3.1.18. Ordre d’intégration.

L’intégrale \(\displaystyle\int_{-1}^2\int_{x^2}^{x+2}\dee{y}\,\dee{x}\) calcule l’aire d’une région du plan \(Oxy\text{.}\) Exprimer la même aire par une intégrale double en inversant l’ordre d’intégration.
L’étape critique lorsqu’on change l’ordre d’intégration est d’illustrer la région dans le plan \(Oxy\text{.}\) À partir de l’intégrale donnée, nous voyons que, sur le domaine d’intégration,
  • \(x\) varie de \(-1\) à \(2\) et que,
  • pour chaque \(x\) fixe, \(y\) varie de la parabole \(y=x^2\) à la droite \(y=x+2\text{.}\)
Cette intégrale itérée correspond aux découpes verticales de la figure de gauche ci-dessous.
Pour inverser l’ordre d’intégration, on doit considérer dans la même région et faire des découpes en rectangles horizontaux, comme illustré dans la figure de droite ci-dessus.
Nous faisons face à un nouveau problème : l’expression de la valeur de \(x\) à l’extrémité gauche d’un rectangle dépend de si son ordonnée \(y\) est inférieure ou supérieure à \(y=1\text{.}\) Plus précisément, nous avons que, sur le domaine d’intégration,
  • \(y\) varie de \(0\) à \(4\) et que,
  • pour chaque \(y\) fixe tel que \(0\leqslant y\leqslant 1\text{,}\) \(x\) varie de \(x=-\sqrt{y}\) à \(x=+\sqrt{y}\) et que,
  • pour chaque \(y\) fixe tel que \(1\leqslant y\leqslant 4\text{,}\) \(x\) varie de \(x=y-2\) à \(x=+\sqrt{y}\text{.}\)
Ainsi,
\begin{gather*} \int_{-1}^2\int_{x^2}^{x+2}\dee{y}\, \dee{x} =\int_0^1\int_{-\sqrt{y}}^{\sqrt{y}}\dee{x}\, \dee{y} + \int_1^4\int_{y-2}^{\sqrt{y}}\dee{x}\, \dee{y}. \end{gather*}
Le but de l’exemple qui suit est de montrer que, malgré que les deux ordres d’intégration doivent donner le même résultat, ils ne sont pas toujours de difficulté semblable.

Exemple 3.1.19.

Évaluer l’intégrale de \(\frac{\sin x}{x}\) sur la région du plan \(Oxy\) qui se situe au-dessus de l’axe des abscisses, à la droite de la droite \(y=x\) et à la gauche de \(x=1\text{.}\)
Solution.
Voici une esquisse du domaine.
Nous allons essayer d’évaluer l’intégrale donnée de deux façons: une fois avec les tranches horizontales (ce qui sera impossible) et une autre fois avec les tranches verticales (la façon la plus simple).
Solution avec des tranches horizontales. Afin de poser l’intégrale en utilisant des tranches horizontales, comme dans la figure ci-dessous, nous observons que, sur notre domaine d’intégration,
  • \(y\) varie de \(0\) à \(1\) et que,
  • pour chaque \(y\) fixe, \(x\) varie de \(x=y\) à \(x=1\text{.}\)
Ainsi, l’intégrale itérée est
\begin{equation*} \int_0^1\int_y^1 \frac{\sin x}{x}\, \dee{x}\ \dee{y}. \end{equation*}
Le problème est que l’intégrande \(\frac{\sin x}{x}\) ne possède pas de primitive élémentaire. Il est impossible d’évaluer \(\int_y^1 \frac{\sin x}{x}\ \dee{x}\) sans faire appel à des méthodes numériques ou à des développements en séries.
Solution avec des tranches verticales.\(\ \ \ \) Afin de poser l’intégrale en utilisant des tranches verticales, comme dans la figure de droite ci-dessus, nous observons que, sur le domaine d’intégration,
  • \(x\) varie de \(0\) à \(1\) et que,
  • pour chaque \(x\) fixe, \(y\) varie de \(0\) à \(y=x\text{.}\)
Donc, l’intégrale itérée est
\begin{equation*} \int_0^1\int_0^x \frac{\sin x}{x}\ \dee{y} \ \dee{x}. \end{equation*}
Cette fois-ci, parce que la variable \(x\) est considérée comme une constante dans l’intégrale intérieure, le calcul de cette intégrale devient simple.
\begin{align*} \int_0^1\int_0^x \frac{\sin x}{x}\ \dee{y} \ \dee{x} &=\int_0^1\frac{\sin x}{x}\left[\int_0^x\dee{y}\right] \dee{x} =\int_0^1 \sin x \ \dee{x} =1-\cos(1) \end{align*}
Voici un dernier exemple afin de rappeler certaines techniques d’intégration de fonctions à une variable.

Exemple 3.1.20.

Calculer le volume du solide situé sous la surface \(z=1-3x^2-2y^2\) et au-dessus du plan \(Oxy\text{.}\)
Solution.
Avant de se précipiter dans l’intégration, il est important de se représenter une esquisse de ce solide et de cette surface. Pour chaque constante \(z_0\text{,}\) la partie de la surface \(z=1-3x^2-2y^2\) qui est projetée dans le plan horizontal \(z=z_0\) est l’ellipse \(3x^2+2y^2=1-z_0\text{.}\) La plus large de ces ellipses est celle qui se trouve dans le plan \(Oxy\) lorsque \(z_0=0\text{,}\) soit l’ellipse d’équation \(3x^2+2y^2=1\text{.}\) Au fur et à mesure que \(z_0\) augmente, l’ellipse se contracte jusqu’au point unique \((0,0,1)\text{,}\) soit lorsque \(z_0=1\text{.}\) Ainsi, la surface est composée d’une superposition d’ellipses, et le solide \(\cV\) est donné par
\begin{equation*} \cV =\Set{(x,y,z)}{3x^2+2y^2\leqslant 1,\ 0\leqslant z\leqslant 1-3x^2-2y^2}. \end{equation*}
Ceci est illustré dans la figure ci-dessous.
La vue d’en haut de la région
\begin{equation*} \cR=\Set{(x,y)}{3x^2+2y^2\leqslant 1} \end{equation*}
est illustrée dans la figure ci-dessous.
En considérant que la relation \(z=1-3x^2-2y^2\) exprimée en fonction de \(x\) est presque identique à celle exprimée en fonction de \(y\) (seulement les coefficients \(2\) et \(3\) sont interchangés), l’utilisation de la méthode des tranches verticales nous amènera à une intégration de même difficulté qu’avec la méthode des tranches horizontales. Nous pouvons donc choisir n’importe laquelle des méthodes; choisissons celle avec les tranches verticales.
L’ellipse dans \(\cR\) coupe l’axe des \(y\) aux points où \(x=0\) et où \(y\) vérifie l’équation \(3(0)^2+2y^2 =1 \text{,}\) donc \(y=\pm\frac{1}{\sqrt{2}}\text{.}\) Ainsi, dans \(\cR\text{,}\) \(-\frac{1}{\sqrt{2}}\leqslant y\leqslant \frac{1}{\sqrt{2}}\) et, pour chaque \(y\) de la région, \(3x^2\leqslant 1-2y^2\text{,}\) ou bien \(-\sqrt{\frac{1-2y^2}{3}}\leqslant x\leqslant \sqrt{\frac{1-2y^2}{3}}\text{.}\) Par conséquent, en utilisant des rectangles verticaux comme dans la figure ci-dessus,
\begin{align*} \text{Volume}(\cV) &=\iint_\cR \big(1-3x^2-2y^2\big)\,\dee{x}\,\dee{y}\\ &=\int_{-\frac{1}{\sqrt{2}}}^{\frac{1}{\sqrt{2}}} \left[\int_{-\sqrt{\frac{1-2y^2}{3}}}^{\sqrt{\frac{1-2y^2}{3}}} \big(1-3x^2-2y^2\big)\ \dee{x}\right] \dee{y}\\ &=4\int_0^{\frac{1}{\sqrt{2}}} \left[\int_0^{\sqrt{\frac{1-2y^2}{3}}} \big(1-3x^2-2y^2\big)\ \dee{x}\right] \dee{y}\\ &=4\int_0^{\frac{1}{\sqrt{2}}} \Big[(1-2y^2)x-x^3\Big]_0^{\sqrt{\frac{1-2y^2}{3}}}\ \dee{y}\\ &=4\int_0^{\frac{1}{\sqrt{2}}} \sqrt{\frac{1-2y^2}{3}}\ \Big[(1-2y^2)-\frac{1-2y^2}{3}\Big]\ \dee{y}\\ &=8\int_0^{\frac{1}{\sqrt{2}}} \left[\frac{1-2y^2}{3}\right]^{3/2}\ \dee{y}. \end{align*}
Nous évaluons cette intégrale en utilisant la substitution trigonométrique et en posant \(2y^2=\sin^2\theta\text{,}\) c’est-à-dire
\begin{equation*} y=\frac{\sin\theta}{\sqrt{2}},\qquad \dee{y}=\frac{\cos\theta}{\sqrt{2}}\,\dee{\theta}. \end{equation*}
On obtient:
\begin{align*} \text{Volume}(\cV) &= 8\int_0^{\frac{\pi}{2}} \overbrace{ \frac{\cos\theta}{\sqrt{2}} }^{\dee{y}} \left[\frac{\cos^2\theta}{3}\right]^{3/2}\ \dee{\theta} =\frac{8}{\sqrt{54}}\int_0^{\frac{\pi}{2}} \cos^4\theta \ \dee{\theta}. \end{align*}
Ensuite, pour intégrer \(\cos^4\theta\text{,}\) nous utilisons l’identité trigonométrique suivante:
\begin{align*} \cos^2\theta &=\frac{\cos(2\theta)+1}{2}\\ \implies \cos^4\theta &= \frac{\big(\cos(2\theta)+1\big)^2}{4} =\frac{\cos^2(2\theta)+2\cos(2\theta)+1}{4}\\ &=\frac{\frac{\cos(4\theta)+1}{2}+2\cos(2\theta)+1}{4}\\ &=\frac{3}{8} +\frac{1}{2}\cos(2\theta) +\frac{1}{8}\cos(4\theta). \end{align*}
Finalement, en considérant que \(\int_0^{\frac{\pi}{2}} \cos(4\theta)\,\dee{\theta} =\int_0^{\frac{\pi}{2}} \cos(2\theta)\,\dee{\theta} =0\text{,}\) on obtient
\begin{gather*} \text{Volume}(\cV) = \frac{8}{\sqrt{54}}\ \frac{3}{8}\ \frac{\pi}{2} =\frac{\pi}{2\sqrt{6}}. \end{gather*}

Sous-section Quelques compléments au sujet de la définition de \(\ \iint_\cR f(x,y)\ \dee{A}\)

En principe, l’intégrale \(\iint_\cR f(x,y)\ \dee{A}\text{,}\)\(\cR\) est une région bornée du plan \(\R^2\text{,}\) est définie comme suit.
  • Subdivisons \(\cR\) en traçant des droites parallèles aux axes de coordonnées.
  • Numérotons les sous-rectangles obtenus de \(1\) à \(n\) (ce sont presque tous des rectangles, sauf ceux qui se trouvent à la frontière de \(\cR\)). Notons que \(n\text{,}\) dans ce contexte, est le nombre total de sous-rectangles, pas le nombre de rectangles verticaux ou horizontaux comme précédemment.
  • Posons \(\De A_i\text{,}\) l’aire du sous-rectangle numéro \(i\text{.}\)
  • Soit \((x_i^*,y_i^*)\) un point d’échantillon sélectionné au hasard dans le sous-rectangle \(i\text{.}\)
  • Considérons la somme \(\sum\limits_{i=1}^n f(x_i^*,y_i^*) \De A_i \text{.}\) Rappelons que la somme tient compte de tous les sous-rectangles dans la région.
Répétons cette procédure en redécoupant \(\cR\) afin d’obtenir des sous-rectangles de plus en plus petits. Si, au fur et à mesure que les dimensions des sous-rectangles tendent vers \(0\text{,}\) cette somme s’approche d’une valeur limite unique (peu importe le choix de la subdivision de la région et des points \((x_i^*,y_i^*)\)), alors nous définissons
\begin{equation*} \iint_\cR f(x,y)\ \dee{A}\ =\lim \sum_{i=1}^n f(x_i^*,y_i^*)\,\De A_i. \end{equation*}
La preuve de ce théorème n’est pas particulièrement difficile, mais elle dépasse le cadre de ce texte. Une conséquence importante de ce théorème est donnée ci-après.
Les hypothèses de ces deux théorèmes peuvent être atténuées légèrement, mais pas trop. Par exemple, il faut que \(f\) « ne soit pas trop discontinue ». En effet, si
\begin{align*} \cR&=\Set{(x,y)}{0\leqslant x\leqslant 1,\ 0\leqslant y\leqslant 1},\\ f(x,y)&=\begin{cases} 1&\text{si } x,y \text{ sont tous deux rationnels} \\ 0&\text{sinon}, \end{cases} \end{align*}
alors l’intégrale \(\iint_\cR f(x,y)\ \dee{A}\) n’existe pas. Effectivement, si tous les \(x_i^*\) et \(y_i^*\) sont choisis rationnels, alors
\begin{equation*} \sum_{i=1}^n f(x_i^*,y_i^*)\,\De A_i =\sum_{i=1}^n \,\De A_i =\text{Aire}(\cR). \end{equation*}
Mais si l’on choisit des nombres irrationnels pour tous les \(x_i^*\) et \(y_i^*\text{,}\) alors
\begin{equation*} \sum_{i=1}^n f(x_i^*,y_i^*)\,\De A_i =\sum_{i=1}^n 0\,\De A_i =0. \end{equation*}
Ainsi, la valeur de la limite de \(\sum\limits_{i=1}^n f(x_i^*,y_i^*)\,\De A_i\) dépend du choix des points d’échantillon \((x_i^*, y_i^*)\text{.}\) L’intégrale \(\iint_\cR f(x,y)\ \dee{A}\) n’existe donc pas.
Voici un exemple encore plus pathologique
 3 
Pour les mathématicien(ne)s, “pathologique” signifie “cool” ou “amusant”.
.

Exemple 3.1.23.

Dans cet exemple, nous atténuons l’une des hypothèses du théorème de Fubini, soit la continuité de la fonction \(f\text{.}\) Nous allons construire un exemple dans lequel les deux intégrales des deux ordres d’intégration existent, mais ne sont pas égales. Plus précisément, choisissons \(\cR=\Set{(x,y)}{0\leqslant x\leqslant 1,\ 0\leqslant y\leqslant 1}\) et définissons une fonction continue sur \(\cR\text{,}\) sauf à l’origine.
Pour commencer, soit \(\d_1,\d_2,\d_3,\ \dots\text{,}\) une suite de nombres réels tels que
\begin{equation*} 1=\d_1 \gt \d_2 \gt \d_3\gt\ \cdots\ \gt\d_n\rightarrow 0. \end{equation*}
Par exemple, \(\d_n=\frac{1}{n}\) ou \(\d_n=\frac{1}{2^{n-1}}\) sont des choix valables. Pour chaque entier positif \(n\text{,}\) posons l’intervalle \(I_n=]\d_{n+1},\d_n]=\Set{t}{\d_{n+1} \lt t\leqslant \d_n}\) et posons \(g_n(t)\) une fonction non négative et continue vérifiant
  • \(g_n(t)=0\) si \(t \not\in I_n\text{,}\) et
  • \(\int_{I_n}g(t)\,dt=1\text{.}\)
Il existe plusieurs exemples d’une telle fonction, par exemple
\begin{equation*} g_n(t)=\left(\frac{2}{\d_n-\d_{n+1}}\right)^2\begin{cases} \d_n-t& \text{si } \half(\d_{n+1}+\d_n)\leqslant t\leqslant \d_n\\ t-\d_{n+1}& \text{si }\d_{n+1}\leqslant t\leqslant \half(\d_{n+1}+\d_n)\\ 0& \text{sinon} \end{cases}. \end{equation*}
Voici un résumé de ce qui a été fait jusqu’à maintenant.
  • On a subdivisé l’intervalle \(0 \lt x\leqslant 1\) en une infinité de sous-intervalles \(I_n\text{.}\) À mesure que \(n\) augmente, les sous-intervalles deviennent de plus en plus petits, et \(I_n\) s’approche de plus en plus près de zéro.
  • On a défini, pour chaque \(n\text{,}\) une fonction non négative et continue \(g_n\) qui s’annule lorsque \(t\) est à l’extérieur de \(I_n\) et dont l’intégrale sur \(I_n\) vaut 1.
Maintenant, définissons l’intégrande \(f(x,y)\) en fonction de ces sous-intervalles \(I_n\) et de ces fonctions \(g_n\text{.}\)
\begin{equation*} f(x,y)=\begin{cases} 0& \text{si } x=0 \\ 0& \text{si } y=0 \\ g_m(x)g_n(y)& \text{si } x\in I_m,\ y\in I_n \text{ et } m=n\\ -g_m(x)g_n(y)& \text{si } x\in I_m,\ y\in I_n \text{ et } m=n+1\\ 0& \text{sinon} \end{cases} \end{equation*}
On peut imaginer \(]0,1]\times]0,1]\) comme l’union de plusieurs petits rectangles de dimensions \(I_m\times I_n\text{,}\) telle qu’illustrée dans la figure ci-dessous. Sur la majorité de ces rectangles, \(f(x,y)\) vaut \(0\text{.}\) Les exceptions sont les rectangles gris foncé de dimensions \(I_n\times I_n\) sur la “diagonale” de la figure ainsi que les rectangles gris pâle de dimensions \(I_{n+1}\times I_n\) à la gauche de la “diagonale”.
Sur chaque rectangle gris foncé, \(f(x,y)\geqslant 0\text{,}\) et le graphe de \(f(x,y)\) est le graphe de \(g_n(x)g_n(y)\) qui ressemble à une pyramide. Sur chaque rectangle gris pâle, \(f(x,y)\leqslant 0\text{,}\) et le graphe de \(f(x,y)\) est celui de \(-g_{n+1}(x)g_n(y)\) qui ressemble plutôt à une pyramide inversée.
Maintenant, fixons un \(y\) tel que \(0\leqslant y\leqslant 1\) et calculons \(\int_0^1 f(x,y)\ \dee{x}\text{.}\) Ainsi, nous intégrons \(f\) le long d’une droite parallèle à l’axe des abscisses. Si \(y=0\text{,}\) alors \(f(x,y)=0\) pour tout \(x\text{,}\) de sorte que \(\int_0^1 f(x,y)\,\dee{x} = 0\text{.}\) Si \(0 \lt y\leqslant 1\text{,}\) alors il existe exactement un entier positif \(n\) tel que \(y\in I_n\) et \(f(x,y)=0\text{,}\) sauf lorsque \(x\) se trouve dans \(I_n\) ou dans \(I_{n+1}\text{.}\) Ainsi, pour \(y\in I_n\text{,}\)
\begin{align*} \int_0^1 f(x,y)\,\dee{x}&=\sum_{m=n,n+1}\int_{I_m} f(x,y)\,\dee{x}\\ &=\int_{I_n} g_n(x)g_n(y)\,\dee{x}-\int_{I_{n+1}} g_{n+1}(x)g_n(y)\,\dee{x}\\ &=g_n(y)\int_{I_n} g_n(x)\,\dee{x}-g_n(y)\int_{I_{n+1}} g_{n+1}(x)\,\dee{x}\\ &=g_n(y)-g_n(y)=0. \end{align*}
Ici, nous avons utilisé à deux reprises le fait que \(\int_{I_m}g(t)\,dt=1\) pour tout \(m\text{.}\) Ainsi, \(\int_0^1 f(x,y)\,\dee{x}=0\) pour tout \(y\text{,}\) et donc \(\int_0^1\Big[\int_0^1 f(x,y)\ \dee{x}\, \dee{y}\Big]=0\text{.}\)
Finalement, nous fixons un \(x\) tel que \(0\leqslant x\leqslant 1\) et nous évaluons \(\int_0^1 f(x,y)\ \dee{y}\text{.}\) Cette fois, nous intégrons \(f\) le long d’une droite parallèle à l’axe des ordonnées. Si \(x=0\text{,}\) alors \(f(x,y)=0\) pour tout \(y\text{,}\) donc \(\int_0^1 f(x,y)\,\dee{y} = 0\text{.}\) Si \(0 \lt x\leqslant 1\text{,}\) alors il existe un unique entier positif \(m\) tel que \(x\in I_m\text{.}\) Si \(m\geqslant 2\text{,}\) alors \(f(x,y)\) est nulle, sauf lorsque \(y\) est dans \(I_m\) ou dans \(I_{m-1}\text{.}\) Mais si \(m=1\text{,}\) alors \(f(x,y)\) est nulle, sauf pour \(y\in I_1\text{.}\) (Voir la figure ci-dessus de nouveau!) Ainsi, pour \(x\in I_m\text{,}\) avec \(m\geqslant 2\text{,}\)
\begin{align*} \int_0^1 f(x,y)\,\dee{y}&=\sum_{n=m,m-1}\int_{I_n} f(x,y)\,\dee{y}\\ &=\int_{I_m} g_m(x)g_m(y)\,\dee{y}-\int_{I_{m-1}} g_{m}(x)g_{m-1}(y)\,\dee{y}\\ &=g_m(x)\int_{I_m} g_m(y)\,\dee{y}-g_m(x)\int_{I_{m-1}} g_{m-1}(y)\,\dee{y}\\ &=g_m(x)-g_m(x)=0. \end{align*}
Mais, pour \(x\in I_1\text{,}\)
\begin{align*} \int_0^1 f(x,y)\,\dee{y}&=\int_{I_1} f(x,y)\,\dee{y} =\int_{I_1} g_1(x)g_1(y)\,\dee{y} =g_1(x)\int_{I_1} g_1(y)\,\dee{y}\\ &=g_1(x). \end{align*}
Ainsi,
\begin{equation*} \int_0^1 f(x,y)\,\dee{y}=\begin{cases} 0&\text{si } x\leqslant \d_2 \\ g_1(x)&\text{si } x\in I_1 \end{cases}, \end{equation*}
de sorte que
\begin{equation*} \int_0^1 \bigg[\int_0^1 f(x,y) \ \dee{y} \bigg] \dee{x} =\int_{I_1}g_1(x)\,\dee{x}=1. \end{equation*}
La conclusion est que, pour la fonction \(f(x,y)\) définie ci-dessus, qui existe pour tout \(0\leqslant x\leqslant 1\text{,}\) \(0\leqslant y\leqslant 1\) et qui est continue partout, sauf en \((0,0)\text{,}\) nous avons
\begin{equation*} \int_0^1\int_0^1 f(x,y)\ \dee{x}\, \dee{y}=0 \qquad\text{et}\qquad \int_0^1\int_0^1 f(x,y)\ \dee{y}\, \dee{x}=1. \end{equation*}

Exercices Exercices

Pour se pratiquer.

1.
Évaluez les intégrales doubles sans calculer en termes d’intégrales itérées. Plutôt, interprétez l’intégrale comme étant, par exemple, une aire ou un volume.
  1. \(\displaystyle \ds\int_{-1}^3\int_{-4}^1 \dee{y}\,\dee{x}\)
  2. \(\displaystyle \ds\int_0^2\int_0^{\sqrt{4-y^2}} \dee{x}\,\dee{y}\)
  3. \(\displaystyle \ds\int_{-3}^3\int_0^{\sqrt{9-y^2}}\sqrt{9-x^2-y^2}\ \dee{x}\,\dee{y}\)
Réponse.
(a) \(20\)
(b) \(\pi\)
(c) \(9\pi\)
2.
Soit \(f(x,y)= 12 x^2y^3\text{.}\) Évaluez
  1. \(\displaystyle \ds\int_0^3 f(x,y)\,\dee{x}\)
  2. \(\displaystyle \ds\int_0^2 f(x,y)\,\dee{y}\)
  3. \(\displaystyle \ds\int_0^2\int_0^3 f(x,y)\,\dee{x}\,\dee{y}\)
  4. \(\displaystyle \ds\int_0^3\int_0^2 f(x,y)\,\dee{y}\,\dee{x}\)
  5. \(\displaystyle \ds\int_0^3\int_0^2 f(x,y)\,\dee{x}\,\dee{y}\)
Réponse.
(a) \(108 y^3\)
(b) \(48x^2\)
(c), (d) \(432\)
(e) \(648\)
3.
Pour chacun des cas suivants, évaluez l’intégrale double donnée en considérant les intégrales comme des intégrales itérées.
  1. \(\displaystyle\iint_R (x^2+y^2)\,\dee{A}\text{,}\)\(R\) est le rectangle \(0\leqslant x\leqslant a,\ 0\leqslant y\leqslant b\text{,}\)\(a \gt 0\) et \(b \gt 0\text{.}\)
  2. \(\displaystyle\iint_T (x-3y)\,\dee{A}\text{,}\)\(T\) est le triangle avec sommets \((0,0),\ (a,0),\ (0,b)\text{.}\)
  3. \(\displaystyle\iint_R xy^2\,\dee{A}\text{,}\)\(R\) est la région finie dans le premier cadran bornée par la courbe \(y=x^2\) et \(x=y^2\text{.}\)
  4. \(\displaystyle\iint_D x\cos y\,\dee{A}\text{,}\)\(D\) est la région finie dans le premier cadran bornée par les axes de coordonnées et la courbe \(y=1-x^2\text{.}\)
  5. \(\displaystyle\iint_R {x\over y}e^y\,\dee{A}\text{,}\)\(R\) est la région \(0\leqslant x\leqslant 1,\ x^2\leqslant y\leqslant x\text{.}\)
  6. \(\displaystyle\iint_T {xy\over 1+x^4}\,\dee{A}\text{,}\)\(T\) est le triangle avec les sommets \((0,0),\ (0,1),\ (1,1)\text{.}\)
Réponse.
(a) \(\frac{1}{3}\big(a^3b+ab^3\big)\)
(b) \(\frac{a^2b}{6}-\frac{ab^2}{2}\)
(c) \(\frac{3}{56}\)
(d) \(\frac{1}{2}(1-\cos 1)\)
(e) \(\frac{1}{2}(e-2)\)
(f) \(\frac{1}{4}\left(\frac{\pi}{4}-\frac{1}{2}\ln 2\right)\)
4.
Combinez la somme des deux intégrales doubles itérées
\begin{equation*} \int_{y=0}^{y=1}\int_{x=0}^{x=y} f(x,y)\ \dee{x}\,\dee{y} +\int_{y=1}^{y=2}\int_{x=0}^{x=2-y} f(x,y)\ \dee{x}\,\dee{y} \end{equation*}
en une seule intégrale double itérée avec l’ordre d’intégration inversé.
Réponse.
\(\int_{x=0}^{x=1}\int_{y=x}^{y=2-x} f(x,y)\ \dee{y}\,\dee{x}\)
5.
Considérez l’intégrale
\begin{equation*} \int_0^1\int_x^1 e^{x/y}\ \dee{y}\ \dee{x} \end{equation*}
  1. Esquissez le domaine d’intégration.
  2. Évaluez l’intégrale en inversant l’ordre d’intégration.
Réponse.
(a)
(b) \(\half(e-1)\)

Pour mieux comprendre.

6.
Pour chacune des intégrales suivantes (i) esquissez le région d’intégration, (ii) écrivez une intégrale double équivalente avec un ordre d’intégration inversé et (iii) évaluez les deux intégrales doubles.
  1. \(\displaystyle \displaystyle\int_0^2\int_1^{e^x}\dee{y}\dee{x}\)
  2. \(\displaystyle \displaystyle\int_0^{\sqrt{2}} \int_{-\sqrt{4-2y^2}}^{\sqrt{4-2y^2}}\dee{x}\dee{y}\)
  3. \(\displaystyle \displaystyle\int_{-2}^1 \int_{x^2+4x}^{3x+2}\dee{y}\dee{x}\)
Réponse.
(a) \(e^2-3\)
(b) \(\frac{8}{3}\)
(c) \(\frac{9}{2}\)
7.
L’integrale \(I\) est définie par
\begin{equation*} I =\iint_R f(x,y)\ \dee{A} = \int_1^{\sqrt{2}} \int_{1/y}^{\sqrt{y}} f(x,y)\,\dee{x}\,\dee{y} +\int_{\sqrt{2}}^4 \int_{y/2}^{\sqrt{y}} f(x,y)\,\dee{x}\,\dee{y} \end{equation*}
  1. Esquissez la région \(R\text{.}\)
  2. Réécrire l’intégrale \(I\) en inversant l’ordre d’intégration.
  3. Calculez l’intégrale \(I\) lorsque \(f(x,y)= x/y\text{.}\)
Réponse.
(a) La région \(R\) est la région ombragée dans la figure
(b) \(I= \int_{1/\sqrt{2}}^1 \int_{1/x}^{2x} f(x,y)\,\dee{y}\,\dee{x} +\int_1^2 \int_{x^2}^{2x} f(x,y)\,\dee{y}\,\dee{x}\)
(c) \(\frac{1}{2}\)
8.
Une région \(E\) dans le plan \(xy\) a la propriété que pour toutes fonctions continues f
\begin{equation*} \iint_E f(x,y)\,\dee{A} = \int_{x=-1}^{x=3}\left[\int_{y=x^2}^{y=2x+3} f(x,y) \dee{y}\right] \dee{x} \end{equation*}
  1. Calculez \(\iint_E x\,\dee{A}\text{.}\)
  2. Esquissez la région E.
  3. Poser \(\iint_E x\,\dee{A}\) comme étant l’intégrale ou la somme d’intégrales dans l’ordre opposé.
Réponse.
(a) \(\frac{32}{3}\)
(b)
(c) \(I = \int_0^1\dee{y} \int_{-\sqrt{y}}^{\sqrt{y}}\dee{x}\ x +\int_1^9\dee{y} \int_{(y-3)/2}^{\sqrt{y}}\dee{x}\ x\)
9.
Calculez l’intégrale:
\begin{equation*} \iint_D \sin(y^2)\ \dee{A} \end{equation*}
\(D\) est la région bornée par \(x + y = 0\text{,}\) \(2x - y = 0\text{,}\) et \(y = 4\text{.}\)
Réponse.
\(\frac{3}{4}\big[1-\cos(16)\big]\)
10.
Considérez l’intégrale
\begin{equation*} I = \int_0^1 \int_{\sqrt{y}}^1 \frac{\sin(\pi x^2)}{x}\ \dee{x}\,\dee{y} \end{equation*}
  1. Esquissez la région d’intégration.
  2. Évaluez \(I\text{.}\)
Réponse.
(a)
(b) \(\frac{1}{\pi}\)
11.
Soit \(I\) l’intégrale double de la fonction \(f(x,y) = y^2 \sin xy\) sur le triangle déterminé par les sommets \((0, 0)\text{,}\) \((0, 1)\) et \((1, 1)\) dans le plan \(xy\text{.}\)
  1. Écrire \(I\) comme une intégrale itérée de deux façons différentes.
  2. Évaluez \(I\text{.}\)
Réponse.
(a) \(I=\int_0^1\dee{x}\int_x^1\dee{y}\ y^2\sin xy =\int_0^1\dee{y}\int_0^y\dee{x}\ y^2\sin xy\)
(b) \(\frac{1-\sin 1}{2}\)
12.
Trouvez le volume \((V)\) du solide borné au-dessus par la surface
\begin{equation*} z = f (x,y) = e^{-x^2}, \end{equation*}
et en-dessous par le plan \(z = 0\) et sur le triangle dans le plan \(xy\) formé par les droites \(x = 1\text{,}\) \(y = 0\) et \(y = x\text{.}\)
Réponse.
\(\frac{1-e^{-1}}{2}\)
13.
Considérez l’intégrale \(\displaystyle I=\int_0^1 \int_y^{2-y}\frac{y}{x} \ \dee{x}\,\dee{y}\text{.}\)
  1. Esquissez la région d’intégration.
  2. Interchanger l’ordre d’intégration.
  3. Évaluez \(I\text{.}\)
Réponse.
(a)
(b) \(I=\int_0^1\dee{x}\int_0^x\dee{y}\ \frac{y}{x} +\int_1^2\dee{x}\int_0^{2-x}\dee{y}\ \frac{y}{x}\)
(c) \(2\ln 2 -1\)
14.
Pour l’intégrale
\begin{equation*} I = \int_0^1 \int_{\sqrt{x}}^1 \sqrt{1+y^3}\ \dee{y}\,\dee{x} \end{equation*}
  1. Esquissez la région d’intégration.
  2. Évaluez \(I\text{.}\)
Réponse.
(a)
(b) \(\frac{2\big(2\sqrt{2}-1\big)}{9}\)
15.
  1. \(D\) est la région bornée par la parabole \(y^2 = x\) et la droite \(y = x - 2\text{.}\) Esquissez \(D\) et évaluez \(J\)
    \begin{equation*} J = \iint_D 3y\ \dee{A} \end{equation*}
  2. Esquissez la région d’intégration et ensuite évaluez l’intégrale \(I\) :
    \begin{equation*} I = \int_0^4 \int_{\frac{1}{2}\sqrt{x}}^1 e^{y^3}\ \dee{y}\,\dee{x} \end{equation*}
Réponse.
(a) \(J=\frac{27}{4}\)
(b) \(I=\frac{4}{3}\big[e-1\big]\)
16.
Considérez l’intégrale itérée
\begin{equation*} \int_{-4}^0\int_{\sqrt{-y}}^2 \cos(x^3)\,\dee{x}\,\dee{y} \end{equation*}
  1. Dessinez la région d’intégration.
  2. Évaluez l’intégrale.
Réponse.
(a)
(b) \(\frac{\sin(8)}{3}\)
17.
  1. Combinez la somme des intégrales itérées
    \begin{equation*} I = \int_0^1\int_{-\sqrt{y}}^{\sqrt{y}} f(x,y)\ \dee{x}\,\dee{y} + \int_1^4\int_{y-2}^{\sqrt{y}} f(x,y)\ \dee{x}\,\dee{y} \end{equation*}
    en une seule intégrale itérée avec l’ordre d’intégration inversé.
  2. Évaluez \(I\) si \(f(x,y)=\frac{e^x}{2-x}\text{.}\)
Réponse.
(a) \(I = \int_{-1}^2\int_{x^2}^{x+2} f(x,y)\ \dee{y}\,\dee{x}\)
(b) \(2e^2 + \frac{1}{e}\)
18.
Soit
\begin{equation*} I=\int_0^4\int_{\sqrt{y}}^{\sqrt{8-y}} f(x,y)\,\dee{x}\,\dee{y} \end{equation*}
  1. Esquissez le domaine d’intégration.
  2. Inverser l’ordre d’intégration.
  3. Évaluez l’integrale pour \(f(x,y)=\frac{1}{(1+y)^2}\text{.}\)
Réponse.
(a)
(b) \(\int_0^2\int_0^{x^2} f(x,y)\,\dee{y}\,\dee{x} +\int_2^{\sqrt{8}}\int_0^{8-x^2} f(x,y)\,\dee{y}\,\dee{x}\)
(c) \(\sqrt{8}-\arctan 2 -\frac{1}{6}\left[\ln\frac{3+\sqrt{8}}{3-\sqrt{8}} -\ln 5\right]\)
19.
Évaluez
\begin{equation*} \int_{-1}^0 \int_{-2}^{2x} e^{y^2}\ \dee{y}\,\dee{x} \end{equation*}
Réponse.
\(\frac{1}{4}\big[e^4-1\big]\)
20.
Soit
\begin{equation*} I = \int_0^2 \int_0^x f(x,y)\ \dee{y}\,\dee{x} + \int_2^6 \int_0^{\sqrt{6-x}} f(x,y)\ \dee{y}\,\dee{x} \end{equation*}
Exprimez \(I\) comme une intégrale où on intègre d’abord en respect à \(x\text{.}\)
Réponse.
\(I = \int_0^2 \int_y^{6-y^2} f(x,y)\ \dee{x}\,\dee{y}\)
21.
Considérez le domaine \(D\) sur l’axe des \(x\) sous la parabole \(y = 1-x^2\) dans le plan \(xy\text{.}\)
  1. Esquissez \(D\text{.}\)
  2. Exprimez
    \begin{equation*} \iint_D f(x,y)\ \dee{A} \end{equation*}
    comme une intégrale itérée correspondante à l’ordre \(\dee{x}\,\dee{y}\text{.}\) Exprimez ensuite cette intégrale comme une intégrale itérée correspondante à l’ordre \(\dee{y}\,\dee{x}\text{.}\)
  3. Calculez l’intégrale dans le cas où \(f(x,y) = e^{x-(x^3/3)}\text{.}\)
Réponse.
(a)
(b) \(\int_0^1 \int_{-\sqrt{1-y}}^{\sqrt{1-y}} f(x,y)\ \dee{x}\,\dee{y}\text{,}\) \(\int_{-1}^1 \int_0^{1-x^2} f(x,y)\ \dee{y}\,\dee{x}\)
(c) \(e^{2/3}-e^{-2/3}\)
22.
Soit \(I=\int_0^1 \int_{x^2}^1 x^3\ \sin(y^3)\ \dee{y}\ \dee{x}\text{.}\)
  1. Esquissez la région d’intégration dans le plan \(xy\text{.}\) Étiquettez votre esquisse suffisamment bien de façon à ce qu’une personne puisse déterminer les bornes de l’intégrale double.
  2. Évaluez \(I\text{.}\)
Réponse.
(a)
(b) \(\frac{1-\cos(1)}{12}\)
23.
Considérez le solide sous la surface \(z=6-xy\text{,}\) borné par les cinq plans \(x=0\text{,}\) \(x=3\text{,}\) \(y=0\text{,}\) \(y=3\text{,}\) \(z=0\text{.}\) Notez qu’aucune partie de ce solide se trouve sous le plan \(xy\text{.}\)
  1. Esquissez la base du solide dans le plan \(xy\text{.}\) Notez que ce n’est pas un carrée!
  2. Calculez le volume du solide.
Réponse.
(a)
(b) \(27+18\ln\frac{3}{2}\approx 34.30\)
24.
Évaluez l’intégrale suivante:
\begin{equation*} \int_{-2}^2\int_{x^2}^4\cos\big(y^{3/2}\big)\ \dee{y}\,\dee{x} \end{equation*}
Réponse.
\(\frac{4}{3}\sin 8\approx 1.319\)
25.
Évaluez l’intégrale suivante:
\begin{equation*} \int_0^9\int_{\sqrt{y}}^3\sin(\pi x^3)\ \dee{x}\dee{y} \end{equation*}
Réponse.
\(\frac{2}{3\pi}\approx 0.212\)
26.
L’intégrale itérée
\begin{equation*} I=\int_0^1\bigg[\int_{-\sqrt{x}}^{\sqrt{x}} \sin\big(y^3-3y\big)\,\dee{y}\bigg] \ \dee{x} \end{equation*}
est égale à \(\iint_R\sin\big(y^3-3y)\ dA\) pour une région \(R\) dans le plan \(xy\text{.}\)
  1. Esquissez la région \(R\text{.}\)
  2. Écrivez l’intégrale \(I\) avec l’ordre d’intégration inversé, et avec les bonnes bornes d’intégrations.
  3. Trouvez \(I\text{.}\)
Réponse.
(a)
(b) \(\ds\int_{-1}^1\bigg[\int_{y^2}^1 \sin\big(y^3-3y\big)\,\dee{x}\bigg]\ \dee{y}\)
(c) \(0\)
27.
Trouvez l’intégrale double de la fonction \(f(x, y) = xy\) sur la région bornée par \(y = x - 1\) et \(y^2 = 2x + 6\text{.}\)
Réponse.
\(36\)
28.
Trouvez le volume du solide à l’intérieur du cylindre \(x^2+2y^2=8\text{,}\) au-dessus du plan \(z=y-4\) et sous le plan \(z=8-x\text{.}\)
Réponse.
\(48\sqrt{2}\,\pi\)

Pour aller plus loin.

29.
Considérez le volume au-dessus du plan \(xy\) qui se trouve à l’intérieur du cylindre circulaire \(x^2+y^2=2y\) et sous la surface \(z=8+2xy\text{.}\)
  1. Exprimez ce volume comme une intégrale double \(I\text{,}\) en posant clairement le domaine sur lequel \(I\) est considérée.
  2. Exprimez en coordonnées cartésiennes, l’intégrale double \(I\) comme une intégrale itérée de deux façons différentes, en indiquant clairement les bornes d’intégrations dans chaque cas.
  3. Quel est ce volume?