Textocfinal II2007

Universidad Mayor de San Simon
Fa ultad de Cien ias y Te nologa
Departamento de Matemati as
Cal ulo III

E ua iones Diferen iales
Hans C. Muller Santa Cruz
Co habamba, 1998.
Contenido
Prefa io : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
E ua iones Diferen iales Ordinarias : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
1.- Vo abulario Basi o : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
2.- E ua iones Diferen iales de Primer Orden : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
3.- Problemas de Existen ia y Uni idad : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
4.- E ua iones de Orden Superior : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
5.- Ejer i ios : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
II.- Sistemas Diferen iales y Apli a iones : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
1.- Con eptos Basi os : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
2.- Apli a iones : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
3.- Sistemas Diferen iales Lineales : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
4.- Estabilidad y Estudio de Puntos Crti os : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
5.- Ejer i ios : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
III.- Complementos : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
1.- E ua iones Diferen iales Exa tas : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
2.- E ua iones a Derivadas Par iales : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
3.- Elementos de Cal ulo Varia ional : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
4.- Ejer i ios : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
5.- Teorema de Hamilton Cayley : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
6.- Endomorsmos Nilpotentes : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
7.- Des omposi ion de Jordan : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
I.-
iii
1
1
4
13
20
33
39
39
44
51
60
66
71
71
77
82
98
103
105
112
Prefa io
El ontenido de la asignatura de Cal ulo III que se imparte en la Fa ultad de Cien ias y Te nologa de la
Universidad Mayor de San Simon, esta orientado al estudio y solu ion de problemas diferen iales. La importan ia de esta materia radi a en el he ho en que las e ua iones diferen iales onstituyen un lenguaje natural
para des ribir la mayor parte de los fenomenos, asimismo onstituye un medio para prede ir omportamientos
de los diferentes pro esos que debera analizar el estudiante en su desempeno profesional ulterior.
El presente texto es el resultado del urso de Cal ulo III impartido durante varios semestres a los
estudiantes de la Carrera de Ingeniera Civil y las Carreras de Fsi a, omo tambien en su momento la Carrera
de Matemati as. Contemplando el umplimiento de los objetivos generales y espe i os de la asignatura, el
texto uenta on tres aptulos.
El primer aptulo presenta las e ua ioenes diferen iales ordinarias, dando el vo abulario basi o y las
te ni as elementales de solu ion para las e ua iones diferen iales mas usuales.
El segundo aptulo pretende mostrar que las e ua iones diferen iales onstituyen una herramienta
poderosa para resolver otro tipos de problemas, que no son diferen iales propiamente. Asimismo se introdu e
el estudio de los sistemas diferen iales bajo una opti a moderna orientada a la implementa ion de metodos
numeri os.
El ter er aptulo es un omplemento del urso donde se abordan temas a un nivel introdu torio.
Para un buen asimila ion de los ono imientos y razonamientos de este texto; las deni iones y on eptos
mas signi ativos estan es ritos en negrillas, estos deberan ser memorizados y manipulados uidamente.
Los resultados mas importantes estan expresados en los teoremas, orolarios y proposi iones, estos deberan
tambien ser memorizados para manejarlos de manera uida. Las demostra iones de este texto deberan ser
trabajadas, on la nalidad de adquirir las diferentes te ni as de demostra ion y resolu ion de problemas
que se emplean en E ua iones Diferen iales. Con nes pedagogi os, en algunos paragrafos se presentan los
resultados fundamentales que seran tratados en el paragrafo en uestion, estos estan es ritos en ara teres
it
ali os.
La pra ti a del urso, es una fuente para pra ti ar los ono imientos adquiridos y as mismo omo un
medio de adquirir ono imientos adi ionales. Por lo tanto, una resolu ion en gran numero de estos ejer i ios,
podra medir el grado de asimila ion del estudiante.
Cap
tulo I
E ua iones Diferen iales Ordinarias
En este primer aptulo, se estudiara las e ua iones diferen iales usuales que uno en uentra en la mayora de
los problemas, donde intervienen e ua iones diferen iales. Para tal efe to, se des ribira el vo abulario basi o,
los metodos de solu ion de las e ua iones mas utilizadas y algunos metodos para aproximar solu iones de
algunos problemas diferen iales.
I.1 Vo abulario B
asi o
Para denir lo que es una e ua ion diferen ial, introduz amos la siguiente nota ion para representar las
derivadas (par iales) de una fun ion. Sea f : U R ! R una fun ion, onsideremos el onjunto de
= (1 ; 2 ; : : : ; n ); onk 0 y k 2 Z; k = 1; : : : ; n:
Denimos
jj f
(x);
D f (x) = 1 2
x x xn
donde x = (
x 1 ; x2 ; : : : ; x n 2 R n
y jj = 1 + 2 + + n.
Remar as
1.- Si n = 1, onvenimos que 2 Z, 0 y
D f (x) = f () (x);
la derivada (ordinaria) de orden de f .

2.- El entero jj es el orden de la derivada (par ial) de f .
Ejemplo
1.- Si n = 2 y 1 = n y 2 = m, se tiene
D f (x) =
n+m f
(x):
xn1 xm
2
Una e ua ion diferen ial es una e ua ion donde apare en una fun ion in ognita y y
algunas de sus derivadas; es de ir, es una expresion
F (x; y; D y; D y; : : : ; D y) = 0
(I.1.1)
donde F es una fun ion ontinua ono ida, x 2 Rn , y(x) la fun ion in ognita.
Deni ion I.1.1.-
I E ua iones Diferen iales Ordinarias
Ejemplo
2.- Para n = 1, la e ua ion y0 = f (x), donde f (x) es una fun ion ono ida es una e ua ion diferen ial; que
por ierto, ya resuelta en Cal ulo I. La solu ion se la ono e on el nombre de primitiva o integral
indenida. Por onsiguiente,
Z
y(x) = f (x) dx:
El estudio de las e ua iones diferen iales pasa por lo tanto, por:

1.- En ontrar una o todas las y on esta propiedad. Es de ir, busqueda de una \primitiva".
2.- >Existe solu iones?
3.- >Se las puede expli itar?
Deni ion I.1.2.- Si la fun ion in ognita y es fun ion de una sola variable, n = 1, la e ua ion diferen ial
es una e ua ion diferen ial ordinaria. Si y es fun ion de varias variables, n > 1, la e ua ion diferen ial es
una e ua ion diferen ial a derivadas par iales.
En lo que sigue este aptulo, solo trataremos e ua iones diferen iales ordinarias, de donde onvenimos
que por e ua ion diferen ial entenderemos e ua ion diferen ial ordinaria.
El Orden de una e ua ion diferen ial
El orden de una e ua ion diferen ial es el orden de deriva ion mas alto de la fun ion
in ognita que apare e en la e ua ion diferen ial.
Deni ion I.1.3.Ejemplos
3.- Las e ua iones algebrai as f (x; y) = 0 son e ua iones diferen iales de orden 0.
4.- y00 + y = 0, es de orden 2.
5.- (y0)3 = y2 es de orden 1 o primer orden.
6.- y00 + y0y = 1 es de orden 2.
Deni ion I.1.4.- Una e ua ion expl ita de orden n, es una e ua ion diferen ial de orden n de la forma
y(n) = f (x; y; y0 ; : : : ; y(n 1) );
(I.1.2)
donde f es una fun ion ontinua.
Solu iones de una E ua ion Diferen ial Ordinaria
Consideramos la e ua ion diferen ial (ordinaria y expl ita)

y(n) = f (x; y; y0 ; : : : ; y(n 1) );
(I.1.2)
donde f : A Rn+1 ! R es ontinua.
Una solu ion o solu ion parti ular de esta e ua ion diferen ial es una fun ion ' : I ! R, on I R
intervalo, n ve es derivable, tal que
'(n) (x) = f (x; '(x); '0 (x); : : : ; '(n 1) (x)) 2 A; 8x 2 I:
La solu ion general de una e ua ion diferen ial es el onjunto de sus solu iones parti ulares.
Remar as
1.- Si bien, la solu ion general de una e ua ion diferen ial es un onjunto de fun iones, es ostumbre denotar
la solu ion general bajo la forma de una fun ion que depende de un parametro o onstante.
2.- La e ua ion diferen ial puede ser dada en forma expl ita, sin que se pueda ha er lo mismo para sus
solu iones. Por ejemplo
2(y + 2)2 ;
y0 =
((y + 2) (x 1))2
si o
I.1 Vo abulario Ba
uyas solu iones estan dadas por
ln(y + 2) = 2 ar tan xy + 21 + C;
C 2 R;
y no puede ser expresado omo omposi ion de fun iones elementales; es de ir, las fun iones estudiadas
en Cal ulo I.
Problemas Diferen iales
En general, las e ua iones diferen iales estan ligadas a iertas ondi iones impuestas en el planteo del modelo
matemati o que des ribira el fenomeno o la a ion analizada. Para ilustrar, la ne esidad de ontar on iertas
ondi iones adi ionales a la e ua ion diferen ial, onsideremos omo ejemplo la ada libre de un objeto.
Supongamos que la gravedad es la uni a fuerza que interviene sobre el objeto, despre iando otras fuerzas.
Por onsiguiente, el movimiento del objeto veri a
a = g;
donde a es la a elera ion y g es la onstante de gravedad sobre la super ie de la Tierra. Tradu ido en forma
de e ua ion diferen ial, se tiene
d2 y
= g;
dt2
de donde, por dos integra iones su esivas, obtenemos
1
y(t) = gt2 + t + d;
(I.1.3)
2
on ; d 2 R onstantes.
Ahora bien, (I.1.3) no des ribe el movimiento del objeto en uestion; es de ir, on solamente la e ua ion
diferen ial no podemos saber la posi ion del objeto en ualquier instante. Requerimos informa ion adi ional,
aparte de la e ua ion diferen ial. En el urso de Fsi a I, este problema fue tratado uando se ono a, la
posi ion ini ial y la velo idad ini ial del objeto. Este es un ejemplo de un:
Deni ion I.1.5.- Un problema diferen ial a valor ini ial o problema diferen ial de Cau hy es darse una
e ua ion diferen ial de orden n
y(n) = f (x; y; y0; : : : ; y(n 1)
y los valores ini iales (x0 ; y0; y1; : : : ; yn 1) 2 Rn+1 ; bus ando una solu ion ' de la e ua ion diferen ial
propuesta, tal que
'(x0 ) = y0 ;
'0 (x0 ) = y1 ;
...
'(n 1) (x0 ) = yn 1 :
Remar a.- Los problemas diferen iales a valores ini iales son muy importantes, pues una vasta antidad de
problemas en me ani a, qumi a y otras ien ias son problemas a valores ini iales. Por otro lado, en Analisis
Numeri o, los metodos de solu ion de e ua iones diferen iales omtemplan este tipo de problema. Por lo
tanto, es importante responder a: >Dada una e ua ion diferen ial, se puede resolver un problema de Cau hy
para esta e ua ion?, >la solu ion de este problema de Cau hy es uni a?
Los problemas diferen iales a valor ini ial no son los uni os problemas diferen iales, existen problemas
diferen iales mas generales, ono idos omo problemas diferen iales on valores en la frontera o valoresn en
los bordes. Como ilustra ion estudiemos el siguiente ejemplo.
Ejemplo
7.- Consideremos una viga de longitud L en posi ion horizontal, sometida a una arga verti al q(x). Por
efe to de esta arga, la viga sufre una deforma ion transversal u(x). Suponiendo la deforma ion peque~na,
u(x) es solu ion de la e ua ion diferen ial
u(4)(x) = q(x):
Ahora bien, dependiendo lo que su ede en las extremidades de la viga, se tiene diferentes problemas
diferen iales on ondi iones en los bordes. En la gura I.1, se tiene las tres variantes mas usuales.
111
000
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
u(x)
111
000
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
u(x)
111
000
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
u(0)=u00 (0)=0
u00 (0)=u00 (1)=0:
u(0)=u0 (0)=0
u0 (0)=u0 (1)=0:
u(x)
u(0)=u0 (0)=0
u00 (1)=u000 (1)=0:
Figura I.1.- Problemas Diferen iales de Vigas
I.2 E ua iones Diferen iales de Primer Orden
Estudiaremos los metodos de resolu ion de las e ua iones de primer orden bajo la forma expl ita, mas
usuales.
E ua iones Lineales
Deni ion I.2.1.- Una e ua ion diferen ial de primer orden lineal es una e ua ion del tipo
y0 = a(x)y + b(x);
(L)
donde a; b : I ! R ontinua e I R intervalo.

Conven ion.- Para no res ribir ada vez, utilizaremos (L) para referirnos a una e ua ion diferen ial lineal
(de primer orden).
Deni ion I.2.2.- Se dira que una e ua ion diferen ial lineal de primer orden es homogenea, si b(x) = 0
8x 2 I en (L). Es de ir
y0 = a(x)y:
(LH)
Solu ion de una e ua ion lineal homogenea
Teorema I.2.3.- La solu ion general de la e ua ion lineal homogenea de primer orden
y0 = a(x)y;
on a : I
! R es ontinua, esta dada por
y(x) = CeA(x)
donde A(x) es una primitiva de a; es de ir, A0 (x) = a(x).
(LH)
(I.2.1)
5
Demostra ion.- Primero mostremos que, toda fun ion de la forma CeA(x) es solu ion; en efe to, derivando
obtenemos
0
0

CeA(x) = C eA(x) = CeA(x) A0 (x) = CeA(x)a(x) = a(x) CeA(x) :
Ahora mostremos que toda solu ion '(x) de (LH) es de la forma CeA(x), donde A es una primitiva de a. En
efe to, si ' es solu ion, onsideremos
(x) = '(x) ;
eA(x)
derivemos , lo ual ondu e a

A(x) '(x)eA(x) A0 (x)
0 (x) = '0 (x)eA(x) '(x)(eA(x) )0 = a(x)'(x)e
= 0:
(eA(x))2
(eA(x))2
Como 0(x) = 0 para todo x 2 I , se dedu e que (x) = C .
Remar as
1.- La hipotesis del teorema pre edente que a(x) es ontinua, asegura la existen ia de una primitiva A(x).
En efe to, el primer teorema fundamental del al ulo lo asegura.
2.- Por el primer teorema fundamental del al ulo, la existen ia de A(x) esta asegurada, lo que no signi a
que A(x) pueda expresarse omo la omposi ion de fun iones elementales. Por ejemplo,
A(x) =
Z x
0
et dt
2
es primitiva de ex2 , pero es imposible expresar A(x) omo la omposi ion de fun iones elementales.
Corolario I.2.4.- Las solu iones parti ulares de una e ua ion diferen ial ordinaria lineal homogenea de
primer orden forman un espa io ve torial real de dimension 1.
Demostra ion.- En el urso de Algebra Lineal, se vio que el onjunto de las fun iones f : I ! R,
por I R es un espa io ve torial para:
denotado
(f + g)(x) = f (x) + g(x);

(f )(x) = f (x); 2 R
Para mostrar que la solu ion general de (LH) es un subespa io ve torial, se debe veri ar que existe al menos
una solu ion de (LH), lo que es ierto porque '(x) = 0, 8x 2 I es solu ion. Luego, dadas dos solu iones ' y
de (LH), enton es ualquier ombina ion lineal de las dos solu iones es solu ion de (LH); es de ir
'; solu iones ) ' + solu ion; 8; 2 R:
En efe to,
(' + )0 (x) = '0 (x) + 0(x)

= (a(x)'(x)) + (a(x) (x))
= a(x)('(x) + (x)) = a(x)(' + (x)):
Por el teorema pre edente eA(x) genera la solu ion general o subespa io ve torial de solu iones parti ulares
de (LH), lo que muestra que la dimension es 1.
Corolario I.2.5.- Un problema de Cau hy, para una e ua ion lineal de primer orden homogenea, y = a(x)y
y (x0 ; y0 ) 2 R2 , admite siempre una solu ion uni a.
Demostra ion.- Mostraremos que el problema diferen ial

y0 = a(x)y;
y(x0 ) = y0
tiene solu ion uni a. En efe to, sea A(x) una primitiva de a(x). (I.2.1) indi a que toda solu ion de (LH) es
de la forma
y(x) = CeA(x):
En parti ular, para el problema a valor ini ial se tiene
y0 = CeA(x0 ) ;
obteniendo de manera uni a
y
C = A(0x0 ) :
e
Ejemplos
1.- La solu ion general de y0 = os(x)y, es
y(x) = Cesin x ;
porque sin0 x = os x.
2.- La solu ion general de y0 = ay, on a 2 R, es
y(x) = Ceax :
3.- La solu ion del problema y0 = a(x)y, y(x0 ) = 0; es

y(x) = 0:
Solu ion de una e ua ion lineal (no homogenea)
Consideremos la e ua ion lineal de primer orden

y0 = a(x)y + b(x)
(L)
donde a(x) y b(x) son fun iones ontinuas denidas sobre un intervalo I . Como ya hemos visto el aso (LH),
supongamos que b(x) no es identi amente nula. A la e ua ion (L), le aso iamos
y0 = a(x)y:
(LH)
Proposi ion I.2.6.- Sea
una solu ion parti ular de (L) y ' ualquier solu ion de (LH), enton es
+'
es una solu ion de (L).
Demostra ion.- En efe to
(' + )0 (x) = '0(x) + 0(x)

= a(x)'(x) + (a(x) (x) + b(x))
= a(x)(' + )(x) + b(x):
Proposi ion I.2.7.- Sea una solu ion parti ular de (L), enton es ualquier solu ion de (L) es de la forma
+ ';
donde ' es solu ion de (LH).
Demostra ion.- Sea una solu ion de (L), su iente ver que
es solu ion de (LH). En efe to
( )0 (x) = 0(x) 0(x)
= (a(x) (x) + b(x)) (a(x)(x) + b(x))
= a(x)( (x) (x)):
Por lo tanto, para ono er la solu ion general de (L), es su iente ono er una solu ion parti ular de (L) y
la solu ion general de (LH), algo que ya sabemos ha er. Este resultado lo expresamos, mediante la siguiente
regla memote ni a.
Solu ion general de (L) = Solu ion general de (LH) + Una solu ion parti ular de (L)
Proposi ion I.2.8.- El problema a valor ini ial
y0 = a(x)y + b(x);
y(x0 ) = y0
tiene solu ion uni a.
Demostra ion.- Ejer i io.

Determina ion de una solu ion parti ular de (L)
La determina ion de la solu ion general de una e ua ion linea de primer orden pasa: por la determina ion de
la solu ion general de la e ua ion lineal homogenea aso iada a la e ua ion en uesti'on (algo que ya sabemos
ha er) y por la determina ion de una solu ion parti ular.
Ahora bien, podemos determinar una solu ion parti ular mediante:
i.- Al tanteo. Este pro edimiento es produ to de la experien ia y la genialidad de la persona que lo
utiliza; sin embargo, en algunos asos se puede en ontrar la solu ion, sobre todo, uando b(x) la parte
no homogenea de (L) es un polinomio, una fun ion sin o os y uando es exponen ial.
Ejemplos
4.- Consideremos la e ua ion y0 = y + x2 . En este ejemplo b(x) = x2 y planteamos omo solu ion
parti ular y(x) = + x + x2 un polinomio del mismo grado que b(x). Remplazando en la
e ua ion obtenemos
+ 2 x = + x + x2 + x2 ;
de donde, planteando las respe tivas e ua iones, obtenemos = 1, = 2 y = 2. Por lo
tanto
y(x) = 2 2x x2
es la solu ion parti ular bus ada.
5.- Consideremos la e ua ion y0 = y + os2x. Podemos intentar on una solu ion de la forma y(x) =
os2x + sin 2x. Remplazando obtenemos
2 sin 2x + 2 os2x = os2x + sin 2x + os2x;
resolviendo las e ua iones algebrai as que se dedu en, se obtiene
1 os2x + 2 sin 2x
y(x) =
5
5
omo solu ion parti ular.

6.- Consideremos la e ua ion y0 = 2y + ex, la solu ion parti ular que estamos bus ando, puede ser de
la forma y(x) = ex. Remplazando obtenemos
ex = 2ex + ex;
de donde = 1 y la solu ion parti ular es
y(x) = ex:
ii)
Cuando no es posible adivinar las solu iones parti ulares, se tiene el

metodo de varia ion de onstantes. Sabiendo que la solu ion general de (LH) es y(x) = CeA(x), donde
A(x) es una primitiva de a(x), la idea del metodo es plantear omo solu ion parti ular de (L)
y(x) = (x)eA(x)
(I.2.2)
donde (x) depende de x. Remplazando en la e ua ion (L), se obtiene
0 (x)eA(x) + (x)eA(x) a(x) = a(x) (x)eA(x) + b(x);
de donde
0 (x) = b(x)e A(x):
(I.2.3)
0
(x) se obtiene de (x) por integra ion.
Ejemplo
7.- Consideremos la e ua ion diferen ial
y 0 = y + ex ;
la solu ion general de la e ua ion lineal homogenea aso iada es
y(x) = Cex :
Planteando la solu ion parti ular bus ada omo y(x) = (x)ex , varia ion de onstantes da
0 (x) = exe x = 1;
de donde (x) = x y la solu ion parti ular sera y(x) = xex. La solu ion general
y(x) = Cex + xex :
Varia ion de Constantes.
Remar a Varia ion de onstantes es un metodo que nun a falla,
E ua iones de Tipo Bernouilli
Son e ua iones diferen iales de primer orden de la forma

y0 = a(x)y + b(x)y
on a 6= 0 y a; b fun iones ontinuas.
(B)
Remar as
1.- La expresion y, on 2 R, esta denida por

y = e ln y ;
(I.2.4)
por onsiguiente, tendra sentido uando y > 0. Cuando es entero o ra ional y puede tener un dominio
de deni ion mayor, por lo que debe tratarse de a uerdo al aso.
9
2.- En la e ua ion de tipo Bernouilli, debera suponerse b(x) es no identi amente nulas, porque sino estaramos en el aso de las e ua iones lineales de primer orden ya vistas.
3.- En la e ua ion de tipo Bernouilli, debera tambien suponerse que 6= 0 y 6= 1, aso ontrario estaramos
en el aso linal ya visto.
Con las remar as he has, supongamos que (x) es una solu ion de (B), planteamos
'(x) = ( (x))1 :
(I.2.5)
Derivando ', se obtiene
'0 (x) = (1 )( (x)) 0 (x):
Como es solu ion de (B), se tiene
'0 (x) = (1 )( (x)) (a(x) (x) + b(x)( (x)) )
= (1 )a(x)( (x))1 + (1 )b(x)
= (1 )a(x)'(x) + (1 )b(x):
Hemos mostrado el:
Teorema I.2.9.- Sean ' y
dos fun iones, si
'(x) = ( (x))1 ; 2 R; 6= 0; 1;
enton es
es solu ion de la e ua ion de tipo Bernouilli
y0 = a(x)y + b(x)y ;
(B)
on a; b ontinuas y no identi amente nulas, si y solamente si ' es solu ion de la e ua ion lineal de primer
orden
z 0 = (1 )a(x)z + (1 )b(x):
Regla Memote ni a.- Cada vez que se en uentre una e ua ion de tipo
z = y1 .
Ejemplo
(L)
Bernouilli, introdu ir la variable
8.- Determinemos la solu ion general de la e ua ion diferen ial

y0 = y + y2 :
Planteamos z = y1 2 = 1=y, de donde zy = 1. Derivando, so obtiene z0y + y0z, multipli ando la e ua ion
por z, se tiene
zy0 = zy + zy2;
remplazando, se obtiene
z 0y = 1 + y;
dividiendo por y, obtenemos la e ua ion lineal para z
z 0 = z 1;
la solu ion general para esta ultima e ua ion es
z (x) = Ce x 1;
10
y por onsiguiente la solu ion general para la e ua ion de Bernouilli es

1 ;
y(x) =
Ce x + 1
teniendo uidado que el denominador no se anule.
Corolario I.2.10.- Un problema diferen ial a valor ini ial para una e ua ion de tipo Bernouilli
y0 = a(x)y + b(x)y ;
y(x0 ) = y0
tiene solu ion uni a, si esta existe.
E ua iones Separables
Las e ua iones diferen iales separables, son e ua iones diferen iales de primer orden, que pueden es ribirse
omo
f (x)
;
(S)
y0 =
g(y)
on f y g fun iones ontinuas y g(y) 6= 0.
Antes de determinar y estudiar las solu iones de la e ua ion (S), analizemos f y sobre todo g. Como
f y g son fun iones ontinuas, por el primer teorema fundamental del Cal ulo Integral, ambas admiten
primitivas, que las denotamos por F (x) y G(y) respe tivamente. Por la deni ion de primitiva, se tiene
F 0 (x) = f (x); G0 (y) = g(y):
Por otro lado, omo g es ontinua y por hipotesis g(y) 6= 0, se tiene, sea g(y) > 0, o sea g(y) < 0. Por
onsiguiente, si g(y) > 0, se tiene que G(y) es estri tamente re iente, y por lo tanto G(y) es inye tiva y
ontinua. En el otro aso, si g(y) < 0, se tiene que G(y) es estri tamente de re iente, y por lo tanto G(y) es
inye tiva y ontinua.
Ahora bien, una fun ion que es ontinua e inye tiva, admite una inversa que es ontinua; en nuestro
analisis la inversa de G, la denotaremos por G 1, de donde
G 1 (G(y)) = y; 8y:
Supongamos que '(x) sea una solu ion de (S), se tiene por onsiguiente
g('(x))'0 (x) = f (x);
on g('(x)) 6= 0. La regla de la adena, da
(G('(x))0 = (F (x))0 ;
de donde
G('(x)) = F (x) + C
'(x) = G 1 (F (x) + C ):
A abamos de mostrar el:
Teorema I.2.11.- Sea la e ua ion de primer orden separable
y0 =
f (x)
;
g(y)
(S)
on f; g ontinuas y g (y ) 6= 0 para todo y ; enton es la solu ion general de (S) esta dada por
y(x) = G
(F (x) + C );
(I.2.6)
11
donde F y G son primitivas de f y g respe tivamente y G
denota la fun ion inversa de G.
Cuando se tiene que resolver una e ua ion separable olo ar en un lado de la

e ua ion todo lo que depende de y e y0, en el otro lado lo que depende de x. Luego integrar.
Regla Memote ni a.Ejemplo
9.- Resolvamos la e ua ion

Esta e ua ion es separable y tenemos
integrando, se obtiene
y0 = x(1 + y2 ):
y0
1 + y2 = x;
ar tan y = 21 x2 + C;
1
y = tan( x2 + C ):
2
Corolario I.2.12.- Un problema diferen ial de Cau hy para una e ua ion diferen ial de tipo separable, on
las misma hipotesis del teorema pre edente, admite una sola solu ion.

Remar as
1.- Es posible en las e ua iones separables (S), debilitar la hipotesis sobre g(y) 6= 0. El pro edimiento de
resolu ion sera el mismo, pero puede su eder que la solu ion que se en uentre no ontemple toda la
solu ion general de la e ua ion en uestion, ni tampo o que los problemas a valores ini iales tengan
solu ion uni a. Este aspe to lo veremos en la siguiente se ion.
2.- En general, la fun ion G 1 no puede ser expresada omo la omposi ion de fun iones elementales, por
lo que resulta mas onveniente expresar la solu ion general impl itamente; es de ir, dejar la solu ion
general omo
G(y) = F (x) + C:
E ua iones de tipo Homogeneo
A no onfundir on las e ua iones lineales homogeneas. Una e ua ion diferen ial de tipo homogenea, es una
e ua ion diferen ial de primer orden que puede es ribirse omo
y
y0 = f ( );
x
donde f es una fun ion ontinua y x 6= 0.

Para resolver este tipo de e ua iones, la substitu ion que salta a la vista onsiste en plantear
De donde zx = y y derivando obtenemos
y
z= :
x
y0 = z 0 x + z = f (z ):
Hemos obtenido por onsiguiente una nueva e ua ion diferen ial dada por
f (z ) z
;
z0 =
que es una e ua ion de tipo separable.
Ejemplo
(H)
12
10.- Resolvamos la e ua ion diferen ial
2y2 :
y 2 + x2
Esta e ua ion es de tipo homogeneo, si la es ribimos de la manera siguiente
2(y=x)2 ;
y0 =
1 + (y=x)2
planteando z = y=x, obtenemos la e ua ion diferen ial separable
2z2 z = z 1 + 2z + z2
xz 0 =
1 + z2
1 + z2
2
= z(1z++z1)2 ;
y0 =
de donde
1 + z2 z0 = 1 ;
x
z (1 + z )2
des omponiendo en fra iones simples o par iales obtenemos

2 z0 = 1 ;
1
z (z + 1)2
x
integrando, se tiene
ln z + 1 +2 z = ln x + C:
Remplazamos z = y=x y efe tuamos opera iones algebrai as, obteniendo
ln y ln x + 2 y +x x = ln x + C;
de donde la solu ion general, esta dada por
ln y + 2 y +x x = C:
Otras E ua iones Diferen iales de Primer Orden
Una forma de resolver e ua iones diferen iales de primer orden, onsiste en ex luir los diferentes tipos de
e ua iones; mas pre isamente lo que se ha e es:
i) Veri ar si es una e ua ion diferen ial lineal, si lo es resolverla, sino.
ii) Veri ar si es una e ua ion diferen ial de tipo Bernouilli, si lo es resolverla, sino.
iii) Veri ar si es una e ua ion separable, si lo es resolverla, sino.
iv) Veri ar si es una e ua ion de tipo homogeneo, si lo es resolverla.
Cuando se han ex luido los uatro asos anteriores, lo que se podra ha er es: intentar una substitu ion
o ambio de variable que sea lo su ientemente evidente para que la nueva e ua ion sea mas simple. La
ele ion de la substitu ion orre ta es fruto de la experien ia, la madures y la genialidad de la persona que
esta resolviendo. Esto se onsigue on la pra ti a.
En el aso en que no se pudiese en ontrar una substitu ion o ambio de variable que permita resolver
la e ua ion se puede intentar aproximar la solu ion mediante algun metodo numeri o. Los rudimentos los
veremos en la siguiente se ion.
A manera de ilustra ion onsideremos los siguientes ejemplos.
Ejemplos
13
I.3 Problemas de Existen ia y Uni idad
11.- La e ua ion diferen ial
y0 = sin(x + y)
no pertene e a ninguno de los tipos de e ua iones estudiadas mas arriba, sin embargo, la substitu ion
que se ve fa ilmente es
z = x + y;
on lo que la e ua ion diferen ial se onvierte en
z 0 = sin z + 1
que es una e ua ion de tipo separable.
12.- En general las e ua iones que se es riben bajo la forma
y0 = f (ax + by + );
on a; b; 2 R y a; b 6= 0, on la substitu ion z = ax + by + , se onvierte en una e ua ion de tipo
separable. >Cual es la nueva e ua ion?
13.- E ua iones diferen iales de primer orden, de la forma

ax + by + r
0
;
y =f
x + dy + s
on ad b 6= 0 y (r; s) 6= 0. Este tipo de e ua iones pueden onvertirse en e ua iones de tipo homogeneo,
mediante las substitu iones
u = x x0
v = y y0
donde (x0 ; y0) es la interse ion de las re tas:
ax + by + r = 0
x + dy + s = 0:
Ejer i io.- Determinar la e ua ion de tipo homogeneo que resulta de estas substitu iones.
La resolu ion de e ua iones diferen iales no tienen mayor utilidad, si no van a ompa~nadas on la solu ion
de un problema diferen ial. En esta se ion veremos la importan ia de las ondi iones o hipotesis que van
on los problemas diferen iales para evitar de ometer errores absurdos.
Comenzemos ilustrando on el siguiente ejemplo, onsideremos el problema diferen ial a valor ini ial
p
y0 = 3 y2
y(a) = 0;
(I.3.1)
p
donde a 2 R. Para ha er un uso mas e iente de exponentes, es ribamos y2=3 en lugar de 3 y2.
Re ordemos en la se ion pre edente, on las hipotesis dadas para ada problema diferen ial a valor
ini ial, la solu ion era uni a. Para el problema planteado deberiamos esperar que la solu ion tambien sea
uni a.
14
Resolvamos (I:3:1), la lasi a ion de e ua iones diferen iales, nos muestra de manera evidente, que
esta es una e ua ion de tipo separable, la resolvemos ingenuamente, sin jarnos en las hipotesis que debe
umplir di ha e ua ion. Se tiene por lo tanto
y0
= 1;
y2=3
3y1=3 = x + C;
1
y = (x + C )3 ;
27
ahora introduz amos la ondi ion ini ial, lo que da C = a, por lo que la solu ion en ontrada es
1
(I.3.2)
y(x) = (x a)3 :
27
La gra a de esta solu ion puede apre iarse en la gura I.2.
Por otro lado, si no hubieramos resuelto ingenuamente, una veri a ion sen illa muestra que
y(x) = 0;
(I.3.3)
es tambien solu ion del problema (I.3.1), on lo que eviden iamos que el problema no tiene solu ion uni a.
Mostraremos a ontinua ion que di ho problema
tiene una innidad de solu iones, en efe to, para b a
8
< = 0 si x b
y(x) = 1
3
:
27 (x b) si x b
es tambien solu ion de (I.3.1). Para a,
8
< 1 si x
y(x) = 27
:
0 si x
es tambien solu ion de (I.3.1). Por ultimo, para8 a b,
1 si x
>
>
>
>
< 27
y(x) = 0 si x b
>
>
>
1 (x b)3 si x b
>
:
27
es tambien solu ion del problema (I.3.1). Estas solu iones pueden apre iarse en la gura I.2.
15
Figura I.2.- Solu iones del problema (I.3.1)
El problema (I.3.1) a primera vista era un problema ompletamente sen illo, sin embargo a abamos de
mostrar que tiene una innidad de solu iones, sin poder elegir aquella que pueda onvenirnos para satisfa er
nuestros requerimientos de solu ion del problema.
Condi iones su ientes para la existen ia y uni idad de solu iones
Deni ion I.3.1.- Se di e que una fun ion f : I U ! R ontinua, on I; U intervalos satisfa e una
ondi ion de Lips hitz si para todo subintervalo errado y a otado J de I , existe una onstante CJ > 0 tal
que 8x 2 J y 8y; z 2 U , se tiene
jf (x; z ) f (x; y)j CJ jz yj :
(I.3.4)
Ejemplos
1.- Consideremos
f (x; y) = a(x)y + b(x);
donde a; b son fun iones ontinuas. f satisfa e una ondi ion de Lips hitz por que
jf (x; z ) f (x; y)j = ja(x)j jz yj max
ja(x)j jy z j :
x2J
2.- Consideremos
f (x; y) = sin y
satisfa e una ondi ion de Lips htiz; en efe to, apli ando el teorema del valor medio, se tiene
jsin y sin z j = j os j jy z j 1 jy z j :
Teorema I.3.2.- Supongamos que f : I U
y0 2 U , enton es el problema a valor ini ial
! R ontinua satisfa e una ondi ion de Lips htiz y x0 2 I ,

y0 = f (x; y);
y(x0 ) = y0
admite una uni a solu ion.
Demostra ion.- Aunque no haremos una demostra ion formal y rigurosa del teorema, daremos un bosquejo
de lo que orresponde a la existen ia de la solu ion. Para tal efe to, deniremos de manera re ursiva la
su esion de fun iones
'0 (x) = y0 ;
Z x
'n+1 (x) = y0 + f (s; 'n (s)) ds:
x0
Se muestra que esta su esion de fun iones es onvergente y su lmite es una fun ion '(x) que es solu ion
del problema a valores ini iales planteado. Lo importante de esta demostra ion y lo que hay que retener es
la su esion de fun iones denidas mas arriba. Esta demostra ion y la onstru ion de la su esion se debe a
Pi ard, matemati o fran es del siglo XIX.
Ejemplo
3.- La fun ion f : R R ! R denida por f (x; y) = y2=3 si bien es ontinua, no satisfa e una ondi ion de
Lips htiz, porque sino el problema a valor ini ial (I.3.1) tendra solu ion uni a.
Veri ar que una fun ion satisfa e una ondi ion de Lips htiz, puede resultar a menudo una tarea
ompli ada. Tenemos un resultado mas manipulable para determinar uando un problema a valor ini ial
16
Corolario I.3.3.- Si f (x; y) es una fun ion ontinua y ontinuamente derivable respe to a y, enton es el
problema a valor ini ial
y0 = f (x; y)
y(x0 ) = yo
admite solu ion uni a.

Corolario I.3.4.- Si f (x; y) es una fun ion ontinua y ontinuamente derivable respe to a y, enton es la
solu ion general de
y0 = f (x; y)
forma una familia uniparametri a de urvas, es de ir la solu ion general puede es ribirse omo
donde ha e el papel de onstante.
F (x; y; ) = 0;
Veri ar que las solu iones generales de las e ua iones diferen iales de primer orden de tipo
lineal, tipo Bernouilli, de tipo separable y de tipo homogeneo pueden expresarse omo
F (x; y; ) = 0:
Ejer i io.-
Solu iones Aproximadas de una E ua ion Diferen ial
En general, es imposible expresar una solu ion parti ular o la solu ion general de una e ua ion diferen ial
omo omposi ion de fun iones elementales.
Por otro lado, en la pra ti a solo se resuelven problemas diferen iales, y la solu ion que a ompa~na el
problema da informa iones ualitativas y en algunos asos uantitativas.
Las e ua iones diferen iales son parte de los modelos matemati os que des riben un fenomeno, por lo
que la des rip ion de ierta manera es par ial, sujeta a errores de aproxima ion, debido a las simpli a iones
que se ha e on el modelo. En onse uen ia es inutil tratar de determinar una solu ion de manera exa ta y
es justi able por la misma razon tratar on aproxima iones de las solu iones.
Itera iones de Pi ard
Ya nos referimos a las itera iones de Pi ard en el bosquejo de demostra ion del teorema I.3.2. Este metodo
esta formulado para resolver problemas a valores ini iales y no en ontrar la solu ion general de una e ua ion
diferen ial.
Consideremos el problema a valor ini ial
y0 = f (x; y);
y(x0 ) = y0
donde f es una fun ion ontinua que satisfa e una ondi ion de Lips htiz. Las itera iones de Pi ard estan
dadas por:
'0 (x) = y0 ;
Z x
'n+1 (x) = y0 + f (s; 'n (s)) ds; k = 0; 1; : : : :
x0
Cuando f satisfa e una ondi ion de Lips hitz, la su esion de fun iones 'n, tiende a una fun ion lmite '(x).
Por lo tanto, se tiene
Z x
0
' (x) = y0 + f (s; '(s)) ds;
x0
se ve inmediatamente que '(x0 ) = y0 y derivando se veri a que ' es solu ion de la e ua ion diferen ial.
Es obvio que utilizando el metod de Pi ard, no se llegara a la solu ion exa ta, porque es humanamente
imposible determinar la totalidad de las fun iones de la su esion. Por onsiguiente, nos ontentamos on
determinar un numero nito de estas fun iones.
Por otro lado debemos remar ar que si estamos utilizando un metodo para aproximar la solu ion del
problema planteado, es porque en ontrar la solu ion de este es sumamente ompli ada, por no de ir imposible.
Ahora bien, el al ulo de las integrales puede ser tambien tan ompli ado omo la resolu ion de la e ua ion
por lo que este metodo tiene mas un ara ter teori o que pra ti o.
17
Ejemplo
4.- En ontremos una aproxima ion de la solu ion del problema a valor ini ial
y0 = y;
y(0) = 1:
En este ejemplo f (x; y) = y, por lo que las itera iones de Pi ard estan dadas por:
'0 (x) = 1;
Z x
'1 (x) = 1 + 1 ds = 1 + x;
0
Z x
1
'2 (x) = 1 + (1 + s) ds = 1 + x + x2 ;
2
0
Z x
1
1
1
'3 (x) = 1 + (1 + s + s2 ) ds = 1 + x + x2 + x3
2
2
3!
0
...
1
1
1
'n (x) = 1 + x + x2 + x3 + + xn :
2 3!
n!
Polgonos de Euler
La utiliza ion del metodo de Pi ard para resolver problemas a valor ini ial en general solo tiene un interes
teori o, por lo que para en ontrar aproxima iones de las solu iones de los problemas a valor ini ial se preere
otras alternativas.
Consideremos la e ua ion diferen ial de primer orden
y0 = f (x; y);
donde f es ontinua. Si '(x) es una solu ion de esta e ua ion, se tiene
'0 (x) = f (x; '(x)):
Por otro lado, se tiene
'(x + h) '(x)
;
'0 (x) = lim
h!0
h
es rito de otra manera, se tiene
'(x + h) = '(x) + h'0 (x) + r(x; h);
donde
r(x; h)
= 0:
lim
h! h
Si despre iamos r(x; h), la expresion
g(h) = '(x) + h'0 (x)
representa la e ua ion parametri a de la re ta tangente al grafo de ' en el punto x.
Por otro lado, '0 (x) = f (x; '(x)), por lo que la e ua ion parametri a de la re ta tangente se onvierte
en
g(h) = '(x) + hf (x; '(x)):
Ahora bien, la re ta tangente de una urva en punto dado es la re ta que mejor aproxima el omportamiento
de la urva. Por onsiguiente si h es bastante peque~no es razonable aproximar '(x + h) pro g(h).
18
De esta manera, podemos denir un metodo de aproxima ion para resolver el problema diferen ial
siguiente
y0 = f (x; y);
y(x0 ) = y0
y queremos determinar y(xf ). Para simpli ar el problema suponemos que xf > x0 . Dividimos el intervalo
[x0 ; xf en
x0 < x 1 < x 2 < < x n = x f ;
planteamos hk = xk xk 1 para k = 1; : : : ; n. Denimos
yk = yk 1 + hk f (xk 1 ; yk 1 ); k = 1; : : : ; n;
y tomamos omo aproxima ion de y(xn) yn.
Ejemplo
5.- Consideremos el problema diferen ial
y0 = y y2 + x;
y(0) = 0;
intervalo [0; 1
y determinemos y(1). Dividiendo el

de manera equidistante on h = 1=10, (n=10), se
obtiene el resultado que puede observarse en la gura I.3, al igual que la solu ion exa ta. La solu ion
obtenida por el metodo de Euler se la remar a por la forma poligonal de su grafo.
Figura I.3.- Ilustra ion del metodo de Euler
19
El metodo de Euler debe utilizarse para obtener aproxima iones de solu iones que nos den mas informa ion ualitativa que uantitativa. Para obtener aproxima iones de mayor pre ision, no es a onsejable
utilizar el metodo de Euler.
Iso linas
El lado dere ho de la e ua ion diferen ial
y0 = f (x; y)
indu e un ampo de pendientes sobre R2 . Remar amos que en el punto (x; y), f (x; y) es la pendiente de la
re ta tangente a la o las solu iones que pasan por (x; y).
Una iso lina es una urva de igual pendiente; es de ir es una urva de e ua ion
f (x; y) = m:
donde m es la pendiente.
La idea del metodo de las iso linas onsiste en trazar iso linas para diferentes valores de pendientes y
luego dibujar las solu iones utilizando las pendientes de ada una de las iso linas.
Ejemplo
6.- Determinemos algunas de las solu iones utilizando iso linas de la e ua ion diferen ial
y0 = y y2 + x:
Las iso linas estan dadas por las e ua iones
x = (y 1=2)2 + m 1=4
En la gura I.4 se tiene las iso linas, on el respe tiva pendiente que representan, en lineas segmentadas.
Los grafos de las solu iones son las lineas ontinuas.
20
Figura I.4.- Ilustra ion del metodo de las Iso linas
I.4 E ua iones de Orden Superior
En las dos se iones pre edentes vimos e ua iones diferen iales y problemas diferen iales rela ionados a
e ua iones de primer orden. En esta se ion estudiaremos las e ua iones de orden superior mas orrientes.
Redu ion del Orden
Consideremos la e ua ion diferen ial de orden n, en forma expl ita,

y(n) = f (x; y; y0 ; : : : ; y(n 1) );
on f ontinua.
Partimos del supuesto que mientras el orden de la e ua ion diferen ial es mayor, la resolu ion es mas
di il. Por onsiguiente, si se puede onvertir la e ua ion a resolver en una de orden mas peque~no, la
resolu ion se habra \simpli ado".
21
La redu ion de orden es posible en los siguientes dos asos:

1.- La e ua ion diferen ial de orden n es de la forma
y(n) = f (x; y0 ; y00 ; ; y(n 1));
(I.4.1)
es de ir, en la e ua ion no interviene la variable y de manera expl ita. Planteando
z (x) = y0 (x);
la e ua ion (I.4.1), se onvierte en la e ua ion
z (n 1) = f (x; z; z 0; : : : ; z (n 2))
(I.4.2)
que es una e ua ion de orden n 1.
2.- La e ua ion diferen ial de orden n es de la forma
y(n) = f (y; y0 ; : : : ; y(n 1) );
(I.4.3)
es de ir, en la e ua ion diferen ial no interviene expl itamente la variable x. Planteamos
y0 = u(y):
A diferen ia del primer aso, en el segundo y0 se expresa en fun ion de y. El siguiente paso es expresar
las otras derivadas de y, en fun ion de y, para lo que derivamos utilizando la regla de la adena.
du dy
y00 = (y0 )0 =
dy dx
= u du
:
dy
Pasamos a la siguiente derivada.

y000 = (y00 )0 =
d du du
(u )
dy dy dy
d u du 3
+ ( dy ) :
= u du
dy dy2
2
Observamos que tanto y00, y000 han bajado de un orden uando se expresa en fun ion de y. Continuando
on el mismo pro edimiento se llega a obtener la e ua ion
u(n 1) = g(y; u; : : : ; u(n 2);
donde las derivadas de u son respe to a y.
22
Ejemplos
y00 = os x:
Mediante la substitu ion z = y0, se obtiene la e ua ion diferen ial de primer orden
z 0 = os x;
uya solu ion general es
z = sin x + C:
Por onsiguiente y satisfa e la e ua ion diferen ial
y0 = sin x + C;
de donde
y(x) = os x + Cx + D:
2.- Resolvamos la e ua ion diferen ial de segundo orden
y00 + y = 0;
en esta e ua ion x no interviene expli tamente, por lo que planteando y0 = u(y), da la e ua ion
diferen ial de primer orden
u0u + y = 0;
que es una e ua ion de tipo separable, la solu ion de esta e ua ion es
u2 + y2 = C:
Suponiendo y 0 y C 0, se tiene
p
y0 = C 2 y2 ;
e ua ion de tipo separable, uya solu ion es
ar sin(y=C ) = x + D;
es de ir
y = C sin(x + D);
que se onvierte en
y(x) = C1 os x + C2 sin x:
E ua iones Diferen iales Lineales de Orden Superior
Una e ua ion diferen ial lineal de orden n es una e ua ion que puede es ribirse de la forma
y(n) + an 1 (x)y(n 1) + + a1 (x)y0 + a0 (x)y = b(x)
on a0; : : : ; an 1 y b fun iones ontinuas sobre un intervalo.
(L);
E ua iones Lineales Homogeneas

Se dira que la e ua ion (L) de orden n, es homogenea si b(x) es identi amente nula, es de ir si la e ua ion
puede expresarse omo
(x)y(n 1) + + a1(x)y0 + a0(x)y = 0:

(LH)
El estudio de las e ua iones lineales homogeneas de orden n sigue el mismo amino que el seguido para las
e ua iones lineales de primer orden.
y(n) + an
23
Proposi ion I.4.1.- La solu ion general de una e ua ion lineal homogenea de orden n, es un subespa io
ve torial de dimension n.
Demostra ion.-Al igual que en aso de primer orden, puede veri arse que la solu ion general de (LH) es
un subespa io ve torial. Probar que la dimension es n, es apa de los objetivos del presente urso.
Al ser la solu ion general de (LH) un espa io ve torial real de dimension n, para determinar ualquier
solu ion de (LH) es su iente ono er n solu iones linealmente independientes.
Deni ion I.4.2.- Se llama sistema fundamental (SF) de la e ua ion diferen ial lineal homogenea de
orden n
y(n) + an 1 (x)y(n 1) + + a1 (x)y0 + a0 (x)y = 0
(LH)
a un onjunto f'1; '2 ; : : : ; 'ng de solu iones de (LH) que sean linealmente independientes.
y(n) + an 1 (x)y(n
y(x0 ) = y0 ;
y0 (x0 ) = y1 ;
..
.
y(n
1)
(x0 ) = yn
1)
+ + a1(x)y0 + a0(x)y = 0;
on x0 en el intervalo de deni ion de las ai , tiene solu ion uni a.
Demostra ion.- Sea '1 ; : : : ; 'n
un sistema fundamental de la e ua ion diferen ial (LH). Por onsiguiente,

ualquier solu ion de (LH) se es ribe de manera uni a omo
y(x) = 1 '1 (x) + 2 '2 (x) + + n 'n (x):
(I.4.4)
Por onsiguiente los i satisfa en el sistema de e ua iones
1 '1 (x0 ) + 2 '2 (x0 ) + +
n 'n (x0 )
= y0
1 '01 (x0 ) + 2 '02 (x0 ) + +
n '0n (x0 )
= y1
...
...
(n 1)
(n 1)
(n 1)
1 '1 (x0 ) + 2 '2 (x0 ) + + n 1 !'n (x0 ) = yn 1
que es rito de manera matri ial
0
'1 (x0 )
'2 (x0 ) 'n (x0 ) 1 0 1 1 0 y0 1
0
'02 (x0 ) '0n (x0 ) C B 2 C B y1 C
B '1 (x0 )
B
B C B C
...
... C
A .. A = .. A;

.
.
(n 1)
(n 1)
(n 1)

y
(x0 ) '2 ({zx0 ) 'n (x0 ) } | {zn } | {zn }
'
| 1
C
Y0
W (x0 )
es de ir
W (x0 )C = Y0 :
((I.4.5)
Ahora bien W (x0 ) es inversible, sino f'1; : : : ; 'ng no sera un sistema fundamental de la e ua ion. Por lo
tanto
C = (W (x0 )) 1 Y0
da la existen ia y la uni idad.
24
Deni ion I.4.4.- Se llama a
0
B
W (x) = B

'1 (x)
'01 (x)
...
(n 1)

'n (x)
'0n (x)
...
1)
1
C
C
A
(x) 'n(n (x)

la matriz wronskiana del sistema f'1; : : : ; 'ng y det W (x) el wronskiano del sistema f'1; : : : ; 'ng.
'1
(x)
'2 (x)
'02 (x)
'2(n
1)
Corolario I.4.5.- Sea f'1 ; : : : ; 'n g n solu iones de la e ua ion (LH) de orden n
y(n) + an 1 (x)y(n 1) + + a1 (x)y0 + a0 (x)y = 0;
enton es f'1 ; : : : ; 'n g es un sistema fundamental de (LH), si y solamente si W (x) es inversible para todo x,
si y solamente si det W (x) 6= 0 para todo x.
Demostra ion.- Conse uen ia dire ta del Algebra Lineal.

Ejemplo
.
3.- Consideremos la e ua ion diferen ial lineal homogenea de segundo orden

y00 + y = 0;
por ejemplo 2, se tiene que os x y sin x son solu iones de la e ua ion en uestion. El wronskiano esta
dado por

os
x
sin
x
det W (x) = det sinx os x = 1;
de donde os x y sin x es un sistema fundamental, por lo que la solu ion general esta dada por
y(x) = 1 os x + 2 sin x:
Determina ion de Sistemas Fundamentales
Habiendo visto que para determinar la solu ion general de una e ua ion diferen ial lineal homogenea es
su iente ono er un sistema fundamental de solu iones, de la misma manera sabiendo re ono er si un
onjunto de solu iones parti ulares de una e ua ion (LH) es un sistema fundamental, el siguiente paso es
ono er los metodos que permitan en ontrar tales sistemas.
Lastimosamente no existe un metodo general para poder determinar sistemas fundamentales de una
e ua ion lineal homogenea. Sin embargo podemos resaltar las siguientes situa iones en las que se puede
determinar un sistema fundamental:
1.- Si la e ua ion es de segundo orden y se ono e una solu ion no nula.
Sea
y00 + p(x)y0 + q(x)y = 0;
y '(x) una solu ion no nula. Suponemos que la otra solu ion del sistema fundamental a determinar es
de la forma
y(x) = (x)'(x);
donde (x) es una fun ion no onstante, aso ontrario tendramos dos solu iones linealmente dependientes. Derivando, se obtiene
y0 (x) = 0 (x)'(x) + (x)'0 (x);
y00 (x) = 00 (x)'(x) + 2 0(x)'0 (x) + (x)'00 (x);
remplazando en la e ua ion diferen ial, se tiene
'(x) 00 (x) + (2'0 (x) + '(x)p(x)) 0 (x) + (x)('00 (x) + p(x)'0 (x) + q(x)'(x)) = 0;
{z
}
|
=0
25
de donde 0(x) es solu ion de la e ua ion de primer orden lineal homogenea, que ya sabemos resolver,
2'0 (x) + '(x)p(x) z(x) = 0;
z 0 (x) +
'(x)
por onsiguiente
0 (x) = eA(x) ;
0
donde A(x) es una primitiva de 2' (x) '+(x')(x)p(x) , por lo tanto
(x) =
Z x
x0
eA(s) ds:
Ejemplo
4.- Consideremos la e ua ion de segundo orden

x2 y00 xy0
3y = 0;
se veri a que y(x) = x3 es una solu ion parti ular no nula de la e ua ion diferen ial de segundo
orden en uestion. Determinemos otra solu ion linealmente independiente. Planteamos
y(x) = (x)x3 ;
derivando y remplazando en lae ua ion diferen ial, se obtiene
x5 00 + (6x4 x4 ) 0 (x) = 0;
es de ir 0(x) es solu ion de la e ua ion diferen ial
5
z 0 = z;
por lo que 0(x) = 1=x5, de donde
= 41 x 4;
La solu ion de la e ua ion diferen ial, esta dada por

y(x) = 1 x3 + 2 :
x
(x) =
2.- La e ua ion diferen ial lineal es a oe ientes onstantes.

Es de ir, la e ua ion (LH), tiene la forma
y(n) + an 1 y(n 1) + + a1 y0 + a0 y = 0:
Para poder en ontrar un sistema fundamental para este tipo de e ua iones, hagamos un peque~no repaso
de variable ompleja.
Un po o de Variable Compleja
Re ordemos que C el plano omplejo es R2 provisto de una adi ion y una multipli a ion que ha en de
C un uerpo onmutativo. Un elemento z 2 C puede es ribirse de la manera siguiente
z = x + iy on x; y 2 R:
Una fun ion f : I R ! C admite la representa ion siguiente
f (x) = u(x) + iv(x);
26
donde u; v : I ! R.
Ahora bien, si u y v son derivables, f es derivable y
f 0 (x) = u0 (x) + v0 (x):
Proposi ion I.4.6.- Sean f; g : I R ! C , a 2 C , enton es, se tiene:
i) (f + g )0 (x) = f 0 (x) + g 0 (x).
ii) (af )0 (x) = af 0 (x).
iii) (f (x)g (x))0 = f 0 (x)g (x) + f (x)g 0 (x).
Remar a.-Las reglas de al ulo para derivadas son las mismas que para las fun iones reales.
Deni ion I.4.7.- Sea a = + i 2 C , x 2 R, se dene
eax = ex( os x + i sin x):
Proposi ion I.4.8.- Se tiene:
i) Dados a; b 2 C , se veri a
ii) Para a 2 C , se veri a
(I.4.6)
(I.4.7)
(I.4.8)
(I.4.9)
e(a+b)x = eax ebx:
(I.4.10)
(eax)0 = aeax:
(I.4.11)
Demostra ion.- El punto i) ejer i io. Demostremos el punto ii). Se tiene

