Notes LM216

Fonctions de plusieurs variables
Notes du Cours LM 216

Albert Cohen
Dans ce cours, on s’intéresse aux fonctions de plusieurs variables
(x1 , · · · , xn ) 7→ f (x1 , · · · , xn ).
Celles-ci interviennent naturellement pour décrire la dépendance de grandeurs

en fonctions de plusieurs paramètres : par exemple la distribution de la
température ou de la pression atmosphérique en fonction de la position dans
l’espace et du temps. Elles peuvent aussi être définies par des expressions
similaires à celles que l’on rencontre pour
√ les fonctions d’une variable : par
2
exemple f (x, y, z) = sin(xy + 1) + z 4 + 3x2 est une fonction des trois
variables x, y et z.
L’objectif principal de ce cours est de généraliser aux fonctions de plusieurs
variables les notions de dérivation et d’intégration qui ont été abordées dans
le cadre des fonctions d’une variable.
Ces notes contiennent la totalité des résultats du cours sous une forme
très condensée. En particulier, la plupart des démonstrations détaillées en
amphi sont omises, et c’est un excellent exercice pour l’étudiant que d’essayer
de les refaire uniquement à partir des indications données dans les notes. La
présentation est orientée sur les fonctions de n variables, où n est un nombre
arbitraire, mais l’étudiant est encouragé à traduire les résultats dans le cas
de fonctions de 2 ou 3 variables afin de s’en forger une intuition visuelle.
Avertissement : ces notes sont régulièrement mises à jour et corrigées.

Toutes les remarques permettant d’améliorer la rédaction peuvent être en-
voyées à l’adresse [email protected].
1 Notions de topologie dans IRn

1.1 Vecteurs et normes dans IRn
Rappelons qu’en dimension 2, on identifie un vecteur x de coordonnées
(x1 , x2 ) avec un point du plan de coordonnées (x1 , x2 ) une fois fixée une
origine. On généralisera ici cette identification en désignant le point ou le
vecteur de coordonnées (x1 , · · · , xn ) par x = (x1 , · · · , xn ) ∈ IRn .
1
La norme euclidienne d’un vecteur x = (x1 , · · · , xn ) ∈ IRn est définie par
q
kxk = x21 + · · · + x2n
Le produit scalaire de x avec y = (y1 , · · · , yn ) est défini par
x · y = x1 y1 + · · · + xn yn .
on a en particulier kxk2 = x · x.
Inégalité de Cauchy-Schwartz : pour tout x et y dans IRn
|x · y| ≤ kxk kyk,
avec égalité si et seulement si x et y sont proportionels.
Inégalité triangulaire : pour tout x et y dans IRn
kx + yk ≤ kxk + kyk.
On rappelle aussi que si θ désigne l’angle entre deux vecteurs x et y non-nuls,

on a
x·y
cos(θ) = .
kxk kyk
D’un point de vue géométrique, la norme euclidienne kxk correspond à
la longueur du vecteur x (ou encore à la distance du point x à l’origine). On
définit ainsi la distance euclidienne entre deux points x et y par
d(x, y) = kx − yk,
c’est à dire la longeur du vecteur reliant le point x au point y.

La norme euclidienne n’est pas l’unique façon de mesurer la taille d’un
vecteur x.
Définition 1.1 On appelle norme sur IRn une application N de IRn dans
IR+ vérifiant les trois propriétés suivantes
1. N (x) = 0 si et seulement si x = 0.
2. N (λx) = |λ|N (x) pour tout x ∈ IRn et λ ∈ IR.
2
3. N (x + y) ≤ N (x) + N (y) pour tout x, y ∈ IRn .
On vérifie facilement que la norme euclidienne vérifie ces trois propriétés.

Voici d’autres exemples de normes : la norme “sup” ou `∞ définie par
kxk∞ := sup |xi |,

i=1,···,n
la norme `p définie pour 1 ≤ p < +∞ par
kxkp := (|x1 |p + · · · + |xn |p )1/p .
On voit que le cas p = 2 correspond à la norme euclidienne.

Etant donné une norme N , un point x et un nombre r ≥ 0 on appelle
boule de centre x et de rayon r pour la norme N l’ensemble des points à
distance au plus r de x, c’est à dire
B(x, r) = {y ; N (y − x) ≤ r}.
On appelle boule unité la boule B(0, 1) de centre 0 et de rayon 1. En

l’absence de précision sur N , il s’agira systématiquement de la norme eucli-
dienne.
Définition 1.2 Deux normes N1 et N2 sont dites “équivalentes” si et seule-

ment si il existe deux constantes c, C > 0 telles que pour tout x ∈ IRn
cN1 (x) ≤ N2 (x) ≤ CN1 (x).
Cette propriété signifie intuitivement qu’un vecteur très petit dans une
norme le sera aussi dans l’autre. On vérifie par exemple aisément (en util-
isant l’inégalité de Cauchy-Schwarz) que
kxk2 ≤ kxk1 ≤ n1/2 kxk2 .
On a aussi
n−1/2 kxk2 ≤ kxk∞ ≤ kxk2 .
On a en fait la propriété fondamentale suivante.
Théorème 1.1 Toutes les normes sur IRn sont équivalentes entre elles.
3
1.2 Quelques notions de topologie
Soit (pk )k>0 une suite de points (ou de vecteurs) de IRn . On dit que cette
suite converge vers une limite p ∈ IRn si et seulement si pour tout ε > 0, il
existe k0 tel que
k ≥ k0 ⇒ kpk − pk ≤ ε.
On note alors limk→+∞ pk = p ou encore parfois pk → p. Dans le cas
contraire on dit que la suite diverge. Dans le cas n = 1, on a kpk − pk =
|pk −p| et on retrouve la définition usuelle de la convergence des suites réelles.
On vérifie aisément que la limite d’une suite si elle existe est unique.
Grâce à l’équivalence des normes, il est possible de remplacer la norme
euclidienne par n’importe quelle autre norme N sans alterer la définition de
la convergence. Il en est de même pour toutes les notions que nous allons
introduire dans cette section.
En particulier en utilisant la norme sup, on voit que si pk = (xk,1 , · · · , xk,n )
et p = (x1 , · · · , xn ), la convergence de la suite (pk )k≥0 vers p est équivalente
à la convergence pour tout i = 1, · · · , n de la suite des i-ème coordonnées
(xk,i )k≥0 vers la coordonnée xi .
Une suite (pk )k>0 est dite “de Cauchy” si et seulement si pour tout ε > 0,
il existe k0 tel que
k, l ≥ k0 ⇒ kpk − pl k ≤ ε.
On vérifie aisément que toute suite convergente est de Cauchy. La réciproque
est aussi vraie : l’espace IRn est complet ce qui signifie que toute suite de
Cauchy est convergente. Cette propriété se déduit aisément de la propriété
analogue connue pour IR en utilisant la norme sup. L’intérêt pratique de
cette propriété est qu’il n’est ainsi pas nécessaire de connaitre la limite d’une
suite pour pouvoir prouver sa convergence.
Définition 1.3 Un ensemble U ⊂ IRn est dit “ouvert” si et seulement si

pour tout x ∈ U il existe r > 0 tel que B(x, r) ⊂ U .
Par convention l’ensemble vide ∅ est ouvert.
Pour un ensemble E quelconque on dit que x est un point intérieur à E

si il existe r > 0 tel que B(x, r) ⊂ E. L’ensemble des points intérieur à E
◦ ◦
est appelé l’intérieur de E et est noté E . On vérifie que E est le plus grand
ouvert contenu dans E et que E est ouvert si et seulement si il est égal à
son intérieur.
4
Définition 1.4 Un ensemble F ⊂ IRn est dit “fermé” si et seulement si son
complémentaire F c = IRn \ F est ouvert.
Théorème 1.2 F est fermé si et seulement toute suite de point de F qui

converge a sa limite contenue dans F .
Pour un ensemble E quelconque on dit que x est une valeur d’adhérence

de E si il existe une suite (xk )k≥0 de points de E qui converge vers x.
L’ensemble des valeurs d’adhérence de E est noté E et est appelé l’adhérence
(ou la fermeture) de E. On vérifie que E est le plus petit fermé contenant
E et que E est fermé si et seulement si E = E.
Quelques propriétés élémentaires:
1. Toute union finie ou infinie d’ouverts est un ouvert.
2. Toute union finie de fermés est un fermé
3. Toute intersection finie ou infinie de fermés est un fermé.
4. Toute intersection finie de d’ouverts est un ouvert
5. Les seuls ensembles à la fois ouverts et fermés sont ∅ et IRN .
6. Un ensemble fini de points de IRn est fermé.

◦
7. On a E ⊂ E ⊂ E.
Définition 1.5 On appelle frontiere d’un ensemble E l’ensemble des points

de son adhérence qui ne sont pas dans son intérieur, c’est à dire
◦ ◦c
∂E := E\ E = E∩ E .
Définition 1.6 Un ensemble E ⊂ IRn est dit borné si et seulement si il

existe R > 0 tel que E ⊂ B(0, R). Un ensemble fermé et borné de IRn est
dit “compact”.
Théorème 1.3 Pour E ⊂ IRn , les trois propriétés suivantes sont équivalentes
1. E est compact
2. Propriété de Bolzano-Weierstrass : toute suite (xk )k≥0 de points de E

admet une sous-suite qui converge vers une limite x ∈ E.
5
3. Propriété de Borel-Lebesgue : de toute famille d’ouverts (Ui )i∈I telle
que E ⊂ ∪i∈I Ui , on peut extraire une sous famille finie (U1 , · · · , Um )
telle que E ⊂ U1 ∪ · · · ∪ Um .
Définition 1.7 Un ensemble E est dit “connexe” si et seulement si il n’existe

aucune paire d’ouverts (U1 , U2 ) tels que E ⊂ U1 ∪ U2 et E ∩ U1 et E ∩ U2
sont disjoints.
Autrement dit, E est connexe si on ne peux pas le séparer en deux partie

disjointes en l’intersectant avec deux ouverts. Dans le cas d’un ensemble
ouvert, cela se traduite par la propriété intuitive suivante.
Théorème 1.4 Un ensemble U est ouvert et connexe si et seulement si

pour tout x et y dans U , il existe une ligne brisée contenue dans U qui les
relie, c’est à dire un ensemble fini de points {p1 , · · · , pm } tel que p1 = x et
pm = y, et tel que les segments d’extremités pi et pi+1 sont tous contenus
dans U .
Définition 1.8 Un ensemble E est dit “convexe” si et seulement si pour

tout x, y ∈ E et t ∈ [0, 1] on a tx + (1 − t)y ∈ E, i.e. le segment d’extremité
x et y est contenu dans E.
On peut vérifier qu’en dimension n = 1 les connexes et les convexes de

