UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS FACULTAD DE CIENCIAS MATEMÁTICAS E.A.P. DE MATEMÁTICA

Aplicaciones del teorema del punto fijo de Banach

TESIS para optar el Título Profesional de Licenciado en Matemática Pura

AUTOR Julio Román Loayza Cerrón

ASESOR Tomás Núñez Lay

Lima – Perú 2006

Aplicaciones del teorema del punto fijo de Banach

Julio Román Loayza Cerrón

Tesis presentada a consideración del Cuerpo Docente de la Facultad de Ciencias Matemáticas de la Universidad Nacional Mayor de San Marcos, como parte de los requisitos para obtener el título profesional de Licenciado en Matemática Pura. Aprobada por:

________________________________ Dr. Víctor Rafael Cabanillas Zannini

________________________________ Dr. Alfonso Pérez Salvatierra

________________________________ Mg. Tomás Núñez Lay

LIMA-PERÚ Diciembre 2006 ii

i A la memoria de LUCIO LOAYZA GALINDO. Mi siempre querido y recordado PADRE.

i A JULIA mi siempre querida MADRE, por sus enseñanzas y constante apoyo.

i A Rosario, Rocío, Aristóteles, Mirko, Rosa mis hermanos por entender lo absorvente que es mi profesión.

i Lucio y Gabriel Alonso la nueva generación de mi familia.

iii

AGRADECIMIENTOS

i Al Prof. Tomás Núñez Lay por asesorarme en la elección, elaboración y desarrollo del tema de mi tesis. Con su apoyo logré participar como expositor en tres eventos, dos nacionales y uno internacional; teniendo la gran responsabilidad, así como gran orgullo de representar a mi universidad y a nuestro país, trascendiendo en mi vida profesional y personal. Una vez más agradezco al Prof. Tomás Núñez Lay por el invalorable apoyo académico que me brindó, pero sobre todo por la enorme paciencia que me tuvo.

i Al Prof. Victor Cabanillas Zannini, así como al Prof. Alfonso Pérez Salvatierra, por sus importantes observaciones luego de revisar mi tesis, contribuyendo así poder hacer algunas importantes precisiones convenientes y en otros casos desarrollar de mejorar manera los capítulos de esta tesis; y en ese camino estar mas conciente de lo que escribía.

i Al Prof. Eugenio Cabanillas Lapa quien de manera muy atenta y cordial aceptó realizar una primera revisión del primer borrador que imprimí de mi tesis; dándome las primeras pautas a seguir. Además por su constante apoyo moral en la culminación de este trabajo.

i A la Prof. Roxana López Cruz por su apoyo como flamante directora de nuestra escuela, en agilizar los documentos, para quedar expedito para la sustentación de mi tesis en esta fecha. Así como su ayuda con el idioma inglés, permitiéndome redactar correctamente el abstract de este trabajo; y apoyo moral permanente.

i Al magíster en matemática pura Walter Huaraca Vargas, reciente título obtenido en la Universidad de Sao Paulo filial Sao Carlos Brasil; amigo sanmarquino, por la ayuda, consejos y apoyo que me dió cuando estuve en esa casa de estudios el verano del 2004, en realidad fue él quien me hizo ver la necesidad de trabajar una tesis, y luego continuar de manera seria y responsable los estudios en matemática pura.

i A José Luis Quispe Castillo amigo sanmarquino, por brindarme algunos apuntes sobre el tema de mi tesis, y sus constantes ánimos a lo largo de la elaboración de mi paper.

i A Jorge Astocondor Félix amigo sanmarquino, por la ayuda que me brindó en elaborar los gráficos que posee este paper, además por su apoyo moral.

i A Orlando Luciano Julca Jesús, Jairo Yamil Esquivel Ortiz primero compañeros de aula en nuestra facultad, luego leales y sinceros amigos, por esa gran amistad , apoyo moral y ánimos para lograr terminar mi tesis.

iv

RESUMEN

APLICACIONES DEL TEOREMA DEL PUNTO FIJO DE BANACH

JULIO ROMÁN LOAYZA CERRÓN

DICIEMBRE 2006

Orientador:

Prof. Tomás Núñez Lay.

Título obtenido:

Licenciado en Matemática.

----------------------------------------------------------------------------------------------------------------Para aplicar el Teorema del Punto Fijo de Banach ( T.P.F.B.), se necesita una aplicación contractiva de un espacio completo en sí mismo; este resultado garantiza la existencia y unicidad de la solución de un problema específico. El teorema nos provee de un método iterativo, para construir la solución aproximada con cierto margen de error previamente fijado. Por lo mencionado, el T.P.F.B. ó método de las aproximaciones sucesivas (M.A.S.) se convierte en una potente herramienta del análisis, lo que quedará evidenciado luego de presentar algunas importantes aplicaciones del T.P.F.B., tales como: 1.- Métodos numéricos. 2.- Ecuaciones integrales. 3.- Ecuaciones diferenciales ordinarias. 4.- Ecuaciones en derivadas parciales. 5.- Problema de Sturm-Liouville. 6.- Programación dinámica. 7.- Dinámica compleja. PALABRAS CLAVES:

i Espacios métricos completos.

i Contracción.

i Puntos fijos.

i Existencia y unicidad. v

ABSTRACT

APPLICATIONS OF THE FIXED POINT BANACH’S THEOREM

JULIO ROMÁN LOAYZA CERRÓN DECEMBER 2006 Management by:

Tomás Núñez Lay.

Obtained Title:

Mathematics License.

----------------------------------------------------------------------------------------------------------------To apply the Fixed Point Banach’s Theorem (F.P.B.T.) , we need a contracting application mapping a complete metric space into itself. The hypothesis guarantees the existences and uniqueness of solution of a specific problem, whose must be planted as a problem to find fixed points. The theorem provides to us with a iterative method to construct the approximated solution with a certain margin of error previously fixed. According to before, the F.P.B.T. or successive approximations method (S.A.M.) becomes a powerful tool into the analysis, it can be demonstrated through the next important applications: 1.- Numerical methods. 2.- Integral equations. 3.- Ordinary differential equations. 4.- Partial differential equations. 5.- Sturm-Liouville Problem. 6.- Dynamic programming. 7.- Complex dynamic. KEY WORDS:

i Complete metric space.

i Contraction.

i Fixed point.

i Existence and uniqueness. vi

INDICE i INTRODUCCIÓN..……………………………………..…………………... 1 i CAPÍTULO I : DEFINICIONES Y TEOREMAS BÁSICOS.

I.1. Punto Fijo..………………………………………………………………. 2 I.2. Contracción….………………………………………………………....... 2 I.3. Iteración..……………………………………………………………........ 6 I.4. Teorema del punto fijo de Banach (T.P.F.B.)….…………….…………. 6 I.5. Interpretación geométrica del T.P.F.B…..……………………………… 23

i CAPÌTULO II : APLICACIONES DEL TEOREMA DEL PUNTO FIJO DE BANACH (T.P.F.B.).

II.1. Aplicación a los sistemas de ecuaciones lineales algebraicas……..…... 30 II.2. Aplicación a las ecuaciones integrales…...................………………….. 38 II.3. Aplicación a las ecuaciones diferenciales ordinarias…….……………. 48 II.4. Aplicación a las ecuaciones diferenciales parciales.....……..…………. 52 II.5. Aplicación al problema de Sturm-Liouville……..……………………...71 II.6. Aplicación a la programación dinámica………..………………........... 79 II.7. Aplicación a la dinámica compleja…..……………………………....... 85

i Bibliografía.………………………………………………………………… 99

INTRODUCCIÓN Presentaremos el Teorema del Punto Fijo de Banach (T.P.F.B.), que es el análogo abstracto del método de las aproximaciones sucesivas de Picard y de otros resultados semejantes. El T.P.F.B. se desarrolla para aplicaciones contractivas definidas en un espacio métrico completo en el mismo; por esto al T.P.F.B. también se le llama teorema de la aplicación contractante (T.A.C.). El T.P.F.B. es importante porque proporciona una técnica que permite la prueba de varios teoremas de existencia y unicidad en diferentes ramas del análisis, así como en otras áreas de la matemática. El teorema además, nos proporciona una forma explícita de construir la solución aproximada mediante un proceso iterativo, que converge hacia la solución exacta; incluso mostraremos un resultado fácilmente deducible, el cual nos provee una fórmula para saber el número de iteraciones suficientes para obtener una aproximación deseada con cierto margen de error previamente establecido. La presente monografía está divida en dos capítulos, en el primero abordamos todo lo concerniente al aspecto teórico dotándolo de ejemplos, los cuales nos permitirán fijar ideas de lo que se hará en adelante. Es así como esto, nos proveerá de suficientes herramientas analíticas para poder mostrar en el segundo capítulo algunas aplicaciones importantes del Teorema del Punto Fijo de Banach (T.P.F.B.).

-1-

CAPÍTULO I DEFINICIONES Y TEOREMAS BÁSICOS I.1. DEFINICIÓN Sea X un conjunto no vacío y T : X → X una aplicación. Un punto x ∈ X se llama punto fijo de T si T ( x) = x i.e. x se mantiene fijo por T . Para cada aplicación T : X → X , podemos definir el conjunto de puntos fijos como

fix(T ) = { x ∈ X : x = T ( x)} . Nótese que el concepto de punto fijo es muy general, ya que no se requiere siquiera que X tenga alguna estructura, tan sólo se requiere un conjunto diferente del vacío. I.1.1. EJEMPLOS I.1.1.1.- La función f :

→

definida por f ( x ) = x tiene dos puntos fijos 0 y 1 . 2

I.1.1.2.- Definimos la función g : todo x ∈

→

por g ( x ) = x + a ; para todo a ∈

con a ≠ 0 y

llamada traslación. Es claro que una traslación no tiene puntos fijos.

I.1.1.3.- La rotación del plano tiene sólo un punto fijo, el centro de rotación. I.1.1.4.- Sea la función π :

( x, y ) ∈ ×

×

→

×

definida por π ( x, y ) = ( x,0) para cualquiera

llamada proyección. Recordemos que cualquier número real x es

identificado como un par ordenado ( x,0) en el plano. Es obvio que la proyección tiene infinitos puntos fijos. I.1.1.5.- Sea h :

− {0} →

− {0} definida por h( x) = x +

1 para todo x ∈ x

− {0} .

Aquí también es claro que h no tiene puntos fijos. I.2. DEFINICIÓN Sean ( X , d X ) y (Y , dY ) espacios métricos. Decimos que una aplicación T : X → Y es una contracción si existe un número real

α con 0 ≤ α < 1 tal que, para todo x, y ∈ X

dY (Tx, Ty ) ≤ α d X ( x, y ) . Decimos que

α es la constante de contracción de T .

Geométricamente, esto significa que cualesquiera que sean los puntos x e y estos tienen imágenes cercanas a estos puntos x e y ; más precisamente, la razón excede a la constante

α que es estrictamente menor que 1 . -2-

dY (Tx, Ty ) no d X ( x, y )

Cuando Y = X , decimos que T es una contracción en X . I.2.1. EJEMPLOS I.2.1.1.- Sea T :

→

definida por

T ( x) = 3 1 + x , x ∈ . Encontrar la solución de

T ( x) = x es equivalente a encontrar las raíces de la ecuación

x3 − x − 1 = 0 , x ∈ . Afirmación: T es una contracción en el subconjunto X = [1,2] porque

α < 1.