(eax)0 = ex( os x + i sin x)0
= ex( os x + i sin x) + ex( sin x + i os x)
= ex( os x + i sin x) + iex( os x + i sin x)
= ( + i)ex( os x + i sin x) = aeax:
Re onsideremos nuevamente la e ua ion diferen ial lineal homogenea de orden n a oe ientes onstantes
y(n) + an 1 y(n 1) + + a1 y0 + a0 y = 0:
y planteamos
L(y) = y(n) + an 1 y(n 1) + + a1 y0 + a0 y:
(LH)
Remar amos que L es una apli a ion lineal y resolver la e ua ion (LH), es en ontrar y(x) tal que
L(y) = 0.
Ahora bien, en lugar de onsiderar solamente solu iones reales, podemos prolongar nuestras solu iones
a solu iones y : R ! C .
Utilizando las reglas de al ulo para derivadas, se tiene para r 2 C
L(erx) = (erx)(n) + an 1 (erx)(n 1) + + a1 (erx)0 + a0 erx
= rn erx + an 1rn 1 erx + + a1rerx + a0 erx
= erx(rn + an 1 rn 1 + + a1r + a0)
= erxl(r):
Deni ion I.4.9.- El polinomio l(r) = rn + an 1rn
+ + a1r + a0 es el polinomio ara tersti o

