École Normale Supérieure – Année 2012-2013

T.D. de Topologie, Analyse et Calcul différentiel

Corrigé de la feuille d’exercices no 12

1. Lemme de Grönwall et variantes

 R R
t t
1. On introduit donc v(t) := exp − a β(s) ds a β(s)u(s) ds, que l’on dérive en temps :

 Z t   Z t Z t
v 0 (t) = β(t)u(t) exp − β(s) ds − β(t) exp − β(s) ds β(s)u(s) ds
a a a
 Z t 
≤ β(t)α(t) exp − β(s) ds .
a

On intègre entre a et t pour obtenir

 Z t Z t Z t  Z s 
exp − β(s) ds β(s)u(s) ds ≤ β(s)α(s) exp − β(τ ) dτ ds.
a a a a

D’où en réutilisant l’hypothèse

Z t
u(t) ≤ α(t) + β(s)u(s) ds
a
Z t Z t  Z s 
≤ α(t) + exp β(s) ds β(s)α(s) exp − β(τ ) dτ ds,
a a a

ce qui donne l’inégalité recherchée.

2. (Une généralisation due à T. Cazenave et A. Haraux). On introduit comme indiqué t. Il est clair que
pour ε > 0, on a t > a. Supposons qu’il ne soit pas égal à b : on voit que pour t ≤ t on a
Z t
w(t) − vw0 +ε (t) ≤ −ε + [f (w(s)) − f (vw0 +ε )(s)] ds,
a

et on trouve une contradiction par la monotonie de f . Ensuite, on fait tendre ε vers 0, on utilise la dépendence
continue des solutions par rapport à la donnée initiale dans Cauchy-Lipschitz et on obtient le résultat.

2. Quelques exemples élémentaires d’utilisation de Cauchy-Lipschitz

1. Soit f : R → R de classe C 1 telle que |f (t) − cos(t)| < 1 (pour tout t ∈ R). On considére le système
différentiel  0
x (t) = f (x(t))
(S)
x(t0 ) = x0