En efecto, observemos que para x, y ∈ X tenemos por el teorema del valor medio que existe c ∈ ]1,2[ tal que

T ( x) − T ( y ) = T ′(c)( x − y )

(I.2.1.1.A)

Aquí,

T ′( x) =

1 3 (1 + x ) 3

2

;

como

1< c < 2 entonces 3

3 4 3 (1 + c ) 9 < < , 3 3 3 2

por consiguiente

1 3 (1 + c ) 3

2


1 entonces

0 0 entonces T1 es creciente en [ a, b ] , luego T1 ( x) ∈ [ a, b ] si x ∈ [ a, b ] . Y por lo tanto

T1 ([ a, b ]) ⊂ [ a, b ] . Además, como

T1′ ( x) = 1 − β T ′( x) ≤ 1 − βα < 1 es decir tenemos

T1′ ( x) ≤ k , x ∈ [ a, b ] donde k = 1 − βα con 0 < k < 1 . Luego por el ejemplo I.2.1.3., tenemos que

T1 : [ a, b ] → [ a, b ] es una contracción en [ a, b ] . Luego por el T.P.F.B., existe un y sólo un x ∈ [ a, b ] punto fijo de T1 . Aún más, dado un punto x0 ∈ [ a, b ] tenemos que la sucesión iterada ( xn )n∈ con xn = T1 xn−1 ; n = 1, 2,… es tal que T1 x → x . n

-14-

I.4.6.7.- La ecuación

s = θ − ε senθ con

0 ≤ θ ≤ 2π , 0 ≤ ε < 1 , 0 < s < 2π

donde

ε es la excentricidad de la orbita de algún satélite, s es la anomalía media o

ángulo que recorre el satélite y θ es anomalía excéntrica; está ecuación nos permite calcular de forma aproximada la distancia de un satélite con respecto al Sol, llamada ecuación de Kepler. Dado s se desea hallar θ tal que

θ − ε senθ = s con este objetivo definimos

T (θ ) = θ − ε senθ − s , y el problema se reduce a hallar los ceros de la función T . En efecto, hacemos

T (0) = − s < 0 , T (2π ) = 2π − s > 0 , entonces Luego tomamos

T (0) < 0 < T (2π ) .

δ∈

tal que

ε < δ < 1 , con esto

T ′(θ ) = 1 − ε cosθ < 1 + δ . Y como tenemos

1 − ε < 1 − ε cosθ = T ′(θ ) , 1 − ε < T ′(θ ) < 1 + δ

entonces

0 < 1 − ε ≤ T ′(θ )
3. I.4.6.9.- Dado el conjunto

C = {( x, y ) ∈

2

y el sistema

-16-

: x ≤ 1, y ≤ 1}

⎧ sen( x) sen( y ) = 2 x ⎨ 2 2 ⎩x + y = 4 y Probaremos que ( 0,0 ) ∈

2

es la única solución en C del sistema dado.

En efecto, observemos que C es completo por ser un conjunto cerrado de y consideremos el campo vectorial f : C →

2

2

y

definido por

⎛ sen( x) sen( y ) x 2 + y 2 ⎞ f ( x, y ) = ⎜ , ⎟. 2 4 ⎝ ⎠ Es inmediato que f (C ) ⊂ C y que f es diferenciable; y la matriz jacobiana viene dada por

⎡cos( x) sen( y ) 1⎢ Df ( x, y ) = 2⎢ ⎢⎣ x

sen( x)cos( y ) ⎤ ⎥ ⎥ ⎥⎦ y

luego, hallamos

∑ ⎣⎡ D i, j

2 1 f ( x) ⎦⎤ = ( cos 2 ( x) sen 2 ( y ) + sen 2 ( x)cos 2 ( y ) + x 2 + y 2 ) 4

j i

≤

1 sen 2 ( y ) + sen 2 ( x) + x 2 + y 2 ) ( 4

≤

1 1 2 sen 2 (1) + 2 ) = ( sen 2 (1) + 1) < 1 . ( 4 2

Pues si para todo x ∈ [ −1,1] hacemos

g ( x) = sen 2 ( x) , entonces g alcanza su máximo en x = 1 . Por lo tanto por el ejemplo I.2.1.5. f es una contracción en C , luego sigue del T.P.F.B. que f tiene un único punto fijo, y en consecuencia ( 0,0 ) es la única solución del sistema considerado en C . I.4.6.10.- Sea A ⊂

n

un abierto convexo y sea f : A → A una aplicación continua.

⎧⎪

Supongamos que f es derivable en A y que sup ⎨

⎪⎩

-17-

∑ ⎡⎣ D i, j

⎫⎪ 2 ⎤ f ( x ) : x A ∈ ⎬ < 1. j i ⎦ ⎪⎭

Entonces f tiene un único punto fijo en A . En efecto, sean x, y ∈ A con x ≠ y aplicando el resultado del ejemplo I.2.1.5. a la función f en el segmento de recta [ x, y ] obtenemos

⎪⎧ f ( y ) − f ( x) 2 ≤ y − x 2 sup ⎨ ⎪⎩

∑ ⎡⎣ D i, j

2 ⎪⎫ f ( x) ⎤⎦ : x ∈ A⎬ . ⎪⎭

j i

Así, denotando

⎧⎪

α = sup ⎨ ⎩⎪

∑ ⎡⎣ D i, j

⎫⎪ 2 ⎤ f ( x ) : x A ∈ ⎬ < 1, j i ⎦ ⎭⎪

tenemos que

f ( y ) − f ( x ) 2 ≤ α y − x 2 ; x, y ∈ A . Por la continuidad de f , también tenemos

f ( y ) − f ( x ) 2 ≤ α y − x 2 ; x, y ∈ A . Entonces por el T.P.F.B. tenemos que existe un único punto fijo de f en A . I.4.6.11.- Demostraremos la existencia de una única función continua y acotada

y : [ 0, +∞[ →

que es solución de la ecuación integral

y (t ) = sen(t ) + ∫ e − s y ( set )ds , t ∈ [ 0, +∞[ . t

2

0

En efecto, consideremos

X = C * ([ 0, +∞[ ;

) = { f : [0, +∞[ →

; f es continua y acotada}

i.e. el espacio de las funciones continuas y acotadas en [ 0,+∞[ , para cualesquiera

f , g ∈ X se define la norma f − g = sup f ( x) − g ( x) . t∈[ 0,+∞[

(

Recordemos que X ,

) es un espacio de Banach.

Para todo x ∈ X y cada t ∈ [ 0, +∞[ , definimos la aplicación T : X → X por

(Tx ) (t ) = sen(t ) + ∫0 e− s x( set )ds . t

-18-

2

Afirmación 1: Tx ∈ X . i).- Tx es continua sobre [ 0,+∞[ pues la función seno es continua, y la función

t

t

∫e 0

− s2

x( set )ds es continua por el teorema fundamental del cálculo.

ii).- Tx es acotada en [ 0,+∞[ . En efecto, para toda x ∈ X y cada t ∈ [ 0, +∞[ , tenemos

(Tx ) (t ) ≤ sen(t ) + ∫0 e− s t

t

2

x( set ) ds

≤ 1 + M ∫ e − s ds ( x es acotada pues x ∈ X ) 2

0

∞

π

≤ 1 + M ∫ e − s ds = 1 + M 2

2

0

.

Afirmación 2: T es una contracción en X . En efecto, para cualesquiera u , v ∈ X , tenemos

(Tu ) (t ) − (Tv ) (t ) = =

(

t

)(

t

sen(t ) + ∫ e − s u ( set )ds − sen(t ) + ∫ e − s v( set )ds 2

0

t

2

0

t t −s ∫0 e ⎡⎣u(se ) − v(se ) ⎤⎦ ds t

2

≤ ∫ e − s u ( set ) − v( set ) ds 2

0

t

≤ sup u ( set ) − v( set ) ∫ e − s ds 0

t∈[ 0,+∞ )

≤ u−v

2

∫

∞

0

π

2

e − s ds =

π

i.e. (Tu ) (t ) − (Tv ) (t ) ≤

2

u−v

u − v ; t ∈ [ 0, +∞[ , u , v ∈ X .

2

Entonces

sup (Tu ) (t ) − (Tv ) (t ) ≤

t∈[ 0, +∞ )

π 2

u − v ; u, v ∈ X ,

por lo tanto

Tu − Tv ≤

π 2

-19-

u − v ; u, v ∈ X .

)

Luego, T es una contracción en el espacio de Banach X = C

*

([0, +∞[ ; ) , y

por lo tanto, por el T.P.F.B. se concluye la prueba. I.4.6.12.- Demostraremos que, para cualquier f ∈ X = C función y ∈ X tal que

y (t ) = f (t ) + ∫ *

([0, +∞[ ; ) existe una única

y ( s cos(t ) )

t

4 + s2

0

En efecto, en X = C

*

ds , t ≥ 0 .

([0, +∞[ ; ) para todo x ∈ X

y todo f ∈ X definimos la

aplicación T : X → X por

(Tx ) (t ) =

f (t ) + ∫

t

x ( s cos(t ) )

0

4 + s2

ds , t ≥ 0 .

Afirmación 1: Tx ∈ X . En efecto, para cada x ∈ X tenemos i).- Tx es continua sobre [ 0,+∞[ , pues f lo es por hipótesis, además la función

t

x( s cos(t ) ) ∫0 4 + s 2 ds es continua por el teorema fundamental del cálculo. t

ii).- Veamos que Tx es acotada en [ 0,+∞[ . Tenemos para todo t ≥ 0

(Tx ) (t ) ≤

f (t ) + ∫

t

0

x ( s cos(t ) ) 4 + s2

ds

ds ( f y x acotadas pues f , x ∈ X ) 0 4 + s2

≤ M + N∫

t

≤ M + N∫

∞

0

ds π ≤ M + N 4 + s2 4

Afirmación 2: T es una contracción en X . En efecto, para u , v ∈ X , tenemos

⎛

(Tu ) (t ) − (Tv ) (t ) = ⎜⎜ f (t ) + ∫0

t

⎝

=

u ( s cos(t ) ) 4 + s2

⎞ ⎛ ⎞ t v ( s cos(t ) ) ds ⎟ − ⎜ f (t ) + ∫ ds ⎟⎟ 2 0 ⎟ ⎜ + 4 s ⎠ ⎝ ⎠

1 ⎡u ( s cos(t ) ) − v ( s cos(t ) ) ⎤⎦ ds 0 4 + s2 ⎣

∫

t

-20-

≤ sup

t∈[ 0,+∞

{) u ( s cos(t ) ) − v ( s cos(t ) ) } ∫ 4 +1 s t

0

i.e. (Tu ) (t ) − (Tv ) (t ) ≤ Por lo tanto

Tu − Tv ≤

π 4

π 4

2

ds =

π 4

u−v .

u − v ; u, v ∈ X , t ≥ 0 .

u − v ; u, v ∈ X .

T es una contracción en X = C * ([ 0, +∞[ ;

) , la conclusión sigue del teorema del

punto fijo de Banach (T.P.F.B.). I.4.6.13.- Observemos que si X ⊂

n

es compacto y la aplicación f : X → X cumple

para arbitrarios x, y ∈ X con x ≠ y , con la condición

f ( x) − f ( y ) < x − y , entonces f posee un único punto fijo en X . En efecto, sea a ∈ X el punto donde la función continua

ϕ:X →

ϕ ( x) = x − f ( x) posee su mínimo c = a − f (a ) .

⋅ Si fuese c ≠ 0 , entonces tendríamos a ≠

f (a ) , de donde

f (a ) − f ( f (a )) < a − f (a) es decir, tenemos

ϕ ( f (a) ) < c ,

lo cual es una contradicción. Por lo tanto

c = 0, implicando que

a = f (a) , i.e. a es un punto fijo de f en X .

⋅ Si fuese a =

f (a ) y b = f (b) con a ≠ b resultaría a − b = f (a) − f (b) < a − b ,

lo cual es absurdo. Por lo tanto a es el único punto fijo de f en X . -21-

dada por

I.1.6.14.- Sean Λ un espacio topológico;

λ0 ∈ Λ , r > 0 y f dada por

f : (t , s, x, λ ) ∈ D = [ a, b ] × [ a, b ] × [ − r , r ] × Λ

f (t , s, x, λ ) ∈

una función continua tal que i). f (t , s,0, λ0 ) = 0 para cualquier s, t ∈ [ a, b ] ;

⎛ ∂f ⎞ ⎟ (t , s,0, λ0 ) = 0 ; s, t ∈ [ a, b ] . ∂ x ⎝ ⎠

⎛ ∂f ⎞ ⎟: D → ∂ x ⎝ ⎠

ii). existe ⎜

continua con ⎜

Entonces existe una vecindad Λ 0 de

λ0 tal que, para todo λ ∈ Λ 0 , la ecuación

integral

y (t ) = ∫ f [t , s, y ( s ), λ ] ds , t ∈ [ a, b ] b

a

tiene única solución uλ ∈ C

([ a, b],[ −r , r ]) , y la aplicación

λ ∈ Λ0

uλ ∈ C ( [ a , b ] , [ − r , r ] )

es continua. En efecto, consideremos

X = C ([ a, b ] , [ − r , r ]) = { f : [ a, b ] → [ − r , r ]; f es continua en [ a, b ] }, para arbitrarios f , g ∈ X definimos

d ( f , g ) = sup f ( z ) − g ( z ) ; z∈[ a ,b ]

recordemos que ( X , d ) es un espacio métrico completo. Por otro lado, dado c < 1 de i) y ii), sigue que existe una vecindad Λ 0 de tal que, para todo

λ0

λ ∈ Λ 0 , s, t ∈ [ a, b ] y cualquier x ∈ [ −r , r ] , tenemos

i’). f ( t , s, y ( s ), λ ) ≤

r , b−a

ii’).

∂f c . ( t , s , x, λ ) ≤ ∂x 2r (b − a )

Dado λ ∈ Λ 0 , u ∈ X y t ∈ [ a, b ] , definimos

(Tλ u ) (t ) = ∫a f [t , s, u ( s), λ ] ds ; b

se sigue que la función Tλ u es continua y de i’), se sigue que Tλ x ∈ X y de ii’), se sigue que la aplicación -22-

T : (λ , u ) ∈ Λ 0 × X

Tλ u ∈ X

es una contracción uniforme, pues dados u , v ∈ X , por la desigualdad del valor medio, tenemos

(Tλ u ) (t ) − (Tλ v ) (t ) ≤ ∫a f [t , s, u ( s), λ ] − f [t , s, v( s), λ ] ds b

b

≤ ∫ u ( s ) − v( s ) a

c ds ≤ c . 2r (b − a )

Por otro lado, es inmediato que, para cualquier

λ ∈ Λ0

u ∈ X , la aplicación

Tλ u ∈ X es continua. El resultado sigue del T.P.F.B.

I.5. INTERPRETACIÓN GEOMÉTRICA DEL T.P.F.B. Para enunciar el teorema del punto fijo de Banach para el caso real, usaremos el ejemplo I.2.1.3. teniendo el siguiente enunciado: con a ≠ b y g : [ a, b ] → [ a, b ] una función diferenciable tal que

“Sean a, b ∈

g ′( x) ≤ k < 1 , x ∈ [ a, b ] . Entonces: _

i).- Existe un único x ∈ [ a, b ] punto fijo de g . _

ii).- Cualquiera sea x0 ∈ [ a, b ] la sucesión iterada xn = g ( xn−1 ) ; n = 1, 2,... converge a x ”. Para interpretar geométricamente el T.P.F.B. para el caso real, observemos que el ítem i) _

_

_

nos dice que x = g ( x) es decir si y1 = x es la función identidad y y2 = g ( x) entonces el _

_

_

_ _

punto fijo x es un punto en el cual y1 = y2 , i.e., ( x, g ( x )) = ( x, x ) es la intersección de las gráficas de y1 y y2 . Para poder observar el comportamiento geométrico del T.P.F.B. en el plano, asumiremos que con la función continua g pueda ocurrir los siguientes dos casos; que ella sea creciente o decreciente. En ambos casos, como se observa en los gráficos I.5.A. y I.5.B., obtendremos una curva llamada quebrada iterada y una espiral hacia adentro (en este caso decimos que _

x es un punto atractor) respectivamente.

-23-

GRÁFICA I.5.A.

GRÁFICA I.5.B. -24-

1.5.1. OBSERVACIÓN 1.5.1.1.- En la observación I.2.3.1. vimos que toda contracción es uniformente continua, y por lo tanto continua. Ahora nos planteamos la siguiente interrogante: ¿Toda función continua es una contracción? Respuesta: ¡NO! + 0

En efecto, definimos la función f :

f ( x) =

→

1 , x∈ x

+ 0

por

+ 0

.

Observemos su gráfica:

de donde es claro que f es continua en Supongamos que f es una contracción en para arbitrarios x, y ∈

+ 0

+ 0. + 0,

entonces por definición tenemos que

se cumple

f ( x) − f ( y ) ≤ α x − y siempre con 0 ≤ α < 1 . En particular reemplazamos x =

1 1 , y = en la desigualdad anterior y obtenemos 3 5

⎛1⎞ f ⎜ ⎟− ⎝3⎠

1 1 ⎛1⎞ f ⎜ ⎟ ≤α − 3 5 ⎝5⎠ 3−5 ≤α -25-

2 15

entonces

α ≥ 15 ¡Absurdo! pues 0 ≤ α < 1 . Por lo tanto f no es una contracción en

+ 0.

1.5.1.2.- En la gráfica I.5.1.2. mostramos la gráfica de g en el caso cuando no se cumple la condición g ′( x) ≤ k < 1 , i.e. g no es una contracción en [ a, b ] . Esto desaparece la conclusión ii), por tanto existe el punto fijo pero no tenemos convergencia de la _

sucesión iterada al punto fijo x . 1.5.1.3.- En la gráfica I.5.1.3. mostramos la gráfica de la función f del ejemplo I.4.6.3., (ver página 13) tenemos el espacio métrico completo pues f no es una contracción en

, pero no tenemos punto fijo

.

1.5.1.4.- De las observaciones I.5.1.2. y I.5.1.3., podemos concluir que el teorema en general sólo se cumple si se tiene las hipótesis dadas, la carencia de una de ellas, colapsa sus conclusiones.

GRÁFICA I.5.1.2. -26-

GRÁFICA I.5.1.3. 1.5.1.5.- Recordemos un resultado que fue enunciado y probado en el capítulo I, en el ejemplo I.4.6.13. (ver página 21), el cual dice: “Si X ⊂

n

es compacto, y la aplicación f : X → X para cualesquiera x, y ∈ X

con x ≠ y , verifica

f ( x) − f ( y ) < x − y

(I.5.1.5.A)

entonces f posee un único punto fijo en X .” Surge la siguiente pregunta: ¿Si X no es compacto y verifica la desigualdad (I.5.1.5.A), f poseerá un único punto fijo en X ?. Respuesta : ¡NO! En efecto; recordemos que

f:

→

no es un compacto. Luego consideremos la función

definida por

f ( x) =

)

(

1 x + 1 + x2 , x ∈ . 2 -27-

Afirmación 1: f verifica la condición (I.5.1.5.A). En efecto; para arbitrarios x, y ∈

f ( x) − f ( y ) =

con x ≠ y , tenemos

) (

(

1 1 x + 1 + x2 − y + 1 + y2 2 2

(

=

1 1 ( x − y) + 2 2

≤

1 x − y + 1 + x2 − 1 + y2 2

1 + x2 − 1 + y2

(

1 ⎛⎜ = x− y + 2⎜ ⎝

(

) )

)

⎛ 1 + x2 + 1 + y 2 ⎞ ⎞ 1+ x − 1+ y ⎜ ⎟⎟ ⎜ 1 + x2 + 1 + y 2 ⎟ ⎟ ⎝ ⎠⎠ 2

2

x− y x+ y 1⎛ = ⎜ x− y + 2 ⎜⎝ 1 + x2 + 1 + y2

)

⎞ ⎟ ⎟ ⎠

⎛ ⎞ x+ y 1 = x − y ⎜1 + ⎟ 2 2 ⎟ ⎜ 2 x y 1 1 + + + ⎝ ⎠ Probaremos para arbitrarios x, y ∈

(I.5.1.5.B)

con x ≠ y , que

x+ y 1 + x2 + 1 + y 2

< 1.

(I.5.1.5.C)

Si aceptamos este resultado, tenemos

x + y < 1 + x2 + 1 + y2 x 2 + 2 xy + y 2 < 1 + x 2 + 2 (1 + x 2 )(1 + y 2 ) xy − 1 < (1 + x 2 )(1 + y 2 )

x 2 y 2 − 2 xy + 1 < 1 + y 2 + x 2 + x 2 y 2 ( x + y)2 > 0 Invirtiendo los pasos, verificamos la desigualdad (I.5.1.5.C). Observemos que para cualesquiera x, y ∈

-28-

con x ≠ y , descartamos que:

x+ y 1 + x2 + 1 + y 2

= 1,

pues si en particular tomamos x = 0 y y = 1 obtendríamos

1 =1 1+ 2 lo cual es un absurdo!. Todo lo anterior lo reemplazamos en (I.5.1.5.B), verificándose la afirmación 1. Afirmación 2: f no posee punto fijo en

.

En efecto; supongamos que f posee un punto fijo en tenemos que existe un w∈

, luego por definición

tal que

f ( w) = w

)

(

1 w + 1 + w2 = w 2 w + 1 + w2 = 2 w 1 + w2 = w

1 + w2 = w2 ¡absurdo! Notemos además que se cumple

f ( w) > w cuya prueba es análoga a la hecha para la desigualdad (I.5.1.5.C). Para fijar ideas, observemos nuevamente la gráfica I.5.1.3. de donde claramente podemos ver que las pruebas que hicimos analíticamente son correctas. Resaltemos que en esta demostración, hemos desarrollado un procedimiento alternativo de prueba a lo mostrado en el ejemplo I.4.6.13.

-29-

CAPÍTULO II APLICACIONES DEL TEOREMA DEL PUNTO FIJO DE BANACH II.1.- APLICACIÓN A LOS SISTEMAS DE ECUACIONES LINEALES ALGEBRAICAS El T.P.F.B. tiene importantes aplicaciones en los métodos de iteración para resolver sistemas de ecuaciones lineales algebraicas y proporciona suficientes condiciones para la convergencia y cota de error. n

Consideremos el espacio euclidiano n-dimensional

y el sistema de ecuaciones lineales

algebraicas n

∑a x ij

j

= yi ; i = 1,2,3,......, n

j=1

que en notación matricial es

y = Ax ,

donde A = ⎡⎣ a ij ⎤⎦ es una matriz n × n , x = [ xi ] y y = [ yi ] son matrices n × 1 . Notemos que el problema de resolver y = Ax se puede reemplazar por el problema equivalente de encontrar los puntos fijos del operador T :

n

→

n

definido por

Tx = ( I − A) x + y ; hacemos

B = ( I − A) = ⎡⎣bij ⎤⎦ la cual es una matriz n × n real fija ( I es la matriz identidad n × n ), es decir, tenemos

Tx = Bx + y . ¿Bajo qué condiciones T es una contracción en de la elección de la métrica en

n

n

?. La respuesta a esta pregunta depende

con la que se va a trabajar. Dediquémonos a encontrar

tales condiciones. Sean

u = ( u1 ,

, un ) , v = ( v1 ,

, vn ) , p = ( p1 ,

, pn ) , q = ( q1 ,

, qn ) ∈

supongamos que p = Tu (i.e. p = Tu = Bu + y ) y q = Tv (i.e. q = Tv = Bv + y ).

i Con la métrica del máximo, tenemos ⎛ n ⎞ ⎛ n ⎞ d ∞ (Tu , Tv ) = d ∞ ( p, q ) = max pi − qi = max ⎜ ∑ bij u j + yi ⎟ − ⎜ ∑ bij v j + yi ⎟ i i ⎝ j =1 ⎠ ⎝ j =1 ⎠

-30-

n

y

= max ∑ bij ( u j − v j ) ≤ max ∑ bij u j − v j n

i

n

i

j =1

(

j =1

)

n ⎛ ⎞ ≤ ⎜ max ∑ bij ⎟ max u j − v j = λ∞ d ∞ (u , v) ⎝ i j =1 ⎠ j n

Entonces, para que T sea una contracción en

debemos exigir que

n

λ∞ = max ∑ bij < 1 . i

(II.1.A)

j =1

i Con la métrica de la suma, tenemos n

d1 (Tu , Tv) = d1 ( p, q ) = ∑ i =1

⎛ n ⎞ ⎛ n ⎞ pi − qi = ∑ ⎜ ∑ bij u j + yi ⎟ − ⎜ ∑ bij v j + yi ⎟ i =1 ⎝ j =1 ⎠ ⎝ j =1 ⎠ n

n n n ⎛ ⎞⎛ n ⎞ = ∑ ∑ bij ( u j − v j ) ≤ ⎜ max ∑ bij ⎟⎜ ∑ u j − v j ⎟ = λ1d1 (u, v) . i =1 j =1 ⎝ i j =1 ⎠⎝ j =1 ⎠ n

Así, el operador T es una contracción en

si

n

λ1 = max ∑ bij < 1 . i

(II.1.B)

j =1

i Por último, considerando la métrica euclidiana, tenemos

[ d 2 (Tu, Tv)] = [ d 2 ( p, q)] 2

2

n

= ∑ ( pi − qi ) i =1

2

⎧⎪⎛ n ⎞ ⎛ n ⎞ ⎫⎪ = ∑ ⎨⎜ ∑ bij x j + yi ⎟ − ⎜ ∑ bij y j + yi ⎟ ⎬ i =1 ⎪ ⎠ ⎝ j =1 ⎠ ⎪⎭ ⎩⎝ j =1 n

2

2

⎧ n ⎫ ⎛ n n 2⎞ 2 = ∑ ⎨∑ bij ( u j − v j ) ⎬ ≤ ⎜ ∑∑ bij ⎟ ⎡⎣ d ( u , v ) ⎤⎦ i =1 ⎩ j =1 ⎭ ⎝ i =1 j =1 ⎠ n

⎛ ⎜ ⎝

i.e. d 2 (Tu , Tv ) ≤ ⎜

⎞ 2 b ⎟ d ( u, v ) ; u, v ∈ ∑∑ ij ⎟ i =1 j =1 ⎠ n

n

n

.

(En la desigualdad anterior hemos utilizado la desigualdad de Schwartz). Análogamente a lo hecho antes, aquí exigimos que n

n

∑∑ b i =1 j =1

2

ij

0 tal que K ( t0 , s, x( s ) ) ≤ N

(

)

para todo ( s, x( s ) ) ∈ [ a, b ] × x [ a, b ] .

-42-

Con estos resultados, tenemos la siguiente desigualdad

(Tx ) (t ) − (Tx ) (t0 ) ≤ de donde es claro que

y (t ) − y (t0 ) + λ b − a

ε λ b−a

+ N t − t0

(Tx ) (t ) → (Tx ) (t0 )

para todo x ∈ X . Por lo tanto, concluimos que la afirmación 1 es válida i.e. Tx ∈ X Afirmación 2: Existe m ≥ 1 tal que T

m

.

es una contracción en X .

En efecto; vamos a demostrar inicialmente que dados u , v ∈ X y t ∈ [a, b] se cumple

(T u ) (t ) − (T v ) (t ) ≤ n

n

λ M n (t − a ) n n

n!

u−v

(II.2.3.1)

La demostración es hecha por inducción matemática.

i Para n = 1 el resultado es inmediato, es decir tenemos

(Tu ) (t ) − (Tv ) (t ) ≤

λ M (t − a) 1!

u−v .

.

i Admitiendo el resultado para n vamos a demostrarlo para n + 1

(T

u ) (t ) − (T n+1v ) (t ) = T (T nu ) (t ) − T (T n v ) (t )

n +1

(

)(

= y (t ) + λ ∫ K ⎡⎣t , s, (T nu ) ( s ) ⎤⎦ ds − y (t ) + λ ∫ K ⎡⎣t , s, (T nv ) ( s ) ⎤⎦ ds a a t

{

t

}

= λ ∫ K ⎡⎣t , s, (T nu ) ( s) ⎤⎦ − K ⎡⎣t , s, (T n v ) ( s ) ⎤⎦ ds a t

t

K ⎡⎣t , s, (T nu ) ( s ) ⎤⎦ − K ⎡⎣t , s, (T nv ) ( s ) ⎤⎦ ds

≤λ

∫

≤λ

∫ M (T u ) (s) − (T v ) ( s) ds

a t

n

n

a

t

λ M n (s − a)n

a

n!

n

≤λ

∫M

=λ

n +1

M

n +1

u − v ds

(t − a) n+1 u−v . (n + 1)!

Luego (II.2.3.1) es cierto, además de (II.2.3.1) tenemos

(T u ) (t ) − (T v ) (t ) ≤ n

n

λ M (t − a) n!

n

d ( u , v ) ; u , v ∈ X y t ∈ [ a , b] ,

-43-

)

pero aun más, como

λ M (t − a)

n

≤

n!

λ M (b − a) n!

n

,

lo que implica

(T u ) (t ) − (T v ) (t ) ≤ n

n

λ M (b − a)

n

d ( u , v ) ; u , v ∈ X y t ∈ [ a , b] ,

n!

entonces

sup (T u ) (t ) − (T v ) (t ) ≤ n

n

λ M (b − a)

n

n!

a ≤t ≤ b

d ( u, v ) ; u, v ∈ X .

Por lo tanto

d (T u , T v ) ≤ n

Para cada n ∈

+

n

λ M (b − a )

n

n!

d ( u, v ) ; u, v ∈ X .

, sea

an =

λ M (b − a)

n

n!

ahora utilizamos el criterio de la razón para sucesiones de números reales:

λ M (b − a) lim n →∞

n +1

λ M (b − a) an+1 (n + 1)! = lim = lim = 0 < 1, n n →∞ n→∞ an n + 1 λ M (b − a) n!

entonces

lim

λ M (b − a )

n →∞

n!

n

= 0,

⎧⎪ λ M (b − a ) n ⎫⎪ por lo tanto ⎨ ⎬ converge, luego sigue por definición de límite de sucesiones n ! ⎩⎪ ⎭⎪n≥1 de números reales que: Dado 0 < ε < 1 , existe n0 ∈

tal que para todo n > n0 tenemos

λ M (b − a) n!

n

< ε < 1.

Con lo que se verifica la afirmación 2. Entonces tomamos m = n0 , para que exista m ≥ 1 tal que T

m

es una contracción en C ([a, b],

) = X . Y la conclusión de este teorema sigue del

corolario 1.4.2. ■ -44-

II.2.4. OBSERVACIÓN Hacemos

λ = 1 en la ecuación integral de Fredholm de segunda especie, y tenemos b

x(t ) = y (t ) + ∫ K (t , s ) x ( s )ds

(II.2.4.1)

a

II.2.5. TEOREMA Supongamos que K : [ a, b ] × [ a, b ] →

sup a ≤t ≤ b

y y : [ a, b ] →

x : [ a, b ] →

{∫

b

a

es una función continua tal que satisface

}

K (t , s) ds < 1

(II.2.5.1)

una función continua. Entonces existe una y sólo una función continua que satisface la ecuación (II.2.4.1).

Demostración: Consideremos el espacio normado X = C ∞

([ a, b], ) el cual con la métrica uniforme

es completo; y definamos la aplicación T : X → X por

x

Tx

(Tx ) (t ) = y (t ) + ∫a K (t , s) x( s)ds . b

Análogamente a lo hecho en la afirmación 1 del Teorema II.2.3. se prueba que Tx ∈ X . Veamos que T es una contracción en X . Para cualesquiera x1 , x2 ∈ X tenemos

Tx1 − Tx2

∞

= sup (Tx1 ) (t ) − (Tx2 ) (t ) a ≤ t ≤b

(

)(

b

b

= sup y (t ) + ∫ K (t , s) x1 ( s )ds − y (t ) + ∫ K (t , s ) x2 ( s )ds a ≤t ≤ b

= sup

a ≤ t ≤b

a

∫

b

a

a

K (t , s) ( x1 ( s ) − x2 ( s ) ) ds

b

≤ sup ∫ K (t , s ) x1 ( s ) − x2 ( s ) ds a ≤ t ≤b a

≤ x1 − x2

b

sup ∫ K (t , s) ds ≤ c x1 − x2 ∞ a ≤t ≤b a

es decir

Tx1 − Tx2

∞

≤ c x1 − x2 ∞ ; x1 , x2 ∈ X -45-

∞

)

siempre que

c = sup

a ≤ t ≤b

Luego T es una contracción en C

{∫

b

a

}

K (t , s) ds < 1 .

([ a, b], ) = X , y la afirmación de este teorema sigue del

T.P.F.B. ■ Dado un x0 ∈ C

([ a, b], ) . El T.P.F.B. nos permite obtener el punto fijo x como un límite x = lim T n x0 .

(II.2.5.2)

n→∞

Es interesante reinterpretar este límite como una serie. Definimos la aplicación

k : C ([ a, b ] ,

) → C ([ a, b], ) por b

kx = ∫ K (t , s ) x( s )ds . a

Esta aplicación es denominada operador de Fredholm, y la función k es llamada el núcleo de K . La ecuación integral (II.2.4.1) puede ser reescrita como

(I − k)x = y

(II.2.5.3)

donde I es el operador identidad. La aplicación T (donde Tx = x ) esta definida por

Tx = y + kx , lo implica que

T n x0 = T n−1 (Tx0 ) = T n−1 ( y + kx0 ) = y + kT n−1 ( x0 ) = y + kT n−2 (Tx0 ) = y + kT n−2 ( y + kx0 ) = y + k ⎡⎣ y + kT n−2 x0 ⎤⎦ = y + k ⎡⎣ y + kT n−3 (Tx0 ) ⎤⎦ = y + k ⎡⎣ y + k ( y +

T n x0 = y + ky +

+ k ( y + kx0 ) ) ⎤⎦

+ k n y + k n+1 x0 .

Usando la ecuación (II.2.5.2), obtenemos n +1

∞

i =0

n =0

x = lim ∑ k i y = ∑ k n y , n→∞

Por otro lado de (II.2.5.3) tenemos

x = (I − k) y , −1

-46-

entonces de las dos igualdades anteriores, conseguimos ∞

(I − k ) = ∑kn . −1

(II.2.5.4)

n =0

Esta serie es llamada serie de Neumann. El uso de sumas parciales de esta serie para aproximar a la inversa se llama aproximación de Born. Explícitamente tenemos

⎛ ∞ ⎞ x(t ) = ⎜ ∑ k n y ⎟ (t ) = ( y + ky + k 2 y + ⎝ n =0 ⎠

) (t ) = y(t ) + ( ky ) (t ) + ( k y ) (t ) + 2

= y (t ) + ∫ K (t , s ) y ( s )ds + ( k (ky ) ) (t ) + b

a

= y (t ) + ∫ K (t , s ) y ( s )ds + ∫ K (t , s ) ( (ky )( s ) ) ds + b

b

a

a

b

b

a

a

b

b

a

a

= y (t ) + ∫ K (t , s ) y ( s )ds + ∫ K (t , s ) x(t ) = y (t ) + ∫ K (t , s ) y ( s )ds + ∫

∫

b

a

( ∫ K (s, r) y(r)dr ) ds + b

a

K (t , s ) K ( s, r ) y (r )drds +

La serie de Neumann es similar a la serie geométrica. De hecho, la igualdad (II.2.5.4) es realmente una serie geométrica que es absolutamente convergente con respecto a determinada norma del operador cuando k

∞

< 1 . Esto explica porque no necesitamos imponer ninguna

condición sobre y para que la ecuación (II.2.5.1) sea una condición que asegure que

( I − k ) es invertible y depende sólo de la función k

-47-

(ver [13]) .

II.3.- APLICACIÓN A LAS ECUACIONES DIFERENCIALES ORDINARIAS II.3.1. DEFINICIÓN Sean Ω ⊂

y f ∈ C (Ω, ) . Decimos que una

×

diferenciable

función continuamente

u : [ c, d ] →

es una solución de la ecuación diferencial ordinaria

y′ = f (t , y )

(II.3.1.1)

si para todo t ∈ [ c, d ] tenemos que

(t , u (t )) ∈ Ω y u′(t ) = f (t , u (t )) . II.3.2. PROPOSICIÓN Sea (t0 , y0 ) ∈Ω . La condición necesaria y suficiente para que u : [ c, d ] → continua con (t , u (t )) ∈ Ω , para todo t ∈ [ c, d ] sea solución de la ecuación diferencial (II.3.1.1) satisfaciendo u (t0 ) = y0 , es que u sea solución continua de la ecuación integral t

y (t ) = y0 + ∫ f ( s, y ( s ))ds

(II.3.2.1)

t0

Demostración: Es inmediato, pues, si f y u son continuas, la función t ∈ [ c, d ] → f (t , u (t )) ∈

,

también lo es. Basta entonces aplicar el Teorema Fundamental del Cálculo. ■ II.3.3. TEOREMA (Teorema de existencia de Cauchy para ecuaciones diferenciales) Sean (t0 , y0 ) ∈

y f : [t0 − a, t0 + a ] × Bb [ y0 ] ⊂

×

continua y lipschitziana en la segunda variable, esto es, existe L ∈

×

→

+

tal que

una función

f (t , z1 ) − f (t , z2 ) ≤ L z1 − z2 ; t ∈ [t0 − a, t0 + a ] , z1 , z2 ∈ Bb [ y0 ]

(II.3.3.1)

Entonces existe a con 0 < a < a tal que la ecuación diferencial *

*

y′ = f (t , y ) * * tiene una única solución u definida en ⎡⎣t0 − a , t0 + a ⎤⎦ satisfaciendo u (t0 ) = y0 .

Demostración: Observemos que f está definida en un dominio compacto, entonces siendo f continua,

f es acotada i.e. existe M > 0 tal que f (t , y ) ≤ M , para ( t , y ) ∈ [t0 − a, t0 + a,] × Bb [ y0 ] . -48-

(Notemos que si M = 0 resulta f = 0 y la afirmación es trivial). Tomamos

⎧ b⎫ a* = min ⎨a, ⎬ . ⎩ M⎭

(II.3.3.2)

Para simplificar la notación, supongamos que tenemos a ≤

b * i.e. a = a . M

Recordemos que X = C ([t0 − a, t0 + a ] , Bb [ y0 ]) dotado de la métrica del supremo *

d (u , v) = u − v =

sup

t∈[t0 − a ,t0 + a ]

{ u (t ) − v(t ) }

es un espacio métrico completo. Sea u ∈ X , definimos Tu : [t0 − a, t0 + a,] →

por

t

(Tu )(t ) = y0 + ∫ f ( s, u ( s ))ds . t0

Es inmediato que Tu es continua. De (II.3.3.2) se tiene que, para todo t con t − t0 ≤ a

(Tu )(t ) − y0 ≤

∫

t

t0

f ( s, u ( s )) ds ≤ M t − t0 ≤ Ma ≤ b ,

entonces Tu ∈ X . Por otro lado, para cualesquiera u , v ∈ X , tenemos

(Tu )(t ) − (Tv)(t ) = ≤ ≤

∫

t

∫

t

∫

t

t0

t0

t0

≤L

t

f ( s, u ( s ))ds − ∫ f ( s, v( s ))ds t0

f ( s, u ( s )) − f ( s, v( s )) ds L u ( s ) − v( s ) ds

sup

t∈[t0 − a ,t0 + a ]

{ u (t ) − v(t ) } t − t

i.e. (Tu )(t ) − (Tv)(t ) ≤ L u − v t − t0 . Además de esto, existe m ∈

+

tal que T

m

0

(II.3.3.3)

es una contracción en X , pues, para arbitrarios

u , v ∈ X tenemos Lm t − t0 (T u )(t ) − (T v)(t ) ≤ m! m

m

-49-

m

u−v .

(II.3.3.4)

En efecto, para arbitrarios u , v ∈ X y para todo t con t − t0 ≤ a , se tiene

(T 2u )(t ) − (T 2v)(t ) = T (Tu )(t ) − T (Tv )(t )

= ≤ ≤

∫ [ f (s,Tu (s)) − f (s,Tv(s))] ds t

t0

∫

t

∫

t

t0

t0

f ( s, Tu ( s )) − f ( s, Tv( s )) ds L Tu ( s ) − Tv( s ) ds

∫

≤ L2 u − v

t

t0

s − t0 ds (viene de II.3.3.3) 2

s − t0 t = L u−v |t0 2 2

t − t0 = L u−v 2

2

2

t − t0 i.e. (T u )(t ) − (T v)(t ) ≤ L 2! 2

2

2

2

u − v ; u, v ∈ X

y, de modo más general

(T m+1u )(t ) − (T m+1v)(t ) ≤ (T mu )(t ) − (T mv)(t ) L

t − t0 ; u, v ∈ X , m +1

de donde resulta (II.3.3.4). Como t − t0 ≤ a tomando el supremo en (II.3.3.4), tenemos

T mu − T m v ≤ Para cada m ∈

+

Lm a m u−v . m!

, sea

um

( La ) = m!

m

,

utilizando el criterio de la razón para sucesiones de números reales

( La ) ( m + 1)! = La lim ⎛ 1 ⎞ = 0 < 1 ( ) m→∞ ⎜ ⎟ m ⎝ m +1⎠ ( La ) m +1

lim

m→∞

um+1 = lim m→∞ um

m! -50-

entonces

lim um

m→∞

( La ) = lim m→∞

m!

m

= 0,

⎛ ( La )m ⎞ por lo tanto ⎜ ⎟ converge, luego sigue por definición de límite de sucesiones de ⎜ m! ⎟ ⎝ ⎠ m≥1 números reales que: Dado 0 < ε < 1 , existe n0 ∈

tal que para todo n > n0 tenemos

( La )

m

m!

Entonces basta tomar m = n0 , para concluir que existe m ≥ 1 tal que T

m

< ε < 1. es una contracción

en C ([t0 − a, t0 + a ] , Bb [ y0 ]) = X , luego por el corolario I.4.2., T tendrá un único punto *

fijo, el cual es la solución de la ecuación integral (II.3.2.1) y por tanto solución de la ecuación diferencial ordinaria (II.3.1.1). ■

-51-

II.4.- APLICACIÓN A LAS ECUACIONES DIFERENCIALES PARCIALES II.4.1. EJEMPLO: Sean f , g :

+

×

funciones continuas y para cada ( x, t ) ∈

→

×

+

considere-

mos la ecuación diferencial en derivadas parciales

∂ 2u ∂ 2u ( x, t ) − 2 ( x, t ) − g ( x, t )u ( x, t ) = f ( x, t ) ∂t 2 ∂x

(II.4.1.1)

con las condiciones iniciales

u ( x,0) = 0 ,

∂u ( x,0) = 0 . ∂t

(II.4.1.2)

Observemos que (II.4.1.1) corresponde al problema de la ecuación de la onda con vibraciones forzadas (semilineal hiperbólica). II.4.2. PROPOSICIÓN

×

La función u :

iniciales (II.4.1.2) con u ,

+

→

es solución de (II.4.1.1), satisfaciendo las condiciones

∂u ∂u ∂ 2u ∂ 2u ∂ 2u , , , , continuas en ∂t ∂x ∂x∂t ∂t 2 ∂x 2

×

+

si sólo si u es

solución de la ecuación integral

u ( x, t ) =

x + ( t −τ ) 1 t d τ [ f (σ ,τ ) + g (σ ,τ )u (σ ,τ )]dσ 2 ∫0 ∫x −( t −τ )

(II.4.2.1)

Demostración: Llamemos F = f + gu , para (II.4.1.1) hacemos el cambio de variable

ξ = x+t , η = x−t , de donde

x=

ξ +η 2

,t=

ξ −η 2

.

∂u ∂u ∂ 2u ∂ 2u ∂ 2u Entonces, ahora las nuevas hipótesis son u , , , , , son funciones ∂ξ ∂η ∂ξ∂η ∂ξ 2 ∂η 2 continuas en

×

+

(Observar que estas implican que

Queremos hallar expresiones para

∂ 2u ∂ 2u ). = ∂ξ∂η ∂η∂ξ

∂ 2u ∂ 2u , para lo cual consideramos el diagrama de ∂t 2 ∂x 2

dependencia: -52-

Luego, con esto hallemos:

i

∂u ∂u ∂ξ ∂u ∂η ∂u ∂u ∂u = + i.e. = − ∂t ∂ξ ∂t ∂η ∂t ∂t ∂ξ ∂η 1

(II.4.2.2)

−1

∂ 2u ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ∂u ⎞ ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ⎞ = ⎜ ⎟= − = − = ⎜ ⎟− ⎜ ⎟ ∂t 2 ∂t ⎝ ∂t ⎠ ∂t ⎜⎝ ∂ξ ∂η ⎟⎠ ∂t ⎜⎝ ∂ξ ⎟⎠ ∂t ⎜⎝ ∂η ⎟⎠ ∂ξ ⎝ ∂t ⎠ ∂η ⎝ ∂t ⎠ (Pues de ambas hipótesis sobre u y sus derivadas parciales, tenemos

∂ 2u ∂ 2u ∂ 2u ∂ 2u y ). = = ∂t ∂ξ ∂ξ∂t ∂t ∂η ∂η∂t =

∂ ⎛ ∂u ∂u ⎞ ∂ ⎛ ∂u ∂u ⎞ ∂ 2u ∂ 2u ∂ 2u ∂ 2u − − − = − − + ∂ξ ⎜⎝ ∂ξ ∂η ⎟⎠ ∂η ⎜⎝ ∂ξ ∂η ⎟⎠ ∂ξ 2 ∂ξ∂η ∂η∂ξ ∂η 2 ∂ 2u ∂ 2u ∂ 2u ∂ 2u = + −2 ∂t 2 ∂ξ 2 ∂η 2 ∂ξ∂η

i

∂u ∂u ∂ξ ∂u ∂η ∂u ∂u ∂u = + i.e. = + ∂x ∂ξ ∂x ∂η ∂x ∂x ∂ξ ∂η 1

(II.4.2.3)

(II.4.2.4)

1

∂ 2u ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ∂u ⎞ ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ⎞ = = ⎜ ⎟= + = + ⎜ ⎟+ ⎜ ⎟ ∂x 2 ∂x ⎝ ∂x ⎠ ∂x ⎜⎝ ∂ξ ∂η ⎟⎠ ∂x ⎜⎝ ∂ξ ⎟⎠ ∂x ⎜⎝ ∂η ⎟⎠ ∂ξ ⎝ ∂x ⎠ ∂η ⎝ ∂x ⎠ (Análogamente tenemos

∂ 2u ∂ 2u ∂ 2u ∂ 2u = = y ). ∂x∂η ∂η∂x ∂x∂ξ ∂ξ∂x

∂ 2u ∂ 2u ∂ 2u ∂ ⎛ ∂u ∂u ⎞ ∂ ⎛ ∂u ∂u ⎞ ∂ 2u = + + + = + + + ∂ξ ⎜⎝ ∂ξ ∂η ⎟⎠ ∂η ⎜⎝ ∂ξ ∂η ⎟⎠ ∂ξ 2 ∂ξ∂η ∂η∂ξ ∂η 2 ∂ 2u ∂ 2u ∂ 2u ∂ 2u . = + +2 ∂x 2 ∂ξ 2 ∂η 2 ∂ξ∂η Luego reemplazamos (II.4.2.3) y (II.4.2.5) y F en (II.4.1.1), obteniendo

-53-

(II.4.2.5)

∂ 2u ⎛ ξ + η ξ − η ⎞ 1 ⎛ ξ +η ξ −η ⎞ , , ⎜ ⎟=− F⎜ ⎟. 2 ⎠ 4 ⎝ 2 2 ⎠ ∂ξ∂η ⎝ 2 Ahora integramos con respecto a ξ

∫γ

ξ

∂ 2u ⎛ ψ + η ψ − η ⎞ 1 ξ ⎛ψ + η ψ −η ⎞ , d F⎜ , ψ = − ⎜ ⎟ ⎟ dψ 2 ⎠ 4 ∫γ ⎝ 2 2 ⎠ ∂ψ∂η ⎝ 2

∂u ⎛ ξ + η ξ − η ⎞ ∂u ⎛ γ + η γ − η ⎞ 1 ξ ⎛ψ + η ψ −η ⎞ , − , = − F⎜ , ⎜ ⎟ ⎜ ⎟ ⎟ dψ . ∂η ⎝ 2 2 ⎠ ∂η ⎝ 2 2 ⎠ 4 ∫γ ⎝ 2 2 ⎠ Aquí hacemos

γ = η y tenemos

∂u ⎛ ξ + η ξ − η ⎞ ∂u 1 ξ ⎛ψ + η ψ −η ⎞ η , = ,0 − F⎜ , ( ) ⎜ ⎟ ⎟ dψ ∂η ⎝ 2 2 ⎠ ∂η 4 ∫γ ⎝ 2 2 ⎠

(II.4.2.6)

Restando (II.4.2.4) de (II.4.2.2), obtenemos

∂u 1 ∂u 1 ∂u = − ∂η 2 ∂x 2 ∂t entonces

∂u 1 ∂u 1 ∂u (η ,0) = (η ,0) − (η ,0) . ∂η 2 ∂x 2 ∂t

Esta última relación la reemplazamos en (II.4.2.6), y conseguimos

∂u ⎛ ξ + η ξ − η ⎞ 1 ∂u 1 ∂u 1 ξ ψ +η ψ −η ⎞ , , (η ,0 ) − (η ,0 ) − ∫γ F ⎛⎜ ⎜ ⎟= ⎟ dψ . ∂η ⎝ 2 2 ⎠ 2 ∂x 2 ∂t 4 2 2 ⎝ ⎠ Integramos nuevamente con respecto a

η de un valor arbitrario de η hasta η = ξ , y

tenemos ξ

∫η

ξ ⎡ 1 ∂u ∂u ⎛ ξ + s ξ − s ⎞ 1 ∂u 1 ξ ξ ⎛ψ + s ψ − s ⎞ ⎤ , ds = s ,0 − s ,0 ds − F⎜ , ( ) ( ) ⎜ ⎟ ⎟ dψ ds ∫η ⎣⎢ 2 ∂x ∂s ⎝ 2 2 ⎠ 2 ∂t 4 ∫η ∫s ⎝ 2 2 ⎠ ⎦⎥

1 ∂u 1 ξ ξ ⎛ψ + s ψ − s ⎞ ⎛ ξ + η ξ − η ⎞ ξ ⎡ 1 ∂u ⎤ u (ξ ,0 ) − u ⎜ , = s ,0 − s ,0 ds − F⎜ , ( ) ( ) ⎟ ⎟ dψ ds 2 ⎠ ∫η ⎣⎢ 2 ∂x 2 ∂t 4 ∫η ∫s ⎝ 2 2 ⎠ ⎝ 2 ⎦⎥ 1 ξ ∂u 1 ξ ∂u 1 ξ ξ ψ +s ψ −s⎞ ⎛ ξ +η ξ −η ⎞ u⎜ , , ( s,0 ) ds + ∫η ( s,0 ) ds + ∫η ∫s F ⎛⎜ ⎟ = u (ξ ,0 ) − ∫η ⎟ dψ ds 2 ⎠ 2 ∂x 2 ∂t 4 2 ⎠ ⎝ 2 ⎝ 2

1 1 1 ξ ∂u 1 ξ ξ ⎛ψ + s ψ − s ⎞ = u (ξ ,0 ) − u (ξ ,0 ) + u (η ,0 ) + ∫ s,0 ) ds ∫ ∫ F ⎜ , ( ⎟ dψ ds 2 2 2 η ∂t 4 η s ⎝ 2 2 ⎠ 1 ξ ∂u 1 ξ ξ ⎛ψ + s ψ − s ⎞ ⎛ ξ +η ξ −η ⎞ 1 ⎡ ⎤ u⎜ , u ,0 u ,0 s ,0 ds F⎜ , ξ η = + + + ( ) ( ) ( ) ⎟ ⎟ dψ ds ⎦ 2 ∫η ∂t 2 ⎠ 2⎣ 4 ∫η ∫s ⎝ 2 2 ⎠ ⎝ 2 -54-

Ahora, aplicamos las condiciones iniciales (II.4.1.2) quedándonos

⎛ ξ + η ξ −η ⎞ 1 ξ ξ ⎛ψ + s ψ − s ⎞ u⎜ , F⎜ , ⎟= ⎟ dψ ds 2 ⎠ 4 ∫η ∫s ⎝ 2 2 ⎠ ⎝ 2

(II.4.2.7)

Para la integral doble, hacemos el cambio de variable

ψ = σ +τ , s = σ −τ luego por el teorema de cambio de variable para integrales dobles tenemos

∂ψ ∂ (ψ , s ) ∂σ ) dσ dτ = drds = det( J (ψ , s )) dσ dτ = det( ∂s ∂ (σ ,τ ) ∂σ 1 1 = dσ dτ = −2 dσ dτ = 2dσ dτ . 1 −1

∂ψ ∂τ dσ dτ ∂s ∂τ

Observemos además que el dominio de integración inicial

η ≤ s ≤ψ ≤ ξ se transforma en

η +τ ≤ σ ≤ ξ −τ , 0 ≤τ ≤

1 (ξ -η ) . 2

Luego de reemplazar en (II.4.2.7), tenemos ξ −η ⎛ ξ + η ξ − η ⎞ 1 2 ξ −τ u⎜ , F (σ ,τ ) dσ dτ ⎟= 2 ⎠ 2 ∫0 ∫η +τ ⎝ 2

ahora tenemos presente que

u ( x, t ) =

ξ = x + t , η = x − t , F = f + gu obteniendo

x + ( t −τ ) 1 t d τ [ f (σ ,τ ) + g (σ ,τ )u (σ ,τ )] dσ , 2 ∫0 ∫x −( t −τ )

y por lo tanto, sigue que u , satisfaciendo (II.4.1.1) y (II.4.1.2), también satisface (II.4.2.1). Recíprocamente, veamos que (II.4.2.1) i.e.

u ( x, t ) =

x + ( t −τ ) 1 t d τ [ f (σ ,τ ) + g (σ ,τ )u (σ ,τ )] dσ 2 ∫0 ∫x −( t −τ )

satisface (II.4.1.1) y las condiciones iniciales (II.4.1.2).

i u ( x,0) = 0

i Probaremos que

∂u ( x,0) = 0 . ∂t -55-

En efecto, sea

u ( x, t ) =

1 t G ( x, t ,τ )dτ 2 ∫0

(II.4.2.A)

donde

G ( x, t ,τ ) = ∫

x +( t −τ )

x −( t −τ )

F (σ ,τ )dσ .

En esta parte de la demostración ignoramos la constante

∂u 1 ( x, t ) = lim [u ( x, t + h) − u ( x, t ) ] h →0 h ∂t

1 ; entonces 2

t 1 t +h = lim ⎡ ∫ G ( x, t + h,τ )dτ − ∫ G ( x, t ,τ )dτ ⎤ ⎥⎦ 0 h →0 h ⎢ ⎣ 0

1 = lim ⎡ h →0 h ⎢ ⎣

( ∫ G( x,t + h,τ )dτ + ∫ t

t +h

0

0

)

G ( x, t + h,τ )dτ − ∫ G ( x, t ,τ )dτ ⎤ ⎥⎦ 0 t

t 1 t 1 t +h = lim ⎡ ∫ G ( x, t + h,τ )dτ − ∫ G ( x, t ,τ )dτ ⎤ + lim ∫ G ( x, t + h,τ )dτ ⎢ 0 0 h →0 h ⎣ ⎦⎥ h→0 h t

1 1 t +h G ( x, t + h,τ ) − G ( x, t ,τ ) ] dτ + lim ∫ G ( x, t + h,τ )dτ [ 0 h →0 h h →0 h t t

= ∫ lim

∂G 1 t +h x, t ,τ ) dτ + lim ∫ G ( x, t + h,τ )dτ . ( 0 ∂t h →0 h t

=∫

t

Para la segunda integral usamos el teorema del valor medio para integrales

1 t +h 1 ( , , ) lim G x t + h τ d τ = G ( x, t , ω ) h = G ( x, t , t ) , h →0 h ∫t h →0 h

lim

en esta última igualdad como t ≤ ω ≤ t + h entonces cuando h → 0 tendremos

ω →t.

Entonces tendremos la siguiente expresión t ∂G ∂u ( x, t ) = ∫ ( x, t ,τ )dτ + G ( x, t , t ) . 0 ∂t ∂t

Tomando t = 0 , tenemos

∂u ( x,0) = 0 . ∂t

Por lo tanto u verifica las condiciones iniciales (II.4.1.2). Ahora mostremos que la ecuación integral (II.4.2.1) verifica la ecuación en derivada parciales (II.4.1.1). En efecto, observemos que -56-

G ( x, t , t ) = 0 entonces t ∂G ∂u ( x, t ) = ∫ ( x, t ,τ )dτ , 0 ∂t ∂t

análogamente encontramos t ∂G ∂u ( x, t ) = ∫ ( x, t ,τ )dτ . 0 ∂x ∂x

A estas dos funciones halladas le volvemos aplicar la definición de derivada parcial con respecto a t y x respectivamente, obteniendo 2 t∂ G ∂ 2u = ( x , t ) ∫0 ∂t 2 ( x, t ,τ )dτ , ∂t 2 2 t∂ G ∂ 2u ( x, t ) = ∫ ( x, t ,τ )dτ . 0 ∂x 2 ∂x 2

Hallemos

∂ 2G ∂ 2G ( x, t ,τ ) y ( x, t ,τ ) . ∂t 2 ∂x 2

i

∂G 1 ( x, t ,τ ) = lim [G ( x, t + h,τ ) − G ( x, t ,τ ) ] h →0 h ∂t x + ( t −τ ) 1 x +( t + h−τ ) = lim ⎡ ∫ F (σ ,τ )dσ − ∫ F (σ ,τ )dσ ⎤ ⎥⎦ τ τ x ( t h ) x ( t ) − + − − − h →0 h ⎢ ⎣

Gráficamente si h > 0 , tenemos:

entonces x + ( t + h −τ ) ∂G 1 x −(t −τ ) ( x, t ,τ ) = lim ⎡ ∫ F (σ ,τ )dσ + ∫ F (σ ,τ )dσ ⎤ ⎢ ⎥⎦ x + ( t −τ ) h→0 h ⎣ x −( t + h −τ ) ∂t x + ( t + h −τ ) 1 x −( t −τ ) ( , ) lim F σ τ d σ + F (σ ,τ )dσ . h→0 h ∫x −( t + h −τ ) h→0 ∫x + ( t −τ )

= lim

Nuevamente a cada uno de estos términos le aplicamos el teorema del valor medio para integrales, obteniendo -57-

∂G ( x, t ,τ ) = F ( x − t + τ ,τ ) + F ( x + t − τ ,τ ) ∂t entonces

∂ 2G ∂F ∂F ( x, t ,τ ) = ( x − t + τ ,τ ) + ( x + t − τ ,τ ) . 2 ∂t ∂t ∂t

i

∂G 1 ( x, t ,τ ) = lim [G ( x + h, t ,τ ) − G ( x, t ,τ ) ] h →0 h ∂x x + ( t −τ ) 1 x + h+( t −τ ) = lim ⎡ ∫ F (σ ,τ )dσ − ∫ F (σ ,τ )dσ ⎤ . ⎥⎦ τ τ x h ( t ) x ( t ) + − − − − h →0 h ⎢ ⎣

Gráficamente h > 0 , tenemos:

Similarmente a lo que hicimos para

∂G ( x, t ,τ ) , llegamos a ∂t

∂G ( x, t ,τ ) = F ( x + t − τ ,τ ) − F ( x − t + τ ,τ ) ∂x entonces

∂ 2G ∂F ∂F = + − − ( x , t , τ ) ( x t τ , τ ) ( x − t + τ ,τ ) . ∂x 2 ∂x ∂x Luego reemplazamos las expresiones obtenidas en: 2 2 t∂ G t∂ G ∂ 2u ∂ 2u ( x, t ) − 2 ( x, t ) = ∫ ( x, t ,τ )dτ − ∫ ( x, t ,τ )dτ 0 ∂t 2 0 ∂x 2 ∂t 2 ∂x t ⎡ ∂F t ⎡ ∂F ∂F ∂F ⎤ ⎤ = ∫ ⎢ ( x − t + τ ,τ ) + ( x + t − τ ,τ ) ⎥ dτ − ∫ ⎢ ( x + t − τ ,τ ) − ( x − t + τ ,τ ) ⎥ dτ 0 ∂t 0 ∂t ∂t ∂t ⎣ ⎦ ⎣ ⎦

= F ( x, t ) + F ( x, t ) t⎧ 1 1 ⎫ − ∫ ⎨lim ⎡⎣ F ( ( x + h) + t − τ ,τ ) − F ( x + t − τ ,τ ) ⎤⎦ − lim ⎡⎣ F ( ( x + h) − t + τ ,τ ) − F ( x − t + τ ,τ ) ⎤⎦ ⎬ dτ 0 h →0 h h →0 h ⎩ ⎭

= 2 F ( x, t ) t⎧ 1 1 ⎫ − ∫ ⎨lim ⎡⎣ F ( x + (t + h) − τ ,τ ) − F ( x + t − τ ,τ ) ⎤⎦ − lim ⎡⎣ F ( x − t + (τ + h),τ ) − F ( x − t + τ ,τ ) ⎤⎦ ⎬ dτ 0 h →0 h h →0 h ⎩ ⎭

-58-

t ⎡ ∂F ∂F ⎤ = 2 F ( x, t ) − ∫ ⎢ ( x + t − τ ,τ ) − ( x − t + τ ,τ ) ⎥ dτ 0 ∂t ∂τ ⎣ ⎦

= 2 F ( x , t ) − F ( x , t ) + F ( x , t ) = 2 F ( x, t ) . Recordemos que ignoramos la constante tenemos

1 que acompaña a la ecuación integral, por lo tanto 2

∂ 2u ∂ 2u ( x, t ) − 2 ( x, t ) = F ( x, t ) . ■ ∂t 2 ∂x Notemos, que en (II.4.2.1), el dominio de integración es el triángulo:

Dx ,t = {(σ ,τ ) ∈ ×

+

: x − σ ≤ t −τ} .

Vamos pues a considerar el problema de existencia y unicidad de soluciones de la ecuación integral (II.4.2.1), es decir de

u ( x, t ) = en

×

+

x + ( t −τ ) 1 t d τ [ f (σ ,τ ) + g (σ ,τ )u (σ ,τ )] dσ 2 ∫0 ∫x −( t −τ )

. Para eso es suficiente probar que, en todo dominio D ⊂

×

+

condición Dx ,t ⊂ D , para todo ( x, t ) ∈ D existe y es única la solución de (II.4.2.1).

-59-

con la

Consideremos el espacio de Banach X = C ( D, ) , y sea u ∈ X ; definimos Tu por

1

x + ( t −τ )

(Tu )( x, t ) = ∫0 dτ ∫x−(t −τ ) [ f (σ ,τ ) + g (σ ,τ )u (σ ,τ )] dσ . 2 t

Afirmación: Tu ∈ X . En efecto, análogamente a lo usado en (II.4.2.A) tenemos que

1

(Tu ) ( x, t ) = ∫0 G ( x, t ,τ )dτ 2 t

donde

G ( x, t ,τ ) = ∫

x + ( t −τ )

x −( t −τ )

Nuevamente volvemos a ignorar la constante

F (σ ,τ )dσ .

1 . Consideremos ( x, t ) , ( x0 , t0 ) ∈ D tal que 2

( x, t ) → ( x0 , t0 ) entonces

( x − x0 ) → 0 y ( t − t0 ) → 0 . Luego, para todo u ∈ X tenemos

(Tu ) ( x, t ) − (Tu ) ( x0 , t0 ) = ∫0 G ( x, t ,τ )dτ − ∫0 G( x0 , t0 ,τ )dτ t

t0

(II.4.2.B)

Además supongamos que t > t0 , con esto la primera integral de (II.4.2.B) se transforma en

∫

t

0

t0

t

0

t0

G ( x, t ,τ )dτ = ∫ G ( x, t0 ,τ )dτ + ∫ G ( x, t0 ,τ )dτ

(II.4.2.C)

Reemplazamos (II.4.2.C) en (II.4.2.B) y agrupando el primer término de (II.4.2.C) con el segundo término de (II.4.2.B), obtenemos

(

(Tu ) ( x, t ) − (Tu ) ( x0 , t0 ) ≤ ∫0 G ( x, t0 ,τ )dτ − ∫0 G ( x, t0 ,τ )dτ ≤∫

t0

0

t0

t0

( G ( x, t0 ,τ ) − G ( x0 , t0 ,τ ) ) dτ + ∫t

t

)

t

+ ∫ G ( x, t ,τ )dτ

G ( x, t ,τ ) dτ .

t0

(II.4.2.D)

0

Ahora calculemos el límite de cada una de las integrales de (II.4.2.D):

i Para la segunda integral, observemos que se está integrando sobre el intervalo [t0 , t ] el cual tiene medida nula, pues ( t − t0 ) → 0 , entonces ella misma tiende a cero.

-60-

i Para la primera integral, hacemos a = t0 − τ y explicitamos G ( x, t0 ,τ ) − G ( x0 , t0 ,τ ) = ∫

x+a

x−a

F (σ ,τ ) dσ − ∫

x0 + a

x0 − a

F (σ ,τ ) dσ .

(II.4.2.E)

Aquí supondremos que x > x0 de donde x + a > x0 + a y x − a > x0 − a . Gráficamente tenemos:

Luego (II.4.2.E) se transforma en

G ( x, t0 ,τ ) − G ( x0 , t0 ,τ ) = ∫

x+a

x0 + a

F (σ ,τ ) dτ − ∫

x −a

x0 − a

F (σ ,τ ) dτ

De forma similar a como tratamos la segunda integral de (II.4.2.D), pero ahora teniendo presente ( x − x0 ) → 0 , concluimos que cada una de las integrales anteriores tiende a cero. Por lo tanto, de (II.4.2.B) para todo u ∈ X tenemos

(Tu ) ( x, t ) − (Tu ) ( x0 , t0 ) → 0 entonces

(Tu ) ( x, t ) → (Tu ) ( x0 , t0 ) . i.e. Tu es continua en D , luego Tu ∈ X . También obtenemos

1

x + ( t −τ )

(Tu )( x, t ) − (Tv )( x, t ) ≤ ∫0 ∫x−(t −τ ) 2 t

g (σ ,τ ) u (σ ,τ ) − v(σ ,τ ) dσ dτ .

Si g (σ ,τ ) ≤ M en el dominio D , entonces

(Tu )( x, t ) − (Tu )( x, t ) ≤

t x + ( t −τ ) 1 M∫ ∫ u (σ ,τ ) − v(σ ,τ ) dσ dτ 0 x −( t −τ ) 2

t x + ( t −τ ) 1 ≤ M sup u (σ ,τ ) − v(σ ,τ ) ∫0 ∫x −( t −τ ) dσ dτ 2 D

1 d (Tu , Tv ) ≤ Mt 2 d ( u , v ) ; u , v ∈ X . 2 Análogamente conseguimos

-61-

(T u ) ( x, t ) − (T v ) ( x, t ) = T (Tu) ( x, t ) − T (Tv )( x, t ) 2

2

≤

1 2

≤

1 t x +(t −τ ) 1 M M τ 2 d ( u , v ) dσ dτ ∫ ∫ 0 ( ) x t − − τ 2 2

t

x + ( t −τ )

0

x −( t −τ )

∫∫

g (σ ,τ ) [ (Tu )(σ ,τ ) − (Tv)(σ ,τ ) ] dσ dτ

t x + ( t −τ ) 1 = M 2d ( u, v ) ∫ ∫ τ 2 dσ dτ τ 0 ( ) x t − − 4

d ( T 2u , T 2 v ) ≤

1 24 M t d ( u, v ) ; u, v ∈ X . 24

Por inducción, resulta

d (T nu , T nv) ≤

1 M n (t * ) 2 n d ( u , v ) (2n)!

donde

t * = sup {( x, t ) ∈ D} . Para cada n ∈

+

, sea

cn =

M (t * ) 2 (2n)!

n

,

ahora utilizamos el criterio de la razón para sucesiones de números reales

M (t * ) 2 lim n →∞

n +1

cn+1 [ 2(n + 1)]! = M (t * )2 lim 1 = 0 < 1 , = lim n n →∞ n→∞ 2( n + 1) cn M (t * ) 2 (2n)!

entonces

lim n →∞

M (t * ) 2

n

n!

= 0,

⎧ M (t * ) 2 ⎫ ⎬ converge, luego sigue por definición de límite de sucesiones de n ! ⎩ ⎭n≥1

por lo tanto ⎨

números reales que: Dado 0 < ε < 1 , existe n0 ∈

tal que para todo n > n0 tenemos -62-

M (t * ) 2 (2n)!

n

< ε < 1.

Entonces tomamos m = n0 , para concluir que existe m ≥ 1 tal que T

C ([a, b], ×

+

)= X .

m

es una contracción en

Así, por el corolario I.4.2. demostramos la existencia y unicidad en

de la solución de (II.4.2.1). ■

II.4.3. EJEMPLO Consideremos ahora la ecuación diferencial parcial

u xy = f ( x, y, u , u x , u y ) . Sobre una curva C , dada por y = ϕ ( x) con

(II.4.3.1)

ϕ ∈ C (1) ( ,

)

y

ϕ ′( x) < 0 , x ∈ . Sean

dadas las funciones continuamente diferenciables u0 y u1 ( ó u2 ). Queremos estudiar la existencia de u ∈ C

(2)

(

) satisfaciendo (II.4.3.1) y

× ,

( ux )|C = u1 ,

u|C = u0 , esto es, tal que, para todo x ∈

[ ó ( u y ) |C = u2 ],

, tenemos

⎛ ∂u ⎞ u ( x,ϕ ( x)) = u0 ( x) , ⎜ ⎟ ( x,ϕ ( x)) = u1 ( x) ⎝ ∂x ⎠

⎛ ∂u ⎞ [ ó ⎜ ⎟ ( x,ϕ ( x)) = u2 ( x) ]. ⎝ ∂y ⎠

Vamos a mostrar inicialmente como podemos reducir el problema a la búsqueda de una solución u de una ecuación análoga a (II.4.3.1) con

u ≡ u x ≡ u y ≡ 0 sobre la recta x + y = 0 . Observemos preliminarmente que, dadas u0 , u1 [respectivamente u2 ], queda determinada u2

[respectivamente u1 ], pues; si u|C = u0 , ( u x ) |C = u1 y si ( u y ) |C = u2 , entonces d u0 ( x,ϕ ( x)) = u x ( x,ϕ ( x)) + u y ( x,ϕ ( x))ϕ ′( x) = u1 ( x,ϕ ( x)) + ϕ ′( x)u2 ( x,ϕ ( x)) dx lo que determina u2 . Observar además que no hay pérdida de generalidad en suponer que la curva C es dada por

x+ y=0 pues siempre podemos hacer un cambio de variables poniendo

ξ = − y , η = ϕ ( x) . Además de eso, si tuviéramos

-63-

(II.4.3.2)

_

u [ x(t ), y (t ) ] = u (t ) , u x [ x(t ), y (t ) ] = p(t ) ,

(II.4.3.3)

u y [ x(t ), y (t ) ] = q(t ) ,

a lo largo de la recta x = t , y = −t , consideremos la ecuación

u xy = 0

(II.4.3.4)

sujeta a esas condiciones. La solución es entonces, _

u1 ( x, y ) =

_

u ( x) + u ( y ) 1 x 1 y + ∫ p (t )dt + ∫ q(t )dt . 2 2 −y 2 −x

Sustrayendo u1 de una solución de (II.4.3.1) que satisfaga también a (II.4.3.3), obtenemos una solución de una nueva ecuación (cuyo segundo miembro, para simplificar, también denotamos por f ).

u xy = f ( x, y, u , u x , u y ) i.e.

(II.4.3.1)

u ≡ u x ≡ u y ≡ 0 sobre la recta x + y = 0 .

(II.4.3.5)

satisfaciendo

El problema de existencia de una solución de los sistemas (II.4.3.1) y (II.4.3.5), que pasamos a considerar. II.4.4. PROPOSICIÓN Supongamos que dado a > 0 , la función f : [ − a, a ] × [ − a, a ] ×

3

→

es continua

y lipschitziana en las tres últimas variables, esto es, existe una constante L > 0 tal que

f ( x, y, u2 , p2 , q2 ) − f ( x, y, u1 , p1 , q1 ) ≤ L ⎣⎡ u2 − u1 + p2 − p1 + q2 − q1 ⎤⎦ para cualesquiera x, y ∈ [ − a, a ] con u1 , u2 , p1 , p2 , q1 , q2 ∈

. Entonces u es solución de

(II.4.3.1) sujeta a las condiciones (II.4.3.5) si y sólo si es solución de

u ( x, y ) = ∫∫ f (ξ ,η , u , uξ , uη )dξ dη D

donde D es el triangulo limitado por las rectas

ξ +η = 0 , ξ = x , η = y , Gráficamente tenemos:

-64-

(II.4.4.1)

con lo que: y

x

−x

−η

x

y

−y

−ξ

u ( x, y ) = ∫ dη ∫ f (ξ ,η , u, uξ , uη )dξ

(II.4.4.2)

= ∫ dξ ∫ f (ξ ,η , u , uξ , uη )dη . Demostración: También en este caso vamos aplicar el corolario 1.4.2. para demostrar la existencia de soluciones de (II.4.3.1) satisfaciendo (II.4.3.5). Vamos a considerar dos casos, expuestos a seguir: Primer caso: Supongamos además que f es lineal en las tres últimas variables i.e. cuando (II.4.3.1) tiene la forma

u xy = g ( x, y )u x + h( x, y )u y + r ( x, y )u + s ( x, y ) ×

donde las funciones g , h, r , s son continuas en Vamos a demostrar que existe u :

2

→

.

solución de (II.4.4.2) y por lo tanto, solución

de (II.4.3.1) satisfaciendo (II.4.3.5). Para eso, es suficiente demostrar que, para todo a > 0 existe una solución de (II.4.2.2) en [ − a, a ] × [ − a, a ] . Sea X = C

(1)

([ −a, a ] × [ −a, a ], ) i.e. X

es el espacio de las funciones continuas

con primeras derivadas parciales continuas. Recordemos que X es un espacio de Banach provisto de la norma

u donde,

X

= u + ux + u y

v = sup { v( x, y ) : x, y ∈ [ − a, a ]} . -65-

Para todo u ∈ X y para cualesquiera x, y ∈ [ − a, a ] , definimos Tu por

(Tu ) ( x, y) = ∫∫ f (ξ ,η , u, uξ , uη )dξ dη .

(II.4.4.3)

(Tu ) x ( x, y ) = ∫− x f ( x,η , u, ux , uη )dη ,

(II.4.4.4)

(Tu ) y ( x, y ) = ∫− y f (ξ , y, u, uξ , uη )dξ .

(II.4.4.5)

D

Luego se cumplen y

x

En efecto, mostraremos (II.4.4.4). De (II.4.4.3) y de la segunda igualdad de (II.4.4.2) tenemos

(Tu ) ( x, y) = ∫− y dξ ∫−ξ f (ξ ,η , u, uξ , uη )dη . x

Sea

y

G (ξ , y, u , uξ , u y ) = ∫ f (ξ ,η , u , uξ , uη )dη . y

−ξ

con esto

(Tu ) ( x, y) = ∫− y G (ξ , y, u, uξ , u y ) dξ . x

Entonces

(Tu ) x ( x, y ) = lim h →0

(Tu ) ( x + h, y ) − (Tu ) ( x, y ) h

x 1 x+h = lim ⎡ ∫ G (ξ , y, u , uξ , u y ) dξ − ∫ G (ξ , y, u , uξ , u y ) dξ ⎤ ⎥⎦ −y h →0 h ⎢ ⎣ −y

1 x+h G (ξ , y, u , uξ , u y ) dξ h→0 h ∫x

= lim

1 = lim G ( w, y, u, uw , u y ) h = G ( x, y, u, u x , u y ) h →0 h para algún w tal que x ≤ w ≤ x + h ; usando el teorema del valor medio para integrales como h → 0 entonces w → x .

(

)

i.e. (Tu ) x ( x, y ) = G x, y, u , u x , u y , por lo tanto

(Tu ) x ( x, y ) = ∫− x f ( x,η , u, ux , uη )dη . y

Luego sigue de (II.4.4.3), (II.4.4.4) y (II.4.4.5) que Tu ∈ X , y para cualesquiera u , v ∈ X , mayoremos: -66-

i (Tu − Tv ) ( x, y ) ≤

∫

y

−x

dη ∫

x

f (ξ ,η , u , uξ , uη ) − f (ξ ,η , v, vξ , vη ) dξ

−η

≤ L ∫ dη ∫ ⎡ u − v + uξ − vξ + uη − vη ⎤ dξ ⎦ −x −η ⎣ y

x

= L u−v

X

∫

y

x

−x

dη ∫ dξ −η

i.e. (Tu − Tv ) ( x, y ) ≤ L u − v

i

( (Tu ) x − (Tv) x ) ( x, y) ≤ ∫− x ∫−η y

x

X

x+ y

2

(II.4.4.6)

f (ξ ,η , u, u x , uη ) − f (ξ ,η , v, vx , vη ) dξ dη

≤ L ∫ ⎡ u − v + u x − vx + uη − vη ⎤ dη ⎦ −x ⎣ y

= L u −v

X

∫

y

−x

dη = L u − v

i.e. ( (Tu ) x − (Tv) x ) ( x, y ) ≤ L u − v

X

X

x+ y

x+ y

(II.4.4.7)

x+ y

(II.4.4.8)

i Análogamente, conseguimos

( (Tu)

y

− (Tv) y ) ( x, y ) ≤ L u − v

X

Para cualesquiera x, y ∈ [ − a, a ] (i.e. x ≤ a , y ≤ b ) tenemos

x + y ≤ x + y ≤ 2a , luego tomando c = max {2a,2} , tendremos

( x + y)

2

2

= x + y = x + y x + y ≤ 2a x + y ≤ c x + y i.e. ( x + y ) ≤ c x + y . 2

y, por tanto, en (II.4.4.6), (II.4.4.7) y (II.4.4.8) podemos tomar como segundo miembro a

cL u − v

X

x+ y .

Entonces tenemos

Tu − Tv

X

= Tu − Tv + (Tu − Tv ) x + (Tu − Tv ) y

≤ 3Lc u − v

X

x+ y .

Similarmente, encontraremos mayoraciones para: -67-

(II.4.4.9)

i (T 2u − T 2v ) ( x, y ) = T (Tu )( x, y ) − T (Tv)( x, y ) ≤

∫

y

−x

dη ∫

x

f (ξ ,η , Tu,(Tu )ξ ,(Tu )η ) − f (ξ ,η , Tv,(Tv)ξ ,(Tv)η ) dξ

−η

≤ L ∫ dη ∫ ⎡ Tu − Tv + (Tu )ξ − (Tv)ξ + (Tu )η − (Tv)η ⎤ dξ ⎦ −x −η ⎣ y

x

≤ 3Lc u − v

∫

y

−x

X

x

dη ∫ ξ + η dξ

(II.4.4.10)

(Aquí utilizamos (II.4.4.9))

−η

Análogamente, llegamos a conseguir las siguientes relaciones:

i

( (T u)

i

( (T u)

x

− (T 2v) x ) ( x, y ) ≤ 3L2c u − v

y

− (T 2v) y ) ( x, y ) ≤ 3L2c u − v

2

2

∫

y

∫

y

−x

X

X

−x

x + η dη

(II.4.4.11)

ξ + y dξ

(II.4.4.12)

Observemos por ejemplo, como obtuvimos

∫

y

−x

⎛ η2 ⎞ x + η dη = ∫ ( x + η ) dη = ⎜ xη + ⎟ −x 2 ⎠ ⎝ y

y2 ⎡ x2 ⎤ = xy + − x (− x) + ⎥ 2 ⎢⎣ 2⎦

y

−x

y 2 x 2 x 2 + 2 xy + y 2 ( x + y ) = xy + + = = 2 2 2 2

2

donde hemos partido del hecho que:

− x < η < y luego x − x < x + η < x + y i.e. 0 < x + η < x + y entonces x + η = x + η . Entonces conseguimos

T 2u − T 2 v

X

= T 2 u − T 2 v + ( T 2u − T 2 v ) + ( T 2u − T 2 v ) ≤ x

y

9 L2c 2 u−v 3!

X

siendo los segundos miembros de (II.4.4.10), (II.4.4.11) y (II.4.4.12) mayorados por

3L2c 2 u − v Por inducción viene que

T u −T v n

Para cada n ∈

+

n

X

( 3Lc ) ≤

2

X

x+ y .

n

n!

u−v

, sea n

3Lc , bn = n! -68-

X

; u, v ∈ X .

x+ y

2

ahora utilizamos el criterio de la razón para sucesiones de números reales n +1

lim n →∞

bn+1 bn

3Lc 1 (n + 1)! = lim = 3Lc lim = 0 m0 tenemos < ε < 1. m!

Entonces tomamos n = m0 para concluir que existe n ≥ 1 tal que T es una contracción en n

C (1) ([ −a, a ] × [ − a, a ] ,

) = X , luego por el corolario I.4.2. T

tiene un único punto fijo en

X , es decir existe una única función u ∈ X la cual satisface (II.4.4.2). Segundo Caso: Supongamos ahora que f esta definida en una vecindad de la recta

x + y = 0 y que la constante de Lipschitz de f (en las tres últimas variables u, p y q ) crezca con u , p , q ; entonces podemos demostrar la existencia y unicidad de la solución de (II.4.2.2) sólo en una vecindad de la recta x + y = 0 . Precisamente nos vamos a restringir a

A = [ −α ,α ] × [ −α ,α ] × Br [ 0] × BR [ 0] × BR [ 0] , donde Br [ 0] , BR [ 0] son bolas cerradas centradas en 0 y de radio r y R respectivamente, además

α es escogido de tal manera que, si u : [ −α ,α ] × [ −α ,α ] →

toma valores en

Br [ 0] y u x , u y en BR [ 0] , entonces lo mismo vale para Tu definida por (II.4.4.3). Es posible escoger

α así, pues si f = M en A ; entonces de (II.4.4.3), (II.4.4.4), y

(II.4.4.5) se sigue que:

-69-

i (Tu ) ( x, y ) =

∫∫ f (ξ ,η , u, uξ , uη )dξ dη

≤M(x + y)

D

≤ ∫∫ f (ξ ,η , u , uξ , uη ) dξ dη ≤ f D

i (Tu ) x ( x, y ) ≤ ∫

y

−x

∫∫ dξ dη = 4M α

2

.

D

f (ξ ,η , u, uξ , uη )dη ≤ M x + y

≤ M ( x + y ) ≤ M (α + α ) = 2 M α . i (Tu ) y ( x, y ) ≤ 2α M . Sea X = C

1

([ −α ,α ] × [ −α ,α ], B [0]) cuyas derivadas asumen valores en B [0] , r

R

los cálculos hechos en el primer caso se aplican. ■

-70-

II.5.- APLICACIÓN AL PROBLEMA DE STURM-LIOUVILLE El problema de Sturm-Liouville consiste en resolver la ecuación diferencial con valores en la frontera

⎧ d 2u ⎪ 2 + λ f ( x, u ) = 0 ⎨ dx ⎪u (0) = u (1) = 0 ⎩

(II.5.1)

la cual presenta una conexión con el problema de una varilla elástica bajo compresión con extremos fijos; aquí u denota el desplazamiento axial y

λ la carga axial. Para valores peque-

λ , la solución es u = 0 (a menos que λ corresponda a un autovalor), el cual corres-

ños de

ponde a la varilla recta. Pero como la carga

λ es aumentada, la varilla será doblada; dando

origen a la solución no nula. Queremos encontrar una ecuación integral equivalente al problema de Sturm-Liouville (II.5.1). Sea u con segunda derivada continua y definamos

∫ ∫ ϕ ( s )dsdy = ∫ ∫ x

0

y

x

0

0

y

0

x

y

0

0

ϕ := −u′′ . Ahora

( −u′′( s) )dsdy = ∫ ( −u′(s ) ) dy

x

x

0

0

= ∫ ( −u′( y ) + u′(0) )dy = ( −u ( y ) + u′(0) y ) = −u ( x) + u′(0) x + u (0) 0

∫ ∫ ϕ ( s )dsdy = −u ( x) + u′(0) x x

i.e.

y

0

0

despejando u ( x) tenemos

u ( x) = xu′(0) − ∫

x

0

∫ ϕ ( s )dsdy . y

0

Integramos por partes la integral doble

u = ∫ ϕ ( s )ds y

dv = dy

0

du = ϕ ( y ) x

∫∫ 0

y

0

(

ϕ ( s )dsdy = y ∫ ϕ ( s )ds y

0

)

v= y x

− ∫ yϕ ( y )dy = x ∫ ϕ ( s )ds − ∫ yϕ ( y )dy

0

x

x

x

0

0

0

= x ∫ ϕ ( y )dy − ∫ yϕ ( y )dy = ∫ ( x − y )ϕ ( y )dy x

x

x

0

0

0

i.e.

x

∫∫ 0

y

0

ϕ ( s )dsdy = ∫ ( x − y )ϕ ( y )dy x

0

-71-

Por lo que

u ( x) = xu′(0) − ∫ ( x − y )ϕ ( y )dy . x

(II.5.2)

0

Tenemos 1 1

∫∫ x

y

ϕ ( s )dsdy = ∫

1 1

1

1

1

( −u′′( s) )dsdy = ∫x ( −u′( s) ) y dy = ∫x ( −u′(1) + u′( y) ) dy x ∫y 1

= ( −u′(1) y + u ( y ) ) = − u′(1) + u (1) − ( −u′(1) x + u ( x) ) x

0

= −u ( x) − u′(1)(1 − x) , entonces sigue

u ( x) = −u′(1)(1 − x) − ∫

1 1

x

∫ ϕ ( s )dsdy . y

Integraremos por partes la integral doble

m = ∫ ϕ ( s )ds y

dn = dy

1

dm = ϕ ( y ) 1 1

∫∫ x

y

1

(

n= y

)

ϕ ( s )dsdy = ∫ − ∫ ϕ ( s )ds dy x

= −∫

1

y

1

∫

y

x 1

⎡

(

ϕ ( s )dsdy = − ⎢ y ∫ ϕ ( s ) ds ⎣

y

1

)

1 ⎤ − ∫ yϕ ( y ) dy ⎥ x x ⎦

1

x 1 1 1 = − ⎡ − x ∫ ϕ ( s ) ds − ∫ yϕ ( y ) dy ⎤ = − ⎡ x ∫ ϕ ( s ) ds − ∫ yϕ ( y ) dy ⎤ ⎢⎣ 1 ⎥⎦ ⎢⎣ x ⎥⎦ x x

1 1 1 = − ⎡ x ∫ ϕ ( y ) dy − ∫ yϕ ( y ) dy ⎤ = − ∫ ( x − y ) ϕ ( y ) dy x x ⎣⎢ x ⎦⎥

i.e.

1 1

1

x

x

∫∫

y

ϕ ( s )dsdy = − ∫ ( x − y )ϕ ( y ) dy

entonces 1

u ( x) = −u′(1)(1 − x) + ∫ ( x − y )ϕ ( y )dy .

(II.5.3)

x

De igual forma 1

1

1

0

0

0

∫

ϕ ( y )dy + u′(1) − u′(0) = ∫ ( −u′′( y ) )dy + u′(1) − u′(0) = ( −u′( y ) ) + u′(1) − u′(0) = −u′(1) + u′(0) + u′(1) − u′(0) = 0 -72-

1

∫ ϕ ( y )dy + u′(1) − u′(0) = 0

i.e.

(II.5.4)

0

para todo x ∈ [ 0,1] .Multiplicando (II.5.2) por (1 − x) , (II.5.3) por x y (II.5.4) por (1 − x ) x y entonces sumando estas tres nuevas ecuaciones y agrupando adecuadamente obtenemos x

1

0

x

u ( x) = ∫ (1 − x ) yϕ ( y )dy + ∫ x (1 − y )ϕ ( y ) dy

(II.5.5)

1

u ( x) = ∫ G ( x, y )ϕ ( y )dy 0

donde G es la función de Green

⎧ x(1 − y ), x ≤ y ⎪ G ( x, y ) = ⎨ ⎪ y (1 − x), y ≤ x ⎩ Pero como

ϕ ( x) = −u′′( x) = λ f ( x, u ( x) ) tenemos 1

u ( x) = λ ∫ G ( x, y ) f ( y, u ( y ) )dy

(II.5.6)

0

2

con G definida como antes (Observemos que G es continua en

).

Recíprocamente, sea ϕ ∈ C [ 0,1] solución de la ecuación integral (II.5.5). Para todo x ∈ [ 0,1] definamos la función x

1

0

x

u ( x) = ∫ (1 − x) yϕ ( y )dy + ∫ x(1 − y )ϕ ( y )dy , entonces u (0) = 0 y u (1) = 0 y por la construcción de la ecuación integral tenemos

−u′′ = ϕ = λ f ( x, u ) . Por consiguiente, el problema con valor en la frontera (II.5.1) y la ecuación integral (II.5.5) son equivalentes. Es decir, resolver el problema de Sturm-Liouville es equivalente a resolver la ecuación integral (II.5.6). Observemos que esta ecuación integral equivalente es la ecuación de Fredholm de primera especie (ver II.2.1.2, página 38) donde a = 0 , b = 1 y el núcleo para cualesquiera x, y ∈ [ 0,1] es K ( x, y ) = G ( x, y ) . Consideremos el espacio de las funciones cuadrado integrable en el dominio [ 0,1] . 2

i.e. L

([0,1]) = { f : [0,1] →

; f medible en [ 0,1] y -73-

∫

1

0

f ( x) dx < ∞ .} 2

2

Recordemos que L

([0,1]) es un espacio de Banach.

Definimos f : [ 0,1] → *

por

f*= f ψ donde ψ : [ 0,1] →

2

con ψ = ( id , u ) , id denota la función identidad y f :

2

→

es

la función dada i.e. tenemos para cada x ∈ [ 0,1]

f * ( x) = ( f ψ ) ( x) = f (ψ ( x)) = f ( x, u ( x) ) . Ahora definamos el operador

T : L2 ([ 0,1]) → L2 ([ 0,1]) u Tu dado por 1

(Tu ) ( x) = λ ∫0 G ( x, y ) f * ( y )dy . Es evidente que un punto fijo de T es solución de la ecuación integral y recíprocamente. 2

Para asegurar que T es una contracción en L

([0,1]) , asumimos que

∂f * ∂f * es acotado i.e. ( x) ≤ M ; 0 ≤ x ≤ 1 , u ∈ L2 ([ 0,1]) . ∂u ∂u Denotemos a

K ( x, y , u ) = G ( x, y ) f * ( y ) ,

(II.5.A)

entonces por el teorema del valor medio, tenemos

K ( x, y , u1 ) − K ( x, y, u2 ) = =

∂K ( x, y, u ) u1 ( y ) − u2 ( y ) ∂u ∂ G ( x, y ) f * ( y ) ) u1 ( y ) − u2 ( y ) ( ∂u

= G ( x, y )

∂f * ( y ) u1 ( y ) − u2 ( y ) ∂u

≤ M G ( x, y ) u1 ( y ) − u2 ( y ) i.e. K ( x, y, u1 ) − K ( x, y, u2 ) ≤ M G ( x, y ) u1 ( y ) − u2 ( y ) para algún u entre u1 y u2 . -74-

(II.5.B)

Sea

μ2 = M 2∫

1 1

∫

0 0

2

G ( x, y ) dxdy .

(II.5.C)

Probaremos que

M2 μ = M ∫ ∫ G ( x, y ) dxdy = . 0 0 90 2

1 1

2

2

En efecto, teniendo en cuenta la función de Green con la que estamos trabajando, calculemos 1 1

∫∫

0 0

1 1

( G ( x, y ) ) 0 ∫0

G ( x, y ) dxdy = ∫ 2

2

1 y 1 2 2 dxdy = ∫ ⎡ ∫ ( G ( x, y ) ) dx + ∫ ( G ( x, y ) ) dx ⎤ dy ⎥⎦ 0⎢ y ⎣ 0

) (∫ = ∫ ( (1 − y ) ∫ x dx + y ∫ (1 − x) dx ) dy

=∫

1

y

0

0

1

x 2 (1 − y ) 2 dx + ∫ y 2 (1 − x) 2 dx dy y

1

2

0

y

2

2

0

⎛ x3 = ∫ ⎜ (1 − y ) 2 0⎜ 3 ⎝ 1

1

2

y

x= y

(1 − x)3 −y 3 2

x =0

(

⎞ ⎟ dy ⎟ x= y ⎠

x =1

)

=

1 1 3 2 3 2 (1 − y ) y + y 1 − y dy ( ) 3 ∫0

=

1 1 4 y − 2 y 3 + y 2 ) dy ( ∫ 3 0 y =1

1 ⎛ y5 y 4 y3 ⎞ 1⎛ 1 1 1⎞ 1 . = ⎜ − + ⎟ = ⎜ − + ⎟= 3⎝ 5 2 3 ⎠ y =0 3 ⎝ 5 2 3 ⎠ 90 i.e.

1 1

∫∫

0 0

2

G ( x, y ) dxdy =

1 . 90 2

Ahora mostraremos que T es una contracción en L

(II.5.D)

([0,1]) para un valor fijo de λ .

En efecto, con (II.5.A) el operador T queda definido por 1

(Tu ) ( x) = λ ∫0 K ( x, y, u ( y ) )dy . Para cualesquiera u1 , u2 ∈ L

2

Tu1 − Tu2

([0,1]) tenemos 2

=

{∫ (Tu ) ( x) − (Tu ) ( x) dx} 1

0

2

1

2

-75-

1 2

.

(II.5.E)

Mayoraremos la función que estamos integrando i.e. 1

1

(Tu1 ) ( x) − (Tu2 ) ( x) = λ ∫0 K ( x, y, u1 ( y ) )dy − λ ∫0 K ( x, y, u1 ( y ) )dy ≤λ

≤λ

∫

K ( x, y, u1 ( y )) − K ( x, y, u2 ( y )) dy

0 1

∫ M G( x, y) u ( y) − u ( y) dy 1

0

≤λ⎡ ⎢⎣ =λ

1

(∫ M 1

2

2

G ( x, y ) dy

0

(∫ M 1

2

0

(usamos (II.5.B))

2

2

G ( x, y ) dy

)

)( ∫ u ( y) − u ( y) dy )⎤⎥⎦

1 2

1

2

1

0

1 2

2

u1 − u2

2

donde, el penúltimo paso lo conseguimos usando la desigualdad de Hölder para integrales. Entonces tenemos

(

1

(Tu1 ) ( x) − (Tu2 ) ( x) ≤ λ ∫0 M

2

2

G ( x, y ) dy

)

1 2

u1 − u2 2 .

Reemplazamos esta desigualdad en (II.5.E)

⎧ ⎪ 1 Tu1 − Tu2 2 ≤ ⎨ ∫ λ 0 ⎪⎩

0

= λ u1 − u2

2

= λ u1 − u2

2

i.e. Tu1 − Tu2 Ahora sea α =

(∫ M 1

2

2

2

G ( x, y ) dy

)

1 2

2

u1 − u2

2

( M ∫ ∫ G( x, y) dydx ) 1 1

2

2

⎫ ⎪ dx ⎬ ⎪⎭

1 2

1 2

0 0

μ

(usamos (II.5.C))

≤ λ μ u1 − u2 2 ; u1 , u2 ∈ L2 ([ 0,1]) .

μ λ > 0 . Notemos que α < 1 sólo si λ < 2

Así, concluimos que T es una contracción en L

1

μ

.

([0,1]) siempre que:

∂f * i). f ∈ L ([ 0,1]) y ( x) ≤ M , para cada 0 ≤ x ≤ 1 y para todo u ∈ L2 ([ 0,1]) . ∂u *

ii).

λ
3 10 .

-78-

II.6. APLICACIÓN A LA PROGRAMACIÓN DINÁMICA La optimización Dinámica como su nombre lo indica, estudia la optimización de sistemas dinámicos, es decir, sistemas que evolucionan en el tiempo. Dado un sistema, se trata de guiar o controlar el sistema de manera óptima mediante el manejo de otras magnitudes llamadas variables de control a lo largo de un horizonte temporal dado, de acuerdo a un objetivo previamente fijado. Consideremos por ejemplo, nuestro propio cuerpo como un sistema dinámico el cual el médico examina tomándonos la presión arterial, la temperatura, analizando los estudios realizados, como son un electrocardiograma, la cantidad de plaquetas, de leucocitos, de glóbulos rojos, el nivel de colesterol, el nivel de enzimas que intervienen en el metabolismo de los aminoácidos (transaminasas), el azúcar de la sangre (glucosa), etc. De esta manera el médico determina el estado en que está nuestro organismo. Si considera que el mismo no es el adecuado, entonces tomará medidas de control para mejorar dicho estado. Entre muchas prescripciones, nos puede recomendar hacer una mayor actividad física, ingerir ciertos medicamentos, seguir una dieta específica de comidas, etc. Llevadas a la práctica estas medidas tenderán a modificar el estado de nuestra salud. Los campos de aplicación de la optimización dinámica no se circunscriben solamente a ciertas ramas de la matemática, como el análisis matemático, la geometría o los sistemas dinámicos, sino que abordan otras áreas del conocimiento, como por ejemplo la ingeniería, la física, la economía, la biología, etc. De manera similar a la búsqueda de máximos y mínimos de una función que modeliza matemáticamente cierto problema, la finalidad de la optimización dinámica radica en determinar la existencia, y eventualmente el cálculo, de los valores de ciertas variables, llamadas variables de estado, consideradas como funciones de otras variables (por ejemplo, el tiempo y/o variables llamadas variables de control), que producen valores óptimos (máximos o mínimos, según sea el caso) de cierta cantidad de estudio llamada funcional objetivo durante cierto intervalo temporal. En cuanto a su tratamiento matemático, los problemas de optimización dinámica pueden resolverse por alguna de estas tres técnicas: 1.- El cálculo de variaciones. 2.- El control óptimo. 3.- La programación dinámica. -79-

Si bien las tres técnicas permiten abordar problemas en tiempo discreto y continuo, los métodos 1 y 2 se utilizan generalmente para el tiempo continuo y el restante para el tiempo discreto. II.6.1. DEFINICIÓN Sean X e Y dos conjuntos arbitrarios. Una correspondencia f de X en Y es una regla que asocia cada elemento x de X a un subconjunto f ( x) de Y . i.e. f : X → P( Y ) donde P( Y ) es el conjunto de partes de Y . Nuestro interés en este trabajo se centra en correspondencias del tipo

f : X → P( X ). Si usamos la notación de funciones, podemos definir una correspondencia X → como una función X

P( X )

{ Y ⊆ P( X )} .

Sean X e Y dos espacios métricos arbitrarios y sea una correspondencia f : X → P( Y ). El grafo de f es el conjunto G f =

∪ ( x, f ( x) ) donde

f ( x) es el conjunto imagen del

x∈X

elemento x (recordemos que f es una correspondencia y por lo tanto la imagen de cualquier punto es un conjunto). Se dice que el grafo es cerrado si el conjunto G f es un conjunto cerrado. El concepto de continuidad se generaliza para correspondencias a partir de los conceptos más débiles de semicontinuidad inferior y semicontinuidad superior. Para nosotros sólo es relevante el último de estos conceptos. II.6.2. DEFINICIÓN Si X ⊆

m

eY⊆

n

, la correspondencia f : X → P( Y ) se llama semicontinua

superior si su grafo es cerrado y para cada punto de X su imagen es un conjunto compacto. Por otro lado, la extensión del concepto de punto fijo a correspondencias es directa. II.6.3. DEFINICIÓN Dado un conjunto X , se dice que un elemento x ∈ X es un punto fijo de la correspondencia f : X → P( X ) si x ∈ f ( x ) . -80-

Obviamente para cada correspondencia f : X → P( X ), podemos definir

fix( f ) = { x ∈ X : x ∈ f ( x)} . Teniendo en mente esta definición y de la definición I.1. del capítulo I, podemos establecer que, tanto para correspondencias como para funciones, se verifica fix( f ) ⊆ X . Es claro que si definimos recursivamente para cualquier i ∈

y x∈ X ,

f (i +1) ( x ) = f ⎡⎣ f ( i ) ( x ) ⎤⎦ con f

(1)

≡ f , esto implica fix ( f ) ⊆ ∩ f (i ) ( X ) . i∈

En otras palabras, los puntos fijos de una función deben pertenecer a la intersección de las imágenes de sus sucesivas iteraciones. II.6.4. APLICACIÓN DEL T.P.F.B. A LA PROGRAMACIÓN DINÁMICA La Teoría de la Programación Dinámica fue elaborada principalmente por Richard Bellman durante la década de 1950, trata acerca de la toma óptima de decisiones en procesos de etapas múltiples, frecuentemente estocásticos. Dada una función continua

g:

n

×

n

→

y una correspondencia nunca vacía

f:

n

→ P(

n

)

semicontinua superiormente, un problema determinístico típico de esta teoría es el de hallar el valor

v ( x0 ) = max ∑ β t −1 g ( xt −1 , xt ) , { xt ∈ f ( xt −1 )}t∈ { xt }

t∈

donde la sucesión { xt }t∈ se elige en el espacio de sucesiones

→

n

y

β ∈ [ 0,1) es un

parámetro que frecuentemente se interpreta como un factor de descuento entre dos tiempos. Siguiendo con esta interpretación de t como un índice de tiempo, esta clase de problema se denomina de horizonte infinito. La estructura recursiva que presenta permite mostrar que la función valor v :

n

→

debe satisfacer la ecuación funcional

v ( xt −1 ) = max { g ( xt −1 , xt ) + β v( xt )} xt ∈ f ( xt −1 )

-81-

(II.6.4.1)

En efecto, en (II.6.4.1) damos valores enteros positivos para t , encontrando para

i t = 1 : v ( x0 ) = max { g ( x0 , x1 ) + β v( x1 )} x1∈ f ( x0 )

i t = 2 : v ( x1 ) = max { g ( x1 , x2 ) + β v( x2 )} x2∈ f ( x1 )

i t = 3 : v ( x2 ) = max { g ( x2 , x3 ) + β v( x3 )} x3∈ f ( x2 )

Reemplazando v ( x2 ) en v ( x1 ) , tenemos

{

}

v ( x1 ) = max g ( x1 , x2 ) + β max { g ( x2 , x3 ) + β v( x3 )} . x2∈ f ( x1 )

x3∈ f ( x2 )

Reemplazando v ( x1 ) en v ( x0 ) , tenemos

{

}

⎡ ⎤ v ( x0 ) = max ⎢ g ( x0 , x1 ) + β max g ( x1 , x2 ) + β v max ( g ( x2 , x3 ) + β v( x3 ) ) ⎥ x1∈ f ( x0 ) ⎣ x2∈ f ( x1 ) x3∈ f ( x2 ) ⎦ = max ⎡⎣ β 0 g ( x0 , x1 ) + β 1 g ( x1 , x2 ) + β 2 g ( x2 , x3 ) + β 3v( x3 ) ⎤⎦ x1∈ f ( x0 ) x2∈ f ( x1 ) x3∈ f ( x2 )

= max

∑β

{ xt }∈ f ( xt −1) t∈

t −1

g ( xt −1 , xt ) = v ( xt −1 )

t∈

(La importancia de la relación (II.6.4.1) es que nos permite maximizar por partes.) Esta ecuación se denomina ecuación de Bellman y se debe satisfacer para todo t ∈

si v

ha de ser una solución. Nótese que la estructura de la ecuación es invariante ante desplazamientos temporales, por lo que se refiere a la determinación de la función de valor, no hay inconveniente en prescindir del índice t . Hacemos: y = xt −1 y x = xt , con esto la ecuación de Bellman es

v( y ) = max [ g ( y, x) + β v( x) ] . x∈ f ( y )

Sea X = C

*

(

n

,

)

el espacio de las funciones continuas y acotadas, el cual dotado de la

métrica del sup es completo. De esta forma asumiendo que el max esta bien definido para cualquier “ y ”, podemos definir unívocamente el operador de Bellman como T : X → X dado por

T ( w ) ( y ) = max [ g ( y, x) + β w( x) ] , w ∈ X . x∈ f ( y )

-82-

En lo que sigue asumiremos que la función g y la correspondencia semicontinua superior f son de tal tipo que max [ g ( y , x) + β w( x) ] está bien definido y es una función continua y x∈ f ( y )

acotada de y . Luego, si fuese posible mostrar que T es una contracción en X , el teorema del punto fijo de Banach garantiza la existencia de una única solución a la ecuación de Bellman. Afirmación: T es una contracción en X . En efecto; sean w1 , w2 ∈ X

T ( w1 )( y ) − T ( w2 )( y ) = max [ g ( y, x) + β w1 ( x) ] − max [ g ( y, x) + β w2 ( x)] x∈ f ( y )

(II.6.4.2)

x∈ f ( y )

Para todo x ∈ f ( y ) , sea

u y ( x) = g ( y, x) + β w1 ( x) . Observemos que u y es continua sobre el compacto f ( y ) , luego u y tiene máximo. Digamos que x ∈ f ( y ) es el punto donde u y alcanza su máximo, entonces tenemos *

max u y ( x) = max [ g ( y, x) + β w1 ( x)] = u y ( x* ) = g ( y, x* ) + β w1 ( x* ) .

x∈ f ( y )

x∈ f ( y )

Análogamente para

v y ( x) = g ( y, x) + β w2 ( x) , x ∈ X ( y )

se tiene que x ∈ f ( y ) es el punto donde v y alcanza su máximo, **

max v y ( x) = max [ g ( y, x) + β w2 ( x)] = u y ( x** ) = g ( y, x** ) + β w2 ( x** ) .

x∈ f ( y )

x∈ f ( y )

Luego reemplazando en (II.6.4.2), tenemos

T ( w1 )( y ) − T ( w2 )( y ) ≤ ( g ( y, x* ) + β w1 ( x* ) ) − ( g ( y, x** ) + β w2 ( x** ) )

(II.6.4.3)

Ahora puede ocurrir que:

i Supongamos que

g ( y, x* ) + β w1 ( x* ) ≥ g ( y, x** ) + β w2 ( x** ) ; como g ( y , x ) + β w2 ( x ) es el máximo de v y ( x ) para todo x ∈ f ( y ) , en particular **

**

para x ∈ f ( y ) tenemos *

g ( y, x* ) + β w1 ( x* ) ≤ g ( y, x** ) + β w2 ( x** ) . -83-

Reemplazando esta última relación en (II.6.4.3) conseguimos

T ( w1 )( y ) − T ( w2 )( y ) ≤ ( g ( y, x* ) + β w1 ( x* ) ) − ( g ( y, x* ) + β w2 ( x* ) ) = β w1 ( x* ) − w2 ( x* ) ≤ sup w1 ( x) − w2 ( x) x

= β d ( w1 , w2 ) . i Si por el contrario, se verificase

g ( y, x* ) + β w1 ( x* ) < g ( y, x** ) + β w2 ( x** ) ; como g ( y , x ) + β w1 ( x ) es el máximo de u y ( x) para todo x ∈ f ( y ) , en particular *

*

para x ∈ f ( y ) , tenemos **

g ( y, x** ) + β w2 ( x** ) ≤ g ( y, x* ) + β w1 ( x* ) . Luego de reemplazar esta desigualdad en (II.6.4.3), tenemos

T ( w1 )( y ) − T ( w2 )( y ) ≤ ( g ( y, x** ) + β w2 ( x** ) ) − ( g ( y, x** ) + β w1 ( x* ) ) = β w2 ( x** ) − w1 ( x** ) ≤ sup w2 ( x) − w1 ( x) x

= β d ( w2 , w1 ) = β d ( w1 , w2 ) De ambos casos al final siempre tenemos

T ( w1 )( y ) − T ( w2 )( y ) ≤ β d ( w1 , w2 ) ; w1 , w2 ∈ X . Luego, el primer miembro es acotado, entonces tiene supremo, y

sup T ( w1 )( y ) − T ( w2 )( y ) ≤ β d ( w1 , w2 ) y

d (Tw1 , Tw2 ) ≤ β d ( w1 , w2 ) ; w1 , w2 ∈ X . Pero como tenemos por hipótesis que 0 ≤ β < 1, T es una contracción en C

*

(

n

,

)= X .

Y por el teorema del punto fijo de Banach (T.P.F.B.), existe una única solución a la ecuación de Bellman. ■

-84-

II.7. APLICACIÓN A LA DINÁMICA COMPLEJA II.7.1. FUNCIONES HOLOMORFAS DE VARIAS VARIABLES COMPLEJAS

=

× ×

entonces

es una coordenada compleja de z . Como z j = x j + iy j ∈

donde

x j = Re( z j ) , y j = Im( z j ) ∈ n

, y si un z = ( z1 , z2 ,..., zn ) ∈

n −veces

cada z j con j ∈

Si z ∈

×

n

n

Recordemos que

, luego tenemos que:

entonces ( x1 , y1, x2 , y2 ,..., xn , yn ) ∈

Para z = ( z1 , z2 ,..., zn ) , w = ( w1 , w2 ,..., wn ) ∈

n

2n

y

i.e. identificamos a

n

con

2n

.

α ∈ , definimos

⎧ z + w = ( z1 + w1 , z2 + w2 ,..., zn + wn ) ⎨ ⎩α z = (α z1 ,α z2 ,...,α zn ) n

con estas operaciones n

Dotaremos a

− espacio vectorial.

es un

de una topología, para esto definimos: n

i). Un polidisco abierto en con r = (r1 , r2 ,..., rn ) ∈

de centro a = ( a1 , a2 ,..., an ) ∈

( )

+ n

n

con r = ( r1 , r2 ,..., rn ) ∈

+ n

y poliradio r ∈ (

)

denotado por Δ ( a, r ) , como el conjunto

Δ ( a, r ) = { z = ( z1 , z2 ,..., zn ) ∈ ii). Un polidisco cerrado en

n

n

: z1 − a1 < r1 , z2 − a2 < r2 ,…, zn − an < rn } .

de centro a = ( a1 , a2 ,..., an ) ∈

( )

n

y poliradio r ∈ (

+ n

denotado por Δ [ a, r ] , como el conjunto

+ n

Δ [ a, r ] = { z = ( z1 , z2 ,..., zn ) ∈

n

: z1 − a1 ≤ r1 , z2 − a2 ≤ r2 ,…, zn − an ≤ rn } .

II.7.1.2. OBSERVACIONES II.7.1.2.1.- Si z = ( z1 , z2 ,..., zn ) ∈

n

la norma de z denotada por z está definida por

z =

2

2

z1 + z2 +

+ zn

es la norma euclidiana. II.7.1.2.2.- Si denotamos:

{ ⎡⎣ a ⎤⎦ = { z ∈

i Drj ( a j ) = z j ∈ : z j − a j < rj , j ∈ i Drj

j

j

: z j − a j ≤ rj , j ∈ -85-

}⊂ }⊂

. .

2

)

Obtenemos respectivamente:

× Drn ( an ) .

i Δ ( a, r ) = Dr1 ( a1 ) × Dr2 ( a2 ) ×

× Drn [ an ] .

i Δ [ a, r ] = Dr1 [ a1 ] × Dr2 [ a2 ] ×

Por ejemplo, gráficamente mostremos al polidisco abierto Δ ( a, r ) ⊂

II.7.1.2.3.-

n

2

:

dotado de la topología cuya base es generada por los polidiscos abiertos es un

espacio topológico equivalente a

2n

dotado de la topología cuya base es

generada por las bolas abiertas i.e. tenemos: n

≈ ↑

2n

topológicamente De esta manera todos los resultados conocidos de la topología de los espacios euclidianos

2n

pueden ser aplicados a

n

.

II.7.2. NOTACIÓN DE MULTI-ÍNDICES DE SCHWARTZ

i Un multi-índice de dimensión n , es una n -upla de números naturales Q = ( q1 , q2 ,…, qn ) ∈ (

i Su norma se define como Q = q1 + q2 + … + qn . -86-

)

+ n . 0

i El factorial de Q se denota por Q ! y se define Q ! = q1 !q2 !… qn !

i Sea U ⊆

n

abierto y f : U ⊆

n

→

, definimos

∂Q f ∂ f . = ∂z Q ∂z1q1 ∂z2q2 ...∂znqn Q

II.7.3. DEFINICIÓN Sea U ⊆

n

abierto y f : U →

. Decimos que f es holomorfa (ó analítica) en

a = (a1 , a2 ,..., an ) ∈U si existe Δ ( a, r ) ⊆ U polidisco abierto tal que ∞

f ( z ) = ∑ cQ ( z − a )Q ; z ∈ Δ ( a, r ) ⊆ U . Q =0

Decimos que f es holomorfa en U si es holomorfa en a , para todo a ∈U . Al conjunto de todas las funciones holomorfas en el abierto U ⊆

n

lo denotaremos por

O (U ) . II.7.3.1. EJEMPLO Las funcionales lineales f :

n

→

definidas por

f ( z1 , z2 ,..., zn ) = a1 z1 + a2 z2 + son funciones holomorfas en

n

+ an zn

.

En general, toda función polinomial es holomorfa en

n

.

II.7.4. OBSERVACIÓN De la teoría elemental de serie de potencias (en una variable en

) tenemos el

siguiente resultado: ∞

“Sea

∑ c ( z − a) n =0

n

n

una serie de potencias convergente en el disco DR ( a ) ⊆

.

Si 0 < r < R entonces la serie converge absolutamente y uniformemente en el disco cerrado

Dr [ a ] ”. Este resultado se generaliza a serie de potencias de varias variables:

-87-

∞

“Sea

∑ c ( z − a)

Q =0

una serie de potencias convergente en el polidisco abierto Δ ( a, R ) con

Q

n

poliradio R = ( R1 , R2 ,

, Rn ) ∈ (

)

+ n 0

, si 0 < rj < R j , para todo 1 ≤ j ≤ n entonces la serie

converge absoluta y uniformemente en el polidisco cerrado Δ [ a, r ] , r = (r1 , r2 ,..., rn ) ”. Por lo tanto, tenemos: “Si U ⊆

n

es un abierto y f : U →

es holomorfa en U entonces f es continua.”

“i.e. O (U ) ⊆ C (U ) ”. II.7.5. DEFINICIÓN Sea U ⊆

n

un abierto. Un campo vectorial holomorfo en U es una función

Z :U → n z Z ( z ) = ( Z1 ( z ), Z 2 ( z ),..., Z n ( z ) ) tal que 1).- Z j : U →

son funciones holomorfas para todo 1 ≤ j ≤ n .

2).- Si z ∈U entonces Z ( z ) ∈

n

es un vector cuyo punto de aplicación es z .

II.7.5.1. NOTACIÓN Denotemos por Χ (U ) al conjunto de todos los campos vectoriales holomorfos. i.e. Χ (U ) = { Z : U →

n

; Z es campo vectorial holomorfo en U }.

II.7.5.2. EJEMPLOS II.7.5.2.1.- Sea Z :

2

→

2

definido por

Z ( z1 , z2 ) = (az1 + bz2 , cz1 + dz2 ) , para cada ( z1 , z2 ) ∈ Sigue que Z ∈ Χ ( II.7.5.2.2.- Sea A = ⎡⎣ aij ⎤⎦ ∈

2

2

donde a , b , c , d ∈

.

) , el cual se denomina campo lineal en

n×n

, y definimos Z :

n

z ⎛

→

2

.

n

.

n

Z ( z ) = Az

n n ⎞ a z , a z ,..., ank zk ⎟ . ∑ ∑ 1k k ∑ 2 k k k =1 k =1 ⎝ k =1 ⎠

i.e. Z ( z ) = ⎜ Sigue que Z ∈ Χ (

n

n

) , el cual se denomina campo lineal de -88-

Note que el ejemplo II.7.5.2.1. es un caso particular del ejemplo II.7.5.2.2., para

⎡a b ⎤ ⎥∈ ⎣c d ⎦

la matriz A = ⎢

II.7.5.2.3.- Sea A = ⎡⎣ aij ⎤⎦ ∈

n×n

2×2

.

, b = [bi ] ∈

n×1

≈

n

definimos Z :

Z ( z ) = Az + b , z ∈

n

n

→

n

por

.

n n ⎛ n ⎞ i.e. Z ( z ) = ⎜ ∑ a1k zk + b1 , ∑ a2 k zk + b2 ,..., ∑ ank zk + bn ⎟ . k =1 k =1 ⎝ k =1 ⎠

Sigue que Z ∈ Χ (

n

) el cual se denomina campo afín en

n

.

II.7.6. DEFINICIÓN Sea U ⊆

n

Z si Z ( z ) = 0 ∈

un abierto y Z ∈ Χ (U ) . Decimos que z ∈U es un punto singular de n

(i.e. z anula el campo). Caso contrario, decimos que z es un punto

regular de Z . Denotaremos por Sing ( Z ) al conjunto de todos los puntos singulares de Z . i.e. Sing ( Z ) = { z ∈U : Z ( z ) = 0} ⊆ U . II.7.6.1. OBSERVACIONES II.7.6.1.1.- Dado Z = ( Z1 , Z 2 ,..., Z n ) ∈ Χ(U ) . Decir que

z = ( z1 , z2 ,..., zn ) ∈ Sing ( Z ) significa

⎧ Z1 ( z1 , z2 ,..., zn ) = 0 ⎪ Z ( z , z ,..., z ) = 0 ⎪ 2 1 2 n ⎨ ⎪ ⎪⎩ Z n ( z1 , z2 ,..., zn ) = 0 II.7.6.1.2.- Si

Z ∈ Χ(

n

) es un campo lineal (ver

θ = (0,0,...,0) ∈ II.7.6.1.3.- Si Z ∈ Χ (

n

n

ejemplo

II.7.5.2.2.) entonces

es un punto singular de Z .

) es un campo afín (ver ejemplo II.7.5.2.3) con b ≠ 0 entonces

θ = (0,0,...,0) ∈

n

es un punto regular de Z . -89-

II.7.6.2. EJEMPLOS 2

II.7.6.2.1.- Sea Z :

→

2

definido por

Z ( z1 , z2 ) = (a, z2 ) para todo ( z1 , z2 ) ∈

2

− {0} . Es claro que Z ∈ Χ(

con a ∈

2

) y que Z no

tiene puntos singulares. 2

II.7.6.2.2.- Sea Z :

→

2

definido por

Z ( z1 , z2 ) = ( z1 − z23 , z1 z2 ) , ( z1 , z2 ) ∈ Entonces Z ∈ Χ

2

.

( ) y (0,0) es el único punto singular de Z . 2

En efecto; sea ( w1 , w2 ) ∈ Sing ( Z ) entonces tenemos

Z ( w1 , w2 ) = 0 , de donde obtenemos

( w1 − w23 , w1w2 ) = (0,0) , esta igualdad se cumple si y sólo si

⎧ w1 − w23 = 0 ⎨ ⎩ w1w2 = 0 donde de la segunda igualdad obtenemos

w1 = 0 ó w2 = 0 lo que implica

w23 = 0 ó w1 = 0 respectivamente, de donde

( w1 , w2 ) = (0,0) . Por lo tanto

Sing ( Z ) = {(0,0)} . II.7.7. DEFINICIÓN Sea U ⊆

n

abierto y Z ∈ Χ (U ) .

i).- La ecuación diferencial ordinaria asociada a Z es dada por

z′ = Z ( z )

-90-

(II.7.7.A)

ii).- Una solución de la E.D.O. (II.7.7.A) es una función holomorfa un disco abierto D ⊆

ϕ : D → U definida en

tal que para todo T ∈ D

ϕ ′(T ) = Z (ϕ (T ) ) . II.7.7.1. OBSERVACIONES II.7.7.1.1.- z′ =

dz ,T∈ . dT

II.7.7.1.2.- Sea Z = ( Z1 , Z 2 ,..., Z n ) ∈ Χ (U ) donde cada Z j : U →

z = ( z1 , z2 ,..., zn ) ∈

con j = 1,2,..., n y

entonces (II.7.7.A) es equivalente a

⎧ z1′ = Z1 ( z1 , z2 ,..., zn ) ⎪ z′ = Z ( z , z ,..., z ) ⎪ 2 2 1 2 n ⎨ ⎪ ⎪⎩ zn′ = Z n ( z1 , z2 ,..., zn ) II.7.7.1.3.- Si

ϕ:D→

n

definido en un disco abierto D ⊆

(II.7.7.1.2.A)

con T ∈

tal que

ϕ (T ) = (ϕ1 (T ),ϕ2 (T ),...,ϕn (T ) ) es solución de (II.7.7.A) ó (II.7.7.1.2.A) entonces

⎧ϕ1′(T ) = Z1 (ϕ1 (T ),ϕ 2 (T ),...,ϕ n (T )) ⎪ϕ ′ (T ) = Z (ϕ (T ),ϕ (T ),...,ϕ (T )) ⎪ 2 n 2 1 2 ; T ∈D ⊆ ⎨ ⎪ ⎩⎪ϕn′ (T ) = Z n (ϕ1 (T ),ϕ 2 (T ),...,ϕ n (T )) Gráficamente tenemos:

-91-

.

II.7.7.1.4.- Si Z ∈ Χ (U ) es lineal entonces la E.D.O. asociada a Z se llama E.D.O. lineal. II.7.7.1.5.- Si Z ∈ Χ (U ) es afín entonces la E.D.O. asociada a Z se llama E.D.O. afín. II.7.8. DEFINICIÓN Sea U ⊆

n

abierto y Z ∈ Χ (U ) , para todo z0 ∈U y cualquier T0 ∈

.

i).- El problema de valor inicial (P.V.I.) asociado a Z es dado por

⎧ z′ = Z ( z ) ⎨ ⎩ z (T0 ) = z0

(II.7.8.1)

(II.7.8.1) es llamado Problema de Cauchy para Campos Vectoriales Holomorfos. ii).- La solución del P.V.I. (II.7.8.1) es una función holomorfa

ϕ : D → U donde D ⊆ ,

es un abierto, tal que para todo T0 ∈ D

⎧ϕ ′ = Z (ϕ (T )), T ∈ D ⎨ ⎩ϕ (T0 ) = z0 II.7.9. PROPOSICIÓN Sea U ⊆

n

abierto, Z ∈ Χ (U ) , para todo z0 ∈U y cualquier T0 ∈

.

Entonces resolver el P.V.I. (II.7.8.1) es equivalente a resolver la ecuación integral T

z (T ) = z0 + ∫ Z ( z (τ ))dτ

(II.7.9.1)

T0

Demostración: Sea

ϕ solución de (II.7.9.1) entonces existe D ⊆ ϕ : D →U T ϕ (T ) -92-

abierto, T0 ∈

tal que

es holomorfa, y

⎧ϕ ′ = Z (ϕ (T )), T ∈ D ⎨ ⎩ϕ (T0 ) = z0 Sea T ∈ D y consideremos el camino recto que une T0 y T :

∫

T

T0

T

ϕ ′(τ )dτ = ∫ Z (ϕ (τ ))dτ T0

entonces

ϕ (T ) − ϕ (T0 ) = ∫ Z (ϕ (τ ) )dτ T

T0

de donde obtenemos

ϕ (T ) = z0 + ∫ Z (ϕ (τ ) )dτ . T

Aquí tomamos T = t0 , obteniendo

Por lo tanto

T0

ϕ (t0 ) = z0

ϕ es solución de (II.7.9.1). El recíproco es evidente. ■

II.7.10. DEFINICIÓN Sea U ⊆

n

abierto, y Z : U →

n

un campo vectorial holomorfo.

II.7.10.1.- Decimos que Z es Lipschitz en U si y sólo si existe c > 0 tal que

Z ( z ) − Z ( w) ≤ c z − w ; z , w ∈U . II.7.10.2.- Decimos que Z es localmente Lipschitz en U si y sólo si para cualquier z0 ∈U existe el polidisco abierto de centro z0 y poliradio r ∈

Z|Δ( z0 ,r ) : Δ ( z0 , r ) →

( )

n +

tal que

n

es Lipschitz en Δ ( z0 , r ) . II.7.10.1. OBSERVACIONES II.7.10.1.1.- Todo campo Lipschitz es localmente Lipschitz. II.7.10.1.2.- Si Z : U →

n

es lipschitziana en U ⊆

n

, entonces el conjunto

⎧⎪ Z ( z ) − Z ( w) ⎫⎪ ∈ ≠ : z , w U , z w ⎨ ⎬⊆ z−w ⎪⎩ ⎪⎭ es acotado superiormente. El supremo de este conjunto es llamado constante de Lipschitz de Z y lo denotamos por Lip ( Z ) . -93-

II.7.10.1.3.- En general se cumple

Z ( z ) − Z ( w) ≤ Lip ( Z ) z − w ; z , w ∈U . II.7.10.1.4.- Si Z : U →

n

es localmente lipschitziana entonces la constante de Lipschitz de

Z|Δ( z0 ,r ) depende del polidisco abierto Δ( z0 , r ) . II.7.11. PROPOSICIÓN Sea U ⊆

n

abierto, Z : U →

n

un campo vectorial holomorfo. Entonces Z es

localmente lipschitziana. Demostración: Sea z0 ∈U entonces existe R = ( R1 , R2 ,..., Rn ) ∈

( )

+ n

tal que Δ [ z0 , R ] ⊂ U .

Como Z es holomorfa en U , entonces Z ′ es continua en U y Δ [ z0 , R ] es un compacto, entonces Z ′ es acotada en Δ [ z0 , R ] i.e. existe M > 0 tal que

Z ′(ω ) ≤ M , ω ∈ Δ [ z0 , R ] . Sean z , w ∈ Δ [ z0 , R ] entonces el segmento de recta [ z , w] ⊂ Δ [ z0 , R ] . Luego definimos

ψ : [ 0,1] → n t ψ (t ) = Z ( (1 − t ) z + tw )

claramente ψ es holomorfa. Además

ψ ′(t ) = Z ′ ( (1 − t ) z + tw ) ( w − z ) , t ∈ [ 0,1] . Entonces 1

Z ( w) − Z ( z ) = ψ (1) − ψ (0) = ∫ ψ ′( s )ds 0

1

Z ( w) − Z ( z ) ≤ ∫ ψ ′( s ) ds 0

1

= ∫ Z ′ ( (1 − s ) z + sw ) ( w − z ) ds 0

1

≤ ∫ Z ′ ( (1 − s ) z + sw ) ( w − z ) ds 0

1

≤ ∫ M ( w − z ) ds = M ( w − z ) 0

Z ( w) − Z ( z ) ≤ M w − z ; z , w ∈ Δ ( z0 , R ) . Por lo tanto Z es localmente lipschitziana. ■ -94-

II.7.12. OBSERVACIONES II.7.12.1.- Sea K ⊆

n

compacto y denotemos

C ( Dα [T0 ] , K ) = { φ : Dα [T0 ] → K ; φ continua}. Para

ϕ ,ψ ∈ C ( Dα [T0 ] , K ) definimos d (ϕ ,ψ ) = max

T ∈Dα [T0 ]

{ ϕ (T ) −ψ (T ) } .

( C ( D [T ], K ) , d ) resulta ser un espacio métrico completo. α

0

II.7.12.2.- Sea K ⊆

n

compacto y denotemos

A ( Dα [T0 ] , K ) = { φ : Dα [T0 ] → K ; φ continua y holomorfa}. Es claro que

A ( Dα [T0 ] , K ) ⊆ C ( Dα [T0 ] , K ) .

( (

) )

( (

) )

Consideremos el espacio métrico A Dα [T0 ] , K , d donde d es la métrica

(

)

de C Dα [T0 ] , K . Probaremos que A Dα [T0 ] , K , d es completo.

(

)

Para esto, es suficiente probar que A Dα [T0 ] , K es cerrado. Recordemos dos resultados de variable compleja: a) Sea U ⊆

abierto,

φk : U →

una función. Decimos que

n

una sucesión de funciones y

φ :U →

n

φk converge a φ uniformemente en las partes

compactas de U lo que denotamos por

φk → φ u.p.c. de U si y sólo si

φk |K → φ |K uniformemente en K , para todo K ⊆ U compacto (ver [24]). b) Sea U ⊆ en U y

abierto, si

φ :U →

n

φk : U →

es tal que

n

es una sucesión de funciones holomorfas

φk → φ u.p.c. de U entonces φ es holomorfa

en U . Teniendo presente estos dos resultados continuemos con la demostración. Sea

φ ∈ A ( Dα [T0 ] , K ) entonces existe (φk ) ⊆ A ( Dα [T0 ] , K ) tal que lim d (φk ,φ ) = 0 . k →∞

-95-

Entonces tenemos Afirmación:

φ ∈ C ( Dα [T0 ] , K ) .

φ es holomorfa en Dα (T0 ) .

En efecto, sea K ′ ⊆ Dα (T0 ) compacto. Como

lim d (φk ,φ ) = 0 , k →∞

tenemos que

φk |K′ → φ |K′ uniformemente en K ′ entonces por a). φk → φ u.p.c de

Dα (T0 ) , luego por b). φ ∈ A ( Dα [T0 ] , K ) .

( (

) )

Por lo tanto, concluimos que A Dα [T0 ] , K , d es un espacio métrico completo. II.7.13. TEOREMA (PICARD) Sea Z : Δ [ z0 , r ] →

n

un campo vectorial Lipschitz en Δ [ z0 , r ] y holomorfo en

Δ ( z0 , r ) entonces para todo T0 ∈

existe una única solución del P.V.I. (I.7.8.1)

⎧ z′ = Z ( z ) ⎨ ⎩ z (T0 ) = z0 definida en el disco Dα [T0 ] , donde

{

r ⎫ ⎧ r1 r2 , ,..., n ⎬ , r = ( r1 , r2 ,..., rn ) ∈ ( N⎭ ⎩N N

α = min ⎨

}

y N = max Z ( z ) : z ∈ Δ [ z0 , r ] .

)

+ n

Demostración: Sabemos por la proposición (II.7.9) que el P.V.I. (II.7.8.1) es equivalente a T

z (T ) = z0 + ∫ Z ( z (τ ))dτ . T0

( (

) )

Consideremos el espacio métrico completo X = A Dα [T0 ] , K , d . Dado

φ ∈ X definimos Fφ : Dα [T0 ] →

n

por T

Fφ (T ) = z0 + ∫ Z (φ (τ ))dτ , T0

claramente Fφ es continua en Dα [T0 ] y holomorfa en Dα (T0 ) . Para z = ( z1 , z2 ,..., zn ) ∈

n

denotemos por

π j ( z ) a la proyección j − ésima definida por

π j ( z ) = z j para todo 1 ≤ j ≤ n . -96-

Ahora hacemos

π j ( Fφ (T )) − z 0j =

∫

T

T0

Z j (φ (τ ))dτ

T

T

T0

T0

≤ ∫ Z j (φ (τ )) dτ ≤ max Z j (φ (τ )) ∫ dτ ≤ N T − T0 ≤ Nα ≤ rj entonces Fφ ∈ X , luego hemos construido

F:X →X φ F (φ ) = Fφ Vamos a probar que F es continua y que existe m ∈

+

tal que F

m

es una contracción

en X . Afirmación:

Lip ( Z ) T − T0 F n (φ1 ) (T ) − F n (φ2 ) (T ) ≤ d (φ1 ,φ2 ) ; n! n

n

(II.7.13.1)

para cada n ≥ 0 , para cualesquiera φ1 ,φ2 ∈ X y para todo T ∈ Dα [T0 ] . En efecto:

i Para n = 0 ¡es evidente!

i Supongamos que la desigualdad se verifica para n = k . Probaremos que se cumple para n = k + 1. Para cualesquiera φ1 ,φ2 ∈ X y cada T ∈ Dα [T0 ] tenemos

(

)

(

π j ( F k +1 (φ1 )(T ) ) − π j ( F k +1 (φ2 )(T ) ) = π j F ( F k (φ1 ) ) (T ) − π j F ( F k (φ2 ) ) (T ) ≤ ∫ Z j ( F k (φ1 )(τ ) ) − Z j ( F k (φ2 )(τ ) ) dτ T

T0 T

≤ ∫ Lip ( Z ) F k (φ1 )(τ ) − F k (φ2 )(τ ) dτ T0

Lip ( Z ) k k ≤ Lip ( Z ) ∫ d (φ1 ,φ2 ) τ − T0 dτ T0 k! T

=

T Lip ( Z ) k +1 k d (φ1 ,φ2 ) ∫ τ − T0 dτ T0 k!

Lip ( Z ) k +1 d (φ1 ,φ2 ) . = (k + 1)! Lo cual prueba la afirmación. -97-

)

Hacemos n = 1 en (II.7.13.1) y tenemos

F (φ1 ) (T ) − F (φ2 ) (T ) ≤ Lip ( Z )α d (φ1 ,φ2 ) entonces

d ( F (φ1 ), F (φ2 ) ) ≤ ( Lip Z α ) d (φ1 ,φ2 )

por lo tanto F es continua. Análogamente tenemos

Lip ( Z ) α F n (φ1 ) (T ) − F n (φ2 ) (T ) ≤ d (φ1 ,φ2 ) . n! n

Para cada m ∈

+

n

, sea

um

( Lip( Z )α ) =

m

m!

,

utilizando el criterio de la razón para sucesiones de números reales tenemos que

lim

m→∞

um+1 = 0 n0 tenemos um < ε < 1 .

Entonces basta tomar m = n0 , para afirmar que F

m

es una contracción en X , luego por el

corolario I.4.2. F tendrá un único punto fijo, el cual es la solución del P.V.I. (II.7.8.1). ■ II.7.15. COROLARIO Sea U ⊆

n

abierto, Z ∈ Χ (U ) , para todo z0 ∈U , para todo T0 ∈

, existe una

única solución del P.V.I.(II.7.9.1) la cual está definida en una vecindad de T0 . Demostración: Como Z ∈ Χ (U ) entonces por proposición (II.7.12.) Z es localmente Lipschitz, luego existe Δ [ z 0 ,r ] ⊆ U tal que Z es Lipschitz en Δ [ z 0 ,r ] . Como Z es holomorfo en Δ ( z 0 ,r ) ⊆ U entonces por el teorema de Picard para campos holomorfos (teorema II.7.14.) se tiene que existe una única solución del P.V.I. (II.7.9.1). ■ Para más referencias sobre la teoría de dinámica compleja puede verse [3] y [15].

-98-

REFERENCIAS BIBLIOGRÁFICAS [1] ARIS R.; “The Optimal Design of Chemical Reactors”, Academic Press Inc, New York, 1961. [2] BELLMAN R.; “Dynamic Programing”; Princeton University Press., New Jersey 1957. [3] BENAZIC R.; “Singularidades de Campos Vectoriales Holomorfos en el Domino de Poincaré”, Pro Matemática, Vol. X, Nº 19-20, 1996. [4] BOYCE W. – DIPRIMA R.; “Ecuaciones Diferenciales y Problemas con Valores en la Frontera”, Instituto Politécnico Rensselear, Troy New York, 1978. [5] BURDEN R. – FAIRES J.; “Análisis Numérico”, Iberoamericana, México D.F, 1985. [6] CAMACHO C. – SAD, P.; “Puntos Singulares de Equaçoes Diferenciais Analiticas”, 16 Coloquio Brasileiro de Matemática, 1987. [7] CERDA E.; “Optimización Dinámica”, Prentice Hall, España, 2001. [8] CONTE S. - BOOR C.; “Análisis Numérico Elemental”, McGraw - Hill S.A., México, 1985. [9] CHAIM H.; “Aplicaçoes da Topologia à Análise”, Instituto de Matemática Pura e Aplicada (IMPA), Sao Paulo, 1976. [10] CHAIM H.; “ Functional Analysis and integro – differential with linear constraints”, Sociedad Brasilera de Matemática, Reuniao de Análise funcional, Campinas, 1974. [11] CHAIM H.; “Análise Funcional e Aplicaçoes (Volume II)”, Publicaçoes ou ediçoes do Instituto de Matemática e Estadística da Universidade de Sao Paulo, Sao Paulo, 1970. [12] CHIANG A.; “Elements of Dynamic Optimization”, McGraw-Hill, 1992. [13] CHUMPITAZ M.; “Análisis Funcional I”, Facultad de Ciencias, Universidad Nacional de Ingeniería. [14] CHURCHILL R. – WARD J.; “Variable Compleja y Aplicaciones”, Mc Graw-Hill, España, 1990. [15] GUNNING R.; “ Introduction to Holomorphic Functions of Several Variables, Vol I Function Theory”, Wardsworth&Brooks/Cole, California, 1990. [16] GUNNING R. – ROSSI H.; “ Analytic Functions of Several Complex Variables”, Prentice Hall, New York, 1965. [17] HASSER B. – SULLIVAN J.; “Análisis Real”, Editorial Trillas México, 1978. [18] HIRSCH M. – SMALE S.; “Differential Equations, Dynamical Systems and Linear algebra”, University of California, Berkeley Academic, INC. New York, 1974.

-99-

[19] KINCAID D. – CHENNEY W.; “Análisis Numérico: Las matemáticas del cálculo científico”, The University of Texas in Austin, Addison – Wesley Iberoamericana Wilmigton, Delaware, E.U.A., 1994. [20] KOLMOGOROV K. – FOMIN S.; “Elements of Theory of Functions and Functional Analysis – Volume 2”, University Wisconsin –Rochester, New York, 1957. [21] KREYSZIG E.; “Introductory Functional Analysis with Applications”, University to Windsor, New York, 1989. [22] LIMA E.; “Curso de Análise-Vol. 2”, Instituto de Matemática Pura e Aplicada (IMPA), Rio de Janeiro, 1976. [23] LIMA E.; “Espaços Métricos”, Instituto de Matemática Pura e Aplicada (IMPA), Rio de Janeiro, 1976. [24] LINS A.; “Funçoes de uma Variável Complexa”, Proyecto Euclides, IMPA, Rio de Janeiro, 1977. [25] MEDEIROS L. – DE ANDRADE G.; “Iniciaçao às Equaçoes Diferenciais Parciais”, Livros Técnicos e Científicos Editora; Rio de Janeiro, 1978. [26] NOLASCO A.; “Análise I”, ICMC, USP SAO CARLOS, Sao Carlos, 2004. [27] NOWOSAD P.; “Introduçao à Análise Funcional”, Universidad Nacional Mayor de San Marcos, Lima, 1967. [28] QUARTERONI A. – SACCO R. – SALERI F.; “Numerical Mathematics”, Springer – Verlag, New York Inc, 2000. [29] REDDY J. – GARVIN C.; “Applied Functional Analysis and Variational Methods in Engineering”, Virginia Polytechnic Institute and State, New York, 1986. [30] SOTOMAYOR J.; “Liçoes de Equaçoes Diferenciais Ordinárias”, Proyecto Euclides, IMPA, Rio de Janeiro, 1979. [31] SHASHKIN YU.; “Lecciones populares de matemáticas: Puntos Fijos”, editorial MIR, Moscú, 1991. [32] VILENKIN N.; “Lecciones populares de matemáticas: Método de aproximaciones sucesivas”, Editorial MIR, Moscú, 1991.

-100.