(PC) de la e ua ion diferen ial lineal homogenea a oe ientes onstantes de orden n
y(n) + an 1 y(n 1) + + a1 y0 + a0 y = 0:
(LHC)
1
27
A abamos de demostrar el:

Teorema I.4.10.- y(x) = erx es solu ion de la e ua ion (LHC) si y solamente si r es una raiz de l().
Remar a.-La e ua ion (LHC) admite al menos una solu ion de la forma y(x) = erx debido al Teorema
Fundamental del Algebra que asegura de al menos una raiz r 2 C del polinomio ara tersti o de (LHC).
El siguiente paso es estudiar las solu iones de la forma erx, teniendo en uenta que nuestro objetivo es
determinar un sistema fundamental
de solu iones \reales" de (LHC). Sea r 2 R, se tiene dos asos:
1) r 2 R. En este aso erx es una fun ion real, por lo que y(x) = erx es una solu ion \real" de (LHC).
2) r = + i on 6= 0. Como el (PC) es a oe ientes reales r = i es tambien raiz del (PC),
por lo que obtenemos las siguientes solu iones \ omplejas"
'(x) = ex ( os x + i sin x)
(x) = ex( os x i sin x):
Como la solu ion general es un subespa io ve torial, se tiene que
1
y1 (x) = ('(x) + (x)) = ex os x;
2
1
y2 (x) = ('(x)
(x)) = ex sin x
2i
son solu iones de (LHC), pero y1(x) y y2(x) son solu iones \reales" de (LHC).
Ejer i io.- Mostrar que y1 (x) y y2 (x) son linealmente independientes. (Utilizar el wronskiano).
En resumen, para ada raiz r 62 R obtenemos un par de solu iones linealmente independientes.
Con lo desarrollado mas arriba estamos en la posibilidad de enun iar nuestro primer resultado.
Teorema I.4.11.- Si el polinomio ara tersti o de la e ua ion diferen ial (LHC) de orden n, tiene sus
n rai es diferentes, enton es se tiene un sistema fundamental para (LHC). Ademas si r 2 R es raiz del
(PC), erx ha e parte del sistema fundamental; si r = + i on 6= 0 ex os x y ex sin x ha en
parte del sistema fundamental.
Demostra ion.- Efe tivamente ada raiz aporta on una
o dos solu iones para onformar el sistema

fundamental, si la raiz es real se tiene una solu ion; si la raiz no es real, esta y su onjugada aportan on
dos solu iones. En resumen se tiene n solu iones no nulas. Utilizando el wronskiano se puede veri ar
que son linealmente independientes.
Ejemplo
5.- Consideremos la e ua ion diferen ial de orden 4.

y(4) y = 0:
El polinomio ara tersti o esta dado por

l() = 4 1 = (2 + 1)( + 1)( 1):
Se dedu e a partir de la fa toriza ion que las rai es son:
1; 1; i; i;
todas diferentes, por lo que el sistema fundamental estara dado por
ex ; e x ; os x; sin x:
La solu ion general de esta e ua ion es
y(x) = 1 ex + 2 e x + 3 os x + 4 sin x:
28
Polinomio ara tersti o on rai es multiples
El teorema pre edente no es apli able uando el polinomio ara tersti o de la e ua ion (LHC) tiene
rai es multiples o multipli idad mayor a 1, ya que faltan solu iones para ompletar un sistema fundamental.
Introdu iendo los operadores o apli a iones lineales
D(y) = y0
I (y) = y:
y utilizando la onven ion Dk (y) = D(Dk 1 (y)) = y(k), el operador L dado mas arriba, se puede es ribir
omo
L = Dn + an 1 Dn 1 + + a1 D + a0 I:
Como l() el polinomio ara tersti o de la e ua ion (LHC), puede fa torizarse en C , sin importar el
orden de los fa tores omo
l() = ( r1 )n1 ( r2 )n2 ( rm )nk ;
donde r1 ; : : : ; rm son las rai es diferentes de l() y los nj son las multipli idades de ada una de las
rai es rj . Observamos que nj 1 y n1 + n2 + + nm = n el orden de la e ua ion (LHC).
La fa toriza ion de l, permite fa torizar por lo tanto L de la manera siguiente
L = (D r1 I )n1 (D r2 I )n2 (D rm )nm :
Proposi ion I.4.12.- Si '(x) es solu ion de (D rj I )nj , donde rj es una raiz del polinomio ara tersti o
de (LHC) y nj es la multipli idad de la raiz, enton es ' tambien es solu ion de (LHC ).
Demostra ion.- Como el orden de los fa tores no importa, L puede es ribirse omo
Y
L = ( (D ri I )ni )(D rj I )nj ;
de donde
i6=j
L('(x)) = (
i6=j
(D
= ( (D
i6=j
=0
ri I )ni )((D rj I )nj ('(x)))

ri I )ni )(0)
Ahora veamos las solu iones on las que aporta (D
rj I )nj
al sistema fundamental.
Proposi ion I.4.13.- Son solu iones linealmente independientes de (D rj I )nj (y) = 0 las fun iones
erx; xerx; : : : ; xnj 1 erx:
Demostra ion.- Por indu ion sobre nj . Para nj = 1, er1 x es solu ion de la e ua ion, la veri a ion
es trivial. Suponemos ierto para nj 1 on nj > 1. Por hipotesis de indu ion, son solu iones
erj ; : : : ; xnj 2 erj x , porque
(D rj I )nj = (D rj I )(D rj I )nj 1:
Solo debemos veri ar que xnj 1erj x es solu ion, en efe to
(D
rj I )nj (xnj 1 erj x ) = (D rj I )nj

= (D rj I )nj
= (D rj I )nj
=0
(D rj )(xnj 1erj x)
(D rj )((nj 1)xnj
1
((nj 1)xnj 2erj x)
1
erj x + rj xnj 1 erj x rj xnj 1 erj x )
29
Remar a.- Si rj es una raiz no real de multipli idad nj , enton es rj es una raiz de multipli idad nj y
ontribuyen ambas rai es on 2nj solu iones al (SF), las uales son
ex os x; xex os x; : : : ; xnj 1 ex os x
ex sin x; xex sin x; : : : ; xnj 1 ex sin x
Finalmente enun iamos:
Teorema I.4.14.- Sean r1 ; : : : ; rm las rai es diferentes de una e ua ion (LHC) de orden n y la fa toriza ion del (PC) esta dada por
l() = ( r1 )n1 ( r2 )n2 ( rm )nk ;
enton es las ontribu iones de solu iones para el sistema fundamental de (LHC) esta dada de la manera
siguiente, si rj 2 R, rj ontribuye on
erx; xerx; : : : ; xnj 1 erx;

si rj
= + i on 6= 0, rj y rj ontribuyen on
ex os x; xex os x; : : : ; xnj
ex sin x; xex sin x; : : : ; xnj
ex os x
1 x
e sin x:
1
Demostra ion.- Un peque~no ejer i io de onteo mostrara que hay exa tamente n solu iones.
skiano veri ara la independen ia lineal.
Ejemplo
El wron
6.- Consideremos la e ua ion diferen ial (LHC)

y(6) 2y(5) + 3y(4) 4y(3) + 3y00 2y0 + y = 0;
el polinomio ara tersti o de la e ua ion esta dado por
l() = 6 25 + 34 43 + 32 2 + 1 = ( 1)2 (2 + 1)2 :
Dedu imos que 1, i y i son rai es del polinomio ara tersti o, las tres de multipli idad 2, por lo
que el sistema fundamental bus ado esta dado por las solu iones:
ex; xex ; os x; sin x; x os x; x sin x:
La solu ion general, sera por onsiguiente
y(x) = 1 ex + 2 xex + 3 os x + 4 sin x + 5 x os x + 6 x sin x:
3.- Si laa e ua ion lineal homogenea es de la forma

xx yn + an 1 xn 1 yn 1 + + a1 xy0 + a0 y = 0;
la substitu ion x = et, onvierte esta e ua ion en una (LHC). Ver ejer i ios.
4.- Utiliza ion de Series Enteras.- Re omendable para determinar la solu ion que satisfage un problema
a valor ini ial. Re ordemos los on eptos ne esarios para desarrollar este pro edimiento.
Deni ion I.4.15.- Sea f : I R una fun ion indenidamente derivable, se di e que f es analiti a en
el punto x0 , si
1
X
f (x) =
n (x x0 )n ; jx x0 j < ;
n=0
30
on 0 < +1.
Damos algunas de las propiedades, sin demostra ion de las fun iones analiti as.
Proposi ion I.4.16.- Supongamos que f; g son analiti as en el punto x0 , enton es

1
X
i)
(f + g)(x) = ( n + dn)(x x0 )n;
n=0
f 0 (x) =
ii)
iv)
n=1
n n (x x0 )n 1 ;
(x0 ) = n! n;
n
1 X
X
k dn
f (x)g(x) =
1
X
1
X
iii)
1
X
(n)
n=0 k=0
donde
f (x) =
n=0
n (x x0 )
g(x) =
n;
n=0
(x
x 0 )n ;
dn (x x0 )n :
Consideremos la e ua ion lineal homogenea de orden n

y(n) + an 1 (x)y(n 1) + + a1 (x)y0 + a0 (x) = 0:
(LH)
x0 es un punto ordinario, si ada una de las fun iones ak (x) es analiti a en x = x0 , aso ontrario se
dira que x0 es un punto singular. Enun iamos sin demostra ion el:
Teorema I.4.17.- Si x0 es un punto ordinario de la e ua ion (LH), enton es ualquier solu ion es
analiti a en x = x0 .
La determina ion de un sistema fundamental para la e ua ion (LH), se la ha e eligiendo un punto x0
que sea ordinario y se onsidera los n problemas a valor ini ial
y(n) + an 1 (x)y(n 1) + + a1 (x)y0 + a0 (x) = 0;
y(k) (x0 ) = kj ; k = 1; : : : ; n 1; (j = 1; : : : ; n)
donde kj = 0 si k 6= j y kj = 1 si j = k. Luego se resuelve ada uno de los problemas planteando
y(x) =
1
X
n=0
n (x x0 )n ;
desarrollando en series de poten ias ada una de las aj (x) y utilizando las reglas de al ulo dadas en la
proposi ion (I.4.16).
Ejemplo
7.- Determinemos la solu ion general de
y00 + xy = 0:
Observamos que la e ua ion es (LH) y los oe ientes de las derivadas de y son analiti as, (a1 (x) = 0,
Elegimos por omodidad x0 = 0 ya que x ya esta
desarrollada en serie de poten ias. Planteando
a0 (x) = x), para ualquier punto de la re ta real.

y(x) =
1
X
n=0
n xn ;
se tiene Y n = y(n)(0)=n! y en parti ular y00(0) = 0y(0) da 2 = 20.

Por otro lado para n 3, se tiene
yn (x) = ( xy(x))(n
2)
nX2
k=0
x(k) y(n k
2)
(x) =
xy(n
2)
(x) (n 2)y(n 3)(x);
31
de donde
n =
(n 2)(n 3)! n
n!
(I.4.12)
Se observa inmediatamente que 3k+2 = 0.

Planteando 0 = y0 y 1 = y1, donde y0 y y1 son los valores ini iales y utilizando la formula (I.4.12),
se obtiene los valores de los restantes oe ientes.
Puntos Singulares Regulares

Deni ion I.4.18.- Se di e que
singular regular si las fun iones
el punto x0 singular de la e ua ion (LH) de orden n es un punto

an k (x x0 )k ; k = 1; : : : ; n
son fun iones anaaliti as.

Lo ideal es trabajar en puntos ordinariso, sin embargo la des rip ion de algunos fenomenos pasa por la
solu ion de problemas diferen iales a valor ini ial, en el que x0 es un punto singular regular.
Para resolver utilizando series enteras, se plantea
y(x) = (x x0 )
1
X
n=0
n (x x0 )
1
X
n=0
n (x x0 )n+
donde 2 R, (se debera determinar su valor) y se desarrolla en serie de poten ias las fun iones
an k (x)(x x0 )k para obtener una e ua ion para y rela iones re ursivas para los k .
Ejemplo
8.- Resolvamos la e ua ion diferen ial
x2 y00 + xy0 m2 y = 0:
Para que el problema tenga sentido, ualquier problema a valor ini ial en x = 0, debe satisfa er la
ondi ion que y0 = 0, dejando libre y1. Res ribiendo la e ua ion en la forma requerida, se tiene
1 m2 y = 0:
y00 + y0
x
x2
Se observa inmediatamente que 1 y m2 onstantes son analiti as.

Planteando
1
X
y(x) =
n xn+ ;
n=0
se obtiene
y0 (x) =
y00 (x) =
1
X
n=0
1
X
n=0
(n + ) nxn+
(n + )(n + 1) nxn+
Remplazando en la e ua ion diferen ial, se tiene

1
X
n=0
(n + )(n + 1) + (n + )
de donde obtenemos las e ua iones
1
X
n=0
(n + )2
(n + )2
= 0;
=0
m2 n xn+
m2 n xn+
m2 n = 0:
32
De estas e ua iones dedu imos que n = 0 para todo n ex epto para solamente uno de los n. Si
suponemos 0 6= 0, se tiene
2 m2 = 0;
por lo que = m o = m. De aqu obtenemos dos solu ioens diferentes:
y(x) = xm ; y(x) = x m :
5.- Realizando una substitu ion onveniente para x. Ver ejer i ios.
Solu ion de E ua iones Lineales (no Homogeneas)

Consideremos la e ua ion diferen ial lineal de orden n
(x)y(n 1) + + a1(x)y0 + a0(x)y = b(x):

(L)
donde las ai (x) y b(x) son fun iones ontinuas denidas sobre un intervalo. Suponemos que b no es identi amente nula, el aso uando b es la fun ion nula ya ha sido estudiado. A la e ua ion (L) le aso iamos la
e ua ion
y(n) + an 1 (x)y(n 1) + + a1 (x)y0 + a0 (x)y = 0:
(LH)
y(n) + an
Proposi ion I.4.19.- Sea
'+
tambien es solu ion de (L).
Demostra ion.- Similar al aso de primer orden.

Proposi ion I.4.20.- Sea una solu ion parti ular de (L), enton es ualquier solu ion de (L) es de la forma
'+ ;

Demostra ion.- Similar al aso de primer orden.
Por lo tanto, para ono er la solu ion general de (L), es su iente ono er una solu ion parti ular de (L)
y la solu ion general de (LH), algo que ya sabemos saber en un gran numero de situa iones. Este resultado
lo expresamos, mediante la regla memote ni a
y(n) + an 1 (x)y(n
y(x0 ) = y0 ;
y0 (x0 ) = y1 ;
..
.
y(n
1)
(x0 ) = yn
1)
+ + a1(x)y0 + a0(x)y = b(x)
Demostra ion.-Ejer i io.

Podemos determinar una solu ion parti ular mediante:
33
i) Al tanteo Es posible adivinar una solu ion en algunos asos, por ejemplo uando b(x) es un polinomio,
se puede intentar la solu ion parti ular on otro polinomio; si b(x) es sin o os se puede intentar on
una expresion que ontenga sin y os; si b(x) es una fun ion exponen ial, la solu ion parti ular debera
ser otra fun ion exponen ial. En todo aso el pro edimiento es el mismo que en el aso de primer orden.
ii) Varia ion de Constantes. Sea f'1; '2; : : : ; 'ng un sistema fundamental de solu iones de (LH). Al
igual que en el aso de primer orden se supone que la solu ion parti ular es de la forma
(x) =
Derivando , se obtiene
n
X
k=1
k (x)'k (x):
0 (x) = X 0 (x)'k (x) + X k (x)'0 (x);

