# Théorème de Cauchy-Lipschitz

En mathématiques et plus précisément en analyse, le théorème de Cauchy-Lipschitz, appelé également théorème de Picard-Lindelöf ou théorème d'existence de Picard, concerne les solutions d'une équation différentielle. Sous certaines hypothèses de régularité de la fonction définissant l'équation, le théorème garantit l'existence d'une solution répondant à une condition initiale dite de Cauchy et l'unicité d'une solution maximale.
Certaines lois physiques, comme le principe fondamental de la dynamique, se traduisent par des équations différentielles vérifiant les hypothèses du théorème. Ce dernier assure alors le caractère déterministe du mécanisme décrit par la loi. Ce déterminisme ne se traduit pas toujours par une possibilité de prédiction, la théorie du chaos montre l'existence de possibles phénomènes fortuits.
Selon les auteurs, le théorème de Cauchy-Lipschitz s'exprime de manière plus ou moins forte. Sous une forme plus élaborée, ce théorème assure que la solution varie continûment si la condition initiale est modifiée, et il en est de même si la fonction définissant l'équation dépend continûment d'un paramètre. Si l'équation est définie par une fonction de classe Cp, la solution est de classe Cp+1. Ce théorème peut encore être généralisé au cas où l'équation différentielle n'est plus à valeurs dans un espace vectoriel, mais dans une variété différentielle.
Une première version est démontrée par Augustin-Louis Cauchy durant la première moitié du XIXe siècle, à l'aide d'une technique d'approximation découverte par Leonhard Euler au siècle précédent. Rudolf Lipschitz généralise l'énoncé en élargissant un peu la classe des équations qui s'y rapportent. Le théorème n'en reste pas moins uniquement un résultat d'existence locale. C'est à la fin de ce siècle que les techniques de démonstration, ainsi que l'énoncé du théorème, sont profondément modifiés. À la suite des travaux de Lazarus Fuchs, les mathématiciens Émile Picard, Paul Painlevé et Henri Poincaré développent une version moderne de l'analyse des équations différentielles. Cette vision permet d'apporter des éléments de réponse sur les solutions maximales, l'unicité et la régularité de la solution. Une version relativement moderne est publiée en 1894 par Ernst Lindelöf. Le théorème se démontre maintenant généralement à l'aide d'un théorème du point fixe et d'une approche topologique, classique en analyse fonctionnelle.


## Préambule
### Énoncés
Dans tout le reste de l'article, 
  
    
      
        E
      
    
    {\displaystyle E}
  
 désigne un espace de Banach, 
  
    
      
        Ω
      
    
    {\displaystyle \Omega }
  
 un ouvert de 
  
    
      
        
          R
        
        ×
        E
      
    
    {\displaystyle \mathbb {R} \times E}
  
 et 
  
    
      
        f
      
    
    {\displaystyle f}
  
 une fonction continue de 
  
    
      
        Ω
      
    
    {\displaystyle \Omega }
  
 dans 
  
    
      
        E
      
    
    {\displaystyle E}
  
. L'objectif est d'étudier l'équation différentielle du premier ordre suivante :

  
    
      
        (
        1
        )
        
        
          
            
              
                d
              
              x
            
            
              
                d
              
              t
            
          
        
        (
        t
        )
        =
        f
        (
        t
        ,
        x
        (
        t
        )
        )
      
    
    {\displaystyle (1)\quad {\frac {\mathrm {d} x}{\mathrm {d} t}}(t)=f(t,x(t))}
  

