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 ${\displaystyle n}$ (ou son hypervolume si ${\displaystyle n>3}$) comme l'intégrale d'une fonction de ${\displaystyle 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).

Soit ${\displaystyle f:D\rightarrow \mathbb {R} }$ une fonction bornée sur une partie bornée ${\displaystyle D}$ de ${\displaystyle \mathbb {R} ^{n}}$. Pour définir (au sens de Riemann) l'intégrabilité de ${\displaystyle f}$ et, le cas échéant, la valeur de son intégrale, notée

${\displaystyle \int \!\dots \!\int _{D}f(x_{1},\dots ,x_{n})\;\mathrm {d} x_{1}\dots \mathrm {d} x_{n}}$

ou plus simplement ${\displaystyle \int f}$, on se ramène d'abord au cas où ${\displaystyle D}$ est un pavé fermé ${\displaystyle P=I_{1}\times \ldots \times I_{n}}$, en choisissant un tel ${\displaystyle P}$ qui contient ${\displaystyle D}$ et en prolongeant ${\displaystyle f}$ par ${\displaystyle 0}$ (on démontre après coup que les définitions suivantes sont indépendantes du choix de ${\displaystyle P}$). Les définitions sont ensuite identiques, mutatis mutandis, à celles dans le cas ${\displaystyle 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^[1]. La première, par exemple, s'écrit :

${\displaystyle \int f=\int _{*}f=\int ^{*}f}$

(c'est-à-dire que ${\displaystyle f}$ est intégrable si et seulement si ses intégrales inférieure, ${\displaystyle \int _{*}f}$, et supérieure, ${\displaystyle \int ^{*}f}$, sont égales, et l'intégrale de ${\displaystyle f}$ est alors par définition cette valeur commune), et la troisième s'écrit :

${\displaystyle \int f=\lim _{\delta \to 0}\sum _{1\leq k\leq N}f(\xi _{k})\operatorname {vol} (P_{k}),}$