k=1
suponiendo
n
X
k=1
se tiene
k=1
k (x)'0k (x) = 0;
n
0 (x) = X 0 (x)'k (x):

k=1
Derivemos 0 y obtenemos
00 (x) = X 0 (x)'0 (x) + X k (x)'00 (x);

k=1
suponiendo
n
X
k=1
se tiene
k=1
0k (x)'0k (x) = 0;
n
00 (x) = X k (x)'0 (x):

k
k=1
Repetimos el pro eso de deriva ion hasta determinar

n
X
y obteniendo por onsiguiente
k=1
(j )
k=1
0 (x)'(n
k
1)
(x) =
(x) +
n
X
k=1
n
X
k=1
k=1
0k (x)'k
(n 1)
(j )
(n)
n
X
k=1
1:
en la respe tiva e ua ion diferen ial, obteniendo
k (x)'k (x) +
(x) +
2:
k (x)'(kj) (x); j = 0; 1; : : : ; n
Reagrupando los terminos, se obtiene

n
X
, suponiendo por lo tanto
0k (x)'(kj) (x) = 0; j = 0; 1; : : : ; n
Finalmente introdu imos los valores de

n
X
(n 1)
nX1
j =0
aj (x)
k (x) 'k (x) +

(n)
n
X
nX1
j =0
k=1
k (x)'k (x)
(j )
= b(x):
aj (x)'k (x)A = b(x)
{z
(j )
=0 solu i
on de (LH)
34
Resumiendo hemos obtenido n e ua iones lineales para 01; 02; : : : ; 0n,

01 (x)'1 (x) + 02 (x)'2 (x) + + 0n (x)'n (x) = 0
01 (x)'01 (x) + 02 (x)'02 (x) + + 0n (x)'0n (x) = 0
...
...
(
n
1)
(
n
1)
(
n
1)
01 (x)'1 (x) + 02 (x)'2 (x) + + 0n (x)'n (x) = b(x)
que es rito de manera matri ial
0
'1 (x)
'2 (x) 'n (x) 1 0 01 (x) 1 0 0 1
0
'02 (x) '0n (x) C B 02 (x) C B 0 C
B '1 (x)
B
B
C B
C
.
... C

A .. A = ... A;
..
.
b(x)
'(n 1) (x) '2(n 1) (x) 'n(n 1) (x)
0 (x)
| 1
{z
} | n{z } | {z }
B (x)
C 0 (x)
W (x)
es de ir
W (x)C 0 (x) = B (x):
Ahora bien, omo f'1; '2 ; : : : ; 'n g es un sistema fundamental de (LH), W (x) es inversible, por lo que
C 0 (x) = (W (x)) 1 B (x);
obtenemos las fun iones k (x) integrando las respe tivas 0k (x).
Ejemplo
9.- Hallemos la solu ion general de
1
sin x :
La e ua ion (LH) aso iada, admite omo sistema fundamental f os x; sin xg. La apli a ion de
varia ion de onstantes ondu e a onsiderar el sistema de e ua iones

os x sin x 01(x) = 0 :
sin x os x 02(x)
1=over sin x
Resolviendo el sistema, por ejemplo por determinantes obtenemos:
os x ;
01 (x) = 1; 02 (x) =
sin x
de donde
1 (x) = x; 2 (x) = ln(sin x):
La solu ion general esta dada por onsiguiente
y(x) = 1 os x + 2 sin x x os x + ln(sin x) sin x:
y00 + y =
I.5 Ejer i ios
1.- En ontrar una solu ion parti ular, luego la solu ion general, de ada una de las e ua iones diferen iales
siguientes:
(a) y00 + y = e22
( ) y00 y = ex
( ) y + y = x
(d) y + 2y = sin x
35
I.5 Ejer i ios
2.- Resolver la e ua ion diferen ial de Bernouilli

(1 x2 )y0 xy = axy2
donde a es un parametro real.
3.- En ontrar la solu ion general de ada una de las e ua iones diferen iales siguientes:
a)
y0 = (x + y)2 ;
b)
(1 + x2)y0 + xy xy2 = 0;
)
y0 + y + (ex + sin x)y3 = 0:
4.- Dar la solu ion general de la e ua ion
y(4)
5.- En ontrar la solu ion general de
2y(3) + 2y(2) 2y0 + y = 0:

y000 y = 2x3 + 7:
6.- Resolver
(1 x2 )y00 xy0 + n2y = 0

sobre el intervalo [ 1; 1 ha iendo la substitu ion x = os t.
7.- Se onsidera la e ua ion diferen ial
x3 y000 + 2x2 y00 xy0 + y = 0:
a) Mostrar que la substitu ion
t = ln x
la transforma en una e ua ion a oe ientes onstantes.
b) Resolver de esta manera la e ua ion propuesta.
8.- a) Mostrar que si G(x) es una fun ion integrable sobre [a; b, enton es 8; 2 R existe una uni a solu ion
f de y00 = G(x) tal que f (a) = y f (b) = .
b) Supongamos que (b a) = 2! n, n un entero; enton es existe una innidad de solu iones f de la
e ua ion y00 = !2y tales que f (a) = f (b) = .
) 8; , existe una uni a solu ion f de y00 = 2 y tal que f (a) = y f (b) = .
9.- Resolver la e ua ion diferen ial
1 :
y00 + y =
sin x
10.- Resolver el problema a valor ini ial
y00 + y(y0 )3 = 0;
1
y( ) = 1;
6
1
y0 ( ) = 2:
6
g00 (x)
2xg0(x) 4g(x) = 0:
36
bus ando una solu ion de la forma
g(x) =
1
X
n=0
an xn :
>Cual es la solu ion mas general de esta forma? 0

>Cual es la solu ion de esta forma on g(0) = 0 y g (0) = 1?
12.- Comprobar que ada una de las siguientes e ua iones es homogenea y resolverla:
(a)
( )
y
x2 y0 = 3(x2 + y2 ) ar tan + xy;
x
y
y
x sin y0 = y sin + x;
x
x
13.- Resolver las siguientes e ua iones diferen iales

y0 = (x + y)2 ;
14.- Resolver las siguientes e ua iones
a)
y0 =
b)
y0 =
y0 =
(b)
(d)
x2 y0 3xy 2y2 = 0
p
xy0 = x2 + y2
y0 = sin2 (x + y
x+y+4
x y 6
x+y 1
;
x + 4y + 2
x+y+4
:
x+y 6
1):
;
15.- Ha iendo el ambio de variable z = y=xn y es ogiendo un valor ade uado de n, mostrar que las e ua iones
diferen iales siguientes pueden transformarse en e ua iones separadas y resolverlas:
1 xy2 ;
y0 =
a)
2x2y
2 + 3xy2 ;
y0 =
b)
4x2y
)
y0 =
y xy2
:
x + x2 y
y0 =
2y + x3 + x tan
x
y
:
x2
17.- Determinar la solu ion general de ada una de las siguientes e ua iones diferen iales:
a)
xy0 + y = x4 y3 ;
b)
xy2 y0 + y3 = x os x;
18.- Una solu ion de y0 sin 2x = 2y + 2 os x permane e a otada uando x ! =2. Hallarla.
19.- Resolver la e ua ion diferen ial xy0 = 2x2y + y ln y, utilizando la substitu ion z = ln y.
20.- Resolver las siguientes e ua iones:
(a) yy00 00 + (y20 )2 = 0; (b) xy2 0000= y0 + (0 y0)3 ; 0 2
( ) y ky = 0; (d) x y = 2xy + (y ) :
37
I.5 Ejer i ios
21.- Hallar la solu ion parti ular espe i ada para ada aso:
a)
(x2 + 2y0)y00 + 2xy0 = 0; y(0) = 1; y0(0) = 0;
1
yy00 = y2 y0 + (y0 )2 ; y(0) = ; y0 (0) = 1;
b)
2
)
y00 = y0 ey ; y(0) = 0; y0 (0) = 2:
22.- Una extension natural a las e ua iones de tipo lineal y tipo Bernouilli es la e ua ion de Ri ati
y0 = a(x) + b(x)y + r(x)y2 :
En general esta e ua ion no se puede resolver por metodos elementales. No obstante si se ono e una
solu ion parti ular y1(x), la solu ion general tiene la forma
y(x) = z (x) + y1 (x);
donde z(x) es la solu ion general de la e ua ion de Bernouilli
z 0 (q(x) + 2r(x)y1 (x))z = r(x)z 2 :
Demostrar esto y al ular la solu ion general de la e ua ion
y0 =
y
+ x3 y 2 x5 ;
x
que tiene omo solu ion parti ular evidente y1(x) = x.

b) Mostrar que la solu ion general tiene la forma de una familia uniparametri a de urvas
f (x) + g(x)
y=
:
F (x) + G(x)
23.- Aprove hando que y = x es solu ion de ada una de las e ua iones que se indi an a ontinua ion, hallar
su solu ion general.
1 y = 0;
x 0
a)
y +
y00
x 1
x 1
b)
x2 y00 + 2xy0 2y = 0
)
x2 y00 x(x + 2)y0 + (x + 2)y = 0:
24.- Comprobar que y(x) = ex es solu ion de
xy00
(2x + 1)y0 + (x + 1)y = 0
y hallar la solu ion general.

25.- Hallar la solu ion general de las siguientes e ua iones diferen iales lineales homogeneas a oe ientes
onstantes
a)
y000 3y00 + 2y0 = 0;
b)
y000 y = 0
)
y000 + y = 0
d)
y(4) + y000 3y00 5y0 2y = 0;
e)
y(5) 6y(4) 8y000 + 48y00 + 16y0 96y = 0:
38
26.- Determinar la solu ion general de ada una de las siguientes e ua iones diferen iales lineales:
(a) y002+ 4y =00 tan 2x;0
(b) y002 + y = se
x tan x;
2
2
( ) (x00 1)0y 2xy + 2y =x (x1 1) (d) x(6)y 2xy0 +1 2y = xe x;
(e) y 3y + 2y = (1 + e ) ;
(f) y y = x 0:
27.- La e ua ion de Chebi hef es
(1 x2 )y00 xy0 + n2y = 0;

on n entero. Veri ar que utilizando desarrollos en series de poten ia, se puede en ontrar una solu ion
polinomial.
28.- Resolver la e ua ion utilizando series alrededor de x = 0,
4x2 y00 8x2y0 + (4x2 + 1)y = 0:
29.- Determine aproximadamente las solu iones, utilizando polgonos de Euler e Iso linas de
a)
y0 = sin y + x
b)
y 0 = y 2 + x2
)
y 0 = y + x2
d)
y0 = tan y:
Cap
tulo II
Sistemas Diferen iales y Apli a iones
II.1 Con eptos B

asi os
Suponemos que el estudiante esta familiarizado on los espa ios ve toriales Rn y los fundamentos elementos
del al ulo de fun iones de varias variables.
Deni ion II.1.1.- Un sistema diferen ial de talla n y orden m es una expresion de la forma
F (t; x; x;
_ : : : ; x(m) ) = 0;
donde
n R n R n ! Rn :
F :R R
{z
}
|
m+1 ve es
ontinua, x(k) 2 Rn k = 0; 1; : : : ; m
Ejemplo
1.- El movimiento balsti o se rige por el sistema de e ua iones diferen iales de orden 2.
x = 0
y = g
donde x denota la omponente horizontal y y la omponente verti al. La variable independiente \t" no
apare e de manera expl ita.
Remar as
1.- Como la imagen de la fun ion F es Rn , el sistema se lo expresa omo

F1 (t; x; x;
_ : : : ; x(m) ) = 0;
F2 (t; x; x;
_ : : : ; x(m) ) = 0;
...
Fn (t; x; x;
_ : : : ; x(m)) = 0;
donde ada Fi : R R(m+1)n ! R es ontinua.
2.- Para evitarse analisis ompli ados en la resolu ion de di has e ua iones, se preere trabajar on sistemas
expli itados, es de ir, sistemas de la forma
x(m) = F (t; x; x;
_ : : : ; x(m 1) ):
Mas todava on sistemas de orden 1.
40
II Sistemas Diferen iales y Apli a iones
Proposi ion II.1.2.- Todo sistema de e ua iones diferen iales de orden m 1, es equivalente a un sistema
de primer orden.
Demostra ion.-Sea
x(m) == F (t; x; x;
_ : : : ; x(m 1));
on x 2 Rn un sistema de orden m y talla n. Introdu iendo las variables
y1 = x
y2 = x_
...
ym = x(m
se obtiene el sistema equivalente
1)
y_1 = y2 ;
y_2 = y3 ;
...
y_m 1 = ym ;
y_m = F (t; y1 ; y2; : : : ; ym );
y denotando y = (y1; : : : ; ym) 2 Rnm se tiene lo que se quiere.
Conven ion.-En lo que sigue el aptulo, solo se onsiderara sistemas de primer orden, a menos que se diga
lo ontrario.
Se dira que una fun ion ' : I R ! Rn ontinuamente derivable es una solu ion del
sistema de primer orden, de talla n,
x_ = f (t; x);
si '_ (t) = f (t; '(t)), 8t 2 I .
Deni ion II.1.3.-
Sistemas Autonomos
Se dira que un sistema diferen ial de talla n es autonomo si se puede es ribir de la

x_ = f (x);
(SDA)
es de ir, la variable independiente t no apare e expl itamente en el sistema diferen ial.
Deni ion II.1.4.-
forma
Proposi ion II.1.5.- Todo sistema diferen ial de primer orden, es equivalente a un sistema diferen ial
autonomo.
Demostra ion.- Sea
x_ = F (t; x);
un sistema diferen ial no autonomo, planteando

obtenemos esl sistema diferen ial autonomo
y deniendo el ve tor X =
x
xn+1
xn+1 = t;
x_ = f (xn+1 ; x)
x_ n+1 = 1;
2 Rn+1 , el sistema pre edente lo es ribimos omo

X_ = F (X ) =
f (xn+1 ; x) :
1
41
Deni ion II.1.6.- Cuando se trabaja on sistemas diferen iales autonomos (SDA), es ostumbre llamar a
Rn espa io fase y las oordnedas de x, xi fase.
Deni ion II.1.7.- Si ' : I ! Rn es una solu ion del (SDA), la imagen de ', '(I ), se llama traye toria
y en algunos asos linea de ujo.
Remar a.- Se designa t la variable independiente ha iendo en general referen ia al tiempo. Una solu ion
des ribira por onsiguiente la ley de movimiento de un objeto, mientras que la traye toria es la traza dejada
por el movimiento del objeto. La uni a informa ion que puede propor ionar una traye toria, es las posi iones por donde ha estado el objeto y no el momento exa to. Sin embargo, a pesar de esta limita ion se
puede obtener informa iones respe to a las solu iones ono iendo las traye torias. En la gra a ion de las
traye torias, es orriente indi ar el sentido del movimiento olo ando una e ha para indi ar. Ver gura II.1.
si os
II.1 Con eptos Ba
xn
x2
x1
Figura II.1.- Traye toria de una solu ion
Consideremos el sistema diferen ial autonomo

x_ = f (x);
y una solu ion ' : I ! Rn . El ve tor tangente, (en fsi a ve tor velo idad), en t = t0 es '0 (t0) 2 Rn . Ahora
bien
'0 (t0 ) = f ('(t0 ));
de donde a la traye toria en el punto x = '(t0 ) se le puede aso iar el ve tor tangente f (x ). Ver gura II.2
x*
Figura II.2.- Ve tor Tangente
Remar amos inmediatamente que el lado dere ho del (SDA) indu e un ampo de ve tores tangentes (a las
traye torias) del Sistema Diferen ial Autonomo, que lo denotamos
v : R n ! Rn
x 7! f (x):
42
Es ostumbre representar un ampo de ve tores en una gra a, aso iando a ada x 2 Rn del espa io de fases
una e ha que orresponde el ve tor de la imagen. Ver gura II.3.
Figura II.3.- Representa ion de un Campo de Ve tores

Tipos de Solu iones
Consideremos el sistema diferen ial autonomo de talla n

x_ = F (x);
y ' : I R ! Rn una solu ion del sistema.

Deni ion II.1.8.- Se dira que ' es una solu ion esta ionaria si '(t) = x 2 Rn para todo t 2 I .
Deni ion II.1.9.- Se dira que ' es una solu ion periodi a si existe T > 0 tal que '(t + T ) = '(t) para
todo T , si es el aso T se llama periodo.
Proposi ion II.1.10.- Sea x_ = F (x) un sistema autonomo, enton es:
a) ' es una solu ion esta ionaria, si y solamente si la traye toria de ' se redu e a un punto en el espa io
de fases, si y solamente si F ('(t)) = 0.
b) ' es una solu ion periodi a no esta ionaria si y solamente si la traye toria de ' es una urva errada y
F no se anula nun a sobre la traye toria.
Remar a.- La proposi ion pre edente para en ontrar solu iones esta ionarias nos da un medio algebrai o,
resolver F (x) = 0 y un medio gra o viendo las traye torias que se redu en a un punto. En las solu iones
esta ionarias, tenemos un medio gra o que el estudio de las traye torias erradas.
Ahora bien, poder de idir si una urva es errada o no, es en general bastante ompli ado si n > 2. En
el plano es posible en general.
Por otro lado, un gran numero de las apli a iones de los sistemas diferen iales autonomos es de talla 2.
Por lo que vale la pena estudiarlos.
Sistemas Autonomos de Talla 2
43
si os
II.1 Con eptos Ba
Aparte de los sistemas autonomos de talla 2

x_
y

=F

x
y

es fre uente estudiar las e ua iones diferen iales de segundo orden de la forma
x = f (x; x_ );
donde x depende de t omo un sistema autonomo de talla 2. Con tal motivo se plantea
x1 = x;
x2 = x_ ;
el sistema diferen ial autonomo equivalente esta dado por
_

x1 =
x2
x2
f (x1 ; x2 ) ;
y las omponentes del plano de fases son x1 = x (la posi ion) y x2 = x_ (la velo idad).
Volvamos al sistema
x_ = F (x);
donde x 2 R2 . F (x) indu e un ampo de ve tores tangentes a las traye torias. Aprove hando la rela ion
dx2 x_ 2
= ;
dx1 x_ 1
vista en Cal ulo I, y si F (x1 ; x2 ) = (F1 (x1 ; x2); F2 (x1 ; x2 )), las traye torias son los grafos de las solu iones
de la e ua ion diferen ial
dx2 F2 (x1 ; x2 )
=
;
dx1 F1 (x1 ; x2 )
o bien los grafos de las solu iones de la e ua ion diferen ial
dx1 F1 (x1 ; x2 )
=
;
dx2 F2 (x1 ; x2 )
A abamos de formular un metodo que nos determina las traye torias, sin ne esidad de ono er las solu iones
del sistema diferen ial autonomo.
Remar a.-Una familia de urvas en R2 pueden ser las traye torias de mu hos sistemas diferen iales, pero
el ampo de pendientes de las urvas es uni o.
Ejemplos
2.- (Predador-Presa).
El omportamiento de dos pobla iones de animales en un ambiente aislado en el

que una de las pobla iones son por ejemplo onejos y la otra pobla ion son lobos se puede modelar por
el sistema diferen ial
x_ = x( y);
y_ = y( x )
donde x(t) represanta la pobla ion de onejos en el instante t e y(t) representa la pobla ion de lobos en
el instante t. Las traye torias satisfa en la e ua ion diferen ial
y( x )
y0 =
x( y)
44
que es una e ua ion diferen ial separable. Resolviendo la e ua ion obtemos omo solu ion general (en
forma impl ita)
x ln x = ln y y + C:
(II.1.1)
Las traye torias del sistema diferen ial son las urvas de nivel de (II.1.1). En la gura II.4 estan
gra adas las urvas de nivel de (II.1.1) para = 1, = 1, = 2 y = 3. Observamos que las urvas
de nivel son urvas erradas, por lo que dedu imos que el omportamiento de estas pobla iones es li o.
Figura II.4.- Traye torias del Sistema Diferen ial
II.2 Apli a iones
Apli a iones Geometri as
Veremos, omo la utiliza ion de on eptos rela ionados a las e ua iones y sistemas diferen iales permite la
resolu ion de problemas rela ionados a familias de urvas, en parti ular urvas uniparametri as, del plano.
Deni ion II.2.1.- Una familia de urvas uniparametri as, es una familia de urvas si existe una fun ion
F : R2 R ! R, de manera que ada urva es el lugar geometri o de los puntos (x; y) que satisfa en la
e ua ion
F (x; y; ) = 0;
dejando 2 R jo. La e ua ion se llama e ua ion general de la familia de urvas y es el parametro.
Damos el siguiente resultado sin demostra ion.
Proposi ion II.2.2.- Una familia de urvas es uniparametri a, si y solamente si existe una fun ion ontinua
f : R R ! R2 , de manera que las fun iones ' : R ! R2 denidas por
' (t) = f (t; );
sean parametriza iones de ada una de las urvas de la familia.
45
II.2 Apli a iones
Ejemplo
1.- Las ir unferen ias entradas en el origen forman una familia uniparametri a. En efe to, la e ua ion
general es
x2 + y2 r2 = 0;
las parametriza iones de las ir unferen ias entradas en el origen estan dadas por
x(t) = r os t
y(t) = r sin t:
La solu ion general de una e ua ion diferen ial de primer orden que satisfa e una ondi ion
de Lips htiz es una familia uniparametri a de urvas. Asimismo, las traye torias pueden de un sistema
diferen ial autonomo forman una familia uniparametri a, si el sistema satisfa e una ondi ion de Lips htiz;
las solu iones del sistema vienen a ser las parametriza iones de las urvas.
Proposi ion II.2.3.- Para toda familia uniparametri a existe una e ua ion diferen ial y0 = f (x; y) uya
Remar a.-
solu i'on general son las urvas de la familia uniparametri a. Ademas tambien existe un sistema diferen ial
autonomo uyas traye torias son las urvas de la familia uniparametri a.
Demostra ion.- Sea F (x; y; )

respe to a x, se obtiene,
= la e ua ion general de la familia de urvas uniparametri a. Derivando

F
F
(
x; y; ) + (x; y; )y0 = 0:
x
y
Por onsiguiente, y0; dejando x; y jos, satisfa en el sistema de e ua iones (algebrai as)
F (x; y; ) = 0
F
(x; y; ) + F
(x; y; )y0 = 0;
x
y
de donde es posible obtener y0 en fun ion de x e y, es de ir
y0 = f (x; y):
Para el sistema, es su iente onsiderar g1(x; y) y g2(x; y) de manera que
g (x; y)
f (x; y) = 2
;
g1 (x; y)
por lo que las urvas son traye torias del sistema
x_ = g1 (x; y)
y_ = g2 (x; y)
Remar a.-Para la existen ia de f (x; y) se debera onsiderar hipotesis suplementarias de F (x; y; ).
Ejemplo
2.- La familia de parabolas de verti e el origen y eje de simetra el eje y, satisfa e la e ua ion general
y = x2 ;
derivando, se obtiene
y0 = 2 x;
46
despejando de la e ua ion general obtenemos la e ua ion diferen ial

2y
y0 = :
x
Condi iones para que una familia de urvas sea uniparametri a
Hemos visto que para que una familia de urvas sea uniparametri a es ne esario la existen ia de una e ua ion
general. Sin embargo en ontrar di ha e ua ion es en general una tarea ompli ada.
Ahora bien, utilizando resultados sobre las traye torias y solu iones de e ua iones y sistemas diferen iales, uno se puede dar uenta si una familia es uniparametri a o no. En efe to, supongamos que tenemos
una familia de urvas uniparametri a y sea
_

