Si la gente no cree que las matemáticas son simples, es solo
porque no se dan cuenta de lo complicado que es la vida.
– John Louis von Neumann

Introducción

¡Hemos llegado al final de la primera unidad de este curso!.

Concluiremos con la demostración de uno de los teoremas más importantes dentro del campo de las ecuaciones diferenciales; el teorema de existencia y unicidad de Picard – Lindelöf. Pero antes, un poco de contexto histórico.

Este resultado fue estudiado y desarrollado entre los años 1820 y 1900 por Cauchy, Liouville, Lipschitz, Picard y Lindelöf. Entre 1820 y 1830 Cauchy probó que si $f=f(x,y)$ es una función continua y existe su derivada parcial ${\displaystyle \frac{df}{dy}}$ continua en cierta región $U\subset {\mathbb{R}}^2$ que contiene el punto $(x_0,y_0)$, entonces existe un intervalo $\delta$ en el que un problema de valor inicial posee una única solución definida en $\delta$.

En 1838, Liouville simplificó la prueba de Cauchy introduciendo el método de las aproximaciones sucesivas, que más tarde continuarían siendo desarrolladas por Picard y que se conocerían como iterantes de Picard.

En 1876, Lipschitz mejoraría el resultado de Cauchy, sustituyendo la condición de que exista la derivada continua de $f$ por una menos fuerte, conocida como condición de Lipschitz.

Posteriormente, todo lo anterior fue ligeramente mejorado y generalizado por Picard (1890) y Lindelöf (1893), siguiendo las mismas ideas dadas por Liouville y Lipschitz.

Actualmente el método y los resultados se les atribuyen a Picard conociéndose como método de las iterantes de Picard y teorema de Picard (o más generalmente, teorema de Picard – Lindelöf).

En las dos últimas entradas hemos presentado una teoría preliminar con todas las herramientas necesarias para demostrar el teorema de Picard – Lindelöf, sin más, demostremos el teorema.

Teorema de existencia y unicidad de Picard – Lindelöf

El resultado global del teorema de existencia y unicidad de Picard – Lindelöf es el siguiente.

Demostración del teorema de Picard – Lindelöf

Sea $\delta =[a,b]$. Como cualquier función $y:\delta \to \mathbb{R}$ tiene su gráfica en $U$ y por hipótesis $f$ es continua en $U$, tenemos, como consecuencia del teorema sobre la ecuación integral, que $y:\delta \to \mathbb{R}$ es solución del PVI si y solo si $y(x)$ es una función continua en $\delta$ y para cada $x\in \delta$ verifica la ecuación integral

$$\begin{array}{}\text{(2)}& y(x)=y_0+{\int }_{x_0}^{x}f(t,y(t))dt\end{array}$$

Necesitamos probar que esta ecuación integral sólo posee una solución continua. Para ello, al ser $U$ una banda vertical, podemos definir sin problema alguno las iterantes de Picard $y_n:\delta \to \mathbb{R}$, las cuales son funciones continuas que verifican $y_n(x_0)=y_0$.

La demostración la dividiremos en tres secciones:

• Primero probaremos que la sucesión de iterantes $\{y_n\}$ converge uniformemente en el intervalo $\delta$ hacia una función continua $y:\delta \to \mathbb{R}$.
• Posteriormente comprobaremos que esta función $y:\delta \to \mathbb{R}$ verifica la ecuación integral ($\text{2}$) y, por tanto, es solución del PVI.
• Finalmente probaremos que el PVI no posee otra solución distinta de $y:\delta \to \mathbb{R}$.

Con los primeros dos puntos estaremos demostrando la existencia de una solución al problema de valor inicial y con el tercer punto estaremos demostrando la unicidad. Es importante mencionar que en cada uno de los tres puntos anteriores, aparte de la continuidad de $f$, haremos uso de manera esencial de la condición de Lipschitz de $f$ respecto de la segunda variable.

$$\begin{array}{}\text{(3)}& |f(x,y_1)-f(x,y_2)|\le L|y_1-y_2|\end{array}$$

para cada par de puntos $(x,y_1),(x,y_2)\in U$. Con $L$ la constante de Lipschitz para $f$ en $U$.

Así mismo, en el primer punto utilizaremos de forma esencial que la convergencia de las iterantes es uniforme, pues no basta con la convergencia puntual.

Demostremos el primer punto.

• Convergencia uniforme de las iterantes de Picard.

Para probar que la sucesión de iterantes $\{y_n\}$ converge uniformemente en el intervalo $\delta$ es conveniente expresarlas de la siguiente forma.

$$\begin{array}{}\text{(4)}& y_n(x)=y_0(x)+\sum _{m=1}^n(y_m(x)-y_{m-1}(x))\end{array}$$

Desglosa la serie anterior para que verifiques la equivalencia.

Fijado $x\in \delta$ es evidente que la sucesión numérica $\{y_n(x)\}$ es convergente en $\mathbb{R}$ si y sólo si la serie numérica $\sum _{m=1}^{\mathrm{\infty }}(y_m(x)-y_{m-1}(x))$ es convergente, para lo cual es suficiente con la convergencia absoluta de la serie para cada $x\in \delta$, es decir

