title: Intégrale multiple
text: # Intégrale multiple
  
  En analyse mathématique, l'intégrale multiple est une forme d'intégrale qui s'applique aux fonctions de plusieurs variables réelles. Les deux principaux outils de calcul sont le changement de variables et le théorème de Fubini. Ce dernier permet de ramener de proche en proche un calcul d'intégrale multiple à des calculs d'intégrales simples, et d'interpréter le « volume » d'un domaine « simple » de dimension n (ou son hypervolume si n > 3) comme l'intégrale d'une fonction de n – 1 variables (Fig. 2), de même que l'intégrale définie d'une fonction continue positive d'une variable est égale à « l'aire sous la courbe » associée (Fig. 1).
  
  ## Définition
  
  Soit f : D → ℝ une fonction bornée sur une partie bornée D de $ℝ^n.$ Pour définir (au sens de Riemann) l'intégrabilité de f et, le cas échéant, la valeur de son intégrale, notée
  
  $\int \!\dots \!\int _{D}f(x_{1},\dots ,x_{n})\;\mathrm {d} x_{1}\dots \mathrm {d} x_{n}$
  
  ou plus simplement $\int f$, on se ramène d'abord au cas où D est un pavé fermé P = $I_1$ × … × $I_n,$ en choisissant un tel P qui contient D et en prolongeant f par 0 (on démontre après coup que les définitions suivantes sont indépendantes du choix de P). Les définitions sont ensuite identiques, mutatis mutandis, à celles dans le cas n = 1 : à l'aide des notions de volume et subdivision de P et de pas, marquage et raffinement d'une telle subdivision, on définit l'intégrabilité et l'intégrale soit en termes d'intégrales inférieure et supérieure, soit en termes de sommes de Darboux inférieure et supérieure, soit en termes de sommes de Riemann, et l'on montre que ces trois définitions sont équivalentes. La première, par exemple, s'écrit :
  
  $\int f=\int _{*}f=\int ^{*}f$
  
  (c'est-à-dire que f est intégrable si et seulement si ses intégrales inférieure, $\int _{*}f$, et supérieure, $\int ^{*}f$, sont égales, et l'intégrale de f est alors par définition cette valeur commune), et la troisième s'écrit :
  
  $\int f=\lim _{\delta \to 0}\sum _{1\leq k\leq N}f(\xi _{k})\operatorname {vol} (P_{k}),$
  
  où chaque $(P_k,$ $ξ_k)_{1≤k≤N}$ est une subdivision marquée de P de pas δ.
  
  Dans les applications pratiques, comme l'ingénierie ou la physique appliquée, on rencontre quasi exclusivement des intégrales simples, doubles ou triples. L'intégrale, sur un domaine D de dimension n = 1, 2 ou 3, de la fonction constante 1, est respectivement la longueur, l'aire ou le volume de D.
  
  ## Propriétés
  
  ### Critère de Lebesgue
  
  Le théorème de Lebesgue-Vitali étend aux fonctions de plusieurs variables le critère de Lebesgue pour l'intégrabilité de Riemann : une fonction bornée sur une partie bornée de $ℝ^n$ est Riemann-intégrable si et seulement si l'ensemble des points de discontinuité de son prolongement par 0 à $ℝ^n$ est négligeable. En particulier :
  
  Théorème — Toute fonction continue et bornée sur une partie cubable de $ℝ^n$ est Riemann-intégrable.
  
  En effet, une partie de $ℝ^n$ est cubable (si et) seulement si elle est bornée et de frontière négligeable. Cette hypothèse de régularité du domaine est indispensable : une fonction continue sur un compact peut ne pas être intégrable (prendre, sur un compact non cubable, une fonction constante non nulle). Mais en pratique, tous les domaines considérés seront cubables, comme intersections finies de domaines « simples », au sens défini ci-dessous.
  
  ### Propriétés élémentaires
  
  Les autres propriétés des intégrales simples — linéarité, monotonie, etc. se généralisent également. On en déduit par exemple que si D est un compact connexe de $ℝ^n,$ cubable mais non négligeable, et f : D → ℝ une fonction continue, alors il existe un point m de D en lequel la « valeur moyenne » ou « moyenne intégrale » de f sur D est atteinte, c'est-à-dire tel que
  
      $f(m)={\frac {1}{\operatorname {vol} (D)}}\int _{D}f.$
  
  Cette propriété généralise le théorème de la moyenne en dimension 1, à cela près qu'en dimension supérieure, on n'est plus assuré de trouver un m intérieur à D.
  
  ### Théorème de Fubini
  
  Le théorème de Fubini pour les intégrales de Riemann est énoncé ci-dessous pour une intégrale double pour simplifier les notations, mais il se généralise en remplaçant ℝ×ℝ par $ℝ^p×ℝ^q.$ De plus, on peut évidemment intervertir les rôles des deux variables (mono- ou multidimensionnelles) x et y.
  
  Théorème — Pour toute application f : [a, b]×[c, d] → ℝ :
  
  $\int _{*}f\leq \int _{*}\left(\int _{*}f(x,y)\,{\rm {d}}y\right)\,{\rm {d}}x\leq \int ^{*}\left(\int ^{*}f(x,y)\,{\rm {d}}y\right)\,{\rm {d}}x\leq \int ^{*}f.$
  
     Corollaire — Lorsque f est Riemann-intégrable, les deux fonctions (intégrales inférieure et supérieure)
     $[a,b]\to \mathbb {R} ,\quad x\mapsto \int _{*}f(x,y)\,{\rm {d}}y\quad {\rm {et}}\quad x\mapsto \int ^{*}f(x,y)\,{\rm {d}}y$
     sont Riemann-intégrables et de même intégrale que f.
  
  L'énoncé de la conclusion se simplifie lorsque de plus, pour tout x ∈ [a, b], la fonction y ↦ f(x, y) est intégrable (sur [c, d]) : alors, son intégrale (fonction de x ∈ [a, b]) est intégrable et
  
  $\iint _{[a,b]\times [c,d]}f(x,y)\,{\rm {d}}x{\rm {d}}y=\int _{a}^{b}\left(\int _{c}^{d}f(x,y)\,{\rm {d}}y\right)\,{\rm {d}}x.$
  
  ### Changement de variables
  
  Soient U et V deux ouverts de $ℝ^n,$ Φ : U → V un $C^1-difféomorphisme$ et f une application de V dans ℝ.
  
  Théorème — Si U et V sont bornés et si f et (f∘Φ) × |det $J_Φ|$ sont Riemann-intégrables, alors
  
  $\int _{V}f=\int _{U}(f\circ \Phi )\times |{\rm {det}}J_{\Phi }|.$
  
  Ces hypothèses sont vérifiées dès que U et V sont cubables, f est continue bornée, et det $J_Φ$ est borné.
  
  ## Méthodes d'intégration
  
  ### Examen direct
  
  **Linéarité**
  
  Par linéarité, pour toute constante c, $\int (cf)=c\int f.$ En particulier, l'intégrale d'une fonction constante égale à c, sur un domaine cubable D de $ℝ^n,$ est le produit de c par le volume de D.
  
  De même, l'intégrale d'une somme de fonctions intégrables peut être transformée en la somme de leurs intégrales : $\int _{D}(f_{1}+f_{2})=\int _{D}f_{1}+\int _{D}f_{2}.$
  
  En particulier, on peut partitionner le domaine d'intégration en un nombre fini de sous-domaines — ce qui revient à écrire f comme la somme des produits de f par les indicatrices de ces sous-domaines — puis appliquer la règle précédente. C'est ce qui remplace, en dimensions supérieures, la relation de Chasles propre à la dimension 1 : $\int _{D_{1}\sqcup D_{2}}f=\int _{D_{1}}f+\int _{D_{2}}f.$
  
  **Exploitation des symétries**
  
  Si la fonction est impaire par rapport à une variable $x_i$ et si le domaine est invariant par la symétrie par rapport à l'hyperplan $x_i$ = 0, l'intégrale s'annule.
  
  Soit la fonction $f(x,y,z)=x\exp(y^{2}+z^{2})$ avec pour domaine d'intégration la boule unité de $ℝ^3.$ Le domaine présente une symétrie par rapport aux trois plans de coordonnées, mais comme la fonction est impaire par rapport à la variable x, il suffit de la symétrie par rapport au plan yOz pour annuler toute l'intégrale.
  
  ### Réduction à des intégrales simples
  
  La réduction en intégrales simples utilise le concept de domaine simple, de façon à exprimer l'intégrale en une composition d'intégrales simples. L'intégration est effectuée de l'intérieur vers l'extérieur, chaque fois par rapport à une seule variable en considérant les autres constantes, de la même façon que pour le calcul de dérivées partielles.
  
  **Domaines simples**
  
  Si A est une partie cubable de $ℝ^{n–1}$ et si deux fonctions g et h, définies sur A, sont Riemann-intégrables — par exemple si elles sont bornées et continues (éventuellement par morceaux) — alors la partie suivante de $ℝ^n,$ appelée un « domaine simple orhogonal à l'hyperplan t = 0 », est cubable (Fig. 3) :
  
  $D=\{(x,t)\in \mathbb {R} ^{n-1}\times \mathbb {R} \mid x\in A\quad {\rm {et}}\quad g(x)<t<h(x)\}.$
  
  On peut bien sûr, dans cette définition et ce théorème, permuter les composantes des n-uplets de $ℝ^n,$ pour choisir laquelle joue le rôle de t.
  
  On peut également, dans la définition de D, remplacer un ou deux < par des ≤, car le graphe d'une fonction intégrable sur A est Jordan-négligeable.
  
  **Théorème**
  
  Si f est une fonction continue bornée sur ce domaine simple D, les hypothèses du critère de Lebesgue sont vérifiées, et le théorème de Fubini devient :
  
  $\int _{D}f=\int _{A}\left(\int _{g(x)}^{h(x)}f(x,t)\,{\rm {d}}t\right)\,{\rm {d}}x.$
  
  **Exemples**
  
  **Réduction d'une intégrale double**
  
  La région $D=\{(x,y)\mid x\geq 0,x^{2}\leq y\leq 1\}$ (Fig. 4) se présente comme le domaine simple (orthogonal à l'axe y = 0), associé aux deux fonctions g(x) = $x^2$ et h(x) = 1, sur l'intervalle A = [0, 1] (on trouve l'extrémité 1, abscisse du point d'intersection des courbes de g et h, en calculant la solution positive de l'équation g(x) = h(x)).
  
  On se propose d'intégrer sur D la fonction f(x, y) = x + y. Appliquons le théorème :
  
      $\iint _{D}(x+y)\,\mathrm {d} x\mathrm {d} y=\int _{0}^{1}\left(\int _{x^{2}}^{1}(x+y)\,\mathrm {d} y\right)\,\mathrm {d} x.$
  
  On calcule donc d'abord, pour x fixé :
  
      $\int _{x^{2}}^{1}(x+y)\,\mathrm {d} y=\left[xy+{\frac {y^{2}}{2}}\right]_{x^{2}}^{1}=x+{\frac {1}{2}}-x^{3}-{\frac {x^{4}}{2}}$
  
  puis l'intégrale définie :
  
      $\int _{0}^{1}\left(x+{\frac {1}{2}}-x^{3}-{\frac {x^{4}}{2}}\right)\;\mathrm {d} x=\ldots ={\frac {13}{20}}.$
  
  Si l'on avait choisi de redéfinir D comme un domaine simple orthogonal à l'axe x = 0, on aurait obtenu les expressions :
  
      $D=\{(x,y)\mid y\leq 1,0\leq x\leq {\sqrt {y}}\},$
      $\iint _{D}(x+y)\;\mathrm {d} x\mathrm {d} y=\int _{0}^{1}\left(\int _{0}^{\sqrt {y}}(x+y)\mathrm {d} x\right)\,\mathrm {d} y,$
  
  conduisant au même résultat.
  
  **Exemples**
  
  **Réduction d'une intégrale triple**
  
  Considérons le domaine T de $ℝ^3$ orthogonal au plan xy (Fig. 5)
  
      $T=\{(x,y,z)\mid (x,y)\in D,{\rm {a}}(x,y)<z<{\rm {b}}(x,y)\}.$
  
  Sa projection sur xy est le domaine D de $ℝ^2$ orthogonal à x défini par
  
      $D=\{(x,y)\in \mathbb {R} ^{2}\mid a<x<b,f(x)<y<g(x)\}.$
      $\iiint _{T}F(x,y,z)\;\mathrm {d} x\mathrm {d} y\mathrm {d} z=\iint _{D}\left(\int _{{\rm {a}}(x,y)}^{{\rm {b}}(x,y)}F(x,y,z)\;\mathrm {d} z\right)\,\mathrm {d} x\mathrm {d} y=\int _{a}^{b}\left(\int _{f(x)}^{g(x)}\left(\int _{{\rm {a}}(x,y)}^{{\rm {b}}(x,y)}F(x,y,z)\;\mathrm {d} z\right)\,{\rm {d}}y\right)\,{\rm {d}}x.$
  
  ### Changements de variables usuels
  
  Lorsque les méthodes précédentes ne sont pas efficaces (domaine difficile à exprimer ou fonction difficile à intégrer), un changement de variables permet parfois de reformuler l'intégrale de façon plus commode.
  
  **Coordonnées polaires**
  
  Dans $ℝ^2,$ si le domaine présente une symétrie circulaire (c.-à-d. décrit un secteur annulaire) et que la fonction a des caractéristiques particulières, on peut appliquer une transformation des coordonnées cartésiennes en coordonnées polaires, pour changer la forme du domaine et faciliter l'intégration.
  
  D'après la définition de ce système de coordonnées, la transformation à effectuer est :
  
  $(x,y)=\Phi (r,\theta )=(r\cos \theta ,r\sin \theta )$
  
  (avec r positif, et θ généralement choisi entre 0 et 2π).
  
  Soit f(x, y) = $x^2$ + $y^2.$ Dans ce cas,
  
      $(f\circ \Phi )(r,\theta )=r^{2}(\cos \theta ^{2}+\sin \theta ^{2})=r^{2}$
  
  (grâce aux identités trigonométriques).
  
  **Coordonnées polaires**
  
  **Calcul du nouveau domaine**
  
  Le calcul du nouveau domaine U = $Φ^{−1}(V)$ s'effectue de façon automatique, en remplaçant x et y en fonction de r et θ dans les conditions sur (x, y) qui définissent l'ancien domaine V.
  
  Si D = {(x, y) | $a^2$ < $x^2$ + $y^2$ < $b^2,$ y > 0} avec b ≥ a ≥ 0 (Fig. 7), T := $Φ^{−1}(D)$ = {(r, θ) | $a^2$ < $r^2$ < $b^2,$ r sinθ > 0} = {(r, θ) | a < r < b, sinθ > 0} = ]a, b[×]0, π[.
  
  **Coordonnées polaires**
  
  **Jacobien, intégration**
  
  La matrice jacobienne de l'application Φ au point (r, θ) est :
  
      $J_{\Phi }(r,\theta )={\begin{pmatrix}\cos \theta &-r\sin \theta \\\sin \theta &r\cos \theta \\\end{pmatrix}}.$
  
  Son déterminant est donc égal à r, et l'on a, sous les hypothèses du théorème de changement de variables :
  
  $\iint _{V}f(x,y)\;\mathrm {d} x\mathrm {d} y=\iint _{U}f(r\cos \theta ,r\sin \theta )r\;\mathrm {d} r\mathrm {d} \theta .$
  
  Soit f(x, y) = y à intégrer sur l'ouvert D de l'exemple précédent, qui correspond via Φ à l'ouvert T (r ∈ ]a, b[ et θ ∈ ]0, π.
  
      $(f\circ \Phi )(r,\theta )=r\sin \theta$ donc
      $\iint _{D}y\;\mathrm {d} x\mathrm {d} y=\iint _{T}r\sin \theta \,r\;\mathrm {d} r\mathrm {d} \theta =\int _{0}^{\pi }\sin \theta \;\mathrm {d} \theta \int _{a}^{b}r^{2}\,\mathrm {d} r=[-\cos \theta ]_{0}^{\pi }\left[{\frac {r^{3}}{3}}\right]_{a}^{b}={\frac {2(b^{3}-a^{3})}{3}}.$
  
  **Coordonnées cylindriques**
  
  Dans $ℝ^3,$ l'intégration sur des domaines ayant pour base une portion de disque peut s'effectuer via un passage en coordonnées cylindriques :
  
  $(x,y,z)=\Phi (r,\theta ,z)=(r\cos \theta ,r\sin \theta ,z).$
  
  Puisque la composante z est inchangée, le calcul du nouveau domaine et du jacobien sont calqués sur ceux du passage en coordonnées polaires, et l'on a :
  
  $\iiint _{V}f(x,y,z)\;\mathrm {d} x\mathrm {d} y\mathrm {d} z=\iiint _{U}f(r\cos \theta ,r\sin \theta ,z)r\;\mathrm {d} r\mathrm {d} \theta \mathrm {d} z.$
  
  Il est conseillé d'utiliser cette méthode dans les cas de domaines cylindriques, coniques, ou tout du moins de régions pour lesquelles il est commode tant de délimiter l'intervalle des z que de transformer la base circulaire et la fonction.
  
  Soit f(x, y) = $x^2$ + $y^2$ + z à intégrer sur V = { (x,y,z) | $x^2$ + $y^2$ ≤ 9, 0 ≤ z ≤ 5 }. Alors,
  
      $(f\circ \Phi )(r,\theta ,z)=r^{2}+z$
  
  et
  
      $U:=\Phi ^{-1}(V)=\{(r,\theta ,z)\mid 0\leq r\leq 3,\ 0\leq \theta \leq 2\pi ,0\leq z\leq 5\}.$
  
  En toute rigueur, U et V ne sont pas ouverts, le jacobien r s'annule en certains points de U, et Φ n'est même pas injective de U vers V. Mais il suffirait, pour remédier à tous ces problèmes, de remplacer le parallélépipède U par son intérieur et le cylindre V par Φ(U), et cela ne modifierait pas les intégrales (puisque les parties fermées qu'on supprime sont négligeables). Ce genre de raccourci, fréquent dans l'application du théorème de changement de variables proprement dit, sera désormais tacite. On obtient donc :
  
      $\iiint _{V}(x^{2}+y^{2}+z)\ \mathrm {d} x\mathrm {d} y\mathrm {d} z=\iiint _{U}(r^{2}+z)r\ \mathrm {d} r\mathrm {d} \theta \mathrm {d} z=\int _{0}^{2\pi }\mathrm {d} \theta \int _{0}^{3}\left(\int _{0}^{5}(r^{3}+zr)\;\mathrm {d} z\right)\,\mathrm {d} r=2\pi \int _{0}^{3}(5r^{3}+{\tfrac {25r}{2}})\mathrm {d} r=2\pi {\frac {5.3^{4}+25.3^{2}}{4}}=315\pi .$
  
  **Coordonnées sphériques**
  
  Certains domaines de $ℝ^3$ présentent une symétrie sphérique, ainsi il est possible d'exprimer les coordonnées de leurs points à l'aide de deux angles θ (longitude) et φ (colatitude) et d'une distance à l'origine ρ. Le passage en coordonnées sphériques est la transformation :
  
  $(x,y,z)=\Phi (\rho ,\theta ,\varphi )=(\rho \cos \theta \sin \varphi ,\rho \sin \theta \sin \varphi ,\rho \cos \varphi )$
  
  (avec ρ positif, θ compris entre 0 et 2π, et φ entre 0 et π) et son jacobien est :
  
      $\det J_{\Phi }(\rho ,\theta ,\varphi )={\begin{vmatrix}\cos \theta \sin \varphi &-\rho \sin \theta \sin \varphi &\rho \cos \theta \cos \varphi \\\sin \theta \sin \varphi &\rho \cos \theta \sin \varphi &\rho \sin \theta \cos \varphi \\\cos \varphi &0&-\rho \sin \varphi \end{vmatrix}}=-\rho ^{2}\sin \varphi .$
  
  Dès lors, on a (en se rappelant qu'il faut prendre la valeur absolue du jacobien) :
  
  $\iiint _{V}f(x,y,z)\;\mathrm {d} x\mathrm {d} y\mathrm {d} z=\iiint _{U}f(\rho \cos \theta \sin \varphi ,\rho \sin \theta \sin \varphi ,\rho \cos \varphi )\,\rho ^{2}\sin \varphi \;\mathrm {d} \rho \mathrm {d} \theta \mathrm {d} \varphi .$
  
  Il est conseillé d'utiliser cette méthode dans le cas de domaines sphériques et de fonctions facilement simplifiables à l'aide des identités trigonométriques, comme dans l'exemple du calcul du volume d'une boule ou l'exemple ci-dessous.
  
  Intégrons f(x, y, z) = $x^2$ + $y^2$ sur V = {(x, y, z) | $x^2$ + $y^2$ + $z^2$ ≤ $R^2}$ (boule de rayon R centrée à l'origine).
  
      $U:=\Phi ^{-1}(V)=\{(\rho ,\theta ,\varphi )\mid 0\leq \rho \leq R,\;0\leq \theta \leq 2\pi ,\;0\leq \varphi \leq \pi \},$
      $(f\circ \Phi )(\rho ,\theta ,\varphi )=\rho ^{2}\cos ^{2}\theta \sin ^{2}\varphi +\rho ^{2}\sin ^{2}\theta \sin ^{2}\varphi =\rho ^{2}\sin ^{2}\varphi ,$
      $\iiint _{U}\rho ^{2}\sin ^{2}\varphi \,\rho ^{2}\sin \varphi \;\mathrm {d} \rho \mathrm {d} \theta \mathrm {d} \varphi =\left[{\frac {\rho ^{5}}{5}}\right]_{0}^{R}[\theta ]_{0}^{2\pi }\int _{0}^{\pi }\sin ^{3}\varphi \,\mathrm {d} \varphi .$
  
  On termine le calcul en utilisant par exemple le changement de variable t = –cos φ :
  
      $\int _{0}^{\pi }\sin ^{3}\varphi \,\mathrm {d} \varphi =\int _{-1}^{1}(1-t^{2})\,\mathrm {d} t=\left[t-{\frac {t^{3}}{3}}\right]_{-1}^{1}={\frac {4}{3}},$
  
  et l'on conclut :
  
      $\iiint _{V}f={\frac {R^{5}}{5}}\times 2\pi \times {\frac {4}{3}}={\frac {8\pi R^{5}}{15}}.$
  
  (Le passage en coordonnées cylindriques mènerait à des calculs un peu plus compliqués.)
  
  ## Exemples d'applications : calculs de volumes
  
  Grâce aux méthodes décrites précédemment, il est possible de démontrer les règles de calcul du volume de certains solides, en calculant l'intégrale triple sur ces solides la fonction constante 1.
  
  Le volume d'un cylindre droit P = D×[0, h] (de base D et de hauteur h) se ramène à l'intégrale double $∬_Dh$ dxdy sur le domaine D du plan xy. On retrouve ainsi, dans le cas particulier d'un cylindre droit, la formule classique : Volume d'un cylindre = aire base × hauteur.
  
  Le volume d'une boule de rayon R est :
  
      $V=\int _{0}^{2\pi }\mathrm {d} \theta \int _{0}^{\pi }\sin \varphi \mathrm {d} \varphi \int _{0}^{R}\rho ^{2}\mathrm {d} \rho =2\pi \left[-\cos \varphi \right]_{0}^{\pi }{\frac {R^{3}}{3}}={\frac {4\pi R^{3}}{3}}.$
  
  Grâce aux formules de réduction, on peut calculer le volume du 3-simplexe (le tétraèdre rectangle dont le sommet est à l'origine et les arêtes génératrices sont placées le long des trois axes et de longueur 1), en considérant, par exemple, l'orthogonalité selon le plan xy et l'axe x :
  
      $V=\int _{0}^{1}\left(\int _{0}^{1-x}\!\left(\int _{0}^{1-x-y}\!\mathrm {d} z\right)\,\mathrm {d} y\right)\,\mathrm {d} x=\int _{0}^{1}\left(\int _{0}^{1-x}\!(1-x-y)\,\mathrm {d} y\right)\,\mathrm {d} x=\int _{0}^{1}{\frac {(1-x)^{2}}{2}}\,\mathrm {d} x={\frac {1}{6}},$
  
  ce qui est conforme à la formule générale : Volume d'un cône = aire base × hauteur / 3, soit ici : (1/2) × 1 / 3.
  
  ## Intégrale multiple impropre
  
  Dans le cas de domaines non bornés ou d'intégrandes non bornés sur le bord du domaine, on parle d'intégrale multiple impropre.
  
  