x = F (x; y) = F1 (x; y)
y
F2 (x; y)
un sistema autonomo uyas traye torias son las urvas de la familia en uestion. En general las fun iones
F (x; y) tienen a lo mas un numero nito de solu iones del sistema F (x; y) = (0; 0) en ada region a otada,
para efe tos pra ti os onsideramos regiones re tangulares de la forma [a; b [ ; d.
Por otro lado observamos que si en el punto (x; y), F (x; y) 6= 0, solo pasa por este punto una traye toria.
Si F (x; y) = 0, (x; y) es un punto esta ionario en que enventualmente pueden on uir y o emerger varias
traye torias. En onse uen ia, si en una region re tangular existen una innidad de puntos en los uales
pasan diferentes urvas, esta familia no es uniparametri a.
Ejemplos
3.- La familia de ir unferen ias del plano, no es una familia uniparmetri a porque por ada punto del plano
pasa una innidad de ir unferen ias.
4.- La familia de ir unferen ias de entro en el eje x y que pasan por el origen, ver gura II.5, puede ser una
familia uniparmetri a, solo existe un punto, el origen, por el ual pasa una innidad de ir unferen ias,
por los otros, pasa exa tamente una ir unferen ia.
Figura II.5.- Familia de Cir unferen ias
47
II.2 Apli a iones
Ahora bien, esta es una familia uniparametri a, uya e ua ion general es fa ilmente dedu ible y es
(x r)2 + y2 r2 = 0;
o bien
x2 + y2 2rx = 0:
Determinemos la e ua ion diferen ial y es sistema diferen ial. Se tiene
x2 + y2
r=
;
x
derivando obtenemos
2
0
0 = 2x + 2yyxx2
despejamos y0, se obtiene
y0 =
de donde el sistema estara dado por
x2
y2
y 2 x2
2xy ;
x_ = 2xy
y_ = y2 x2
Remar a.- Es importante en ontrar un sistema diferen ial autonomo uyas traye torias sean las urvas de
una familia de urvas uniparametri a, porque nos permite determinar un ampo de ve tores tangentes. Y
manipular ve tores (tangentes) es mu ho mas sen illo que manipular las urvas.
Familias de urvas que forman un angulo dado on una familia de urvas dada
Consideremos el siguiente problema: \Dada una familia de urvas uniparametri a, en ontrar una familia de
urvas que forme un angulo ". Interpretemos lo que signi a este problema. El on epto de angulo es un
on epto para re tas, que se generaliza a las urvas, deniendo el angulo entre dos urvas, omo el angulo
de los ve tores tangentes, ver gura II.6.
D
u
v
C
Figura II.6.- Angulo entre dos urvas
Por onsiguiente, determinamos primero el ampo de ve tores tangentes a la familia de urvas dada, ya
sabemos omo, sea

v
1 (x; y )
v(x; y) = v (x; y)
2
el ampo de ve tores dados. Obtenemos el ampo u(x; y) de la familia de urvas que forman el angulo dado,
rotando por un angulo , el ampo v(x; y), es de ir

u1 (x; y) = os
sin v1(x; y) :
u2 (x; y)
sin os v2(x; y)
El siguiente paso es onsiderar el sistema autonomo
_

x = u1 (x; y)
y
u2 (x; y)
48
y determinar las traye torias, que ya sabemos ha er.
Ejemplo
5.- Consideremos nuevamente la familia de ir unferen ias de entro en el eje x y que pasan por el origen y
en ontremos la familia de urvas ortogonales a esta familia de ir unferen ias. En el ejemplo 4, hemos
en ontrado un ampo de ve tores tangentes, dado por

2
xy
v(x; y) =
:
y 2 x2
Apli ando una rota ion de 90 grados obtenemos

2xy = x2 y2 ;
u(x; y) = 01 01
y 2 x2
2xy
de donde, las urvas ortogonales, son traye torias del sistema diferen ial
_ 2

x = x y2 :
y
2xy
Determinamos las traye torias, resolviendo la e ua ion diferen ial
2xy ;
y0 =
x2 y2
uya solu ion general es
x2 + y2 y = 0;
que es la e ua ion general de una familia de ir unferen ias de entro el eje y y que pasan por el origen.
Figura II.7.- Curvas ortogonales
49
II.2 Apli a iones
Problemas de Perse u ion
Los problemas de perse u ion, son interesantes porque onstituyen un buen ejer i io mental, desde la omprension del fenomeno, la formula ion matemati a del modelo y la misma solu ion.
Un ejemplo trabajado
En un momento dado, un perro se en uentra a una distan ia a a la dere ha de un gato, al que omienza
a perseguirlo. El gato se es apa a un arbol que se en uentra a una distan ia h en dire ion de su frente.
Suponiendo que el gato orre on una rapidez u inferior a la rapidez v del perro, se tiene las siguientes
preguntas.
a) >Cual debe ser la distan ia maxima para que el gato no sea atrapado por el perro?
En aso que el perro atrape al gato,
b) >Cual sera el tiempo que trans urre para que el gato sea atrapado?
) >Que longitud habra re orrido el perro para atrapar al gato?
Para resolver el problema, ha emos una representa ion de las posi iones del perro y el gato en el plano
artesiano. Por onsiguiente, la posi ion del gato y el perro resp tivamente en el momento ini ial, estan
dadas por (0; 0) y (a; 0) respe tivamente. Las oordenadas del arbol son (0; h).
El siguiente paso es modelar la perse u ion, para eso, denotamos por g : [0; +1) ! R2 y p : [0; +1) !
2
R las leyes de movimiento del gato y el perro respe tivamente. Puesto que el gato se dirige al
arbol, la ley
de moviento para el gato sera

g(t) = 0 :
ut
La parte mas interesante del problema, es la formula ion matemati a de la ley de movimiento del perro. Para
tal efe to, debemos determinar la manera omo persigue el perro al gato. Institivamente en ada instante el
perro bus ara el amino mas orto que le permita atrapar al gato. Intuitavemente sabemos que el amino mas
orto que une dos puntos es el segmento que une estos dos puntos. Por onsiguiente la velo idad del perro,
en el instante t tendra la misma dire ion del ve tor uyo origen es la posi ion del perro y su extremidad es
la posi ion del gato en el instante t. Ver gura II.8. Por lo tanto, la velo idad del perro sera
v
!
p_(t) =
p(t)g(t):
!

p(t)g (t)
g(t)
p(t)
Figura II.8.- Movimiento del perro
Expresando en oordenadas, obtenemos el sistema diferen ial de primer orden, que des ribe el movimiento
del perro,
_

v
x =p
x :
y
x2 + (ut y)2 ut y
Para responder las interrogantes planteadas, podemos resolver el sistema, lo que por el momento no sabemos
todava ha er, o bien determinar la traye toria del perro.
Observamos que el sistema obtenido no es autonomo, porque apare e expli itamente la variable t; sin
embargo podemos suponer que t depende de x, lo que da
y0 =
ut y
;
x
50
obteniendo as, la e ua ion de primer orden
xy0 = y ut:
Esta e ua ion no la podemos resolver tal omo esta planteada, por que tenemos 2 fun iones in ognitas.
Derivamos para obtener
xy00 = ut0;
pero t0 = 1=x_ , de donde
xy00 = ut0
= x_u
p
2
2
= u x + (ut y)
vp
x
u x2 + x2 y02
;
v
x
=
planteando r = u=v y viendo la gura II.8 es razonable suponer que x 0, por lo que la traye toria es la
solu ion de la e ua ion diferen ial de segundo orden
p
xy00 = r 1 + y02 ;
que satisfa e los valores ini iales
y(a) = 0; y0 (a) = 0:
Planteando z = y0, obtenemos la e ua ion diferen ial de tipo separable
p 02
1+z ;
z0 = r
x
que al resolverla, se obtiene
z + 1 + z 2 = Cxr :
Ahora bien, z(a) = y0(a) = 0, de donde C = a r . Despejando z de la expresion pre edente, se tiene
z=
Integrando, se obtiene
1
2
a r xr

ar x r :
1 a r xr+1 ar x1 r + D:
2 r+1
1 r
Utilizando el he ho que y(a) = 0, se tiene
ar
D=
1 r2 :
En la gura II.9, se tiene gra ada la traye toria para a = 100, u = 1 y v = 2.
y=
51
II.3 Sistemas Diferen iales Lineales
Figura II.9.- Traye toria del Perro
Obseervamos que para x = 0, se tiene
y=D=
auv
:
v u2
2
Por lo tanto, para que el gato no sea atrapado, la distan ia maxima a la que se debe en ontrar el arbol debe
ser
auv
:
v 2 u2
Si el gato es atrapado, lo sera en el punto (0; auv=(v2 u2)), de donde el tiempo de perse usion sera
av
T = 2 2;
v u
y la longitud que habra re orrido el perro sera
L=
av2
:
v u2
2
Un sistema diferen ial lineal de talla n, es un sistema diferen ial que puede expresarse omo
x_ 1 = a11 (t)x1 + a12 (t)x2 + + a1n (t)xn + b1(t)
x_ 2 = a21 (t)x1 + a22 (t)x2 + + a2n (t)xn + b2(t)
...
x_ n = an1 (t)x1 + an2 (t)x2 + + ann (t)xn + bn (t)
52
donde aij ; bi : I ! R ontinuas.
Por uestiones de omodidad en el manejo de smbolos, el mismo sistema se puede expresar de manera
matri ial, omo
x_ = A(t)x + b(t)
(L)
donde
0
1
0
a11 (t) a12 (t) a1n (t) 1
0
b
x1 1
1 (t)
B a21 (t) a22 (t) a2n (t) C
B . C
.
x = ... A ; A(t) = B
... C
..
A ; b(t) = . A :
.
xn
bn (t)
an1 (t) an2 (t) ann (t)
Se dira que un sistema diferen ial lineal es homogeneo si b(t) = 0; es de ir si el sistema es de la forma
x_ = A(t)x:
(LH)
Sistemas Diferen iales Homogeneos

Teorema II.3.1.- La solu ion general del sistema diferen ial linal homogeneo de talla n
x_ = A(t)x;
donde A(t) es una matriz uyos ooe ientes son fun iones ontinuas, es un subespa io ve torial de dimension
n.
Demostra ion.- Ejer i io para la parte de subespa io ve torial, a eptamos que la dimension es n.
Por onsiguiente, para ono er la solu ion general de un sistema diferen ial lineal homogeneo, es su iente ono er n solu iones linealmente independientes. A una tal familia f'1; : : : ; 'ng de solu iones linealmente independientes, se llama sistema fundamental de solu iones.
Conse uen ia del Algebra Lineal, se tiene la siguiente proposi ion que da riterios para de idir si un
sistema de solu iones es un sistema fundamental.
Proposi ion II.3.2.- Sea f'1 ; '2 ; : : : ; 'n g un sistema de solu iones del sistema diferen ial lineal
x_ = A(t)x;
enton es f'1; '2; : : : ; 'ng es un sistema fundamental, si y solamente si
R(t) = ( '1 (t) '2 (t) 'n (t) )
es una matriz inversible, si y solamente si
det R(t) 6= 0:
Deni ion II.3.3.- Cuando f'1 ; '2 ; : : : ; 'n g es un sistema fundamental del sistema
x_ = A(t)x;
la matriz R(t) se llama matriz resolvente del sistema diferen ial lineal homogeneo.
Proposi ion II.3.4.- Si R(t) es una matriz resolvente del sistema diferen ial lineal de talla n
x_ = A(t)x;
enton es R(t) veri a:
i) R(t) es inversible.
ii) Toda solu ion '(t) del sistema se es ribe de manera uni a omo
'(t) = R(t)C;
53
donde C 2 Rn .
iii) R(t) veri a
R_ (t) = A(t)R(t):
Demostra ion.- El punto i) es onse uen ia de su deni ion.

ii) Las olumnas de R(t) forman un sistema fundamental de solu iones, del sistema diferen ial lineal. Denotemos por '1; '2 ; : : : ; 'n las olumnas de R(t). Por lo tanto, toda solu ion ' se es ribe de manera uni a
omo
'(t) = 1 '1 (t) + 2 '2 (t) + + n 'n (t);
introdu iendo la multipli a ion matri ial, se obtiene
0
'(t) = ( '1 (t) '2 (t)