IR sont exactement les intervalles (de taille finie ou infinie). En revanche, en
dimension n > 1 tout convexe est connexe mais la réciproque est fausse.
On vérifie aussi qu’une intersection finie ou infinie de compacts est un
compact et de même pour les convexes.
2 Fonctions continues
2.1 Définition et propriétés fondamentales
On va s’intéresser à des fonctions de plusieurs variables et à valeurs réelles
f : IRn → IR : pour x = (x1 , · · · , xn ) ∈ IRn , f (x) = f (x1 , · · · , xn ) ∈ IR.
On étend à ce cadre la notion de domaine de définition connue pour les
fonctions d’une pvariable. Par exemple le domaine de définition de la fonction
f (x, y) = log( 3 − x2 − y 2 ) est l’ensemble D = {(x, √ y) ; x2 + y 2 < 3} c’est
à dire l’intérieur de la boule de centre 0 et de rayon 3.
On s’intéressera aussi parfois à la notion plus générale de fonctions de n
variables et à valeurs vectorielles f : IRn → IRm : pour x = (x1 , · · · , xn ) ∈ IRn
6
f (x) = (f1 (x), · · · , fm (x)) ∈ IRm . Notons que chacune des fonctions fi est
une fonction de n variables et à valeur réelle. On dit parfois que f est un
champ de vecteurs à m composantes définis sur IRn .
Si E ⊂ IRn on appelle image de E par f l’ensemble f (E) des images par
f des points de E. On a donc
f (E) = {f (x) ; x ∈ E}.
On dit que la fonction f est bornée sur E si f (E) est un ensemble borné,
autrement dit, il existe R > 0 tel que kf (x)k ≤ R pour tout x ∈ E.
Si F ⊂ IRm , on appelle image réciproque de F par f l’ensemble f −1 (F )
des antécédents par f des points de F . On a donc
f −1 (F ) = {x ∈ IRn ; f(x) ∈ F}.
Une fonction d’une variable et à valeur réelle est décrite par son graphe,
qui est le sous-ensemble de IR2 défini par
Gf = {(x, f (x)) ; x ∈ IR},
et que l’on représente par la courbe du plan d’équation y = f (x). Dans le

cas d’une fonction de deux variables à valeur réelle on définit de même le
graphe
Gf = {(x, y, f (x, y)) ; (x, y) ∈ IR2 },
que l’on peut représenter comme une surface d’équation z = f (x, y) dans
l’espace à 3 dimensions. Pour λ ∈ IR, on appelle courbe de niveau λ de la
fonction f l’ensemble Cλ des points du plan dont l’image par f vaut λ c’est
à dire
Cλ := {(x, y) ; f (x, y) = λ} = f −1 ({λ}).
Les notions de graphe et de courbes de niveau s’étendent de manière évidentes
au cas des fonctions de n variables.
Définition 2.1 Une fonction f : IRn → IR est continue au point x =

(x1 , · · · , xn ) ∈ IRn si et seulement si pour tout ε > 0, il existe α > 0 tel
que pour tout y = (y1 , · · · , yn ) ∈ IRn ,
ky − xk ≤ α ⇒ |f (y) − f (x)| ≤ ε.
Une fonction f est continue sur un ensemble E si et seulement si elle est

continue en tout point de cet ensemble. L’ensemble des fonctions continues
sur E est noté C(E) ou C 0 (E) .
7
Dans le cas de fonctions f : IRn → IRm à valeur vectorielle, il faut
remplacer |f (y) − f (x)| par kf (y) − f (x)k dans la définition ci-dessus. Grâce
à l’équivalence des normes, il est possible de remplacer la norme euclidienne
par n’importe quelle autre norme N sans alterer la définition de la continuité.
On remarque que si f est continue au point x alors elle est toujours bornée
sur une boule B(x, r) pour r suffisament petit.
Une manière équivalente d’exprimer la continuité au point x est de dire
que pour tout suite (xk ) telle que xk → x, on a f (xk ) → f (x). On peut aussi
démontrer que f est continue sur IRn si et seulement si l’image réciproque
de tout ensemble ouvert est un ensemble ouvert.
On démontre enfin que l’image d’un connexe par une fonction continue
est aussi connexe, et que l’image d’un compact par une fonction continue est
aussi compact. Cette dernière propriété entraine qu’une fonction continue
est toujours bornée sur un compact.
Définition 2.2 Une fonction f : IRn → IRm est lipschtizienne si et seule-

ment si pour tout x, y ∈ IRn on a
kf (x) − f (y)k ≤ Ckx − yk,
où C est une constante indépendante de x et y.
On dit parfois que f est C-lipschitzienne. On voit immédiatement que

toute fonction lipschitzienne est continue. En revanche la réciproque est
fausse. Un exemple de fonction 1-lipschitzienne est l’application x 7→ kxk
(cela se déduit aisément de l’inégalité triangulaire).
Si f : IRn → IRm est une application linéaire, alors en introduisant la
base cannonique (e1 , · · · , en ) de IRn et en notant M = max kf (ei )k, on voit
que pour tout x = (x1 , · · · , xn ) ∈ IRn on a
n
X n
X
kf (x)k = k xi f (ei )k ≤ M |xi | ≤ Bkxk,
i=1 i=1
√
avec B = M n. En remplacant x par x − y on voit ainsi que toute ap-
plication linéaire est lipschitzienne et donc continue. On appelle norme de
l’application linéaire f la quantité
kf (x)k
|||f ||| := sup .
x6=0 kxk
On vérifie que |||f ||| est la plus petite constante C telle que kf (x)k ≤ Ckxk
pour tout x.
8
2.2 Opérations sur les fonctions continues
Voici trois propriétés dont les démonstrations sont élémentaires :
1. Une combinaison linéaire ou une multiplication de deux (ou d’un nom-
bre fini de) fonctions continues en un point x est continue en x.
2. Si f : IRn → IRm est continue au point x et g : IRm → IRp est continue

au point f (x), alors g ◦ f : IRn → IRp est continue au point x.
3. La fonction f : IRn → IR qui à x = (x1 , · · · , xn ) associe la i-ème

coordonnée f (x) = xi est continue. On l’appelle la i-ème application
coordonnée
Ces propriétés permettent de construire de nombreuses fonctions contin-
ues en utilisant les résultats connus sur la continuité
p des fonctions d’une seule
variable. Par exemple, la fonction f (x, y) = log( 3 − x2 − y 2 ) est continue
sur son domaine de définition, par composition des applications coordonnées,
carré, racine, logarithme. De même, on voit que si f : IRn → IR est continue
au point x et f (x) 6= 0, alors 1/f est aussi continue au point x.
Une classe importante de fonction continue est celle des polynômes. On
rappelle qu’un polynôme d’une variable est une fonction de la forme
p(x) = a0 + a1 x + · · · + am xm ,
c’est à dire une combinaison linéaire d’un nombre fini de fonctions puissance
xk pour k entier positif. L’entier m est appelé “degré” du polynôme p. Le
terme ak xk est appelé monôme de degré k. Dans le cas de fonctions à
n variable, un polynôme est une combinaison linéaire d’un nombre fini de
fonctions du type xk11 xk22 · · · xknn où k1 , · · · , kn sont des entiers positifs. Le
degré de xk11 · · · xknn est par définition k1 + · · · + kn . Le degré de p est celui du
monôme de plus haut degré présent dans p. En notant k = (k1 , · · · , kn ) ∈
INn , on voit ainsi qu’un polynôme de degré m est de la forme
ak xk11 · · · xknn .
X
p(x) = p(x1 , · · · , xn ) =
k1 +···+kn ≤m
Une manière naturelle d’étudier les variations d’une fonction f de n

variables consiste à fixer n − 1 d’entre elles et à étudier alors la fonction par
rapport à l’unique variable restante. Ainsi si x = (x1 , · · · , xn ) est un point
fixé, on définit la i-ème application partielle au point x par
fi,x (t) = f (x1 , · · · , xi−1 , t, xi+1 , · · · , xn ).
9
Il est facile de vérifier que la continuité de f en x entraine celle de l’application
partielle fi,x au point t = xi . En revanche, la réciproque est fausse : par ex-
emple la fonction définie par f (x, y) = x4xy +y 4
si (x, y) 6= (0, 0) et f (0, 0) = 0
admet des applications partielles continues en tout point, mais cependant
cette fonction n’est pas continue au point (0, 0).
2.3 Fonctions continues et ensembles compacts

Une fonction est dite uniformément continue sur un ensemble E si et seule-
ment si pour tout ε > 0, il existe α > 0 tel que pour tout x, y ∈ E
kx − yk ≤ α ⇒ |f (x) − f (y)| ≤ ε,
(kf (x) − f (y)k dans le cas des fonctions à valeurs vectorielles). Le point
important est que α dépend seulement de ε, pas de x et y. Ainsi une
fonction uniformément continue est continue mais la réciproque est fausse
en générale. On a cependant le résultat suivant que l’on peut démontrer à
l’aide de la propriété de Borel-Lebesgue.
Théorème 2.1 Si E est un ensemble compact, toute fonction continue sur
E est aussi uniformément continue.
Une fonction f à valeur réele est dite minorée si et seulement si il existe
M tel que f (x) ≥ M pour tout x ∈ E. Dans ce cas, l’ensemble f (E) admet
un infimum qui est aussi le plus grand M ayant la propriété précèdente.
Cette valeur est notée inf x∈E f (x) et est caractérisée par les propriétés
1. Pour tout x ∈ E on a f (x) ≥ inf x∈E f (x)
2. Il existe une suite (xk ) de points de E telle que f (xk ) → inf x∈E f (x).
Dans le cas où f n’est pas minorée, on pose inf x∈E f (x) = −∞. Dans le cas
où l’infimum est fini, il peut-être atteint, c’est à dire qu’il existe x0 ∈ E tel
que f (x0 ) = inf x∈E f (x), ou ne pas être atteint. Si il est atteint on parle de
minimum de f sur E et on note alors f (x0 ) = minx∈E f (x). On définit de
manière similaire le suprémum supx∈E f (x) sur E (fini ou égal à +∞) que
l’on appelle maximum maxx∈E f (x) dans le cas ou il est atteint. Le théorème
suivant peut se démontrer à l’aide de la propriété de Bolzano-Weierstrass.
Théorème 2.2 Si E est un ensemble compact, toute fonction continue sur
E admet un minimum et un maximum.
On peut utiliser ce résultat afin de prouver l’équivalence de toutes les
normes définies sur IRn , en montrant leur équivalence avec la norme sup.
10
3 Différentiabilité et dérivées partielles
3.1 Définitions et propriétés élémentaires
Définition 3.1 Soit f une fonction à valeur réelle définie sur un ouvert
U ⊂ IRn . On dit que f admet une dérivée partielle au point x = (x1 , · · · , xn ) ∈
U par rapport à la i-ème coordonnée si la i-ème application partielle fi,x (t)
associée au point x est dérivable au point t = xi . On note
∂f 0 f (x1 , · · · , xi−1 , xi + h, xi+1 , · · · , xn ) − f (x)
(x) = fi,x (xi ) = lim .
∂xi h→0 h
∂f
Le calcul de ∂xi consiste donc à ne dériver l’expression de f que par
∂f
rapport à la variable xi . Notons que les dérivées partielles ∂x i
sont aussi
des fonctions de n variables à valeurs dans IR. Par exemple si f (x, y, z) =
−2x cos(y), on a ∂f ∂f ∂f
∂x (x, y, z) = −2 cos(y), ∂y (x, y, z) = 2x sin(y) et ∂z (x, y, z) =
0.
Définition 3.2 Soit f une fonction à valeur réelle définie sur un ouvert
U ⊂ IRn et soit v ∈ IRn un vecteur fixé. La dérivée de f en x ∈ U suivant
le vecteur v est la valeur de la dérivée de la fonction d’une variable fv (t) =
f (x + tv) en t = 0 si elle existe, c’est à dire limt→0 f (x+tv)−f
t
(x)
.
Intuitivement cette derniere quantité mesure les variations de f lorsqu’on