avec la condition de Cauchy C : x(t0) = x0, où le couple (t0, x0) est un élément de Ω. (On ne se préoccupe que des équations différentielles d'ordre 1 parce que celles d'ordre n s'y ramènent.) On peut remarquer que si f est une fonction de classe Cp, la dérivée de x l'est aussi, et les solutions sont donc de classe Cp+1.
Il existe plusieurs manières d'exprimer le théorème de Cauchy-Lipschitz, dont la forme suivante :

Si la fonction f est localement lipschitzienne par rapport à la deuxième variable, il existe une et une seule solution maximale à l'équation (1) respectant la condition de Cauchy C. Son intervalle de définition est ouvert.
On déduit de ce théorème les corollaires suivants :

Sous les mêmes conditions que le théorème, toute solution de l'équation différentielle (1) respectant la condition C et définie sur un intervalle est une restriction de la solution maximale.
Les graphes des solutions maximales forment une partition de Ω.
Lorsque l'ouvert Ω est un produit I × E, où I est un intervalle ouvert de R, le théorème de Cauchy-Lipschitz global apporte un complément :

Si f est lipschitzienne par rapport à la deuxième variable (localement par rapport à la première seulement), alors toute solution maximale est globale (c'est-à-dire définie sur I tout entier).
L'expression théorème de Cauchy-Lipschitz est aussi utilisée pour désigner d'autres résultats plus évolués. Il est possible de considérer non pas seulement la fonction solution s mais aussi la fonction qui au couple (x, t) associe l'image de t par la solution égale à x à un instant initial t0. On obtient une fonction appelée flot. Si la fonction f est de classe Cp, le flot l'est aussi. Il est aussi possible d'étudier la régularité des solutions si la fonction f dépend d'un paramètre.
D'autres résultats, qui ne portent pas le nom de théorème de Cauchy-Lipschitz sont présentés dans le paragraphe : Généralisations.


### Approche intuitive
L'idée à la base du théorème est la suivante : Certaines équations différentielles possèdent des solutions stationnaires. Par exemple, l'équation 
  
    
      
        d
        y
        
          /
        
        d
        t
        =
        a
        y
      
    
    {\displaystyle dy/dt=ay}
  
 possède la solution stationnaire 
  
    
      
        y
        (
        t
        )
        =
        0
      
    
    {\displaystyle y(t)=0}
  
 pour la condition initiale 
  
    
      
        y
        (
        0
        )
        =
        0
      
    
    {\displaystyle y(0)=0}
  
. Pour une autre condition initiale 
  
    
      
        y
        (
        0
        )
        =
        
          y
          
            0
          
        
        ≠
        0
      
    
    {\displaystyle y(0)=y_{0}\neq 0}
  
, la solution stationnaire est atteinte au bout d'un temps infini, et la solution est unique. Cependant, si la solution stationnaire peut être atteinte en un temps fini, l'unicité est violée. Considérons par exemple l'équation

  
    
      
        d
        y
        
          /
        
        d
        t
        =
        a
        
          y
          
            2
            
              /
            
            3
          
        
        .
      
    
    {\displaystyle dy/dt=ay^{2/3}.}
  

Pour la condition initiale 
  
    
      
        y
        (
        0
        )
        =
        0
      
    
    {\displaystyle y(0)=0}
  
, nous pouvons avoir la solution 
  
    
      
        y
        (
        t
        )
        =
        0
      
    
    {\displaystyle y(t)=0}
  
 ou la solution

  
    
      
        y
        (
        t
        )
        =
        
          
            {
            
              
                
                  (
                  a
                  t
                  
                    /
                  
                  3
                  
                    )
                    
                      3
                    
                  
                
                
                  t
                  <
                  0
                
              
              
                
                  0
                
                
                  t
                  ≥
                  0.
                
              
            
            
          
        
      
    
    {\displaystyle y(t)={\begin{cases}(at/3)^{3}&t<0\\0&t\geq 0.\end{cases}}}
  

On peut noter que la fonction 
  
    
      
        f
        (
        y
        )
        =
        
          y
          
            2
            
              /
            
            3
          
        
      
    
    {\displaystyle f(y)=y^{2/3}}
  
 a une pente infinie en 
  
    
      
        y
        =
        0
      
    
    {\displaystyle y=0}
  
 et ne respecte pas la condition de continuité de Lipschitz. La condition de continuité de Lipschitz élimine ce type d'équations différentielles.


### Déterminisme et chaos
Pour le mathématicien Vladimir Arnold, « Les équations différentielles sont le pivot de la conception scientifique du monde ». En physique, les équations différentielles vérifient de manière très générale les hypothèses du théorème de l'article, il en est ainsi de celle qui régit la loi de la gravitation universelle ou l'interaction électromagnétique. Ces équations sont intrinsèquement déterministes et c'est l'une des conséquences du théorème de Cauchy-Lipschitz. Ce déterminisme, c'est-à-dire la capacité de la physique à prévoir, par exemple la trajectoire des planètes, était connu bien avant le théorème. La célèbre citation de Voltaire l'atteste : « Mais en apercevant l’ordre, l'artifice prodigieux, les lois mécaniques et géométriques qui règnent dans l’univers […], je suis saisi d’admiration et de respect. ». Il en déduisait que « le Chaos est précisément l'opposé de toutes les lois de la nature ». Curieusement, les analyses fines de la fin du XIXe siècle démontrent la pertinence du caractère déterministe de ces lois physiques et réduisent cependant cette idée de Voltaire à l'état de chimère.

Les travaux mathématiques de cette époque montrent effectivement qu'une trajectoire, suivant une loi physique de cette nature, est unique et parfaitement déterminée. En revanche, les études de Poincaré sur la stabilité du système solaire mettent en évidence la faute de raisonnement. Ce mathématicien précise : 
« Une cause très petite, qui nous échappe, détermine un effet considérable que nous ne pouvons pas ne pas voir, et alors nous disons que cet effet est dû au hasard. Si nous connaissions exactement les lois de la nature et la situation de l'univers à l'instant initial, nous pourrions prédire exactement la situation de ce même univers à un instant ultérieur. Mais, lors même que les lois naturelles n'auraient plus de secret pour nous, nous ne pourrons connaître la situation initiale qu'approximativement. Si cela nous permet de prévoir la situation ultérieure avec la même approximation, c'est tout ce qu'il nous faut, nous disons que le phénomène a été prévu, qu'il est régi par des lois ; mais il n'en est pas toujours ainsi, il peut arriver que de petites différences dans les conditions initiales en engendrent de très grandes dans les phénomènes finaux ; une petite erreur sur les premières produirait une erreur énorme sur les derniers. La prédiction devient impossible et nous avons le phénomène fortuit. »

Le théorème de Cauchy-Lipschitz indique bien qu'une prévision parfaite est possible, mais uniquement sous réserve de connaître parfaitement la condition initiale. Certains systèmes dynamiques sont rapidement imprévisibles pour la raison qu'indique Poincaré. On utilise maintenant le terme de chaos pour décrire cette situation. Elle ne se produit pas systématiquement, le théorème de Poincaré-Bendixson précise un contexte où elle ne peut avoir lieu, mais les hypothèses sont restrictives et le résultat peu généralisable.


## Vocabulaire
Plusieurs termes ne sont pas nécessairement intuitifs, ce paragraphe en propose les définitions.


### Vocabulaire générique aux équations différentielles
Une partie du vocabulaire utilisée pour exprimer le théorème est générique aux équations différentielles ou fonctionnelles. Le terme espace de Banach désigne un espace vectoriel normé et complet, c'est-à-dire que pour la métrique associée à la norme, toute suite de Cauchy converge. Un exemple simple est un espace vectoriel réel normé de dimension finie. Pour une compréhension plus facile, le lecteur peut imaginer que E désigne l'ensemble des nombres réels. Dans l'article, tous les espaces de Banach considérés sont des espaces réels.

Une solution, ou courbe intégrale de l'équation (1) est une fonction d'un intervalle de R dans E dont le graphe est inclus dans Ω, et qui est solution de l'équation (1).


### Équation différentielle autonome du premier ordre
Il existe un cas particulier, celui où la fonction f est définie sur un ouvert de E et non pas sur un ouvert de R×E. L'équation (1) s'écrit alors

  
    
      
        
          
            
              x
              ˙
            
          
        
        (
        t
        )
        =
        f
        (
        x
        (
        t
        )
        )
        .
      
    
    {\displaystyle {\dot {x}}(t)=f(x(t)).}
  

Ω désigne un ouvert de E et l'équation est dite autonome.


### Vocabulaire spécifique
D'autres termes sont plus spécifiquement utilisés dans le contexte du théorème de Cauchy-Lipschitz. L'équation différentielle (1) possède généralement plusieurs solutions. Pour cette raison, on ajoute parfois une condition particulière :

Une condition de Cauchy C est un couple (t0, x0) élément de Ω. Une fonction s définie sur un intervalle de R et à valeurs dans E solution de l'équation (1) vérifie la condition de Cauchy C si l'intervalle de définition de s contient t0 et s(t0) = x0.
Remarque :  À la place de l'expression condition de Cauchy, on parle parfois de condition initiale. Les deux expressions sont synonymes.
Résoudre le problème de Cauchy consiste à trouver une solution de l'équation (1) vérifiant la condition de Cauchy C. Cette définition est un peu intuitive. Si l'équation différentielle modélise un courant parcourant une étendue d'eau, c'est-à-dire qu'à l'instant t et au point x le courant est égal à f(t,x), il est possible de poser dans l'eau un bouchon à l'instant t0 au point x0, la trajectoire du bouchon est la solution de l'équation. Il existe bien une solution pour chaque couple (t0, x0).
Une condition de Cauchy ne suffit pas à rendre la solution unique. Considérons l'ensemble S des solutions de l'équation (1) définies sur un intervalle. Cet ensemble est muni d'une relation d'ordre (partiel). Une solution s1 de S est plus petite qu'une solution s2 de S lorsque s2 est un prolongement de s1.

Une courbe intégrale de S est dite maximale si elle est maximale pour la relation d'ordre ci-dessus, c'est-à-dire qu'elle ne peut pas être prolongée en une solution définie sur un intervalle strictement plus grand.
Lorsque f ne remplit pas les conditions du théorème de Cauchy-Lipschitz, en utilisant une forme ou une autre de l'axiome du choix, on montre que toute solution dans S est majorée par une solution maximale (de). Par conséquent, tout problème de Cauchy possède une solution maximale, et s'il n'en possède qu'une alors elle est maximum (parmi les solutions de ce problème). L'article sur le théorème de Cauchy-Peano-Arzelà donne des exemples montrant que cette unicité n'est généralement pas garantie.
Lorsque tout problème de Cauchy associé à f possède une unique solution locale, on montre directement (sans axiome du choix) qu'il possède une solution maximum. Pour cela, il suffit que f soit suffisamment régulière :

La fonction f est dite :
lipschitzienne par rapport à la deuxième variable sur une partie W de Ω s'il existe une constante k telle que :
  
    
      
        (
        t
        ,
        
          x
          
            1
          
        
        )
        ,
        (
        t
        ,
        
          x
          
            2
          
        
        )
        ∈
        W
        ⇒
        ‖
        f
        (
        t
        ,
        
          x
          
            1
          
        
        )
        −
        f
        (
        t
        ,
        
          x
          
            2
          
        
        )
        ‖
        ≤
        k
        ‖
        
          x
          
            1
          
        
        −
        
          x
          
            2
          
        
        ‖
      
    
    {\displaystyle (t,x_{1}),(t,x_{2})\in W\Rightarrow \|f(t,x_{1})-f(t,x_{2})\|\leq k\|x_{1}-x_{2}\|}
  
 ;
localement lipschitzienne par rapport à la deuxième variable si tout point de Ω possède un voisinage sur lequel f est lipschitzienne par rapport à la deuxième variable ;
Lorsque l'ouvert Ω est un produit I × V, où I est un intervalle ouvert de R et V un ouvert de E, une solution de (1) est dite globale si elle est définie sur I tout entier. Une solution globale est évidemment maximale. Une hypothèse supplémentaire assure la réciproque :

La fonction f est dite lipschitzienne par rapport à la deuxième variable localement seulement par rapport à la première variable si tout point de I possède un voisinage J tel que, sur J×V, f soit lipschitzienne par rapport à la deuxième variable.


## Exemples et usages
### Champ constant
Dans cet exemple E est égal au plan réel R2. L'équation différentielle est autonome ; elle s'écrit x' = f(x) où f est la fonction constante de valeur le vecteur v de coordonnées (1, 0), sur Ω = ]–3, 3[2. Ce domaine de f est représenté par le carré bleu de la figure de droite. Une méthode de représentation de la fonction f consiste à dessiner les vecteurs f(x), où x est un point du domaine, en plaçant leurs origines au point x.
On recherche les solutions vérifiant la condition de Cauchy (0, x0), ici les coordonnées de x0 sont (–2, 2). Une solution s est de la forme :

  
    
      
        s
        (
        t
        )
        =
        
          x
          
            0
          
        
        +
        t
        v
        .
      
    
    {\displaystyle s(t)=x_{0}+tv.}
  

La solution ne peut quitter le domaine Ω, en conséquence, la variable t prend nécessairement ses valeurs dans ]–1, 5[. La fonction f est bien continue et k-lipschitzienne en x, avec dans ce cas particulier k = 0. Le théorème de Cauchy-Lipschitz garantit qu'il n'existe qu'une unique courbe intégrale vérifiant une condition de Cauchy précise, comme ici celle donnée par le couple (0, x0). Le graphe de cette solution est représenté en vert sur la figure.


### Équation linéaire
Le théorème de Cauchy-Lipschitz global s'applique : un problème de Cauchy associé à une équation différentielle linéaire possède une solution globale, dont toute autre solution est restriction.


### Étude d'une courbe intégrale
Le théorème de Cauchy-Lipschitz peut être vu comme un outil permettant l'étude d'une courbe intégrale. Illustrons-le par l'exemple d'une fonction logistique, définie par le problème de Cauchy suivant :

  
    
      
        
          
            
              
                d
              
              x
            
            
              
                d
              
              t
            
          
        
        =
        x
        (
        1
        −
        x
        )
        
        
          
            avec
          
        
        
        x
        (
        0
        )
        =
        
          
            1
            2
          
        
        .
      
    
    {\displaystyle {\frac {\mathrm {d} x}{\mathrm {d} t}}=x(1-x)\quad {\textrm {avec}}\quad x(0)={\frac {1}{2}}.}
  

L'équation différentielle s'écrit aussi x' = f(x) avec f(x) = x(1 – x). La fonction f est une fonction polynomiale de degré 2 de zéros 0 et 1 et positive entre ses racines. Le théorème de Cauchy-Lipschitz s'applique car f est continûment dérivable donc localement lipschitzienne : pour toute condition initiale, le problème de Cauchy possède une unique solution maximale. En particulier, toute solution de l'équation différentielle qui prend la valeur 0 ou 1 est constamment égale à cette valeur sur son intervalle de définition.
Soient s la solution maximale du problème de Cauchy ci-dessus et ]a, b[ son intervalle de définition (avec –∞ ≤ a < 0 < b ≤ +∞). D'après le théorème des valeurs intermédiaires, s(]a, b[) est un intervalle. Il contient s(0) = 1/2 donc, d'après la remarque précédente, ne contient ni 0 ni 1, si bien qu'il est inclus dans ]0, 1[.
Le fait que s prenne ses valeurs entre 0 et 1 montre que sa dérivée est strictement positive. L'application s est donc strictement croissante sur ]a, b[. Elle possède donc en b une limite à gauche c ∈ [1/2, 1] et sa dérivée a pour limite c(1 – c), ce qui prouve que b = +∞ (sans quoi, on pourrait prolonger s en une solution définie sur ]a, b], ce qui contredirait la maximalité de s). Comme s est bornée, la seule valeur limite possible en +∞ de sa dérivée est 0, ce qui montre que c = 1. Le même raisonnement montre que a = –∞ (le domaine de définition de s est donc égal à R) et que la limite de s en –∞ est 0
Enfin, on remarque que la fonction qui à t associe 1 – s(–t) est une solution du même problème de Cauchy donc est égale à s, ce qui signifie que la fonction s – 1/2 est impaire. L'équation montre que le point d'abscisse 0 est l'unique point d'inflexion de la courbe et que la dérivée de s en 0 est égale à 1/4, ce qui permet d'établir son graphe, représenté à droite.
Dans le cas de l'exemple choisi, on peut résoudre l'équation différentielle et utiliser les méthodes classiques pour l'étude de la courbe intégrale, mais ce n'est pas toujours possible. Une équation différentielle n'a pas nécessairement des solutions s'exprimant sous la forme d'une expression algébrique construite à l'aide des fonctions élémentaires.


### Théorème de Poincaré-Bendixson
L'usage du théorème de Cauchy-Lipschitz ne se limite pas à la résolution pratique des équations. Il sert aussi d'outil théorique, par exemple pour mieux comprendre le comportement qualitatif d'une équation différentielle. On peut considérer le cas d'une équation différentielle autonome dans R2 ayant une solution p périodique. La courbe intégrale est un lacet simple, c'est-à-dire que son graphe forme une boucle sans point double. Soit maintenant une condition de Cauchy C correspondant à un point à l'intérieur de la boucle. Le théorème de Cauchy-Lipschitz indique que la courbe intégrale maximale s vérifiant C ne pourra jamais traverser la boucle, si la fonction f définissant l'équation est localement lipschitzienne. La courbe s est donc bornée et si l'on suppose que le domaine de f contient l'intérieur de la boucle, la courbe s ne s'approche jamais trop du bord du domaine. Ceci suffit à démontrer que le domaine de définition de s est R tout entier.
Le théorème de Poincaré-Bendixson permet d'aller plus loin. Il indique que soit la courbe s est convergente, soit son comportement s'approche de plus en plus d'une fonction périodique. Cette configuration interdit les trajectoires chaotiques.
L'usage fait ici du théorème permet de comprendre qualitativement le comportement d'une courbe intégrale. Pour l'étude d'équations différentielles plus complexes, par exemple certains systèmes dynamiques, cette approche est indispensable. Il n'est en effet plus toujours possible de résoudre explicitement ces équations, ni même d'approcher leurs solutions sur des longues périodes de temps.


## Fragments d'histoire
### Origines
L'origine de la question traitée par le théorème est ancienne, elle porte initialement le nom de « problème inverse des tangentes ». À chaque point de l'espace, une droite est associée, le problème à résoudre est de trouver la courbe ayant, comme tangente en chaque point, une de ces droites, ce qui correspond en termes modernes à une équation différentielle autonome. Kepler est un initiateur de cette question[source insuffisante] dans une étude sur la contenance d'un tonneau de vin en 1615. Si cette question est abordée par des mathématiciens comme Descartes, Fermat ou Roberval, qui développent cette approche dans des cas particuliers durant le XVIIe siècle, le progrès essentiel est l'œuvre de Newton et Leibniz avec la découverte du calcul infinitésimal. Newton cherche surtout à obtenir un résultat à l'aide d'une série, Leibnitz recherche aussi des solutions exactes, sous forme de primitives de fonctions connues.
Le siècle suivant est l'objet d'une systématisation de l'étude. Dans un premier temps « Des trésors d'ingéniosité ont été dépensés pour ramener à des quadratures d'innombrables équations différentielles particulières et comme l'écrit Paul Painlevé : la vague s'arrêta quand tout ce qui était intégrable, dans les problèmes naturels fut intégré. ». Euler, en 1768 étudie la manière d'approcher une solution. Il étudie le cas particulier x'(t) = f(t) et cherche une solution sur un intervalle [a, b]. Pour ce faire, il partitionne l'intervalle à l'aide d'une suite a0 = a, a1, …, an = b et, si c est élément de l'intervalle [ai, ai+1], il propose l'approximation suivante :

  
    
      
        x
        (
        c
        )
        ≈
        x
        (
        a
        )
        +
        
          ∑
          
            j
            =
            0
          
          
            i
          
        
        (
        
          a
          
            j
            +
            1
          
        
        −
        
          a
          
            j
          
        
        )
        f
        (
        
          a
          
            j
          
        
        )
        .
      
    
    {\displaystyle x(c)\approx x(a)+\sum _{j=0}^{i}(a_{j+1}-a_{j})f(a_{j}).}
  

Dans le cas plus général de l'équation x'(t) = f(t, x(t)), il utilise la même méthode, qui donne :

  
    
      
        
          x
          
            0
          
        
        =
        x
        (
        a
        )
        ,
        
        
          x
          
            1
          
        
        ≈
        
          x
          
            0
          
        
        +
        (
        
          a
          
            1
          
        
        −
        
          a
          
            0
          
        
        )
        f
        (
        
          a
          
            1
          
        
        ,
        
          x
          
            0
          
        
        )
        ,
        
        ⋯
        ,
        
        x
        (
        
          a
          
            i
          
        
        )
        ≈
        
          x
          
            0
          
        
        +
        
          ∑
          
            j
            =
            0
          
          
            i
          
        
        (
        
          a
          
            j
            +
            1
          
        
        −
        
          a
          
            j
          
        
        )
        f
        (
        
          a
          
            j
          
        
        ,
        
          x
          
            j
            −
            1
          
        
        )
        .
      
    
    {\displaystyle x_{0}=x(a),\;x_{1}\approx x_{0}+(a_{1}-a_{0})f(a_{1},x_{0}),\;\cdots ,\;x(a_{i})\approx x_{0}+\sum _{j=0}^{i}(a_{j+1}-a_{j})f(a_{j},x_{j-1}).}
  

Euler ne se pose pas la question de la convergence si le découpage est de plus en plus fin.


### Les apports de Cauchy et de Lipschitz
Augustin Louis Cauchy (1789 - 1857) établit les premiers résultats généraux. Dans un premier temps, il précise sa méthode : « Dans mes leçons données à l'École Polytechnique, comme dans la plupart des ouvrages ou mémoires que j'ai publiés sur le calcul intégral, j'ai cru devoir renverser cet ordre et placer en premier lieu la recherche, non pas des intégrales générales, mais des particulières ; en sorte que la détermination des constantes ou des fonctions arbitraires [pour l'équation aux dérivées partielles] ne fût plus séparée de la recherche des intégrales. » Ce que Cauchy décrit ici est la démarche formalisée par le problème de Cauchy.
À l'aide de la formalisation de la notion de limite par Bolzano, son approche lui permet d'aller plus loin. La question x' = f(x) sur l'intervalle [a, b], si on lui ajoute la condition de Cauchy x(a) = x0, change de nature. Ce n'est plus la recherche d'une primitive de f, mais le calcul d'une intégrale. Il montre que, avec les notations du paragraphe précédent, si ai + 1 – ai tend vers 0 et si f est continue, l'approximation converge. Plus précisément, il applique cette démarche à l'équation x' (t) = f(t, x(t)) dans le cas où f, ainsi que sa différentielle partielle par rapport à la seconde variable, sont continues et bornées : avec les mêmes conditions, on obtient encore une convergence. C'est la première version du théorème, démontrée par Cauchy dès 1820. En 1835, sa méthode est généralisée aux fonctions holomorphes.
Rudolf Lipschitz (1832 - 1903), sans manifestement connaître la teneur des travaux de Cauchy, démontre un théorème « essentiellement équivalent » au théorème local de son prédécesseur, en introduisant la condition qui porte maintenant son nom. Cauchy utilisait cette condition en la déduisant de son hypothèse sur la différentielle partielle de f, à l'aide de l'une de ses découvertes, qu'il affectionne particulièrement : le théorème des accroissements finis.


### Le formalisme moderne
La fin du XIXe siècle voit une profonde mutation du théorème, tant dans la manière de le démontrer que dans l'enrichissement de son contenu, nécessaire pour mieux comprendre l'équation différentielle.
À la suite des travaux de Fuchs, les objectifs deviennent plus ambitieux. Paul Painlevé et Émile Picard s'intéressent au cas général des équations différentielles de premier et second ordre et donc à leurs singularités. Les objectifs d'Henri Poincaré sont encore plus généraux. À la suite d'une étude sur la stabilité du système solaire, il cherche à établir une théorie générale qui prendra le nom de système dynamique. Devant l'impossibilité d'établir des solutions explicites, Poincaré fonde les bases d'une théorie qualitative. Les objectifs sont, dans un premier temps l'étude des singularités, puis le comportement asymptotique (c'est-à-dire l'étude du comportement une fois le système stabilisé) et la sensibilité à la condition initiale.
Pour atteindre ces nouveaux objectifs, la formulation par Cauchy ou Lipschitz du théorème devient insuffisante. On cherche maintenant des éléments de réponse globaux sur les courbes intégrales et non plus uniquement un résultat local. Les questions sur la régularité de la solution deviennent essentielles. Enfin, on cherche à déterminer la nature de la modification de la courbe intégrale en fonction d'une modification de la condition initiale ou d'un paramètre de l'équation. Les méthodes pour y arriver diffèrent radicalement de l'approche de Cauchy qui étudiait le comportement limite de la fonction polygonale imaginée par Euler. Elles se rapprochent de l'analyse fonctionnelle, le contexte de l'étude étant maintenant un espace de fonctions disposant de propriétés géométriques. Le théorème à la source de la démonstration est celui du point fixe. Sa forme initiale est l'œuvre de Picard. Ce théorème est maintenant vu comme une propriété générale d'espaces vectoriels particuliers, formalisés par Stefan Banach. Son application au théorème de l'article est l'œuvre du mathématicien finlandais Ernst Lindelöf en 1894. Pour cette raison, le théorème nommé en France d'après Cauchy et Lipschitz prend le nom, en anglais, de Picard–Lindelöf theorem.


## Généralisations
Il existe de nombreuses généralisations du théorème de Cauchy-Lipschitz. Les techniques utilisées pour sa démonstration permettent d'aller plus loin dans l'analyse des solutions d'une équation différentielle. Une rapide analyse montre qu'une solution est continûment dérivable, mais rien n'indique ce qui se produit si la fonction f est plus régulière. La branche des mathématiques étudiant un système dynamique se pose deux questions essentielles : quel est le comportement d'une courbe intégrale si t s'approche des bornes du domaine de définition ? Et quelle est la sensibilité à la condition initiale ? c'est-à-dire que se passe-t-il si x0 subit une petite modification.
Enfin, l'article suppose que l'ensemble E est un espace de Banach, ce qui n'est pas l'unique cas étudié. Dans certaines situations, il est utile de considérer E comme une variété différentielle. Les théorèmes se transposent aisément dans ce nouvel univers.


### Flot
L'un des objectifs de l'étude des systèmes dynamiques est celle de la sensibilité à la condition de Cauchy. Le vocabulaire et la représentation géométrique sont un peu différents de ce qui a été utilisé jusqu'à présent. Pour en comprendre l'origine, le plus simple est d'imaginer que Ω est un plan d'eau et que R représente le temps. Le plan d'eau est agité par un courant, représenté par la fonction f, appelé champ de vecteurs. En dimension 2, on représente ce champ de vecteurs en associant à certains points x de E une représentation graphique du vecteur f(x) (si le champ de vecteurs ne dépend pas du temps t), à l'image de la figure de droite. Une courbe intégrale satisfaisant à la condition de Cauchy C peut s'imaginer comme la trajectoire d'un bouchon placé dans l'eau à l'instant t0 à la position x0. Pour connaître d'un seul coup toutes les solutions de l'équation différentielle, il suffit de connaître le mouvement de la surface de l'eau, appelé flot, coulée ou encore courant.
Avec ce concept, le théorème de Cauchy-Lipschitz prend une nouvelle forme :

Si f est continue et localement lipschitzienne par rapport à la deuxième variable, le flot est continu, dérivable par rapport à la première variable et localement lipschitzien par rapport à la deuxième.
Il est possible d'aller plus loin si f est plus régulière.

Si f est de classe Cp, le flot l'est aussi.
Cette forme du théorème est plus forte que la précédente. Elle montre la régularité des solutions sous l'action d'une petite modification de la condition de Cauchy. Cet énoncé garantit l'absence d'une certaine forme de chaos. Les bifurcations sont impossibles et, sous réserve que t ne grandisse pas trop, les trajectoires restent proches si la condition de Cauchy est peu modifiée. Néanmoins, si le champ de vecteurs évolue avec le temps ou si E est de dimension strictement supérieure à deux, une autre forme de chaos peut s'installer.
Les formes précises des théorèmes et les démonstrations sont données dans l'article détaillé.


### Théorème de Cauchy-Peano-Arzelà
À la fin du XIXe siècle, Giuseppe Peano, un mathématicien italien, s'interroge sur les généralisations possibles du théorème de Cauchy-Lipschitz sous des hypothèses plus faibles : que se passe-t-il si la fonction f, définissant l'équation (1), reste continue mais n'est pas localement-lipschitzienne par rapport à la deuxième variable ?
Une situation physique correspondant à ce cas est une bille roulant sans frottement sur l'arête d'un toit en V inversé aux branches de pente constante. Notons x la coordonnée transversale (perpendiculairement à l'arête) et y la coordonnée longitudinale (parallèlement à l'arête). En appliquant les lois de la mécanique, on obtient une équation différentielle autonome du premier ordre dans R2, si l'on ne considère que les solutions qui penchent vers les valeurs positives de x :

  
    
      
        
          
            {
            
              
                
                  
                    
                      
                        
                          d
                        
                        x
                      
                      
                        
                          d
                        
                        t
                      
                    
                  
                  =
                  λ
                  
                    
                      x
                    
                  
                
              
              
                
                  
                    
                      
                        
                          d
                        
                        y
                      
                      
                        
                          d
                        
                        t
                      
                    
                  
                  =
                  μ
                  .
                
              
            
            
          
        
      
    
    {\displaystyle {\begin{cases}{\frac {\mathrm {d} x}{\mathrm {d} t}}=\lambda {\sqrt {x}}\\{\frac {\mathrm {d} y}{\mathrm {d} t}}=\mu .\end{cases}}}
  

Si t0 désigne un réel positif, toutes les solutions suivantes vérifient l'équation précédente, pour la même condition de Cauchy, à savoir que le solide est à l'instant 0, à la position (0,0) :

  
    
      
        
          
            {
            
              
                
                  x
                  =
                  0
                  
                     si 
                  
                  t
                  ≤
                  
                    t
                    
                      0
                    
                  
                  
                  ,
                  
                  x
                  =
                  
                    
                      
                        λ
                        
                          2
                        
                      
                      4
                    
                  
                  (
                  t
                  −
                  
                    t
                    
                      0
                    
                  
                  
                    )
                    
                      2
                    
                  
                  
                     sinon
                  
                
              
              
                
                  y
                  =
                  μ
                  t
                  .
                
              
            
            
          
        
      
    
    {\displaystyle {\begin{cases}x=0{\text{ si }}t\leq t_{0}\,,\;x={\frac {\lambda ^{2}}{4}}(t-t_{0})^{2}{\text{ sinon}}\\y=\mu t.\end{cases}}}
  

D'autres exemples sont donnés dans l'article détaillé.
Ici, on remarque que si l'unicité n'est plus valable, l'existence reste vraie. De manière plus générale, si la fonction f de l'équation (1) est continue et bornée et si E est de dimension finie, l'existence d'une solution est garantie. Ce résultat est appelé théorème de Cauchy-Peano-Arzelà.
Quand on lui adjoint le critère d'unicité d'Osgood, on obtient la même conclusion que Cauchy-Lipschitz, sous des hypothèses plus faibles.


### Extension aux équations aux dérivées partielles
La généralisation du théorème pour les équations aux dérivées partielles suppose des hypothèses plus fortes pour un résultat plus faible. La fonction f définissant l'équation doit être analytique ainsi que les conditions aux limites qui remplacent celle de Cauchy. Aucune information n'est fournie dans le cas d'une modification des conditions aux limites ou d'un paramètre. Ce théorème est appelé : théorème de Cauchy-Kowalevski.


## Démonstrations
Les preuves du théorème utilisent de façon fondamentale un théorème d'analyse fonctionnelle : ou bien le théorème d'Ascoli, ou bien le théorème du point fixe de Banach. C'est le second qui sera invoqué ici.


### Préambule
Pour comprendre le mécanisme de la preuve, reprenons l'exemple de l'équation différentielle définissant la courbe logistique, avec la condition de Cauchy C :

  
    
      
        
          
            
              
                d
              
              x
            
            
              
                d
              
              t
            
          
        
        =
        x
        −
        
          x
          
            2
          
        
        
        
          avec la condition 
        
        C
         
        :
        
        x
        (
        0
        )
        =
        
          
            1
            2
          
        
      
    
    {\displaystyle {\frac {\mathrm {d} x}{\mathrm {d} t}}=x-x^{2}\quad {\text{avec la condition }}C~:\quad x(0)={\frac {1}{2}}}
  
.
On définit une suite de fonctions polynomiales (un) par récurrence :

  
    
      
        ∀
        t
        ∈
        
          R
        
        
        
          u
          
            0
          
        
        (
        t
        )
        =
        
          
            1
            2
          
        
        
        
          et
        
        
        
          u
          
            n
            +
            1
          
        
        (
        t
        )
        =
        
          
            1
            2
          
        
        +
        
          ∫
          
            0
          
          
            t
          
        
        [
        
          u
          
            n
          
        
        (
        τ
        )
        −
        
          u
          
            n
          
          
            2
          
        
        (
        τ
        )
        ]
        
        
          d
        
        τ
        .
      
    
    {\displaystyle \forall t\in \mathbb {R} \quad u_{0}(t)={\frac {1}{2}}\quad {\text{et}}\quad u_{n+1}(t)={\frac {1}{2}}+\int _{0}^{t}[u_{n}(\tau )-u_{n}^{2}(\tau )]\,\mathrm {d} \tau .}
  

Sur l'intervalle [–5/2, 5/2], la suite (un) converge uniformément. Sa limite est un point fixe de la fonction φ qui à une fonction continue f de [–5/2, 5/2] dans R associe la fonction φf :

  
    
      
        
          φ
          
            f
          
        
        (
        t
        )
        =
        
          
            1
            2
          
        
        +
        
          ∫
          
            0
          
          
            t
          
        
        [
        f
        (
        τ
        )
        −
        
          f
          
            2
          
        
        (
        τ
        )
        ]
        
        
          d
        
        τ
      
    
    {\displaystyle \varphi _{f}(t)={\frac {1}{2}}+\int _{0}^{t}[f(\tau )-f^{2}(\tau )]\,\mathrm {d} \tau }
  
.

En dérivant l'égalité φu = u, on vérifie que u est bien une solution de l'équation différentielle étudiée et par construction u(0) = 1/2. Cette démarche est celle de la démonstration de Lindelöf. Il constate que si l'intervalle est bien choisi, la fonction φ vérifie la condition de Lipschitz avec un coefficient strictement plus petit que 1, elle est donc contractante, ce qui permet de faire usage du théorème du point fixe. Comme une application contractante n'admet qu'un unique point fixe, l'unicité d'une solution locale est démontrée.
Cette méthode permet de trouver localement une solution. En revanche, pour la valeur 3, la suite (un) diverge. Cependant rien n'empêche de répéter la même démarche avec les deux conditions de Cauchy (5/2, u(5/2)) et (–5/2, u(–5/2)), il devient ainsi possible de prolonger la solution. Quitte à réitérer la démarche une infinité de fois, on finit par obtenir une solution maximale.
L'intérêt de la démarche est surtout théorique. Dans le cas particulier de l'exemple, il est aisé d'intégrer directement l'équation différentielle. Dans le cas général, il existe des méthodes plus rapides pour obtenir une approximation de la solution, comme celle d'Euler décrite dans la partie « Fragments d'histoire » ou encore celle de Runge-Kutta. En revanche, il est possible de démontrer dans le cadre général du théorème de l'article qu'une fonction construite de la même manière que φ est contractante, ce qui montre l'existence et l'unicité d'un point fixe, solution de l'équation différentielle. Cette démarche, caractéristique de l'analyse fonctionnelle, permet donc de démontrer un résultat plus fort que celui qu'avait démontré Cauchy.


### Solution locale
L'objectif de ce paragraphe est d'expliquer la méthode permettant de montrer que l'équation (1) admet localement une unique solution satisfaisant la condition de Cauchy C. Par translations, on se ramène sans peine au cas (t0, x0) = (0, 0). La condition C s'écrit alors : x(0) = 0.
On suppose que f est localement lipschitzienne par rapport à la seconde variable et continue. Comme Ω est un ouvert, il existe un réel strictement positif a tel que le produit du segment [–a, a] par la boule fermée B de centre 0 et de rayon a soit inclus dans Ω et tel que sur ce produit, f soit k-lipschitzienne par rapport à la deuxième variable pour un certain k, et de norme majorée par un réel m. Soit enfin b ≤ a un réel strictement positif, assez petit pour que bm ≤ a et bk < 1. Ce réel b permet de définir un espace F de fonctions sur lequel on construit une application Φ satisfaisant au théorème du point fixe :

F est l'espace des applications continues de [–b, b] dans B, muni de la norme de la convergence uniforme. Il est complet car B l'est
Φ est l'application qui à toute fonction u de F associe la fonction Φu définie par intégration (à valeurs dans E) :
  
    
      
        ∀
        t
        ∈
        [
        −
        b
        ,
        b
        ]
        
        
          Φ
          
            u
          
        
        (
        t
        )
        =
        
          ∫
          
            0
          
          
            t
          
        
        f
        (
        τ
        ,
        u
        (
        τ
        )
        )
        
          d
        
        τ
        .
      
    
    {\displaystyle \forall t\in [-b,b]\quad \Phi _{u}(t)=\int _{0}^{t}f(\tau ,u(\tau ))\mathrm {d} \tau .}
  

On démontre alors que la fonction Φ est à valeurs dans F, et qu'elle est bk-lipschitzienne donc contractante, puisque bk < 1.
Les hypothèses du théorème du point fixe étant réunies, on en déduit que Φ admet un unique point fixe dans F. D'après le théorème fondamental de l'analyse (étendu aux espaces de Banach), cela revient à dire que le problème de Cauchy admet sur [–b, b] une unique solution à valeurs dans B.
On montrerait de même que pour tout ε de [0, b], il existe sur [–ε, 0] et sur [0, ε], une unique solution à valeurs dans B du problème de Cauchy. Pour finir de démontrer l'unicité locale à gauche et à droite, on montre de plus qu'une solution sur [–ε, 0] ou [0, ε] du problème de Cauchy est nécessairement à valeurs dans B.
Tout ceci se résume en un premier théorème :

Une méthode pratique pour trouver le point fixe est de construire une suite (un) qui vérifie la relation de récurrence : un+1 = Φun, la suite converge nécessairement vers le point fixe. Cette technique est celle utilisée dans le préambule.


### Solution maximale
À l'exception de l'unicité de la solution locale, le paragraphe précédent résume le résultat démontré à l'époque de Cauchy. Il n'offre aucune information sur l'existence ou l'unicité d'une solution maximale vérifiant la condition de Cauchy C. Ces résultats, plus tardifs, sont les deux points suivants, qui se déduisent respectivement de l'existence et de l'unicité du théorème local ci-dessus :

Autrement dit, (cf. § « Vocabulaire spécifique ») : ce problème possède une solution non seulement maximale mais maximum, c'est-à-dire dont toute solution définie sur un intervalle est une restriction. Ainsi, grâce à l'unicité locale, on a redémontré directement à peu de frais que toute solution possède un prolongement en une solution maximale.


### Équation différentielle non autonome
Signalons accessoirement qu'à l'aide d'un jeu d'écriture, il est possible de généraliser le cas particulier des équations autonomes aux équations dépendantes du temps. Soit

  
    
      
        g
        :
        Ω
        →
        
          R
        
        ×
        E
        ,
         
        y
        ↦
        (
        1
        ,
        f
        (
        y
        )
        )
      
    
    {\displaystyle g:\Omega \to \mathbb {R} \times E,\ y\mapsto (1,f(y))}
  
.
En désignant par y0 le point (t0, x0) de Ω, on considère le problème de Cauchy suivant :

  
    
      
        (
        2
        )
        
        
          y
          ′
        
        =
        g
        (
        y
        )
        
        
          avec la condition de Cauchy 
        
        
          C
          
            2
          
        
        
        y
        (
        
          t
          
            0
          
        
        )
        =
        
          y
          
            0
          
        
        .
      
    
    {\displaystyle (2)\quad y'=g(y)\quad {\text{avec la condition de Cauchy }}C_{2}\quad y(t_{0})=y_{0}.}
  

La fonction g est localement lipschitzienne dès que f l'est. Cependant, le théorème général ne se déduit pas immédiatement de cette remarque car son énoncé ne suppose pas que f est localement lipschitzienne, mais seulement qu'elle est continue et localement lipschitzienne par rapport à sa deuxième variable.


## Solutions maximales et solutions globales
On suppose ici que l'ouvert Ω est de la forme I × E, où I est un intervalle ouvert de R. Toute solution globale (c'est-à-dire définie sur I tout entier) de l'équation (1) est évidemment maximale, mais la réciproque est fausse en général, comme le montre l'exemple de l'équation x' = x2. Divers théorèmes « d'échappement » ou « d'explosion », parfois joints au lemme de Grönwall, donnent des conditions suffisantes pour une telle réciproque, le plus connu étant le lemme des bouts, ou théorème d'explosion en temps fini. L'énoncé suivant, qui suffit par exemple dans la théorie des équations différentielles linéaires, se démontre directement :


## Notes et références
### Notes
### Références
### Ouvrages cités
Marcel Berger et Bernard Gostiaux, Géométrie différentielle : variétés, courbes et surfaces [détail des éditions] 
Srishti D. Chatterji, Cours d'analyse, vol. 3 : Équations différentielles ordinaires et aux dérivées partielles, PPUR, 1998 (lire en ligne)
Jean-Pierre Demailly, Analyse numérique et équations différentielles [détail des éditions]
(en) A. N. Kolmogorov et A. P. Yushkevich, Mathematics of the 19th Century, vol. 3 : Function Theory According to Chebyshev: Constructive Function Theory, Ordinary Differential Equations, Calculus of Variations, Theory of Finite Differences, Birkhäuser, 1998 (lire en ligne)
(en) S. Lang, Analysis, vol. II, Addison-Wesley, 1969
François Laudenbach, « Équations différentielles », dans Nicole Berline et Claude Sabbah, Aspects des systèmes dynamiques, Éditions de l’École polytechnique, 2009 (ISBN 978-2-73021560-2), p. 1-16
D. Leborgne, Calcul différentiel et géométrie, PUF, 1982 (ISBN 978-2-13-037495-4)
P. Malliavin, Géométrie différentielle intrinsèque, Hermann, 1972 (ISBN 978-2-7056-5696-6)
Jean Mawhin, « Problème de Cauchy pour les équations différentielles et théorie de l'intégration : influences mutuelles », Cahiers du séminaire d'histoire des mathématiques, vol. 9,‎ 1988 (lire en ligne)
Dany-Jack Mercier, L'épreuve d'exposé au CAPES mathématiques, vol. 2, Publibook, 2006 (ISBN 978-2-74833001-4)
F. Paulin, Topologie, analyse et calcul différentiel, École Normale Supérieure, 2008 (lire en ligne), chap. 7.5 (« Théorie de Cauchy-Lipschitz »), p. 255-272


## Voir aussi
### Articles connexes
Critère de Nagumo
Théorèmes d'unicité 


### Bibliographie
J.-M. Arnaudiès, Équations différentielles de fonctions de variable réelle ou complexe, Ellipses, 2000 (ISBN 978-2-7298-0045-1)
S. Lang, Analyse Réelle, Paris, InterEditions, 1977, 230 p. (ISBN 978-2-7296-0059-4)
Lev Pontriaguine, Équations différentielles ordinaires, Moscou, Éditions Mir, 1969
Nicolas Rouche et Jean Mawhin, Équations différentielles ordinaires (vol. 1 : Théorie générale, vol. 2 : Stabilité et solutions périodiques), Masson, Paris, 1973


### Liens externes
#### Histoire
A. L. Cauchy, Leçons de calcul différentiel et intégral, 1844, 26e leçon, p. 385-396


#### Mathématiques
A. Popier et O. Wintenberger, Équations différentielle, École Polytechnique 
N. Petit et P. Rouchon, Théorème de Cauchy, École des Mines de Paris, p. 203 
R. Rolland, Équations différentielles ordinaires, 

 Portail de l'analyse