|
{z
R(t)
1 1
B 2 C
B C
'n }) .. A :
| {z }
iii) Se tiene
R_ (t) = ( '_ 1 (t) '_ 2 (t) 'n (x) )

= ( A(t)'1 (t) A(t)'2 (t) A(t)'n (t) )
= A(t) ( '1(t) '2(t) 'n(t) )
= A(t)R(t):
Corolario II.3.5.- El problema a valor ini ial

x_ = A(t)x
x(t0 ) = x
Demostra ion.- Sea R(t) la resolvente del sistema, se tiene

x = R(t0 )C;
y omo R(t0) es inversible, C existe y es uni o.
Sistemas Diferen iales no Homogeneos
Consideremos nuevamente el sistema lineal
x_ = A(t)x + b(t);
donde x 2 Rn , A y b ontinuas. A este sistema le aso iamos el sistema linal homogeneo

x_ = A(t)x
Al igual que las e ua iones diferen iales lineales del aptulo I, se tiene:
Proposi ion II.3.6.- Sea
+'
(L)
(LH ):
54
es una solu ion de (L).

Proposi ion II.3.7.- Sea una solu ion parti ular de (L), enton es ualquier solu ion de (L) es de la forma
+ ';
Por lo tanto, para ono er la solu ion general de (L), es su iente ono er una solu ion parti ular de (L) y
la solu ion general de (LH). Este resultado lo expresamos, mediante la siguiente regla memote ni a.
Proposi ion II.3.8.- El problema a valor ini ial
y0 = a(t)y + b(t);
y(t0 ) = y0

Al igual que en el aptulo pre edente, se puede determinar una solu ion parti ular al tanteo, o bien utilizando
el metodo de varia ion de onstantes.
Sea R(t) la resolvente del sistema (LH), y suponemos que la solu ion parti ular bus ada es de la forma
(t) = R(t)C (t);
derivando, se tiene
_ (t) = R_ (t)C (t) + R(t)C_ (t) = A(t)R(t)C (t) + R(t)C_ (t) = A(t) (t) + b(t);
de donde
R(t)C_ (t) = b(t);
por lo que
Z t
C (t) = R 1 (s)b(s) ds:
t0
Ejemplo
1.- Consideremos el sistema diferen ial de talla 2
x_ = y + sin t
y_ = x:
Expresando de manera matri ial, se tiene

_

x = 0 1
x + sin t :
y
1 0 y
0
Se puede veri ar que

t sin t
R(t) = os
sin t os t
55
es una matriz resolvente del sistema lineal homogeneo aso iado. Apli ando varia ion de onstantes, se
obtiene

os t sin t _1 = sin t :
0
sin t os t _2
Resolviendo el sistema lineal para _1 y _2 , se tiene
_1 = sin t os t; 2 = sin2 t;
lo que da
1
1 1
1 (t) = sin2 t; 2 (t) = t + sin 2t:
2
2 4
Por lo tanto, la solu ion general sera
1 1
1
x(t) = 1 os t + 2 sin t + sin2 t os t ( t
2
2 4 sin2t) sin t;
1
1 1 sin2t) os t:
y(t) = 1 sin t 2 sin t + sin3 t + ( t
2
2 4
Remar a.- Para poder realizar el metodo de varia ion de onstantes hemos supuesto la existen ia de la
matriz resolvente, por onsiguiente la existen ia de un sistema fundamental de solu iones. Sin embargo,
uando la matriz A(t) ontinua es arbitraria, no existe un metodo que permita determinar la resolvente.
Sistemas Lineales a Coe ientes Constantes
Consideremos el sistema lineal de talla n

donde
x_ = Ax;
a11 a12
B a21
B
A=
...
an1 an2
a1n 1
ann
C
C;
A
aij 2 R:
El aso no homogeneo ha sido ya visto mas arriba.

Cuando n = 1, se tiene la e ua ion lineal de primer orden x_ = ax, uya solu ion general es
x(t) = eta = eta ;
por lo que R(t) = eta. Por otro lado, re ordemos que
x
+
2 3! + :
La deni ion mediante una serie de la fun ion ex, permite generalizar la fun ion exponen ial para matri es
de la siguiente manera.
Deni ion II.3.9.- Sea A una matriz n n a oe ientes reales, se dene
ex = 1 + x +
x2
1
2
1
3!
eA = I + A + A 2 + A 3 + =
1
X
1 Ak :
k=0 k !
Aunque el al ulo de eA parez a muy ompli ado, existen tipos de matri es en las que el al ulo es
inmediato y sen illo.
Ejemplos
56
2.- Sea D una matriz diagonal
0
B
D=B

1
0
1
C
... C
A
0 0 n
la utiliza ion de la deni ion de la exponen ial de una matriz, permite llegar a
1
0
e1 0

C
B 0 e0
eD = B
... C
A

0 0 en
que es una matriz diagonal, uyos oe ientes de la diagonal son las exponen iales de los oe ientes de
la matriz diagonal D.
3.- Si A es una matriz diagonalizable; es de ir que existe una matriz T inversible tal que
A = T DT 1;
donde D es una matriz diagonal. Se tiene
1
eA = I + A + A2 +
2
= I + T DT 1 + 21 T D T| {z1T} DT 1 + dots
=I
1
1
1
= T IT + T DT + 2 T D2T 1 +

1
2
= T I + D + 2D + T 1
= T eD T 1 :
Por el ejemplo 3, observamos que el al ulo de la exponen ial de una matriz diagonalizable es sen illo a
ondi ion de ono er la matriz diagonal a la que es semejante.
Repaso de Algebra Lineal

Sea A una matriz de n n,
supongamos que sea diagonalizable, es de ir que existen D diagonal y T

A = T DT 1;
denotando T = ( T1 T2 Tn ) por sus olumnas, se tiene
0
1
1 0
B 0 2
C
A ( T1 T2 Tn ) = ( T1 T2 Tn ) B
... C

A;
0 0 n
( AT1 AT2 ATn ) = ( 1T1 2 T2 n ) :
Como las T es inversible, las olumnas Tk son ve tores no nulos, de donde los Tk son ve tores propios respe to
a k valor propio de A. Mas todavia los Tk forman una base de Rn . Hemos mostrado:
inversibles tales que
Proposi ion II.3.10.- Una matriz A es diagonalizable, si y solamente si de los ve tores propios de A se
puede formar una base de Rn .
Remar a.- La determina ion
de valores propios de una matriz A pasa por la determina ion de las rai es
del polinomio ara tersti o de A
PA (x) = det(xI A):
57
La Resolvente de un Sistema Lineal a Coe ientes Constantes

Teorema II.3.11.- Consideremos el sistema lineal homogeneo a oe ientes onstantes
x_ = Ax:
Enton es
R(t) = etA
es una resolvente del sistema; es de ir la solu ion general es
x(t) = etAC:
Demostra ion.- Es su iente ver que
y que R(t) = etA es inversible. En efe to
R_ (t) = A(t)R(t);
1 1
1 tk Ak = X
(tk )Ak
k
!
k
!
k=0
k=0
1
X
1
X
k k 1 k
t A
k
k=1 !
1
X
= A (k 1 1)! tk 1 Ak 1
k=1
1
X
= A k1! tk Ak = AetA:
k=0
Si AB = BA se tiene eA+B = eAeB , demostra ion que dejamos omo ejer i io. Por lo tanto
I = e0 = etA tA = etAe tA ;
de donde etA es inversible.
Ejemplo
4.- Consideremos el sistema lineal a oe ientes onstantes

x_ = x + 2y
y_ = x + 4y:
Es rito de manera matri ial, se tiene
_

x = 1 2
x
y
1 4 y :
Determinemos los valores propios de la matriz, el polinomio ara tersti o es
p() = 2 5 + 6 = ( 2)( 3);
Por lo tanto, 2 y 3 son valores propios de la matriz aso iada al sistema. Determinemos los ve tores
propios. Para = 2, se tiene la e ua ion
x + 2y = 0
58
planteando y = 1, se tiene x = 2. Para = 3, se tiene la e ua ion

2x + 2y = 0;
que da fa ilmente x = 1, y = 1. Por lo tanto

2
1
T= 1 1
y
T

= 11 21 :
La solu ion general del sistema esta dada por

x(t) = 2 1
e2t 0
1 1 1
y(t)
1 1 0 e3t
1 2 2
que desarrollando, se obtiene
x(t) = 2( 1 2 )e2t + (2 2 1 )e3t ;
y(t) = ( 1 2 )e2t + (2 2 1 )e3t :
Variante del Metodo de la Exponen ial
Hemos visto que la determina ion de la exponen ial de una matriz (diagonalizable) pasa por la determina ion
de los valores propios y la determina ion de la matriz inversible T . Ahora bien, al ular ve tores propios no
es una tarea sen illa y simple, lo mismo que invertir una matriz. Por otro lado, las matri es no siempre son
diagonalizables, sin que eso no signique no exista la exponen ial.
A ontinua ion formularemos un pro edimiento que nos permita al ular la solu ion general de un
sistema lineal homogeneo a oe ientes onstantes. Consideremos el sistema
x_ = Ax;
Por el momento supondremos que A es diagonalizable en los reales y que los valores propios son diferentes.
Por onsiguiente
0
1
1 0
B 0 0
C
1
A=TB
... C

AT :
0 0 n
La solu ion general del sistema estara dada por
0 1 t
1
e
0 t
B 0
C
e2
C T 1 C:
x(t) = T B
...

A

t
n
0 0 e
Ahora bien, T 1C sigue siendo un ve tor onstante, que lo denotamos una vez mas C ; por onsiguiente
0
1 e1 t 1
x(t) = T 2 e2 t A ;
... ent
n
59
por otro lado, los oe ientes de T pueden ser vistos omo onstantes, de donde la solu ion puede es ribirse
omo
0
1 t
2 t
n t 1
11 e
x(t) =
B 21 e1 t
B

...
12 e
22 e2 t
n1 e1 t n2 e2 t
1n e
2n en t C
C
nn en t
Remar amos que la solu ion esta expresada en fun ion de n2 onstantes, omo la solu ion general es de
dimension n, solo se requiere n onstantes, por lo que n2 n deben ser expresadas en fun ion de n onstantes.
Esto se obtiene remplazando en el sistema y obteniendo e ua iones lineales que deberan ser resueltas.
Como la determina ion de los i pasa por la solu ion del polinomio ara tersti o de la matriz A y las
rai es no siempre son diferentes y no siempre son reales, podemos anar nuestro resultado sobre la solu ion
de la siguiente manera.
Teorema II.3.12.- Consideremos el sistema diferen ial lineal a oe ientes onstantes
x_ = Ax:
Sean r1 ; : : : ; rm las rai es diferentes del polinomio ara tersti o de la matriz A y uya fa toriza ion esta
dada por
PA () = ( r1 )n1 ( r2 )n2 ( rm )nk ;

y onsideremos la familia de fun iones i : R ! R, i = 1; : : : ; n onformada por:
si rj 2 R, rj ontribuye a la familia on
erx; xerx; : : : ; xnj 1 erx;
si rj
= + i on 6= 0, rj y rj ontribuye on
ex os x; xex os x; : : : ; xnj
ex sin x; xex sin x; : : : ; xnj
ex os x
1 x
e sin x:
1
Enton es la solu ion general puede ser expresada omo
11 1 (t) 12 2 (t)
B 21 1 (t) 22 2 (t)
0
x(t) = B

donde n2
..
.
n1 1 (t) n2 2 (t)
1n n (t) 1
2n n (t) C
C
nn n (t)
n de las ij dependen de n ij .
Ejemplo
5.- Consideremos el sistema
x_ = x + 2y
y_ = 2x y:
La matriz aso iada al sistema es por onsiguiente

A=
uyo polinomio ara tersti o es
1 2 ;
2 1
PA () = 2 + 2 + 1 + 4 = ( + 1)2 + 22 ;
por lo tanto, los valores propios son:

1 =
1 + 2i;
2 =
1 2i;
60
por el teorema pre edente, la solu ion general tendra la forma

x(t) = 11 e t os(2t) + 12 e t sin(2t)
y(t) = 21 e t os(2t) + 22 e t sin(2t):
Planteando, por ejemplo, 11 = 1 y 12 = 2, remplazando en la primera e ua ion se obtiene
( 1 + 2 2)e t os(2t) + ( 2 1 2)e t sin(2t) = ( 1 + 2 21)e t os(2t) + ( 2 + 2 22)e t sin(2t):
De donde, 21 = 2 y 22 = 1 y la solu ion general sera
x(t) = 1 e t os(2t) + 2 e t sin(2t)
y(t) = 2 e t os(2t) 1 e t sin(2t):
Un Metodo Alternativo
Tambien podemos determinar la solu ion general de un sistema diferen ial lineal a oe ientes sin tener que
utilizar la exponen ial de una matriz, onvirtiendo el sistema en una e ua ion diferen ial ordinaria lineal
respe to a una de las variables del sistema. Re ordemos la nota ion opera ional del primer aptulo
Dxj = x;
_ Ixj = xj
por lo tanto un sistema lineal a oe ientes puede expresarse omo
(a11D b11I )x1 + (a12 D b12I )x2 + + (a1nD b1nI )xn = 0
(a21D b21I )x1 + (a22 D b22I )x2 + + (a2nD b2nI )xn = 0
...
(an1D bn1I )x1 + (an2 D bn2I )x2 + + (ann D bnnI )xn = 0:
Luego se apli a el pro edimiento de elimina ion de Gauss, para obtener una e ua ion diferen ial para ada
una de las variables.
Ejemplo
6.- Consideremos el sistema diferen ial de talla 3
x_ = y + z =
y_ = x + z
z_ = x + y;
que la es ribimos en la nota ion opera ional omo
=0
=0
+
= 0:
Hagamos desapare er la variable x de las dos ultimas e ua iones, para eso se debe realizar
Primera e ua ion + D(Segunda e ua ion)
Primera e ua ion + D(Ter era e ua ion)
obteniendo de esta manera
(D2 I )y
(D + I )z = 0
(D + I )y + (D2 I )z = 0
eliminamos ahora la variable y, multipli ando por (D I ) la ultima e ua ion. De esta manera
(D + I )(D 2I )Dz = 0;
Dx
Ix
Ix
Iy
Dy
Iy
Iz
Iz
Dz
II.4 Estabilidad y Estudio de Puntos Cr

ti os
por lo que
61
z (t) = 31 e t + 32 e2t + 33 1:
Por simetra del problema, se obtiene las mismas e ua iones para x e y, de donde
x(t) = 11 e t + 12 e2t + 13 1;
y(t) = 21 e t + 22 e2t + 23 1:
Remplazando en las e ua iones del sistema, se obtiene:
11 e t + 2 12 e2t = ( 21 + 31 )e t + ( 22 + 32 )e2t + ( 23 + 33 );
21 e t + 2 22 e2t = ( 11 + 31 )e t + ( 12 + 32 )e2t + ( 13 + 33 );
31 e t + 2 32 e2t = ( 11 + 21 )e t + ( 12 + 22 )e2t + ( 13 + 13 );
Resolviendo las e ua iones lineales que resultan, se tiene
13 = 23 = 33 = 0; 31 = ( 11 + 21 ); 32 = 22 = 12 :
Planteando 11 = 1, 21 = 2 y 12 = 3, la solu ion general del sistema esta dada por
x(t) = 1 e t + 3 e2t ;
y(t) = 2 e t + 3 e2t ;
z (t) = ( 1 + 2 )e t + 2 e2t :
ti os
En esta se ion estudiaremos los on eptos de estabilidad, para tal efe to, onsideremos el sistema diferen ial
de primer orden y talla n
x_ = f (t; x):
Introdu imos la nota ion
x(t; t0 ; x0 );
que indi a que es la solu ion al problema a valor ini ial
x_ = f (t; x);
x(t0 ) = x0 :
Con la nota ion introdu ida podemos plantear las situa iones on las que se puede onfrontar en la solu ion
de e ua iones y o sistemas diferen iales, ademas de los problemas diferen iales a valor ini ial.
Usualmente los sistemas diferen iales, al igual que las e ua iones diferen iales son el lenguaje matemati o
para des ribir fenomenos, pero tambien en la mayora de los asos son simpli a iones de la realidad. Por
onsiguiente, la formula ion diferen ial (simplista), la podemos es ribir
x_ = f (t; x);
y la situa ion real
x_ = f (t; x) + g(t; x):
En general no se puede expresar g(t; x).
62
La otra situa ion, es que se debe resolver a problemas a valor ini ial
x_ = f (t; x)
x(t0 ) = x0 :
Sin embargo, por la pre ision de las medidas y otros fa tores, uno resuelve
x_ = f (t; x)
x(t0 ) = x0 + x0 :
Proposi ion II.4.1.- La solu ion del problema
x_ = f (t; x) + g(t; x); x(t0 ) = x0 ;
esta, para jj peque~no, dada por
donde
Ejemplo
x (ta) = x0 (t) + x1 (t) + ;
x_ 0 (t) = f (t; x0 (t)); x0 (t) = x0

f
x_ 1 (t) = (t; x0 (t)) x1 (t) + g(t; x0 (t)); x1 (t0 ) = 0:
x
1.- Consideremos el problema a valor ini ial para la e ua ion de Bernouilli

x_ = x + x2 ; x(0) = 1:
El teorema pre edente nos permite des omponer la solu ion de este problema en dos problemas a valor
ini ial:
x_ 0 = x0 ; x0 (0) = 1;
uya solu ion es x0(t) = et y el problema
x_ 1 = 1x1 + (e2 t); x1 (0) = 0;
la solu ion de este problema, sera
x1 (t) = et + e2t ;
por lo que
x(t) = et + ( et + e2t) + :
Resolviendo la e ua ion de Bernouilli, se tiene
1
et
x(t) =
=
t
(1 + )e (1 + ) et ;
se dedu e que la primera solu ion es mas simple de al ularla.
Deni ion II.4.2.- Una solu ion x(t; t0 ; x0 ) del sistema
x_ = f (t; x);
se dira que es estable en el sentido de Liapunov, si 8 > 0, 9 > 0, tal que
kx0 k < y t t0 ) ky(t; t0 ; x0 + x0 ) y(t; t0 ; x0 )k < :
63

ti os
Se dira que es asintoti amente, si es estable y ademas existe 0 tal que

kx0 k < 0 ) t!lim
ky(t; t0 ; x0 + x0 ) y(t; t0 ; x0 )k = 0:
+1
La solu ion es inestable si no es estable.
Ejemplos
2.- Consideremos el problema a valor ini ial
y0 = y
y(0) = y0 :
La solu ion dex esta e ua ion es y(x) = y0ex. Ahora bien, si < 0, se tiene ex ! 0 uando x ! +1
y si > 0, e ! +1 uando x ! +1. Por lo tanto, si < 0 toda solu ion y(x; 0; y0) es estable y
asintoti amente estable. Si > 0, toda solu ion sera inestable. Ver gura II.10.
Figura II.10.- Estabilidad de Solu iones
y0 = y(y
1):
Esta e ua ion tiene dos solu iones esta ionarias:

y(x) = 0; y(x) = 1:
y(x) = 0 es una solu ion asintoti amente estable, mientras que y(x) = 1 es inestable; ver gura II.11.
Figura II.11.- Estabilidad de Solu iones
Estabilidad de las Solu iones Esta ionarias
La resolu ion de sistemas diferen iales por medio de metodos analiti os es posible para un numero redu ido
de tipos de sistemas diferen iales. Para ontrarestar esta limita ion, hemos introdu ido metodos que nos
permiten estudiar las traye torias de di hos sistemas y en onse uen ia obtener informa ion respe to a las
solu iones. Un tipo importante de solu ion son las solu iones esta ionarias, porque en mu has apli a iones
se bus an solu iones esta ionarias.
Nuestro objetivo sera estudiar las solu iones de sistemas autonomos de talla 2 en las proximidades de
las solu iones esta ionarias y lasi ar los tipos de solu iones.
Consideremos el sistema autonomo de talla 2
_

x = f (x; y)
y
g(x; y)
64
Sea (x ; y un punto en que la solu ion es esta ionaria y suponiendo que f y g sean lo su ientemente
derivables, se tiene, ver urso Cal ulo II,
f !
f

f (x; y) f (x ; y) + x (x ; y ) y (x ; y )
x x ;
g
g
g(x; y)
g(x ; y )
y y
x (x ; y ) y (x ; y )
|
{z
}
=0
uando (x; y) es lo bastante proximo a (x ; y).

Por onsiguiente, uando se este er a del punto rti o (x ; y), planteando u = x
obtenemos el sistema que aproxima en las proximidades de (x ; y).
f !
f
_
u ;
u = x (x ; y ) y (x ; y )
g (x ; y ) g (x ; y )
v
v
x
y
x , v
=y
y ,
que no es nada mas que un sistema lineal a oefe ientes onstantes.

En onse uen ia, estudiar las solu iones alrededor de las solu iones esta ionarias, es estudiar los sistemas
de la forma
_

u :
u = a b
v
v
d
|
{z
La solu ion del sistema esta rela ionada a los valores propios de la matriz A. Sean 1 y 2 los valores propios
de A. Se tiene los siguientes asos:

ti os
Figura II.12.- Clasi a ion de puntos rti os
1.- 1 < 2 < 0. En un sistema ade uado de oordenadas , , se tiene

(t) = 1 e1 t ;
(t) = 2 e2 t :
65
66
Se observa que uando t ! +1, tanto (t) y (t) tienden a (0; 0). En este aso, se dira que (x ; y) es
un punto de equilibrio estable. Ver gura II.12.
2.- 0 < 1 < 2. En un sistema ade uado de oordenadas , , se tiene
(t) = 1 e1 t ;
(t) = 2 e2 t :
Se observa que uando t ! +1, tanto e1t , omo e2t divergen. En este aso, se dira que (x ; y) es un
punto de equilibrio inestable. Ver gura II.12.
3.- 01 < 0 < 2 . En un sistema ade uado de oordenadas , , se tiene
(t) = 1 e1 t ;
(t) = 2 e2 t :
Se observa que uando t ! +1, e1t tiende a 0, mientras que e2 t diverge. Se dira que (x ; y) es un
punto silla o umbre. Ver gura II.12.
4.- Supongamos que 1 = + i on 6= 0, si < 0, la solu ion en un sistema ade uado de oordenadas
y , sera
= 1 et os(t)
= 2 et sin(t)
uando t ! +1, (t) y (t) tienden a 0. La traye toria de este tipo de solu iones son espirales que
onvergen al entro de la espiral. (x ; y ) se llama fo o absorvente. Ver gura II.12.
5.- Supongamos que 1 = + i on 6= 0, si > 0, la solu ion en un sistema ade uado de oordenadas
y , sera
= 1 et os(t)
= 2 et sin(t)
uando t ! +1, (t) y (t) divergen. La traye toria de este tipo de solu iones son espirales que divergen
del entro de la espiral. (x ; y) se llama fo o repelente. Ver gura II.12
6.- Si 1 = i on 6= 0, las solu iones en un sistema ade uado de oordenadas seran
= 1 os(t)
= 2 sin(t)
En este aso, todas las traye torias son erradas. Por lo que (x ; y), se llama fo o li o. Ver gura
II.12.
Remar a.- Si 1 = 2 = , tenemos los siguientes es enario:
i) < 0 y la solu ion es un aso lmite de 1) y 4), en ambos asos la solu ion tiende al punto de equilibrio,
por lo que (x ; y) es un punto de equilibrio degenerado Ver gura II.12.
ii) > 0, la solu ion es un aso lmite de 2) y 5), se dira que (x ; y) es un punto de equilibrio inestable.
Ver gura II.12.
Ejemplos
4.- Consideremos el sistema diferen ial
x_ =x(y 1)
y_ =(x y)(x + y
3)
67
Son puntos rti os de este sistema (0; 0), (1; 1), (0; 3) (2; 1). En la gura II.13 se observan las traye torias
en el re tangulo [ 1; 4 [ 1; 4
II.5 Ejer i ios
Figura II.13.-Puntos Crti os
Se observa que los puntos (0; 0), (0; 3) y (2; 1) son puntos sillas y que el punto (1; 1) es un fo o repelente.
II.5 Ejer i ios
1.- Dados los sistemas diferen iales, determinar que solu iones son esta ionarias, uales son periodi as y
uales no son periodi as
(
x_ = x2 y2
a)
y_ = x2 + y2 ;
(
x_ = x + y 1
b)
y_ = x y + 1
2.- Gra ar el ampo de ve tores dado por los sistemas de la pregunta 1.
3.- Consideremos la e ua ion del pendulo simple
+ !2 sin = 0:
a) Convertir en un sistema diferen ial de orden 1. Cuales son las omponentes del espa io de fa es.
b) Determinar las traye torias de este sistema diferen ial.
68
) Suponiendo que el pendulo en el instante 0 satisfa e = 0, lasique las solu iones (esta ionarias,
periodi as y no periodi as), de a uerdo al valor que toma la velo idad angular ini ial.
4.- Determine un ampo de ve tores tangentes, siempre que sea posible y en uentre las familias de urvas
que forman un angulo re to on la familia de:
a) Re tas que pasan por el origen.
b) Cir unferen ias del plano.
) Cir unferen ias de entro la re ta y = x y que pasan por el punto (1; 1).
d) Las hiperbolas de entro el origen y asintotas las re tas y = x y y = x.
e) Las parabolas de eje de simetria la re ta y = 2x + 1 y verti e (3; 1).
f) Las parabolas de dire triz la re ta y + x = 0 y fo o sobre la re ta x y = 0.
g) Las ir unfereren ias de entro el origen.
5.- Determine la familia de urvas que forman un angulo de 45, siempre y uando sea posible, de las
familias del ejer i io 4.
6.- La vida media del Carbono 14 es de 5568 a~nos. En un gramo de arbon de madera re ientemente
produ ido hay en promedio 6.08 desintegra iones debidas al Carbono 14 por minuto. En un trozo
en ontrado en unas grutas en Fran ia, se midio 0.97 desintegra iones por mnuto y por gramo de arbon.
Estimar la fe ha probable de esta muestra.
7.- >Para que valores de , existe una solu ion (no nula) al problema
x + 2 x = 0;
x(0) = 0; x_ (1) = 0?
8.- En ontrar las solu iones generales de las e ua iones siguientes:

x_ + x = 0; x_ + x = 0; 2 x + x = 0; x + 2 x = 0:
Dis utir el aso en que tiende ha ia ero.
9.- Resolver las e ua iones:
x_ + x tx2 = 0; y x_ = ax bx2 donde a; b > 0:
10.- a) Mostrar que toda solu ion de x + x + x3 = 0 es periodi a.
b) Lo mismo para x + x3 = 0 y x + 1 +x x2 = 0:
11.- Estudiar la perio idad de las solu iones de x_ + x 2x3 = 0, suponiendo que
1
x_ (0)2 + x(0)2 x(0)4 6= :
4
12.- Dar la solu ion general de la e ua ion x_ = Ax, donde
0
2 1 11
A = 1 2 1A:
1 1 2
13.- Dar la solu ion general de la e ua ion
x(4) + 2x(2) + x = 0:
69
14.- Un punto P es arrastrado por el plano xy mediante una uerda P T de longitud a. Si T arran a del
origne y se mueve a lo largo del eje y positivo, y si P arran a del punto (a; 0). >Cual es la traye toria
de P? Esta urva se llama tra triz
15.- Un rio bastante an ho, sigue su ause siguiendo la dire ion norte sur on una velo idad a. Un bote en
el rio se hunde, el o upante de la bar a divisa un islote siguiendo la dire ion este a una distan ia b;
en la dire ion oeste a una distan ia > b se en uentra una de las orillas. Si la persona nada on una
velo idad d, indique a donde debe dirigirse para estar a salvo. Dis uta las diferentes posibilidades
16.- >Que urva situada por en ima del eje x tiene la propiedad de que la longitud del ar o que une ualesquiera dos puntos sobre ella es propor ional al area bajo di ho ar o?
17.- Determinar las solu iones generales de los sistemas:
II.5 Ejer i ios
x_ = x + 3y
;
y_ = 3x + y
x_ = x + 2y
y_ = 3x + 2y
18.- Determinar la solu ion general del sistema