$$\begin{array}{}\text{(5)}& \sum _{m=1}^{\mathrm{\infty }}|y_m(x)-y_{m-1}(x)|<\mathrm{\infty }\end{array}$$

Si la serie ($\text{5}$) fuese convergente para cada $x\in \delta$, entonces se tendría que la sucesión de iterantes converge puntualmente en $\delta$, sin embargo no es suficiente con la convergencia puntual; necesitamos algo más fuerte, como lo es la convergencia uniforme.

Para probar que la serie funcional $\sum _{m=1}^{\mathrm{\infty }}(y_m(x)-y_{m-1}(x))$ converge uniformemente en $\delta$ y, por tanto, la sucesión de iterantes, vamos a usar el criterio mayorante de Weierstrass para lo cual necesitamos probar que existen unas constantes $M_m\in {\mathbb{R}}^{+}$, tales que

$$\begin{array}{}\text{(6)}& |y_m(x)-y_{m-1}(x)|\le M_m\end{array}$$

para cada $x\in \delta$, cada $m=1,2,3\cdots$, y $\sum _{m=1}^{\mathrm{\infty }}{M}_m<\mathrm{\infty }$.

Vamos a comenzar con los casos $m=1$ y $m=2$, es decir, vamos a hallar las constantes $M_1$ y $M_2$, tales que

$$|y_1(x)-y_0(x)|\le M_{1}y|y_2(x)-y_1(x)|\le M_2$$

y con estos resultados intentaremos encontrar una relación de recurrencia para las constantes $M_m$ para luego corroborar que $\sum _{m=1}^{\mathrm{\infty }}{M}_m<\mathrm{\infty }$ y de esta manera probar la convergencia ($\text{5}$).

Partiendo de la ecuación de las iterantes de Picard

$$\begin{array}{}\text{(7)}& y_m(x)=y_0+{\int }_{x_0}^{x}f(t,y_{m-1}(t))dt;y_0=y_0(x)\end{array}$$

las primeras iterantes son

$$y_1(x)=y_0(x)+{\int }_{x_0}^{x}f(t,y_0(t))dty{y}_2(x)=y_0(x)+{\int }_{x_0}^{x}f(t,y_1(t))dt$$

de donde,

$$|y_1(x)-y_0(x)|=|{\int }_{x_0}^{x}f(t,y_0)dt|$$

y

$$|y_2(x)-y_1(x)|=|{\int }_{x_0}^{x}f(t,y_1(t))-f(t,y_0(t))dt|$$

Al momento de estimar $|y_1(x)-y_0(x)|$ necesitamos hacer la siguiente consideración. La función $f$ es continua en $U$ y, por tanto, la función

$$g:\delta \to \mathbb{R};x\to g(x)=f(x,y_0(x))$$

es continua en $\delta$. Como $\delta$ es compacto, la función $g(x)$ está acotada en $\delta$, es decir, existe una constante $H>0$, tal que

$$|g(x)|=|f(x,y_0(x))|\le H$$

para cada $x\in \delta$ y, por tanto, se verifica lo siguiente.

$$\begin{array}{rl}|y_1(x)-y_0(x)|& =|{\int }_{x_0}^{x}f(t,y_0(t))dt|\\ & \le {\int }_{x_0}^x|f(t,y_0(t))|dt\\ & \le {\int }_{x_0}^{x}Hdt\\ & =H|x-x_0|\end{array}$$

esto es,

$$\begin{array}{}\text{(8)}& |y_1(x)-y_0(x)|\le H|x-x_0|\end{array}$$

Si consideramos todo el intervalo $\delta =[a,b]$ podríamos obtener finalmente la estimación

$$\begin{array}{}\text{(9)}& |y_1(x)-y_0(x)|\le H(b-a)=M_1\end{array}$$

Para poder estimar adecuadamente $|y_2(x)-y_1(x)|$ consideremos el resultado ($\text{8}$), además de las siguientes dos desigualdades.

$$\begin{array}{rl}|y_2(x)-y_1(x)|& =|{\int }_{x_0}^{x}f(t,y_1(t))-f(t,y_0(t))dt|\\ \text{(10)}& & \le {\int }_{x_0}^x|f(t,y_1(t))-f(t,y_0(t))|dt\end{array}$$

y la condición de Lipschitz

$$\begin{array}{}\text{(11)}& |f(t,y_1(t))-f(t,y_0)|\le L|y_1(t)-y_0(t)|\end{array}$$

Supongamos que $x>x_0$. Usando ($\text{10}$) y ($\text{11}$), además del resultado ($\text{8}$), se tiene

