# Picard's existence theorem, successive approximations and the global solution

Picard’s existence theorem states that if $U$ is an open subset of $\mathbb{R}^2$ and $f$ is a continuous function on $U$ that is Lipschitz continuous with respect to the second variable then there is some $h>0$ such that the initial value problem

$$y’=f(x,y), \\y(a)=b$$

has a unique solution on an interval $(a-h,a+h)$ in $U$.

Now assume that I wanted to solve $$y’=y, \\ y(0)=1$$ using the fixed point theorem. Solving the above problem is equivalent to solving
$$\phi'(x)=\int_0^x\phi(t)dt+1$$
for $\phi$.

Now consider the map $F:C([-r,r])\rightarrow C([-r,r])$ given by

$$(F(\psi))(x)=\int_0^x\psi(t)dt+1$$
where $r>0$.

Since $C([-r,r])$ is complete, if we only proved that $F$ is a contraction then we could use the fixed point theorem to find the solution. But $F$ is a contraction only if $2r<1$. So choose $r<\frac{1}{2}$.

Using successive approximations we find out that $y=e^x$ is the solution on an interval $(-\frac{1}{2},\frac{1}{2})$. But obviously this is not the maximal open interval on which $y=e^x$ is our solution. I wonder how we can proceed from here to conclude that $y=e^x$ is the solution on whole $\mathbb{R}$.

#### Solutions Collecting From Web of "Picard's existence theorem, successive approximations and the global solution"

The global solution of an initial value problem (IVP)
$$x’=f(t,x),\quad x(t_0)=x_0\tag{1}$$
is not obtained via Picard’s theorem, but through analytic continuation. The grand picture is as follows:

We are given an open domain $\Omega\subset{\mathbb R}\times{\mathbb R}^n$ and a continuous function $$f:\quad \Omega\to{\mathbb R}^n,\qquad(t,x)\mapsto f(t,x)$$
($n=1$ in the sequel) which is locally Lipschitz continuous with respect to the second variable. This means that for every point $(t_0,x_0)\in\Omega$ there is a rectangular window $W$ with center $(t_0,x_0)$ and a constant $C$ such that $$|f(t,x)-f(t,x’)|\leq C|x-x’|\qquad\bigl((t,x),(t,x’)\in W\bigr)\ .$$
An $f\in C^1(\Omega)$ automatically fulfills these conditions. Picard’s theorem then guarantees that for every $(t_0,x_0)$ there is a window $$[t_0-h,t_0+h]\times [x_0-q,x_0+q]\subset W$$ and a micro-solution $$\phi_0:\quad[t_0-h,t_0+h]\to [x_0-q,x_0+q] \tag{2}$$ of $(1)$.

The micro-solution $(2)$ is continuous up to $t_1:=t_0+h$, and even satisfies $\phi_0′(t_1-)=f\bigl(t_1,\phi_0(t_1)\bigr)$. The uniqueness part of Picard’s theorem then allows to conclude that immediately to the left of $t_1$ this $\phi_0$ coincides with the micro-solution $\phi_1$ of the IVP
$$x’=f(t,x),\quad x(t_1)=\phi_0(t_1)\ .$$
Note that $\phi_1$ is defined up to the point $t_2:=t_1+h’$ for some $h’>0$, hence extends $\phi_0$ further to the right. Concatenating $\phi_0$ and $\phi_1$, and proceeding in this way forever (and similarly to the left of $t_0$) we arrive at the global solution $$\phi_*:\quad I_* \to{\mathbb R}$$ of the IVP $(1)$. The interval $I_*$ on which $\phi_*$ is defined may be infinite, and in general depends on the initial point $(t_0,x_0)$. About this $\phi_*$ one can say the following:

(i) Global uniqueness: Any solution $\phi$ of $(1)$, defined in some open interval $I\subset{\mathbb R}$, is part of $\phi_*$. This means that $I\subset I_*$, and that $\phi(t)=\phi_*(t)$ for all $t\in I$.

(ii) The graph of $\phi_*$ extends to infinity, or to the boundary of $\Omega$. More precisely: Given any compact set $K\subset\Omega$ there are points $\bigl(t,\phi_*(t)\bigr)\notin K$.

These things are proven in detail in books titled “Theory of differential equations”.

First I prove a lemma :

If $F$ has one iterated which is a contraction, when $F$ has a unique fixed point.

Proof : If $F^n$ is a contraction then it has a unique fixed point $\psi$. The equalities $$F^n\circ F(\psi)=F\circ F^n(\psi)=F(\psi)$$ show that, by unicity of fixed point, $F(\psi)=\psi,$ and so existence of fixed point for $F.$ For unicity, we see that the existence a new fixed point $\varphi$ implies that $$F^n(\varphi)=F\circ\dots\circ F\circ F(\varphi)=F\circ\dots\circ F(\varphi)=\dots=\varphi$$ and so $\varphi=\psi.$ $\square$