se déplace autour de x dans la direction du vecteur u à la “vitesse” kuk. On
∂f
voit en particulier que ∂xi
(x) correspond à la dérivée en x suivant le vecteur
de base ei (avec e1 = (1, 0, · · · , 0), e2 = (0, 1, 0, · · · , 0),...,en = (0, · · · , 0, 1)).
Définition 3.3 Soit f une fonction définie sur un ouvert U ⊂ IRn et à

valeur dans IRm . On dit que f est différentiable en un point x ∈ U si et
seulement si il existe une application linéaire de IRn dans IRm notée dfx telle
que pour tout h ∈ IRn suffisement petit
f (x + h) = f (x) + dfx (h) + khkε(h),
avec limh→0 ε(h) = 0. L’application dfx est appelée différentielle de f au

point x ou encore application linéaire tangente de f au point x.
La formule ci-dessus n’est rien d’autre qu’un développement limité de f

à l’ordre 1 au voisinage de x. Dans le cas d’une fonction d’une variable à
valeur réelle, cette définition redonne la notion classique de dérivabilité, et
11
on a alors dfx (h) = f 0 (x)h. Clairement, si f est différentiable au point x,
elle est continue en ce point.
Dans le cas d’une fonction de n variables et à valeur dans IR, dfx est
une forme linéaire de IRn dans IR. Pour h = (h1 , · · · , hn ) ∈ IRn on a donc
∂f
dfx (h) = a1 h1 + · · · + an hn . On vérifie aisément que ai = ∂xi
(x). Le vecteur
∂f ∂f
∇f (x) := ( (x), · · · , (x)),
∂x1 ∂xn
est appelé le gradient de f au point x. On a donc
n
X ∂f
f (x + h) = f (x) + ∇f (x) · h + khkε(h) = f (x) + (x)hi + khkε(h).
i=1
∂xi
On déduit en particulier que la dérivée de f en x suivant le vecteur v est

donnée par le produit scalaire ∇f (x) · v.
Ainsi une fonction différentiable en un point x admet des dérivées par-
tielles et plus généralement des dérivées suivant tout vecteur en ce point.
La réciproque est fausse: il existe des fonctions qui admettent des dérivées
suivant tout vecteur en un point mais ne sont pas même continues en ce
point. Un exemple est le comportement au point (0, 0) de la fonction définie
par f (x, y) = y 2 /x si x 6= 0 et f (0, y) = 0.
Dans le cas d’une fonction f de n variables à valeur dans IRm , on vérifie
que dfx est la matrice m × n dont l’élément de ligne i et de colonne j est
∂fi
donné par ∂x j
(x). Cette matrice est appelée matrice jacobienne de f au
point x. On peut écrire pour chaque composante de f
n
X ∂fi
fi (x + h) = fi (x) + (x)hi + khkεi (h).
j=1
∂xj
3.2 Opérations sur les dérivées, fonctions C 1

Les propriétés suivantes se prouvent de façon élémentaire :
1. Si f est constante, alors dfx est l’application nulle en tout point x. Si
f est linéaire alors dfx = f en tout point.
2. Une combinaison linéaire de deux (ou d’un nombre fini de) fonctions
différentiables en un point est différentiable, et on a d(αf + βg)x =
αdfx + βdgx , ce qui s’écrit aussi ∇(αf + βg) = α∇f + β∇g ou encore
∂(αf +βg) ∂f ∂g
∂xi = α ∂x i
+ β ∂xi
.
12
3. Si f et g sont deux fonctions à valeur réelle différentiables au point x,
alors f g l’est aussi et on a la formule de Leibniz ∇(f g) = f ∇g + g∇f .
4. Si f est différentiable au point x et g est différentiable au point f (x)

alors g ◦ f est différentiable en x avec d(g ◦ f )x = dgf (x) ◦ dfx .
La derniere propriété signifie que la matrice jacobienne de g ◦ f en x est
obtenu en effectuant le produit de la matrice jacobienne de g en f (x) par la
matrice jacobienne de f en x (on rappelle que le produit des matrices n’est
pas commutatif). Dans le cas où f : IRn → IR et g : IR → IR, on obtient
∇(g ◦ f )(x) = g 0 (f (x))∇f (x). Un exemple d’application est le suivant : si
f est différentiable et non nulle au point x, g = 1/f l’est aussi et
∇f (x)
∇g(x) = − .
f 2 (x)
Un cas plus délicat est celui où f : IRm → IRn et g : IRn → R. Le théorème
de composition donne alors ∇(g ◦ f )(x) = (dfx )T ∇g(f (x)). Ceci se traduit
plus intuitivement par la règle
n
∂g ◦ f X ∂fj ∂g
(x) = (x) (f (x)).
∂xi j=1
∂xi ∂xj
Une application de la propriété de composition est la suivante : si f est

une fonction de IRn dans IR qui est homogène de degré α, c’est à dire telle
que f (λx) = λα f (x) pour tout x ∈ IRn et λ ≥ 0, et si f est différentiable au
point x = (x1 , · · · , xn ), on a alors l’identité d’Euler
n
X ∂f
xi (x) = αf (x).
i=1
∂xi
Les résultats qui suivent généralisent des propriétés bien connue pour les
fonctions d’une variable.
Théorème 3.1 Si une fonction f à valeur réelle est différentiable en tout

point du segment S d’extremités x et y, on a
|f (x) − f (y)| ≤ Ckx − yk,
avec C = supz∈S ||∇f (z)||. Par conséquent si f est différentiable sur un ou-
vert convexe U et si C = supz∈U ||∇f (z)|| < +∞, alors f est C-lipschitzienne
sur U .
13
Dans le cas des fonctions à valeur vectorielles, on a le même résultat en
rempacant |f (x) − f (y)| par kf (x) − f (y)k et ||∇f (z)|| par |||dfz |||.
Théorème 3.2 Si U est un ouvert connexe et ∇f (x) est nul en tout point
x ∈ U alors f est constante sur U .
Dans le cas des fonctions à valeur vectorielles, on a le même résultat

en rempacant ∇f (x) par dfx . La nécessité de supposer U connexe dans ce
dernier théorème vient du fait que sinon, on peut le séparer en deux partie
disjointes U ∩ U1 et U ∩ U2 et considérer une fonction f qui prend deux
valeurs constantes différentes sur chaque partie.
Théorème 3.3 Si f est une fonction définie sur IRn ou plus généralement
∂f
sur un convexe et à valeur dans IR, et si la i-ème dérivée partielle ∂xi
est
partout nulle, alors f (x) est une fonction des variables xj pour j 6= i, i.e.
ne dépend pas de i.
On a vu qu’une fonction dont les dérivées partielles sont partout définies

n’est pas nécessairement différentiable. La situation est différente lorsque
ces dérivées partielles sont continues.
Définition 3.4 Une fonction f définie sur un ouvert U de IRn et à valeur

∂f
dans IR est dite de classe C 1 sur U si et seulement ses dérivées partielles ∂xi
sont des fonctions continues sur U . On note C 1 (U ) l’ensemble des fonctions
de classe C 1 sur U .
Cette définition se généralise aux fonctions vectorielles en écrivant que

toutes les composantes f1 , · · · , fm ont cette propriété.
Théorème 3.4 Une fonction C 1 sur un ouvert U de IRn est différentiable

en tout point de U et sa différentielle est une fonction continue.
Définition 3.5 Une fonction f de classe C 1 définie sur un ouvert U de IRn

à valeur dans IRm est un C 1 -difféomorphisme si et seulement si il existe
un ouvert V de IRm tel que f est bijective de U dans V , et si la bijection
réciproque f −1 est aussi de classe C 1 .
Un C 1 -difféomorphisme est parfois appelé changement de variable. Le

théorème de composition montre que pour tout x ∈ U on a alors
dfx ◦ dff−1 −1
(x) = I et dff (x) ◦ dfx = I,
14
où I est l’application identité. Conséquence : dfx est inversible en tout point
x et
dff−1 −1
(x) = (dfx ) .
Ceci nous montre aussi qu’on a nécessairement m = n : un C 1 -difféomorphisme

ne peut exister qu’entre espaces de même dimension.
Définition 3.6 Pour une fonction f de classe C 1 définie sur un ouvert U

de IRn à valeur dans IRn , on appelle “jacobien de f au point x”, la valeur ab-
solue du déterminant de la matrice jacobienne dfx , et on note cette quantité
Jf (x).
Dans le cas d’un C 1 -difféomorphisme, Jf (x) ne s’annulle pas puisque dfx

est inversible en tout point x
3.3 Dérivées partielles d’ordre supérieur, fonctions C k

Définition 3.7 Une fonction f définie sur un ouvert U ⊂ IRn et à valeur
dans IR est de classe C 2 si et seulement si elle est de classe C 1 et si ses
dérivées partielles sont aussi de classe C 1 . On note
∂f
∂2f ∂ ∂xj
=
∂xi ∂xj ∂xi
c’est à dire la i-ème dérivée partielle de la j-ème dérivée partielle de f .

2
Lorsque i = j on utilise la notation simplifiée ∂∂xf2 .
i
Un résultat fondamental est la relation de Schwarz qui affirme que l’ordre

des variables peut-être échangé dans le calcul des dérivées partielles sans
altérer le résultat.
Théorème 3.5 Si f est de classe C 2 sur un ouvert U ⊂ IRn , on a en tout

point de U
∂2f ∂2f
= .
∂xi ∂xj ∂xj ∂xi
Ce résultat nous montre que la matrice des dérivées partielles d’ordre 2

au point x
∂2f
d2 fx = ( (x))i,j=1,···,n
∂xi ∂xj
15
est une matrice symmétrique. On l’appelle la matrice hessienne de f au
point x.
Plus généralement, si les propriétés de f le permettent, on introduit par
récurrence les dérivées partielles d’ordre k,
k−1
∂kf ∂ ∂xi ∂ ···∂x

f
i1
k−1
= ,
∂xik · · · ∂xi1 ∂xik
où i1 , · · · , ik sont des indices choisis dans {1, · · · , n} qui correspondent aux
variables par rapports auxquelles on dérive successivement.
Définition 3.8 Une fonction f définie sur un ouvert U ⊂ IRn et à valeur

dans IR est de classe C k si et seulement si toutes ses dérivées partielles
d’ordre au plus k sont définies et continues sur U . On dit que f est C ∞ si
et seulement si elle est de classe C k pour tout k > 0.
Cette définition se généralise aux fonctions vectorielles en écrivant que