$$\begin{array}{rl}|y_2(x)-y_1(x)|& \le {\int }_{x_0}^x|f(t,y_1(t))-f(t,y_0(t))|dt\\ & \le {\int }_{x_0}^{x}L|y_1(t)-y_0(t)|dt\\ & \le LH{\int }_{x_0}^x|t-x_0|dt\\ & =LH{\int }_{x_0}^x(t-x_0)dt\\ & =LH{\displaystyle \frac{(x-x_0{)}^2}{2}}\end{array}$$

Por otro lado, para $x<x_0$, se tiene

$$\begin{array}{rl}|y_2(x)-y_1(x)|& \le {\int }_x^{x_0}|f(t,y_1(t))-f(t,y_0(t))|dt\\ & \le L{\int }_x^{x_0}|y_1(t)-y_0(t)|dt\\ & \le LH{\int }_x^{x_0}|t-x_0|dt\\ & =LH{\int }_x^{x_0}(x_0-t)dt\\ & =LH{\displaystyle \frac{(x_0-x{)}^2}{2}}\end{array}$$

De ambos resultados, podemos afirmar que para cada $x\in \delta$

$$\begin{array}{}\text{(12)}& |y_2(x)-y_1(x)|\le HL{\displaystyle \frac{|x-x_0{|}^2}{2!}}\end{array}$$

La desigualdad ($\text{8}$) la podemos escribir de forma similar a ($\text{12}$) de la siguiente forma.

$$|y_1(x)-y_0(x)|\le H{L}^0{\displaystyle \frac{|x-x_0{|}^1}{1!}}$$

De estas dos relaciones establecemos una relación de recurrencia que vamos a probar por inducción sobre $m$. Proponemos que para cada $m=1,2,3,\cdots$, y para cada $x\in \delta$, se cumple

$$\begin{array}{}\text{(13)}& |y_m(x)-y_{m-1}(x)|\le H{L}^{m-1}{\displaystyle \frac{|x-x_0{|}^m}{m!}}\end{array}$$

La desigualdad ha sido probada anteriormente para $m=1$ y $m=2$. Supongamos que es cierta para $m$ y vamos a probar que es válida para $m+1$ siguiendo el mismo razonamiento que en la obtención del caso $m=2$. Vamos a mostrar el caso $x>x_0$, pero la prueba es similar para el caso $x<x_0$.

Si $x>x_0$, se tiene

$$\begin{array}{rl}|y_{m+1}(x)-y_m(x)|& \le {\int }_{x_0}^x|f(t,y_m(t))-f(t,y_{m-1}(t))|dt\\ & \le L{\int }_{x_0}^x|y_m(t)-y_{m-1}(t)|dt\\ & \le {\displaystyle \frac{H{L}^m}{m!}}{\int }_{x_0}^x(t-x_0{)}^{m}dt\\ & =H{L}^m{\displaystyle \frac{(x-x_0{)}^{m+1}}{(m+1)!}}\end{array}$$

De forma similar, si $x<x_0$, se tiene

$$|y_{m+1}(x)-y_m(x)|\le H{L}^m{\displaystyle \frac{(x_0-x{)}^{m+1}}{(m+1)!}}$$

De ambos resultados concluimos que

$$|y_{m+1}(x)-y_m(x)|\le H{L}^m{\displaystyle \frac{|x-x_0{|}^{m+1}}{(m+1)!}}$$

Es Importante hacer énfasis que en este desarrollo ha sido fundamental que las iterantes $\{y_n\}$ tengan sus gráficas en una región $U$ donde $f$ es lipschitziana.

De lo obtenido anteriormente, y considerando el intervalo completo $\delta =[a,b]$, obtenemos finalmente la siguiente desigualdad.

$$\begin{array}{}\text{(14)}& |y_m(x)-y_{m-1}(x)|\le {\displaystyle \frac{H}{L}}{\displaystyle \frac{(L(b-a){)}^m}{m!}}=M_m\end{array}$$

para cada $x\in \delta$ y cada $m=1,2,3,\cdots$. Como el intervalo $\delta$ es acotado, entonces $M_m\in {\mathbb{R}}^{+}$ y sabemos que

$$\begin{array}{}\text{(15)}& \sum _{m=1}^{\mathrm{\infty }}{\displaystyle \frac{(L(b-a){)}^m}{m!}}=e^{L(b-a)}-1<\mathrm{\infty }\end{array}$$

En definitiva,

$$\sum _{m=1}^{\mathrm{\infty }}{M}_m<\mathrm{\infty }$$

es decir la serie es convergente. Con esto queda probada la condición ($\text{5}$) y debido a que la prueba se hizo utilizando el criterio mayorante de Weierstrass concluimos que se trata de una convergencia uniforme de las iterantes de Picard en el intervalo $\delta$ hacia una función $y:\delta \to \mathbb{R}$.

Es bien conocido que si una sucesión $y_n:\delta \to \mathbb{R}$, $n=1,2,3,\cdots$, de funciones continuas sobre $\delta$ que convergen uniformemente en $\delta$ hacia una función $y:\delta \to \mathbb{R}$, la función límite uniforme $y$ también es continua en $\delta$.