(
19.- Hallar los puntos rti os de:

a)
b)
x_ = x + 2y + t 1
:
y_ = 3x + 2y 5t 2:
x + x_ (x3 + x2 2x) = 0
(
x_ = y2 5x + 6;
:
y_ = x y
20.- Para ada uno de los siguientes sistemas no lineales:

i) Hallar sus puntos rti os;
ii) Gra ar algunas traye torias;
(
x_ = y(x2 + 1)
a)
y_ = 2xy2 ;
(
x_ = y(x2 + 1)
b)
y_ = x(x2 + 1):
21.- Dibujar el diagrama de fases de la e ua ion x = 2x3 y veri ar que tiene en el origen un punto rti o
en el origen.
22.- Mostrar que (0; 0) es un punto rti o asintoti amente estable para ada uno de estos sistemas:
(
23.- La e ua ion de Van der Pol es
x_ = 3x3 y
;
y_ = x5 2y3
x + (x2
x_ 2x + xy2
:
y_ = x2 y2 y3
1)x_ + x = 0:
Investigar las propiedades de estabilidad del punto rti o x = 0 y x_ = 0 en el espa io de fases para los
asos > 0 y < 0.
70
24.- La mayora de los muelles en la realidad son no lineales. Un muelle no lineal se di e que es duro o blando
segun que la magnitud de la fuerza de re upera ion rez a mas o menos rapidamente que una fun ion
lineal del desplazamiento. La e ua ion
x_ + kx + ax3 = 0; k > 0;
des ribe el movimiento de un muelle duro si a > 0 y de uno blando si a < 0. Dibujar las traye torias en
el plano de fases y estudiar el punto rti o.
25.- Sea f : R ! R ontinuamente derivable y sea la e ua ion diferen ial
y00 = f (y; y0):
Si f (0; 0) = 0, mostrar que toda solu ion no nula de esta e ua ion tiene eros simples, si es que los tiene.
Ejemplos y00 + y = 0, y00 + sin y = 0.
26.- Mostrar que si A(x) es una matriz antisymetri a, enton es, se puede en ontrar un sistema fundamental,
respe to a la ual la resolvente es una matriz ortogonal.
27.- Sea
y0 = A(x)y
donde A(x) es una matriz de talla n. Mostrar que si se ono e una solu ion no trivial '(x) del sistema,
se puede redu ir el sistema a un problema analogo pero donde la matriz es de talla n 1.
Si se supone que ' 6= 0, bus ar una solu i
on y (x) = u'(x) + z donde u : R ! R y donde
0 1
0
B z2 C
C
z=B
.. A :
zn
28.- Utilizando el ejer i io 27, determinar un sistema fundamental de solu iones del sistema
0

y1 = 1=x
1 y1
y2
1=x2 2=x y2
sabiendo que
2
y(x) = x x
es una solu ion.
29.- En base al ejer i io pre edente, resolver el problema a valor ini ial

1
=x
1
0
2
y = 1=x2 2=x y + ( 1= rx ) ; y(1) = 00
30.- Consideremos la e ua ion del pendulo
y00 + sin y = 0; y(0) = ; y0 (0) = 0;
donde la amplitud se supone peque~na. Mostrar que la solu ion puede ser es rita bajo la forma
y(x) = y1(x) + 2 y2 (x) + :
Cal ular y1(x) y y2(x).
31.- Cal ular las solu iones esta ionarias de
y1 = y1 y22 2y2
y2 = y1 y12 y2
y mostrar que la solu ion nula es estable.
32.- Estudiar la estabilidad de los 4 puntos rti os del ejemplo 4 de la se ion II.4.
Cap
tulo III
Complementos
III.1 E ua iones Diferen iales Exa tas
A menudo se en uentran expresiones de la forma

a(x; y) dx + b(x; y) dy = 0;
(DE)
o
a(x; y; z ) dx + b(x; y; z ) dy + (x; y; z ) dz = 0;
(DE)
o mas generalmente
a1 (x) dx1 + a2 (x) dx2 + + an (x) dxn = 0;
(DE)
on x 2 Rn .
Limitaremos el estudio de este tipo de e ua iones al aso en que n 3, por omodidad supondremos
n = 3. El primer paso para determinar el signi ado de este tipo de e ua iones es saber lo que es el smbolo
dxi . Por lo he ho hasta ahora, las solu iones de los problemas diferen iales provenientes de e ua iones y o
sistemas pueden representarse mediante urvas. Por onsiguiente, si : R ! R, on
1
0
x1 (t)
C
(t) = B
x2 (t) A
...; x (t)
n
es una parametriza ion de una urva, se tiene por deni ion

dxi (t) = x_ i (t) dt;
donde dt : R ! R es la apli a ion lineal dada por
dt : R ! R
t 7! t:
Geometri amente dxi , llamada diferen ial de dxi mide la varia ion de xi (t) uando se modi a t.
Regresando a la e ua ion (DE), si x(t) es una parametriza ion que satisfa e la e ua ion (DE), se tiene
(a1(x(t))x_ 1 (t) + + an(x(t))x_ n (t)) dt = 0:
De donde, ualquier solu ion x(t) de (DE) satisfa e
a1 (x(t))x_ 1 (t) + + an (x(t))x_ n (t) = 0:
(III.1.1)
72
III Complementos
Remar a.- En el aso n = 2, utilizando el he ho que

dx2
dx1
= xx__ 2 ;
1
la e ua ion (DE), es equivalente a la e ua ion diferen ial ordinaria de primer orden

a1 (x1 ; x2 )
dx2
=
:
dx1
a2 (x1 ; x2 )
Por otro lado, las solu iones de (DE) forman sub onjuntos de Rn , para n = 2, se tiene urvas y para
n = 3 seran super ies. Estos onjuntos se representan en forma de e ua iones algebrai as omo
g(x; y) = ; o bien g(x; y; z ) = ;
donde es una onstante.
Ahora bien, si (x(t); y(t); z(t)) es la parametriza ion de una urva que esta dentro una super ie de
e ua ion g(x; y; z) = , se tiene
f (t) = g(x(t); y(t); z (t)) =
lo que signi a que f_(t) = 0, para todo t, apli ando la regla de la adena se tiene
f_(t) =
g
g
g
x_ (t) + y_ (t) + z_ (t) = 0:
x
y
z
(III.1.2)
Dedu imos que g(x; y; z) = C es solu ion de la e ua ion

g
g
g
dx + dy + dz = 0:
x
y
z
Por onsiguiente, si la e ua ion diferen ial (DE)

a(x; y; z ) dx + b(x; y; z ) dy + (x; y; z ) dz = 0;
es tal que existe una fun ion g : R3 ! R diferen iable es tal que
0
1
a(x; y; z )
grad g(x; y; ) = b(x; y; ) A = a(x; y; z)i + b(x; y; z)j + (x; y; z)k
(x; y; z )
se dira que es exa ta, g se llamara primitiva y la solu ion general del problema sera g(x; y; z) = .
(DE)
Ejemplos
y dx + x dy = 0;
Observamos que g(x; y) = xy es una primitiva de la e ua ion y por onsiguiente la solu ion general es
xy = C .
x dx + y dy + z dz = 0:
La fun ion
1 1 1
g(x; y; z ) = x2 + y2 + z 2
2 2 2
73
es una primitiva de la e ua ion y por onsiguiente la solu ion general es

1 x2 + 1 y2 + 1 z2 = C:
2 2 2
3.- Consideremos
y dx x dy = 0;
y queremos en ontrar una fun ion g(x; y) tal que
g
x
= y;
g
y
x;
dedu imos que tal g debera ser de la forma g(x; y) = yx + (y) si tomamos en uenta la derivada
respe to a x. Derivando g respe to a y obtenemos
g
y
= x + 0 (y) =
x;
lo que es imposible, porque se tendra 0 (y) = 2x y depende solamente de y no de x. Por onsiguiente

esta e ua ion no admite primitiva, lo que no signi a que no tenga solu ion.
Condi iones ne esarias y su ientes para la existen ia de primitivas
Si g : R3 ! R es dos ve es ontinuamente diferen iable, se tiene

rot(grad g(x; y; z)) = 0;
la veri a ion es sen illa y dejamos omo ejer i io. A ontinua ion enun iamos el siguiente teorema sin
demostra ion.
Teorema III.1.1.- La e ua ion diferen ial
a(x; y; z ) dx + b(x; y; z ) dy + (x; y; z ) dz
admite primitiva, si y solamente si
rot(a(x; y; z); b(x; y; z); (x; y; z)) = 0:

Corolario III.1.2.- La e ua ion diferen ial
a(x; y) dx + (x; y) dy = 0

admite primitiva, si y solamente si
a(x; y)
y
x; y)
= b(x
:
Utilizando el riterio del teorema y o orolario se puede determinar si una e ua ion (DE) admite o no
primitiva.
Ejemplo
4.- La e ua ion
y dx x dy = 0;
74
III Complementos
no admite primitiva por que
y
y
= 1 6= 1 = xx :
Determina ion de primitiva
Una vez que se ha determinado que la e ua ion diferen ial

admite primitiva, el problema onsiste en determinar la primitiva, en algunos asos la primitiva salta a la
vista, en otros no es tan evidente. Por onsiguiente es importante ono er un metodo que nos permita
determinarla.
Repasando Cal ulo II y el urso de Analisis Ve torial, la integral de lnea
Z
a(x; y; z ) dx + b(x; y; z ) dy + (x; y; z ) dz;
donde C es un ar o de urva orientado, en el aso en que

rot(a(x; y; z)i + b(x; y; )j + (x; y; z)k) = 0;
solo depende de las extremidades de C y no de la misma urva C . Ver gura III.1
B
A
Figura III.1.- Integrales de Lnea
Por onsiguiente es onsistente introdu ir la nota ion

Z B
a(x; y; z ) dx + b(x; y; z ) dy + (x; y; z ) dz =
en el aso en que rot(a(x; y; z)i + b(x; y; )j + (x; y; z)k) = 0.

Eligiendo un punto A 2 R3 y dejandolo jo, planteamos
g(x; y; z ) =
Z (x;y;z)
a(u; v; w) du + b(u; v; w) dv + (u; v; w) dw:
Se muestra que grad g(x; y; z) = a(x; y; z)i + b(x; y; z)j + (x; y; z)k.
Una vez que sabemos, omo determinar una primitiva para una e ua ion diferen ial exa ta, por medio
de una integral de lnea, el problema se resume a en ontrar un buen punto origen de las urvas y sobre todo
las buenas urvas.
Ejemplo.-
75
(x2 + y2) dx + 2xy dy = 0

Veriquemos que la e ua ion admite primitiva:
x2 +2
= 2y; 2xy = 2y:
y
x
Ahora elegimos omo origen de integra ion el punto (0; 0), luego integramos siguiendo la traye toria que
une (0; 0) on (0; y) y (0; y) on (x; y), lo que da
Z x
1
g(x; y) = (s2 + y2 ) ds = x3 + xy2 :
3
0
Por lo tanto, la solu ion general sera
1 x3 + xy2 = C:
3
Cuando la e ua ion diferen ia (DE) admite primitiva, la resolu ion se onvierte en un problema de resolver
una integral de lnea. Cuando la e ua ion no admite primitiva, se tiene varias op iones para en ontrar la
solu ion general de la e ua ion diferen ial (DE):
1.- Convertir en una e ua ion diferen ial ordinaria, esto es posible uando la e ua ion depende solamente
de 2 variables:
a(x; y) dx + b(x; y) dy = 0;
se onvierte en
a(x; y)
y0 =
:
b(x; y)
Ejemplo
6.- Resolvamos
Es equivalente a resolver
y dx x dy = 0:
uya solu ion general es y(x) = Celn x = Cx.
y0 = y
x
2.- Fa tores Integrantes
Si la e ua ion
a(x; y; z ) dx + b(x; y; z ) dy + (x; y; z ) dz = 0
no admite primitiva, se puede multipli ar di ha e ua ion por una fun ion (x; y; z) no identi amente
nula de manera que
(x; y; z )a(x; y; z ) dx + (x; y; z )b(x; y; z ) dy + (x; y; z ) (x; y; z ) dz = 0
admita una primitiva. La fun ion se llama fa tor integrante. En general determinar un ade uado
puede resultar una tarea mu ho mas di il que resolver la misma e ua ion, pues el fa tor integrante es
solu ion de una e ua ion a derivadas par iales, que en prin ipio es mas di il resolverla que una e ua ion
diferen ial ordinaria.
Ejer i io.- Determinar la e ua ion a derivadas par iales en que el fa tor integrante es solu ion.
Sin embargo, existen algunos asos en que determinar el fa tor integrante salta a la vista, por ejemplo
si se puede separar la e ua ion, es de ir onvirtiendola en una e ua ion de la forma
p(x) dx + q(y) dy + r(z ) dz = 0:
76
III Complementos
Por otro lado, en el aso bidimensional
a(x; y) dx + b(x; y) dy = 0;
3.-
el fa tor intgrante satisfa e la e ua ion a derivadas par iales

a b
= x ;
y
que desarrollandola se tiene

a b
= b

a :
y x
x
y
Ahora bien, si por ejemplo

a b
= m(x)
y x
dependo solamente de x y b lo mismo, se puede suponer que dependera solamente de x, lo que da la
siguiente e ua ion diferen ial de primer orden lineal homogenea
b(x)0 = m(x)
que la sabemos resolver muy bien.
7.- Consideremos
yx2 dx x dy = 0
Veri amos que esta e ua ion no admite primitiva, por lo que es solu ion de la e ua ion diferen ial
(x2 + 1) = x0 ;
por lo tanto, resolviendo la e ua ion, obtenemos
2
e x =2
(x) =
x
La nueva e ua ion estara dada por
2
2
yxe x =2 dx e x =2 dy = 0;
obteniendo omo primitiva
2
g(x; y) = e x =2 y;
lo que da omo solu ion general
2
e x =2 y = C:
Utilizando reglas de al ulo de los diferen iales, la idea es luego de manipular la e ua ion diferen ial
original, llegar a una expresion de la forma
dg = 0;
lo que impli a que g(x; y; z) = . A manera de ilustra ion re ordemos que
d(u + v) = du + dv;
d(uv) = v du + u dv;
v du u dv
d(u=v) =
;
v2
d( os u) = sin u du;
Todas las reglas y formulas de deriva ion aprendidas en al ulo I, son validas para los diferen iales.
8.- Consideremos nuevamente la e ua ion
y dx x dy = 0;
dividiendola por y2, se obtiene
y dx x dy
= d(x=y) = 0;
y2
lo que da x=y = , es de ir y = x.
77
III.2 E ua iones a Derivadas Par iales
E ua iones A Derivadas Par iales de Primer Orden y Lineales
Son e ua iones de la forma
a(x; y)
u
+ b(x; y) u
x
y
= (x; y; u);
(EDP)
donde a; b; son fun iones ontinuas y u es la fun ion in ognita.

Suponemos que ono emos u y sea r : R ! R2 una urva parametrizada, es de ir r(t) = (x(t); y(t)).
Planteamos f = u r. La derivada de f esta dada por
f_(t) =
u
u
x_ + y:_
x
y
Comparando el lado izquierdo de la e ua ion (EDP) on esta ultima identidad obtenemos

(
x_ = a(x; y)
y_ = b(x; y)
(EC)
que las llamamos e ua iones ara tersti as de la e ua ion diferen ial a derivadas par iales de primer orden
lineal. Por otro lado, el lado dere ho de las e ua iones ara tersti as dan un ampo ve torial de ve tores tangentes llamado ve tores ara tersti os o dire iones ara tersti as. Las solu iones de las e ua iones
ara tersti as son urvas llamadas urvas ara tersti as.
Ahora bien, para proseguir on la solu ion del problema, ne esitamos una urva C de ondi iones ini iales.
Para asegurar existen ia y uni idad de las urvas ara tersti as las dire iones ara tersti as no deben ser
tangentes a la urva de ondi iones ini iales. Ver gura III.2.
Figura III.2.- Curvas, Dire iones Cara tersti as y Curva de Condi ion Ini ial
Consideremos la urva ara tersit a que pasa por (x0 ; y0) 2 C , es de ir la solu ion del problema a valor
ini ial
x_ = a(x; y)
y_ = b(x; y)
x(0) = x0 ; y(0) = y0 ;
78
Supongamos que ono emos tambien u sobre la urva de ondi iones ini iales, enton e el lado dere ho de la
e ua ion (EDP) se onvierte en
f_ = (x(t); y(t); f ) = g(t; f )
f (0) = u(x0 ; y0)
que es un problema a valor ini ial aso iado a una e ua ion diferen ial de primer orden.
Por onsiguiente hemos mostrado omo determinar la solu ion de un problema a valores ini iales aso iado a una e ua ion diferen ial a derivadas par iales de primer orden y lineal, utilizando el metodo de las
ara tersti as. Remar amos que u se determina ono iendo a que urva ara tersti a pertene e ada punto
(x; y).
III Complementos
Ejemplos
1.- Consideremos el problema
u
+ x u
x
y
u(x; 0) = 2:
y
u + u2
Los ve tores ara tersti os estan dados por el ampo

y
x
y la urva de ondi ion ini ial es el eje x. Observamos que las dire iones ara tersti as no son tangentes
a la urva de ondi iones ini iales, por lo que podemos ontinuar. Ahora onsideremos el problema a
valor ini ial
x_ = y
y_ = x
x(0) = x0 ; y(0) = 0:
Resolviendo este problema obtenemos la urva ara tersti a que pasa por (x0 ; 0)
1 1
x(t) = x0 ( et + e t ) = x0 osh t
2 2
1 t
1
y(t) = x0 ( et
2 2 e ) = x0 sinh t
El lado dere ho de la e ua ion a derivadas par iales se onvierte en
f_ = f + f 2
f (0) = 2
que es un problema a valor ini ial aso iado a una e ua ion de tipo de Bernouilli, uya solu ion es
2 :
f (t) =
1 et
El siguiente paso es determinar u(x; y). Sea (x; y). Tenemos
x = x0 osh t
y = x0 sinh t
despejando x0 y t obtenemos
t = tanh 1 (y=x):
Lo que da omo solu ion
2
u(x; y) =
tanh 1 (y=x)
2 e
79
La veri a ion de que u(x; y) dado una linea mas arriba, la dejamos omo ejer i io.
E ua iones a Derivadas Par iales de Segundo Orden
Las e ua iones a derivadas par iales de segundo orden lineales des riben fenomenos fsi os de diferente
naturaleza omo la propaga ion de ondas, las vibra iones, preoblemas de difusion y problemas de indole
esta ionario. A diferen ia del metodo de las ara tersti as empleado para resolver e ua iones a derivadas
par iales de primer orden, las e ua iones
La e ua ion de la Cuerda Vibrante.
Consideremos una uerda de longitud L sometida a una tension T .
Ver gura III.3. La uerda es sometida

a una deforma ion o impulso transversal lo que o asiona un movimiento transversal de la uerda denotado
por u(x; t).
Observamos que existen dos puntos (las extremidades) que permane en inmoviles durante todo el tiempo,
esto lo denotamos por
u(0; t) = u(L; t) = 0
(CB);
que llamamos ondi iones de borde (CB). Asimismo en el instante t = 0, la uerda tiene una posi ion, omo
tambien una velo idad ini ial,
u(x; 0) = f (x);
(CI)
u
(
x; 0) = g(x);
t
que son las ondi iones ini iales.
La e ua ion diferen ial que des ribe el movimiento de la uerda en el trans urso del tiempo esta dada
por
2
2u
2 u
=
v
;
(EDP)
t2
x2
p
donde v es la velo idad de propaga ion de las ondas, (v = T=), donde es la densidad lineal de la uerda.
111
000
000
111
000
111
000
111
000
111
111
000
000
111
000
111
000
111
000
111
u(x,y)
L
Figura III.3.- La Cuerda Vibrante
Para resolver este problema, apli amos el metodo de separa ion de variables que onsiste en plantear u(x; t) =
X (x)T (t). Por las ondi iones de borde (CB), dedu imos que
X (0) = X (L);
(CBX)
si queremos que T (t) no sea identi amente nula.
Introdu imos las fun iones X (x) y T (t) en la e ua ion (EDP), por lo que obtenemos
T(t)X (x) = v2 X 00 (x)T (t);
luego
T(t)
X 00 (x)
=
v2
:
T (t)
X (x)
Por onsiguiente, ada miembro de la ultima e ua ion debe ser obligatoriamente onstante; para ver eso,
jamos un t y movemos los x. Hemos obtenido dos e ua iones diferen iales:
X 00 (x)
T(t)
=
v2 ;
= :
T (t)
X (x)
80
Resolvemos primero, la e ua ion diferen ial para X (x). Para tenemos dos posibilidades:
i) = 2 0, en este aso la solu ion tendra la forma
X (x) = ex + e x :
Las ondi iones de borde (CBX) dan = = 0, por lo que si 0, la uni a solu ion posible es la nula.
Remar amos que nos interesa en ontrar solu iones no nulas.
ii) = 2, en este aso la solu ion sera de la forma
X (x) = os(x) + sin(x);
La primera ondi ion (CBS), X (0) = 0 ondu e a que = 0, por lo que nos queda omo solu ion
X (x) = sin(x):
Nos interesa que la solu ion que en ontremos no sea identi amente nula, por que exigimos omo ondi ion
sin(L) = 0:
Esto es posible, si y solamente si

