Kurs:Analysis (Osnabrück 2021-2023)/Teil II/Vorlesung 56/latex

\setcounter{section}{56}






\zwischenueberschrift{Differential- und Integralgleichungen}

Mit dem Begriff des Integrals einer Kurve kann man Differentialgleichungen auch als Integralgleichungen schreiben.




\inputfaktbeweis
{Differentialgleichung/Stetig/Integralgleichung/Fakt}
{Lemma}
{}
{

\faktsituation {Es sei $V$ ein \definitionsverweis {endlichdimensionaler}{}{} \definitionsverweis {reeller Vektorraum}{}{,}
\mavergleichskette
{\vergleichskette
{I }
{ \subseteq }{\R }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} ein \definitionsverweis {reelles Intervall}{}{,}
\mavergleichskette
{\vergleichskette
{U }
{ \subseteq }{V }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} eine \definitionsverweis {offene Menge}{}{} und \maabbeledisp {f} {I\times U} {V } {(t,v)} {f(t,v) } {,} ein stetiges \definitionsverweis {Vektorfeld}{}{} auf $U$. Es sei
\mavergleichskette
{\vergleichskette
{ (t_0,w) }
{ \in }{I \times U }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} vorgegeben.}
\faktfolgerung {Dann ist eine \definitionsverweis {stetige Abbildung}{}{} \maabbeledisp {v} {J} {U } {t} {v(t) } {,} auf einem \definitionsverweis {Intervall}{}{}
\mavergleichskette
{\vergleichskette
{J }
{ \subseteq }{I }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} mit
\mavergleichskette
{\vergleichskette
{t_0 }
{ \in }{J }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} genau dann eine \definitionsverweis {Lösung des Anfangswertproblems}{}{  \zusatzklammer {insbesondere muss $v$ differenzierbar sein} {} {}}
\mathdisp {v'=f(t,v) \text{ und } v(t_0)=w} { , }
wenn $v$ die Integralgleichung
\mavergleichskettedisp
{\vergleichskette
{ v(t) }
{ =} { w + \int_{ t_0 }^{ t } f(s,v(s)) \, d s }
{ } { }
{ } { }
{ } { }
} {}{}{} erfüllt.}
\faktzusatz {}
\faktzusatz {}

}
{

\teilbeweis {}{}{}
{Es sei die Integralbedingung erfüllt. Dann ist
\mavergleichskettedisp
{\vergleichskette
{ v(t_0) }
{ =} { w }
{ } { }
{ } { }
{ } { }
} {}{}{} und aufgrund des Hauptsatzes der Infinitesimalrechnung gilt
\mavergleichskette
{\vergleichskette
{ v'(t) }
{ = }{ f(t,v(t)) }
{ }{ }
{ }{ }
{ }{ }
} {}{}{.} Insbesondere sichert die Integralbedingung, dass $v$ \definitionsverweis {differenzierbar}{}{} ist.}
{}\teilbeweis {}{}{}
{Wenn umgekehrt $v$ eine Lösung des Anfangswertproblems ist, so ist
\mavergleichskette
{\vergleichskette
{ v'(s) }
{ = }{ f(s,v(s)) }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} und daher
\mavergleichskettedisp
{\vergleichskette
{ w + \int_{ t_0 }^{ t } f(s,v(s)) \, d s }
{ =} { w + \int_{ t_0 }^{ t } v'(s) \, d s }
{ =} { w+ v(t)-v(t_0) }
{ =} { v(t) }
{ } { }
} {}{}{.}}
{}

}







\zwischenueberschrift{Der Satz von Picard-Lindelöf}

Wir kommen nun zum wichtigsten Existenz- und Eindeutigkeitssatz für die Lösungen von gewöhnlichen Differentialgleichungen.




\inputfaktbeweis
{Picard Lindelöf/Lokal Lipschitz/Lokale Existenz und Eindeutigkeit/Fakt}
{Satz}
{}
{

\faktsituation {Es sei $V$ ein \definitionsverweis {endlichdimensionaler}{}{} \definitionsverweis {reeller Vektorraum}{}{,}
\mavergleichskette
{\vergleichskette
{I }
{ \subseteq }{\R }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} ein \definitionsverweis {reelles Intervall}{}{,}
\mavergleichskette
{\vergleichskette
{U }
{ \subseteq }{V }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} eine \definitionsverweis {offene Menge}{}{} und \maabbeledisp {f} {I\times U} {V } {(t,v)} {f(t,v) } {,} ein \definitionsverweis {Vektorfeld}{}{} auf $U$.}
\faktvoraussetzung {Es sei vorausgesetzt, dass dieses Vektorfeld \definitionsverweis {stetig}{}{} sei und \definitionsverweis {lokal einer Lipschitz-Bedingung}{}{} genüge.}
\faktfolgerung {Dann gibt es zu jedem
\mavergleichskette
{\vergleichskette
{ (t_0,w) }
{ \in }{ I \times U }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} ein \definitionsverweis {offenes Intervall}{}{}
\mathbed {J} {mit}
{t_0 \in J \subseteq I} {}
{} {} {} {} derart, dass auf diesem Intervall eine eindeutige \definitionsverweis {Lösung für das Anfangswertproblem}{}{}
\mathdisp {v'=f(t,v) \text{ und } v(t_0)=w} { }
existiert.}
\faktzusatz {}
\faktzusatz {}

}
{

Nach Lemma 56.1 ist eine \definitionsverweis {stetige Abbildung}{}{} \maabbdisp {v} {J} {V } {} genau dann eine \definitionsverweis {Lösung des Anfangswertproblems}{}{,} wenn $v$ die Integralgleichung
\mavergleichskettedisp
{\vergleichskette
{ v(t) }
{ =} { w + \int_{ t_0 }^{ t } f(s,v(s)) \, d s }
{ } { }
{ } { }
{ } { }
} {}{}{} erfüllt. Wir wollen die Existenz und Eindeutigkeit einer Lösung für diese Integralgleichung unter Verwendung des Banachschen Fixpunktsatzes dadurch erweisen, dass wir für die Abbildung \zusatzklammer {man spricht von einem \stichwort {Funktional} {}} {} {}
\mathdisp {\psi \longmapsto (t \mapsto w + \int_{ t_0 }^{ t } f(s, \psi (s)) \, d s )} { }
einen Fixpunkt finden. Hierbei stehen links und rechts Abbildungen in $t$ \zusatzklammer {aus einem gewissen Teilintervall von $I$ mit Werten in $V$} {} {.} Die Fixpunkteigenschaft
\mavergleichskette
{\vergleichskette
{ H(\psi) }
{ = }{ \psi }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} bedeutet gerade, dass
\mavergleichskette
{\vergleichskette
{ \psi(t) }
{ = }{ w + \int_{t_0}^t f(s, \psi(s)) ds }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} ist. Um den Fixpunktsatz anwenden zu können müssen wir ein Definitionsintervall festlegen, und eine Metrik auf dem Abbildungsraum nach $V$ definieren, diesen metrischen Raum dann als \definitionsverweis {vollständig}{}{} und das Funktional als \definitionsverweis {stark kontrahierend}{}{} nachweisen. \teilbeweis {}{}{}
{Aufgrund der Voraussetzung über die lokale Lipschitz-Bedingung gibt es eine offene Umgebung
\mavergleichskettedisp
{\vergleichskette
{ (t_0,w) }
{ \in} { J' \times U { \left( w,\epsilon \right) } }
{ \subseteq} { I \times U }
{ } { }
{ } { }
} {}{}{} und ein
\mavergleichskette
{\vergleichskette
{ L }
{ \in }{ \R_{\geq 0} }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} mit
\mavergleichskettedisp
{\vergleichskette
{ \Vert {f(t,v)-f(t, \tilde{v}) } \Vert }
{ \leq} { L \Vert {v - \tilde{v}} \Vert }
{ } { }
{ } { }
{ } { }
} {}{}{} für alle \mathkor {} {t \in J'} {und} {v,\tilde{v} \in U { \left( w,\epsilon \right) }} {.} Durch Verkleinern der Radien können wir annehmen, dass der \definitionsverweis {Abschluss}{}{} von
\mathl{J' \times U { \left( w,\epsilon \right) }}{,} also das Produkt des abgeschlossenen Intervalls mit der abgeschlossenen Kugel, ebenfalls in
\mathl{I \times U}{} liegt. Aufgrund von Satz 36.12 gibt es ein
\mavergleichskette
{\vergleichskette
{ M }
{ \in }{ \R_+ }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} mit
\mathdisp {\Vert {f(t,v)} \Vert \leq M \text { für alle } (t,v) \in J' \times U { \left( w,\epsilon \right) }} { }
\zusatzklammer {da diese Beschränktheit auf dem Abschluss gilt} {} {.} Wir ersetzen nun $J'$ durch ein kleineres Intervall
\mavergleichskettedisp
{\vergleichskette
{ J }
{ =} { [t_0- \delta,t_0+ \delta ] }
{ \subseteq} { J' }
{ } { }
{ } { }
} {}{}{} mit
\mavergleichskette
{\vergleichskette
{ \delta }
{ > }{ 0 }
{ }{ }
{ }{ }
{ }{ }
} {}{}{,}
\mavergleichskette
{\vergleichskette
{ \delta }
{ \leq }{ \epsilon/M }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} und
\mavergleichskette
{\vergleichskette
{ \delta }
{ \leq }{ 1/(2L) }
{ }{ }
{ }{ }
{ }{ }
} {}{}{.}}
{} \teilbeweis {}{}{}
{Wir betrachten nun die Menge der \definitionsverweis {stetigen Abbildungen}{}{}
\mavergleichskettealignhandlinks
{\vergleichskettealignhandlinks
{C }
{ =} { { \left\{ \psi:J \rightarrow V \mid \psi \text{ stetig} , \, \Vert {\psi(t)- w} \Vert \leq \epsilon \text { für alle } t \in J \right\} } }
{ =} { { \left\{ \psi:J \rightarrow V \mid \psi \text{ stetig} , \, \Vert {\psi- w} \Vert \leq \epsilon \right\} } }
{ } { }
{ } { }
} {} {}{.} Dabei wird also $C$ mit der \definitionsverweis {Maximumsnorm}{}{} auf $J$ versehen. Dieser Raum ist nach Satz 55.9 und nach Aufgabe 36.22 wieder ein vollständiger metrischer Raum.}
{} \teilbeweis {}{}{}
{Wir betrachten nun auf diesem konstruierten Intervall $J$ bzw. der zugehörigen Menge $C$ die Abbildung \maabbeledisp {H} {C} {C } {\psi} {H(\psi) = (t \mapsto w + \int_{ t_0 }^{ t } f(s, \psi (s)) \, d s ) } {.} Dazu müssen wir zunächst zeigen, dass
\mathl{H(\psi)}{} wieder zu $C$ gehört. Für
\mavergleichskette
{\vergleichskette
{ t }
{ \in }{ J }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} ist aber nach Satz 39.1
\mavergleichskettealign
{\vergleichskettealign
{ \Vert { H( \psi)(t) - w} \Vert }
{ =} { \Vert { \int_{ t_0 }^{ t } f(s,\psi(s)) \, d s } \Vert }
{ \leq} { {{op:Betrag \int_{ t_0 }^{ t } \Vert {f(s,\psi(s))} \Vert \, d s }} }
{ \leq} { \betrag { t-t_0 } M }
{ \leq} { \frac{ \epsilon}{M} M }
} {
\vergleichskettefortsetzungalign
{ =} { \epsilon }
{ } {}
{ } {}
{ } {}
} {}{,} und
\mathl{H(\psi)}{} ist stetig, da es durch ein Integral definiert wird.}
{} \teilbeweis {}{}{}
{Zum Nachweis der Kontraktionseigenschaft seien
\mavergleichskette
{\vergleichskette
{ \psi_1, \psi_2 }
{ \in }{ C }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} gegeben. Für ein
\mavergleichskette
{\vergleichskette
{ t }
{ \in }{ J }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} ist
\mavergleichskettealignhandlinks
{\vergleichskettealignhandlinks
{ \Vert { H(\psi_1)(t)- H(\psi_2)(t) } \Vert }
{ =} { \Vert { \int_{ t_0 }^{ t } f(s,\psi_1(s)) \, d s - \int_{ t_0 }^{ t } f(s,\psi_2(s)) \, d s } \Vert }
{ =} { \Vert { \int_{ t_0 }^{ t } (f(s,\psi_1(s)) - f(s,\psi_2(s))) \, d s } \Vert }
{ \leq} { \betrag { \int_{ t_0 }^{ t } \Vert { f(s,\psi_1(s)) - f(s,\psi_2(s))} \Vert \, d s } }
{ \leq} { \betrag { \int_{ t_0 }^{ t } L \Vert {\psi_1(s) - \psi_2(s)} \Vert \, d s } }
} {
\vergleichskettefortsetzungalign
{ =} { L \cdot \betrag { \int_{ t_0 }^{ t } \Vert {\psi_1(s) - \psi_2(s)} \Vert \, d s } }
{ \leq} { L \cdot \betrag { \int_{ t_0 }^{ t } \Vert {\psi_1 - \psi_2} \Vert \, d s } }
{ \leq} { L \betrag { t-t_0 } \cdot \Vert {\psi_1 - \psi_2} \Vert }
{ \leq} { \frac{1}{2} \Vert {\psi_1- \psi_2} \Vert }
} {}{.} Da dies für jedes
\mavergleichskette
{\vergleichskette
{ t }
{ \in }{ J }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} gilt, folgt aus dieser Abschätzung direkt
\mavergleichskettedisp
{\vergleichskette
{ \Vert { H(\psi_1)- H(\psi_2) } \Vert }
{ \leq} {\frac{1}{2} \Vert {\psi_1- \psi_2} \Vert }
{ } { }
{ } { }
{ } { }
} {}{}{,} d.h. es liegt eine \definitionsverweis {starke Kontraktion}{}{} vor. Daher besitzt $H$ ein eindeutiges Fixelement
\mavergleichskette
{\vergleichskette
{ \psi }
{ \in }{ C }
{ }{ }
{ }{ }
{ }{ }
} {}{}{,} und diese Abbildung löst die Differentialgleichung. Dies gilt dann erst recht auf jedem offenen Teilintervall von $J$.}
{} \teilbeweis {}{}{}
{Damit haben wir insbesondere bewiesen, dass es in $C$ nur eine Lösung geben kann, wir wollen aber generell auf dem Intervall $J$ Eindeutigkeit erhalten. Für eine Lösung \maabb {v} {J} {V } {} gilt aber wegen der Integralbeziehung wieder
\mavergleichskettedisp
{\vergleichskette
{ v(t) }
{ =} { w + \int_{t_0}^t f(s,v(s)) ds }
{ } { }
{ } { }
{ } { }
} {}{}{} und die gleichen Abschätzungen wie weiter oben zeigen, dass die Lösung zu $C$ gehören muss.}
{}

}







\zwischenueberschrift{Die Picard-Lindelöf-Iteration}

Der Beweis des Satzes von Picard-Lindelöf ist prinzipiell konstruktiv. Darauf beruht die \stichwort {Picard-Lindelöf-Iteration} {,} mit der man Lösungen approximieren kann. Die Güte der Approximationen wird dabei durch geeignete Normen auf Funktionenräumen gemessen, was wir nicht ausführen.






\inputbemerkung
{}
{

Es sei $V$ ein \definitionsverweis {endlichdimensionaler}{}{} \definitionsverweis {reeller Vektorraum}{}{,}
\mavergleichskette
{\vergleichskette
{I }
{ \subseteq }{\R }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} ein \definitionsverweis {reelles Intervall}{}{,}
\mavergleichskette
{\vergleichskette
{U }
{ \subseteq }{V }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} eine \definitionsverweis {offene Menge}{}{} und \maabbeledisp {F} {I\times U} {V } {(t,v)} {F(t,v) } {,} ein \definitionsverweis {Vektorfeld}{}{} auf $U$. Es sei
\mavergleichskette
{\vergleichskette
{ (t_0,w) }
{ \in }{ I \times U }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} eine Anfangsbedingung. Es sei vorausgesetzt, dass dieses Vektorfeld \definitionsverweis {stetig}{}{} sei und \definitionsverweis {lokal einer Lipschitz-Bedingung}{}{} genüge. In der \stichwort {Picard-Lindelöf-Iteration} {} definiert man iterativ eine \definitionsverweis {Folge}{}{} von Funktionen \maabbdisp {\varphi_n} {I} {V } {} durch
\mavergleichskette
{\vergleichskette
{ \varphi_0 }
{ = }{ w }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} \zusatzklammer {dies ist also die konstante Funktion mit dem Wert $w$} {} {} und durch
\mavergleichskettedisp
{\vergleichskette
{\varphi_{n+1} (t) }
{ =} { w + \int_{t_0}^t F(s, \varphi_n(s))\,ds }
{ } { }
{ } { }
{ } { }
} {}{}{.} Dann gibt es ein Teilintervall
\mavergleichskette
{\vergleichskette
{ {]a,b[} }
{ \subseteq }{ I }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} mit
\mavergleichskette
{\vergleichskette
{ t_0 }
{ \in }{ {]a,b[} }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} derart, dass für
\mavergleichskette
{\vergleichskette
{ t }
{ \in }{ {]a,b[} }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} die Folge
\mathl{\varphi_n(t)}{} gegen einen Punkt
\mathl{\varphi(t)}{} konvergiert, wobei gleichmäßige Konvergenz vorliegt. Die Grenzfunktion $\varphi$ ist dann eine Lösung des Anfangswertproblems
\mathdisp {v'=f(t,v) \text{ und } v(t_0)=w} { . }
Bei einer \definitionsverweis {linearen Differentialgleichung}{}{} mit stetigen Koeffizientenfunktionen konvergiert dieses Verfahren auf ganz $I$.

}






\inputbemerkung
{}
{

Zu einem ortsunabhängigen Vektorfeld \maabbeledisp {F} {I \times \R^n} { \R^n } {(t,v)} {F(t,v) = F(t) } {,} und der Anfangsbedingung
\mavergleichskette
{\vergleichskette
{ v(t_0) }
{ = }{ w }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} führt die erste Picard-Lindelöf-Iteration auf
\mavergleichskettedisp
{\vergleichskette
{ \varphi_{1} (t) }
{ =} { w + \int_{t_0}^t F(s, w) ds }
{ =} { w + \int_{t_0}^t F(s) ds }
{ =} { w + G(t) }
{ } { }
} {}{}{,} wobei
\mathl{G(t)}{} eine Stammkurve zu
\mathl{F(t)}{} mit
\mavergleichskette
{\vergleichskette
{ G(t_0) }
{ = }{ 0 }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} sei. Die erste Iteration liefert hier also direkt die Lösung. Die kontrahierende Abbildung im Beweis zu Satz 56.2 ist in dieser Situation konstant.

}

Wir wenden dieses approximative Verfahren auf eine Differentialgleichung mit getrennten Variablen an, für die wir die Lösung schon kennen \zusatzklammer {siehe Aufgabe 30.11} {} {.}




\inputbeispiel{}
{

Wir wenden die \definitionsverweis {Picard-Lindelöf-Iteration}{}{} auf die Differentialgleichung
\mavergleichskettedisp
{\vergleichskette
{ y' }
{ =} {F(t,y) }
{ =} {ty }
{ } { }
{ } { }
} {}{}{} mit der Anfangsbedingung
\mavergleichskettedisp
{\vergleichskette
{ y(0) }
{ =} {1 }
{ } { }
{ } { }
{ } { }
} {}{}{} an \zusatzklammer {die Lösung ist \mathlk{e^{ { \frac{ 1 }{ 2 } } t^2}}{}} {} {.} Daher ist
\mavergleichskette
{\vergleichskette
{ \varphi_0 }
{ = }{ 1 }
{ }{ }
{ }{ }
{ }{ }
} {}{}{.} Die erste Iteration liefert
\mavergleichskettedisp
{\vergleichskette
{ \varphi_1 (t) }
{ =} { 1 + \int_0^t s ds }
{ =} { 1 + { \frac{ 1 }{ 2 } } t^2 }
{ } { }
{ } { }
} {}{}{.} Die zweite Iteration liefert
\mavergleichskettealign
{\vergleichskettealign
{ \varphi_2 (t) }
{ =} { 1 + \int_0^t F { \left( s, \varphi_1(s) \right) } ds }
{ =} { 1 + \int_0^t F { \left( s, 1 + { \frac{ 1 }{ 2 } } s^2 \right) } ds }
{ =} { 1 + \int_0^t s+ { \frac{ 1 }{ 2 } } s^3 ds }
{ =} { 1 + { \frac{ 1 }{ 2 } } t^2 + { \frac{ 1 }{ 8 } } t^4 }
} {} {}{.} Die dritte Iteration liefert
\mavergleichskettealign
{\vergleichskettealign
{ \varphi_3 (t) }
{ =} { 1 + \int_0^t F { \left( s, \varphi_2(s) \right) } ds }
{ =} { 1 + \int_0^t F { \left( s,1 + { \frac{ 1 }{ 2 } } s^2 + { \frac{ 1 }{ 8 } } s^4 \right) } ds }
{ =} { 1 + \int_0^t s+ { \frac{ 1 }{ 2 } } s^3 + { \frac{ 1 }{ 8 } } s^5 ds }
{ =} { 1 + { \frac{ 1 }{ 2 } } t^2 + { \frac{ 1 }{ 8 } } t^4 + { \frac{ 1 }{ 48 } } t^6 }
} {} {}{.} Dabei stimmt die $i$-te Iteration mit der Taylor-Entwicklung der Ordnung $2i$ der Lösung überein.


}






\inputbemerkung
{}
{

Es sei
\mavergleichskettedisp
{\vergleichskette
{v' }
{ =} {Mv }
{ } { }
{ } { }
{ } { }
} {}{}{} eine \definitionsverweis {lineare Differentialgleichung mit konstanten Koeffizienten}{}{} auf dem $\R^n$ und es sei eine Anfangsbedingung
\mavergleichskette
{\vergleichskette
{v(0) }
{ = }{w }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} gegeben. Wir behaupten, dass die $n$-te \definitionsverweis {Picard-Lindelöf-Iteration}{}{} gleich
\mavergleichskettedisp
{\vergleichskette
{ \varphi_n(t) }
{ =} { \sum_{k = 0}^n { \frac{ 1 }{ k! } } M ^{\circ k} (w) t^k }
{ } { }
{ } { }
{ } { }
} {}{}{} ist, wobei
\mathl{M ^{\circ k}}{} die $k$-fache Potenz der Matrix bezeichnet. Diese Aussage zeigen wir durch Induktion nach $n$. Für
\mavergleichskette
{\vergleichskette
{n }
{ = }{0 }
{ }{ }
{ }{ }
{ }{ }
} {}{}{} steht rechts einfach die konstante Kurve $w$. Es sei die Aussage nun für $n$ schon bewiesen. Dann ist
\mavergleichskettealign
{\vergleichskettealign
{\varphi_{n+1} (t) }
{ =} { w + \int_0^t M( \varphi_n(s) ) ds }
{ =} { w+ \int_0^t M { \left( \sum_{k = 0}^n { \frac{ 1 }{ k! } } M ^{\circ k} (w) s^k \right) } ds }
{ =} { w + \sum_{k = 0}^n { \frac{ 1 }{ k! } } \int_0^t s^k M { \left( M ^{\circ k} (w) \right) } ds }
{ =} {w + \sum_{k = 0}^n { \frac{ 1 }{ k! } } \int_0^t s^k M ^{\circ k+1} (w) ds }
} {
\vergleichskettefortsetzungalign
{ =} {w + \sum_{k = 0}^n { \frac{ 1 }{ k! } } { \frac{ 1 }{ (k+1) } } t^{k+1} M ^{\circ k+1} (w) }
{ =} {w + \sum_{k = 0}^n { \frac{ 1 }{ (k+1)! } } t^{k+1} M ^{\circ k+1} (w) }
{ =} {w + \sum_{\ell = 1}^{n+1} { \frac{ 1 }{ \ell ! } } t^\ell M ^{\circ \ell } (w) }
{ =} { \sum_{\ell = 0}^{n+1} { \frac{ 1 }{ \ell ! } } t^\ell M ^{\circ \ell } (w) }
} {}{} und die Aussage ist auch für
\mathl{n+1}{} richtig. Diese Approximationen sind die Anfangsglieder in der \anfuehrung{Exponentialreihe in dem Ausdruck}{} $t M$. Man kann zeigen, dass diese Exponentialreihe auf $\R$ konvergiert und in der Tat die Lösung des Anfangswertproblems ist \zusatzklammer {der Satz von Picard-Lindelöf sichert nur die Konvergenz auf einer Intervallumgebung} {} {.}

}




\inputbeispiel{}
{

Wir wenden die \definitionsverweis {Picard-Lindelöf-Iteration}{}{} auf das Anfangswertproblem
\mavergleichskettedisp
{\vergleichskette
{v(0) }
{ =} {(0,0) }
{ } { }
{ } { }
{ } { }
} {}{}{} zum Vektorfeld \maabbeledisp {F} {\R \times \R^2} { \R^2 } {(t,x,y)} { (t^2y-xy^2, x^3-y+t) } {} an. Es ist
\mavergleichskettedisp
{\vergleichskette
{\varphi_0 (t) }
{ =} { (0,0) }
{ } { }
{ } { }
{ } { }
} {}{}{.} Daher ist
\mavergleichskettedisp
{\vergleichskette
{\varphi_1(t) }
{ =} { \int_0^t (0,s) ds }
{ =} { \left( 0 , \, { \frac{ 1 }{ 2 } } t^2 \right) }
{ } { }
{ } { }
} {}{}{.} Es ist
\mavergleichskettedisp
{\vergleichskette
{F(s, \varphi_1(s)) }
{ =} { \left( s^2 { \frac{ 1 }{ 2 } } s^2 , \, - { \frac{ 1 }{ 2 } } s^2 +s \right) }
{ =} { \left( { \frac{ 1 }{ 2 } } s^4 , \, - { \frac{ 1 }{ 2 } } s^2 +s \right) }
{ } { }
{ } { }
} {}{}{} und daher
\mavergleichskettedisp
{\vergleichskette
{ \varphi_2(t) }
{ =} { \int_0^t \left( { \frac{ 1 }{ 2 } } s^4 , \, - { \frac{ 1 }{ 2 } } s^2 +s \right) ds }
{ =} { \left( { \frac{ 1 }{ 10 } } t^5 , \, - { \frac{ 1 }{ 6 } } t^3 + { \frac{ 1 }{ 2 } } t^2 \right) }
{ } { }
{ } { }
} {}{}{.} Wegen
\mavergleichskettealignhandlinks
{\vergleichskettealignhandlinks
{F(s, \varphi_2 (s)) }
{ =} { \left( s^2 { \left( - { \frac{ 1 }{ 6 } } s^3 + { \frac{ 1 }{ 2 } } s^2 \right) } - { \frac{ 1 }{ 10 } } s^5 { \left( - { \frac{ 1 }{ 6 } } s^3 + { \frac{ 1 }{ 2 } } s^2 \right) }^2 , \, { \left( { \frac{ 1 }{ 10 } } t^5 \right) }^3 - { \left( - { \frac{ 1 }{ 6 } } s^3 + { \frac{ 1 }{ 2 } } s^2 \right) } +s \right) }
{ =} { \left( - { \frac{ 1 }{ 6 } } s^5 + { \frac{ 1 }{ 2 } } s^4 - { \frac{ 1 }{ 10 } } s^5 { \left( { \frac{ 1 }{ 36 } } s^6 - { \frac{ 1 }{ 6 } } s^5 + { \frac{ 1 }{ 4 } } s^4 \right) } , \, { \frac{ 1 }{ 1000 } } s^{15} + { \frac{ 1 }{ 6 } } s^3 - { \frac{ 1 }{ 2 } } s^2 + s \right) }
{ =} { \left( - { \frac{ 1 }{ 360 } } s^{11} + { \frac{ 1 }{ 60 } } s^{10} - { \frac{ 1 }{ 40 } } s^{9} - { \frac{ 1 }{ 6 } } s^5 + { \frac{ 1 }{ 2 } } s^4 , \, { \frac{ 1 }{ 1000 } } s^{15} + { \frac{ 1 }{ 6 } } s^3 - { \frac{ 1 }{ 2 } } s^2 + s \right) }
{ } { }
} {} {}{} und daher
\mavergleichskettealignhandlinks
{\vergleichskettealignhandlinks
{ \varphi_3(t) }
{ =} { \int_0^t \left( - { \frac{ 1 }{ 360 } } s^{11} + { \frac{ 1 }{ 60 } } s^{10} - { \frac{ 1 }{ 40 } } s^{9} - { \frac{ 1 }{ 6 } } s^5 + { \frac{ 1 }{ 2 } } s^4 , \, { \frac{ 1 }{ 1000 } } s^{15} + { \frac{ 1 }{ 6 } } s^3 - { \frac{ 1 }{ 2 } } s^2 + s \right) ds }
{ =} { \left( - { \frac{ 1 }{ 4320 } } t^{12} + { \frac{ 1 }{ 660 } } t^{11} - { \frac{ 1 }{ 400 } } t^{10} - { \frac{ 1 }{ 36 } } t^6 + { \frac{ 1 }{ 10 } } t^5 , \, { \frac{ 1 }{ 16000 } } t^{16} + { \frac{ 1 }{ 24 } } t^4 - { \frac{ 1 }{ 6 } } t^3 + { \frac{ 1 }{ 2 } } t^2 \right) }
{ } { }
{ } { }
} {} {}{.}


}