Queda así demostrado el primer punto de la prueba. Ahora verifiquemos que la función límite uniforme $y:\delta \to \mathbb{R}$ verifica la ecuación integral ($\text{2}$) siendo la solución al problema de valor inicial.

• La existencia de la solución.

Sea $y:\delta \to \mathbb{R}$ la función obtenida anteriormente como límite uniforme de las iterantes de Picard $\{y_n\}$. La convergencia uniforme de $\{y_n\}$ hacia $y(x)$ en el intervalo $\delta$ significa que dado cualquier $\hat{\epsilon }>0$ existe un natural $N=N(\hat{\epsilon })$, tal que para cada $n>N$ y cada $x\in \delta$

$$\begin{array}{}\text{(16)}& |y_n(x)-y(x)|<\hat{\epsilon }\end{array}$$

Sabemos que la convergencia uniforme implica la convergencia puntual (pero no al revés), de manera que para cada $x\in \delta$ se cumple que

$$\begin{array}{}\text{(17)}& \underset{n\to \mathrm{\infty }}{lim}{y}_n(x)=y(x)\end{array}$$

Fijemos un $x\in \delta$. De acuerdo a ($\text{17}$) y usando ($\text{7}$), se tiene

$$y(x)=\underset{n\to \mathrm{\infty }}{lim}{y}_{n+1}(x)=y_0+\underset{n\to \mathrm{\infty }}{lim}{\int }_{x_0}^{x}f(t,y_n(t))dt$$

Por otro lado, sabemos que la función solución que satisface el PVI satisface también la ecuación integral ($\text{2}$),

$$y(x)=y_0+{\int }_{x_0}^{x}f(t,y(t))dt$$

Nuestro objetivo es probar que

$$\begin{array}{}\text{(18)}& \underset{n\to \mathrm{\infty }}{lim}{\int }_{x_0}^{x}f(t,y_n(t))dt={\int }_{x_0}^{x}f(t,y(t))dt\end{array}$$

Pues de esta forma la función límite uniforme $y$ verificaría la ecuación integral y, por tanto sería solución del PVI en el intervalo $\delta$, quedando así probada la existencia de la solución.

Demostrar la relación ($\text{18}$) es equivalente a probar que $\mathrm{\forall }\epsilon >0$ existe $N=N(\epsilon )\in \mathbb{N}$, tal que para cada $n>N$ y cada $x\in \delta$

$$\begin{array}{}\text{(19)}& |{\int }_{x_0}^{x}f(t,y_n(t))dt-{\int }_{x_0}^{x}f(t,y(t))dt|<\epsilon \end{array}$$

Para probar la relación ($\text{19}$) de nuevo haremos uso de la condición de Lipschitz ($\text{3}$) y de la convergencia uniforme de las iterantes hacia $y$ en $\delta$ ($\text{16}$).

$$\begin{array}{rl}|{\int }_{x_0}^{x}f(t,y_n(t))dt-{\int }_{x_0}^{x}f(t,y(t))dt|& \le {\int }_{x_0}^x|f(t,y_n(t))-f(t,y(t))|dt\\ & \le {\int }_a^b|f(t,y_n(t))-f(t,y(t))|dt\\ & \le L{\int }_a^b|y_n(t)-y(t)|dt\end{array}$$

Dado $\epsilon >0$, definimos

$$\begin{array}{}\text{(20)}& \hat{\epsilon }={\displaystyle \frac{\epsilon }{L(b-a)}}\end{array}$$

Con esto, la desigualdad ($\text{16}$) se puede escribir como

$$\begin{array}{}\text{(21)}& |y_n(t)-y(t)|<{\displaystyle \frac{\epsilon }{L(b-a)}}\end{array}$$

Usando esta desigualdad notamos que, para cada $n>N$ y cada $x\in \delta$

$$\begin{array}{rl}|{\int }_{x_0}^{x}f(t,y_n(t))dt-{\int }_{x_0}^{x}f(t,y(t))dt|& \le L{\int }_a^b|y_n(t)-y(t)|dt\\ & \le {\displaystyle \frac{L}{L(b-a)}}{\int }_a^b\epsilon dt\\ & ={\displaystyle \frac{L}{L(b-a)}}\epsilon (b-a)\\ & =\epsilon \end{array}$$

Por lo tanto, $\mathrm{\forall }\epsilon >0$ existe $N=N(\epsilon )\in \mathbb{N}$, tal que para cada $n>N$ y cada $x\in \delta$

$$|{\int }_{x_0}^{x}f(t,y_n(t))dt-{\int }_{x_0}^{x}f(t,y(t))dt|<\epsilon$$

lo que confirma la relación ($\text{18}$) que es lo que queríamos demostrar. De esta forma queda demostrada la existencia de la solución $y:\delta \to \mathbb{R}$ para el problema de valor inicial. Finalmente demostremos la unicidad de esta solución.

• Demostración de la unicidad.