= k ; k 2 Z:
L
Sin demostra ion a eptaremos que las solu iones que nos interesan son uando k = k=L on k 2 N.
La e ua ion diferen ial para T (t) es por onsiguiente
T_ = v2 2k T;
lo que da omo solu ion general
T (t) = (ak os(vk t) + bk sin(vk t));
de donde
u(x; t) = (ak os(vk t) + bk sin(vk t))(sin(k x):
Ahora bien, omo la e ua ion es lineal, ualquier suma de solu iones, es tambien solu ion, por lo tanto
III Complementos
u(x; t) =
1
X
k=1
(ak os(vk t) + bk sin(vk t))(sin(k x):
El siguiente paso es determinar los valores que toman ak y bk , para tal efe to, nos referimos a las ondi iones
ini ales (CI). Obtenemos
1
X
f (x) = u(x; 0) = ak sin(k x);
k=1
g(x) =
1
X
u
(
x; 0) = vk bk sin(k x):
t
k=1
Proposi ion III.2.1.- Se tiene:

Z L
0
= k=L.
donde k
8
<
0 si n = m;
sin(m x) sin(n x) dx = : L
2 si n 6= m:
81
Utilizando esta ultima proposi ion obtenemos

Z L
0
de donde
L
f (x) sin(n x)dx = an ;
an =
Para
bn =
Z L
Lvn
2 Z L f (x) sin(n x) dx:
g(x) sin(n x) dx =
Por lo tanto el problema tiene solu ion uni a.
vk
Z L
0
g(x) sin(n x) dx:
Ejemplo
2.- Consideremos una uerda de longitud , on v = 1, satisfa iendo las ondi iones ini iales siguientes:
u(x; 0) = 0;
(
1 si x 2 [=2 1=2; =2 + 1=2;
u
(
x; 0) =
t
0 sino.
Utilizando el metodo de separa ion de variables, obtenemos que k = k, para todo k. Llegamos a la
solu ion general de la forma
u(x; t) =
1
X
n=1
(an os(nt) + bn sin(nt)) sin(nx):
Puesto que u(x; 0) = 0, dedu imos inmediatamente que los an = 0. Cal ulemos los bn.
Z =2+1=2
2
sin(nx) dx = 2 ( os(n(=2 1=2)) os(n(=2 + 1=2))) :
b =
n
n =2 1=2
n2
En la gurar III.4 podemos observar la solu ion del movimiento de esta uerda en la gra a de u(x; t)
Figura III.4.- Movimiento de una uerda vibrante
82
III Complementos
La e ua ion del Calor
Consideremos el problema diferen ial

on ondi iones de borde
ondi iones ini iales
u
t
u
= x
2
2
(EDP)
u(0; t) = u(L; t) = 0
(CB)
u(x; 0) = f (x):
(CI)
Al igual que en el problema de la uerda vibrante, suponemos que

u(x; t) = X (x)T (t);
de la misma manera obtenemos las ondi iones de borde para X (x):
X (0) = X (L) = 0:
Remplazamos u(x; t) = X (x)T (t) en la e ua ion (EDP) y obtenemos
T_ (t)X (x) = T (t)X 00(x);
lo que da
X 00 (x)
T_ (t)
=

:
T (t)
X (x)
El mismo razonamiento que en el aso de la uerda vibrante ondu e a que
X (x) = sin(k x);
donde k = k=L.
Por lo tanto T es solu ion
T_ = 2k T;
de donde
2
T (t) = ak e k t
La solu ion general de la e ua ion diferen ial on las ondi iones de borde dadas es
u(x; t) =
1
X
k=1
2
ak e k t sin(k x):
Los ak estan dados por las ondi iones ini iales (CI), se obtienen de la misma manera que en el aso de la
uerda vibrante. Por lo tanto,
2 Z L f (x) sin(k x) dx:
ak =
L
0
III.3 Elementos de C
al ulo Varia ional
Mu hos problemas ingenieriles estan rela ionados a problemas de optimiza ion. Con el al ulo varia ional
estos problemas se pueden onvertir en problemas diferen iales. Para motivar formularemos algunos ejemplos
donde intervienen problemas de optimiza ion.
83
l ulo Varia ional

III.3 Elementos de Ca
Ejemplos
1.- Se quiere onstruir un tobogan que tiene una altura de a metros y una longitud horizontal de b metros.
Para el dise~no de di ho tobogan se supone que no existe fri ion entre la super ie del tobogan y el
objeto que se desliza. Ademas se pide que el tiempo de deslizamiento sea mnimo.
Para la solu ion del problema, podemos aso iar la super ie deslizable al grafo de una fun ion y(x) on
y(0) = a y y(b) = 0, ver gura III.5. Por otro lado apli ando los prin ipios y leyes de la fsi a, se tiene
que el objeto a deslizarse satisfa e
v2 (x) = 2g(a y(x));
donde v es la velo idad que tiene el objeto en el punto (x; y(x)).
Por otro lado, se tiene
ds
v= ;
dt
donde s representa la longitud de la urva. Re ordando al ulo I, sabemos que
de donde
ds
dx
= 1 + y02;
Z b ds
Z b p
02
dt
=dx
Tiempo = dx dt = ds dx = p 1 + y dx:
2g(a y)
0
0
0
dt
Por onsiguiente, el problema onsiste en en ontrar una fun ion y(x) on ondi iones de borde y(0) = a
e y(b) = 0 tal que
Z b p
1 + y02 dx ! min :
p
(a y)
0
Z b
2.- Sean a < b A > 0 y B > 0, se bus a una fun ion y(x) ontinuamente diferen iable on y(a) = A y
y(b) = B tal que el area de la super ie de revolu ion generada por y(x) sea minimal. Formulemos
matemati amente este problema. El area de la super ie de revolu ion generada por la urva y(x) esta
dada por
Z b
p
Area = 2y 1 + y02 dx;
a
Por onsiguiente el problema varia ional onsiste en en ontrar una fun ion y(x) on y(a) = A y y(b) = B,
tal que
Z b p
y 1 + y02 dx ! min :
a
En base a los dos ejemplos presentados, podemos enun iar la lase de problemas varia ionales que
resolveremos en este urso.
El problema, \En ontrar y(x) on y(a) = A y y(b) = B tal que
Z b
L(x; y(x); y0 (x)) dx
(PV)
donde L : R3 ! R es ontinua, es un problema varia ional. La fun ion L se llama fun ion de Lagrange o
Lagrangiano.
El siguiente paso, es onvertir el problema varia ional en un problema diferen ial, (estamos en el urso
de E ua iones Diferen iales). Para tal efe to, suponemos que y(x) es solu ion de (PV) y onsideramos una
fun ion h : [a; b ! R ontinuamente diferen iable on
h(a) = h(b) = 0;
84
Consideramos la fun ion g : R ! R, dada por
g() =
Z b
III Complementos
L(x; y(x) + h(x); y0 (x) + h0 (x)) dx:
Se puede mostrar que g es una fun ion diferen iable y sobre todo que g tiene un mnimo en = 0. Por
Cal ulo I g0(0) = 0. Por otro lado, se tiene

Z b
L(x; y(x); y0 (x))
L(x; y(x); y0 (x)) 0
g0 (0) =
h(x) +
h (x) dx
y0
y
Integrando por partes, el segundo termino de la integral se tiene

Z b
L(a; y(a); y0(a))
L(b; y(b); y0(b))
L(x; y(x); y0 (x))
h(b)
h(a)
h(x) dx +
g0 (0) =
a
y
y0
Por lo tanto
Z b
y0
=0
(x; y(x); y0 (x)) 0 h(x) dx:
Z b
L
{z
L(x; y(x); y0 (x))

h(x)
y
y0

L(x; y(x); y0 (x)) 0
dx = 0:
y0
Esta e ua ion es valida para ualquier fun ion h : [a; b ! R ontinuamente diferen iable, a ondi ion que
h(a) = h(b) = 0, lo que impli a, (sin demostra ion),

L(x; y(x); y0 (x))
L(x; y(x); y0 (x)) 0
= 0;
y
y0
si y(x) es solu ion de (PV). Hemos mostrado el teorema

Teorema III.3.1.- Una ondi ion ne esaria para que y(x) sea solu ion del problema varia ional \En ontrar
y(x) on y(a) = A e y(b) = B tal que
Z b
L(x; y(x); y0 (x)) dx;
a
es que y (x) sea solu ion de la e ua ion diferen ial de Eular-Lagrange

L(x; y; y0 )
L(x; y; y0 ) 0
= 0:
(EL)
h(x)
y
y0
Desarrollando la e ua ion de Euler-Lagrange, obtenemos la e ua ion

(EL)
Ly Lxy0 Lyy0 y0 Ly0y0 y00 = 0
que es una e ua ion diferen ial de segundo orden, donde Ly denota la derivada par ial de L respe to a y,
et .
Ahora bien, la e ua ion (EL) puede ser simpli ada, si se tiene informa ion adi ional respe to a L, por
ejemplo:
1.- L no depende expl itamente de y, en este aso Ly = 0 y por lo tanto
Ly0 (x; y0 ) = C:
2.- L no depende expl itamente de y0, en este aso L0y = 0 de donde
Ly (x; y) = 0;
85
l ulo Varia ional

es una e ua ion algebrai a, aso no interesante y en general sin solu ion.

3.- L no depende expl itamente de x, obtenemos por onsiguiente
Ly Lyy0 y0 Ly0 y0 y00 = 0;
multipli amos por y0,
y0 Ly (y0 )2 Lyy0 y0 y00 Ly0 y0 = 0;
efe tuando las agrupa iones y manipula iones algebrai as, se tiene
(y0 Ly + y00 L0y ) (y00L0y + (y0)2 Lyy0 + y0y00Ly0y0 ) = L0 (y0 L0y )0 = 0;
de donde
L y0 L0y = C:
(EL)
Remar a.- Una gran antidad de problemas varia ionales, presentan lagrangianos independientes de x, por
lo que es onveniente memorizar la ter era version de la e ua ion (EL).
Resolvamos los ejemplos pre edentes.
Ejemplos
3.- En ontrar y(x) on y(0) = a e y(b) = 0 tal que

Z b
1 + y02 dx ! min :
a y
Apli amos la e ua ion de Euler-Lagrange para este problema, se obtiene

s
1 + y02
a y
ha emos manipula iones algebra as y tenemos

por lo tanto
y02
(1 + y02)(a y)
1
(1 + y02 )(a
y)
= C;
= C;
(1 + y02 )(a y) = C:
Para resolver esta e ua ion de primer orden que no esta dada de manera expl ita, podemos despejar
y0 o bien parametrizar y, x respe to a donde es el angulo de las tangentes de la urva. Se tiene
y0 = tan , de donde
se 2 (a y) = C;
lo que da
y() = a + C os2 :
Por otro lado,
dx dy
os sin = 2C os2 ;
= d = 2C tan
d y0

integrando, obtenemos
1
x() = D C ( + sin 2):
2
Planteando # = 2, se tiene nalmente
x(#) = D C (# + sin #);
y(#) = a + C (1 + os #);
86
III Complementos
El siguiente paso es determinar C y D. La ondi ion y = a uando x = 0, la es ribimos

x(#0 ) = D C (#0 + sin #0 ) = 0;
y(#0 ) = a + C (1 + os #0 ) = a:
La segunda e ua ion es equivalente a es ribir C (1+ os #0) = 0, omo queremos solu iones no onstantes,
#0 = (la ele ion de no es arbitraria, proviene del gra o de la solu ion ver gura III.5. La
primera e ua ion dara D = C . Por onsiguiente las parametriza iones de la urva seran
x(#) = C ( + # + sin #);
y(#) = a + C (1 + os #);
La ondi ion y = 0, uando x = a, le es ribimos
x(#f ) = C ( + #f + sin #f ) = b;
y(#f ) = a + C (1 + os #f ) = 0;
Dividiendo estas dos ultimas e ua iones, se obtiene la ondi ion para #f
+ #f + sin #f b
1 + os #f = a :
Determinando #f se obtiene el valor de C . En la gura III.5 se observa la solu ion del problema on
a = 12, b = 15.
Figura III.5.- Traye toria mas rapida
Se podra expli itar la rela ion entre x e y, ver ejer i ios. Es importante observar que las solu iones que
en ontremos seran validas mientras que y(x) > 0.
87
l ulo Varia ional

4.- Resolvamos el ejemplo 2, es de ir: En ontrar y(x) 0 on y(a) = A > 0 y y(b) = B > 0, tal que
Z b p
1 + y02 dx ! min :
Dedu imos las e ua iones de Euler Lagrange para este problema, se obtiene
y
yy02
1 + y02
luego
1 + y02
= C;
= C:
1 + y02
Con el mismo razonamiento que el ejemplo 3, planteamos y0 = tan , lo que da
p
1 + tan2
de donde
= C;
y = C jse j :
Analizando el gra o del problema, ver gura III.6, =2 =2.
de donde
dx
d
Por otro lado,
tan
= C se tan
= C se ;
x = C ln(se + tan ) + D:
El siguiente paso, es eliminar de x e y, se tiene:
p
x D = C ln(y=C + (y=C )2 1);

p
x D
e C = y=C + (y=C )2 1);
(y=C )2 2(y=C )e x CD + e2 x CD = (y=C )2 1; y = C e

x D
):
y = C osh(
C
x D
C
+e
2
x D
C
Los valores de C y D estan determinados por las ondi iones de borde; es de ir:
a D
) = A;
C
b D
) = B:
C osh(
C
C osh(
En la gura III.6 se tiene la solu ion del problema para a = 0, b = 5, A = 5 y B = 5.
88
III Complementos
Figura III.6.- Curva generadora de super ie de area mnima

Problemas Isoperimetri os
Con fre uen ia se debe resolver problemas de optimiza ion en las que las solu iones deben satisfa er una
ondi ion adi ional, omo por ejemplo: que tenga ierta longitud dada de antemano, que el area de la
super ie de revolu ion generada por la urva tenga ierto valor, et .
Son problemas que pueden formularse omo: \En ontrar una fun ion y(x) on y(a) = A e y(b) tal que
Z b
Z b
F (x; y; y0 ) dx ! min;
G(x; y; y0 ) dx = 0;
donde F y G son fun iones ontinuas. F sera la fun ion objetivo y G es una fun ion restri ion. Remar amos
que G no es identi amente nula, solamente la integral se anula. Utilizando la idea de los multipli adores de
Lagrange, visto en Cal ulo II, el problema original se lo replantea, deniendo
L(x; y; y0 ; ) = F (x; y; y0 ) G(x; y; y0 )
y luego resolviendo: \En ontrar y(x) diferen iable on y(a) = A y y(b) = B tal que
Z b
L(x; y; y0 ; ) dx ! min :
Ejemplo
5.- Determinar la altura mnima que deben tener dos torres de alta tension sabiendo que la distan ia que
separa las torres es 100 m, la longitud del able es de 105m y que por uestiones de seguridad el punto
mas bajo del tendido debe estar a 3 m sobre la tierra.
Para resolver este problema, onsideramos la familia de fun iones y(x) on y( 50) = y(50) = 0 uya
longitud es 105; esta restri ion la es ribimos omo
Z 50 p
50
1 + y02 1:05
dx = 0:
89
La forma del able esta sujeta al prin ipio fsi o de menor energa; es de ir un uerpo al anza su estado
nal uando se en uentra en un nivel de mnima energa. Como la uni a energia que se puede onsiderar
en este problema es la poten ial, esta es igual a
l ulo Varia ional

Z b p
EP = y 1 + y02 dx:
a
donde es la densidad lineal del able. Por onsiguiente el problema onsistira en en ontrar una fun ion
y(x) on y( 50) = y(50) = 0 tal que
Z b p
Z 50 p
50
1 + y02 dx ! min;
1 + y02 1:05
dx = 0;
de donde la fun ion a apli ar las e ua iones de Euler-Lagrange es

p
L(x; y; y0 ; ) = y 1 + y02
p
1 + y02 1:05
Como L no depende de x, apli amos Euler-Lagrange en su variante respe tiva

!
p
p

yy02
y0 2
02 1 + y02 1:05 = :
p
p
y

1
+
y
1 + y02
1 + y02
Realizamos las simpli a iones algebra as que orresponden y se obtiene
y+
p
= ;
1 + y02
planteamos y0 = tan , de donde
y() = + se :
Por otro lado
dx se tan
= tan = se ;
d
por lo tanto
x() = ln(se + tan ) + d;
por simetra del problema dedu imos que d = 0. El siguiente paso es expresar y en fun ion de x, se
tiene:
ex= = se + tan ; (ex= se )2 = tan2 = se 2 1;
e2x= 2ex= se = 1;
x=
x=
se = e +2e = osh(x= );
de donde
y(x) = + osh(x= ):
La ondi ion de borde y(50) = 0, da + osh(50= ) = 0. La longitud del able se tradu e
Z 50 p
50
1 + y02 dx =
Z 50
50
se dx =
Z 50
50
osh(x= ) dx = 2 sinh(50= ) = 105:
De la ulima e ua ion obtenemos el valor de = 91:96 y la otra ondi ion da = 105.
90
III Complementos
Tenemos y(0) = 13:, de donde la altura de las torres respetando los tres metros reglamentarios, debe
ser por lo menos 16 m. En la gura III.7 se observa la forma que adquiere el able.
Figura III.7.- Cable entre dos torres

Problemas Varia ionales a Varias Variables
Los problemas varia ionales vistos mas arriba se referan a en ontrar una fun ion y que depende de una
sola variable x. Ahora bien, existe una antidad de problemas donde las fun iones in ognitas tienen varias
omponentes y dependen de una sola variable independiente, que en general es el tiempo t.
Los problemas varia ionales que estudiaremos en lo que sigue, on iernen a urvas, es de ir: \en ontrar
una apli a ion : [a; b ! Rn ( (t) = x(t)) on x(a) = A 2 Rn y x(b) = B 2 Rn tal que
Z b
L(t; x; x_ ) dt ! min;
on L : R Rn Rn ! R ontinuamente diferen iable."

De la misma manera on se pro edio en el aso de una sola variable, se plantea
g() =
Z b
L(t; x(t) + h(t); x_ + h_ (t));
donde x(t) es supuesta la solu ion del problema y h(t) ontinuamente diferen ible que satisfa e h(a) = 0 =
Se tiene que g tiene un mnimo en = 0, de donde g0(0) = 0. Finalmente obtenemos las formulas de
Euler-Lagrange, dejando el detalle de la determina ion al estudiante,
h(b).
sistema diferen ial de orden 2 y talla n.

Remar as.-
dLx_ 1
dt
dLx_ 2
dt
Lx1 = 0;
dLx_ n
dt
Lxn = 0;
Lx2 = 0;
...
1.- En el aso en que n > 1 no se puede simpli ar las e ua iones de Euler-Lagrange de la misma manera
que se hizo mas arriba.
2.- Para los problemas isoparametri os, la utiliza ion de los multipli adores de Lagrange es valida omo
antes.
91
l ulo Varia ional

Ejemplos
6.- Determinar la urva de longitud mnima que une dos puntos A y B del espa io. Por onsiguiente, el problema es determinar la urva (x(t); y(t); z(t)), 0 t 1 on (x(0); y(0); z(0)) = A y (x(1); y(1); z(1)) = B
tal que
Z 1p
x_ 2 + y_ 2 + z_ 2 dt:
p
L(t; x; y; z; x_ ; y;_ z_ ) = x_ + y_ + y_ , omo no depende ni de x, ni de y, ni de z , se tiene Lx = 0, Ly = 0

y Lz = 0 por lo que Lx_ = 1, Ly_ = 2, Lz_ = 3 , es de ir
2
x_
x_ 2 + y_ 2 + z_ 2
y_
p
x_ 2 + y_ 2 + z_ 2
z_
p
2
x_ + y_ 2 + z_ 2
p
= 1 ;
= 1 ;
= 1 :
Sin perder generalidad, podemos suponer que x_ 2 + y_ 2 + z_ 2 es onstante, ya que para toda urva existe
una parametriza ion de tales ara tersti as. Por onsiguiente
x_ = 1
y_ = 2
z_ = 3
Resolviendo este sistema, se obtiene
x = 1 t + d1
y = 2 t + d2
z = 3 t + d3
que es la parametriza ion de una re ta, los valores de i y di estan determinados por A y B.
7.- Determinar la urva que esta por en ima del eje x, que une los punto ( 1; 0) y (1; 0) de manera que el
area limitada por la urva y el eje x sea A y la longitud de la urva sea mnima.
Por onsiguiente tenemos un problema isoperimetri o. La formula de Green permite evaluar el area
en errada por la urva. Para tal efe to, el sentido de la urva debe ser ontrario a las mane illas del
reloj, por lo que debemos en ontrar una urva : [0; 1 ! R2 on (0) = (1; 0) y (1) = (0; 1). El area
de la region a hurada sera por lo tanto
Z
1 y dx 1 x dx = Z 1 1 (yx_ xy_ ) dt:
Area =
2
0 2
2
La longitud de la urva esta dada por Z 1
p
x_ 2 + y_ 2 dt:
0
92
III Complementos
De donde el problema a resolver, sera

Z 1 p
0
x_ 2 + y_ 2 (yx_ xy_ ) dt:
Apli ando las e ua iones de Euler Lagrange, obtenemos el sistema

!
x_
y =dt = y;_
x_ 2 + y_ 2
!
y_
d p 2 2 + x =dt = x;
_
x_ + y_
d
Como en el ejemplo pre edente, suponemos que x_ 2 + y_ 2 = , de donde

x
= 2y_

y
= 2x_

La onstante la in orporamos al lado dere ho de ada una de las e ua iones, obteniendo

d(x_ y)
= 0;
dt
d(y_ + x)
dt
de donde
= 0;
x_ = y + d1
y_ = x + d2
La traye toria de esta e ua ion diferen ial es una ir unferen ia, la veri a ion la dejamos omo ejer i io.
Por lo tanto la urva que une los dos puntos es un ar o de ir unferen ia. La determina ion de di ha
ir unferen ia es un simple ejer i io de geometra.
Problemas Geodesi os
Mas arriba hemos abordado problemas varia ionales en las que las fun iones a en ontrar no tenan restri ion
alguna, o se ono ia algun valor de esta determinada por medio de una integral. Ahora bien, en algunas
situa iones es ne esario que las urvas se en uentren dentro un onjunto dado omo ser una super ie.
La formula ion de este tipo de problemas esta dada por, en ontrar una fun tion : [0; 1 ! Rn on
(0) = A y (0) = B tal que
Z 1
F (t; x; x_ ) dt ! min
0
x2M
Rn :
es en general una super ie si n = 3.

Por onsiguiente, nos limitaremos al aso en que M = fx 2 Rn jg(x) = 0g, donde g : Rn ! R es ontinuamente diferen iable. Utilizando multipli adores de Lagrange, el problema se onvierte en: \En ontrar
: [0; 1 ! Rn on x(0) = A y x(1) = B tal que
Z 1
0
remar ando que depende de x".
(F (t; x; x_ )
(x)g(x)) dt ! min
93
l ulo Varia ional

Las e ua iones de Euler-Lagrange son por lo tanto:

dFx_ 1
dt
dFx_ 2
dt
(Fx1
(Fx2
dFx_ n
dt
(Fxn
(x)gx1 ) = 0
(x)gx2 ) = 0
...
(x)gxn ) = 0
Ahora bien, en los problemas usuales F depende solamente de x_ , por lo que las e ua iones on esta supo ion
son iguales a:
dFx_ 1
+ (x)gx1 = 0
dt
dFx_ 2
+ (x)gx2 = 0
dt
...
dFx_ n
dt
+ (x)gxn = 0:
Volvamos a nuestro problema de geodesi as.

Deni ion III.3.2.- Sea M R un onjunto (super ie), A; B 2 M . Una geodesi a en M que une A y B
es una urva en M que une A y B de longitud minimal.
En el aso de las geodesi as, se tiene
q
F (x_ ) = x_ 21 + + x_ 2n :
p
Ahora bien, podemos elegir una parametriza ion de manera que x_ 21 + + x_ 2n sea onstante, por lo que
las e ua iones de Euler-Lagrange daran
x1 = (x)gx1 ;
...
xn = (x)gxn ;
Por lo tanto, dedu imos que x y grad g son linealmente dependientes en todo x 2 M .
Ejemplo
8.- Determinemos las geodesi as de una ir unferen ia. Sin perder generalidad podemos suponer que la
esfera esta entrada en el origen. Por lo tanto
g(x; y; z ) R2 = 0:
Ahora bien por lo he ho anteriormente (x; y; z) y grad g(x; y; z) = (2x; 2y; 2z) son linealmente independientes y esto su ede si y solamente si
(x; y; z) grad g(x; y; z) = 0:
Obtenemos las siguientes e ua iones
yz zy = 0;
xz zx = 0;
yx xy = 0:
94
III Complementos
Luego se tiene
de donde
yz zy = 0 = (yz

_ yz_ );
xz zx = 0 = (xz
_ xz_ );
xz xy = 0 = (yx
_ yx_ );
yz
_ yz_ = 1
xz
_ xz_ = 2
yx
_ yx_ = 3
Multipli ando la primera e ua ion por x, la segunda por y, la ter era por z, luego adi ionando se
tiene
1 x 2 y + 3 z = 0:
Por lo tanto, las geodesi as son las inter iones de la esfera on planos que pasan por el origen.
III.4 Ejer i ios
1.- Resolver las siguientes e ua iones diferen iales:

a) x dx + y dy = 0;
b) x dx y dy = 0;
) y dx + x dy = 0;
d) y dx x dy = 0;
e) x dy y dx (1 x2 ) dx = 0;
f) (x2 y) dx x dy = 0;
g) y(x 2y) dx x2 dy = 0;
h) (x + y os x) dx + sin y dy = 0.
2.- Resolver
y dx x dy
(x + y)2 + dy = dx
utilizando fa tor integrante y luego mediante manipula iones on los diferen iales.
3.- Resolver
4y2 2x2 dx + 8y2 x2 dy = 0
4xy2 x3
4y3 x2 y
utilizando si es ne esario un fa tor integrante y luego resolviendo omo una e ua ion diferen ial ordinaria.
M (x; y) dx + N (x; y) dy = 0:
Probar que si (My Nx)=(Ny Mx) es una fun ion g(z) donde z = xy, enton es
= eG(z)
es un fa tor integrante, donde G0(z) = g(z).
b) >Que forma tendra el fa tor integrante si z = x + y.
5.- Resolver las siguientes e ua iones:
a) (y ln y 2xy) dx + (x + y) dx = 0;
b) ex dx + (ex ot y + 2y s y) dx = 0;
95
III.4 Ejer i ios
) x dy + x dx = pxy dy.
6.- Resolver las siguientes e ua iones on derivadas par iales de primer orden:
a) u
+ u = x2 + y2 on u(x; 0) = sin x;
x y
u
b) u
= u2 on u(x; 0) = ex;
x y
1 u = 0 on u(0; y) = y.
) u
x x2 y
u
+
= 1 + u2 on u(x; 0) = 1.
d) u
x y
7.- Resolver la e ua ion de la uerda vibrante uya e ua ion es tu2 =
dadas por: u(0; t) = u(; t) = 0 y ondi iones ini iales dadas por:
u(x; 0) = sin x + 3 sin 3x
u
(x; 0) = 0
x
2
2u
x2
on ondi iones de frontera
u
8.- Resolver la e ua ion de difusion uya e ua ion es u
on ondi iones de frontera dadas por:
=
t
x2
u(0; t) = u(; t) = 0 y ondi iones ini iales dadas por:
u(x; 0) = sin2x + 3 sin5x
2
2u 2u
+ y2 on ondi iones de
9.- Resolver la e ua ion de la membrana uadrada uya e ua ion es tu2 = x
2
frontera dadas por: u(0; y; t) = u(; y; t) = u(x; 0; t) = u(x; ; t) = 0 y ondi iones ini iales dadas por:
u(x; y; 0) = sin 2x sin y + 3 sin4x sin3y
u
(x; y; 0) = 0
x
10.- Resolver la e ua ion de difusion uya e ua ion es
2
u
t
= 4u;
on ondi iones de frontera dadas por

u = 0 sobre los bordes del uadradado de lado :
Condi iones ini iales dadas por:
u(x; y; 0) = sin x sin y + sin 2x sin 3y:
u
11.- Resolver la e ua ion de difusion uya e ua ion es u
on ondi iones de frontera dadas por:
=
t
x2
u(0; t) = 0; u(; t) = 100 y ondi iones ini iales dadas por:
u(x; 0) = 0
12.- Hallar la urva on y( 1) = y(1) = 1 que genera una super ie de revolu ion de area minimal.
b) Sabiendo que el area de la super ie es minimal y el volumen es igual a 4.
13.- Hallar la e ua ion que des ribe una uerda de longitud 3 sometida a la a ion de la gravedad, sabiendo
que sus extremidades tienen las oordenadas ( 1; 0) y (1; 0).
14.- Mostrar que la urva que une dos puntos del plano, es un segmento de re ta.
15.- DeterminarRlas1 p
fun iones extremales siguientes:
a) 1(y) = R 1 y(1 + y02) dx, on y( 1) = y(1) = b > 0;
b) 2(y) = ab 1+y0y22 dx on y(a) = A, y(b) = B;
2
96
III Complementos
) 3(y) = RRab (xp y)2 dx on y(a) = A y y(b) = B;

d) 4(y) = R ab p
x 1 + y02 dx on y(a) = A y y(b) = B ;
b
e) 5(y) = a x(1 + y02 dx on y(0) = A y y(b) = B, b > 0.
16.- Si P y Q son dos puntos de un plano, la longitud de una urva que los une es, en oordenadas polares,
Z Qp
dr2 + r2 d2
Hallar la e ua ion polar de una re ta ha iendo que esta integral sea mnima.
a) on omo variable independiente;
b) on r omo la variable independiente.

Textocfinal II2007

Cargado por

Copyright:

Formatos disponibles

Textocfinal II2007

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Textocfinal II2007

Cargado por

Copyright:

Formatos disponibles

Universidad Mayor de San Simon

Fa ultad de Cien ias y Te nologa

Cal ulo III

la derivada (ordinaria) de orden de f .

I E ua iones Diferen iales Ordinarias

El estudio de las e ua iones diferen iales pasa por lo tanto, por:

Consideramos la e ua ion diferen ial (ordinaria y expl ita)

uyas solu iones estan dadas por

Problemas Diferen iales

I E ua iones Diferen iales Ordinarias

Figura I.1.- Problemas Diferen iales de Vigas

I.2 E ua iones Diferen iales de Primer Orden

donde a; b : I ! R ontinua e I  R intervalo.

! R es ontinua, esta dada por

I.2 E ua iones Diferen iales de Primer Orden

derivemos , lo ual ondu e a

(f + g)(x) = f (x) + g(x);

( ' + )0 (x) = '0 (x) + 0(x)

I E ua iones Diferen iales Ordinarias

Demostra ion.- Mostraremos que el problema diferen ial

1.- La solu ion general de y0 = os(x)y, es

3.- La solu ion del problema y0 = a(x)y, y(x0 ) = 0; es

Consideremos la e ua ion lineal de primer orden

Proposi ion I.2.6.- Sea

Demostra ion.- En efe to

(' + )0 (x) = '0(x) + 0(x)

I.2 E ua iones Diferen iales de Primer Orden

tiene solu ion uni a.

Demostra ion.- Ejer i io.

I E ua iones Diferen iales Ordinarias

omo solu ion parti ular.

Cuando no es posible adivinar las solu iones parti ulares, se tiene el

Varia ion de Constantes.

Remar a Varia ion de onstantes es un metodo que nun a falla,

E ua iones de Tipo Bernouilli

Son e ua iones diferen iales de primer orden de la forma

1.- La expresion y , on 2 R, esta de nida por

Teorema I.2.9.- Sean ' y

dos fun iones, si

es solu ion de la e ua ion de tipo Bernouilli

Bernouilli, introdu ir la variable

8.- Determinemos la solu ion general de la e ua ion diferen ial

I E ua iones Diferen iales Ordinarias

y por onsiguiente la solu ion general para la e ua ion de Bernouilli es

tiene solu ion uni a, si esta existe.

Demostra ion.- Ejer i io.

I.2 E ua iones Diferen iales de Primer Orden

donde F y G son primitivas de f y g respe tivamente y G

denota la fun ion inversa de G.

Cuando se tiene que resolver una e ua ion separable olo ar en un lado de la

Regla Memote ni a.Ejemplo

9.- Resolvamos la e ua ion

Demostra ion.- Ejer i io.

E ua iones de tipo Homogeneo

donde f es una fun ion ontinua y x 6= 0.

I E ua iones Diferen iales Ordinarias

I.3 Problemas de Existen ia y Uni idad

11.- La e ua ion diferen ial

I E ua iones Diferen iales Ordinarias

Universidad Mayor de San Simon

Fa ultad de Cien ias y Te nologa

Cal ulo III

Consideramos la e ua ion diferen ial (ordinaria y expl ita)

uyas solu iones estan dadas por

donde a; b : I ! R ontinua e I R intervalo.

! R es ontinua, esta dada por

(' + )0 (x) = '0 (x) + 0(x)

Demostra ion.- Mostraremos que el problema diferen ial

1.- La solu ion general de y0 = os(x)y, es

3.- La solu ion del problema y0 = a(x)y, y(x0 ) = 0; es

Consideremos la e ua ion lineal de primer orden

Proposi ion I.2.6.- Sea

Demostra ion.- En efe to

tiene solu ion uni a.

Demostra ion.- Ejer i io.

omo solu ion parti ular.

Varia ion de Constantes.

Remar a Varia ion de onstantes es un metodo que nun a falla,

1.- La expresion y, on 2 R, esta denida por

es solu ion de la e ua ion de tipo Bernouilli

8.- Determinemos la solu ion general de la e ua ion diferen ial

y por onsiguiente la solu ion general para la e ua ion de Bernouilli es

tiene solu ion uni a, si esta existe.

Demostra ion.- Ejer i io.

denota la fun ion inversa de G.

Cuando se tiene que resolver una e ua ion separable olo ar en un lado de la

Regla Memote ni a.Ejemplo

9.- Resolvamos la e ua ion

Demostra ion.- Ejer i io.

E ua iones de tipo Homogeneo

donde f es una fun ion ontinua y x 6= 0.

11.- La e ua ion diferen ial

! R ontinua satisfa e una ondi ion de Lips htiz y x0 2 I ,

admite una uni a solu ion.

admite solu ion uni a.

Solu iones Aproximadas de una E ua ion Diferen ial

Figura I.3.- Ilustra ion del metodo de Euler

El lado dere ho de la e ua ion diferen ial

Figura I.4.- Ilustra ion del metodo de las Iso linas

Consideremos la e ua ion diferen ial de orden n, en forma expl ita,

La redu ion de orden es posible en los siguientes dos asos:

1.- Resolvamos la e ua ion

on x0 en el intervalo de deni ion de las ai , tiene solu ion uni a.

Demostra ion.- Sea '1 ; : : : ; 'n

un sistema fundamental de la e ua ion diferen ial (LH). Por onsiguiente,

Deni ion I.4.4.- Se llama a

Demostra ion.- Conse uen ia dire ta del Algebra Lineal.

3.- Consideremos la e ua ion diferen ial lineal homogenea de segundo orden

Determina ion de Sistemas Fundamentales

4.- Consideremos la e ua ion de segundo orden

2.- La e ua ion diferen ial lineal es a oe ientes onstantes.

Demostra ion.- El punto i) ejer i io. Demostremos el punto ii). Se tiene

+ + a1r + a0 es el polinomio ara tersti o

Demostra ion.- Efe tivamente ada raiz aporta on una

5.- Consideremos la e ua ion diferen ial de orden 4.

El polinomio ara tersti o esta dado por

Polinomio ara tersti o on rai es multiples

skiano veri ara la independen ia lineal.

3.- Si laa e ua ion lineal homogenea es de la forma

Proposi ion I.4.16.- Supongamos que f; g son analiti as en el punto x0 , enton es

Consideremos la e ua ion lineal homogenea de orden n

7.- Determinemos la solu ion general de