où chaque ${\displaystyle (P_{k},\xi _{k})_{1\leq k\leq 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 ${\displaystyle D}$ de dimension ${\displaystyle n=1}$, ${\displaystyle 2}$ ou ${\displaystyle 3}$, de la fonction constante ${\displaystyle 1}$, est respectivement la longueur, l'aire ou le volume de ${\displaystyle D}$.

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 ${\displaystyle \mathbb {R} ^{n}}$ est Riemann-intégrable si et seulement si l'ensemble des points de discontinuité de son prolongement par ${\displaystyle 0}$ à ${\displaystyle \mathbb {R} ^{n}}$ est négligeable^[2]. En particulier :

En effet, une partie de ${\displaystyle \mathbb {R} ^{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.

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 ${\displaystyle D}$ est un compact connexe de ${\displaystyle \mathbb {R} ^{n}}$, cubable mais non négligeable, et ${\displaystyle f:D\rightarrow \mathbb {R} }$ une fonction continue, alors il existe un point ${\displaystyle m}$ de ${\displaystyle D}$ en lequel la « valeur moyenne » ou « moyenne intégrale » de ${\displaystyle f}$ sur ${\displaystyle D}$ est atteinte, c'est-à-dire tel que

${\displaystyle 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 ${\displaystyle 1}$, à cela près qu'en dimension supérieure, on n'est plus assuré de trouver un ${\displaystyle m}$ intérieur à ${\displaystyle D}$.

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 ${\displaystyle \mathbb {R} \times \mathbb {R} }$ par ${\displaystyle \mathbb {R} ^{p}\times \mathbb {R} ^{q}}$. De plus, on peut évidemment intervertir les rôles des deux variables (mono- ou multidimensionnelles) ${\displaystyle x}$ et ${\displaystyle y}$.

Théorème^[3] — Pour toute application ${\displaystyle f:[a,b]\times [c,d]\rightarrow \mathbb {R} }$ :

${\displaystyle \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^[4] — Lorsque ${\displaystyle f}$ est Riemann-intégrable, les deux fonctions (intégrales inférieure et supérieure)

${\displaystyle [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 ${\displaystyle f}$.

L'énoncé de la conclusion se simplifie^[5] lorsque de plus, pour tout ${\displaystyle x\in [a,b]}$, la fonction ${\displaystyle y\mapsto f(x,y)}$ est intégrable (sur ${\displaystyle [c,d]}$) : alors, son intégrale (fonction de ${\displaystyle x\in [a,b]}$) est intégrable et

${\displaystyle \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.}$

Soient ${\displaystyle U}$ et ${\displaystyle V}$ deux ouverts de ${\displaystyle \mathbb {R} ^{n}}$, ${\displaystyle \Phi :U\rightarrow V}$ un C¹-difféomorphisme et ${\displaystyle f}$ une application de ${\displaystyle V}$ dans ${\displaystyle \mathbb {R} }$.

Théorème^[6] — Si ${\displaystyle U}$ et ${\displaystyle V}$ sont bornés et si ${\displaystyle f}$ et (f∘Φ) × |det J_Φ| sont Riemann-intégrables, alors

${\displaystyle \int _{V}f=\int _{U}(f\circ \Phi )\times {|}{\det }J_{\Phi }{|}.}$

Ces hypothèses sont vérifiées dès que ${\displaystyle U}$ et ${\displaystyle V}$ sont cubables, ${\displaystyle f}$ est continue bornée, et ${\displaystyle \det J_{\Phi }}$ est borné^[7].

Par linéarité, pour toute constante ${\displaystyle c}$, ${\displaystyle \int (cf)=c\int f.}$ En particulier, l'intégrale d'une fonction constante égale à ${\displaystyle c}$, sur un domaine cubable ${\displaystyle D}$ de ${\displaystyle \mathbb {R} ^{n}}$, est le produit de ${\displaystyle c}$ par le volume de ${\displaystyle 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 : ${\displaystyle \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 ${\displaystyle f}$ comme la somme des produits de ${\displaystyle 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 ${\displaystyle 1}$ : ${\displaystyle \int _{D_{1}\sqcup D_{2}}f=\int _{D_{1}}f+\int _{D_{2}}f.}$

Si la fonction est impaire par rapport à une variable ${\displaystyle x_{i}}$ et si le domaine est invariant par la symétrie par rapport à l'hyperplan ${\displaystyle x_{i}=0}$, l'intégrale s'annule.

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.

Si ${\displaystyle A}$ est une partie cubable de ${\displaystyle \mathbb {R} ^{n-1}}$ et si deux fonctions ${\displaystyle g}$ et ${\displaystyle h}$, définies sur ${\displaystyle A}$, sont Riemann-intégrables — par exemple si elles sont bornées et continues (éventuellement par morceaux) — alors la partie suivante de ${\displaystyle \mathbb {R} ^{n}}$, appelée un « domaine simple orhogonal à l'hyperplan ${\displaystyle t=0}$ », est cubable (Fig. 3) :

${\displaystyle 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 ${\displaystyle \mathbb {R} ^{n}}$, pour choisir laquelle joue le rôle de ${\displaystyle t}$.

On peut également, dans la définition de ${\displaystyle D}$, remplacer un ou deux ${\displaystyle <}$ par des ${\displaystyle \leq }$, car le graphe d'une fonction intégrable sur ${\displaystyle A}$ est Jordan-négligeable.

Si ${\displaystyle f}$ est une fonction continue bornée sur ce domaine simple ${\displaystyle D}$, les hypothèses du critère de Lebesgue sont vérifiées, et le théorème de Fubini devient :

${\displaystyle \int _{D}f=\int _{A}\left(\int _{g(x)}^{h(x)}f(x,t)\,{\rm {d}}t\right)\,{\rm {d}}x.}$

La région ${\displaystyle 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 ${\displaystyle y=0}$), associé aux deux fonctions ${\displaystyle g(x)=x^{2}}$ et ${\displaystyle h(x)=1}$, sur l'intervalle ${\displaystyle A=[0,1]}$ (on trouve l'extrémité ${\displaystyle 1}$, abscisse du point d'intersection des courbes de ${\displaystyle g}$ et ${\displaystyle h}$, en calculant la solution positive de l'équation ${\displaystyle g(x)=h(x)}$).

On se propose d'intégrer sur ${\displaystyle D}$ la fonction ${\displaystyle f(x,y)=x+y}$. Appliquons le théorème :

${\displaystyle \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 ${\displaystyle x}$ fixé :

${\displaystyle \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 :

${\displaystyle \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 ${\displaystyle D}$ comme un domaine simple orthogonal à l'axe ${\displaystyle x=0}$, on aurait obtenu les expressions :

${\displaystyle D=\{(x,y)\mid y\leq 1,0\leq x\leq {\sqrt {y}}\},}$

${\displaystyle \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.

Considérons le domaine ${\displaystyle T}$ de ${\displaystyle \mathbb {R} ^{3}}$ orthogonal au plan ${\displaystyle xy}$ (Fig. 5)

${\displaystyle T=\{(x,y,z)\mid (x,y)\in D,{\rm {a}}(x,y)<z<{\rm {b}}(x,y)\}.}$

Sa projection sur ${\displaystyle xy}$ est le domaine ${\displaystyle D}$ de ${\displaystyle \mathbb {R} ^{2}}$ orthogonal à ${\displaystyle x}$ défini par^[8]

${\displaystyle D=\{(x,y)\in \mathbb {R} ^{2}\mid a<x<b,f(x)<y<g(x)\}.}$

${\displaystyle \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.}$

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.

Dans ${\displaystyle \mathbb {R} ^{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 :

${\displaystyle (x,y)=\Phi (r,\theta )=(r\cos \theta ,r\sin \theta )}$

(avec ${\displaystyle r}$ positif, et ${\displaystyle \theta }$ généralement choisi entre ${\displaystyle 0}$ et ${\displaystyle 2\pi }$).

Le calcul du nouveau domaine U = Φ⁻¹(V) s'effectue de façon automatique, en remplaçant ${\displaystyle x}$ et ${\displaystyle y}$ en fonction de ${\displaystyle r}$ et ${\displaystyle \theta }$ dans les conditions sur ${\displaystyle (x,y)}$ qui définissent l'ancien domaine ${\displaystyle V}$.

La matrice jacobienne de l'application ${\displaystyle \Phi }$ au point ${\displaystyle (r,\theta )}$ est :

${\displaystyle J_{\Phi }(r,\theta )={\begin{pmatrix}\cos \theta &-r\sin \theta \\\sin \theta &r\cos \theta \\\end{pmatrix}}.}$

Son déterminant est donc égal à ${\displaystyle r}$, et l'on a, sous les hypothèses du théorème de changement de variables :

${\displaystyle \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 .}$

Dans ${\displaystyle \mathbb {R} ^{3}}$, l'intégration sur des domaines ayant pour base une portion de disque peut s'effectuer via un passage en coordonnées cylindriques :

${\displaystyle (x,y,z)=\Phi (r,\theta ,z)=(r\cos \theta ,r\sin \theta ,z).}$

Puisque la composante ${\displaystyle 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 :

${\displaystyle \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 ${\displaystyle z}$ que de transformer la base circulaire et la fonction.

Certains domaines de ${\displaystyle \mathbb {R} ^{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 ${\displaystyle \theta }$ (longitude) et ${\displaystyle \varphi }$ (colatitude) et d'une distance à l'origine ${\displaystyle \rho }$. Le passage en coordonnées sphériques est la transformation :

${\displaystyle (x,y,z)=\Phi (\rho ,\theta ,\varphi )=(\rho \cos \theta \sin \varphi ,\rho \sin \theta \sin \varphi ,\rho \cos \varphi )}$

(avec ${\displaystyle \rho }$ positif, ${\displaystyle \theta }$ compris entre ${\displaystyle 0}$ et ${\displaystyle 2\pi }$, et ${\displaystyle \varphi }$ entre ${\displaystyle 0}$ et ${\displaystyle \pi }$) et son jacobien est :

${\displaystyle \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) :

${\displaystyle \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.

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 ${\displaystyle 1}$.

Cette section est vide, insuffisamment détaillée ou incomplète. Votre aide est la bienvenue ! Comment faire ?

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.

Jacques Douchet et Bruno Zwahlen, Calcul différentiel et intégral : Fonctions réelles de plusieurs variables réelles, Lausanne, PPUR, 2004, 2^e éd. (1^re éd. 1998), 172 p. (ISBN 2-88074-257-9)

v · m Analyse vectorielle
Objets d'étude	Champ de vecteurs Champ scalaire Équation aux dérivées partielles Équation de Laplace Équation de Poisson
Opérateurs	Nabla Gradient Rotationnel Rotationnel du rotationnel Divergence Laplacien scalaire Bilaplacien Laplacien vectoriel D'alembertien Advection
Théorèmes	Théorème de Green Théorème de Stokes Théorème de Helmholtz-Hodge Théorème de la divergence Théorème du gradient Théorème du rotationnel

• Portail de l'analyse