Con los dos puntos anteriores estamos convencidos de la existencia de una solución $y:\delta \to \mathbb{R}$ que satisface el problema de valor inicial ($\text{1}$), así como la ecuación integral ($\text{2}$). La prueba de la unicidad se basa en la suposición de la existencia de otra solución $\hat{y}:\delta \to \mathbb{R}$ que igualmente cumple con los dos puntos anteriores y el objetivo será demostrar que $\hat{y}(x)=y(x)$.

De tarea moral demostrarás que la solución $\hat{y}(x)$ es también una función límite uniforme de las mismas iterantes de Picard para cada $x\in \delta$, esto es

$$\begin{array}{}\text{(22)}& \underset{n\to \mathrm{\infty }}{lim}{y}_n(x)=\hat{y}(x)\end{array}$$

o, lo que es equivalente, mostrar que

$$\begin{array}{}\text{(23)}& \underset{n\to \mathrm{\infty }}{lim}|\hat{y}(x)-y_n(x)|=0\end{array}$$

y por la ecuación ($\text{17}$) concluir que $\hat{y}(x)=y(x)$.

En esta situación se procede de manera muy similar a la prueba del primer punto en el que debemos encontrar una relación de recurrencia que acote a la cantidad $|\hat{y}(x)-y_n(x)|$ para cada $x\in \delta$ de la siguiente manera

$$\begin{array}{}\text{(24)}& 0\le |\hat{y}(x)-y_n(x)|\le B_n\end{array}$$

y si se prueba que

$$\underset{n\to \mathrm{\infty }}{lim}{B}_n=0$$

entonces quedará probada la relación ($\text{23}$).

A continuación te damos algunos hints y resultados que deberás obtener a lo largo de tu demostración.

Estudia lo que sucede con $n=1$ y $n=2$ y con los resultados encuentra la relación de recurrencia general para cada $n\in \mathbb{N}$, para ello considera la máxima distancia entre $\hat{y}$ y $y_0$, esto es

$$A=\underset{x\in \delta }{max}|\hat{y}-y_0|$$

El máximo $A\in {\mathbb{R}}^{+}$ esta asegurado gracias a la continuidad de la función $\hat{y}$ en el intervalo compacto $\delta$. Como la gráfica de la función $\hat{y}(x)$ está contenida en $U$ y $f=f(x,\hat{y})$ es una función lipschitziana, demuestra que para cada $x\in \delta$

$$\begin{array}{}\text{(25)}& |\hat{y}(x)-y_1(x)|\le L{\int }_{x_0}^x|\hat{y}(t)-y_0(t)|dt\le AL|x-x_0|\end{array}$$

Usando este resultado demuestra que

$$\begin{array}{}\text{(26)}& |\hat{y}(x)-y_2(x)|\le A{L}^2{\int }_{x_0}^x|t-x_0|dt\le A{L}^2{\displaystyle \frac{|x-x_0{|}^2}{2!}}\end{array}$$

Demuestra por inducción que en general, para cada $n=1,2,3,\cdots ,$ y $x\in \delta$

$$\begin{array}{}\text{(27)}& |\hat{y}(x)-y_n(x)|\le A{L}^n{\displaystyle \frac{|x-x_0|}{n!}}\end{array}$$

Este resultado te permite concluir que para cada $x\in \delta =[a,b]$

$$\begin{array}{}\text{(28)}& |\hat{y}(x)-y_n(x)|\le A{\displaystyle \frac{(L(b-a){)}^n}{n!}}=B_n\end{array}$$

Prueba que

$$\begin{array}{}\text{(29)}& \underset{n\to \mathrm{\infty }}{lim}{B}_n=\underset{n\to \mathrm{\infty }}{lim}{\displaystyle \frac{(L(b-a){)}^n}{n!}}=0\end{array}$$

Así finalmente queda demostrada la relación ($\text{23}$) y por lo tanto $\hat{y}(x)=y(x)$ para cada $x\in \delta$.

Realizar este ejercicio te servirá para consolidar mucho mejor lo que hemos realizado a lo largo de la demostración. Sin embargo, la demostración de la unicidad puede ser mucho más simple si aplicamos el lema de Gronwall. Demostremos la unicidad por esta opción.

Sean $y:\delta \to \mathbb{R}$ y $\hat{y}:\delta \to \mathbb{R}$ soluciones al PVI ($\text{1}$) y la ecuación integral ($\text{2}$).

$$y(x)=y_0+{\int }_{x_0}^{x}f(t,y(t))dty\hat{y}(x)=y_0+{\int }_{x_0}^{x}f(t,\hat{y}(t))dt$$

Restemos las dos ecuaciones anteriores y consideramos su valor absoluto.

$$\begin{array}{rl}|y(x)-\hat{y}(x)|& =|{\int }_{x_0}^{x}f(t,y(t))dt-{\int }_{x_0}^{x}f(t,\hat{y}(t))dt|\\ & =|{\int }_{x_0}^{x}f(t,y(t))-f(t,\hat{y}(t))dt|\\ & \le {\int }_{x_0}^x|f(t,y(t))-f(t,\hat{y}(t))|dt\end{array}$$