toutes les composantes f1 , · · · , fm ont cette propriété.
En appliquant itérativement la relation de Schwarz, on voit que l’ordre
d’apparition des variables xi dans la dérivation peut être permuté. On
emploie donc la notation
∂kf
k1
,
∂x1 · · · ∂xknn
pour désigner la dérivée partielle d’ordre k faisant intervenir ki fois la vari-
able xi avec k1 + · · · + kn = k.
Les propriétés suivantes se démontrent par récurrence sur k :
1. Une combinaison linéaire de deux (ou d’un nombre fini de) fonctions
C k est aussi de classe C k .
2. Si f et g sont deux fonctions à valeur réelle et de classe C k , alors f g

l’est aussi.
3. Si f est de classe C k d’un ouvert U ⊂ IRn dans IRm et si g est de classe

C k d’un ouvert V ⊂ IRm contenant f (U ) dans IRp , alors g ◦ f est de
classe C k de U dans IRp .
4. Si f est un C 1 -difféomophisme et est de classe C k , la réciproque f −1

est aussi de classe C k .
16
3.4 Points critiques et extremas
Définition 3.9 Soit f une fonction C 1 définie sur un ouvert U ⊂ IRn et à
valeur réelle. On dit que x ∈ U est un point critique de f si et seulement si
∇f (x) = 0.
Cette définition s’étend aux fonctions à valeur vectorielle en remplacant

∇f (x) par dfx .
Définition 3.10 Soit f une fonction définie sur un ensemble E ⊂ IRn et à

valeur réelle. On dit que f admet un minimum global au point x ∈ E si et
seulement si f (x) ≤ f (y) pour tout y ∈ E. On dit que f admet un minimum
local au point x ∈ E si et seulement si il existe un ouvert U contenant x tel
que f (x) ≤ f (y) pour tout y ∈ U ∩ E.
On définit de la même maniere la notion de maximum global ou local.

On utilise la dénomination d’extremum pour désigner sans distinction un
maximum ou minimum. On a déjà signalé l’existence d’un maximum et
minimum global sur un compact. Par ailleurs si f admet un extremum local
en x, alors pour tout u ∈ IRn la fonction d’une variable gu (t) = f (x + tu)
admet un extremum local en t = 0 et par conséquent g 0 (0) = ∇f (x) · u = 0.
Comme u est arbitraire, ceci entraine le résultat suivant.
Théorème 3.6 Soit f une fonction C 1 définie sur un ouvert U ⊂ IRn et à

valeur dans IR. Si f admet un extremum local en un point x ∈ U alors x
est un point critique de f .
Etre un point critique est donc une condition nécessaire pour être un
extremum, elle n’est cependant pas suffisante : considérer par exemple la
fonction f (x, y) = x2 − y 2 au point (0, 0).
Afin d’identifier plus précisément les extrema parmis les points critiques,
on est amené à introduire un développement limité au deuxième ordre qui
fait intervenir la matrice hessienne d2 fx .

valeur réelle. Pour x ∈ U et pour tout h ∈ IRn suffisament petit, on a
1
f (x + h) = f (x) + ∇f (x) · h + d2 fx h · h + khk2 ε(h),
2
où limh→0 ε(h) = 0.
17
Le terme du deuxième ordre s’écrit plus précisément en fonction des
dérivées partielles
n X
n
X ∂2f
d2 fx h · h = (x)hi hj .
i=1 j=1
∂xi ∂xj
C’est la forme quadratique associée à la matrice d2 fx .

On rappelle qu’une matrice symmétrique se diagonalise dans une base or-
thonormée et que toutes ses valeurs propres sont réelles. En notant (f1 , · · · , fn )
et (λ1 , · · · , λn ) la base orthonormée et les valeurs propres pour la matrice
d2 fx , et en décomposant h suivant cette base selon h = ni=1 ai fi on obtient
P
l’expression plus simple

n
X
d2 fx h · h = λi a2i .
i=1
On peut à partir de cette remarque obtenir une condition nécessaire pour

avoir un extremum au point x qui précise celle de point critique énoncée
précedemment.

valeur réelle. Si f admet un minimum local au point x alors les λi sont tous
positifs ou nuls. Si f admet un maximum local en x alors les alors les λi
sont tous négatifs ou nuls.
On peut aussi obtenir une condition suffisante pour avoir un extremum

au point x.
Théorème 3.9 Soit f une fonction C 2 définie sur un ouvert U ⊂ IRn et

à valeur réelle. Soit x ∈ U un point critique de f . Si les λi sont tous
strictement positifs alors f admet un minimum local au point x. Si les λi
sont tous strictement négatifs alors f admet un maximum local au point x
Dans le cas n = 2 la matrice hessienne est de taille 2 × 2 et il est partic-

ulierement facile de vérifier le signe des λi en remarquant que le déterminant
Jf (x) = det(d2 fx ) vaut λ1 λ2 et que la trace tr(d2 fx ) vaut λ1 + λ2 . On voit
donc que si x est un point critique et si Jf (x) > 0, alors f admet un ex-
tremum au point x. On voit aussi que cet extremum est un maximum si
tr(d2 fx ) < 0 et un minimum si tr(d2 fx ) > 0.
18
3.5 Formes différentielles
Considérons une fonction f de classe C 1 de IRn dans IR. On sait que
n
X ∂f
dfx (h) = ∇f (x) · h = (x)hi .
i=1
∂xi
Par conséquent, en notant dxi l’application coordonnée dxi (h) = hi , on peut

écrire n
X ∂f
dfx = (x)dxi
i=1
∂x i
Cette notation a une autre interpretation plus intuitive où dfx représente
la variation au premier ordre de f lorsque l’on effectue un déplacement
infinitésimal de (dx1 , · · · , dxn ) autour du point x.
De façon plus générale, étant donné un champ de vecteur g = (g1 , · · · , gn )
de classe C 1 de IRn dans IRn , on définit la 1-forme différentielle ω associée à
ce champ de vecteur par la formule
ω(x) = g1 (x)dx1 + · · · + gn (x)dxn .
D’un point de vue rigoureux ω est une application de IRn dans l’ensemble
des formes linéaires sur IRn . Intuitivement, on peut à nouveau penser à ω(x)
comme une combinaison linéaire de variations infinitésimales (dx1 , · · · , dxn ).
Nous utiliserons ultérieurement cette notion dans la définition des intégrales
curvilignes.
Définition 3.11 La forme différentielle ω est dite exacte si et seulement

si il existe une fonction f telle que ω(x) = dfx pour tout x, c’est à dire
∂f
gi (x) = ∂x i
(x) pour tout i et pour tout x.
Définition 3.12 La forme différentielle ω est dite fermée si et seulement

∂g ∂gi
si pour tout i, j dans {1, · · · , n} on a ∂xji = ∂xj
.
D’après la relation de Schwarz, on voit que toute forme différentielle

exacte est fermée. Un résultat remarquable est que la réciproque est vraie.
Théorème 3.10 Une forme différentielle ω définie sur IRn est exacte si et
seulement si elle est fermée.
19
4 Résultats plus avancés
4.1 Equations aux dérivées partielles
De nombreux phénomènes physiques se modélisent par des équations faisant
intervenir les dérivées partielles de la fonction solution de l’équation. De
telles équations sont appelées équations aux dérivées partielles et souvent
notées EDP.
L’exemple le plus simple d’EDP est l’équation de transport: étant donné
un réel c 6= 0, on cherche u(x, t) de classe C 1 de IR2 dans IR qui vérifie
∂u ∂u
+c = 0.
∂t ∂x
On peut faire le changement de variable y = x − ct et z = x + ct qui est un
C 1 difféomorphisme de IR2 et on écrit donc u(x, t) = v(y, z). On voit alors
que l’équation se re-écrit
∂v
= 0,
∂z
Ce qui montre que u est uniquement fonction de y = x−ct. On dit aussi que
u est constante le long des droites caractéristiques d’équation x = x0 + ct,
pour tout x0 ∈ IR.
En imposant une condition initiale au temps t = 0,
u(x, 0) = u0 (x),
où u0 est une fonction fixée, on voit ainsi que la solution de l’équation
transport s’écrit
u(x, t) = u0 (x − ct),
ce qui correspond au “déplacement” de la fonction u0 à la vitesse c le long
de l’axe des x.
Voici d’autres exemples d’équations aux dérivées partielles:
1. L’équation de la chaleur: étant donné a > 0 et u0 , on cherche u(x, t)
de classe C 2 de IR × IR+ dans IR qui vérifie
∂u ∂2u
− a 2 = 0, t > 0 et u(x, 0) = u0 (x).
∂t ∂x
On pourra essayer de montrer que la solution est donnée par la formule
y2
Z
−1/2
u(x, t) = (2πat) u0 (x − y)e− 2at dy.
IR
20
2. L’équation des ondes: étant donné c ∈ IR on cherche u(x, t) de classe
C 2 de IR2 dans IR qui vérifie
∂2u 2
2∂ u
− c = 0.
∂t2 ∂x2
On pourra à nouveau utiliser le changement de variable (x, t) → (y, z)
pour décrire la solution à partir des conditions initiales u0 (x) = u(x, 0)
et u1 (x) = ∂u
∂t (x, 0).
3. L’équation de Burgers: on cherche u(x, t) de classe C 1 de IR × IR+ dans

IR qui vérifie
∂u ∂u
+u = 0, t > 0 et u(x, 0) = u0 (x).
∂t ∂x
Il s’agit d’une équation non-linéaire à cause de la présence du facteur u
qui signifie que la vitesse est proportionelle à l’amplitude de la solution.
On pourra essayer de montrer que u est constante le long des droites
caractéristiques d’équation x = x0 + u0 (x0 )t pour tout x0 ∈ IR, et de
comprendre pourquoi les solutions de cette équation ne peuvent pas
rester de classe C 1 (ni même continues) au bout d’un temps critique
T > 0.
4.2 Inversion locale et fonctions implicites

Un résultat classique pour les fonctions d’une variable affirme que si f est C 1
et f 0 (x) 6= 0, alors f est strictement monotone sur un intervalle I contenant
x et est une bijection entre les intervalles I et f (I). De plus la bijection
réciproque est aussi de classe C 1 .
Afin d’établir un résultat du même type pour les fonctions de plusieurs
variables, on utilise un résultat simple mais fondamental : le théorème du
point fixe de Picard. Rappelons que x est un point fixe d’une fonction f si
et seulement si f (x) = x.
Théorème 4.1 Soit E ⊂ IRn et f : E → E une fonction C-Lipschitzienne

avec C < 1. Alors f admet un unique point fixe x ∈ E.
Ce théorème se démontre en remarquant que toute suite (xn )n≥0 définie

par xn = f (xn−1 ) avec x0 ∈ E vérifie kxn+1 − xn k ≤ C n kx1 − x0 k ce qui
entraine qu’elle est de Cauchy. Sa limite est nécessairement un point fixe de
f et son unicité découle de la propriéé kf (x) − f (y)k ≤ Ckx − yk.
21
Soit à présent f de classe C 1 sur un ouvert U ⊂ IRn et à valeur dans
n
IR et telle que dfx est inversible pour un point x ∈ U . On peut étudier
l’inversibilité de f au voisinage de x en considérant z ∈ IRn tel que kz −
f (x)k ≤ ε, et en cherchant une solution y = x + h à l’équation z = f (y).
En utilisant la différentiabilité de f en x, on peut re-écrire cette équation
suivant
z = f (x) + dfx (h) + φ(h),
où φ est une fonction C 1 telle que φ(0) = 0 et dφ0 = 0. La solution h est
solution du problème de point fixe
h = dfx−1 (f (x) − z) + dfx−1 (φ(h)) = F (h).
On vérifie alors que pour ε suffisament petit, il existe α > 0 tel que F
est C-Lipschitzienne de B(0, α) dans B(0, α) avec C < 1. En appliquant
théorème du point fixe, on en déduit l’existence d’un unique antécédent y
de z dans B(x, α). Ce raisonnement nous conduit au théorème d’inversion
locale énoncé ci-dessous.
Théorème 4.2 Soit f de classe C 1 sur un ouvert U ⊂ IRn et à valeur