1. Comme f est de classe C 1 , elle est localement lipchitzienne et l’existence d’une unique solution maxi-
male résulte du Théorème de Cauchy-Lipschitz. Il reste à voir que cette solution maximale est définie
sur tout R. Or, pour tout t ∈ R on a |f (t)| ≤ 2, et donc par le critère d’explosion en temps fini, les
solutions maximales sont définies sur R.
2. On remarque que −1 + cos(x) ≤ f (x) ≤ 1 + cos(x), d’où f (2kπ) ≥ 0 et f ((2k + 1)π) ≤ 0 et l’existence
de tk suit du Théorème des valeurs intermédiaires.
 3. Soit k tel que x0 ∈ [tk , tk+1 [. Si x0 = tk , la solution est stationnaire. Si, pour tout t, x(t) ∈]tk , tk+1 [,
alors, x(t) est bornée. Supposons qu’il existe t tel que x(t) ∈]t / k , tk +1[. Comme x0 ∈ [tk , tk+1 [, par
le théorème des valeurs intermédiaires, il existe t tel que x(t) = tk ou x(t) = tk+1 . Mais alors, les
solutions constantes x = ti sont les uniques solutions de (S) qui prennent la valeur ti en un point (par
Cauchy-Lipschitz). Par conséquent, une telle fonction x est constante donc bornée!

2. Soient f, g : [0, 1] × R → R deux fonctions continues vérifiant f (t, x) ≤ g(t, x) pour tout (t, x) ∈ [0, 1] × R.
On fixe t0 ∈ [0, 1[ et a ∈ R et on considère les systèmes
 0  0
x (t) = f (t, x(t)) y (t) = g(t, y(t))
(Sf ) = et (Sg ) =
x(t0 ) = a y(t0 ) = a.

1. Comme x et y sont des solutions (C 1 ) de (Sf ) et (Sg ), il existe un intervalle (semi-)ouvert [t0 , t0 + ε[
contenant t sur lequel elles sont toutes deux définies. Donc pour tout t ∈ [t0 , t0 + ε[, on a
Z t Z t
x(t) − y(t) = x0 (s) − y 0 (s) ds = f (s, x(s)) − g(s, y(s)) ds (1)
t0 t0
Z t



= f (s, x(s)) − g(s, x(s)) + g(s, x(s)) − g(s, y(s)) ds (2)
t0

Comme f − g est continue sur le compact [0, 1] × K, où K = x([0, ε/2]) ∪ y([0, ε/2]) est bien compact
(par union d’images continue d’un compact dans un séparé), f − g y a un maximum −α < 0. De
plus, g est uniformément continue sur [0, 1] × K, donc il existe η > 0 tel que, pour tout |x − y| ≤ η et
s ∈ [0, 1], on ait |g(s, x(s)) − g(s, y(s))| ≤ α/2. Comme x(t0 ) = y(t0 ), il existe un intervalle [t0 , t0 + δ]
(on peut choisir δ < ε/2) sur lequel |x(t) − a| ≤ η/2 et |y(t) − a| ≤ η/2, en particulier |x(t) − y(t)| ≤ η.
On déduit alors de (2) que, pour tout t ∈ [t0 , t0 + δ[, on a
x(t) − y(t) ≤ −(t − t0 )α + (t − t0 )α/2 < 0.

2. On considère J = {t ∈ [t0 , 1], pour tout s ∈ [t0 , t], x(t) ≤ y(t)}. On note t+ = supt∈J (t). Il est clair
que J est un intervalle d’intérieur non-vide et contenant t0 (d’après la question précédente). Il est donc
de la forme [t0 , t+ ] ou [t0 , t+ [. De plus, par continuité, on a clairement par continuité que x(t+ ) ≤ y(t+ ).
Il reste donc à montrer que t+ = 1. Sinon, si x(t+ ) < y(t+ ), par continuité, ceci est encore vrai dans
un voisinage ouvert de t+ ce qui contredirait sa maximalité. Par conséquent, x(t+ ) = y(t+ ) et alors la
question précédente permet uen nouvelle fois de contredire sa maximalité.

3. Soit f : Rn → R de classe C 2 telle que lim f (x) = +∞.

kxk→∞

1. Déjà, f C 2 , donc X est C 1 et par Cauchy-Lipschitz, on a l’existence d’une solution maximale sur un
intervalle de la forme ]a, b[. On veut montrer que b = +∞. Montrons que t 7→ f (x(t)) est décroissante.
On a  
n
∂ X ∂f 0 X  ∂f 
(f (x(t))) = xi (t) = − x0i (t) ≤ 0

∂t ∂xi  ∂xi


i=1

car x est solution de (SX ). Il suit que f (x(t)) ≤ f (x0 ) pour tout t ≥ t0 . Mais alors, par le critère
d’explosion, si b < +∞, alors |x(t)| → +∞ quand t tend vers b, et donc f (x(t)) est non borné quand t
tend vers b ce qui contredit f (x(t)) ≤ f (x0 ). Il suit que b = +∞.
2. Pour f (x) = x4 /4, l’équation différentielle devient x0 (t) = −x3 (t). On choisit t0 = 0 et x(0) 6= 0.
Notons que, par unicité de la solution, ceci implique que x(t) n’est jamasi nul sur l’intervalle [a, +∞[
car la fonction nulle est solution (la seule constante) du système. En particulier, on peut donc diviser
l’équation différentielle par x(t) sur ]b, +∞[ ce qui donne:
x0 (t) 1 1
=1 ⇐⇒ − = t − 0.
x3 (t) 2x2 (t) 2x(0)2

1
x0
Ceci force que t > − 2x(0) 2 et x(t) = p et donc b = − 2x1 2 .
1 + 2x20 t 0
3. Équations linéaires à coefficients constants
1. On connaı̂t la solution explicite de l’équation ẋ = Ax, c’est x(t) = etA x0 . Calculons etA . Soit λ1 , ..., λp les
valeurs propres de A, de multiplicités
Lrespectives m1 , ..., mp . On décompose Cn en somme directe des sous-
p
espaces caractéristiques de A: C = j=1 Ker(A−λj Id)mj . Soit Aj la restriction de A à Ker(A−λj Id)mj ,
n
m
on a Aj = λj Id +Nj , avec Nj nilpotent (Nj j = 0). Les sous-espaces caractéristiques sont stables par etA , et
k
on a sur Ker(A − λj Id)mj , etA = etλj etNj = etλj (Id +tNj + ... + tk! Njk ) pour un certain entier k (dépendant
de j). Ce calcul va permettre de répondre à la question.
(i) Une condition nécessaire et suffisante pour que toutes les solutions de ẋ = Ax tendent vers 0 en +∞
est que les valeurs propres de A soient de partie réelle strictement négative. En effet en prenant pour x0 un
vecteur propre, on voit que la condition est nécessaire. Montrons qu’elle est suffisante. Soit x0Pune donnée
intiale quelconque, on décompose x0 sur la somme directe: x0 = x1 + ... + xp , alors etA x0 = j etAj xj , il
tk
suffit donc de prouver que les etAj tendent vers 0. Or on a ketAj k ≤ et Re(λj ) (1 + tkNj k + ... + k
k! kNj k) →0
quand t → +∞.
(ii) La condition nécessaire et suffisante pour que les solutions soient bornées sur R+ est que Re λj ≤ 0
pour tout j, et Nj = 0 si Re λj = 0 (i.e. sur les sous-espaces caractéristiques associés aux valeurs propres
imaginaires pures, la matrice est diagonalisable). La condition est suffisante, en effet si Re(λj ) < 0, on
vient de voir que etAj tend vers 0, donc est borné, et si Re λj = 0, on a etAj = etλj Id qui est borné
aussi, donc etA est borné sur R+ . La condition est nécessaire: si on avait Re(λj ) > 0, alors tout vecteur
propre x0 associé à λj donnerait une solution non bornée, et si Re(λj ) = 0 mais que Nj 6= 0, on prend
un x0 dans Ker(A − λj Id)mj tel que Njk x0 6= 0 (ici k est le plus grand entier tel que Njk 6= 0). Alors
tk tk tk
etA x0 = etλj (x0 + tNj x0 + ... + k
k! Nj x0 ), d’où ketA x0 k = kx0 + ... + k
k! Nj x0 k ∼ k
k! kNj x0 k → +∞.
(iii) Enfin, la condition pour que les solutions soient bornées sur R est que A et −A vérifient la condition
précédente, ce qui donne: le spectre de A est imaginaire pur et A est diagonalisable.

2. La condition cherchée est ∀j, Re(λj ) > 0. C’est nécessaire, en effet si on prend b = 0, on a une solution
bornée évidente (la solution nulle). On veut que ce soit la seule, i.e. pour tout x0 6= 0, etA x0 doit être non
borné. En particulier pour x0 vecteur propre, on voit que les valeurs propres doivent avoir une partie réelle
> 0.
Montrons que cette condition est suffisante. La solution de l’équation ẋ = Ax + b, avec x(0) = x0 ,
s’obtient par la méthode de la variation des constantes: on pose x(t) = etA y(t), et on forme une équation
Rt Rt
sur y. Il vient ẏ = e−tA b, d’où y(t) = 0 e−sA b(s) ds, puis x(t) = etA x0 + 0 e(t−s)A b(s) ds. Il en résulte
facilement qu’il y a au plus une solution bornée: si x0 6= x̃0 , alors x(t) − x̃(t) = etA (x0 − x̃0 ) n’est pas
borné, autrement dit la différence entre deux solutions distinctes n’est Rjamais bornée. Montrons qu’il existe
t
une solution bornée. Réécrivons x(t) sous la forme x(t) = etA (x0 + 0 e−sA b(s) ds). On devine alors en
R∞
faisant tendre t vers +∞ que la solution bornée va être obtenue en choisissant x0 = − 0 e−sA b(s) ds.
Cette intégrale est bien convergente, en effet il existe des constante C et α > 0 telles que ke−sA k ≤ Ce−sα .
Pour le prouver, il suffit d’obtenir une majoration du même type sur chaque sous-espace caractéristique:
k
ke−sAj k ≤ e−s Re λj (1 + ... + sk! kNj kk ) ≤ Cj e−sα en prenant α < min Re(λj ). La solution x(t) associée à ce
R ∞ (t−s)A R∞
x0 est donnée par x(t) = t e b(s) ds. Elle est effectivement bornée: kx(t)k ≤ t ke(t−s)A kkbk∞ ds ≤
R ∞ −uα
Ckbk∞ 0 e du = Cα kbk∞ . Ceci achève la preuve.

4. Une étude de systèmes différentiels

2
1. On note X = (x, y). L’équation est de la forme Ẋ = V (X), avec V (X) = (λ−e|X| )X, qui est de classe C 1 ,
donc localement lipschitzienne. Le théorème de Cauchy-Lipschitz fournit alors l’existence locale et l’unicité
des solutions pour toute donnée initiale. Cherchons les solutions constantes : ce sont les X(t) = (x0 , y0 ) tels
que V (x0 , y0 ) = 0, ce qui donne (x0 , y0 ) = (0, 0) ou x20 + y02 = ln λ.
2. Soit u(t) = x(t) + iy(t) = r(t)eiθ(t) . Pour une solution u non indentiquement nulle, on a u(t) 6= 0 pour
tout t (d’après Cauchy-Lipschitz), alors le théorème de relèvement garantit l’existence d’une détermination
2
de classe C 1 de la fonction θ. On a u̇ = ẋ + iẏ = (λ − er )u d’une part, et d’autre part u̇ = ṙeiθ + irθ̇eiθ ,
2
d’où r + irθ̇ = (λ − er )r. On en déduit que θ̇ = 0, ainsi θ est constant, les trajectoires sont donc rectilignes.
2
On en déduit aussi que ṙ = (λ − er )r. Ainsi si r2 < ln λ, alors ṙ > 0, r est croissant, et si r2 > ln λ alors
ṙ < 0, r est décroissant. Le signe de ṙ ne peut varier au cours du temps, sinon la trajectoire rencontrerait
une solution constante de module ln λ.
√
3. Lorsque r02 < ln λ, la solution considérée dans R2 est confinée dans le cercle de rayon √ ln λ, donc n’explose
pas : elle est globale. Une autre façon de le dire est que le champ est complet sur B(0, ln λ).
Lorsque r02 > ln λ, r est décroissant donc pour t ≥ 0 on a r(t) ≤ r0 , la solution est donc bien définie sur
2
R+ . Pour t ≤ 0, on a r(t) ≥ r0 . Soit C > 0 tel que ∀x > r0 , ex − λ ≥ Cx. On a pour tout t ≤ 0, ṙ ≤ −Cr2 ,
ou encore − rṙ2 ≥ C. On intègre cette relation entre t < 0 et 0, il vient r(t) ≥ r−11+Ct (t ≤ 0). Il y a donc
0
explosion à une date −T ∗ ∈ ]−(Cr0 )−1 , 0[.
2
4. Les limites possibles pour r en ±∞ sont nécessairement des zéros de la fonction g(r) = r(λ − er ). Notons
l = lim r(t) en +∞. Si par l’absurde g(l) 6= 0, alors on a par exemple g(l) > 0, et g est > 0 au voisinage de
l: il existe  > 0, δ > 0 tels que g(r) ≥ δ si |r − l| ≤ . Pour t assez grand (t > t0 ), on a |r(t) − l| ≤ , d’où
ṙ = g(r) ≥ δ. En intégrant, il vient r(t) ≥ r(t0 ) + δ(t − t0 ), ce qui contredit manifestement le fait que r(t)
converge vers l. Ainsi on a bien g(l) = 0.
Lorsque 0 < r02 < λ, la fonction r(t) est croissante. Par √
monotonie, elle admet des limites quand t → ±∞.
Les seules limites possibles sont alors: lim√ t→+∞ r(t) = ln λ et limt→−∞ r(t) = 0. Dans le second cas,
r02 > ln λ, on a de même limt→+∞ r(t) = ln λ.
√
5. La trajectoire est soit un segment ouvert joignant l’origine à un point du cercle de rayon ln λ, soit une
demi-droite radiale issue d’un point de ce même cercle.

5. Une étude de systèmes différentiels II

1. On a:

• L’équation est donnée par une fonction F (x, y) de classe C 1 , on a donc existence locale et unicité
des solutions. La fonction E est une intégrale première: pour tout solution x(t), y(t), la fonction
E(x(t), y(t)) est constante. Pour le prouver, il suffit de dériver et de remplacer à l’aide de l’équation.
On peut aussi constater que l’équation se réécrit ẋ = ∂E ∂E
∂y , ẏ = − ∂x . On dit de ce système qu’il est
hamiltonien, de fonction de Hamilton (ou hamiltonien) E.
√ √
• On a E(x, y) = √ x2 + (x2 + y)2 . On en déduit que x2 ≤ E et que |y| ≤ E + x2 ≤ E + E, d’où
|x| + |y| ≤ E + 2 E, et il en résulte que E(x, y) → +∞ quand k(x, y)k → ∞.
Puisque E est une intégrale première, toute solution du système est contenue dans une ligne de niveau
de E. Ces lignes de niveau sont bornées, donc toute solution est bornée. Il en découle que les solutions
sont globales.
• Déterminons les√lignes de niveau de E. La courbe E(x, y)√= C est la réunion
√ du graphe de la fonction
f1 (x) = −x2 + C − x2 et du graphe de f2 (x) = −x2 − C − x2 (|x| ≤ C). C’est donc une √ courbe
fermée simple
√ (on pourra représenter graphiquement les lignes de niveau de E). Soit A = ( C, −C)
et B = (− C, −C) les points de la courbe d’abscisse maximale et minimale respectivement.
Supposons par exemple que (x0 , y0 ) est un point du graphe de f1 . Tant p que la solution parcourt le
graphe de f1 , on a y(t) = f1 (x(t)), d’où ẋ(t) = 2x(t)2 + 2f1 (x(t)) = 2 C − x(t)2 . On intègre cette
Rt R x(t) du
équation à variables séparées: t = 0 √ẋ(s) ds 2 = x0 2√C−u 2
, à l’aide du changement de variable
2 C−x(t)
u = x(t). Clairement x(t)√est croissante, et la solution va atteindre le point A en un temps fini, à une
R C du
date t1 donnée par t1 = x0 2√C−u 2
< ∞ (sinon la précédente relation fournit une contradiction pour
 tout t > t1 ). Une fois arrivée
√ au point A, la solution parcourt le graphe de f2 en rétrogradant, et on a
alors ẋ = 2x2 +2f2 (x) = −2 C − x2 , tant que l’on n’a pas atteint le point B. Ceci s’intègre de la même
R √C du
R √C du
√
façon en t − t1 = x(t) 2√C−u2 , et on voit que l’on atteint B à la date t2 = t1 + − C 2√C−u . Ensuite,
2

on est à nouveau sur le graphe de f1 , et par le même calcul que précédemment, on arrive en A à la date
R √C du
R √C du
R1 dv
t3 = t2 + −√C 2√C−u2 . Ainsi la solution est périodique, de période T = −√C √C−u2 = −1 √1−v2 = π.

d
2. On réutilise la fonction E de la question précédente. On a cette fois dt E(x(t), y(t)) = −2ax(t)2 . Supposons
d
par exemple a > 0. Il vient dt E(x(t), y(t)) ≤ 0, donc E décroı̂t. Alors pour tout t positif, la solution
est contenue dans le domaine borné {(x, y) | E(x, y) ≤ E(x(0), y(0))}, donc la solution existe pour tout
d
t > 0. Pour les temps négatifs, on remarque que dt E(x, y) ≥ −2aE(x, y). Cette inégalité est du type
2at d
Gronwall : on pose φ(t) = E(x(t), y(t))e , il vient dt φ ≥ 0, d’où φ(t) ≤ φ(0) pour t < 0, c’est-à-dire
−2at
E(x, y) ≤ E(x(0), y(0))e . Il en résulte que la solution est bornée sur tout intervalle borné, donc elle est
globale (il ne peut y avoir d’explosion en temps fini).

3. ...

6. Equation de transport
On introduira X(t, x) la solution au système différentiel X 0 (t, x) = A(t, X(t, x)), avec X(0) = x. On calculera
alors ∂t f (t, X(t, x)).
On en déduira que f (t, x) = f0 (Y (0, x)), en notant Y (t, x) la solution à

d
Y (s, x) = A(t, Y (s, x)), Y (t, x) = x.
ds

Autrement dit, la donnée initiale f0 est transportée le long des lignes de champ de A.

7. Théorème de Hadamard
1. f −1 est continue à valeur dans un espace séparé, donc transforme en compact, ce qui prouve que f est
propre. Il est {evident que i) implique que df(x) est un isomorphisme pour tout x.
Par le théorème de l’image ouverte (ou d’inversion locale), f (Rn ) est un ouvert du connexe Rn . Il suffit
de montrer qu’il est fermé pour conclure à la surjectivité. Soit (f (xn ))n∈N uen suite de f (Rn ) qui converge
vers y ∈ Rn . Il suffit de montrer que y ∈ f (Rn ). La réunion K = {f (xn ), n ∈ N} ∪ {y} est un compact,
donc (xn )n∈N est une suite à valeur dans le compact f −1 (K). Elle admet une sous-suite (xφ(n) )n∈N qui
converge vers x ∈ f −1 (K). Par continuité, f (xφ(n) ) → f (x) et f (xφ(n) ) → y par définition de y. Il suit que
y = f (x) ∈ f (Rn ).

2. Soit x ∈ S (donc f (x) = f (z)), comme df(x) est inversible, il existe un voisinage de x sur lequel f est un
difféomorphisme, et donc x est l’unique solution de f (y) = f (z) dans un voisinage de x. Ceci assure que S
est discret. Comme f est propre, S est aussi compact, on en déduit qu’il est fini.

3. On note {z1 , . . . , zm } l’ensemble S (c’est à
dire les solutions de f (x) = f (z)). Soit g : Rn → Rn la fonction
définie par g(x) = (df(x) )−1 · (f (x) − f (z)) qui est de classe C 1 par hypothèse. En particulier le système
différentiel 
ẋ(t) = −g(x(t))
(Sz ) =
x(0) = x0
admet une solution maximale définie sur un intervalle [0, b[. Notons que les solutions constantes de (Sz ) sont
exactement données par les points zi de S.
 Montrons que b = +∞. Soit F (x) = hf (x) − f (z); f (x) − f (z)i. Notons que f propre implique que F est
propre. On dérive la fonction t 7→ F (x(t)) (où x(t) est solution de (Sz )):

∂

F (x(t)) = 2hf (x) − f (z); dfx(t) · ẋ(t)i
∂t
= −2hf (x) − f (z); f dfx(t) · dfx−1 ((x) − f (z)) i

= −2F (x(t)).

Il suit que F (x(t)) = F (x0 ) exp(−2t). On en déduit que x(t) ∈ F −1 [0, F (x0 )] qui est un compact car F est
propre. On en conclut que x(t) est bornée pour t ≥ 0, donc sa solution maximale est définie sur [0, +∞[.
Pour tout i = 1 . . . m, il existe une boule ouverte B(zi , δi ) de centre i, telle que la restriction de f à cette
boule soit un difféomorphisme sur un voisinage de f (zi ) = f (z) dans Rn . On peut restreindre le voisinage
d’arrivée à une boule ouverte B(f (zi ), ε). En faisant le changement de variable y(t) = f (x(t)) − f (z) (qui
ets bien aussi un difféomorphisme) au voisinage de zi (c’est à dire pour tout y(0) ∈ B(0, ε)), le système
(Sz ) devient ẏ(t) = −y(t) dans un voisinage de zi , dont la solution est y(t) = y(0) exp(−t). On voit que
cette solution reste dans B(0, ε) pour t ≥ 0, et donc est définie sur [0, +∞[ et de plus elle tend vers 0 quand
t → +∞. Par conséquent, pour x0 ans B(zi , δi ), on a x(t) = h(y(0) exp(−t)) qui tend vers zi quand t → +∞.
On vient de montrer que les solutions constantes x(t) = zi sont des points critiques stables.
On montre maintenant que pour tout x0 , la solution x(t) tend vers un des zi lorsque t → +∞. On note
A(zi ) = {x0 ∈ Rn , x(t) → zi } (notons que par unicité des solutions maximales, et comme celles ci sont
t→+∞
définies sur [0, +∞[, A(zi ) est bien défini et non-vide vu ci-dessus). De fait, on a vu que x([0, +∞[) est
inclus dans un compact C. On en déduit qu’il existe une suite x(φ(n))n∈N qui converge dans C. De plus on
a vu que F (x(t)) = F (x0 ) exp(−2t) −→ 0. Donc lim F (x(φ(n))) = 0 d’où x(φ(n)) converge vers un des
t→+∞
zi . Mais alors il existe un n tel que x(φ(n)) ⊂ B(zi , δi ) et d’après le paragraphe ci-dessus, et l’unicité des
solutions maximales associées au système différentiel pour la valeur initiale x(φ(n)) en φ(n), on en déduit
que x(t + φ(n)) −→ zi , d’où x(t) −→ zi . Par conséquent, on a bien que pour tout x0 , x(t) converge vers
t→+∞ t→+∞
un des zi lorsque t tend vers +∞.
Par ailleurs, les ensembles A(zi ) sont ouverts. En effet, si x0 ∈ A(zi ), il existe t0 tel que x(t0 ) ∈ B(zi , δi /2).
Mais par continuité des solutions du système diffŕentiel par rapport aux conditions initiales, il existe une
boule ouverte B(x0 , η) telle que pour tout y0 ∈ B(x0 , η), la solution maximale y(t) associée à la conditon
initiale y(0) = y0 vérifie |y(t0 ) − x(t0 )| < δi /2. Par conséquent, y(t0 ) ∈ B(zi , δi ), et donc, les arguments
précédents assurent que y(t) −→ zi ; c’est à dire y0 ∈ A(zi ). On a bien montré que A(zi ) est ouvert.
t→+∞
Sm
Finalement, on a montré que Rn = i=1 A(zi ) où les A(zi ) sont des ouverts disjoints (car S est dis-
cret). Par connexité de Rn , m = 1 et donc f est injective. Par le théorème d’inversion globale, f est un
difféomorphisme global !