Como $f$ es una función lipschitziana respecto de la segunda variable, entonces

$${\int }_{x_0}^x|f(t,y(t))-f(t,\hat{y}(t))|dt\le {\int }_{x_0}^{x}L|y(t)-\hat{y}(t)|dt=L{\int }_{x_0}^x|y(t)-\hat{y}(t)|dt$$

es decir,

$$\begin{array}{}\text{(30)}& |y(x)-\hat{y}(x)|\le L{\int }_{x_0}^x|y(t)-\hat{y}(t)|dt\end{array}$$

Para que este resultado nos sea más familiar definamos lo siguiente.

$$h(x)=|y(x)-\hat{y}(x)|y\alpha =0,\beta =L$$

Usando esto reescribimos a la ecuación ($\text{30}$) como

$$\begin{array}{}\text{(31)}& 0\le h(x)\le \alpha +\beta {\int }_{x_0}^{x}h(t)dt\end{array}$$

Estamos en las condiciones del lema de Gronwall, pero en el caso especial en el que $\alpha =0$, así que aplicando el corolario del lema de Gronwall podemos concluir que para cada $x\in \delta$

$$h(x)=|y(x)-\hat{y}(x)|=0$$

lo que significa que $\mathrm{\forall }x\in \delta$, $y(x)=\hat{y}(x)$, es decir, la solución al problema de valor inicial es única.

Con esto quedan demostrados los tres puntos de la prueba, por lo tanto concluimos que el problema de valor inicial ($\text{1}$) posee una única solución en $\delta =[a,b]$ y además, las iterantes de Picard asociadas al PVI convergen uniformemente en $\delta$ hacia la solución $y:\delta \to \mathbb{R}$ del PVI.

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

¡Listo!. Hemos demostrado el teorema de existencia y unicidad de Picard – Lindelöf.

Apliquemos este resultado al caso de las ecuaciones diferenciales lineales.

Existencia y unicidad en ecuaciones lineales

Apliquemos el teorema de Picard – Lindelöf al caso de las ecuaciones diferenciales lineales de primer orden.

$$\begin{array}{}\text{(32)}& {\displaystyle \frac{dy}{dx}}=P(x)y+Q(x);y(x_0)=y_0\end{array}$$

Donde las funciones $P(x)$ y $Q(x)$ son continuas en un intervalo compacto $\delta =[a,b]$, $x_0\in \delta$ y $y_0\in \mathbb{R}$.

En este caso $U=\delta \times \mathbb{R}$ y

$$f(x,y)=P(x)y+Q(x)$$

Notemos que se verifica lo siguiente:

• $U$ es una banda vertical de base compacta, pues $\delta$ es un intervalo compacto.

• Como $P(x)$ y $Q(x)$ son continuas en $\delta$, entonces $f$ es continua en $U$.

• Como $P(x)$ es continua en el intervalo $\delta$ y éste es compacto, entonces la función $P(x)$ es acotada, así que podemos fijar $L>0$, tal que $|P(x)|<L$ para todo $x\in \delta$. Considerando esto tenemos que

$$\begin{array}{rl}|f(x,y_1)-f(x,y_2)|& =|(P(x)y_1+Q(x))-(P(x)y_2+Q(x))|\\ & =|P(x)y_1-P(x)y_2|\\ & =|P(x)||y_1-y_2|\\ & \le L|y_1-y_2|\end{array}$$

esto es

$$\begin{array}{}\text{(33)}& |f(x,y_1)-f(x,y_2)|\le L|y_1-y_2|\end{array}$$

para cada par de puntos $(x,y_1),(x,y_2)\in U$ es decir, $f$ es una función lipschitziana.

Por tanto, se cumplen las condiciones del teorema de existencia y unicidad global. En consecuencia ratificamos el resultado visto anteriormente en el que se asegura que cualquier problema de valor inicial asociado a una ecuación lineal posee solución única en el intervalo $\delta$. Además, ahora podemos afirmar que las iterantes de Picard asociadas convergen uniformemente hacia la solución del PVI.

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Un resultado importante que debemos revisar es que si dos problemas de valor inicial tienen valores iniciales muy cercanos entre sí, entonces las soluciones a cada PVI serán funciones muy próximas. A esto le llamamos dependencia continua de las soluciones respecto a condiciones iniciales. Revisemos este resultado. En la demostración será de uso esencial el lema de Gronwall.

Dependencia continua de la condición inicial

Demostración: Como $y(x)$ y $\hat{y}(x)$ son solución de sus respectivos PVI, entonces cada solución verifica una ecuación integral.

$$y(x)=y_0+{\int }_{x_0}^{x}f(t,y(t))dty\hat{y}(x)={\hat{y}}_0+{\int }_{x_0}^{x}f(t,\hat{y}(t))dt$$

Vemos que