dans IRn et soit x ∈ U telle que dfx est inversible. Alors il existe deux
ouverts, U 0 ⊂ U et V tels que x ∈ U 0 et f (x) ∈ V et tels que f est un C 1
difféomorphisme de U 0 dans V .
Il existe une version de ce théorème qui permet d’établir la propriété

de C 1 difféomorphisme sur la totalité du domaine U . C’est le théorème
d’inversion globale.
Théorème 4.3 Soit f de classe C 1 sur un ouvert U ⊂ IRn et à valeur

dans IRn , et telle que f est injective sur U . Alors f (U ) est un ouvert V de
IRn . Dans ces conditions f est un C 1 difféomorphisme de U dans V si et
seulement si dfx est inversible en tout point x de U .
Une application du théorème d’inversion locale concerne le problème

suivant : pour f fonction C 1 de deux variables, on considère l’équation
f (x, y) = 0 et on cherche à comprendre si cette équation est en un certain
sens équivalente à y = g(x) où g est une fonction d’une variable. L’exemple
de la fonction f (x, y) = x2 + y 2 − 1 nous montre que ceci n’est possible que
localement : certaines portions du cercle unité s’identifient au graphe de la
fonction y = (1 − x2 )1/2 , d’autres à celui de la fonction y = −(1 − x2 )1/2 .
22
D’autre portions, telles qu’au voisinage du point (1, 0) ne peuvent s’identifier
à un graphe. Le théorème des fonctions implicite donne un résultat général
allant dans ce sens.
Théorème 4.4 Soit f une fonction C 1 de IR2 dans IR et (x0 , y0 ) un couple

tel que f (x0 , y0 ) = 0 et tel que ∂f
∂y (x0 , y0 ) 6= 0. Alors il existe deux intervalles
ouverts I et J tels que x0 ∈ I et y0 ∈ J, et une fonction g de classe C 1 de I
dans J telle que pour tout (x, y) ∈ I × J on ait
f (x, y) = 0 ⇔ y = g(x).
De plus on a
∂f
(x, g(x))
g 0 (x) = − ∂f
∂x
.
∂y (x, g(x))
La preuve de ce théorème peut s’effectuer en considérant la fonction

h(x, y) = (x, f (x, y)) et en lui appliquant le théorème d’inversion locale qui
permet de définir (x, g(x)) comme l’antécédent de (x, 0).
Le théorème des fonctions implicite se généralise aux fonctions de plus
de deux variables : si f est C 1 de IRn × IRm dans IRm avec f (x0 , y0 ) = 0 et si
sa différentielle par rapport à la variable y est inversible en (x0 , y0 ), alors il
existe deux ouverts U et V tels que x0 ∈ U et y0 ∈ V , et une fonction g de
classe C 1 de U dans V telle que pour tout (x, y) ∈ U × V on ait f (x, y) = 0
si et seulement si y = g(x).
4.3 Extremas liés

On s’est déjà intéressé à l’étude des extremas d’une fonction f définie sur
un ouvert de IRn . Dans de nombreuses applications, il est fréquent que
l’on cherche à minimiser ou maximiser f sur un sous-ensemble restreint ce
qui correspond à imposer une contrainte sur la variable x. Par exemple,
minimiser une fonction f sur la sphère unité signifie que l’on impose sur x
la contrainte x21 + · · · + x2n − 1 = 0.
De façon générale, étant données deux fonctions f et g de classe C 1
définies sur un ouvert U de IRn et à valeur dans IR, on considère l’ensemble
E := {x ∈ U ; g(x) = 0},
et on s’intéresse aux extrémas de f sur M . Le théorème des extrémas

liés permet de préciser les propriétés de tels extrémas et d’en faciliter la
recherche.
23
Théorème 4.5 Soit x ∈ E tel que ∇g(x) 6= 0 et tel que f restreinte à E
admet un extremum local en x. Alors, il existe un réel λ appelé multiplicateur
de Lagrange tel que
∇f (x) = λ∇g(x).
Esquissons la preuve: pour u ∈ IRn et v ∈ IRn tel que ∇g(x) · v 6= 0, on

définit une fonction de IR2 dans IR
F (s, t) = g(x + su + tv).
On s’intéresse aux (s, t) tels que F (s, t) = 0 c’est à dire x + su + tv ∈ E.

Le théorème des fonctions implicites montre qu’il existe une fonction h telle
que t = h(s). On voit que si x est un extremum local de f sur E, alors la
fonction d’une variable G(s) = f (x + su + h(s)v) admet un extremum local
en s = 0. On a donc
∇g(x) · u
0 = G0 (0) = ∇f (x) · u + h0 (0)∇f (x) · v = ∇f (x) · u − ∇f (x) · v,
∇g(x) · v
Ce qui montre que ∇f (x) · u = λ∇g(x) · u pour tout u avec λ = ∇f (x)·v

∇g(x)·v , et
donc ∇f (x) = λ∇g(x).
Voici finalement une généralisation dans le cas où l’ensemble E est défini
par plusieurs contraintes associées à des fonctions (g1 , · · · , gm ), c’est à dire
E := {x ∈ U ; g1 (x) = · · · = gm (x) = 0},
Théorème 4.6 Soit x ∈ E tel que les vecteurs (∇g1 (x), · · · , ∇gn (x)) sont
indépendants et tel que f restreinte à E admet un extremum local en x.
Alors, il existe des réels λi appelés multiplicateurs de Lagrange tel que
∇f (x) = λ1 ∇g1 (x) + · · · + λn ∇gn (x).
5 Construction de l’intégrale multiple

5.1 Intégrale sur un rectangle
Pour une fonction f d’une variable définie sur un intervalle [a, b], rappelons
que l’intégrale ab f (t)dt est définie par un passage à la limite des sommes de
R
Riemann
m−1
X
Σ= (xk+1 − xk )f (yk ),
k=0
24
(avec a = x0 < x1 < · · · < xm−1 < xm = b et yk ∈ [xk , xk+1 ]) lorsque
maxk |xk+1 − xk | tend vers 0.
Il est aisé de reproduire cette construction pour une fonction f de n
variables à valeur dans IR lorsque le domaine d’intégration est de forme
rectangulaire, c’est à dire
R = [a1 , b1 ] × · · · × [an , bn ].
Par abus de notation, on appellera “rectangle” un tel domaine, quelquesoit

la dimension n considérée. On notera
|R| = vol(R) = (b1 − a1 ) · · · (bn − an ),
le volume de R.
On peut dans ce cas définir une subdivision de R à partir de subdivisions
de chaque intervalle [ai , bi ]: si pour tout i, on a la subdivision
ai = x0,i < x1,i · · · xmi −1,i < xmi ,1 = bi ,
alors l’ensemble des rectangles
[xk1 ,1 , xk1 +1,1 ] × · · · × [xkn ,n , xkn+1 ,n ],
lorsque (k1 , · · · , kn ) parcourt {0, · · · , m1 − 1} × · · · × {0, · · · , mn − 1} forme

une partition de R (à ceci près que certains de ces rectangles ont des faces
communes). Afin de simplifier, on notera
S = (Rj )j=1,···,M ,
où M = m1 m2 · · · mn est le nombre total de rectangles présent dans S.

Appelons S la subdivision de R constituée par cet ensemble de rectangle
et désignons par
δ(S) := max{|xki +1,i − xki ,i | ; ki = 0, · · · , mi − 1, i = 1, · · · , n}
le pas de la subdivision qui correspond au plus grand côté sur l’ensemble des
rectangles Rj de S.
Comme dans le cas unidimensionel, on dit que S2 est plus fine que S1 si
tout rectangle de S1 est l’union de plusieurs rectangles de S2 . Si S1 et S2
sont deux subdivisions de R alors on peut construire une subdivision S3 =
S3 (S1 , S2 ) plus fine que S1 et S2 en considérant pour chaque i = 1, · · · , n
l’union des subdivisions d’une variable associées à S1 et S2 .
25
Pour S = (Rj )j=1,···,M une subdivision de R et yj ∈ Rj , on définit la
somme de Riemann
M
X
Σ = Σ(S) = f (yj )|Rj |.
j=1
Remarquons que Σ dépend aussi du choix des yj dans chaque rectangle Rj .

On peut par exemple décider de prendre pour yj le centre de Rj .
Définition 5.1 Une fonction f définie sur R est dite intégrable (au sens de
Riemann) sur R si et seulement si pour tout ε > 0, il existe α > 0 tel que
pour toutes subdivisions S1 et S2 de pas inférieur à α, on a
|Σ(S1 ) − Σ(S2 )| ≤ ε,
pour toutes les sommes de Riemann associées à S1 et S2 . En conséquence

il existe un nombre I, tel que
|I − Σ(S)| ≤ ε.
pour tout S tel que δ(S) ≤ α, c’est à dire I = limδ(S)=0 Σ(S). On dit que I
est l’intégrale de f sur R et on la note
Z
I= f (x)dx.
R
Théorème 5.1 Si f est une fonction continue sur R, alors elle est intégrable
sur R.
La preuve de ce théorème s’effectue comme dans le cas des fonctions

d’une variable en utilisant la subdivision S3 (S1 , S2 ) et la continuité uniforme
de f sur le compact R, qui permet de prouver que |Σ(Si ) − Σ(S3 )| ≤ ε/2
pour i = 1, 2 et α > 0 suffisament petit.
On utilise parfois une autre notation pour l’intégrale multiple :
Z Z
I= f (x1 , x2 )dx1 dx2 ,
R
dans le cas des intégrales à 2 variables,

Z Z Z
I= f (x1 , x2 , x3 )dx1 dx2 dx3 ,
R
dans le cas des intégrales à 3 variables, etc. Cette notation sera justifie par
une technique de calcul des intégrales multiples donnée en section 6.
26
5.2 Intégrale sur un domaine quarrable
On souhaite à présent étendre la notion d’intégrale multiple à des domaines
plus généraux que des rectangles. Afin de réaliser cet objectif, nous al-
lons nous limiter à une classe restreinte, néammoins suffisament riche, de
domaines.
Soit E un ensemble borné de IRn . Cela entraine en particulier que E est
contenu dans un rectangle suffisament grand:
E ⊂ R = [a1 , b1 ] × · · · × [an , bn ].
Considérons à présent une subdivision S de R. On lui associe le pavage

extérieur et intérieur de E par S qui sont les sous-ensembles de S définis
par
Sext (E) = {Rj ∈ S ; Rj ∩ E 6= ∅},
et ◦
Sint (E) = {Rj ∈ S ; Rj ⊂E }.
On a clairement
∪Rj ∈Sint (E) Rj ⊂ E ⊂ ∪Rj ∈Sext (E) Rj
On remarque aussi que la différence entre le pavage extérieur et intérieur