Now let $r\in\mathbb{R}^+,$ and let’s show that one of the iterated of $F$ is a contraction on $C([0,r]).$ We have that :
$$|(F(\psi)-F(\varphi))(t)|=|\int_0^t\psi(s)-\varphi(s)|\mathrm{d}s\leq t\,||\psi-\varphi||_\infty$$ and so
$$||F(\psi)-F(\varphi)||_{\infty}\leq r\cdot||\psi-\varphi||_{\infty}$$ so $F$ is $r-$lipschitz. Now see that
$$|(F^2(\psi)-F^2(\varphi))(t)|\leq\int_0^t|F(\psi(s))-F(\varphi(s))|\mathrm{d}s\leq\int_0^t s||\psi-\varphi||_\infty\mathrm{d}s=\frac{t^2}{2}||\psi-\varphi||_{\infty}.$$

One can see that $$||F^n(\psi)-F^n(\varphi)||_\infty\leq\frac{r^n}{n!}||\psi-\varphi||_\infty$$ and as it exists $n\in\mathbb{N}$ such that $n!>r^n$ we establish our result.

For ODE defined on a stripe $(a,b)×\Bbb R^n$ with a global Lipschitz constant $L$ there, one can use the modified supremum norm
$$\|y\|_L=\sup e^{-2L·|t-t_0|}\|y(t)\|$$
where the initial condition is at $t_0\in (a,b)$. One can then show that in this norm the Picard iteration is contractive for continuous functions over the whole of $(a,b)$, giving the maximally possible interval at once.

Existence is usually related to compactness in some or the other way. In your case you’d try to estimate the derivate of the solution on the interval for which you have a solution, $(s,t)$, say, and try to show that the solution is bounded and has bounded derivates on $[s,t]$. You then continue to solve in a neighbourhood of these points and extend your solution. In case you can generalize this procedure independently of $s, t$ then you can conclude that the set on which the solution exists is closed. Picard’s theorem tells you it’s open, so it’s connected. Picard also gives you uniqueness, as a nice addon.

Now the only subset of $\mathbb{R}$ which is closed and open is $\mathbb{R}$ itself, so you are done (in this case. Of course proving compactness prove is often not possible).

$$\phi'(x) = \phi(x),\;\;\; \phi(0)=1.$$
You can see that $\phi'(x)=\phi(x),\; \phi(a)=C$ has a unique solution on any interval $[a,a+h]$ for sufficiently small $h$. Because of this, if you solve $\phi_1’=\phi_1,\;\;\phi(a_1)=C$ on $[a_1,b_1]$, then you can next solve $\phi_2’=\phi_2,\;\;\; \phi_2(a_2)=\phi_1(b_1)$, glue the two functions together on $[a_1,b_2]$ and conclude that
$$\phi(x) = \int_{a_1}^{x}\phi'(t)dt+C,\;\;\; a_1 \le x \le b_2.$$
If you let $\mathscr{C}=\{ y \ge a_1 : \mbox{integral equation solution extends to } [a,y] \}$, then you can show that $\mathscr{C}$ is both open and closed in $[a,\infty)$, which leads to the conclusion that $\mathscr{C}=[a_1,\infty)$. Showing the set is open uses the above extension argument; you can use a continuity argument to extend solutions on $[a,y)$ to $[a,y]$ in order show that $\mathscr{C}$ is closed. Arguing in this way leads to global existence and uniqueness, knowing only that you have a local Lipschitz condition.
Picard iteration has an added feature that can be characterized by showing that an operator $T^n$ is contractive for some $n \ge 1$. For example, if you iterate
$$\phi(x)=\int_{a}^{x}\phi(t_1)dt_1+A \\ = \int_{a}^{x}\int_{a}^{t_1}\phi(t_2)dt_2dt_1+A(x-a)+A,$$
then you can stop after $n$ iterations so long as the following integral operator is contractive:
$$L_n\phi = \int_{a}^{x}\int_{a}^{t_1}\int_{a}^{t_{2}}\cdots\int_{a}^{t_{n-1}}\phi(t_n)dt_n dt_{n-1}\cdots dt_1.$$
For $a \le x \le b$,
$$\|L_n\phi\| \le \frac{(b-a)^{n}}{n!}\|\phi\|$$
For fixed $b$ and large enough $n$, the norm coefficient $\frac{(b-a)^{n}}{n!}$ can be made arbitrarily small because it is the general term of a convergent exponential series.