$$\begin{array}{rl}|y(x)-\hat{y}(x)|& =|y_0-{\hat{y}}_0+{\int }_{x_0}^{x}f(t,y(t))dt-{\int }_{x_0}^{x}f(t,\hat{y}(t))dt|\\ & \le |y_0-{\hat{y}}_0|+|{\int }_{x_0}^{x}f(t,y(t))-f(t,\hat{y}(t))dt|\\ & \le |y_0-{\hat{y}}_0|+{\int }_{x_0}^x|f(t,y(t))-f(t,\hat{y}(t))|dt\end{array}$$

Sabemos que $f$ es una función lipschitziana respecto de la segunda variable en $R$ de manera que

$$\begin{array}{}\text{(35)}& |f(x,y(x))-f(x,\hat{y}(x))|\le L|y(x)-\hat{y}(x)|\end{array}$$

con $L$ la constante de Lipschitz para $f$ en $R$. Entonces,

$$\begin{array}{}\text{(36)}& |y(x)-\hat{y}(x)|\le |y_0-{\hat{y}}_0|+L{\int }_{x_0}^x|y(t)-\hat{y}(t)|dt\end{array}$$

Definamos

$$0<g(x)=|y(x)-\hat{y}(x)|,\alpha =|y_0-{\hat{y}}_0|y\beta =L$$

Con esto la desigualdad ($\text{36}$) la podemos reescribir como

$$0<g(x)\le \alpha +\beta {\int }_{x_0}^{x}g(t)dt$$

Ahora podemos aplicar el lema de Gronwall.

$$g(x)\le \alpha e^{\beta (x-x_0)}$$

es decir,

$$|y(x)-\hat{y}(x)|\le |y_0-{\hat{y}}_0|e^{L(x-x_0)}$$

Que es lo que queríamos demostrar.

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

De este resultado observamos que si

$$y_0\to {\hat{y}}_0\Rightarrow |y_0-{\hat{y}}_0|\to 0$$

Entonces las soluciones de los PVI serán funciones muy próximas

$$|y(x)-\hat{y}(x)|\to 0\Rightarrow y(x)\to \hat{y}(x)$$

Para concluir la entrada hagamos un breve comentario sobre el resultado local del teorema de Picard y realicemos unos ejemplos al respecto.

Teorema de existencia y unicidad local

Recordemos que el resultado local del teorema de existencia y unicidad de Picard – Lindelöf establece lo siguiente.

La demostración a este teorema corresponde a una adaptación de la demostración vista para el caso global, teniendo en cuenta que las gráficas de las iterantes de Picard, así como la de cualquier posible solución, definidas en el intervalo $\delta =[x_0-h,x_0+h]$, están dentro del rectángulo $R$ donde la función $f$ es continua y lipschitziana respecto de la segunda variable. Los pasos claves a seguir y las técnicas son prácticamente una repetición de lo visto anteriormente cambiando la banda vertical $U=[a,b]\times \mathbb{R}$ por el rectángulo

$$R=\{(x,y)\in \mathbb{R}\mid |x-x_0|\le a,|y-y_0|\le b,a,b\in \mathbb{R}\}$$

Para conocer sobre los detalles puedes revisar la demostración del teorema local en los videos de este mismo curso.

Finalmente, resolvamos algunos ejemplos.

Ejemplo: Mostrar que el problema de valor inicial

$${\displaystyle \frac{dy}{dx}}={\sin }^2(x-y);y(0)=0$$

posee una única solución definida en $\mathbb{R}$.

Solución: En este caso tenemos la función $f:U=\mathbb{R}\times \mathbb{R}\to \mathbb{R}$ definida como

$$f(x,y)={\sin }^2(x-y)$$

Es claro que $f$ es continua en la región $U=\mathbb{R}\times \mathbb{R}$. La función derivada parcial ${\displaystyle \frac{\partial f}{\partial y}}:U\to \mathbb{R}$ está dada como

$${\displaystyle \frac{\partial f}{\partial y}}=-2\sin (x-y)\cos (x-y)$$

Como

$$|\sin (x-y)|\le 1y|\cos (x-y)|\le 1$$

para todo $(x,y)\in U$, entonces

$$|{\displaystyle \frac{\partial f}{\partial y}}(x,y)|\le 2$$

para todo $(x,y)\in U$. En consecuencia $f$ es una función lipschitziana en $U$ respecto de la segunda variable. Con esto hemos probado que se satisfacen las hipótesis del teorema de existencia y unicidad global por lo que podemos asegurar que el PVI posee una única solución definida en $\mathbb{R}$.

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Calcular las iterantes no siempre será sencillo. En el ejemplo anterior las iterantes pueden no ser fácil de desarrollar, pero debido a que satisface el teorema de Picard – Lindelöf podemos asegurar que dichas iterantes van a converger a la solución del PVI.

Ejemplo: Mostrar que aplicando el resultado global del teorema de Picard – Lindelöf al problema de valor inicial

$${\displaystyle \frac{dy}{dx}}=y^2;y(0)=1$$

no es posible asegurar la existencia y unicidad de la solución.

Solución: La función $f:U=\mathbb{R}\times \mathbb{R}\to \mathbb{R}$ definida como

$$f(x,y)=y^2$$