correspond aux rectangles qui touchent la frontiere ∂E de E :
Sext (E) \ Sint (E) = {Rj ∈ S ; Rj ∩ ∂E 6= ∅}.
Définition 5.2 Un ensemble E est dit quarrable si et seulement si la pro-

priété suivante est vérifiée: pour tout ε > 0, il existe α > 0 tel que pour tout
S tel que δ(S) ≤ α on a
X X
|Rj | − |Rj | ≤ ε.
Rj ∈Sext (E) Rj ∈Sint (E)
Ceci équivaut à l’existence d’un un nombre A tel que

X X
lim |Rj | = lim |Rj | = A.
δ(S)→0 δ(S)→0
Rj ∈Sext (E) Rj ∈Sint (E)
Le nombre A est appelé volume ou mesure de E (ou aire en dimension

n = 2) et noté
A = |E| = vol(E).
27
Cette définition nous montre en particulier que E est quarrable si et seule-
ment si sa frontiere est quarrable et est de volume nul:
|∂E| = 0.
Les propriétés suivantes se démontrent assez facilement:

1. Un rectangle R est quarrable et son aire définie ci-dessus est bien égale
au produit de ses côtés.
2. Si E et F sont quarrables alors E ∩ F et E ∪ F le sont aussi et on a
|E ∪ F | + |E ∩ F | = |E| + |F |.
3. Si E et F sont quarrables alors E \ F l’est aussi et on a
|E \ F | + |E ∩ F | = |E|.
4. L’aire est invariante par translation et par rotation (on le démontre

d’abord sur un rectangle puis sur un ensemble général).
On peut à présent définir la notion de somme de Riemman pour une
fonction f définie sur E : étant donné une subdivision S du rectangle R, on
définit X
Σ = Σ(S) = f (yj )|Rj |.
Rj ∈S(E)
où yj ∈ Rj et S(E) est un ensemble vérifiant Sint (E) ⊂ S(E) ⊂ Sext (E).
Définition 5.3 Soit E un ensemble quarrable fermé. Une fonction f définie

sur E est dite intégrable sur E si et seulement si pour tout ε > 0, il existe
α > 0 tel que pour toutes subdivision S1 et S2 de pas inférieur à α, on a
|Σ(S1 ) − Σ(S2 )| ≤ ε,
pour toutes les sommes de Riemann associées à S1 et S2 . En conséquence

il existe un nombre I, tel que
|I − Σ(S)| ≤ ε.
pour tout S tel que δ(S) ≤ α, c’est à dire I = limδ(S)=0 Σ(S). On dit que I
est l’intégrale de f sur E et on la note
Z
I= f (x)dx.
E
28
Théorème 5.2 Si f est une fonction continue sur E, alors elle est intégrable
sur E.
La preuve de ce théorème ressemble à celle dans le cas du rectangle, mais

il faut en plus utiliser le fait que l’ensemble est quarrable. Les propriétés
suivantes se démontrent facilement:
R
1. L’intégrale de la constante 1 donne l’aire : E 1dx = |E|
2. L’intégrale d’une fonction sur un domaine de volume nul est toujours

nulle.
R R R
3. Linéarité de l’intégrale : E (αf +βg)(x)dx =α E f (x)dx+β E g(x)dx.
4. Additivité par rapport

R
aux ensembles
R
: si E et
R
F sont quarrables et
|E ∩ F | = 0 alors E∪F f (x)dx = E f (x)dx + F f (x)dx.
R R
5. Si f et |f | sont intégrables, on a | E f (x)dx| ≤ E |f (x)|dx.
f (x)g(x)dx|2 ≤ ( |f (x)|2 dx)( |g(x)|2 dx).

R R R
6. Cauchy-Schwarz : | E E E
R R
7. Monotonie : si f ≤ g sur E on a E f (x)dx ≤ E g(x)dx.
De la derniere propriété, on déduit que si f est continue sur E on a

Z
|E| min f (x) ≤ f (x)dx ≤ |E| max f (x),
x∈E E x∈E
ce qui entraine la formule de la moyenne en supposant l’ensemble E connexe

et d’aire non-nulle : il existe x ∈ E tel que
1
Z
f (x) = f (x)dx.
|E| E
Comme dans le cas des fonctions d’une variable il est possible d’étendre
la notion d’intégrale à des classes de fonctions plus générales que les fonctions
continues, en particulier les fonctions continues par morceaux.
Théorème 5.3 Si un ensemble quarrable E peut s’écrire
E = ∪m
i=1 Ei ,
ou chaque Ei est quarrable et Ei ∩ Ej = ∅, et si f est une fonction uni-

formément continue sur chaque partie Ei , alors f est intégrable sur E.
29
On voit en particulier que si f est continue sur un fermé quarrable E et
que l’on définit l’extension de f sur le rectangle R par
f˜(x) = f (x) si x ∈ E, f˜(x) = 0, si x ∈

/ E,
alors on a Z Z
f˜(x)dx = f (x)dx.
R E
En particulier, si on définit la fonction indicatrice 1E de E par
1E (x) = 1 si x ∈ E, 1E (x) = 0, si x ∈
/ E,
on a Z
1E (x)dx = |E|.
R
On peut remarquer au passage qu’un ensemble est quarrable si et seulement
si son indicatrice est intégrable.
Il est temps de donner une classe importante d’ensembles quarrables.
Théorème 5.4 Soit E ⊂ IRn un ensemble dont la frontiere ∂E est contenue

dans une union finie de boules Bi = B(xi , ri ), i = 1, · · · , m, avec la propriété
suivante : pour chaque i, il existe un repere orthonormé de coordonnées
(α1 , · · · , αn ) centré autour du point xi et une fonction continue fi de IRn−1
dans IR , tels que ∂E ∩ Bi s’identifie au graphe de fi :
∂E ∩ Bi = {(α1 , · · · , αn ) ; αn = fi (α1 , · · · , αn−1 ) et α12 + · · · + αn−1

2
≤ ri2 }.
Alors l’ensemble E est quarrable.
Ce théorème se démontre en prouvant que ∂E est quarrable avec |∂E| =

0. Entre autres conséquences, on voit qu’un polygone (ensemble dont la
frontiere est une ligne brisée fermée) est quarrable, que l’ensemble des points
compris entre les graphes de deux fonctions continues f ≤ g restreintes à un
ensemble quarrable est quarrable, etc.
Il est aussi facile d’exhiber des ensembles non-quarrables : par exemple
l’ensemble E des points d’un rectangle R dont les coordonnées sont des
nombres rationels. On voit aussi que la fonction indicatrice de E n’est pas
intégrable au sens de Riemann sur R (elle l’est avec une définition plus fine
de l’intégrale, due à Lebesgue).
30
6 Calcul des intégrales multiples
6.1 Le théorème de Fubini
Il est bien connu que le calcul des intégrales des fonctions d’une variable est
lié au calcul des primitives : si f est telle que f = F 0 , on a
Z b
f (t)dt = F (b) − F (a) := [F (x)]ba
a
Le calcul pratique des intégrales multiple se heurte au fait que la notion

de primitive est uniquement bien définie pour les fonctions d’une variable.
Le théorème de Fubini offre un moyen de ramener le calcul des intégrales
multiples à celui des intégrales des fonctions d’une variable. Une version
élémentaire de ce théorème s’énonce pour les fonctions de deux variables
définies sur un rectangle.
ThéorèmeR 6.1 Soit f une fonction continue de R = [a, b] × [c, d] dans IR.
Si F (x) = cd f (x, y)dy alors F est continue et on a
Z Z Z b
f (x, y)dx dy = F (x)dx.
R a
Ce résultat peut aussi s’écrire sous la forme

Z Z Z b Z d
f (x, y)dx dy = ( f (x, y)dy)dx.
R a c
En échangeant le rôle des variables x et y, on voit que l’on a aussi

Z Z Z d Z b
f (x, y)dx dy = ( f (x, y)dx)dy,
R c a
autrement dit l’ordre des variables d’intégration n’importe pas.

Il existe de multiples formes plus générales du théorème de Fubini. On
peut tout d’abord considérer des fonctions de n variables continues sur un
rectangle R = [a1 , b1 ] × · · · × [an , bn ]. On a alors
Z Z b1 Z b2 Z bn−1 Z bn
f (x)dx = ( ··· ( f (x)dxn )dxn−1 · · · dx2 )dx1 ,
R a1 a2 an−1 an
et l’ordre d’intégration peut être arbitrairement permuté sans changer le

résultat.
31
On peut ensuite considérer des fonctions de n variables continues sur un
ensemble quarrable E contenu dans R = [a1 , b1 ] × · · · × [an , bn ]. Divisons
l’ensemble des variables en deux parties {1, · · · , m} et {m + 1, · · · , n}. et no-
tons x = (x1 , · · · , xn ) := (y, z) avec y = (x1 , · · · , xm ) et z = (xm+1 , · · · , xn ),
et f (y, z) := f (x). On note T = [a1 , b1 ] × · · · × [am , bm ], et
Ey := {z ∈ IRn−m (y, z) ∈ E},
la section de l’ensemble E dont les m premieres coordonnées sont données
par y. En supposant que les ensembles Ey sont quarrables, on peut alors
montrer que la fonction F définie par
Z
F (y) = f (y, z)dz
Ey
est intégrable au sens de Riemann, et que l’on a

Z Z Z
f (x)dx = ( f (y, z)dz)dy.
E T Ey
Par exemple, dans le cas des fonctions de 3 variables et avec m = 1, cela

revient à “couper le domaine en tranches” et à calculer sur chaque tranche
une intégrale de deux variables que l’on integre ensuite par rapport à la
variable restante. A nouveau, il est possible de changer l’ordre des variables
dans l’intégration. Il est intéressant de noter que pour un domaine général,
la continuité de f n’entraine pas nécessairement celle de F (y) : considérer
par exemple une fonction constante sur un domaine en forme de L consitué
de 3 carrés adjacents.
Enfin, il est possible d’appliquer le théorème de Fubini à des classes de
fonctions plus générales que les fonctions continues, par exemples des fonc-
tions continues par morceaux du type considéré dans la section précèdente.
6.2 Changements de variables

Soit φ un C 1 difféomorphisme envoyant un ensemble F sur un E = φ(F ).
On vérifie aisément que F est fermé (resp. ouvert) si et seulement si E
est fermé (resp. ouvert). On peut aussi montrer que E est quarrable si et
seulement si F est quarrable. Si f est une fonction intégrable sur E, on va
chercher à relier l’intégrale de f sur E à celle de f ◦ φ sur F .
Considérons tout d’abord le cas particulièr où φ est une application
affine, c’est à dire
φ(x) = Ax + b,
avec b ∈ IRn et A une matrice n × n inversible.
32
Théorème 6.2 Lorsque φ est une application affine du type ci-dessus, on
a Z Z
f (x)dx = |det(A)| (f ◦ φ)(y)dy.
E F
La preuve de ce théorème s’effectue en remarquant que si S est une

subdivision en rectangles utilisée pour approcher l’intégrale sur E par la
somme de Riemann X
Σ= f (yj )|Rj |,
Rj ∈S(E)
alors on obtient un pavage S(F ) par les parallellépipedes Tj = φ−1 (Rj ) tels
que |Rj | = |det(A)| |Tj |, et on peut écrire
X
Σ = |det(A)| (f ◦ φ)(xj )|Tj |,
Tj ∈S(F )
avec xj = φ−1 (yj ) ∈ Tj .

Intuitivement l’élément infinitésimal rectangulaire dy est transformé par
φ en un élément infinitésimal dx avec |dx| = |det(A)| |dy|.
Dans le cas où φ est un C 1 -difféomorphisme plus général, lorsque l’on se
place au voisinage de y, φ est proche de l’application affine
φ̃(z) = φ(y) + dφy (z − y),
ce qui signifie (toujours intuitivement) que l’élément infinitésimal rectangu-

laire dy au voisinage de y est transformé par φ en un élément infinitésimal
dx avec |dx| = Jφ (y) |dy|, où Jφ (y) = |det(dφy )| est le Jacobien de φ en y.
Ceci justifie le résultat suivant dont la preuve est difficile.
Théorème 6.3 Lorsque φ est un C 1 difféomorphisme général on a

Z Z
f (x)dx = (f ◦ φ)(y)Jφ (y)dy.
E F
Voici l’application de ce théorème pour trois changements de variables

très classiques :
1. Coordonnées polaires (r, θ) dans IR2 : (x, y) = φ(r, θ) = (r cos θ, r sin θ)

et Jφ (r, θ) = r. Si f (x, y) = g(r, θ) avec g = f ◦ φ, on peut écrire
Z Z
f (x, y)dx dy = g(r, θ)rdr dθ
E F
33
2. Coordonnées cylindriques (r, θ, z) dans IR3 : (x, y, z) = φ(r, θ, z) =
(r cos θ, r sin θ, z) et Jφ (r, θ) = r. Si f (x, y, z) = g(r, θ, z) avec g = f ◦φ,
on peut écrire
Z Z
f (x, y, z)dx dy dz = g(r, θ, z)rdr dθ dz
E F
3. Coordonnées sphèriques (r, θ, ϕ) dans IR3 : (x, y, z) = φ(r, θ, ϕ) =

(r sin θ cos ϕ, r sin θ sin ϕ, r cos θ) et Jφ (r, θ) = r2 | sin θ|. Si f (x, y, z) =
g(r, θ, ϕ) avec g = f ◦ φ, on peut écrire
Z Z
f (x, y, z)dx dy dz = g(r, θ, ϕ)r2 | sin θ|dr dθ dϕ
E F
6.3 Intégrales généralisées

Il est aussi possible d’étendre aux fonctions de plusieurs variables la notion
d’intégrale généralisée.
Définition 6.1 Soit E un ensemble borné ou non. On dit que f est locale-
ment intégrable sur E si et seulement si elle est intégrable sur tout fermé
borné quarrable contenu dans E.
Définition 6.2 Soit E un ensemble borné ou non. Une fonction f est dite
“intégrable sur E” si et seulement si f et |f | sont localement intégrables sur
E et si pour tout ε > 0, il existe un fermé borné F ⊂ E tel que pour tout
fermé borné quarrable G ⊂ E \ F on a
Z
|f (x)|dx ≤ ε.
G
On voit que cette propriété entraine l’existence d’un nombre I tel que
pour tout ε > 0, il existe un fermé borné F ⊂ E tel que pour tout fermé
borné quarrable G vérifiant F ⊂ G ⊂ E on a
Z
|I − f (x)dx| ≤ ε.
G
R
Le nombre I est appelé l’intégrale de f sur E et est noté I = E f (x)dx.
On peut ainsi définir l’intégrale généralisée sur un domaine non-borné un
ensemble fermé privé d’un point, de sa frontière, etc. La définition ci-dessus
équivaut à la propriété qu’il existe une constante C < +∞ telle que
Z
|f (x)|dx ≤ C,
G
34
pour tout ensemble G fermé borné quarrable contenu dans E. Cette dernière
propriété est souvent plus facile à vérifier car il suffit qu’elle soit valable pour
une famille F de fermés bornés telle que tout fermé borné quarrable contenu
dans E soit contenu dans un G ∈ F. Par exemple, pour que f admette une
intégrale généralisée sur IRn , il suffit qu’il existe C < +∞ telle que
Z
|f (x)|dx ≤ C,
B(0,R)
pour tout R > 0. Les propriétés classiques des intégrales énoncées dans les
sections précèdentes s’étendent aux intégrales généralisées. Il est en partic-
ulier possible d’appliquer les techniques de calcul basées sur le changement
de variable et le théorème de Fubini pour prouver l’existence d’intégrales
généralisées et les calculer.
2 2
A titre d’exemple, on peut considérer le calcul de I = IR2 e−(x +y ) dxdy.
R
En se restreignant sur B(0, R) et en effectuant le changement de variable en

coordonnées polaires, on obtient
Z Z 2π Z R Z R
−(x2 +y 2 ) −r2 2 2
e dxdy = e r dr dθ = 2π e−r r dr = π[−e−r ]R
0,
B(0,R) 0 0 0
d’où I = π. D’autre part le théorème de Fubini nous montre que

Z Z Z
2 2 2
I= e−y ( e−x dx)dy = ( e−x dx)2 .
IR IR IR
Ce calcul permet ainsi d’obtenir l’intégrale de la gaussienne pour laquelle on
2 √
ne dispose pas de primitive explicite : IR e−x dx = π.
R
Teminons par un résultat simple mais important : le théorème de majo-

ration.
Théorème 6.4 Soient f et g définies sur un ensemble E tel que f et |f |