es continua en $U=\mathbb{R}\times \mathbb{R}$, sin embargo su derivada parcial

$${\displaystyle \frac{\partial f}{\partial y}}=2y$$

no está acotada en $U$ por lo que $f$ no es una función lipschitziana en $U$.

Una observación más es que la solución al PVI dada por

$$y(x)={\displaystyle \frac{1}{1-x}}$$

no está definida en $\mathbb{R}$ si no en el intervalo $(-\mathrm{\infty },1)$.

En definitiva, como no se cumple la tercera condición del teorema global, entonces no podemos asegurar nada sobre la existencia y unicidad de la solución del PVI.

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Veamos ahora la importancia del resultado local del teorema de existencia y unicidad de Picard – Lindelöf . Resolvamos de nuevo el ejemplo anterior, pero ahora considerando una región $R$ alrededor del punto dado por la condición inicial.

Ejemplo: Mostrar que aplicando el resultado local del teorema de Picard – Lindelöf, el problema de valor inicial

$${\displaystyle \frac{dy}{dx}}=y^2;y(0)=1$$

posee una única solución. Encontrar el intervalo de existencia y unicidad.

Solución: Es claro que la función

$$f(x,y)=y^2$$

es continua en ${\mathbb{R}}^2$ por lo que $f$ será una función lipschitziana en cualquier conjunto $R$ convexo y compacto. Consideremos el rectángulo centrado en el valor inicial $(0,1)$ de dimensiones $a=2$ y $b=1$, es decir

$$R=\{(x,y)\in \mathbb{R}\mid |x|\le 2,|y-1|\le 1\}=[-2,2]\times [0,2]$$

En la región $R$ la función $f$ si es lipschitziana y continua por lo que se satisfacen las condiciones del teorema local de existencia y unicidad de Picard. Este teorema nos dice que existe una única solución definida en el intervalo $\delta =[x_0-h,x_0+h]$ donde

$$h=min\{a,{\displaystyle \frac{b}{M}}\}yM\ge \underset{(x,y)\in R}{max}|f(x,y)|$$

En este caso, como el máximo valor que puede tomar $y$ en el rectángulo $R$ es $y=2$, entonces

$$M=\underset{(x,y)\in R}{max}|f(x,y)|=\underset{(x,y)\in R}{max}|y^2|=4$$

Usando este resultado, se tiene

$$h=min\{a,{\displaystyle \frac{b}{M}}\}=min\{2,{\displaystyle \frac{1}{4}}\}={\displaystyle \frac{1}{4}}$$

Por lo tanto, podemos asegurar la existencia y unicidad de la solución del PVI en el intervalo

$$\delta =[x_0-h,x_0+h]=[-{\displaystyle \frac{1}{4}},{\displaystyle \frac{1}{4}}]$$

Además podemos asegurar que las iterantes de Picard convergen uniformemente en el intervalo $\delta$ hacia la solución única del PVI. A saber, convergen a

$$y(x)={\displaystyle \frac{1}{1-x}}$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Con esto concluimos la primera unidad del curso.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

1. Completar la demostración de la unicidad de la solución a un problema de valor inicial que cumple con las hipótesis del teorema global de existencia y unicidad. Recuerda que el objetivo es demostrar que $$\underset{n\to \mathrm{\infty }}{lim}|\hat{y}(x)-y_n(x)|=0$$

2. Comprobar que el problema de valor inicial $${\displaystyle \frac{dy}{dx}}={\displaystyle \frac{y}{x}};y(0)=0$$ posee infinitas soluciones en cualquier intervalo $\delta$ en el que $0\in \delta$.
   ¿Porqué no contradice esto al teorema de existencia y unicidad local?.

3. Determinar, por el método de iterantes de Picard, la solución del siguiente problema de valor inicial: $${\displaystyle \frac{dy}{dx}}=2y(1+x);y(-1)=1$$

4. Comprobar que el mayor intervalo que proporciona el teorema local de existencia y unicidad de Picard para el problema de valor inicial $${\displaystyle \frac{dy}{dx}}=1+x^2;y(0)=0$$ donde se asegura la existencia y unicidad de la solución, es $\delta =[-{\displaystyle \frac{1}{2}},{\displaystyle \frac{1}{2}}]$.

Más adelante…

Con la demostración del teorema de existencia y unicidad de Picard – Lindelöf justificamos la teoría realizada a lo largo de esta primera unidad.

En la siguiente entrada comenzaremos la unidad 2 del curso. En dicha unidad estudiaremos las ecuaciones diferenciales de orden mayor a uno, en particular estudiaremos con mayor detalle las ecuaciones diferenciales de segundo orden.

Entradas relacionadas

• Página principal del curso: Ecuaciones Diferenciales I
• Entrada anterior del curso: Teorema de Existencia y Unicidad – Iterantes de Picard y Convergencia
• Siguiente entrada del curso: Ecuaciones diferenciales de orden superior

• Video relacionado al tema: Teorema de Existencia y Unicidad de Picard

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»