sont localement intégrables sur E et |f (x)| ≤ g(x) pour tout x ∈ E. Si g est
intégrable sur E, alors f l’est aussi.
7 Intégrales curvilignes
7.1 Arcs orientés et courbes fermées orientées
De façon générale, une courbe du plan est définie par une fonction continue
γ = (γx , γy ) : IR → IR2 .
35
Les coordonnées γx et γy sont donc des fonctions de IR dans IR. Afin de
définir la notion d’intégrale curviligne, nous allons tout d’abord nous re-
streindre à certains types de courbes.
Définition 7.1 Un arc orienté Γ est défini par la donnée d’une fonction γ
continue et injective d’un intervalle [a, b] ⊂ IR dans IR2 . On note
Γ = {γ(t) ; t ∈ [a, b]},
On dit que l’arc est C 1 si γ est C 1 sur ]a, b[.
On dit que γ(t) est un paramètrage de Γ. Dans le cas d’un arc C 1 on

définit la vitesse du paramétrage γ au point γ(t) de la courbe par le vecteur
γ 0 (t) = (γx0 (t), γy0 (t)).
Lorsqu’il est non nul, le vecteur γ 0 (t) est un vecteur tangent à l’arc, de
longueur kγ 0 (t)k. Il faut cependant remarquer qu’aux points où γ 0 s’annule,
la courbe Γ n’admet pas nécessairement une tangente bien définie en tout
point même si γ est C 1 : considérer par exemple les courbes γ(t) = (t3 , t2 )
ou γ(t) = (t3 , |t|3 ) au point (0, 0).
L’injectivité entraine que la courbe Γ ne se croise pas elle-même. No-
tons que plusieurs paramètrages différents peuvent donner la même courbe.
L’orientation de l’arc Γ précise le sens dans lequel on le parcours, c’est à
dire de x0 = γ(a) vers x1 = γ(b): on voit ainsi que si on définit µ(t) =
γ(b + t(a − b)) pour t ∈ [0, 1] on obtient un paramètrage du même arc
mais orienté dans le sens inverse. Le résultat facile suivant montre comment
obtenir différents paramètrages du même arc avec la même orientation.
Théorème 7.1 Deux paramétrages γ sur [a, b] et µ sur [c, d] définissent le

même arc Γ avec la même orientation si et seulement si il existe ϕ une fonc-
tion continue et strictement croissante de [c, d] dans [a, b], telle que ϕ(c) = a
et ϕ(d) = b et µ = γ ◦ ϕ.
Définition 7.2 La longueur de l’arc orienté Γ est définie par l’intégrale de

la norme du vecteur vitesse
Z b
L= kγ 0 (t)kdt.
a
36
Il est facile de voir que cette définition est indépendante du paramétrage:
si µ = γ ◦ ϕ, on a
Z b Z d Z d
kγ 0 (t)kdt = kγ 0 (ϕ(s))kϕ0 (s)ds = kµ0 (s)kds.
a c c
On vérifie aussi que la longueur est invariante par changement d’orientation

du paramétrage.
Définition 7.3 Le paramétrage γ sur [a, b] est dit “normal” pour l’arc Γ si
et seulement si kγ 0 (t)k = 1 pour tout t ∈]a, b[. On a alors L = b − a.
Définition 7.4 Une courbe fermée orientée Γ est définie par la donnée
d’une fonction γ continue d’un intervalle [a, b] ⊂ IR dans IR2 , telle que
γ est injective sur [a, b[ et γ(a) = γ(b). On note
Γ = {γ(t) ; t ∈ [a, b]},
La fonction γ peut être étendue par périodicité à une fonction continue

définie sur tout IR de période b − a. On dit que la courbe est C 1 si γ ainsi
prolongée est C 1 sur IR.
A nouveau, plusieurs paramétrages différents peuvent aboutir à une

même courbe fermée. On note d’autre part que puisque la courbe est
fermée, on peut la parcourir à partir de n’importe quel point: par exemple
les fonctions γ(t) = (sin(t), cos(t)) et µ(t) = (sin(t + π/3), cos(t + π/3)) pour
t ∈ [0, 2π] sont deux paramétrages du cercle unité parcouru dans le même
sens. Remarquons aussi que selon notre définition ν(t) = (sin(2t), cos(2t))
ne donne par un paramétrage du cercle sur [0, 2π] car l’injectivité n’est pas
vérifiée, mais donne un paramétrage sur [0, π]. Par analogie avec le cas
des arcs orienté, on introduit la définition suivante pour l’orientation d’une
courbe fermée.
Définition 7.5 Un paramétrage µ de Γ défini sur [c, d] à la même orienta-

tion que γ si et seulement si il existe une fonction ϕ continue et strictement
croissante sur [c, d] telle que ϕ(d) − ϕ(c) = b − a et telle que µ = γ ◦ ϕ.
Deux sens de parcours sont uniquement possibles: celui “des aiguilles

d’une montre” dit sens indirect (ou négatif), et le sens inverse dit sens direct
(ou positif). Une façon de déterminer l’orientation du paramétrage est la
37
suivante: pour tout y ∈ IR tel qu’il existe des points d’ordonnée y contenus
dans la courbe Γ on considère
α(y) = min{x ; (x, y) ∈ Γ} et β(y) = max{x ; (x, y) ∈ Γ},
c’est à dire l’abscisse minimale et maximale des points de la courbe ayant

l’ordonnée y. On regarde alors la deuxième coordonée γy0 (t) du vecteur
vitesse aux points (α(y), y) et (β(y), y). On peut montre que deux cas sont
alors possibles:
1. Pour tout y, on a γy0 (t) ≤ 0 aux points (α(y), y) et γy0 (t) ≥ 0 aux points
(β(y), y) : la courbe est orientée dans le sens direct.
2. Pour tout y, on a γy0 (t) ≥ 0 aux points (α(y), y) et γy0 (t) ≤ 0 aux points
(β(y), y) : la courbe est orientée dans le sens indirect.
7.2 Intégrale curviligne et formule de Green-Riemann

Définition 7.6 Soit Γ un arc orienté ou une courbe fermée orientée de
paramétrage γ sur un intervalle [a, b], et soit g une fonction continue de IR2
dans IR. On pose
Z Z b Z Z b
g(x, y)dx = g(γ(t))γx0 (t)dt et g(x, y)dy = g(γ(t))γy0 (t)dt.
Γ a Γ a
En utilisant
R
la formule deRchangement de variable, on vérifie aisément que
les quantités Γ g(x, y)dx et Γ g(x, y)dy sont invariantes par changement de
paramétrage préservant l’orientation. En revanche, elle changent de signe
si on change d’orientation. Intuitivement dx = γx0 (t)dt et dy = γy0 (t)dt
correspondent à la mesure des déplacement infinitésimaux dans les directions
x et y lorsque l’on se déplace sur la courbe en faisant varier le paramètre t
de dt.
Il est aussi possible de définir une intégrale où le déplacement infinitésimal
est mesuré le long de la courbe.
Définition 7.7 Soit Γ un arc orienté ou une courbe fermée orientée de

paramétrage γ sur un intervalle [a, b], et soit soit g une fonction continue
de IR2 dans IR. On appelle intégrale de g sur Γ la quantité
Z Z b
g(x, y)ds = g(γ(t))kγ 0 (t)kdt.
Γ a
38
On peut vérifier que cette dernière quantité ne dépend pas de l’orientation
du paramétrage. On retrouve propriétés classiques des intégrales:
R R R
1. Linéarité: Γ αf + βg = α Γf +β Γg pour les intégrales par rapport
à dx, dy ou ds.
R R
2. Monotonie: f ≥ 0 ⇒ Γ fR = Γ 0 seulement pour l’intégrale par rap-
port à ds. En particulier Γ 1ds donne la longueur de Γ.
R R
3. RAdditivité: si on coupe Γ en deux parties (Γ1 , Γ2 ), on a Γ f = Γ1 f +
Γ2 f . Dans le cas des intégrales par rapport à dx ou dy, il faut supposer
que Γ1 et Γ2 ont la même orientation que Γ.
Définition 7.8 Soit Γ un arc orienté ou une courbe fermée orientée, de

paramétrage γ sur un intervalle [a, b], et soit g = (g1 , g2 ) un champ de
vecteur continu de IR2 dans IR2 . On appelle “circulation de g” le long de Γ
la quantité
Z Z b
g1 (x, y)dx + g2 (x, y)dy = g(γ(t)) · γ 0 (t)dt.
Γ a
Cette quantité est aussi appelée circulation de la 1-forme différentielle

g1 (x, y)dx + g2 (x, y)dy. On vérifie que c’est aussi l’intégrale sur Γ de la
fonction g(x, y) · τ (x, y) où τ est le vecteur unitaire tangent à la courbe
orienté dans la direction choisie.
Nous allons à présent établir un lien entre les intégrales curvilignes et les
intégrales multiples.
Théorème 7.2 Soit Ω un ensemble fermé dont la frontière est une courbe
Γ fermée et orientée dans le sens direct. Soit g une fonction de classe C 1
sur un ouvert contenant Ω. On a alors
∂g
Z Z
(x, y)dx dy = g(x, y)dy
Ω ∂x Γ
et
∂g
Z Z
(x, y)dx dy = − g(x, y)dx
Ω ∂y Γ
Il est important de noter ici la disymmétrie des rôles des variables x et

y par la présence du signe − dans la deuxième formule. Dans le cas d’une
courbe orientée dans le sens indirect le signe − est alors dans la première
formule. On peut facilement démontrer ces formules dans le cas particulier
39
où les sections horizontales et verticales de Ω sont des intervalles fermés.
Pour la première formule par exemple, on considère les intervalles
[α(y), β(y)] = {x ; (x, y) ∈ Ω}, (7.1)
où y varie dans l’intervalle [y0 , y1 ] en dehors duquel les sections horizontales
de Γ sont vides. On applique alors le théorème de Fubini qui donne
R β(y)
∂g
= yy01 ( α(y) ∂x
∂g
R R
Ω ∂x (x, y)dx dy (x, y)dx)dy
= y0 g(β(y), y))dy − yy01 g(α(y), y))dy
R y1 R
R R
= Γβ g(x, y)dy − Γα g(x, y)dy,
où Γα et Γβ sont les deux parties de Γ définies respectivement par les
parmétrage γα (y) = (α(y), y) et γβ (y) = (β(y), y) pour y ∈ [y0 , 1]. En
remarquant que Γβ est orientée dans le même sens que Γ et Γα dans le sens
inverse, et en appliquant la propriété d’additivité, on a
Z Z Z
g(x, y)dy − g(x, y)dy = g(x, y)dy,
Γβ Γα Γ
ce qui prouve la premiere formule. La preuve de la deuxième est similaire.

La formule de Green-Riemann synthétise les deux formules de la façon
suivante: si g = (g1 , g2 ) est une fonction de classe C 1 sur un ouvert contenant
Ω dans IR2 , on a
∂g2 ∂g1
Z Z
g1 (x, y)dx + g2 (x, y)dy = ( (x, y) − (x, y))dxdy.
Γ Ω ∂x ∂y
Une conséquence immédiate de cette formule est la suivante : la circulation
sur Γ d’une 1-forme différentielle fermée sur Ω est nulle.
Si on applique cette formule à la fonction g̃ = (−g2 , g1 ), on obtient
∂g1 ∂g2
Z Z
g1 (x, y)dx − g2 (x, y)dy = (
(x, y) + (x, y))dxdy.
Γ Ω ∂x ∂y
Le membre de gauche correspond à l’intégrale sur Γ du produit scalaire
g̃(x, y) · τ (x, y) qui est égal à g(x, y) · n(x, y) où n(x, y) désigne le vecteur
unitaire normal à Γ au point (x, y) et pointant vers l’extérieur de Γ. Cette
quantité est appelée flux du champ de vecteur g sur Γ. La quantité intégrée
dans le membre de droite est appellée divergence de g, plus généralement
∂gi
pour une fonction g = (g1 , · · · , gn ) de IRn dans IRn par div(g) = ni=1 ∂x
P
i
.
On a ainsi établi la formule de la divergence qui exprime que le flux de g sur
Γ est égal à l’intégrale de sa divergence sur Ω:
Z Z
g(x, y) · n(x, y)ds = div(g)(x, y)dx dy.
Γ Ω
40

Notes LM216

Transféré par

Droits d'auteur :

Formats disponibles

Notes LM216

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Notes LM216

Transféré par

Droits d'auteur :

Formats disponibles

Fonctions de plusieurs variables

Notes du Cours LM 216

Dans ce cours, on s’intéresse aux fonctions de plusieurs variables

Celles-ci interviennent naturellement pour décrire la dépendance de grandeurs

Avertissement : ces notes sont régulièrement mises à jour et corrigées.

1 Notions de topologie dans IRn

Le produit scalaire de x avec y = (y1 , · · · , yn ) est défini par

Inégalité de Cauchy-Schwartz : pour tout x et y dans IRn

avec égalité si et seulement si x et y sont proportionels.

Inégalité triangulaire : pour tout x et y dans IRn

On rappelle aussi que si θ désigne l’angle entre deux vecteurs x et y non-nuls,

c’est à dire la longeur du vecteur reliant le point x au point y.

2. N (λx) = |λ|N (x) pour tout x ∈ IRn et λ ∈ IR.

On vérifie facilement que la norme euclidienne vérifie ces trois propriétés.

kxk∞ := sup |xi |,

la norme `p définie pour 1 ≤ p < +∞ par

kxkp := (|x1 |p + · · · + |xn |p )1/p .

On voit que le cas p = 2 correspond à la norme euclidienne.

On appelle boule unité la boule B(0, 1) de centre 0 et de rayon 1. En

Définition 1.2 Deux normes N1 et N2 sont dites “équivalentes” si et seule-

cN1 (x) ≤ N2 (x) ≤ CN1 (x).

kxk2 ≤ kxk1 ≤ n1/2 kxk2 .

Définition 1.3 Un ensemble U ⊂ IRn est dit “ouvert” si et seulement si

Par convention l’ensemble vide ∅ est ouvert.

Pour un ensemble E quelconque on dit que x est un point intérieur à E

Théorème 1.2 F est fermé si et seulement toute suite de point de F qui

Pour un ensemble E quelconque on dit que x est une valeur d’adhérence

Quelques propriétés élémentaires:

1. Toute union finie ou infinie d’ouverts est un ouvert.

2. Toute union finie de fermés est un fermé

3. Toute intersection finie ou infinie de fermés est un fermé.

4. Toute intersection finie de d’ouverts est un ouvert

5. Les seuls ensembles à la fois ouverts et fermés sont ∅ et IRN .

6. Un ensemble fini de points de IRn est fermé.

Définition 1.5 On appelle frontiere d’un ensemble E l’ensemble des points

Définition 1.6 Un ensemble E ⊂ IRn est dit borné si et seulement si il

2. Propriété de Bolzano-Weierstrass : toute suite (xk )k≥0 de points de E

Définition 1.7 Un ensemble E est dit “connexe” si et seulement si il n’existe

Autrement dit, E est connexe si on ne peux pas le séparer en deux partie

Théorème 1.4 Un ensemble U est ouvert et connexe si et seulement si

Définition 1.8 Un ensemble E est dit “convexe” si et seulement si pour

On peut vérifier qu’en dimension n = 1 les connexes et les convexes de

f (E) = {f (x) ; x ∈ E}.

f −1 (F ) = {x ∈ IRn ; f(x) ∈ F}.

Gf = {(x, f (x)) ; x ∈ IR},

et que l’on représente par la courbe du plan d’équation y = f (x). Dans le

Définition 2.1 Une fonction f : IRn → IR est continue au point x =

Une fonction f est continue sur un ensemble E si et seulement si elle est

Définition 2.2 Une fonction f : IRn → IRm est lipschtizienne si et seule-

On dit parfois que f est C-lipschitzienne. On voit immédiatement que

2. Si f : IRn → IRm est continue au point x et g : IRm → IRp est continue

3. La fonction f : IRn → IR qui à x = (x1 , · · · , xn ) associe la i-ème

Une manière naturelle d’étudier les variations d’une fonction f de n

fi,x (t) = f (x1 , · · · , xi−1 , t, xi+1 , · · · , xn ).

2.3 Fonctions continues et ensembles compacts

Intuitivement cette derniere quantité mesure les variations de f lorsqu’on

Définition 3.3 Soit f une fonction définie sur un ouvert U ⊂ IRn et à

f (x + h) = f (x) + dfx (h) + khkε(h),

avec limh→0 ε(h) = 0. L’application dfx est appelée différentielle de f au