EDO Notas de Clase V4
EDO Notas de Clase V4
EDO Notas de Clase V4
ECUACIONES DIFERENCIALES
ORDINARIAS
Tcnicas de Resolucin
Edixon M. Rojas
[email protected]
Universidad Nacional de Colombia
Bogot-Colombia
Copyright
c 2016 Luz Marina Moya & Edixon Rojas
N OTAS DE CURSO
1 Nociones Fundamentales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1 Definiciones y Ejemplos 8
1.2 Solucin de una Ecuacin Diferencial 13
1.2.1 Solucin General y Particular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.2.2 Problema de Valor Inicial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Nociones Fundamentales
Definicin 1.1 Un objeto es una unidad elemental sobre la que se pueden hacer observa-
ciones y cuya estructura interna no existe o se puede ignorar. Un sistema es una coleccin
de objetos relacionados entre si. Una descripcin es una representacin de un fenmeno
por medio de un lenguaje, en nuestro caso, lenguaje matemtico, explicando sus distintas
etapas, partes y/o cualidades. Un modelo es una descripcin matemtica de un sistema.
En las ciencias los modelos matemticos son usados para buscar comprender fenmenos
naturales, sistemas reales biolgicos o fsicos, a travs de la identificacin de las variables
que intervienen en la situacin o problema y a continuacin cuantificando, generalmente
mediante el uso de ecuaciones, las relaciones existentes entre las variables identificadas.
Este procedimiento permite predecir el comportamiento futuro aproximado del sistema bajo
estudio y, de ser posible, la manipulacin del mismo.
Algunos problemas que son de especial inters en reas de las ciencias e ingeniera, son
problemas dinmicos, que involucran variables independientes que cambian respecto a otras
variables dependientes y es precisamente el concepto de derivada, lo que permite describir
cuantitativamente las relaciones entre las variables involucradas en la descripcin del pro-
blema en estudio. Como consecuencia, las ecuaciones diferenciales aparecen constituyendo
ejemplos de modelos matemticos.
En apariencia las ecuaciones diferenciales parecen ser algo nuevo para nosotros, sin embargo
ya nos hemos encontrado con el problema de resolver ecuaciones diferenciales en reas
como la fsica o de la misma matemtica. Por ejemplo, al determinar la posicin de un
cuerpo en movimiento o en la bsqueda de las antiderivadas de una funcin dada.
8 Nociones Fundamentales
que involucra una funcin desconocida f en las variables x1 , . . . , xt y algunas de sus deriva-
das, donde F denota una funcin de n = n1 + + nt variables.
A continuacin presentaremos algunos ejemplos de ecuaciones diferenciales que aparecen
en algunas reas de las ciencias.
F = ma.
donde la fuerza F es una funcin que depende del tiempo t, de la posicin x(t) y de la
velocidad dx/dt.
Ejemplo 1.5 Las ecuaciones diferenciales dadas en (1.2), (1.3) y (1.6) asociadas a los
modelos descritos en los Ejemplos 1.1, 1.2 y 1.3, son ejemplos de ecuaciones diferenciales
ordinarias.
x x
+ = 3,
y t
Ejemplo 1.7 Otro ejemplo de ecuacin diferencial parcial viene dado por
2u 2
2 u u
2
= c 2
p qu,
t x t
donde c es una constante. Esta ecuacin es llamada ecuacin del telgrafo, dado que
apareci por primer vz al tratar de determinar la distribucin de la corriente y el voltaje a
lo largo de las lineas terrestres de los telgrafos.
Ejemplo 1.8 El flujo de calor en un slido est gobernado por la ecuacin diferencial
parcial
T T T
x + y = c ,
x x y y t
donde x y y son las conductividades trmicas del solido en las direcciones x y y, T es la
temperatura, la densidad del slido y c una constante.
Ejemplo 1.11 a) x3 y000 + 2x2 y00 + 3xy0 + 4x = 0 es una ecuacin diferencial ordinaria
de tercer orden y primer grado.
b) (y00 )2 = sin(x) + yex es una ecuacin diferencial ordinaria de segundo orden y de
grado 2.
Obs Anlogo al caso de las ecuaciones algebraicas, en general, no tiene sentido hablar del
grado de una EDO. Por ejemplo, si en una ecuacin diferencial, alguna de las derivadas
de la funcin incgnita aparece con un exponente que no es un nmero natural o
tambien cuando aparecen haciendo parte del argumento de una funcin trascendente,
como por ejemplo sin, ln, exp o cualquier otra de este estilo.
Ejemplo 1.12 La EDO x3 y000 + 2x2 y00 + 3xy0 + 4x = 0 es una ecuacin diferencial ordinaria
de tercer orden y grado uno.
Ejemplo 1.13 La ecuacin diferencial (y00 )2 = sin(x) + yex es una ecuacin diferencial
ordinaria cuadrtica de segundo orden.
d2 y
dy
Ejemplo 1.14 La ecuacin diferencial e + dx
dx2 = sin(x), es un ejemplo de ecuacin
diferencial de segundo orden no lineal. Ntese que este es un ejemplo en el cual no tiene
sentido hablar del grado de la ecuacin.
12 Nociones Fundamentales
Ejemplo 1.15 La ecuacin de Van der Pol, es una ecuacin diferencial ordinaria cuadrtica
de segundo orden
d2x dx
2
(1 x2 ) + x = 0,
dt dt
donde es una constante, usada para describir el comportamiento de un oscilador no
conservativo con amortiguamiento no lineal.
u 3u u
+ 6 f (t)u = g(t, u),
t x x
llamada ecuacin de Korteweig-de Vries (KdV), usada para modelar el comportamiento de
ondas en aguas poco profundas.
dy d3y
Ejemplo 1.17 Las ecuaciones diferenciales dx = xy2 y dx3
+ y2 = 0, tambin son ejemplos
de ecuaciones diferenciales ordinarias no lineales.
Ejercicios
1. Clasificar las siguientes ecuaciones diferenciales de acuerdo a su tipo, indicar las
variables dependiente e independiente y proporcionar el orden de la ecuacin. Si la
ecuacin diferencial es ordinaria, determinar si es lineal o no lineal.
a) 3y000 4xy0 + x3 y 4 = 0.
b) ln(x) 3xy00 4xy = sin(x).
c) 3xy00 3y0 + 4y y2 = 0.
d 2 y 3 dy
d) dx 2 cos(x) + 4xy dx + y2 ex = 0.
3 3
e) 3 xy3 z3y = 0.
f) ln(x) + 6xy00 4xy = cos(x).
d3y dy 4
g) x dx 3 ( dx ) + y = 0.
h) (sin( ))y00 (cos( ))y0 = 2.
Obs Al usar la expresin resolver una ecuacin diferencial o hallar la solucin de una
ecuacin diferencial queremos significar hallar una funcin la cual es una solucin
de la ecuacin diferencial. De forma anloga, cuando nos referimos a cierta ecuacin
como solucin de una ecuacin diferencial, queremos decir que la funcin definida por
la ecuacin es solucin. As, si la ecuacin no define una funcin, entonces esta no es
una solucin de una ecuacin diferencial a pesar de que siguiendo el procedimiento
formal podamos comprobar que la ecuacin satisface la ecuacin diferencial.
p
Ejemplo 1.19 Comprobar que la frmula y = (1 + x2 ) satisface la ecuacin diferencial
x + yy0 = 0, (1.11)
sin embargo no es una solucin de esta.
p x
Solucin. Derivando y = (1 + x2 ) tenemos que y0 = .
Sustituyendo en la
(1+x2 )
p
ecuacin diferencial (1.11) vemos que se satisface la identidad, sin embargo y = (1 + x2 )
no define una funcin y por lo tanto no es una solucin de la ecuacin diferencial (1.11).
Ejemplo 1.21 Verificar que la funcin f (x) = x + ex , con dominio toda la recta real R,
define una solucin explcita de la ecuacin diferencial y0 + y = x + 1.
Solucin. Calculando la primera derivada de f (x) y obtenemos f 0 (x) = 1 ex la cual est
definida en todo R. Luego
f 0 (x) + f (x) = 1 ex + x + ex = 1 + x.
Es decir, se satisface la igualdad para cualquier x R.
Ejemplo 1.22 Verificar que para cualquier par de nmeros reales c1 y c2 , la funcin
(x) = c1 e2x + c2 e3x
00 (x) 5 0 (x) + 6 (x) = (4c1 e2x + 9c2 e3x ) 5(2c1 e2x + 3c2 e3x ) + 6(c1 e2x + c2 e3x )
= e2x (4c1 10c1 + 6c1 ) + e3x (9c2 15c2 + 6c2 )
= 0.
Ejemplo 1.23 Verificar que la funcin y = tan(x) x, definida para todo x 6= (2n + 1) 2 ,
n = 0, 1, 2, . . . define una solucin explcita de la ecuacin diferencial y0 = (x + y)2 .
Solucin. Calculando la primera derivada de y obtenemos y0 = sec2 (x) 1 = tan2 (x).
Sustituyendo y y y0 obtenemos la identidad
Se tiene as la igualdad para cualquier x en cada uno de los intervalos en donde la funcin y
junto con su derivada estn definidas.
Notemos que la ecuacin diferencial est definida en todo R, sin embargo la solucin
y = tan(x) x solo est definida en ciertos intervalos contenidos en R.
Obs El ejemplo anterior muestra que el dominio de definicin de una ecuacin diferencial
puede ser mayor al dominio de una solucin. Tambin es posible que una funcin que
est definida en un intervalo sea solucin de una ecuacin diferencial solo en una parte
de este intervalo. Por ejemplo, la funcin f (x) = |x| est definida en todo x R y no
es diferenciable en x = 0. Esta funcin satisface la ecuacin diferencial f 0 (x) = 1, en
el intervalo x > 0 y f 0 (x) = 1 en el intervalo x < 0. Sin embargo f (x) no es solucin
de ninguna ecuacin diferencial en cualquier intervalo que contenga al punto x = 0.
Definicin 1.8 Solucin Implcita. Una relacin g(x, y) = 0 define una solucin
implcita de la ecuacin diferencial ordinaria F(x, y, y0 , y00 , , y(n) ) = 0 en el intervalo I, si
se satisfacen las siguientes condiciones:
1. La relacin g(x, y) = 0 define implcitamente a y como funcin de x sobre el intervalo
I; es decir, existe una funcin (x) definida sobre I, tal que para todo x I se verifica
g(x, (x)) = 0.
para todo x I.
16 Nociones Fundamentales
Obs Ntese que g(x, y) = 0 es una ecuacin y una ecuacin nunca es una solucin de una
ecuacin diferencial ya que solo las funciones pueden ser soluciones de estas. Lo que
queremos significar al decir que la relacin g(x, y) = 0 define una solucin implcita de
una ecuacin diferencial es que la funcin y = (x) definida por la relacin g(x, y) = 0
es la solucin.
Obs El procedimiento estndar para mostrar que una relacin es una solucin de una
ecuacin diferencial es derivar implcitamente la relacin y sustituirla en la ecuacin.
Sin embargo, se debe tener cuidado con el dominio de la solucin al realizar este
procedimiento ya que puede ocurrir, por ejemplo, que afirmemos que x2 + y2 = 0
es una solucin implcita de x + yy0 = 0, pues al derivar implcitamente x2 + y2 = 0
y sustituir en x + yy0 = 0 obtenemos la identidad, pero x2 + y2 = 0 no define a y
implcitamente como una funcin de x en ningn intervalo ya que solo el punto (0, 0)
satisface esta frmula. As, concluir que x2 + y2 = 0 es una solucin implcita de
x + yy0 = 0 no tiene sentido.
2
Ejemplo 1.26 Demostrar que la familia de funciones y = ke2x definidas en todo la recta
real, hacen parte de la solucin general de la ecuacin diferencial
dy
= 4xy.
dx
2
Solucin. Derivando y obtenemos y0 = 4kxe2x , es decir, y0 = 4xy.
2 /2
Ejemplo 1.27 Demostrar que la funcin y = 3 3ex es una solucin particular de la
ecuacin diferencial
dy
+ xy = 3x.
dx
18 Nociones Fundamentales
2 /2
Solucin. Derivando y obtenemos y0 = 3xex . Luego,
dy x2 x2 x2 x2
+ xy = 3xe 2 + x(3 3e 2 ) = 3xe 2 + 3x 3xe 2 = 3x.
dx
Obs El uso de los trminos solucin general y solucin singular como aqu estn definidos
es un poco controversial dado que para algunos autores estos no son del todo correctos.
Por ejemplo, consideremos la ecuacin diferencial de primer orden y0 = 2y3/2 . La
1
familia de funciones de la forma y = (x+c) 2 forman parte de la solucin general de
c 2
Pero, por otro lado, cualquier funcin de la forma y = (cx+1) 2 tambin es solucin de la
Ejemplo 1.30 Verificar que y = xex + cx, donde c es una constante arbitraria, proporciona
una familia a un parmetro de soluciones explcitas de la ecuacin
dy y
= xex .
dx x
Graficar varias de estas curvas solucin y determinar la solucin particular que satisfaga la
condicin inicial y(1) = e 1.
Solucin. Derivando y = xex + cx obtenemos y0 = ex + xex + c. Ahora,
dy y (xex + cx)
= ex + xex + c = ex + xex + c ex c = xex .
dx x x
En la Figura 1.1 se bosquejan algunas curvas solucin de y = xex + cx.
Por ltimo, se determina c de modo que satisfaga la condicin inicial y(1) = e 1. Al hacer
x = 1 y y = e 1 en la ecuacin y = xex + cx se tiene
e 1 = e + c,
Ejercicios
1. Comprobar que la funcin o familia de funciones indicada es una solucin explicita o
implcita de la ecuacin diferencial dada.
a) y0 3y = 0; y = 2e3x .
b) y00 y0 + y = 0; y = c1 ex/2 cos(x 3/2) + c2 ex/2 sin(x 3/2).
c) y0 cos(t) + y sin(t) = 0; y = sin(t) + cos(t).
d) xydx + (x2 /2 + y)dy = 0; x2 y + y2 = k.
e) y00 + y = cot(x); y = c1 cos(x) + c2 sin(x) + sin(x) ln[(1 cos(x))/ sin(x)].
f ) y00 + y0 12y = 0; y = c1 e4x + c2 e
3x .
3
g) y000 3y00 + 3y0 y ex = 0; y = ex c1 + c2 x + c3 x2 + x6 .
h) y00 + 2y0 + y = 0; ey = c1 x + c2 .
i) y00 + (y0 )2 = 0; y = ln |x + c1 | + c2 .
000 00 2x 2 x 1 9 cos(2x)7 sin(2x)
x
j) y y e sin (x) = 0; y = c1 + c2 e + c3 e + 12 + 520 e2x .
k) x3 y000 + 2x2 y00 xy0 + y = 12x2 ; y = c1 x + c2 x ln |x| + 4x2 , x > 0.
l) y2 1 (2y + xy)y0 = 0; y2 1 = (x + 2)2 .
dy
m) exy + eyx dx = 0; e2y + e2x = 1.
n) xy0 = 2y; y = x2 .
) y00 = 1 2 ; y = arcsin(x) + 1 x2 .
1x
o) y0 = 13 y; y = 2ex/3 .
dr
p) cos( ) d 2r sin( ) = 0; r = a sec2 ( ).
2. Demostrar que la ecuacin diferencial
dy
+ |y| + 3 = 0
dx
no tiene soluciones.
Como vimos, una ecuacin de diferencial de primer orden es una relacin de la forma
Cuando esta relacin es lineal en las variables t, y e y0 ; es decir, que se pueda reescribir
como
a1 (t)y0 + a2 (t)y + g(t) = 0, a1 (t) 6 0
se dice que la ecuacin diferencial de primer orden es lineal. En caso contrario se dice no
lineal.
En este captulo se describen varios mtodos de resolucin para ciertos tipos de ecuaciones
diferenciales de primer orden. De manera ms precisa, se estudiarn tcnicas de resolucin
de ecuaciones de primer orden no lineales llamadas ecuaciones de variable separable
y ecuaciones exactas. Se mostrarn estrategias para hallar algunas funciones auxiliares
llamadas factores integrantes las cuales permitirn reducir algunas ecuaciones no lineales a
ecuaciones exactas. Tambin se resolvern ecuaciones diferenciales lineales y finalmente,
mediante sustituciones convenientes, se podrn reducir algunas ecuaciones diferenciales
tales como las llamadas ecuaciones Homogneas, de Bernoulli, de Ricatti, con coeficien-
tes lineales a ecuaciones lineales y a otras ecuaciones estudiadas previamente. Primero
presentaremos el teorema de existencia y unicidad de soluciones de ecuaciones de primer
orden.
Teorema 2.1 Existencia y unicidad de soluciones para EDO de primer orden
f
orden. Sean las funciones f y y continuas en algn rectngulo a < t < b, c < y < d
que contiene al punto (t0 , y0 ). Entonces en algn intervalo t0 h < t < t0 + h contenido en
a < t < b existe una nica solucin y = (t) del problema de valor inicial
dy
= r(x, y)
dx
dy
= f (x) g(y)
dx
a veces por comodidad escribimos
dy f (x)
=
dx p(y)
con g(y) = 1/p(y). Para resolver esta ecuacin diferencial se separan las funciones f y g
como
dy
p(y) = f (x)
dx
y se integra ambos lados con respecto a x
Z Z
dy
p(y) dx = f (x)dx.
dx
Dado que y depende de x, de la regla de la cadena tenemos entonces que el diferencial
dy
dy = dx dx, por lo tanto la ecuacin anterior queda en la forma
Z Z
p(y)dy = f (x)dx.
x2
Figura 2.1: Familia de curvas de y = 2 + 2x + c.
En trminos geomtricos esto significa que las soluciones de la ecuacin diferencial son
aquellas funciones y cuya distancia en cada punto x de su dominio a la funcin f (x) 1
x4 x4
es igual a la funcin e 4 ec con c R arbitrario. Estas soluciones son y(x) = 1 e 4 k,
donde por simplicidad en la presentacin hemos definido k = ec > 0. Ahora bien, estas
dos soluciones son parte del conjunto solucin general de las soluciones de la ecuacin
x4
diferencial de la forma: y(x) = 1 + e 4 k, donde se admiten valores negativos para k. Ahora,
usando la condicin inicial tenemos
3 = y(0) = 1 + ke0 = 1 + k.
De este modo tenemos que la solucin del problema de valor inicial es:
x4
y(x) = 1 + 2e 4 .
Ejercicios
1. Encontrar la solucin general de las siguientes ecuaciones diferenciales de variables
separables.
dy
a) dx = xy .
b) (1 y2 )dx xydy = 0.
c) (y + 1)dx + (y 1)(1 + x2 )dy = 0.
d) (2x + 1)y0 + y2 = 0.
e) y0 x3 = x3 y.
f ) xy2 dy
+ (x2 + 1)dx = 0.
g) xy2 + 1 + x2 y0 = 0.
h) xy2 dx + xdy = 0.
dy
i) dx = 1 + x + y + xy.
2.1 Ecuaciones Diferenciales No Lineales 25
dy
j) x2 dx = 1 x2 + y2 x2 y2 .
2. Resuelva el problema de valor inicial respectivo.
dy
a) dx = 4(x2 + 1); y(0) = 1.
dy y2 1
b) dx = x2 1 ; y(2) = 2.
dy
c) x2 dx = y xy; y(1) = 1.
dy 2 2
d) dx = 3x (y + 1); y(0) = 1.
dy
e) 2y dx = 2x ; y(5) = 2.
x 16
dy
f) tan(x) dx = y; y 2 = 2 .
dy 2 2 2
g) dx = 2xy + 3x y ; y(1) = 1.
3. Hallar todas las funciones f (x) tales que su derivada es el cuadrado de la funcin.
es llamada exacta si existe una funcin F(x, y) de dos variables con derivadas parciales
continuas hasta de segundo orden en un dominio , tal que:
F F
= M(x, y), = N(x, y).
x y
M N
= .
y x
F F
= M(x, y), = N(x, y).
x y
2F M 2F N
= , = .
y x y x y x
M N
= . (2.2)
y x
existencia de una funcin F(x, y) tal que dF = Mdx + Ndy. Si esta funcin existe,
debe verificar que
F F
= M(x, y), = N(x, y).
x y
Integrando la primera ecuacin con respecto a x, y manteniendo a y como constante,
se obtiene que
Z
F(x, y) = M(x, y)dx + g(y) (2.3)
Para determinar g(y) es importante que el lado derecho de la ecuacin (2.4) sea una
funcin que depende solamente de y. Para establecer este hecho derivamos (2.4) con
respecto a x, obteniendo as
d 0 N(x, y) M(x, y)
g (y) = .
dx x y
Por la ecuacin (2.2) esta expresin es cero, por lo tanto el lado derecho de (2.4) no
depende de x, integrando (2.4) y sustituyendo g(y) en (2.3), se obtiene la solucin de
la ecuacin (2.1)
Z Z R
( M(x, y)dx)
k = F(x, y) = M(x, y)dx + N(x, y) dy.
y
Este es el mtodo que se emplea para determinar F(x, y) y es indiferente empezar el clculo
de esta manera o bien con: Z
N(x, y)dy + g(x),
p(y)dy f (x)dx = 0
M N
= ex + xex =
y x
se tiene que la ecuacin es exacta. Para determinar F(x, y), comenzamos integrando con
respecto a x
F
=M(x, y) = ex y + xex y
x Z
F(x, y) = (ex y + xex y)dx + g(y)
F(x, y) =xex y + g(y). (2.6)
F
=xex + g0 (y)
y
xex + 2 =xex + g0 (y)
2 =g0 (y).
As, 2 = g0 (y), y g(y) = 2y. La constante de integracin se puede omitir ya que cualquier
solucin de la ecuacin diferencial se satisface. Por tanto, sustituyendo g(y) en (2.6) tenemos
xex y + 2y = k.
2.1 Ecuaciones Diferenciales No Lineales 29
Comparando (tomando los trminos comunes una vez y los no comunes), obtenemos
F(x, y) = xyex + 2y. Y la solucin de la ecuacin est dada de manera explcita por
k =xex y + 2y
k
y= x ,
xe + 2
con la condicin inicial x = 0, tenemos y = 1 por lo que k = 2. De esta manera la
solucin es:
2
y= x .
xe + 2
M N
y
y x
y encontramos que
M N
= sin(xy) + xy cos(xy) =
y x
luego la ecuacin es exacta. Para determinar F(x, y), comenzamos integrando M con respecto
ax
F
=M(x, y) = y sin(xy) + x2
x
30 Ecuaciones Diferenciales de Primer Orden
Z
F(x, y) = (y sin(xy) + x2 )dx + g(y)
x3
F(x, y) = cos(xy) + + g(y).
3
Ahora derivamos parcialmente F con respecto a y y sustituimos x sin(xy) + y2 en lugar de
N(x, y) ya que F/ y = N(x, y):
F
=x sin(xy) + g0 (y)
y
x sin(xy) + y2 =x sin(xy) + g0 (y)
y2 =g0 (y).
As,
Z
y3
g(y) = y2 dy = .
3
Por tanto tenemos que,
x3 y3
F(x, y) = cos(xy) + +
3 3
es decir,
x3 y3
k = cos(xy) +
+ .
3 3
Ahora usemos la condicin inicial y(0) = 1. Tenemos
03 (1)3
k = cos(0(1)) + +
3 3
1 4
k =1 = .
3 3
As la solucin del problema de valor inicial viene dada en forma implcita por
4 x3 y3
= cos(xy) + + .
3 3 3
Ejercicios
1. Determinar si la ecuacin es exacta. Si lo es, resulvala.
a) x2 dy + 2xydx = x2 dx.
b) (y2 2x)dx + 2xydy 2= 0.
1 2
c) x + 2y x dx + (2yx cos(y))dy = 0.
d) (y2 sin(x))dx + 1x xy dy = 0.
e) (tan(y) 2)dx + x sec2 (y) 1y dy = 0.
2.1 Ecuaciones Diferenciales No Lineales 31
2x 3y2 2y x2 1
f ) y x4 dx + x3 y2 + y dy = 0.
f (x)dx + g(y)dy = 0
es exacta.
3. Resuelva el problema de valor inicial respectivo.
a) (4x3 y2 6x2 y 2x 3)dx + (2x4 y 2x3 )dy = 0, y(1) = 3.
b) (4y cos(x)+4 cos(x) sin(x)+sec2 (x))dx+(4y4 sin(x))dy = 0, y 4 = 0.
c) (y3 x)ex dx + 3y2 (ex + y)dy = 0, y(0) = 0.
d) (sin(x) y sin(x))dx + (cos(x) + y)dy = 0, y(0) = 1.
4. Determinar la funcin ms general M(x, y) de modo que la ecuacin sea exacta:
2 x
a) M(x, y)dx + sec (y) y2 dy = 0.
b) M(x, y)dx + (x3 y3 )dy = 0.
c) M(x, y)dx + (e2x ey sin(x))dy = 0.
5. Determinar la funcin ms general N(x, y) de modo que la ecuacin sea exacta:
a) (x2 y3 3xy + 2y2 )dx + N(x, y)dy = 0.
b) (y sin(x) + x sin(y))dx + N(x, y)dy.
6. Encuentre condiciones para las constantes A, B, C y D tales que la ecuacin
sea exacta.
7. Demostrar que la ecuacin y0 + y = 0 es exacta si la multiplicamos por ex y hallar la
solucin general.
Sea,
una ecuacin diferencial que no es exacta. Si multiplicamos esta ecuacin por una funcin
(x, y) adecuada que haga que la ecuacin resultante
sea exacta, se denominar a esta funcin factor de integracin para la ecuacin. Para de-
terminar la funcin (x, y) verificamos el criterio de exactitud en la ecuacin transformada:
(x, y)M(x, y) (x, y)N(x, y)
=
y x
al usar la regla del producto para derivadas esto se reduce a la ecuacin
N M
M N = . (2.8)
y x x y
Por simplicidad, consideraremos que depende de una sola variable. En este caso, supone-
mos que depende de x, es decir = (x), de tal forma que la ecuacin (2.8) se reduce a la
ecuacin de variable separable
" M N #
y x
=
x N
donde
M
y Nx
N
2.1 Ecuaciones Diferenciales No Lineales 33
Ahora, supongamos que la ecuacin (2.7) tiene un factor integrante que solo depende de y,
es decir = (y), en este caso la ecuacin (2.8) se reduce a la ecuacin separable
" N M #
x y
=
y M
donde
N
x My
M
solo depende de y. Entonces en este caso el factor integrante para (2.7) es:
Z N M
!
x y
(y) = exp dy .
M
obtenemos una funcin que solo depende de x, de modo que el factor integrante es dado por
la frmula
Z M N
!
y x
(x) = exp dx ,
N
34 Ecuaciones Diferenciales de Primer Orden
Z
!
1
(x) = exp dx = x.
x
M N
y
y x
y encontramos que
M N
= 4xy + 2x = .
y x
As tenemos en efecto que la ecuacin es exacta. Para determinar F(x, y) comenzamos
integrando M con respecto a x
F
=M(x, y) = 2xy2 + 2xy + 4x3
x Z
F(x, y) = (2xy2 + 2xy + 4x3 )dx + g(y)
F(x, y) =x2 y2 + x2 y + x4 + g(y).
F
=2x2 y + x2 + g0 (y)
y
2x y + x2 =2x2 y + x2 + g0 (y)
2
0 =g0 (y).
As tenemos que
g(y) = c.
Por tanto,
F(x, y) = x2 y2 + x2 y + x4 .
Es decir, la solucin de la ecuacin diferencial viene dada de manera implcita por la relacin
k = x2 y2 + x2 y + x4 .
2.1 Ecuaciones Diferenciales No Lineales 35
Ejercicios
1. Encontrar el factor integrante y resolver las siguientes ecuaciones diferenciales:
a) 2xy2 dx + 3x2 ydy = 0.
b) (x y)dx + xdy = 0.
c) (2y3 + 6xy2 )dx + (3xy2 + 4x2 y)dy = 0.
d) (y2 + 2xy)dx x2 dy = 0.
e) (x2 y + 4xy + 2y)dx + (x2 + x)dy = 0.
f ) (y ln |y| + yex )dx + (x + y cos(y))dy = 0.
2. Demostrar que si ( N/ x M/ y)/(xM yN) solo depende del producto xy, es
decir
( N/ x M/ y)
= H(xy)
(xM yN)
entonces la ecuacin M(x, y)dx + N(x, y)dy = 0 tiene un factor integrante de la forma
(xy). Proporcione la frmula general para (xy).
3. Encontrar el factor integrante (xy) y resolver las siguientes ecuaciones diferenciales:
a) ydx + (x 3x2 y2 )dy = 0.
b) ydx + (x 3x3 y2 )dy = 0.
c) y(x2 y2 + xy)dx + x(x2 y2 1)dy = 0.
4. Suponga que a, b, c y d son constantes tales que cb ad 6= 0 y sean p, q, r y s nmeros
reales arbitrarios. Demostrar que
Sea M(x, y) = (x)P(x)y (x)Q(x) y N(x, y) = (x). Para que esta ecuacin sea exacta
se debe tener que
M N
=
y x
(x)P(x) = 0 (x)
1
d(x) =P(x)dx.
(x)
Al integrar ambos lados tenemos el factor integrante
R
P(x)dx
(x) = e . (2.13)
2.2 Ecuaciones Diferenciales Lineales 37
dy d(x)
(x) + y = (x)Q(x)
dx dx
d (x)y
= (x)Q(x)
dx
R
integrando ambos lados de la ecuacin tenemos (x)y = (x)Q(x)dx + c. Despejando
y(x) obtenemos,
Z
1
y(x) = (x)Q(x)dx + c (2.14)
(x)
Donde R
P(x)dx
(x) = e
es el factor integrante de la ecuacin lineal. Ntese, por qu la necesidad de que P y Q sean
continuas (condicin suficiente para que las integrales existan).
Solucin. Para poder escribir la forma general de la ecuacin diferencial lineal dividimos
por x, ya que x 6= 0, y obtenemos:
dy y 1
+2 = 4
dx x x
dy 2 1
+ y = 4.
dx x x
Donde vemos que
2 1
P(x) = y Q(x) = 4 ,
x x
38 Ecuaciones Diferenciales de Primer Orden
Ejercicios
1. Resolver los siguientes problemas de valor inicial:
a) y0 2y = 3e2x , y(0) = 0.
b) xy0 + 2y = 3x, y(1) = 5.
c) 3xy0 + y = 12x, y(1) = 1.
d) xy0 3y = x3 , y(1) = 10.
e) y0 = (1 y) cos(x), y(0) = 1.
f ) y0 = 1 + x + y + xy, y(0) = 0.
2. Resuelva las ecuaciones diferenciales considerando a y como la variable independiente
en lugar de x.
dy
a) (1 4xy2 ) dx = y3 .
dy
b) (x + yey ) dx = 1.
dy
c) (1 + 2xy) dx = 1 + y2 .
3. Sea y0 una solucin distinta de cero de y0 + p(x)y = 0 y y1 es una solucin de
y0 + p(x)y = q(x).
a) Demostrar por sustitucin que y1 + cy0 es una solucin de y0 + p(x)y = q(x)
para todo nmero c.
b) Demostrar que toda solucin y de y0 + p(x)y = q(x) puede ser escrita como
y = y1 + cy0 para algn c.
4. Sean y1 y y2 soluciones distintas de y0 + p(x)y = q(x), y sea y0 = y1 y2 .
a) Demostrar que y0 es solucin de y0 + p(x)y = 0.
b) Escriba toda solucin de y0 + p(x)y = q(x) en trminos de y1 y y2 .
5. Verificar que las funciones dadas en cada numeral son soluciones de la ecuacin
diferencial y escriba la solucin general por el mtodo del Ejercicio 4.
a) 1 ex , 1 + 2ex ; y0 + y = 1.
b) x 1/x, x + 1/x; y0 + (1/x)y = 2.
c) x, x2 ; (x2 x)y0 + (1 2x)y = x2 .
2.3 Soluciones de Algunas EDO Mediante Sustituciones 39
5. a) y = 1 + cex .
b) y = x + c/x.
c) y = c + c(x2 x).
es homognea.
Solucin. Escribimos la ecuacin en la forma
dy xy2 y3
=
dx 2x2 y
dy xy2 y3
= 2 2
dx 2x y 2x y
dy y y2
= 2
dx 2x 2x
dy 1 y 1 y 2
= .
dx 2x 2 x
Vemos que el lado derecho solo depende de xy , es decir la ecuacin es homognea.
Para resolver una ecuacin homognea hacemos la sustitucin u = xy de tal manera que si
dy
derivamos con respecto a x obtenemos dx = u + x du
dx , entonces sustituimos en la ecuacin
(2.18) para obtener
du
u+x = r(1, u).
dx
y
La nueva ecuacin es de variables separables x y u, una vez resuelta reemplazamos u por x
para dejar la solucin en trminos de las variables originales.
Solucin. Como vimos en el ejemplo anterior esta ecuacin es homognea, ahora hacemos
la sustitucin u = xy y dx
dy
= u + x du
dx .
2.3 Soluciones de Algunas EDO Mediante Sustituciones 41
As:
dy 1 y 1 y 2
=
dx 2 x 2 x
du 1 1
u + x = u u2 .
dx 2 2
La ecuacin anterior es de variable separable. Por tanto obtenemos,
Z Z
2 dx
2
du =
u+u x
2 ln |u| ln |u + 1| = ln |x| + c.
y
Ahora reemplazando u = x tenemos que
y y + x
2 ln + 2 ln = ln |x| + c
x x
es la solucin en forma implcita de la ecuacin.
Ejercicios
1. En las siguientes funciones, determinar si la funcin dada es homognea. Si lo es,
indique su orden de homogeneidad.
4
a) f (x, y) = 3x2 y x3 + 2yx .
2x2 y+xy2
b) f (x, y) = y .
3 2 2
c) f (x, y) = 4xy +2x
xy
y
.
d) f (x, y) = (x + 3y 2)2 .
2. Resolver las siguientes ecuaciones diferenciales.
a) y0 = y+xx .
0 2yx
b) y = y .
y
c) y0 = x+y .
x2 +y2
d) y0 = xy .
e) y0 = x
x+y .
2x2 +3y2 +y
f) y0 = x .
g) (y + xy + x )dx + (2x2 + 3xy + y2 )dy = 0.
2 2
dy y(2x3 y3 )
=
dx x(2y3 x3 )
y 3+ 2x2 +3y2
f)
3
c = 0.
x x
g) 4xy + y2 + x2 = c.
x
h) 2(y+x) ln |y + x| = 21 .
2.3 Soluciones de Algunas EDO Mediante Sustituciones 43
3. a) 21 ln(x2 + y2 ) + arctan y/x = c.
b) 23 ln |y| 13 ln | 3x + y| = c.
c) ln |x+y|3
2
ln |xy|
= c.
2 3
7. a) 3x y + y = c.
b) y = ce2 1x/y .
c) y = 2x arctan(cx).
2
d) xy2 + ln |xy| = 0.
e) y sin xy = 0.
f ) 2ex/y + ln |y| = c.
g) ln(x2 ) ey/x sin xy + cos xy = c.
8. a) y2 = x2 + x.
b) ln |x| + ey/x = 1.
c) x = ex/y1 .
dy a1 x + b1 y a1 + b1 xy
= =
dx a2 x + b2 y a2 + b2 xy
dv a1 (u + h) + b1 (v + k) + c1
=
du a2 (u + h) + b2 (v + k) + c2
dv a1 u + b1 v + (a1 h + b1 k + c1 )
= . (2.20)
du a2 u + b2 v + (a2 h + b2 k + c2 )
44 Ecuaciones Diferenciales de Primer Orden
Como a1 b2 6= a2 b1 , el sistema
(
a1 h + b1 k + c1 = 0
(2.21)
a2 h + b2 k + c2 = 0
tiene solucin. As, para h y k soluciones del sistema (2.21) se tiene que la ecuacin
(2.20) es homognea.
du 3u + 4 du 3u
1 = , o =
dx 2u 7 dx 2u + 7
la cual es una ecuacin de variables separables. Integrando con respecto a x obtenemos
Z Z
2u + 7 du
dx = dx
3 u dx
2u 13 ln |3 u| =x + c.
De esta forma, la solucin general en forma implcita de la ecuacin est dada por
x + c = 2(x + y) 13 ln |3 (x + y)|.
dz 3u zu
z+u =
du u + zu
dz 3 z
z+u =
du 1 + z
dz 3 2z z2
u =
du 1+z
la cual es una ecuacin de variable separable. Integramos con respecto a u para resolverla.
Z Z
z + 1 dz du
2
du =
3 2z z du u
Z Z
z+1 du
2
dz =
z + 2z 3 u
1 2
ln z + 2z 3 = ln |u| + c1
2
2
lnz + 2z 3 = 2 ln |u| + c1
2
lnz + 2z 3 = ln |u2 | + c1
z2 + 2z 3 = u2 ec1 = cu2 .
Ejercicios
1. Hallar la solucin general de las siguientes ecuaciones diferenciales con coeficientes
lineales.
a) (x + 2y 4)dx (2x 4y)dy = 0.
b) (3x + 2y + 1)dx (3x + 2y 1)dy = 0.
c) (x + y 1)dx + (2x + 2y 3)dy = 0.
d) (x + y)dx + (2x + 2y 1)dy = 0.
e) (x + 2y)dx + (3x + 6y + 3)dy = 0.
2. Hallar la solucin del problema de valor inicial.
a) (x + y)dx + (3x + 3y 4)dy = 0; y(1) = 0.
b) (x + 7)dx + (2x + y + 3)dy = 0, y(0) = 1.
c) (x + y + 2)dx (x y 4)dy = 0; y(1) = 0.
sustituyendo u = y1n ,
du
+ (1 n)P(x)u = (1 n)Q(x)
dx
la cual es una ecuacin diferencial lineal.
Ejemplo 2.12 Hallar la solucin general de cada una de las siguientes ecuaciones diferen-
ciales:
dy
1. dx y = ex y2 .
dy 1
2. dx y = ex y 2 .
du
y2 y = ex y2 . (2.25)
dx
du
+ y1 = ex , (2.26)
dx
en (2.26) sustituimos u = y1 ,
du
+ u = ex , (2.27)
dx
como la ecuacin diferencial es lineal, el factor integrante para la ecuacin (2.27) es
R
1dx
(x) = e = ex .
En (2.28) reemplazamos u = y1
ex + 2cex
y1 =
2
2
y= x .
e + 2cex
48 Ecuaciones Diferenciales de Primer Orden
dy
2. Para esta ecuacin de Bernoulli n = 12 , entonces u = y1/2 , por tanto dx = 2y1/2 du
dx .
1/2
Sustituyendo esto y la igualdad u = y en la ecuacin original tenemos
du
2y1/2 y =ex y1/2
dx
du 1 1/2 1 x
y = e
dx 2 2
du 1 1
u = ex .
dx 2 2
Como la ecuacin diferencial es lineal, el factor integrante es
R 1 1
(x) = e 2 dx = e 2 x .
Utilizando la frmula (2.14) tenemos
Z
1
x 1
x 1 x
u(x) =e 2 e 2 e dx + c ,
2
1
1
u(x) =e 2 x e 2 x + c . (2.29)
1
En (2.29) sustituimos u = y 2 . Es decir,
1 1
y 2 =ex ce 2 x .
Ejercicios
1. Hallar la solucin general de las siguientes ecuaciones de Bernoulli
a) xy0 + y = y2 ln |x|.
b) y0 + y = xy3 .
dy
c) (1 x3 ) dx 2(1 + x)y = y5/2 .
d) xy0 y(2y ln |x| 1) = 0.
e) x2 (x 1)y0 y2 x(x 2)y = 0.
2. Hallar la solucin de los siguientes problemas de valor inicial
2
a) y0 + 1x y = yx ; y(1) = 1.
b) 2 cos(x)dy = (y sin(x) y3 )dx; y(0) = 1.
Ejemplo 2.13 Hallar una familia mono paramtrica de soluciones de la ecuacin dife-
rencial
dy 4 1
= 2 y + y2
dx x x
2
en donde y1 = x es una solucin conocida de la ecuacin.
du ds
= u2 .
dx dx
Reemplazando en la ecuacin (2.33) nos queda
2 ds 3
u u = u2 , (2.34)
dx x
En (2.35) sustituimos s = u1 :
ds 3
+ s = 1. (2.36)
dx x
Como la ecuacin diferencial es lineal, el factor integrante para la ecuacin (2.36) es
R 3
(x) = e x dx = e3 ln |x| = x3 .
x4 + c
u1 =
4x3
2.3 Soluciones de Algunas EDO Mediante Sustituciones 51
4x3
u= . (2.38)
x4 + c
Sustituyendo u = y 2x en (2.38) y despejando la funcin y, tenemos que la solucin de la
ecuacin diferencial viene dada por
4x3 2
y=u= 4
+ .
x + c x
Ejercicios
1. Sea la ecuacin de Riccati y0 = f0 (x) + f1 (x)y + f2 (x)y2 , con f2 (x) 6= 0. Si y1 es una
solucin particular de la ecuacin, Demostrar que la sustitucin y = y1 + 1u trasforma
esta ecuacin en una ecuacin lineal de la forma u0 + ( f1 (x) + 2 f2 (x)y1 )u = f2 (x).
2. Dada la solucin particular y1 hallar la solucin general de las ecuaciones de Ricatti
dadas.
a) y0 = x3 + 2x y 1x y2 ; y1 (x) = x2 .
b) y0 = 2 tan(x) sec(x) y2 sin(x); y1 (x) = sec(x).
c) y0 = x12 xy y2 ; y1 (x) = 1x .
2
d) y0 = 1 + xy xy2 ; y1 (x) = x.
y = g(y0 )x + f (y0 ),
dx g0 (p)x + f 0 (p)
=
dp p g(p)
dx 0
g (p) f 0 (p)
x=
d p p g(p) p g(p)
la cual es una ecuacin lineal con solucin
Z
1
x(p) = (p)Q(p)d p + c = h(p, c)
(p)
y = x(y0 + 3) 2(y0 )2 .
y = x(p + 3) 2(p)2
dp dp
y0 = p + 3 + x 4p
dx dx
dp dp
p = p + 3 + x 4p
dx dx
dp
3 = (x 4p),
dx
2.3 Soluciones de Algunas EDO Mediante Sustituciones 53
Ejercicios
1. Hallar las ecuaciones paramtricas de la soluciones de las ecuaciones diferenciales
siguientes
a) y = x(y0 )2 + ln |y0 |.
b) y = x(y 0 2 0 3
)0 +1(y 0) 2. 2 0 3/2
c) y = x y + 2 (y ) 3 (y ) .
d) y = (cos(y0 ) + y0 )x + 12 sin(y0 ) cos(y0 ) + 21 y0 .
54 Ecuaciones Diferenciales de Primer Orden
Ejercicios
1. Hallar las soluciones de las ecuaciones diferenciales siguientes
a) y0 = (x y)2 .
b) y0 = sin(x y + 1).
dy
c) dx = (x + y 1)2 .
dy
d) dx = tan2(x + y).
e) y0 = 2 + y 2x + 3.
dy
f ) dx = 1xy
x+y .
dy
g) dx = 1 + eyx+5 .
cualquiera con una rapidez proporcional al nmero de personas presentes en dicho instante.
Si la poblacin se duplica en 5 aos. Cunto demorar en triplicarse?
Solucin. De acuerdo con el enunciado del problema, un modelo para estudiar este fen-
meno es el siguiente:
dP
= kP, (3.2)
dt
Otra aplicacin modelada mediante ecuaciones autnomas son los llamados problemas
financieros de inters: Consideremos una cantidad A a ser invertida al % anual. El capital
P al final del ao ser
P =A 1 + si el inters es compuesto es anualmente,
100
2
P =A 1 + si el inters es compuesto es semestralmente,
2 100
4
P =A 1 + si el inters es compuesto es trimestralmente,
4 100
12
P =A 1 + si el inters es compuesto es mensualmente.
12 100
En general, el capital P al final del ao ser
m
P = A 1+ .
m 100
Si la tasa de inters es r por ciento anual compuesto m veces por ao, al final de n aos el
capital ser:
h r m in
P = A 1+ .
m
Si el nmero m incrementa sin lmite entonces
h
r m in r m/r nr
P = lm A 1 + = lm A 1 + = Aenr .
m m m m
Finalmente, reemplazando n por t obtenemos que
P = Aert
Ejercicios
1. Suponga que la tasa de crecimiento por da de una poblacin es 0.7944 por miembro.
La poblacin inicial es de 2 miembros, hallar el tamao de la poblacin al final de 5
das.
2. En 1990 el Departamento de Recursos Naturales liber 2.000 ejemplares de una
especie de pez en un lago. En 1997, la poblacin de estos peces en el lago se estim
en 5.000. Use la ley de Malthus para el crecimiento de poblaciones y estime la
poblacin de estos peces en el lago en el ao 2010.
3. Suponga que una poblacin tiene 50 miembros en t = 0 y 75 al final de 100 das.
Hallar la poblacin al final de 150 das.
4. En cualquier tiempo t la cantidad de bacterias en un cultivo crece a razn proporcional
al nmero de bacterias presentes. Al cabo de dos horas se observa que hay 600
individuos. Despus de 15 horas hay 3.000 especmenes. Cul era la cantidad inicial
de bacterias?
5. En 1990 la poblacin de lagartos en los terrenos del Centro Espacial Kennedy se
estim en 350. En 2000, la poblacin haba aumentado hasta un estimado de 1.300.
Use la ley de Malthus para el crecimiento de poblaciones y estime la poblacin de
lagartos en dichos terrenos en el ao 2010.
6. Suponga que la tasa de crecimiento diaria de una poblacin dada es 0.21 por miembro.
Si el tamao de la poblacin en un da particular es 100, hallar el tamao de la
poblacin 7 das despus.
7. La poblacin de una comunidad crece a razn proporcional a la poblacin en cualquier
momento t. Su poblacin inicial es de 6.500 y aumenta 15 % en 10 aos. Cul ser
la poblacin en 30 aos?
8. Suponga que una poblacin tiene 39 miembros en t = 8 y 60 miembros en t = 12.
Cual fue el tamao de la poblacin en t = 0?
9. A qu tasa de inters se duplicarn $300 en 8 aos?
10. Suponga que una poblacin tiene 15 miembros en t = 4 y 10 miembros en t = 12.
Cual fue el tamao de la poblacin en t = 0?
11. Cual ser el rdito de $1.000 al 4,5 % de inters luego de 20 aos?
12. En 1986 la poblacin de la tierra era de 5 109 . Use un modelo de crecimiento
exponencial con la tasa de crecimiento de la poblacin del 2 % por ao observado en
1986 para predecir el tamao de la poblacin de la tierra en el ao 2000.
13. Cuanto dinero se debe depositar en un banco al 5 % de inters para poder retirar
$3.600 anuales y que el capital sea consumido al final de este tiempo si el dinero se
62 Algunas Aplicaciones de EDO de Primer Orden
retira continuamente desde la fecha del deposito una cantidad de 3.600/365 al da?
14. Una poblacin de bacterias es inoculada en un placa de Petri a una densidad de 10/ml.
La densidad de las bacteria se duplica en 20 horas. Suponga que esta situacin es
descrita mediante la ecuacin diferencial
dx
= Cx,
dt
donde x es la densidad de las bacterias y C es constante. Cuanto tiempo se requiere
para que la densidad se incremente en 8 veces su valor original?
15. Suponga que una poblacin tiene una tasa de crecimiento constante de r por miembro,
por unidad de tiempo y que el tamao inicial de la poblacin es x0 . Demostrar que el
tamao de la poblacin en un tiempo t es x0 er(tt0 ) .
16. Suponga que una poblacin tiene una tasa de crecimiento de r por miembro, por cada
unidad de tiempo, con r > 0. Demostrar que el tiempo requerido para que la poblacin
duplique su tamao inicial es ln(2)/r.
17. Suponga que una poblacin tiene una tasa de crecimiento de r por miembro, por
unidad de tiempo con r < 0. Demostrar que el tiempo requerido para que la poblacin
decrezca a la mitad de su tamao inicial es ln(2)/r.
18. Suponga que una poblacin tiene una tasa de crecimiento r(t) que depende de la
poblacin en cada instante de tiempo t, as que el tamao x(t) de la poblacin es
gobernado por el problema de valor inicial
dx
= r(t)x, x(0) = x0 .
dt
Demostrar que el tamao deR la poblacin en cualquier instante de tiempo t viene dado
t
por la expresin x(t) = x0 e 0 r(s)ds .
o equivalentemente,
NP
= ekt+C1 = ekt eC1 = C2 ekt .
N P
Debemos considerar dos situaciones:
Caso I: Si N P > 0, entonces
NP
= C2 ekt
N P
NP = NC2 ekt PC2 ekt
P N +C2 ekt = NC2 ekt
NC2 ekt N
P= = N kt
N +C2 ekt C2 e +1
N
P=
NC3 ekt + 1
donde C3 = C12 .
Caso II: Si N P < 0, entonces
NP
= C2 ekt
N P
NP = NC2 ekt + PC2 ekt
P N C2 ekt = NC2 ekt
NC2 ekt N
P= kt
= N kt
N C2 e C e 1
2
N
P=
1 NC3 ekt
donde, otra vez, hemos hecho C3 = C12 .
dP
Si la poblacin inicial es P(0) = P0 < N entonces la poblacin crecer (pues dt es positivo
suponiendo k > 0) hasta el soporte:
N
lm P = lm = N.
t t 1 NC3 ekt
Ejemplo 3.3 En cierta ciudad hay 4.800 personas susceptibles de contraer una cierta enfer-
medad contagiosa. La tasa de crecimiento de una epidemia es conjuntamente proporcional
al nmero de personas infectadas y al nmero de personas que todava no se han infectado.
Inicialmente contraen la enfermedad 300 personas y 10 das despus se han infectado 1.200.
1. Cuntas personas se habrn infectado al cabo de 20 das?
2. Cundo se propagar con mayor rapidez la enfermedad? (es decir, al cabo de cuntos
das se habr infectado la mitad de las personas susceptibles de contraer la enferme-
dad).
Solucin. En este caso N = 4,800 y la ecuacin est dada por
dP
=k1 P(4.800 P)
dt
P
=4.800k1 P 1 .
4.800
Ntese que se trata de una ecuacin logstica con k =4.800k1 y soporte N =4.800, adems
P(0) = 300 <4.800= N de tal forma que la solucin est dada por
N
P(t) =
NC3 ekt + 1
4.800
= .
4.800C3 ekt + 1
4.800
P(t) = .
15ekt + 1
Se sabe tambin que P(10) =1.200, lo que permite encontrar el valor del parmetro k.
4.800
P(10) = = 1.200
15ek10 + 1
4
=1
15ek10 + 1
1
ek10 =
5
1 1
k = ln 0,16094.
10 5
4.800
P(t) = .
15e0,16094t +1
4.800
P(20) = 3.000.
15e0,16094(20) + 1
P(t) =2.400
4.800 4. 800
0,16094t
= = 2. 400
15e +1 2
2
=1
15e0,16094t + 1
15e0,16094t + 1 =2
1
e0,16094t =
15
1 1
t = ln 16,8 das.
0,16094 15
3.1 Modelos de Crecimiento y Decrecimiento 67
Ejercicios
1. Las observaciones sobre el crecimiento de tumores en animales indican que el tamao
y(t) del tumor en el tiempo t puede describirse mediante la ecuacin diferencial
dy y
= ky ln ,
dt a
N
P=
1 + ekt
10. Suponga que una poblacin satisface una ecuacin diferencial de tipo logstica pero
con tasa de crecimiento variable que depende de t:
dP P
= k(t)P 1 , P(0) = P0 .
dt N
Muestre que la solucin viene dada por
N p0
P(t) = Rt .
P0 + (N P0 )e 0 k(s)ds
11. El rea de la biologa que trata con las tasas de crecimiento comparativas de diferentes
cantidades biolgicas se llama alometra. Sea x(t) la tasa de crecimiento de una
cantidad biolgica A y y(t) la tasa de crecimiento de una cantidad biolgica B.
Suponga que ambas tasas de crecimiento son proporcionales
1 dy 1 dx
= .
y dt x dt
Demostrar que la siguiente ley de alometra se satisface y = Kx , donde K es cons-
tante.
12. La tasa metablica (TM) se define como la tasa de produccin de calor por unidad de
tiempo. La energa de los alimentos requerida en reposo y ayuno (la tasa metablica
basal) a travs de una ventana de tiempo T est dada por
TM T.
Ejemplo 3.5 En una cueva de Sudfrica se hall un crneo humano junto con los restos de
una hoguera. Los arquelogos creen que la edad del crneo es igual a la edad de la hoguera.
Se ha determinado que solo queda 2 % de la cantidad original de carbono 14 en los restos
de madera en la hoguera. Estime la edad del crneo.
70 Algunas Aplicaciones de EDO de Primer Orden
Solucin. Nuevamente, se tiene que A(t) = A0 ekt donde A(0) = A0 cantidad inicial. Cuando
t =5.600 aos, A(5.600) = A0 /2, de lo cual es posible determinar el valor de k como sigue:
A0
=A0 e5.600k
2
1
5.600k = ln = ln(2)
2
ln(2)
k = = 0,00012378.
5.600
Por lo tanto, A(t) = A0 e0,00012378t . Si queda el 2 % de la cantidad inicial entonces
2
A0 = A0 e0,00012378t
100
despejando t resulta
2
ln = 0,00012378t
100
ln 2/100
t = 31.604,6 31.605 aos.
0,00012378
De esta forma determinamos que la edad del crneo es de aproximadamente 31.605 aos.
Ejercicios
1. Se sabe que un material radioactivo se desintegra con una rapidez proporcional a
la cantidad presente en cualquier instante. Si inicialmente hay 150 mg de material
y, despus de tres aos, se observa que el 10 % de la masa original se desintegr,
determinar:
a) Una expresin para la masa al momento t.
b) El tiempo necesario para que se desintegre el 20 % de la masa original.
2. Un material radioactivo se desintegra tres quintos en 2.000 aos. Determinar su vida
media.
3. En los problemas a) y b) suponga que la razn de decaimiento de una sustancia
radiactiva es proporcional a la cantidad de sustancia presente. La vida media de una
sustancia radiactiva es el tiempo que tarda en desintegrarse la mitad de la sustancia.
a) Si en un principio se tiene 80 g de una sustancia radiactiva y despus de 2 das
solo restan 20 g, qu porcentaje de la cantidad original quedar despus de 5
das?
b) Si en un principio se tienen 400 g de una sustancia radiactiva y despus de 10
aos restan 100 g. Cunto tiempo deber transcurrir para que solo queden 5 g?
4. El 209 Pb, isotopo radioactivo del plomo, se desintegra con una rapidez proporcional
a la cantidad presente en cualquier tiempo t y tiene una vida media de 3,3 horas. Si
al principio haba 0,5 gramos de plomo, cunto tiempo debe transcurrir para que se
desintegre 95 %?
3.2 Problemas de Mezclas 71
dA
= (Rapidez de entrada) (Rapidez de salida)
dt
dA
= R1 R2 (3.4)
dt
donde
Rapidez de entrada = (flujo de entrada) (concentracin).
Flujo de salida
3
= 25 + ce 125 t .
A(0) = 25 + c = 0
c = 25.
Luego,
3
A(t) = 25 25e 125 t .
La cual es la cantidad de medicamento en el rgano en cada instante t. Para determinar la
concentracin, se divide la cantidad entre el volumen. Es decir,
3
A(t) 25 25e 125 t
C(t) = =
125 125
3
125 t
1 1e
C(t) = . (3.6)
5
La ecuacin (3.6) es la concentracin de medicamento en el rgano en cualquier instante t.
Para determinar en que momento t0 la concentracin es 0, 1 g/cm3 reemplazamos en (3.6)
3
1 e 125 t0
0,1 =
5
125 ln |0,5|
t0 = 28,88 29.
3
Ejemplo 3.7 Una solucin de cido ntrico entra a una razn constante de 6 litros/minutos
en un tanque de gran tamao que en un principio contena 200 litros de una solucin de
cido ntrico al 0,5 %. La solucin dentro del tanque se mantiene bien revuelta y sale del
tanque a razn de 8 litros/minutos. Si la solucin que entra en el tanque tiene cido ntrico
al 20 %, determinar el volumen de cido ntrico en el tanque despus de t minutos.
Solucin. Sea V (t) el volumen de cido ntrico en el tanque en cada instante de tiempo t.
En problemas de esta clase la rapidez neta con que V (t) cambia est dada por
dV
= R1 R2 (3.7)
dt
donde R1 es la rapidez de entrada y R2 la rapidez de salida.
Primero debemos determinar la rapidez con que el cido ntrico entra al tanque, sabemos
que la solucin entra a una razn constante de 6 litros/minuto. La solucin que entra en el
74 Algunas Aplicaciones de EDO de Primer Orden
20
tanque tiene cido ntrico al 20 %, la concentracin es 100 , concluimos que la rapidez de
entrada en el tanque es:
litros 20 6 litros
R1 = 6 = .
minutos 100 5 minutos
Ahora debemos determinar la rapidez de salida. La diferencia entre la razn de flujo de
entrada y la razn de flujo de salida es 6 8 = 2 litros/minuto, de modo que el volumen
de fluido en el tanque despus de t minutos es (200 2t) litros. Por lo tanto la rapidez con
que el cido ntrico sale del tanque es
litros V (t) 4V (t) litros
R2 = 8 = .
minutos 200 2t 100 t minutos
En un principio el tanque contena 200 litros de una solucin de cido ntrico al 0,5 %,
0,5
de modo que V (0) = 200 litros, es decir 100 200 = 1 litros, al sustituir R1 y R2 en (3.7)
tenemos:
dV 6 4V
= , V (0) = 1litro
dt 5 100 t
dV 4V 6
+ = . (3.8)
dt 100 t 5
4
La ecuacin (3.8) es lineal, con P(t) = 100t y Q(t) = 56 , entonces podemos hallar el factor
integrante sustituyendo en la frmula (2.13),
R 4
(t) = e 100t dt = e4 ln(100t) = (100 t)4 .
Luego,
V (t) = 0,4(100 t) (3,9 104 )(100 t)4 .
El cual es el volumen de cido ntrico en el tanque despus de t minutos.
3.2 Problemas de Mezclas 75
Ejercicios
1. Un tanque de 400 galones contiene inicialmente 200 galones de agua que contiene 2
partes por 1.000 millones en peso de una dioxina que es una carcingena extrema-
damente potente. Suponga que el agua que contiene 5 partes por 1.000 millones de
dioxina fluye hacia arriba del tanque a razn de 4 galones por minuto. El agua en el
tanque se mantiene bien mezclada y se retiran dos galones por minuto por el fondo
del tanque. Cunta dioxina se encuentra en el tanque cuando est lleno?
2. Una cubeta de 5 galones est llena de agua pura. Suponga que empezamos a aadir
sal a la cubeta a razn de 14 de libra por minuto. Adems, abrimos el grifo de manera
que salga 12 galn por minuto de la cubeta y agregamos agua pura para mantener
llena la cubeta. Si la solucin de agua salada est siempre bien mezclada. Cul es la
cantidad de sal en la cubeta despus de:
a) 1 minuto?
b) 10 minutos?
c) 60 minutos?
3. A un tanque parcialmente lleno con 100 galones de salmuera, con 10 lb de sal disuelta,
le entra salmuera con 1/2 lb de sal por galn a razn de 6 gal/min. El contenido del
tanque est bien mezclado y de l sale a razn de 4 gal/min de solucin. Calcular la
cantidad de libras de sal que hay en el tanque a los 30 minutos.
4. Un depsito contiene 50 litros de una solucin compuesta por 90 % de agua y 10 %
de alcohol. Se vierte en el depsito, a razn de 4 litros/minuto, una segunda solucin
que contiene 50 % de agua y 50 % de alcohol. Al mismo tiempo se vaca el depsito a
razn de 5 litros/minuto. Suponiendo que la solucin se agita constantemente, calcular
la cantidad de alcohol que queda despus de 10 minutos.
5. Una piscina cuyo volumen es de 10.000 galones contiene agua con cloro al 0,01 %.
A partir del instante t = 0, se bombea agua del servicio publico con cloro al 0,001 %
hacia la piscina, a razn de 5 galones/minuto. El agua sale de la piscina con la misma
razn. Cul es el porcentaje de cloro en la piscina despus de 1 hora? En qu
momento el agua de la piscina tendr 0,002 % de cloro?
6. La sangre conduce un medicamento a un rgano a razn de 3 cm3 /seg y sale con la
misma razn. El rgano tiene un volumen lquido de 125 cm3 . Si la concentracin
del medicamento en la sangre que entra al rgano es de 0,2 g/cm3 . Cul es la
concentracin del medicamento en el rgano en el instante t, si inicialmente no
haba rastros de dicho medicamento? En qu momento llegar la concentracin del
medicamento en el rgano a 0,1 g/cm3 ?
4. 26,55 lts.
5. 210,72 min.
6. 28,81 seg.
Ejemplo 3.8 Un termmetro se saca de una habitacin, en donde la temperatura del aire
es de 700 F, al exterior, en donde la temperatura es de 100 F. Despus de 1/2 minuto el
termmetro marca 500 F. Cunto marca el termmetro cuando t = 1 minuto? Cunto
demorar el termmetro en alcanzar los 150 F?
Solucin. Es claro que Tm = 100 F, reemplazando en la ecuacin (3.9) se tiene:
dT
= k(T 10) (3.10)
dt
T (0) = 700 F
1
T = 500 F.
2
Se debe resolver la ecuacin (3.10) con condiciones iniciales, la cual es lineal y tambien
de variables separables. Empleando la tcnica de resolucin de ecuaciones de variables
separables obtenemos:
dT
=kdt
T 10
Z Z
dT
= kdt
T 10
ln | T 10 |=kt + c1
T 10 =c2 ekt , c2 R.
Entonces T = 10+c2 ekt . Cuando t = 0, se tiene que T = 700 F, de modo que 70 = 10+c2 es
decir c2 = 60 y, por lo tanto, T = 10 + 60ekt . Ahora podemos hallar k ya que T (1/2) = 500 F
1
50 =10 + 60ek 2
1 40
e2k =
60
3.3 Problemas de Temperatura: La Ley de Enfriamiento de Newton 77
2
k =2 ln 0,81093.
3
En consecuencia,
T (t) = 10 + 60e0,81093t .
Para t = 1 tenemos T (1) = 10 + 60e0,810931 = 36,66 370 F, el termmetro marca 370 F
cuando t = 1 minuto. Ahora debemos hallar cunto se demora para alcanzar 150 F,
15 =10 + 60e0,81093t
5
=e0,81093t
60
ln 1/12
t= = 3,06427 3,1minutos.
0,81093
As, para alcanzar 150 F se demora 3,1 minutos.
1
y K es la constante de tiempo para el edificio (sin calefaccin o aire acondicionado).
Ejemplo 3.9 En una calurosa maana de sbado, cuando las personas trabajan dentro del
edificio, el aire acondicionado mantiene la temperatura interior en 240 C. A medioda, el aire
acondicionado se apaga y las personas se van a casa. La temperatura exterior es constante e
igual a 350 C durante el resto de la tarde. Si la constante de tiempo del edificio es de 4 horas,
cul ser la temperatura dentro del edificio a las 2 : 00 p.m.? Y a las 6 : 00 p.m.? En qu
momento llegar la temperatura interior del edificio a 270 C?
1
Solucin. Se tiene M(t) = 350 C, T (0) = 240 C, H(t) = 0, U(t) = 0, K = 4, por lo tanto
tenemos que K = 14 . Hallemos T (2). Reemplazando en la ecuacin
Z
Kt Kt
T (t) = e e [KM(t) + H(t) +U(t)]dt + c
obtenemos,
Z
14 t 1 1
T (t) = e e 35dt + c4t
4
1
h 1 i
4t t
T (t) = e e 35 + c
4
1
T (t) = 35 + ce 4 t .
Por lo tanto,
K 2
B w
c = 16 B0 + .
K 2
1+ w
Es decir,
2 !
B cos(wt) + Kw sin(wt) B Kw Kt
T (t) = B0 + 16 B0 + 2 e .
w 2 K
1+ t 1+ w
Como H0 = 0 y B0 = M0 cuando no hay razn de calentamiento adicional dentro del edificio
H0 = 0 la temperatura promedio interior B0 es igual a la temperatura promedio exterior
M0 . Por lo tanto, B0 = 240 C. Como el trmino exponencial se ha extinguido, entonces c = 0.
Luego, reemplazando tenemos
B Kw sin(wt) + Kw cos(wt)
T (t) = 24 2 .
1 + Kw
M(t) = M0 + B cos(wt)
M(t) = 24 B cos t .
12
Para t = 0 y M(t) = 16 se tiene B = 8. Por lo tanto,
8 Kw sin(wt) + Kw cos(wt)
T (t) = 24 2 .
1 + Kw
8(4)2
T (0) = 24 2
16,50 C
1 + (4)
1
para K = 5, se tiene K = 15 . La temperatura mnima es
4 2
8
T (0) = 24 5
2 = 20,80 C.
1 + 45
1
Para hallar la temperatura mxima t = 12, K = 1, se tiene K = 1
8(4)2
T (12) = 24 + 2
31,50 C
1 + (4)
3.3 Problemas de Temperatura: La Ley de Enfriamiento de Newton 81
1 1
para K = 5, se tiene K = 5 la temperatura mxima es
4 2
8
T (12) = 24 + 5
2 = 27,10 C.
1 + 45
Ejemplo 3.11 Durante el verano la temperatura dentro de una camioneta llega a 550 C,
mientras que en el exterior es constante e igual a 350 C. Cuando la conductora entra a la
camioneta, enciende el aire acondicionado con el termostato en 160 C. Si la constante de
tiempo para la camioneta es K1 = 2 horas y para la camioneta con el aire acondicionado es
1 1 0
K1 = 3 hora. En qu momento llegar la temperatura dentro de la camioneta a los 27 C?
K1 = K + Ku
1
3 = + Ku
2
5
Ku = .
2
Reemplazando en (3.13) obtenemos
dT 1 5
+ 3T = 35 + 16
dt 2 2
dT 115
+ 3T =
dt 2
82 Algunas Aplicaciones de EDO de Primer Orden
la cual es una ecuacin diferencial lineal. Su factor integrante es (t) = e3t , por lo tanto la
solucin es
Z
3t 3t 115
T =e e dt + c
2
3t 3t 115
T =e e dt + c
6
115
T= + ce3t .
6
Como T (0) = 55, tenemos
115
55 = +c
6
215
c=
6
115 215 3t
T (t) = + e ,
6 6
para hallar en qu momento t0 la temperatura llegar a 270 C, reemplazamos en la ltima
ecuacin y despejamos t0 .
115 215 3t0
27 = + e
6 6
t0 = 0,5068 horas
t0 = 30,4 minutos.
Ejercicios
1. Un termmetro se lleva del interior de una habitacin al exterior, donde la temperatura
del aire es 600 F. Despus de un minuto, el termmetro indica 580 F; cinco minutos
despus marca 320 F. Cul era la temperatura interior?
2. Si una barra metlica pequea, cuya temperatura inicial es de 250 F, se deja caer
en un recipiente con agua hirviente cuya temperatura es de 2120 F. Cunto tiempo
tardar en alcanzar 1000 F si se sabe que la temperatura aument 30 F en un segundo?
Cunto tiempo tardar en llegar a 1050 F?
3. Una taza de caf caliente, inicialmente a 1000C, se enfra hasta 850C en 5 minutos
al estar en un cuarto con temperatura de 210C. Use solo la ley de enfriamiento de
Newton y determinar el momento en que la temperatura del caf estar a unos 550C.
4. Una cerveza fra, inicialmente a 330 F, se calienta hasta 450 F en 3 minutos estando
en un cuarto con temperatura 650 F. Qu tan caliente estar la cerveza si se deja ah
durante 10 minutos?
5. Una taza de t caliente est inicialmente a 1900 F y se deja en un cuarto que tiene
una temperatura ambiente de 750 F. Suponga que a partir del tiempo t = 0 se enfra a
razn de 250 F por minuto.
3.3 Problemas de Temperatura: La Ley de Enfriamiento de Newton 83
1. 71, 5160 F.
2. a) Para llegar a 1000 F tarda 31,7 seg.
b) Para llegar a 1050 F tarda 34,526 seg.
3. 20,02 min.
4. 58,30 F.
5. a) T = 75 + 115e0,24512t .
b) 3, 82 min.
6. a) T = 72 + 108e0,26329t .
b) 5, 13 min.
c) Muri a las 4 a.m. con 21 minutos.
84 Algunas Aplicaciones de EDO de Primer Orden
p(t) = mv(t)
b
F2 = bv(t)
t=0
x(t)
b
F1 = mg
b t
Solucin. Hay dos fuerzas actuando sobre el objeto: una fuerza constante debida al empuje
hacia abajo de la gravedad y una fuerza debida a la resistencia del aire que es proporcional
a la velocidad del objeto y actun en forma opuesta al movimiento del objeto. Por lo tanto,
el movimiento del objeto se realizar a lo largo de un eje vertical. En este eje, elegimos el
origen como el punto donde el objeto fue lanzado inicialmente y definimos x(t) como la
distancia que ha caido el objeto hasta el instante t.
Las fuerzas que actan sobre el objeto a lo largo de este eje se pueden expresar de la
siguiente manera:
La fuerza debida a la gravedad
F1 = mg
donde g es la aceleracin de la gravedad en la superficie de la tierra.
La fuerza producto de la resistencia del aire
F2 = bv(t)
F = F1 + F2 = mg bv(t). (3.15)
dx
v(t) = = 4,905 + 45,095e2t
dt
donde, integrando respecto a t concluimos
Puesto que el objeto se liber a 500 m sobre el suelo, podemos determinar el momento en
que el objeto toca el suelo haciendo x(t) = 500 y despejando t. As escribimos
Solucin. Solo nos interesa el momento en que el paracaidista toca el suelo, no el lugar. As,
nicamente consideramos la componente vertical de su descenso. Para esto, necesitamos
usar dos ecuaciones, una para describir el movimiento antes de abrir el paracadas y la otra
para aplicar despus de abrirlo.
m
Antes de abrirse: m = 75 kg, v0 = v(0) = 0, b = b1 = 30 N.seg/m, g = 9,81 seg 2 y
Ejercicios
1. Un objeto de 400 libras se libera desde el reposo a 500 pies sobre el suelo y se le
permite caer bajo la influencia de la gravedad. Suponiendo que la fuerza en libras
debida a la resistencia del aire es 10v, donde v es la velocidad del objeto en pies/seg,
determinar la ecuacin de movimiento del objeto. En qu momento tocar el objeto
el suelo?
2. Un objeto de masa 8 kg recibe una velocidad inicial hacia abajo de 20 m/seg y luego
se le permite caer bajo la influencia de la gravedad. Suponga que la fuerza en Newtons
debida a la resistencia del aire es 16v, donde v es la velocidad del objeto en m/seg.
Determinar la ecuacin de movimiento del objeto. Si el objeto est inicialmente a 100
m sobre el suelo, determinar el momento en que el objeto golpear el suelo.
3. Un objeto con masa de 2 kg se lanza desde el reposo de una plataforma a 30 m sobre el
agua y se deja caer bajo la influencia de la gravedad. Despus de que el objeto golpea
el agua, comienza a hundirse, con la gravedad jalndolo hacia abajo y una fuerza de
flotacin empujndolo hacia arriba. Suponga que la fuerza de gravedad es constante,
que no hay cambios en el momento del objeto al golpear el agua, que la fuerza de
flotacin es 1/2 del peso (peso = mg), y que la fuerza debida a la resistencia del aire o
del agua es proporcional a la velocidad del objeto, con constante de proporcionalidad
b1 = 10N.seg/m en el aire y b2 = 30N.seg/m en el agua. Determinar la ecuacin de
movimiento del objeto. Cul es la velocidad del objeto al minuto despus de ser
arrojado?
4. Un paracaidista cuya masa es de 100 kg se arroja de un helicptero que vuela a
3.000 m sobre el suelo y cae bajo la influencia de la gravedad. Suponga que la fuerza
debida a la resistencia del aire es proporcional a la velocidad del paracaidista, con
la constante de proporcionalidad b3 = 30 N.seg/m cuando el paracadas est cerrado
y b4 = 50 N.seg/m cuando se abre. Si el paracadas no se abre hasta 30 segundos
despus de que el paracaidista sale del helicptero, despus de cuntos segundos
llegar l al suelo? Si el paracadas no se abre hasta 1 minuto despus de que el
paracaidista sale del helicptero, despus de cuntos segundos llegar l al suelo?
5. Un objeto con masa de 100 kg se lanza desde el reposo de una lancha hacia el agua y se
deja hundir. Aunque la gravedad jala el objeto hacia abajo, una fuerza de flotacin de
1/40 veces el peso del objeto lo empuja hacia arriba (peso = mg). Si suponemos que
la resistencia del agua ejerce sobre el objeto una fuerza proporcional a la velocidad
del objeto, con constante de proporcionalidad 10 N.seg/m. Determinar la ecuacin de
movimiento del objeto. Despus de cuantos segundos ocurrir que la velocidad del
objeto es igual a 70 m/seg?
6. Un objeto de masa m se libera desde el reposo y cae bajo la influencia de la gravedad.
Si la magnitud de la fuerza debida a la resistencia del aire es bvn , donde b y n son
constantes positivas, determinar la velocidad lmite del objeto (suponiendo que este
lmite existe). (Sugerencia: Justifique que la existencia de una velocidad lmite (finita)
implica dvdt 0 cuando t +).
7. Cuando la velocidad v de un objeto es muy grande, la magnitud de la fuerza debida a
la resistencia del aire es proporcional a v2 y la fuerza acta en direccin opuesta al
3.5 Trayectorias Ortogonales 89
movimiento del objeto. Un proyectil con masa de 3 kg se lanza desde el suelo hacia
arriba y con una velocidad inicial de 500 m/s. Si la magnitud de la fuerza debida a
la resistencia del aire es 0,1v2 , en qu momento alcanzar el proyectil su mxima
altura sobre el suelo?. Cul es esa mxima altura?
8. Si un objeto cae desde un aeroplano, la velocidad hacia abajo despus de x seg es
aproximadamente y0 = gk (1 ekx ) en donde g = 9,8 segm
2 y k = 0, 2seg
1 . y(x) es la
distancia que ha cado, por lo tanto y(0) = 0. Si este objeto cae desde una altura de
5.000 m. Cunto le tomar llegar a la tierra?
1. 13,74 seg.
2. 488,45 seg.
3. La ecuacin del movimiento es y(t) = 1,96t + 0,392e2t 0,392 y la velocidad
despus de un minuto es 0,327 m/seg.
4. El paracaidista llega al suelo, si el paracadas abre luego de 30 seg., en 424,28 seg.
5. Ecuacin de movimiento x(t) = 96t + 960e0,1t 960. en t = 13, 06 el cuerpo alcanza
la
qvelocidad de 70 m/seg.
n mg
6. b .
Dos curvas se cortan en ngulo recto si las respectivas rectas tangentes, en los puntos de
interseccin, son rectas perpendiculares. En esta seccin utilizaremos ecuaciones diferencia-
les de primer orden para establecer cuando dos familias de curvas son ortogonales; es decir,
cuando cada curva de una familia corta a las curvas pertenecientes a la otra familia en un
ngulo recto.
Definicin 3.1 Cuando todas las curvas de una familia F(x, y, c1 ) = 0 cortan ortogo-
nalmente a todas las curvas de otra familia G(x, y, c2 ) = 0, se dice que las familias son
ortogonales. Esto significa que cada curva de la familia F(x, y, c1 ) = 0 corta en ngulo recto
a toda curva de la otra familia G(x, y, c2 ) = 0.
Ejemplo 3.14 En la Figura 3.3 se ve que la familia de rectas que pasan por el origen
y = c1 x y la familia de crculos concntricos con centro en el origen x2 + y2 = c2 son
trayectorias ortogonales.
90 Algunas Aplicaciones de EDO de Primer Orden
Para encontrar las trayectorias ortogonales de una familia de curvas F(x, y, c1 ) = 0, se halla
dy
de esta ecuacin (por diferenciacin implcita) la derivada dx , que usualmente se puede
expresar de la forma siguente:
dy
= r(x, y)
dx
r(x, y) = F/ x
F/ y , en todo punto (x, y) de F(x, y, c1 ) = 0 donde F
y 6
= 0 .
dy x
=
dx y ln(y)
dy
y ln(y) = x
Z
dx Z
dy
y ln(y) dx = xdx
dx
y ln(y) y2
2 x2
= + c2
2 4 2
2 2 2
2y ln(y) y + 2x = 4c2
2y2 ln(y) y2 + 2x2 = c3
y = cxn , n R. (3.20)
dy
= yxn nxn1 = nx1 y.
dx
dy x
La ecuacin diferencial de la familia de curvas es dx = ny . Esta ltima ecuacin se resuelve
por el mtodo de variable separables, de donde tenemos
dy
ny =x
Z
dx Z
dy
n dx =
y xdx
dx
y2 x2
n = + c.
2 2
ny2 2
La trayectorias ortogonales vienen dada por la frmula 2 + x2 = c.
92 Algunas Aplicaciones de EDO de Primer Orden
Ejercicios
1. Encuentre las trayectorias ortogonales a la familia:
a) x2 + 2xy y2 = k.
b) x2 + (y k)2 = k2 .
c) y = x 1 + c1 ex .
d) y2 = 4px.
e) xy = kx 1.
1
f ) y = x+c 1
.
k
g) y = x , para x > 0, c > 0.
h) xy = c.
i) ex cos(y) = k.
2
j) sin(y) = kex .
2. Determinar el valor de a para que las familias de curvas y3 = C1 x y x2 + ay2 = C2
sean ortogonales.
3. Hallar las trayectorias ortogonales da cada una de las siguientes familias de curvas
a) Una familia de lineas rectas que cortan el origen.
b) Una familia de crculos con radio variable, centradas en el eje x y pasando a
travs del origen.
c) Una familia de elipses con centros en el origen y vrtices en (1, 0).
j) x = cos2 (y).
3. a) x2 + y2 = c.
b) x2 + y2 = cy.
2 2
c) x2 = cex +y .
Captulo 4
Como vimos anteriormente una ecuacin diferencial de orden n lineal es una ecuacin de la
forma:
dny d n1 y dy
an (x) n + an1 (x) n1 + + a1 (x) + a0 (x)y = g(x).
dx dx dx
Obs Desde el punto de vista geomtrico, la unicidad en el teorema permite afirmar que dos
curvas solucin que pasan por el mismo punto deben ser idnticas en todo intervalo
(a, b).
Solucin. Primero notemos que a2 (x) 1 y a0 (x) 1 son funciones continuas, luego
del Teorema de Existencia y Unicidad (Teorema 4.1), la ecuacin y00 y = 0 tiene solucin.
Consideremos y(x) = ex , luego y00 (x) = ex y as tenemos que y(x) = ex es una solucin
de la ecuacin lo que comprueba 1. Para y(x) = ex , tenemos que y00 (x) = ex y as esta
funcin tambin es una solucin de la ecuacin, de donde tenemos 2. Para demostrar 3.
notemos que para y(x) = 1 ex + 2 ex tambin se satisface que y(x) = y00 (x) por lo cual se
tiene que y00 y = 0, es decir, se satisface la ecuacin para esta funcin.
El ejemplo anterior muestra que dada una solucin de y(x) de una ecuacin diferencial
homognea, entonces todo mltiplo de ella y(x) tambin lo es. Ms an, dadas dos
soluciones diferentes y1 (x) y y2 (x), toda combinacin lineal de ellas tambin es solucin
de la ecuacin. Lo anterior nos lleva a preguntar cuando dos (o ms) soluciones de una
ecuacin diferencial son realmente diferentes, para ello usaremos las siguientes nociones.
4.1 Teora Preliminar 97
Ejemplo 4.2 Sea f1 (x) = x2 , f2 (x) = x|x| y consideremos un intervalo abierto I que
contenga al punto x = 0. El Wronskiano de { f1 , f2 } es
2
x x|x|
W [ f1 , f2 ](x) = = x3 |x|0 + x2 |x| 2x2 |x|.
2x x|x|0 + |x|
c1 x2 + c2 |x| = 0.
dny d n1 y dy
an (x) n
+ a n1 (x) n1
+ + a1 (x) + a0 (x)y = 0
dx dx dx
es linealmente independiente en un intervalo I si, y slo si, el Wronskiano W [y1 , y2 , . . . , yn ](x)
no es idnticamente cero en el intervalo I.
Para y1 y y2 soluciones de la ecuacin diferencial lineal de orden 2 homognea,
dny d n1 y dy
an (x) n
+ a n1 (x) n1
+ + a1 (x) + a0 (x)y = 0
dx dx dx
en un intervalo I = [a, b], entonces el Wronskiano W [y1 , y2 , . . . , yn ](x) o es idnticamente
cero en el intervalo I o no es cero para cualquier x I.
Definicin 4.3 Conjunto fundamental de soluciones. Las soluciones y1 , y2 , . . . , yn
de la ecuacin diferencial
dny d n1 y dy
an (x) n + an1 (x) n1 + + a1 (x) + a0 (x)y = 0,
dx dx dx
forman un conjunto fundamental de soluciones en I si
W [y1 , y2 , , yn ](x) 6= 0.
Donde 2e2x es diferente de cero en todo el conjunto de los nmeros reales, por lo tanto del
corolario anterior {y1 , y2 } es un conjunto linealmente independiente.
Ahora consideraremos el caso de diferentes soluciones de una ecuacin diferencial con dis-
tintos coeficientes no homogneos. En este caso tenemos el siguiente importante resultado.
Teorema 4.7 Principio de Superposicin. Sea y1 una solucin de la ecuacin dife-
rencial
dny d n1 y dy
an (x) n
+ a n1 (x) n1
+ + a1 (x) + a0 (x)y = g1 (x)
dx dx dx
y sea y2 una solucin de la ecuacin diferencial
dny d n1 y dy
an (x) n
+ a n1 (x) n1
+ + a1 (x) + a0 (x)y = g2 (x),
dx dx dx
100 Ecuaciones Diferenciales Lineales de Orden Superior
dny d n1 y dy
an (x) n
+ an1 (x) n1
+ + a1 (x) + a0 (x)y = g(x)
dx dx dx
se puede escribir como
dny d n1 y dy
an (x) n
+ a n1 (x) n1
+ + a1 (x) + a0 (x)y = 0
dx dx dx
y Y es una solucin particular de la ecuacin no homognea.
As, para obtener la solucin de una ecuacin diferencial lineal no homognea
dny d n1 y dy
an (x) n
+ an1 (x) n1
+ + a1 (x) + a0 (x)y = g(x)
dx dx dx
se siguen los siguientes dos pasos:
Encontramos la solucin general yh de la ecuacin homognea
dny d n1 y dy
an (x) n
+ a n1 (x) n1
+ + a1 (x) + a0 (x)y = 0 y despus
dx dx dx
encontramos una solucin particular y p de la ecuacin no homognea
dny d n1 y dy
an (x) n
+ an1 (x) n1
+ + a1 (x) + a0 (x)y = g(x).
dx dx dx
De donde se tiene entonces que la solucin general yg de la ecuacin no homognea
viene dada por
yg = yh + y p .
Como 2e2x es diferente de cero para todos los reales, entonces y1 , y2 son linealmente
independientes.
2. Para determinar una solucin general de la ecuacin diferencial, en primer lugar
sustituimos y1 y y2 en y00 2y0 + 5y = 0 para determinar si son soluciones particulares
de la ecuacin. De donde obtenemos las identidades
3ex cos(2x) 4ex sin(2x) 2ex cos(2x) + 4ex sin(2x) + 5ex cos(2x) =0
4ex cos(2x) 3ex sin(2x) 4ex cos(2x) 2ex sin(2x) + 5ex sin(2x) =0.
3. Con las condiciones iniciales y(0) = 2, y0 (0) = 0 podemos determinar una solucin
particular, sustituyendo en la ecuacin diferencial y en su derivada tenemos:
2 =c1
0 =c1 + 2c2 .
Ejercicios
1. En los problemas a) a d), cada familia de funciones es la solucin general de la
ecuacin diferencial en el intervalo indicado. Determinar un miembro de la familia
que sea solucin del problema del valor inicial.
a) y(x) = c1 ex + c2 ex ; (, ); y00 y = 0, y(0) = 0, y0 (0) = 1.
b) y(x) = c1 e4x + c2 ex ; (, ); y00 3y0 4y = 0, y(0) = 1, y0 (0) = 2.
c) y(x) = c1 x + c2 x ln(x); (0, ); x2 y00 xy0 + y = 0, y(1) = 3, y0 (1) = 1.
d) y(x) = c1 + c2 cos(x) + c3 sin(x); (, ); y000 + y0 = 0, y() = 0, y0 () = 2,
y00 () = 1.
2. Si y(x) = c1 + c2 x2 es una familia biparamtrica de soluciones de xy00 y0 = 0 en el
intervalo (, ), demostrar que las constantes c1 y c2 no se pueden determinar de tal
manera que un miembro de la familia satisfaga las condiciones y(0) = 0, y0 (0) = 1.
Explique por qu esto no contradice el Teorema 4.1.
102 Ecuaciones Diferenciales Lineales de Orden Superior
siguientes intervalos?
a) [0, ).
b) (, 0).
c) (, ).
d) Calcular el Wronskiano W [y1 , y2 ](x) en el intervalo (, ).
12. Encontrar las regiones donde el Teorema de existencia y unicidad garantiza la unicidad
para el problema de Cauchy:
y00 y0 + y = sin(x),
y00 y0 + y = e2x ,
dny d n1 y dy
an (x) n + an1 (x) n1 + + a1 (x) + a0 (x)y = 0
dx dx dx
es un espacio vectorial.
16. Demostrar que si y p (x) es una solucin de an (x)y(n) + + a0 (x)y = g(x), entonces
Ay p (x), con A constante, es una solucin de an (x)y(n) + + a0 (x)y = Ag(x).
17. Si y p (x) = u(x) + iv(x), aqu i = 1, es una solucin particular de an (x)y(n) + +
a0 (x)y = R(x) + iS(x) donde ai (x) son funciones reales, entonces
a) La parte real de y p (x); esto es u(x), es una solucin de an (x)y(n) + + a0 (x)y =
R(x).
b) La parte imaginaria de y p (x); esto es v(x), es una solucin de an (x)y(n) + +
a0 (x)y = S(x).
9. b) Linealmente independiente.
c) Linealmente independiente.
d) Linealmente dependiente.
10. a) y(x) = 3x2 + x1 .
c) y(x) = 5et1 t 2 2t 2.
11. b) No son linealmente independientes.
d) 0.
13. b) y(x) = cos(x) e2x .
c) y(x) = 4 cos(x) + 6e2x .
Caso II: Las races del polinomio son reales pero algunas con multiplicidad mayor que uno.
Primero supongamos que tenemos una ecuacin de orden dos ay00 + by0 + cy = 0
cuyo polinomio caracterstico tiene dos races iguales, entonces y1 (x) = emx es una
solucin, para encontrar la otraR utilizamos el mtodo de reduccin de orden, con
Z p(x)dx
e
la frmula, y2 (x) = y1 (x) dx en este caso p(x) = ba y m = b
2a ya que
[y1 (x)]2
b2 4ac = 0. Sustituyendo en la ecuacin obtenemos,
R b
Z
mx e a dx
y2 (x) = e b 2 dx
e 2a x
calculando la integral obtenemos y2 (x) = xemx . Estas dos soluciones son linealmente
independientes y la solucin general de la ecuacin es
m1 =a + ib,
m2 =a ib.
Las soluciones fundamentales son y1 (x) = e(a+ib)x , y2 (x) = e(aib)x . Puede tomarse
la expresin de estas exponenciales segn la frmula de Euler:
1. y000 6y00 y0 + 6y = 0.
2. 4y00 + 20y0 + 25y = 0.
3. y000 y00 + y0 + 3y = 0.
Solucin. 1. Si tratamos de hallar la solucin de la forma y(x) = emx , entonces esto nos
conduce a hallar las races de la ecuacin auxiliar
m3 6m2 m + 6 = 0,
y(x) = c1 e6x + c2 ex + c3 ex .
2. Si tratamos de hallar la solucin de la forma y(x) = emx , esto se reduce a hallar las
races de la ecuacin auxiliar
4m2 + 20m + 25 = 0.
m3 m2 + m + 3 = 0,
m1 = 1
2 4 12 2 8
m2 = = = 1 i 2.
2 2
Aqu = 1, = 2, y as la solucin general es
y(x) = c1 ex + c2 ex cos( 2x) + c3 ex sin( 2x).
4.3 Ecuaciones Diferenciales no Homogneas 107
Ejercicios
1. Resolver las siguientes ecuaciones diferenciales homogneas:
a) 3y00 + 11y0 7y = 0.
b) y000 + y00 5y0 + 3y = 0.
c) y000 + y00 + 3y0 5y = 0.
d) 4y00 + 4y0 + 6y = 0.
e) 4y00 4y0 + 26y = 0.
f ) y00 + 2y0 + 2y = 0; y(0) = 2, y0 (0) = 1.
g) y000 3y00 + 3y0 y = 0.
h) y000 y00 + y0 + 3y = 0.
i) y000 + 2y00 5y0 6y = 0; y(0) = 1, y0 (0) = 0, y00 (0) = 0.
j) y000 4y00 + 7y0 6y = 0; y(0) = y0 (0) = 0, y00 (0) = 1.
k) y00 + 16y = 0; y(0) = 2, y0 (0) = 2.
l) y(4) + y000 + 3y00 = 0.
m) y(4) 2y00 + y = 0.
n) y00 + y0 + 2y = 0; y(0) = y0 (0) = 0.
) y00 2y0 + y = 0; y(0) = 1, y0 (0) = 2.
o) y00 2y0 + 2y = 0; y(0) = 1, y0 () = 1.
2. Examinar la ecuacin de segundo orden con coeficientes constantes y00 + by0 + cy = 0.
Si y(x) es una solucin de la ecuacin describa qu condiciones deben satisfacer b y
c para que lmx y(x) = 0.
3. Para ver el efecto de cambiar el parmetro b en el problema con valores iniciales
g) y(x) = c1 ex + c2 xex + c3 x2 ex .
i) y(x) = 15 e3x + 15 e2x + ex .
k) y(x) = 12 sin(4x) + 2 cos(4x).
m) y(x) = c1 ex + c2 xex + c3 ex + c4 xex .
n) y(x) = 0.
y(x) = ex cos(2x) + (e 21) ex sin(2x).
o)
y00 + 2y0 + y = ex ,
dado que g(x) = ex , entonces la funcin candidata a ser solucin particular de la ecua-
cin tendr la forma y p (x) = ex . Ahora debemos determinar el coeficiente , para ello
derivamos esta funcin y sustituimos en la ecuacin, as tenemos que
ex + 2ex + ex = ex ,
g(x) yp (x)
1. pn (x) = an xn + + a
1 x + a0
s s n
x Pn (x) = x (An x + + A1 x + A0 )
2. aex xs Aex
3. a cos( x) + b sin( x) xs (A cos( x) + B sin( x))
4. pn (x)ex xs (An xn + + A1 x + A0 )ex
5. pn (x) cos( x) + qm (x) sin( x), don- xs {PN (x) cos( x) + QN (x) sin( x)}, don-
de, qm (x) = bm xm + + b1 x + b0 . de QN (x) = BN xN + + B1 x + B0 y N =
max(m, n).
6. aex cos( x) + bex sin( x) xs (Aex cos( x) + Bex sin( x))
7. ex [pn (x) cos( x) + qm (x) sin( x)] xs {ex (PN (x) cos( x) + QN (x) sin( x))},
y N = max(m, n).
Tabla 4.1: Coeficiente no homogneo para el mtodo de los coeficientes indeterminados y
la correspondiente forma de la solucin particular. Aqu el entero no negativo s es el menor
entero tal que ningn trmino de la solucin y p sea solucin de la ecuacin homognea.
4.3 Ecuaciones Diferenciales no Homogneas 109
Simplificando tenemos,
1
y = c1 e2x + c2 xe2x + x3 e2x .
6
2. Primero resolvemos la ecuacin homognea y00 + 5y0 + 6y = 0 y obtenemos como
solucin
yh (x) = c1 e3x + c2 e2x .
En este caso g(x) = sin(x) cos(2x) que no es de las formas dadas en la Tabla
4.1, pero es la suma de dos de tales trminos. Entonces utilizamos el principio de
superposicin y consideramos por separado las ecuaciones
Simplificando tenemos
y concluimos que
1 1
A= y B= .
10 10
1 1
Por lo tanto, y p1 (x) = 10 cos(x) + 10 sin(x). Ahora resolvemos la ecuacin (4.3)
para la cual g2 (x) = cos(2x) es decir
Simplificando tenemos,
y concluimos que
1 5
C= y D= .
52 52
1 5
Por lo tanto y p2 (x) = 52 cos(2x) 52 sin(2x), es decir la solucin general es
1 1 1 5
y(x) = C1 e3x +C2 e2x cos(x) + sin(x) cos(2x) sin(2x).
10 10 52 52
4.3 Ecuaciones Diferenciales no Homogneas 111
Ejercicios
1. Determinar una solucin general de la ecuacin diferencial dada.
a) y00 8y0 + 16y = 6xe4x + 2 + 16x + 16x2 .
b) y00 + 4y0 + 3y = ex (2 + 8x); y(0) = 1, y0 (0) = 2.
c) y00 y0 2y = cos(x) sin(2x); y(0) = 20 7
, y0 (0) = 15 .
d) y00 4y0 + 4y = x2 e2x e2x .
e) y00 3y0 + 2y = e3x (1 + x).
f ) y00 + 2y0 + y = x2 + 1 ex ; y(0) = 0, y0 (0) = 2.
g) y00 + 6y0 + 13y = e3x cos(2x).
h) y000 + 8y00 = 6x2 + 9x + 2
i) y00 + y0 + 41 y = ex (sin(3x) cos(3x)).
j) y00 + 5y0 6y = 10e2x ; y(0) = y0 (0) = 1.
k) y000 2y00 4y0 + 8y = 6xe2x .
l) y00 + 3y0 2y = e2x ((4 + 20x) cos(3x) + (26 32x) sin(3x)).
m) y(4) y00 = 4x + 2xex .
n) y000 + 8y = 2x 5 + 8e2x ; y(0) = 5, y0 (0) = 3, y00 (0) = 4.
) y00 4y0 + 4y = 6e2x + 25 sin(x); y(0) = 5, y0 (0) = 3.
o) y00 4y0 + 4y = e2x (1 + x) + e2x (cos(x) sin(x)) + 3e3x + 1 + x.
2. Mtodo del anulador. El mtodo del anulador permite justificar el mtodo de
coeficientes indeterminados. Sabemos que una ecuacin diferencial se puede es-
cribir de la forma an Dn y + an1 Dn1 y + + a1 Dy + a0 y = g(x) en donde Dk =
dk y
dxk
, k = 1, 2, , n. Cuando nos convenga, representaremos tambin la ecuacin de
la forma L(x) = g(x), donde L representa el operador diferencial lineal de orden n:
L = an Dn + an1 Dn1 + + a1 D + a0 . Se dice que L es un anulador si L( f (x)) = 0.
Ejemplo Dn anula a cada una de las funciones 1, x, x2 , , xn1 , (D )n anula las
funciones ex , xex , x2 ex , , xn1 ex , [D2 2D + ( 2 + 2 )]n anula las funciones
ex cos( x), xex cos( x), x2 ex cos( x), , xn1 ex cos( x), ex sin( x),
xex sin( x), x2 ex sin( x), , xn1 ex sin( x).
Para ilustrar considere la ecuacin:
c1 e3x + c2 ex + c3 ex (4.6)
es una solucin general de (4.5) y concluya que existe una solucin particular
de (4.4) de la forma (4.6).
112 Ecuaciones Diferenciales Lineales de Orden Superior
dny d n1 y dy
an (x) n
+ an1 (x) n1
+ + a1 (x) + a0 (x)y = f (x)
dx dx dx
y una solucin no trivial y1 (x) = g(x), el cambio de variable y(x) = g(x)u(x) transforma la
ecuacin diferencial dada en una ecuacin lineal de orden n 1 en la variable w = du dx .
Primero veamos el caso n = 2 y f (x) 0.
Sea,
d2y dy
a2 (x) + a 1 (x) + a0 (x)y = 0 (4.7)
dx2 dx
y una solucin no nula y1 (x) = g(x). Tenemos y2 (x) = g(x)u(x) la segunda solucin enton-
ces:
Podemos ver que el primer corchete toma el valor cero por ser g(x) solucin de la ecuacin
homognea, entonces nos queda
a2 (x)gu00 + 2a2 (x)g0 u0 + a1 (x)gu0 = 0.
Ahora llamamos w = u0 , la ecuacin queda en la forma
a2 (x)gw0 + 2a2 (x)g0 w + a1 (x)gw = 0.
Esta ecuacin es de variables, as integrando respecto a x
w0 2a2 (x)g0 a1 (x)
= +
w a2 (x)g a2 (x)
Z Z Z
1 dw 2g0 a1 (x)
dx = dx dx
w dx g a2 (x)
Z
2 a1 (x)
ln |w| = ln(g ) dx
a2 (x)
Z
a1 (x)
dx
w =g2 e a2 (x) .
De manera similar, el cambio y = uy1 nos permite hallar la solucin general de una ecuacin
no homognea
dny d n1 y dy
an (x) n
+ an1 (x) n1
+ + a1 (x) + a0 (x)y = f (x)
dx dx dx
si conocemos una solucin y1 de la correspondiente ecuacin homognea
dny d n1 y dy
an (x) n
+ a n1 (x) n1
+ + a1 (x) + a0 (x)y = 0.
dx dx dx
Ejemplo 4.8 Hallar la solucin general de la ecuacin
Solucin. Sea y(x) = u(x)ex , luego tenemos que y0 = u0 ex + uex y y00 = u00 ex + 2u0 ex + uex .
Reemplazamos esto en el lado izquierdo de la ecuacin y obtenemos
xy00 (2x + 1)y0 + (x + 1)y =x(u00 ex + 2u0 ex + uex ) (2x + 1)(u0 ex + uex )
+ (x + 1)(uex )
=(xu00 u0 )ex .
(xu00 u0 )ex = x2 .
(xz0 z)ex = x2 .
Solucin. Sea y(x) = u(x)y1 (x) = u(x)x. Tenemos que y0 = u0 x + u y y00 = u00 x + 2u0 . Re-
emplazamos en la ecuacin
Ejercicios
1. Se da una ecuacin diferencial y una solucin no trivial f . Determinar una segunda
solucin linealmente independiente.
a) x2 y00 xy0 + 2y = 0; f (x) = x sin(ln |x|).
b) y00 25y = 0; f (x) = e5x .
c) x2 y00 + 2xy0 6y = 0; f (x) = x2 .
d) xy00 + (1 2x)y0 + (x 1)y = 0; x > 0; f (x) = ex .
e) (1 2x x2 )y00 + 2(1 + x)y0 2y = 0; f (x) = x + 1.
f ) 4x2 sin(x)y00 4x(x cos(x) + sin(x))y0 + (2x cos(x) + 3 sin(x))y = 0; f (x) =
x1/2 .
116 Ecuaciones Diferenciales Lineales de Orden Superior
donde es un parmetro: Use la frmula para reduccin de orden para obtener una
representacin integral de una segunda solucin linealmente independiente de la
ecuacin de Legendre para el valor dado de y la solucin correspondiente de f (x).
a) = 1, f (x) = x.
b) = 2, f (x) = 3x2 1.
c) = 3, f (x) = 5x3 3x.
3. La ecuacin xy000 + (1 x)y00 + xy0 y = 0, tiene a f (x) = x como solucin. Use la
sustitucin y(x) = v(x) f (x) para reducir esta ecuacin de tercer orden a una ecuacin
lineal homognea de segundo orden en la variable w = v0 .
4. Forma normal. Muestre que la sustitucin y(x) = u(x)v(x), donde
Z
1
v(x) = exp p(x)dx ,
2
digamos
yh = c1 y1 + c2 y2 + + cn yn .
El mtodo de variacin de parmetros propone una solucin particular de la forma
de donde y1 (x) = cos(x), y2 (x) = sin(x), W1 (x) = (1)21W [y2 ](x) = sin(x) y
W2 (x) = (1)22W [y1 ](x) = cos(x). El Wronskiano de {y1 , y2 } es:
cos(x) sin(x)
W [y1 , y2 ](x) = = cos2 (x) + sin2 (x) = 1.
sin(x) cos(x)
Por lo tanto,
Z Z
tan2 (x) sin(x) tan2 (x) cos(x)
u1 (x) = dx, u2 (x) = dx
1 1
(cos2 (x) + 1)
u1 (x) = , u2 (x) = ln | sec(x) + tan(x)| sin(x)
cos(x)
es decir,
(cos2 (x) + 1)
y p (x) = cos(x) + (ln | sec(x) + tan(x)| sin(x)) sin(x)
cos(x)
y p (x) = sin(x)(ln | sec(x) + tan(x)|) 2
Ahora
y p (x) = u1 (x)e2x + u2 (x)xe2x
con y1 (x) = e2x , y2 (x) = xe2x , W1 (x) = (1)21W [y2 ](x) = xe2x y W2 (x) =
(1)22W [y1 ](x) = e2x . El Wronskiano de {y1 , y2 } es:
2x
e xe2x
W [y1 , y2 ](x) = = e4x 2xe4x + 2xe4x = e4x .
2e2x e 2xe
2x 2x
Por lo tanto,
Z Z
e2x ln |x|xe2x e2x ln |x|e2x
u1 (x) = 4x
dx, u 2 (x) = dx
Ze Z e4x
u1 (x) = x ln |x|dx, u2 (x) = ln |x|dx
2x2 ln |x| + x2
u1 (x) = , u2 (x) = x ln |x| x.
4
Es decir,
2x2 ln |x| + x2 2x
y p (x) = e + (x ln |x| x)xe2x
4
4.3 Ecuaciones Diferenciales no Homogneas 119
e2x x2 (2 ln |x| 3)
y p (x) =
4
luego, la solucin general de la ecuacin es:
e2x x2 (2 ln |x| 3)
y(x) = c1 e2x + c2 xe2x + .
4
y000 2y00 y0 + 2y = 0
2x
e ex
W2 (x) = (1)32W [y1 , y3 ](x) = 2x x = e3x ,
2e e
2x
e ex
W3 (x) = (1) 33
W [y1 , y2 ](x) = 2x = 3ex .
2e e
x
El Wronskiano de {y1 , y2 , y3 } es
2x
e x x
2x e x ex
W [y1 , y2 , y3 ](x) = 2e e e = 6e2x .
4e2x ex ex
Ahora calculamos
Z
2e3x 1
u1 (x) = 2x
dx = ex ,
6e 3
Z 3x 3x
e e 1
u2 (x) = 2x
dx = e4x ,
6e 24
Z x 3x
3e e 1
u3 (x) = 2x
dx = e2x .
6e 4
120 Ecuaciones Diferenciales Lineales de Orden Superior
es
2 19 3 1
y(x) = e2x + ex + ex + e3x .
3 24 4 8
Ejercicios
1. Resolver cada una de las siguientes ecuaciones diferenciales:
a) y00 + y = tan(x).
b) y00 y = cosh(x).
2x
c) y00 4y = ex .
d) y00 + 3y0 + 2y = sin(ex ).
e) y000 + y0 = tan(x).
2
f ) y00 2y0 + y = 14x 3 ex .
4ex
g) y00 y = 1e 2x .
h) y y + y y = ex sin(x).
000 00 0
+ x sin(x)
4 (sin(x) cos(x) x).
3. a) y(x) = c1 x1/2 ex + c2 x1/2 ex + ln |x|x1/2 ex .
b) y(x) = c1 x2 + c2 x3 x 61 .
d) y(x) = ex csc(x) x
+ xe xcos(x)
cot(x).
f) y(x) = c1 x + c2 xe + 4 xex .
5. a = 1, y(x) = x3 + c1 x2 + c2 .
m1 =a + ib
m2 =a ib.
x3 y000 6y = 0, x > 0.
Solucin. Si tratamos de hallar la solucin de la forma y(x) = xm , entonces esto nos conduce
a hallar las races de la ecuacin auxiliar
m(m 1)(m 2) 6 =0
m3 3m2 + 2m 6 =0.
4.3 Ecuaciones Diferenciales no Homogneas 125
dny n1 d
n1 y dy
an x n n
+ a n1 x n1
+ + a1 x + a0 y = g(x), x>0
dx dx dx
esta ecuacin se puede reducir a una ecuacin con coeficientes constantes mediante el
cambio de variable x = et , esto nos permite hallar la solucin general de esta ecuacin sin
necesidad de conocer previamente alguna otra solucin.
x2 y00 + 3xy0 + y = x1 .
t = ln(x)
dt 1
=
dx x
dy dy dt dy 1
= =
dx dt dx dt x
d 2 y 1 d 2 y dt 1 dy 1 d 2 y dy
= 2 = 2 .
dx2 x dt 2 dx x dt x dt 2 dt
dy dy
x = (4.9)
dx dt
2 d 2 y dy
2d y
x = . (4.10)
dx2 dt 2 dt
d 2 y dy dy
+ 3 + y(t) = et
dt 2 dt dt
d2y dy
2
+ 2 + y(t) = et ,
dt dt
esta ltima ecuacin es lineal de segundo orden con coeficientes constantes no homognea
la cual sabemos resolver. Primero resolvemos la ecuacin homognea
y00 + 2y0 + y = 0
por lo tanto,
Z Z
et tet et et
u1 (t) = dt, u2 (t) = dt
e2tZ Z e
2t
u1 (t) = t, u2 (t) = dt
t 2
u1 (t) = , u2 (t) = t
2
es decir,
t 2 t 2 t
y p (t) = e +t e
2
t 2 et
y p (t) = .
2
As, la solucin general de la ecuacin con coeficientes constantes es
t t t 2 et
y(t) = c1 e + c2te + .
2
Sustituimos x = et y se obtiene la solucin general de la ecuacin de Cauchy-Euler
ln2 (x)x1
y(x) = c1 x1 + c2 ln(x)x1 + .
2
d3y d2y dy d 2 y dy dy
3 + 2 + + =t
dt 3 dt 2 dt dt 2 dt dt
d3y d2y dy
4 + 4 =t, (4.12)
dt 3 dt 2 dt
la cual es una ecuacin no homognea con coeficientes constantes. Para hallar su solucin
general primero hallemos la solucin de la ecuacin homognea. La ecuacin auxiliar de la
ecuacin homognea es
as, el conjunto fundamental de soluciones es {y1 (t) = 1, y2 (t) = e2t , y3 (t) = te2t }. Hallare-
mos una solucin particular de la ecuacin no homognea usando el mtodo de variacin de
parmetros. Calculemos el Wronskiano
1 e2t te2t
W [y1 , y2 , y3 ](t) = 0 2e (2t + 1)e = 4e4t .
2t 2t
0 4e2t (4t + 4)e2t
1 ln(x) 1
y0 (x) =2c2 x + c3 (x + 2x ln(x)) + + ,
4 x
4x
00 1 1 ln(x) 1
y (x) =2c2 + c3 (3 + 2 ln(x)) + ,
4 x2 4x2
0 = y(1) =c1 + c2
1
1 = y0 (1) =2c2 + c3 +
4
1 = y00 (1) =2c2 + 3c3
cuya solucin es c1 = 23/8, c2 = 23/8 y c3 = 9/8. El problema de valor inicial tiene por
solucin
23 23 2 9 2 1 1
y(x) = x + x ln(x) + ln2 (x) + ln(x).
8 8 8 8 4
Ejercicios
1. En los ejercicios a) y b) disee una modificacin del mtodo de ecuaciones de Cauchy-
Euler para determinar una solucin general de la ecuacin dada.
a) (x 2)2 y00 (x) 7(x 2)y0 (x) + 7y(x) = 0, x > 2.
b) (x + 1)2 y00 (x) + 10(x + 1)y0 (x) + 14y(x) = 0, x > 1.
2. Determinar una solucin general de la ecuacin diferencial dada (considere x > 0
cuando sea necesario)
a) x2 y00 + 7xy0 7y = 0.
b) y00 1x y0 + x52 y = 0.
c) x2 y00 3xy0 + 6y = 0.
d) x2 y00 3xy0 2y = 0.
e) x2 y00 + 5xy0 + 4y = 0.
4.3 Ecuaciones Diferenciales no Homogneas 129
f ) x2 y00 + 3xy0 + 5y = 0.
g) x2 y00 + 9xy0 + 17y = 0.
h) x3 y000 + xy0 y = 0.
i) x2 y00 5xy0 + 8y = 0; y(2) = 32; y0 (2) = 0.
j) y000 = 24(x+y)
x3
.
k) x y x y + 2xy0 2y = x3 .
3 000 2 00
1 5
x2 y00 + xy0 + (x2 )y = x 2 , x>0
4
1
tiene dos soluciones linealmente independientes, y1 (x) = x 2 cos(x),
1
y2 (x) = x 2 sin(x). Determinar una solucin general de la ecuacin no homognea.
d2x dx
2
+ 2r + 2 x = f (t),
dt dt
con f una funcin continua en un intervalo I, r y constantes positivas.
d2x
+ 2x = 0
dt 2
donde es una constante positiva y x describe la posicin de la partcula como una funcin
del tiempo t. Las races de la ecuacin auxiliar asociada m2 + = 0 son r1,2 = i y por lo
tanto la solucin de la ecuacin diferencial es
As, la solucin de la ecuacin del movimiento armnico simple puede ser escrita en
cualquiera de las siguientes formas:
Ejemplo 4.16 Una partcula ejecuta un movimiento armnico simple, la frecuencia natural
de este movimiento es 4 rad/seg. Si la partcula comienza desde el punto de equilibrio con
una velocidad de 4 mts/seg hallar:
1. La ecuacin del movimiento de la partcula.
2. La amplitud del movimiento.
3. El ngulo fase.
4. El periodo del movimiento.
5. La frecuencia del movimiento en ciclos por segundo.
0 = A cos( ), 4 = 4A sin( ).
Ejemplo 4.17 Una partcula que se mueve en linea recta es atrada haca el origen por
una fuerza F. Si la fuerza de atraccin es proporcional a la distancia x de la partcula desde
el origen, demostrar que la partcula ejecuta un movimiento armnico simple. Describa el
movimiento.
F = kx
A P2 o P1 A
d2x
F =m = kx,
dt 2
d2x k
2
= x.
dt m
4.4 Aplicaciones: El Oscilador Armnico 133
x(t) = A sin(t + ).
F = k`
k` = mg.
134 Ecuaciones Diferenciales Lineales de Orden Superior
0 0 0
y=0 y=0 posicion de equilibrio
del resorte
+ y
Por la segunda ley de Newton, la fuerza total del sistema es igual a la masa por su aceleracin
la cual es una ecuacin diferencial de segundo orden, pues la aceleracin es la segunda
derivada de la posicin y con respeto al tiempo. Resulta (tomando la direccin positiva haca
abajo) que el desplazamiento y(t) del resorte sigue la ecuacin diferencial
d2y
m =mg k(` + y)
dt 2
= ky.
O, equivalentemente,
d2y k
+ y = 0. (4.15)
dt 2 m
Escribimos la ecuacin auxiliar asociada a (4.15) que es z2 + mk z = 0, cuyas races son
q
complejas conjugadas mk i. Por lo tanto la solucin general de (4.15) es:
r r
k k
y(t) =c1 cos t + c2 sin t , (4.16)
m m
r
k
=A sin t + . (4.17)
m
q
c1
Con amplitud de movimiento A = c21 + c22 , ngulo fase dado por tan( ) = c2 . El movi-
2 k/m
miento es peridico con periodo y frecuencia natural 2 .
k/m
4.4 Aplicaciones: El Oscilador Armnico 135
Ejemplo 4.18 Una masa de 3 Kg est unida a un resorte con rigidez k = 48. La masa se
1
desplaza mts a la izquierda del punto de equilibrio y recibe una velocidad de 2 mts/seg
2
hacia la derecha. La fuerza de amortiguamiento es despreciable. Determinar la ecuacin
de movimiento de la masa junto con su amplitud, periodo y frecuencia. Cunto tiempo
despus de su liberacin la masa pasa por su posicin de equilibrio?
Solucin. Tenemos un caso de vibracin libre no amortiguada, la ecuacin de movimiento
que usamos es (4.15). La frecuencia angular es
r r
k 48 radianes
= = =4 ,
m 3 seg
sustituyendo en (4.16) obtenemos
y(t) = c1 cos(4t) + c2 sin(4t).
Ahora usamos las condiciones iniciales y(0) = 12 mts, y0 (0) = 2 mts/seg para hallar c1 y
c2 . Reemplazamos las condiciones iniciales en la ltima ecuacin y obtenemos
1
=c1
2
2 =4c2
de modo que 21 = c1 y c2 = 12 . Por lo tanto la ecuacin de movimiento de la masa es
1 1
y(t) = cos(4t) + sin(4t).
2 2
q q 2
2 2
Como A = c21 + c22 , tenemos que A = 21 + 12 , luego la amplitud es A = 2 mts.
Para hallar el ngulo fase tomamos tan( ) = cc21 , reemplazando tenemos tan( ) = 1/2
1/2 = 1,
2
de donde = arctan(1) = 4 radianes que es el ngulo fase. El periodo es 4 = 2 seg y
la frecuencia
4 2 ciclos
= = .
2 2 seg
Por ltimo, para determinar el momento en que la masa pasa por su posicin de equilibrio
hacemos y = 0. Reemplazando en la ecuacin tenemos,
1 1 2
y(t) = cos(4t) + sin(4t) = sin(4t + )
2 2 2
2
0= sin(4t + ).
2
Ahora resolvemos la ecuacin la cual se satisface si 4t + = n o t = n+/4
4 , donde n es
un entero. Si tomamos que es la primera vez que la masa pasa por la posicin de equilibrio
obtenemos
1 + /4 5
t= = seg.
4 16
136 Ecuaciones Diferenciales Lineales de Orden Superior
tenemos
dy
y(t) = cos( + 4t), v(t) = = 4 sin( + 4t).
dt
Cuando y = 12 obtenemos de la primera ecuacin
1
= cos( + 4t), + 4t = arc cos(1/2).
2
Luego, cuando y = 1
2, sin( + 4t) = 23 .Finalmente, concluimos que
v(t) = 2 3,
la cual es la velocidad del cuerpo cuando pasa por el punto y = 12 . El signo + indica que el
cuerpo se mueve haca abajo mientras que el signo indica que lo hace haca arriba.
Ejercicios
1. Una partcula ejecuta un movimiento armnico simple con 2 por segundo. Si la
partcula comienza desde la posicin x = 4 mts con una velocidad de 4 mts/seg
hallar:
a) La ecuacin del movimiento.
b) La amplitud del movimiento.
c) La frecuencia del movimiento.
d) El ngulo fase.
e) El tiempo t cuando la partcula atraviesa por primera vez la posicin de equili-
brio.
2. Una partcula ejecuta un movimiento armnico simple. En t = 0 su velocidad es 0 y
est a 5 mts de la posicin de equilibrio. En t = 14 su velocidad de nuevo es 0 y su
posicin una vez ms es 5 mts de la posicin de equilibrio.
a) Hallar la posicin y velocidad como funcin del tiempo.
b) Hallar la frecuencia y amplitud.
c) Cundo y con qu velocidad la partcula atraviesa la posicin de equilibrio?
3. Una partcula ejecuta un movimiento armnico simple y al final de cada 34 de segundo
pasa por la posicin de equilibrio con una velocidad de 8 mts/seg.
4.4 Aplicaciones: El Oscilador Armnico 137
9. a) 0,16
mts, -5,37 mts/seg.
2
b) 32 seg, -5,66 mts/seg.
c) 2 6 mts/seg.
d) 0,31 seg, -0,47 mts.
2
e) 8 seg, 8 2 rad/seg, 12 mts.
As, el grfico que representa a la solucin x(t) cruza por el origen solo una vez,
adems x(t) 0 cuando t . Calculando dx dt tenemos
dx
= c1 et + c2 et , , < 0.
dt
Esta ecuacin tiene la misma forma que x(t) y por lo tanto solo tiene un punto t tal
que x0 (t ) = 0, es decir la curva determinada por x(t) tiene a lo ms un mximo y un
mnimo. Esto significa que el movimiento no es oscilatorio y termina con el tiempo.
En este caso la fuerza de amortiguamiento r es mayor que la fuerza de restitucin y
as evita la oscilacin. El sistema es llamado sobreamortiguado.
Caso II: r2 2 = 0. Las dos races z1 y z2 son iguales. Esto quiere decir que la fuerza de
amortiguamiento r es igual a la causada por la fuerza de restitucin . La solucin
del sistema es:
Dado que r > 0, entonces ert , tert 0 cuando t y, como en el caso anterior,
existe un nico valor t tal que x0 (t ) = 0, as el movimiento no es oscilatorio y
termina con el tiempo. Decimos que el movimiento es crticamente amortiguado
ya que una pequea disminucin de la fuerza de amortiguamiento provocara un
movimiento oscilatorio.
2 2
r < 0. En este caso las races son imaginarias conjugadas z1,2 = r
Caso III:
2 r2 i. La ecuacin resultante del movimiento es:
p p
rt
x(t) = e c1 cos 2 r2t + c2 sin 2 r2t ,
140 Ecuaciones Diferenciales Lineales de Orden Superior
d2x dx
m 2
+ b + kx = 0.
dt dt
b k
Dividiendo por m, 2 = m y haciendo w2 = m tenemos que
d2x dx
+ 2 + w2 x = 0.
dt 2 dt
Ejemplo 4.20 El movimiento de un sistema masa-resorte con amortiguamiento est
descrito por
1
y(0) = , y0 (0) = 2.
y00 (t) + by0 (t) + 64y(t) = 0,
2
Determinar la ecuacin de movimiento y bosqueje la grfica para b = 0, 10, 16 y 20.
m2 + bm + 64 = 0,
b b2 256 b b2 256
m= = .
2 2 2
1. b = 0. En este caso las raices de la ecuacin son complejas (puras)
16i
m1 = = 8i
2
m2 = 8i.
Ahora usamos las condiciones iniciales y(0) = 1/2, y0 (0) = 2 para hallar c1 y c2 ,
reemplazamos las condiciones iniciales en la ltima ecuacin y obtenemos
c1 = 1,
c2 = 0.
y(t) = cos(8t)
q
donde A = c21 + c22 = 1 es el factor de amortiguamiento.
142 Ecuaciones Diferenciales Lineales de Orden Superior
Figura 4.7: y(t) = sin 8t + 2 .
78 5t
Figura 4.8: y(t) = 39 e (sin( 39t) + ).
3. b = 16. En este caso las races son reales y se repiten m1 = m2 = 8. La ecuacin del
movimiento (crticamente amortiguado) es
1
y(t) = e8t 2te8t .
2
El grfico que describe el movimiento de la masa se da en la Figura 4.9:
Solucin. La segunda ley de Newton, debido al factor de resistencia, en este caso queda
como
d2y 1 dy
m 2 = ky .
dt 2 dt
Como 2 libras estiran el resorte 32 pulgadas (1 pie=12 pulgadas), entonces 32 pulgadas = 83
pies. As, dado que ky = mg tenemos
8
k = 2, g = 32 pies/seg2
3
4.4 Aplicaciones: El Oscilador Armnico 145
ky 1
m= = .
g 16
1 d 2 y 1 dy 3 d2y dy
+ + y = 0, + 8 + 12y = 0
16 dt 2 2 dt 4 dt 2 dt
cuya solucin es
Las condiciones iniciales son y(0) = 1, y0 (0) = 0. Usando las condiciones iniciales nos
queda el sistema
1 =c1 + c2
0 = 2c2 6c2
3
cuya solucin es c1 = 2 y c2 = 21 . Luego, ecuacin del movimiento del objeto es
3 1
y(t) = e2t e6t .
2 2
3
Ejemplo 4.22 Si el coeficiente de resistencia en el ejemplo anterior es 8 en vez de 12 ,
hallar:
1. La ecuacin del movimiento del sistema.
2. El factor de amortiguamiento.
3. La amplitud amortiguada del movimiento.
4. El periodo amortiguado del movimiento.
5. La frecuencia amortiguada del movimiento.
6. La constante del tiempo.
1 d 2 y 3 dy 3 d2y dy
+ + y = 0, + 6 + 12y = 0.
16 dt 2 8 dt 4 dt 2 dt
La solucin es
y(t) = Ae3t sin( 3t + )
y la velocidad es
y0 (t) = 3Ae3t sin( 3t + ) + 3Ae3t cos( 3t + ).
146 Ecuaciones Diferenciales Lineales de Orden Superior
1 =A sin( )
0 = 3A sin( ) + 3A cos( ).
1
Luego, A = sin( ) y
3
0 = 3 + 3 cot( ), cot( ) = = 3,
3
de donde tenemos
7 1
= , sin( ) = .
6 6 2
Si suponemos sin( ) = 21 , tenemos que A = 2 y la ecuacin del movimiento queda
3t
y(t) = 2e sin 3t + .
6
Ejercicios
1. Una partcula se mueve en linea recta segn la ley
d2x dx
2
+ 2r + x = 0,
dt dt
donde r es constante y x es el desplazamiento de la partcula desde su punto de
equilibrio.
a) Para qu valores de r el sistema ser subamortiguado, crticamente amortiguado,
sobreamortiguado?
b) Compruebe sus respuesta resolviendo la ecuacin con r = 1, 12 , 12 , 1, 2.
c) Para cual valor de r el movimiento ser oscilatorio con un periodo amortiguado
de 3?
d) Existe un valor de r que haga el periodo amortiguado menor que 2?
2. Una partcula ejecuta un movimiento armnico amortiguado. En 10 segundos el factor
de amortiguamiento decrece en un 80 %. Su periodo de amortiguamiento es 2 seg.
Hallar la ecuacin diferencial del movimiento.
4.4 Aplicaciones: El Oscilador Armnico 147
Z
T f (x) = K(t, s) f (t)dt.
I
siempre que esta integral impropia converja. Ntese la transformada de Laplace aplicada a
una funcin es de nuevo una funcin, a saber, L( f ) F.
Solucin.
Z Z N
L( f ) = 1est dt = lm 1est dt
0 N 0
t=N
1 st 1 sN 1
= lm e = lm e +
N s t=0 N s s
1
= , siempre que s > 0.
s
Ejemplo 5.3 Determinar la transformada de Laplace de la funcin dada por f (t) = cos(2t).
Solucin.
Z
L(cos(2t)) = cos(2t)est dt
0
Z N
= lm cos(2t)est dt
N 0
2 sin(2t)est s cos(2t)est N
= lm
N s2 + 4 s2 + 4 0
2 sin(2N)esN s cos(2N)esN s
= lm + 2 .
N s2 + 4 s2 + 4 s +4
2 sin(2N)esN sN
Como s2 +4
s cos(2N)e
s2 +4
0 cuando s > 0 y N , concluimos que
s
F(s) = 2 , para s > 0.
s +4
5.1 Transformada de Laplace: Definicin y Propiedades 151
Solucin.
Z 3 Z
L( f (t)) = e2t est dt + 1est dt
0 3
Z 3 Z N
2t st
= e e dt + lm est dt
0 N 0
t(2s) 3 st N
e
= + lm e
2 s 0 N s 3
e3(2s) 1 esN es3
= + lm .
2s 2 s N s s
eNs
Como s 0 cuando s > 0 y N , obtenemos:
1 e3(2s) e3s
F(s) = + , para s > 2.
s2 s
Solucin.
Z st Z
L [t] = [t]e dt = [t]est dt
0 1
Z n+1 Z n+1
st
= [t]e dt = nest dt
n=1 n n=1 n
1 st t=n+1 1
=n e = n esn es(n+1)
n=1 s t=n s n=1
1 sn s 1 es
=
s n=1
ne (1 e ) =
s n=1 n(es)n
1 es es 1
= 2 = s
, si s > 0.
s 1 es s(e 1)
Por tanto,
1
L [t] = .
s(es 1)
152 Transformada de Laplace y Problemas de Cauchy
Ntese cmo la transformada de Laplace suaviza una funcin como parte entera.
L(k f1 + f2 ) = kL( f1 ) + L( f2 ).
En otras palabras, f es de orden exponencial si | f | se puede acotar por encima por una
funcin de la forma Met con M > 0 a partir de un nmero positivo T en adelante.
5.1 Transformada de Laplace: Definicin y Propiedades 153
1
Figura 5.1: Las funciones f (t) = sin(t) + 12 y g(t) = e 5 t . La funcin f (t) es de orden
exponencial 1/5 para todo t 2.
Ejemplo 5.7 Cualquier funcin de la forma f (t) = Met claramente es de orden exponen-
cial con t > 0.
Ejemplo 5.8 Toda funcin f acotada con f (t) M es de orden exponencial, tomando
= 0. En particular, las funciones sin(x) y cos(x) son de orden exponencial.
tn
lm = 0.
t et
n
Esto significa que dado > 0 existe N > 0 tal que si T > N entonces ett < , en particular
n
tomando = 1, existir N0 > 0 tal que ett < 1 si, y slo si, t n < et siempre que
t > N0 . Como consecuencia de esto, puede verse que toda funcin polinmica es de orden
exponencial.
Ejemplo 5.10 La funcin parte entera es de orden exponencial. Esto se debe al hecho que
[t] t + 1, para t > 0
pero la funcin del lado derecho es de orden exponencial (pues es un polinomio), luego la
funcin parte entera es de orden exponencial.
2
Ejemplo 5.11 La funcin definida por f (t) = et , no es de orden exponencial puesto que
para cualquier > 0, la funcin
t2
e
= et 2 t
et
no es acotada.
154 Transformada de Laplace y Problemas de Cauchy
La siguiente proposicin nos da un criterio sencillo para saber cundo una funcin F no
puede ser una transformada de Laplace.
Proposicin 5.3 Si lm F(s) 6= 0, entonces F no es la transformada de Laplace de ninguna
s
funcin de orden exponencial.
(n + 1) = n!.
Ejemplo 5.14
( + 1)
L(t ) = , para > 1.
s+1
Solucin. Por definicin tenemos
Z
L(t ) = t est dt.
0
En particular tenemos,
n!
L(t n ) = , n N.
sn+1
( 1
t 2 , si t > 0
h(t) = ,
0, si t = 0
1 ( 12 )
L(t 2 ) = 1 .
s2
Adems, dado que la funcin gamma satisface la siguiente propiedad de reflexin: ()(1
) = sin() concluimos que
1 1 2 1
1 = = = .
2 2 2 sin( 2 )
1
Por lo tanto, ( 21 ) = y as L(t 2 ) =
1 .
s2
s4
L(e4t cos(t)) =
(s 4)2 + 1
5.1 Transformada de Laplace: Definicin y Propiedades 157
Ejemplo 5.18 Hallar la transformada de Laplace de la funcin parte decimal definida por
f (t) = t [t], conocida comnmente como diente de sierra.
Solucin. Aplicando la linealidad de la transformada de Laplace tenemos:
L t [t] =L(t) L([t])
1 1
= 2 s .
s s(e 1)
Ejemplo 5.19 Con los resultados que aparecen en la Tabla 5.1 calcular L sin2 (at) .
Teorema 5.4 Si la transformada de laplace L( f (t)) = F(s) existe para s > , entonces
Ejemplo 5.20 Como sabemos la funcin definida como f (t) = 1 tiene transformada 1s .
1
Por tanto, la funcin g definida como g(t) = eat 1 tiene transformada G(s) = F(s a) = sa ,
como se haba visto antes.
Ejemplo 5.21 Como sabemos la funcin definida como f (t) = cos(bt) tiene transformada
s
F(s) = s2 +b 2 . Por tanto, la funcin g definida como g(t) = eat cos(bt) tiene transformada
sa
G(s) = F(s a) = (sa) 2 +b2 , como se haba visto antes.
158 Transformada de Laplace y Problemas de Cauchy
s
Ejemplo 5.22 Encontrar una funcin cuya transformada de Laplace sea G(s) = s2 +s+1
.
Solucin.
1 1
G(s) = = .
s2 + s + 1 1 2 3
s+ 2 +4
q
Vemos que corresponde a un caso particular del ejemplo anterior con a = 12 y b = 3
4 =
3 2 3 12 t 3
2 , luego una funcin puede ser g(t) = 3 e sin 2 t .
El siguiente teorema nos dice que la transformada de Laplace es inyectiva, lo que permite
hablar de la transformada inversa, y como veremos, tambin es lineal. Este hecho es clave
para resolver algunos problemas de ecuaciones diferenciales con valor inicial en el origen
como se haba mencionado al comienzo del captulo.
Teorema 5.5 Sean f y g funciones continuas sobre [0, ) de orden exponencial tales
que la transformada de Laplace existe para cada una y cuyas transformadas son iguales,
L( f ) = L(g). Entonces f (t) = g(t) para todo t 0.
Definicin 5.4 Se define la transformada inversa de Laplace como L1 (F(s)), de modo
que
L( f (t)) = F(s), entonces L1 (F(s)) = f (t).
Como es natural, mientras conozcamos transformadas de muchas funciones, ser ms fcil
encontrar la transformada inversa de alguna funcin.
Teorema 5.6 Suponga que L1 (F1 ) y L1 (F2 ) existen y son funciones continuas en [0, )
y sea k una constante, entonces
donde Q(s) y P(s) son polinomios en la variable s (con grado de P mayor que el de Q),
mediante el uso de fracciones parciales, para ello procederemos de la siguiente manera:
Z Descomponemos F(s) en fracciones simples, es decir:
Encontramos la races de P(s). Tomamos una de sus races reales con multi-
plicidad m, y a ib a una de las races complejas y su conjugada (en caso de
haber). De esta forma podemos expresar P(s) = (s )m [(s a)2 + b2 ].
Y as escribimos F(s) como:
7s2 + 23s + 30
F(s) = .
(s 2)(s2 + 2s + 5)
Solucin. Primero hallamos la descomposicin en fracciones parciales de F(s). Obser-
vamos que el factor cuadrtico s2 + 2s + 5 tiene una raz imaginaria, con su correspon-
diente conjugada, que es 1 2i de donde podemos expresar este factor cuadrtico como:
s2 + 2s + 5 = (s + 1)2 + 22 , por lo tanto su desarrollo en fracciones parciales queda de la
forma
7s2 + 23s + 30 A B(s + 1) + 2C
2
= + . (5.1)
(s 2)(s + 2s + 5) (s 2) (s + 1)2 + 4
Al multiplicar ambos lados por el comn denominador obtenemos
28 + 46 + 30 =A(4 + 4 + 5)
104 =A13
A =8.
7 23 + 30 =8(1 2 + 5) + 2C(3)
14 =32 6C
C =3.
1 =A(9 + 6 + 2)
1 =A17
1
A= .
17
1
Con s = 1 y A = 17 , sustituyendo en (5.4) tenemos
1
1= (1 2 + 2) +C(4)
17
5.1 Transformada de Laplace: Definicin y Propiedades 161
1
1= 4C
17
4
C = .
17
4 1
Finalmente, si s = 0, C = 17 yA= 17 sustituimos en (5.4) y obtenemos el valor
1 4
1= (2) + B(3) (3)
17 17
1
1 = 4C
17
2 12
1 = 3B +
17 17
1
B = .
17
1 4 1
Dado que B = 17 , C= 17 yA= 17 , sustituyendo en (5.3) se tiene la siguiente descom-
posicin
1 1/17 1/17(s + 1) 4/17
= +
(s 3)(s2 + 2s + 2) (s 3) (s + 1)2 + 1
con esto podemos hallar la transformada inversa de Laplace de F(s).
1 1
L =
(s 3)(s2 + 2s + 2)
1 1 1 1(s + 1) 1 4
=L +L +L
17(s 3) 17 (s + 1)2 + 1 17 (s + 1)2 + 4
1 1 1 1 1 (s + 1) 4 1 1
= L L L
17 (s 3) 17 (s + 1)2 + 1 17 (s + 1)2 + 1
1 1 4
= e3t et cos(t) et sin(t).
17 17 17
Obs El lector habr notado que nuestro inters siempre estuvo en fracciones propias. Esto
se justifica por lo siguiente. Si F no fuera una fraccin propia, lms F(s) 6= 0, luego
F no sera la transformada de Laplace de ninguna funcin de orden exponencial (que
son las de nuestro inters).
a t
1
Figura 5.2: Funcin escaln unitaria con desplazamiento a > 0.
= f (t)es(a+t) dt
0
as
as Z
e L f (t) =e f (t)est dt
0
as
Z
e L f (x a) = f (x a)esx dx
a
Z a Z
= 0 f (x a)esx dx + 1 f (x a)esx dx
Z0 a
sx
= U(x a) f (x a)e dx = L U(x a) f (x a) .
0
5.1 Transformada de Laplace: Definicin y Propiedades 163
Ejemplo 5.25 Expresar la funcin dada mediante funciones escaln unitario y calcular su
transformada.
1.
0, 0t <1
2, 1t <2
g(t) =
1, 2t <3
3, t 3.
g(t)
t
1 1 2 3 4
1
Solucin. 1. La grfica de la funcin g(t) es:
g(t)
t
1 1 2 3 4 5 6
1
g(t + a) = g(t 2 + 2) = t
5.1 Transformada de Laplace: Definicin y Propiedades 165
g(t + a) = g(t 3 + 3) = t
la transformada de Laplace de g es
1
L g(t + a) = L(t) = 2
s
y por la frmula (5.5) tenemos
1
L (t 3)U(t 3) = e3s 2 .
s
Finalmente, la transformada de Laplace de g es
L (t 1)U(t 1) 2L (t 2)U(t 2) + L (t 3)U(t 3) =
1s 2s 3s 1
(e 2e + e ) 2 .
s
se3s
s2 + 4s + 5
como el producto eas F(s). Para esto hacemos eas = e3s y
s
F(s) = .
s2 + 4s + 5
As, a = 3 y
1 s
f (t) = L = e2t cos(t) 2e2t sin(t).
s2 + 4s + 5
166 Transformada de Laplace y Problemas de Cauchy
e3s (s 5)
(s + 1)(s + 2)
f (t)
a 2a 3a 4a
t
1
5.1 Transformada de Laplace: Definicin y Propiedades 167
Algunas aplicaciones tales como circuitos elctricos y sistemas mecnicos, entre otras, estn
sometidas a unas fuerzas de gran magnitud que solamente acta durante un intervalo de
tiempo muy corto. Tales fuerzas son, por ejemplo, una tensin elctrica en el caso de un
circuito elctrico, una descarga elctrica, el golpe de un martillo, el golpe a una pelota de
tenis, etc. La funcin impulso unitario puede servir como un modelo para estas fuerzas.
Definicin 5.7 Impulso Unitario. La funcin impulso unitario : [0, ) R se
define por la frmula
0, si 0 t t0
1
(t t0 ) = 2 , si t0 < t t0 +
0, si t > t0 +
con > 0 y t0 > 0.
El rea bajo la funcin impulso unitario es:
Z
(t t0 )dt = 1.
0
1
2
Z
(t t0 )dt
0
t0 t0 t0 +
En la prctica es conveniente trabajar con otro tipo de impulso llamado funcin generalizada
de Dirac.
Definicin 5.8 Funcin generalizada delta de Dirac. La funcin generalizada delta
de Dirac est dada por
lm (t t0 ) = (t t0 ).
0
t0 t0 t0 +
0
De donde se tiene que la funcin generalizada delta de Dirac satisface las siguientes
propiedades:
(1) (
, si t = t0
(t t0 ) =
0, si t 6= t0 .
(2) Z
(t t0 )dt = 1.
Ms an, dado que (t t0 ) 6= 0 para t suficientemente cerca de t0 se cumple que
Z
(t t0 )dt = 1, para cada a < t0 .
a
Note que ninguna funcin satisface estas propiedades y por tanto (t t0 ) no es una funcin
en el sentido clsico. Sin embargo, a esta se le puede hallar su transformada de Laplace para
t0 > 0.
Z Z
st
L ( (t t0 )) = e (t t0 )dt = est lm (t t0 )dt.
0 0 0
Como (t t0 ) 6= 0 solo para t0 < t < t0 + , entonces est (t t0 ) 6= 0 solo si
t0 < t < t0 + . As,
Z
Z
st est (t t0 )dt, t0 < t < t0 +
e (t t0 )dt = 0 .
0 0, en otro caso
5.1 Transformada de Laplace: Definicin y Propiedades 169
Por lo tanto,
Z Z Z
st st0 st0
L ( (t t0 )) = lm e (t t0 )dt = e (t t0 )dt = e (t t0 )dt = est0 .
0 0 0 0
Ejercicios
1. Use la definicin para determinar la transformada de Laplace de la funcin dada.
a) f (t) = et sin(2t).
b) f (t) = te3t .
c) f (t) = [t]eat .
d)
(
0, 0 t < 2
f (t) =
t, t 2.
e)
(
sin(t), 0 t <
f (t) =
0, t .
f)
(
e2t , 0 t < 3
f (t) =
1, t 3.
f (t)
t
1 1 2 3
1
b)
(
0, 0 t 3/2
f (t) =
sin(t), t > 3/2.
c)
1, 0 t < 4
f (t) = 0, 4 t < 5
1, t 5.
5.1 Transformada de Laplace: Definicin y Propiedades 171
d) f con grfica
a b
1
e) f con grfica
1 1 2 3
1
4. Encuentre la transformada de Laplace.
a) (t 1)U(t 1).
b) e2t U(t 2).
c) (3t + 1)U(t 1).
d) (t t0 ).
e) f (t) (t t0 ), f una funcin dada.
5. Determinar L( f (t))(s), donde la funcin peridica queda descrita mediante su grfica.
a 2a 3a 4a
1
6. Demostrar las frmulas (5.6), (5.7) y (5.8).
7. Demostrar la frmula (5.9)
8. Demostrar el Teorema 5.1.
9. Demostrar el Teorema 5.2.
10. Demostrar la Proposicin 5.3.
Teorema 5.8 Sea f (t) continua por partes en [0, ) y de orden exponencial . Entonces
para s > ,
dnF
L(t n f (t)) = (1)n (s).
dsn
5.1 Transformada de Laplace: Definicin y Propiedades 173
Ejemplo 5.29 Calcular la transformada de la funcin definida por g(t) = t 2 [t] donde [t]
denota la parte entera de t.
n
Ejemplo 5.30 Verificar que si L(t n f (t)) = (1)n ddsFn (s), su transformada inversa es
1 dnF
L (s) = (1)nt n f (t) (5.10)
dsn
donde f (t) = L1 (F). Usar esta ecuacin para calcular la transformada inversa de
s + 2
F(s) = ln .
s 5
Solucin. Utilizando la ecuacin (5.10), primero tenemos que derivar F(s) = ln s+2
s5 , por
lo tanto tenemos:
dF 7
= (5.11)
ds (s + 2)(s 5)
como derivamos una vez, n = 1 en la ecuacin (5.10), as, sustituyendo en la ecuacin (5.10)
obtenemos
1 dF
L =(t) f (t)
ds
7
L1 =(t) f (t).
(s + 2)(s 5)
7 A B
= + . (5.12)
(s + 2)(s 5) (s + 2) (s 5)
La convolucin de dos funciones se puede ver como una operacin algebraica entre ellas la
cual satisface las siguientes propiedades:
1. f g = g f (conmutativa).
2. f (g + h) = ( f g) + ( f h) (distributiva respecto a la suma).
3. ( f g) h = f (g h) (distributiva respecto a la convolucin).
4. f 0 = 0 (existencia del anulador para la convolucin).
La trasformada de Laplace convierte la convolucin de funciones en producto como se
muestra en el siguiente resultado.
Teorema 5.9 Teorema de Convolucin. Sean f y g funciones continuas por partes
en [0, ) y de orden exponencial . Sea F(s) = L( f (t))(s) y G(s) = L(g(t))(s). Entonces,
L( f g) = F(s)G(s)
o equivalentemente,
L1 (F(s)G(s)) = ( f g)(t).
Ejercicios
dnF
1. Sea L1 dsn (t) = (t)n f (t) donde f (t) = L1 (F). Use esta ecuacin para calcular
L1 (F).
176 Transformada de Laplace y Problemas de Cauchy
a) F(s) = ln s+2 .
s5
b) F(s) = ln s3 .
s4
2
c) F(s) = ln ss2+29+1
.
d) F(s) = arctan 1s .
n
2. Use L t n f (t) = (1)n ddsFn (s) para demostrar que
a) L t 2 y0 (t) (s) = sY 00 (s) + 2Y 0 (s), donde Y (s) = L (y(t)) (s).
b) L t 2 y0 (t) (s) = s2Y 00 (s) + 24sY 0 (s) + 2Y (s), donde Y (s) = L (y(t)) (s).
3. Usando la convolucin hallar la transformada inversa de Laplace de
s
a) F(s) = (s2 +1)(s2) .
1
b) F(s) = (s1)(s 2 +1) .
40 + 156 + 28 32 =D5(3)
D =8,
178 Transformada de Laplace y Problemas de Cauchy
32 =10B + 8 40
B =0.
52 =4A + 8 64
A =1.
Ejemplo 5.34 Resolver el problema con valor inicial mediante la transformada de Laplace
g(t) =t 4, a=2
5.2 Problemas de Cauchy con Condicin Inicial en el Origen 179
g(t + 2) =t 2
1 2
L (g(t + 2)) =L(t 2) = L(t) 2L(1) = ,
s2 s
por lo tanto, L((t 4)U(t 2)) = e2s s12 2s .
En la ecuacin (5.18) reemplazamos las condiciones iniciales y L[(t 4)U(t 2)],
2 1 2s 1 2
s L(y) 1 + L(y) = 2 e
s s2 s
2 1 2s 1 2
L(y)(s + 1) = 2 + 1 e
s s2 s
1 1 2s 1 2
y(t) = L e
s2 s2 s
y(t) = t + (4 t 2 cos(t 2) + sin(t 2))U(t 2).
e s
Para resolver el trmino 1 utilizamos la ecuacin (5.6), primero expresamos s2 +2s+1
como el producto eas F(s). Para esto hacemos eas = es y F(s) = (s+1) 1
2 . As
a = y f (t) = L1 (s+1)
1
2 = et t. Por lo tanto,
1 es
L = f (t )U(t )
s2 + 2s + 1
=U(t ) e(t) (t ) .
5.2 Problemas de Cauchy con Condicin Inicial en el Origen 181
Ahora aplicamos transformada de Laplace a ambos lados, reemplazamos con cada una de
las frmulas obteniendo as:
Ejercicios
1. Use la transformada de Laplace para resolver el problema de valor inicial respectivo.
a) y0 + 6y = e4t , y(0) = 2.
b) y00 + 2y0 + y = 6 sin(t) 4 cos(t), y(0) = 1, y0 (0) = 1.
c) y00 + y = sin(2t), y(0) = 0, y0 (0) = 1.
d) y00 + 3y0 + 2y = et , y(0) = 0, y0 (0) = 1.
e) y00 + 4y = 8 sin(2t) + 9 cos(t), y(0) = 1, y0 (0) = 0.
182 Transformada de Laplace y Problemas de Cauchy
2. Hallar la solucin
Z de las ecuaciones integro-diferenciales
t
a) y0 y + y(v)etv (t v)dv = et , y(0) = 1
0 Z t+1
e2 2
b) y00 + 4 y +U(t + 1) y0 (v) (t v 1)e(tv) dv = (t 2), y(0) = y0 (0) = 0.
Z t 0
c) y0 + 4 y(v)dv = t sin(t), y(0) = 2.
Z t0
d) y0 = y(v) cos(t v)dv, y(0) = 1.
0
3. Sea L (t f (t)) = dF
ds donde F(s) = L( f ). Use esta ecuacin para resolver el problema
de valor inicial
) y(t) = 14 14 cos(2t) 1 1 1
2 sin(2t) 4 (t 1)U(t 1) + 4 cos(2t 2)U(t 1).
3
2 e 2 t sin 3
2. a) y(t) = 3 2 t .
e
b) y(t) = U(t 2)(t 2)e 2 (t2) .
d) y(t) = 14 13 cos(t) + 12
25
cos(2t).
Captulo 6
n
n=0 cn (x a) es convergente
Definicin 6.2 Convergencia. Una serie de potencias
en un valor determinado si su sucesin de sumas parciales converge; es decir, si existe
N
lm SN (x) = lm cn (x a)n .
N N
n=0
c0 + c1 + c2 + + cn +
converge absolutamente si
cn+1
cn = L < 1.
En este caso el radio de convergencia es dado por R = 1/L. Si R 6= 0 o R 6= , el intervalo
de convergencia puede o no incluir los extremos a R y a + R.
Teorema 6.1 Si una serie de potencias converge en un intervalo |x a| < R con R > 0,
entonces la serie define una funcin f (x) continua en todo el intervalo |x a| < R.
El teorema anterior establece que si una serie de potencia converge entonces esta define una
funcin continua, sin embargo conocer tal funcin no es fcil, de hecho muchas series de
potencias convergentes no definen funciones elementales.
n n
Sean f (x) = n=0 an (x a) y g(x) = n=0 bn (x a) dos funciones defi-
Teorema 6.2
nida mediante series de potencias convergente en un intervalo |x a| < R. Entonces
1. f (x) = g(x) si, slo si, an = bn , para todo n = 0, . . . .
2. f (x) g(x) = n n n
n=0 an (x a) n=0 bn (x a) = n=0 (an bn )(x a) .
n n n
3. f (x)g(x) = n=0 an (x a) n=0 bn (x a) = n=0 cn (x a) , donde cn = a0 bn +
a1 bn1 + + an b0 .
f (x) n
4. Si g(x) 6= 0 en |x a| < R, entonces g(x) n=0 dn (x a) , donde en la mayora
=
de los casos los coeficientes se pueden obtener con mayor facilidad al igualar los
coeficientes en la relacin equivalente:
!
n
n n n
an(x a) = dn(x a) bn(x a) = dk bnk (x a)n .
n=0 n=0 n=0 n=0 k=0
El siguiente resultado muestra que una funcin representada por una serie de potencias tiene
derivada de todos los ordenes en el intervalo de convergencia de la serie. Estas derivadas se
calculan derivando trmino a trmino.
Teorema 6.3 Una serie de potencias
f (x) = an(x a)n
n=0
con radio de convergencia R > 0 tiene derivada de todos los ordenes en su intervalo de
convergencia y
0
f (x) = an n(x a)n1
n=1
f 00 (x) = an n(n 1)(x a)n2
n=2
..
.
f (k) (x) = an n(n 1)(n 2) (n k + 1)(x a)nk .
n=k
Adems todas estas series tienen el mismo radio de convergencia R.
Obs Note que la serie de la j-sima derivada de una funcin f (x) expresada en serie de
potencias puede comenzar en el ndice de sumatoria n = 0, o cualquier valor de n < j
puesto que los correspondientes coeficientes son cero y por tanto no aportan valor a la
sumatoria.
188 Soluciones de EDO Mediante Series
escriba la funcin x f 00 (x) f (x) de modo que el trmino general sea un mltiplo constante
de (x a)n .
Luego,
x f 00 (x) f (x) = x an n(n 1)(x a)n2 an (x a)n
n=2 n=0
=(x a + a) an n(n 1)(x a)n2 an (x a)n
n=2 n=0
=(x a) an n(n 1)(x a)n2 + a an n(n 1)(x a)n2 an (x a)n
n=2 n=2 n=0
= an n(n 1)(x a)n1 + an n(n 1)a(x a)n2 an (x a)n .
n=2 n=2 n=0
entonces
x f 00 (x) f (x) = an+1 (n + 1)n(x a)n + an+2 (n + 2)(n + 1)a(x a)n
n=0 n=0
an (x a)n
n=0
= [an+2 (n + 2)(n + 1)a + an+1 (n + 1)n an ](x a)n .
n=0
6.1 Series de Potencias 189
Balanceando los ndices de las sumatorias, podemos reescribir las series anteriores como
n1 n
ann(n 1)(x + 2) = an+1 (n + 1)n(x + 2) = an+1(n + 1)n(x + 2)n
n=2 n=1 n=0
ann(n 1)(x + 2)n2 = an+2(n + 2)(n + 1)(x + 2)n.
n=2 n=0
De forma similar,
3(x + 2) an n(x + 2)n1 = 3 an n(x + 2)n = 3 an n(x + 2)n .
n=1 n=1 n=0
Concluimos entonces,
(x 1) f 00 (x) 3(x + 2) f 0 (x) + 3 f (x)
= an+1 (n + 1)n(x + 2) 3 an+2 (n + 2)(n + 1)(x + 2) 3 an n(x + 2)n
n n
n=0 n=0 n=0
+ 3 an (x + 2)n
n=0
= (an+1 n(n + 1) 3an+2 (n + 2)(n + 1) 3nan + 3an ) (x + 2)n .
n=0
190 Soluciones de EDO Mediante Series
f 00 (x0 ) f (n)
f (x) = f (x0 ) + f 0 (x0 )(x x0 ) + (x x0 )2 + + (x x0 )n + , |x x0 | < R
2! n!
es llamada una expansin de Taylor de f (x) en potencias de (x x0 ). Si x0 = 0, la serie es
llamada expansin en serie de Maclaurin de f (x).
Teorema 6.4 Taylor. Si una funcin f (x) tiene derivada de todos los ordenes en un
intervalo |x x0 | < R, entonces
f 00 (x0 ) f (n)
f (x) = f (x0 ) + f 0 (x0 )(x x0 ) + (x x0 )2 + + (x x0n ) + , |x x0 | < R,
2! n!
es decir, la serie infinita converge a f (x).
Definicin 6.6 Funcin Analtica. Una funcin f (x) es analtica en un punto x0 si se
puede expresar como una expansin de Taylor en potencias de x x0 para cualquier x en
una vecindad de x0 .
Veamos la expansin en series de Maclaurin de algunas funciones elementales:
xn
ex = , converge en < x < .
n=0 n!
x2n+1
sin(x) = (1)n , converge en < x < .
n=0 (2n + 1)!
x2n
cos(x) = (1)n , converge en < x < .
n=0 (2n)!
x2n+1
sinh(x) = , converge en < x < .
n=0 (2n + 1)!
x2n
cosh(x) = , converge en < x < .
n=0 (2n)!
1
= xn , converge en 1 < x < 1.
1 x n=0
6.1 Series de Potencias 191
xn
ln(1 + x) = (1)n+1 , converge en 1 < x 1.
n=1 n
x2n+1
tan1 (x) = (1)n , converge en 1 x 1.
n=0 2n + 1
1 3 5 (2n + 1) x2n+1
sin1 (x) = , converge en 1 x 1.
n=0 2 4 6 2n 2n + 1
r(r 1)x2 r(r 1)(r 2)x3
(1 + x)r =1 + rx + + + , converge en 1 < x < 1.
2! 3!
Ejercicios
1. Determinar si cada una de las siguientes series es convergente o divergente.
n
a) n=1 (4n+3)(4n1) .
n+1
b) n=1 2n .
2+(1)n
c) n=1 2n .
ln(n)
d) n=2 nn+1 .
2
n cos (n/3)
e) n=1 2n .
1/n
f ) n=1 (n 1) .
n
n3 ( 2+(1)n )n
g) n=1 3n .
2. Para cada una de las siguientes series hallar el radio de convergencia R. Si R > 0,
hallar el intervalo de convergencia.
(1)n n
a) n=0 2n n (x 1) .
n
b) n7 n
n=0 (1) n! x .
n n
c) n=0 2 n(xn 2) .
3
d) n=0 4n+1 (n+1)2 (x + 7)n .
3. Para y(x) = n
n=0 an x , hallar la serie de potencias resultante
00
a) (x 3)y + xy 2y. 0
ai (x) g(x)
donde Pi (x) = an (x) , i = 0, . . . , n 1 y Q(x) = an (x) .
Definicin 6.7 Se dice que un punto x0 es punto ordinario de la ecuacin (6.1) si Pi (x)
i = 0, . . . n 1 y Q(x) dados en (6.2), son funciones analticas en x0 . Se dice que un punto
que no es ordinario es punto singular de la ecuacin.
y00 xy = 0 (6.5)
Para sumar dos series, es necesario que ambos ndices de suma comiencen en el
mismo nmero y que las potencias de estas series comiencen con la misma potencia,
para lograr esto hacemos lo siguiente: Para la primera serie hacemos el cambio en
el ndice k = n 2, para la segunda k = n y para la tercera tomamos k = n y as
balanceamos los exponentes de la variable x. Por lo tanto, el lado derecho se convierte
en:
(k + 2)(k + 1)ck+2xk 2 ck kxk + ck xk .
k=0 k=1
| {z } | {z } k=0
| {z }
k=n2 k=n k=n
Ahora queremos balancear los ndices de la sumas, como la primera serie comienza
en 0 y la segunda en 1. Si se escribe el primer trmino de la primera serie fuera de la
notacin sigma, y lo mismo hacemos con la tercera serie obtenemos:
2 1 c2 + (k + 2)(k + 1)ck+2 xk 2 ck kxk + c0 + ck xk = 0
k=1 k=1 k=1
como la ecuacin (6.6) es igual a cero es necesario que el coeficiente de cada potencia
de x sea igual a cero, esto es, que 2c2 + c0 = 0 (es el coeficiente de x0 ), y tenemos:
c0
c2 =
2
c1
para k = 1, c3 = 3!
3c2
para k = 2, c4 = 43 , pero como conocemos c2 sustituimos y obtenemos
3c0
c4 =
4!
5c3
para k = 3, c5 = 45 , pero como conocemos c3 sustituimos y obtenemos
c1
c5 =
5!
7c4
para k = 4, c6 = 56 , pero como conocemos c4 sustituimos y obtenemos
21c0
c6 =
6!
9c5
para k = 5, c7 = 67 , pero como conocemos c5 sustituimos y obtenemos
45c1
c7 = .
7!
Y as sucesivamente. Por tanto,
c0 c1 3c0 4 5c1 5 21c0 6 45c1 7
y =c0 + c1 x x2 + x3 x + x x + x +
2! 3! 4! 5! 6! 7!
x2 3x4 21x6 x3 5x5 45x7
=c0 1 + c1 x + + + + .
2! 4! 6! 3! 5! 7!
y00 + sin(x)y
n2
x2n+1
n
n
= cnn(n 1)x + sin(x) = (1) (2n + 1)! cnx =0
n=2 n=0 n=0
= c2 2 + 6c3 x + 12c4 x2 + 20c5 x3 + 30c6 x4 + +
x3 x5 x7
+ x + + c0 + c1 x + c2 x2 + = 0
3! 5! 7!
6.2 Soluciones en Series de Potencias 195
c0
= 2c2 + x(6c3 + c0 ) + x2 (12c4 + c1 ) + x3 20c5 + c2 +
3!
c1
+ x4 30c6 + c3 + = 0.
3!
Por consiguiente,
2c2 = 0, 6c3 + c0 = 0, 20c5 + c2 + c3!0 = 0, 30c6 + c3 + c3!1 = 0 y as sucesiva-
mente. Esto da como resultado c2 = 0, c3 = c60 , c4 = 12
c1 c0
, c5 = 120 , . Al
agrupar trminos llegamos a la solucin general
c0 c1 c0 5
y =c0 + c1 x x3 x4 + x +
6 12 120
x3 x5 x4
=c0 1 + + + c1 x + .
6 120 12
Ejemplo 6.7 Hallar la solucin en serie de potencias del problema de valor inicial
Solucin. Dado que las funciones x 1, x2 y x son polinomios, entonces son funciones
analticas en todo punto. Luego del Teorema de Existencia y Unicidad, podemos hallar una
solucin en series de potencias del problema de valor inicial dado. Como conocemos el
valor de la eventual solucin y su derivada en 0, entonces busquemos una solucin de la
forma
y00 (0) 2 y000 (0) 3
y(x) = y(0) + y0 (0)x + x + x + . (6.8)
2! 3!
Reemplazando las condiciones iniciales, x = 0, y(0) = 1, y0 (0) = 1 en (6.7) obtenemos
De esta forma tenemos los coeficientes y(0), y0 (0) y y00 (0) de la solucin (6.8). Para hallar
los otros coeficientes derivamos (6.7) y usamos las condiciones iniciales. Esto es:
evaluando en y(0) = y0 (0) = y00 (0) = 1, obtenemos y000 (0) = 3 y y(4) (0) = 3. Sustituyendo
en la solucin (6.8) tenemos
x2 x3 x4
y(x) = 1 + x + + + + , < x < .
2 2 8
196 Soluciones de EDO Mediante Series
y00 + y = 0.
ahora balanceamos los ndices de las series, hacemos k = n 2 para la primera serie y k = n
para la segunda. De donde obtenemos
ak+2(k + 1)(k + 2)(x x0)k + ak (x x0)k =0,
k=0 k=0
[ak+2(k + 1)(k + 2) + ak ] (x x0)k =0.
k=0
ak+2 (k + 1)(k + 2) + ak = 0,
Este tipo de ecuaciones aparece en diferentes aplicaciones con x0 = 0, e incluye como casos
particulares a la ecuacin de Legrendre (6.3), la ecuacin de Airy (6.5), la ecuacin de
Chebyshev:
y la ecuacin de Hermite:
y definamos
s
b j = 1, s < r, no importa los b j .
j=r
Note que la ecuacin (6.12) es un caso particular de la ecuacin (6.9). El siguiente teorema
muestra las recurrencias para el caso general de la ecuacin (6.9).
Los coeficientes {an } de cualquier solucin y(x) = n
Teorema 6.6 n=0 an (x x0 ) de la
ecuacin diferencial
(1 + (x x0 )2 )y00 + (x x0 )y0 + y = 0
Ejemplo 6.10 Determinar la solucin en serie de potencias del problema de valor inicial
1 2
x + 2x + 3 y00 + 2(x + 1)y0 y = 0, y(1) = 1, y0 (1) = 1. (6.14)
2
Solucin. Completando cuadrados reescribimos
1 2 1 2
1 2
x + 2x + 3 = (x + 1) + 2 = 1 + (x + 1)
2 2 2
por tanto, la ecuacin (6.14) queda como
1
1 + (x + 1) y00 + 2(x + 1)y0 y = 0
2
2
que es de la forma general (6.9) con x0 = 1, = 1/2, = 2 y = 1. Del teorema
anterior podemos calcular los coeficientes de la solucin mediante el polinomio
1
p(n) = n(n 1) + 2n 1.
2
Los coeficientes se calculan de la siguiente manera:
1
2 (2m)(2m 1) + 4m 1 2m2 + 2m 1
a2m+2 = a2m = a2m
(2m + 2)(2m + 1) (2m + 2)(2m + 1)
6.2 Soluciones en Series de Potencias 201
1
2 (2m + 1)(2m + 1 1) + 2(2m + 1) 1
a2m+3 = a2m+1
(2m + 2)(2m + 3)
2m2 + 5m + 1
a2m+1 .
(2m + 2)(2m + 3)
1 1
a2 = a0 =
12 12
2+21 3 1 13
a4 = a2 = =
34 43 12 1234
24+41 11 13 1 3 11
a6 = a4 = =
56 56 1234 6!
18 + 6 1 23 1 3 11 1 3 11 23
a8 = a6 = =
78 78 6! 8!
..
.
m
m j=1
(2( j 1)( j + 2) + 3)
a2m =(1) .
(2m)!
1 1
a3 = (1) =
23 23
2+5+1 8 1 18
a5 = a3 = =
45 45 23 5!
24+52+1 19 18 1 8 19
a7 = a5 = =
67 67 5! 7!
2
3 2 + 15 + 1 34 1 8 19 1 8 19 34
a9 = a7 = =
89 89 7! 9!
..
.
mj=1 (2 j2 + 5 j + 1)
a2m+1 =(1)m .
(2m + 1)!
Ejercicios
1. Determinar los puntos singulares de la ecuacin diferencial dada.
a) (x + 1)y00 x2 y0 + 3y = 0
b) (x2 + x)y00 + 3y0 6xy = 0
c) sin(x)y00 + cos(x)y = 0
2. Determinar al menos los primeros cuatro trminos no nulos en un desarrollo en serie
de potencias en torno de x = 0 para una solucin general de la ecuacin diferencial
dada.
a) y00 xy = 0 (Ecuacin de Airy)
b) y00 + (x 1)y0 + y = 0
c) y00 2y0 + y = 0
d) (2x 3)y00 xy0 + y = 0
3. Determinar un desarrollo en series de potencias en torno de x = 0 para una solucin
general de la ecuacin diferencial dada. Su respuesta debe incluir una frmula general
para los coeficientes.
a) y0 2xy = 0
b) y00 xy0 + 4y = 0
c) y00 x2 y0 xy = 0
d) (x2 + 1)y00 xy0 + y = 0
4. Determinar al menos los primeros cuatro trminos no nulos en un desarrollo en serie
de potencias en torno de x = 0 para una solucin general de la ecuacin diferencial
dada.
a) y0 + sin(x)y = 0; y(0) = 1
b) y0 ex y = 0; y(0) = 1
c) (x2 + 1)y00 ex y0 + y = 0; y(0) = 1, y0 (0) = 1
d) y0 xy = sin(x)
e) y00 2xy0 + 3y = x2
f ) (1 + x2 )y00 xy0 + y = ex
g) y0 + sin(x)y = cos(x)
5. Hallar la serie de potencias en x para la solucin general de:
a) (1 + x2 )y00 + 6xy0 + 6y = 0.
b) (1 + x2 )y00 9xy0 + 16y = 0.
c) (1 + x2 )y00 + 7xy0 3y = 0.
d) (1 + x2 )y00 + 3xy0 2y = 0.
6.3 Soluciones en Series Infinitas en Puntos Singulares Regulares 203
m 1 m1 2 2m+1 .
m=0 (1) (2m+1)! j=0 ((2 j) 3)x
a1
n1 2 1
y(x) = a0 n k=0 8k 12k+ 2 (x 3)2n +
6. a) n=0 (1) (2n)!
nk=0 8k2 4k+ 29
a1
n=0 (2n+1)! (x 3)2n+1 .
2 5
m
k=0 k + 6 k+2
c) y(x) = a0 n=0 (2n+2)! (x 1)2n+2 +
m 2 11 5
a1 n+1 k=0 k + 6 k+ 3 (x 1)2n+3 .
n=0 (1) (2n+3)!
Definicin 6.8 Se dice que un punto singular x0 es punto singular regular de la ecuacin
(6.15) si p(x) = (x x0 )P(x) y q(x) = (x x0 )2 Q(x) son funciones analticas en x0 . Se dice
que un punto que no es singular regular es punto singular irregular de la ecuacin.
1 x2
p(x) = , q(x) = .
x1 x1
Ambas funciones son analticas en x0 = 0,
1
= (1 + x + x2 + )
x1
x2
=x2 + x3 + x4 + , |x| < 1.
x1
En el caso x0 = 1 se tiene
1
p(x) = , q(x) = (x 1)
x
las cuales son analticas en x0 = 1:
1
=1 (x 1) + (x 1)2 (x 1)3 +
x
(x 1) =1 x + 0 + .
2. Reescribamos la ecuacin,
1 1
y00 + y0 y = 0.
x2 (x 1)2 (x 1)2
De donde tenemos
1 1
p(x) = (x x0 ) , q(x) = (x x0 )2 .
x2 (x 1)2 (x 1)2
Para x0 = 0
1 x2
p(x) = , q(x) = .
x(x 1)2 (x 1)2
1
p(x) = , q(x) = 1
x2 (x 1)
Queremos hallar una solucin en serie de la ecuacin (6.15) en torno a un punto singular
regular x0 , el mtodo que mostraremos es vlido en una vecindad del punto x0 y la serie que
se obtiene, llamada serie de Frobenius, viene dada en la forma
Notemos que si x0 es un punto singular regular de (6.15), entonces P(x) tiene un factor
(x x0 ) en su denominador y/o Q(x) tiene a (x x0 )2 en su denominador, en cualquier caso
multiplicando (6.15) por (x x0 )2 tenemos que la ecuacin se transforma en
en la cual p(x) y q(x) son analticas en x0 . Veamos ahora el teorema que garantiza la
existencia de soluciones en forma de serie de Frobenius de una ecuacin diferencial lineal
con coeficientes variables.
206 Soluciones de EDO Mediante Series
Teorema 6.7 Frobenius. Sea x0 un punto singular regular de (6.17). Entonces (6.17)
tiene al menos una solucin en forma de serie de Frobenius de la forma (6.16). Esto es
vlido en el intervalo comn de convergencia de p(x) y q(x) de (6.17), excepto quiz en
x = x0 ; es decir, si cada expansin en serie de Taylor de p(x) y q(x) es vlida en el intervalo
|x x0 | < r, entonces existe al menos una solucin en serie de Frobenius en |x x0 | < r,
excepto quiz en x = x0 .
Note que el teorema anterior no garantiza la existencia de dos soluciones linealmente inde-
pendientes de la forma (6.16) para la ecuacin diferencial (6.17). Para hallar las soluciones
de (6.16), sea y(x) una solucin en serie de Frobenius de la ecuacin (6.17). Haciendo una
traslacin de eje podemos considerar x0 = 0. De esta manera la ecuacin (6.17) queda en la
forma
p(x) =b0 + b1 x + b2 x2 +
q(x) =c0 + c1 x + c2 x2 + .
La funcin y(x) es una solucin de la ecuacin (6.18) si la trasforma en una igualdad. As,
evaluando en (6.18), se debe satisfacer:
x2 s(s 1)a0 xs2 + (s + 1)sa1 xs1 + (s + 2)(s + 1)a2 x2 +
+ x b0 + b1 x + b2 x2 + sa0 xs1 + (s + 1)a1 xs + (s + 2)a2 xs+1 +
+ c0 + c1 x + c2 x2 + a0 xs + a1 xs+1 + a2 xs+2 + = 0.
Del lgebra de las series (Teorema 6.2), agrupando trminos semejantes obtenemos,
s(s 1) + b0 s + c0 = 0. (6.20)
Esta ecuacin cuadrtica es llamada ecuacin indicial y sus dos soluciones s1 y s2 satisfacen
uno de los siguientes casos.
Caso I: Las soluciones son diferentes y su diferencia no es un nmero entero: s1 6= s2 y
s1 s2
/Z
Caso II: Las soluciones son diferentes y su diferencia es un nmero entero: s1 6= s2 y
s1 s2 Z,
Caso III: Las soluciones son iguales: s1 = s2 .
Ahora veremos la manera de hallar las soluciones en serie de Frobenius de la ecuacin
diferencial (6.18) para cada uno de estos casos.
Ejercicios
1. Hallar los puntos singulares regulares de las siguientes ecuaciones diferenciales
a) (x2 1)y00 3(x + 1)y0 + (x2 + x)y = 0.
b) (x 1)3 x2 y00 + 8(x 1)xy0 + 2y = 0.
c) (x + 3)2 x4 y00 3(x + 3)xy0 + y = 0.
d) (x + 1)2 y00 + xy0 (x 1)y = 0.
2. Defina punto singular regular para una ecuacin diferencial lineal de tercer orden
puesto que los coeficientes bk y ck (k = 0, . . . ) de las funciones analticas p(x) y q(x) res-
pectivamente, son dados. Para este s1 y los valores a1 , a2 , a3 , . . . , cada coeficiente de (6.19)
es cero. Es decir, para este valor s = s1 y los correspondientes coeficientes a1 , a2 , a3 , . . . , la
funcin y(x) es una solucin de la ecuacin diferencial (6.18).
La segunda solucin se obtiene de manera similar sustituyendo s2 en (6.19). Sin embargo,
como veremos en las Secciones 6.3.2 y 6.3.3, no toda ecuacin de la forma (6.18) tiene dos
soluciones en serie linealmente independiente. En caso de existir dos soluciones linealmente
independientes, el intervalo de existencia de las dos soluciones viene dado en el siguiente
resultado.
Teorema 6.8 Las dos soluciones en serie de Frobenius de la ecuacin diferencial (6.18)
son linealmente independiente. Cada solucin existe en el intervalo comn de convergencia
de p(x) y q(x) excepto quizs en x = 0.
x ex
y00 + y0
+ y = 0.
x2 (1 + x2 ) x2
En este caso
x 1 ex
p(x) = x = , q(x) = x2 = ex .
x (1 + x ) 1 + x2
2 2 x2
Tenemos que x = 0 es el punto singular regular. Usaremos (6.19) para hallar la solucin en
serie de Frobenius. Para esto necesitaremos los coeficientes de p(x) y q(x)
1 1
p(x) = + x, b0 = , b1 = 1, b2 = b3 = = 0,
2 2
1 1
q(x) = x2 , c0 = , c2 = 1, c1 = c3 = c4 = = 0.
2 2
La ecuacin indicial (6.20) queda como
1 1
s(s 1) + s = 0
2 2
cuyas races son s1 = 1 y s2 = 1/2. Como s1 s2
/ Z, podemos sustituir s1 = 1 en (6.19).
s+n
Hagamos cero al coeficiente que acompaa a x , teniendo en cuenta que
b2 , b3 , . . . , c1 , c3 , c4 , = 0.
Esto es
1 1
an (1 + n)n + (1 n) + nan1 + an2 (0 1) + a0 (0) = 0.
2 2
4 9 9
7a2 = 2a1 + a0 = a0 + a0 = a0 , a2 = a0 ,
5 5 35
27 27 2 82
a3 = 3a2 + a1 = a0 a0 , a3 = a0 .
2 35 35 945
As tenemos que
2 9 2 82 3
y1 (x) = a0 x 1 x + x x +
5 35 945
210 Soluciones de EDO Mediante Series
son los cuatro primeros trminos de una solucin en serie de Frobenius de la ecuacin
diferencial dada.
Para hallar la segunda solucin en serie, evaluamos s2 = 1/2 como antes. En este caso el
trmino que acompaa a xs+n queda en la forma
1 3 1 1 1 3
an n n + n + an1 n + an2 (0 1) = 0
2 2 2 2 2 2
de la cual se obtiene la relacin de recurrencia
2n2 3n 3
an = n an1 + an2 , n 2.
2 2
3xy00 + (1 + x)y0 + y = 0.
Como s1 s2
/ Z, podemos sustituir s1 = 0 en (6.19) para hacer cero el coeficiente que
s+n
acompaa a x , o podemos derivar y sustituir en la ecuacin. Hagamos lo segundo.
y(x) = an x n , y0 (x) = annxn1, y00 (x) = ann(n 1)xn2.
n=0 n=1 n=2
3x an n(n 1)xn2 + (1 + x) an nxn1 + an xn = 0
n=2 n=1 n=0
n1 n1
3ann(n 1)x + an nx + an nx + an xn = 0
n
n=2 n=1 n=1 n=0
3an+1n(n + 1)xn + an+1(n + 1)xn + annxn + anxn = 0
n=1 n=0 n=0 n=0
(3an+1(n + 1)n + an+1(n + 1) + ann + an)xn = 0.
n=0
De donde tenemos que la ecuacin anterior se convierte en una identidad para todo x si se
satisface la recurrencia
Esto es,
Ejemplo 6.15 Hallar los 6 primeros trminos de las dos soluciones en serie de Frobenius
linealmente independientes de la ecuacin diferencial
2x2 (x2 + x + 1)y00 + x(11x2 + 11x + 9)y0 + (7x2 + 10x + 6)y = 0.
Solucin. Definamos los polinomios como en (6.22).
p0 (r) =2r(r 1) + 9r + 6 = (2r + 3)(r + 2)
p1 (r) =2r(r 1) + 11r + 10 = (2r + 5)(r + 2)
p2 (r) =2r(r 1) + 11r + 7 = (2r + 7)(r + 1).
Las races de la ecuacin indicial p0 (r) = 0 son r1 = 3/2 y r2 = 2. Note que r1 r2 =
1
2
/ Z. Por lo tanto, las soluciones de la ecuacin diferencial son
y1 (x) =x3/2 an(3/2)x2
n=0
y2 (x) =x2 an(2)xn.
n=0
a1 (2) =0
..
.
(n 1)an1 (2) + (n 3)an2 (2)
an (2) = , n 2.
n
Calculando el resto de trminos usando esta recurrencia, las dos soluciones de la ecuacin
son:
3/2 1 2 2 5 3 7 4 76 5
y1 (x) =x 1 x+ x x + x + x +
3 5 21 135 1155
2 1 2 1 3 1 4 1 5
y2 (x) =x 1+ x x + x + x + .
2 3 8 30
a0 (r) =1
..
.
p1 (n + r 1)
an (r) = an1 (r), n 1.
p0 (n + r)
Por lo que podemos hallar explcitamente los coeficientes an (r) dependiendo solamente
de an1 (r). Ms an, se tiene que
n
p1 ( j + r 1)
an (r) = (1)n .
j=1 p0 ( j + r)
a0 (r) =1
a1 (r) =0
..
.
p2 (n + r 2)
an (r) = an2 (r), n 2.
p0 (n + r)
216 Soluciones de EDO Mediante Series
donde,
a0 (r) =1
p2 (2n + r 2)
a2n (r) = a2n2 , n 1
p0 (2n + r)
n
p2 (2 j + r 2)
=(1)n .
j=1 p0 (2 j + r)
Ejercicios
1. Verificar que x = 0 es un punto singular regular de cada ecuacin diferencial y que
las races de la ecuacin indicial no difieren por un entero. Use el mtodo de Frobe-
nius para hallar las dos soluciones linealmente independientes y el correspondiente
intervalo de convergencia.
a) 2xy00 y0 + 2y = 0.
b) 3xy00 + (2 x)y0 + y = 0.
c) 2x2 y00 + 3xy0 + (2x 1)y = 0.
d) 2(x2 + x3 )y00 + (x 5x2)y0 + y = 0.
e) x2 y00 + (x 3/4)y = 0.
2. Hallar las dos soluciones en serie de Frobenius linealmente independientes de cada
ecuacin diferencial. Calcular a0 , . . . , an para n al menos 7 en cada solucin.
a) 2x2 (x2 + x + 1)y00 + x(5x2 + 3x + 3)y0 y = 0.
b) x2 (x + 8)y00 + x(3x + 2)y0 + (x + 1)y = 0.
c) 4x2 y00 + x(4x2 + 2x + 7)y0 (7x2 4x + 1)y = 0.
d) 2x2 (3x + 2)y00 + x(11x + 4)y0 + (1 x)y = 0.
e) 18x2 (x + 1)y00 + 3x(x2 + 11x + 5)y0 + (5x2 + 2x 1)y = 0.
f ) x2 (x + 6)y00 + x(4x + 11)y0 + (2x + 1)y = 0.
g) 10x2 (2x2 + x + 1)y00 + x(13 + 13x + 66x2 )y0 (1 + 4x + 10x2 )y = 0.
3. Hallar las dos soluciones en serie de Frobenius linealmente independientes de cada
ecuacin diferencial. Dar la forma explcita de los coeficientes en cada una de las
soluciones
a) 8x2 (1 x2 )y00 + 2x(1 13x2 )y0 + (1 9x2 )y = 0.
b) x2 (x + 3)y00 + x(4x + 5)y0 (2x + 1)y = 0.
c) 2x2 (x2 + 2)y00 x(7x2 + 12)y0 + (3x2 + 7)y = 0.
d) x2 (4x + 3)y00 + x(18x + 5)y0 + (12x 1)y = 0.
e) x2 (x2 2)y00 + x(x2 + 3)y0 y = 0.
f ) 2x2 (x + 3)y00 + x(5x + 1)y0 + (x + 1)y = 0.
g) x2 (x2 + 3)y00 + (2 x2 )y0 8y = 0.
6.3 Soluciones en Series Infinitas en Puntos Singulares Regulares 217
s1 y s2 = s1 + N, N Z+ .
Luego, se tiene que s1 + N es raz de la ecuacin indicial y por tanto satisface que
Comparando esto con el trmino que acompaa a an en la ecuacin (6.19) vemos que estos
trminos son iguales. Esto quiere decir que si usamos la raz s1 para hallar los coeficientes
que acompaan a xk en la representacin de la solucin en serie de Frobenius, no podemos
determinar todos los valores de estos coeficientes ya que el trmino aN es cero como
acabamos de ver. Por lo tanto, si las races de la ecuacin indicial difieren por un nmero
entero positivo, consideraremos dos casos por separado.
yg (x) = xs (a0 + a1 x + x2 x2 + ).
En esta situacin la raz s1 + N de la ecuacin indicial determina los valores de los coefi-
cientes ak de la solucin de la ecuacin diferencial en trminos de a0 , mientras que la raz
s1 determina dos conjuntos de valores para los coeficientes ak ; uno en trminos de a0 y otro
en trminos de aN . Sin embargo, la solucin en serie de Frobenius que se obtiene mediante
la raz s1 + N en trminos de a0 no es linealmente independiente con la solucin obtenida
en trminos de a0 a travs de la raz s1 . As, en este caso, la raz s1 genera dos soluciones
linealmente independientes cuya combinacin lineal es la solucin general de la ecuacin
diferencial.
La ecuacin indicial es
1
s(s 1)
=0
4
cuyas races son s1 = 1/2 y s2 = 1/2 y la diferencia entre ellas es N = 1. Usando la menor
de las races, s1 = 1/2 en el trmino que acompaa a xs1 +N , esto es:
De aqu tenemos que el coeficiente que acompaa a a1 es cero, as como el resto de trminos
en esta ecuacin (es decir, el que acompaa a a0 ), esto quiere decir que
Simplificando tenemos,
1 3 1 1
an +n +n + +n + an1 (0) + an2 (1) = 0
2 2 2 4
an (n2 n) = an2
6.3 Soluciones en Series Infinitas en Puntos Singulares Regulares 219
an2
an = , n 2.
n2 n
Hallemos los primeros trminos de esta recurrencia. Note que en lo que sigue a0 y a1 son
arbitrarios.
1
a2 = a0
2
1
a3 = a1
6
a2 a0
a4 = =
12 24
a3 a1
a5 = =
20 120
a4 a0
a6 = =
30 720
a5 a1
a7 = = .
42 5040
De aqu tenemos que los primeros siete trminos de las soluciones linealmente independien-
tes son:
1/2 1 2 1 4 1 6
y1 (x) =a0 x 1 x + x x +
2 24 720
1/2 1 3 1 5 1 7
y2 (x) =a1 x x x + x x + .
6 120 5040
La solucin general es
1/2 1 2 1 4 1 6
yg (x) =a0 x 1 x + x x +
2 24 720
1/2 1 2 1 4 1 6
+ a1 x 1 x + x x + .
6 120 5040
Dado que p(x) y q(x) son analticas en todo R, estas series convergen para todo x quizs en
x = 0. La primera serie converge en todo R menos en 0, mientras que la segunda converge
en todo R. Por lo tanto, la solucin general converge para 0 < |x| < . Note que para x < 0
la serie es compleja as que podemos escoger su parte imaginaria tomando a0 = ic con c un
nmero real.
El Coeficiente de aN es Cero y los Restantes Coeficientes de xs1 +N no lo son
En este caso slo la raz s1 + N de la ecuacin indicial determina un conjunto de valores para
los coeficientes ak en trminos de a0 para la solucin en serie de Frobenius de la ecuacin
diferencial
y1 (x) = xs (a0 + a1 x + a2 x2 + ).
Una segunda solucin linealmente independiente en este caso viene dada por
y2 (x) = u(x) N y1 (x) ln(x), x > 0, (6.24)
220 Soluciones de EDO Mediante Series
donde N es la diferencia entre las dos races de la ecuacin indicial, y1 (x) es una solucin en
serie de Frobenius que se obtiene mediante la raz s1 + N y u(x) es una serie de Frobenius
u(x) = xs1 (0 + 1 x + 2 x2 + ).
Sustituyendo estas funciones en la ecuacin diferencial tenemos que y2 (x) es una solucin
si se satisface que
xy00 y = 0.
p(x) =0 + 0 , bn = 0, n 0, x R
1
q(x) = x2 = 0 x + 0x2 + , c0 = 0, c1 = 1, cn = 0, n 2.
x
El punto singular regular es x = 0. La ecuacin indicial es s(s 1) = 0 con races s1 = 0
y s2 = 1, con diferencia N = 1. Usamos la raz s1 = 0 en el trmino que acompaa a x0+1 .
Esto es,
(a1 (0) + a0 (1))x.
As, el trmino que acompaa a a1 es cero pero no as el que acompaa a a0 . Por lo tanto
s1 = 0 no proporciona una solucin de la ecuacin. Para hallar una solucin en serie de
6.3 Soluciones en Series Infinitas en Puntos Singulares Regulares 221
Frobenius usemos la raz s2 = 1. Veamos la recurrencia que anula al trmino que acompaa
a x1+n .
Ejercicios
1. Verificar que x = 0 es un punto singular regular de cada ecuacin diferencial y que las
races de la ecuacin indicial difieren por un entero. Use el mtodo de Frobenius para
hallar al menos una solucin en serie de Frobenius y el correspondiente intervalo de
convergencia.
a) x2 y00 + xy0 + (x2 1)y = 0.
b) xy00 + xy0 + (x2 1/4)y = 0.
c) y00 + 3x y0 2y = 0.
d) xy00 + (1 x)y0 y = 0.
e) xy00 + y0 + y = 0.
222 Soluciones de EDO Mediante Series
Ejemplo 6.18 Hallar una solucin en serie de Frobenius para la ecuacin de Bessel de
ndice cero
x2 y00 + xy0 + x2 y = 0.
Calculemos ahora las soluciones de ecuaciones diferenciales de la forma
en el caso en que p0 (r) = a0 r(r 1) + b0 r + c0 = 0 tiene una raz repetida r1 . Sabemos que
en este caso la ecuacin (6.26) tiene una solucin de la forma
y1 (x) = xr1 anxn.
n=0
Sin embargo la otra solucin viene dada por una solucin logartmica. Para hallar esta use-
mos los polinomios (6.22) y los coeficientes (6.23). Tenemos que las soluciones linealmente
independientes de (6.26) son
y1 (x) =xr1 anxn
n=0
y2 (x) =y1 (x) ln(x) + xr1 a0n(r1)xn, x>0
n=1
dan (r)
donde a0n (r1 ) = dr |r=r1 .
224 Soluciones de EDO Mediante Series
Ejemplo 6.19 Hallar las dos soluciones linealmente independientes de la ecuacin dife-
rencial
2x2 (x + 2)y00 + 5x2 y0 + (x + 1)y = 0.
Dar las frmulas explcitas de los coeficientes en las soluciones.
Tenemos que r1 = 12 es una raz repetida del polinomio indicial p0 (r), por lo tanto las
soluciones de la ecuacin vienen dadas por
1/2 1 n
y1 (x) =x an x
n=0 2
1/2 0 1 n
y2 (x) =y1 (x) ln(x) + x an x .
n=1 2
a0 (r) =1
p1 (n + r 1)
an (r) = an1 (r)
p0 (n + r)
(n + r)(2n + 2r 1)
= an1 (r)
(2n + 2r 1)2
n+r
= an1 (r), n 0.
2n + 2r 1
De aqu se obtiene
n
j+r
an (r) = (1)n , n 0.
j=1 2 j + 2r 1
Evaluamos en r = 1/2
n j + 21 n
2j+1 nj=1 (2 j + 1)
an (1/2) = (1)n = (1)n = (1)n , n 0.
j=1 2j j=1 4 j 4n n!
Para calcular la segunda solucin calculamos a0n (1/2), n = 1, 2, . . . . Para esto usemos
diferenciacin logartmica
n
| j + r|
|an (r)| = , n1
j=1 |2 j + 2r 1|
n
ln |an (r)| = (ln | j + r| ln |2 j + 2r 1|) .
j=1
Derivamos respecto a r,
n
a0n (r) 1 2
=
an (r) j=1 j + r 2 j + 2r 1
n
0 1 2
an (r) =an (r) .
j=1 j + r 2 j + 2r 1
Evaluamos en r = 1/2
!
nj=1 (2 j + 1) n
1 1
a0n (1/2) =(1)n 1
4n n! j=1 j+ 2 j
n
n j=1
(2 j + 1) n 1
= (1) .
4n n! j=1 j(2 j + 1)
Obs En los casos en que la solucin y1 (x)0 sea una suma finita, el usar diferenciacin
logartmica para hallar los coeficientes an (r) en la segunda solucin y2 (x) trae algunas
dificultades. Veamos como solventar esto en el siguiente ejemplo.
Ejemplo 6.20 Hallar las dos soluciones linealmente independientes de la ecuacin dife-
rencial
x2 y00 x(5 x)y0 + (9 4x)y = 0.
p2 (r) =0.
Tenemos que r1 = 3 es una raz de multiplicidad 2 del polinomio indicial p0 (r). As, la
primera solucin es
y1 (x) = x2 an(3)xn.
n=0
Donde de (6.23)
a0 (r) =1
p1 (n + r 1)
an (r) = an1
p0 (n + r)
n+r5
= an1 (r)
(n + r 3)2
n
j+r5
=(1)n 2
, n 1. (6.27)
j=1 ( j + r 3)
Evaluamos en r = 3
n
j2
an (3) = (1)n .
j=1 j2
De aqu obtenemos
a1 (3) = 1, an (3) = 0, n 2.
Por lo tanto la solucin y1 (x) es
y1 (x) = x3 (1 + x).
de forma que
n
ln |an (r)| = (ln | j + r 5| 2 ln | j + r 3|) .
j=1
Derivamos respecto a r,
n
a0n (r) 1 2
=
an (r) j=1 j+r5 j+r3
6.3 Soluciones en Series Infinitas en Puntos Singulares Regulares 227
n
0 1 2
an (r) =an (r)
j=1 j + r 5 j+r3
Si evaluamos r = 3 en esta expresin tenemos que j = 2 nos da una indeterminacin. Por esto,
en estos casos podemos resolver este inconveniente haciendo el siguiente procedimiento:
De (6.27) para n = 1 tenemos
r4
a1 (r) = ,
(r 2)2
luego
r6
a01 (r) =
(r 2)2
a01 (3) = 3
y para n 2
nj=3 ( j + r 5)
an (r) = (1)n (r 4)(r 3) = (r 3)cn (r)
nj=1 ( j + r 3)2
donde
n nj=3 ( j + r 5)
cn (r) = (1) (r 4) n , n 2.
j=1 ( j + r 3)2
De esta manera
a0n (r) = cn (r) + (r 3)c0n (r), n2
lo que implica que a0n (3) = cn (3) si n 2.
nj=3 ( j 2)
a0n (3) =cn (3) = (1)n
nj=1 j2
(1)n+1 1 2 3
= 2 2 2
22 3 4 5
(1) n+1
= , n 2.
n(n 1)n!
228 Soluciones de EDO Mediante Series
Ejercicios
1. Verificar que x = 0 es un punto singular regular de cada ecuacin diferencial y que
la ecuacin indicial tiene una nica raz. Use el mtodo de Frobenius para hallar la
nica solucin en serie de Frobenius y el correspondiente intervalo de convergencia.
a) x2 y00 3xy0 + 4(x + 1)y = 0.
b) xy00 + (1 x)y0 + 12 y = 0.
2. Hallar las dos soluciones linealmente independientes de la ecuacin diferencial. Dar
las frmulas explcitas de las soluciones
a) 4x2 y00 + (4x + 1)y = 0.
b) x2 y00 + x(x2 1)y0 + (x2 + 1)y = 0.
c) x2 (2 x2 )y00 2(1 + 2x2 )y0 + (2 2x)y = 0.
d) 36x2 (2x + 1)y00 + 24x(9x + 1)y0 (70x 1)y = 0.
e) x2 (x2 + 1)y00 + 3x(x2 1)y0 + 4y = 0.
f ) x(x + 1)y00 + (1 x)y0 + y = 0.
g) 2x2 (x + 1)y00 x(6 x)y0 + (8 x)y = 0.
3. Una ecuacin diferencial de la forma
x(1 x)y00 + (c (a + b + 1)x)y0 aby = 0, a, b, c constantes
es llamada ecuacin de Gauss o ecuacin hipergeomtrica. Verificar que
a) x = 0 es un punto singular regular.
b) Las races de la ecuacin indicial son s1 = 0 y s2 = 1 c.
c)
(n + s + a)(n + s + b)
an+1 = an , n 0.
(n + s + 1)(n + s + c)
donde a0 , a1 , . . . son los coeficientes de la serie de Frobenius xs (a0 + a1 x +
a2 x2 + ), a0 6= 0.
d) Las soluciones para s1 = 0 y s2 = 1 c son respectivamente
ab a(a + 1)b(b + 1) x2
y1 (x) =1 + x+ +
c c(c + 1) 2!
a(a + 1) (a + n 1)b(b + 1) (b + n 1) xn
+ + ,
c(c + 1) (c + n 1) n!
c 6= 0, 1, 2, 3, . . .
1c (a c + 1)(b c + 1)
y2 (x) =x 1+ x + , c 6= 2, 3, 4, . . . .
2c
e) Ambas soluciones son iguales si c = 1.
f ) Si c no es entero la solucin general es
y1 + y2
y converge si |x| < 1, excepto quizs en x = 0. Si c es entero solo una de estas
funciones es una solucin.
6.3 Soluciones en Series Infinitas en Puntos Singulares Regulares 229
x = x1 , x0 = x2 , x00 = x3 , y = y1 , y01 = y2 ,
x10 =x2
x20 =x3
x30 = f1 (t, x1 , x2 , x3 )
y01 =y2
y02 = f2 (t, y1 , y2 ).
y = y1 , y0 = y2 , y00 = y3 , y(3) = y4 .
y01 =y2
y02 =y3
y03 =y4
y04 =et 3y4 + 2y3 y2 + 5y1 .
Cuando las funciones fi (t, y1 , . . . , yn ), i = 1, . . . , n son lineales, el sistema es llamado sistema
de n ecuaciones diferenciales lineales de primer orden.
Definicin 7.3 Un sistema de n ecuaciones es llamado lineal de primer orden si es de la
forma
o, como
y0 = A(t)y + Q(t)
donde
y1 f11 (t) f12 (t) fn1 (t) q1 (t)
.. .. .. .. ,
y = . , A(t) = . .
...
. y Q(t) = ... .
yn fn1 (t) fn2 (t) fnn (t) qn (t)
Obs Similar a lo que ocurre con las ecuaciones diferenciales, en general no existe un
mtodo para hallar la solucin de un sistema de ecuaciones diferenciales en trmino
de funciones elementales. Sin embargo, podemos usar las tcnicas de resolucin para
ecuaciones diferenciales dadas en los captulos anteriores para hallar las soluciones de
estos sistemas.
234 Sistemas de Ecuaciones Diferenciales Lineales
Luego, 6t
2 4 c1 e + c2 e2t 6c1 e6t 2c2 e2t
= .
4 2 c1 e6t c2 e2t 6c1 e6t + 2c2 e2t
Es decir, y es solucin del sistema de ecuaciones diferenciales.
siempre que este determinante sea diferente de cero. En este caso el sistema es llamado no
degenerado, en caso contrario se llama degenerado.
equivalentemente,
1
y02 = t c1 et + 1.
3
c2 4
c2 2c3 3 = 1, c3 = .
2
De esta manera tenemos que la segunda solucin del sistema es
1 c2 4
y2 (t) = t 2 + c1 et + t +
6 2
y la solucin general del sistema es el conjunto {y1 , y2 }.
d 2 y2 dy1
2
4y2 + = 0, y2 (0) = 0, y02 (0) = 1
dt dt
dy2 d 2 y1
2y1 4 + 2 = 0, y1 (0) = 1, y01 (0) = 2.
dt dt
Solucin. Aplicamos la transformada de Laplace a cada ecuacin, usamos sus propiedades
de linealidad y el Teorema 5.7. El sistema queda equivalentemente como,
Ahora multiplicamos la primera ecuacin del sistema anterior por 4s y la segunda por s2 4
y sumamos, de donde obtenemos la siguiente ecuacin:
Ejercicios
1. Reescribir el sistema en forma matricial y verificar que la funcin vectorial dada
satisface el sistema para cualquier valor de las constantes
c1 yc2 .
1 2 3t
a) y01 = 2y1 2y2 , y02 = 5y1 + y2 , y = c1 e4t + c2 e .
1 5
5 2t 2
b) y01 = 4y1 10y2 , y02 = 3y1 + 7y2 , y = c1 e + c2 et .
3 1
0 0 2 cos(t)
c) y1 = 3y1 +2y2 +32t, y2 = 5y1 +3y2 +63t, y = c1 +
3 cos(t) sin(t)
2 sin(t) 1
c2 + .
3 sin(t) + cos(t) t
d) y01 = 6y 2t
1 3y2 + 14e + 12e
t 0
, y2 =y1 2y2 + 7e2t 12et ,
3 5t 1 3t e2t + 3et
y = c1 e + c2 e + .
1 1 2e2t 3et
e) y01 = 2y2 + 2y
0
3 , y2 = 2y1+ 2y
0
3 , y3 = 2y1 +
2y2 ,
1 0 1
y = c1 0 e + c2 1 e + c3 1 e4t .
2t 2t
1 1 1
0
f ) y1 = y 0 0 = 2y + 2y y ,
1 + 2y2 + 2y3 , y2= 2y1 y2 + 2y 3 , y3 1 2 3
1 0 1
y = c1 0 e3t + c2 1 e3t + c3 1 e3t .
1 1 1
2. Reescribir la ecuacin diferencial lineal de orden n en un sistema de ecuaciones n n
equivalente
an (x)yn + an1 (x)yn1 + + a0 (x)y = g(x).
3. Resolver los siguientes sistemas de ecuaciones
a) dy 2 dy1
dt = y2 , dt = y1 .
2
b) dy 2 dy1
dt = y2t, dt = y1t.
2
c) dy dy2 dy1
dt = 2t, dt = 3y3 + 2t, dt = y3 + 4y2 + t.
3
d) dy dy1
dt = y2 + sin(t), dt = t y1 .
2
d 2 y2 dy2 2
e) dt 2
+ 4y2 = 3 sin(t), dt ddty21 + y = 2 cos(t).
d 2 y2 2
f) dt 2
+ y2 ddty21 y1 = cos(2t), 2 dy 2 dy1
dt dt y1 = 0.
7.2 Sistemas Homogneos: Teora Preliminar 239
b) 3 dy t dy1
dt + 3y2 + 2y1 = e , 4y2 3 dt + 3y1 = 3t, y2 (0) = 1, t1 (0) = 1.
2
2
c) ddty22 dy dy2 dy1 t 0
dt = 1 t, dt + dt = 4e + y2 , y2 (0) = 0, y1 (0) = 0, y2 (0) = 1.
1
d) dy dy1 dy1
dt y2 2y1 dt = 0, y2 dt = 15 cos(t)U(t ), y2 (0) = x0 , y1 (0) = y0 .
2
e) dy dy1
dt y2 + y1 = 2 sin(t)(1 U(t )), 2y2 dt y1 = 0, y2 (0) = y1 (0) = 0.
2
y1 (0) = 0.
Teorema 7.3 Sea A(t) una matriz n n continua en un intervalo (a, b). Un conjunto
{y1 , . . . , yn } de n soluciones de
y0 = A(t)y,
es un conjunto fundamental de soluciones si y slo si es linealmente independiente en (a, b).
Sean {y1 , . . . , yn } n soluciones de un sistema de ecuaciones diferenciales lineales homog-
neas, digamos
y11 y12 y1n
.. .. ..
y1 = . , y2 = . , . . . , yn = . .
yn1 yn2 ynn
Construimos la matriz Y (t) cuyas columnas se forman con estas funciones; es decir,
y11 (t) y12 (t) yn1 (t)
.. ,
Y (t) = ... ..
.
..
. . (7.3)
yn1 (t) yn2 (t) ynn (t)
y0 Y (t)y.
Dado que toda ecuacin diferencial de orden n se reduce a (o se puede escribir equivalen-
temente como) un sistema de ecuaciones de primer orden, esta definicin de Wronskiano
es equivalente a la dada en el Capitulo 4. Podemos calcular este Wronskiano mediante la
frmula de Abel que en este caso viene dada por el siguiente resultado.
Teorema 7.4 Frmula de Abel. Sea A(t) una matriz n n continua en un intervalo
(a, b) y sean y1 , . . . , yn soluciones en (a, b) del sistema
y0 = A(t)y.
Note de la frmula de Abel que W [y1 , . . . , yn ](t) no tiene ceros en (a, b) o W [y1 , . . . , yn ](t)
0 en (a, b).
7.2 Sistemas Homogneos: Teora Preliminar 241
Teorema 7.5 Sea A(t) una matriz n n continua en un intervalo (a, b) y sean y1 , . . . , yn
soluciones en (a, b) del sistema
y0 = A(t)y.
Entonces las siguientes afirmaciones son equivalentes.
1. La solucin general en (a, b) del sistema sistema es yg = an yn + + a1 y1 , con
a1 , . . . , an constantes.
2. {y1 , . . . , yn } es un conjunto fundamental de soluciones del sistema.
3. {y1 , . . . , yn } es linealmente independiente.
4. El Wronskiano de {y1 , . . . , yn } es diferente de cero en algn punto de (a, b).
5. El Wronskiano de {y1 , . . . , yn } es diferente de cero en todo punto de (a, b).
La matriz Y (t) dada en (7.3) es llamada una matriz fundamental para y0 = A(t)y si una
(y por tanto todas) afirmacin del teorema anterior es vlida para las columnas de Y . En
este caso la solucin general del sistema se escribe como yg = Y c donde c es un vector
constante.
Ejercicios
1. Sea V el conjunto de soluciones y = (y1 (t), y2 (y), y3 (t))T del sistema
0 1 0
y0 = 0 0 1 y.
6 11 6
Demostrar que
1 1 1
y1 (t) = 1 et ,
y2 (t) = 2 e2t ,
y3 (t) = 3 e3t ,
1 4 9
forman un conjunto fundamental de soluciones del sistema y0 = Ay. As, la solucin general
del sistema es
yg (t) = c1 v1 e1t + c2 v2 e2t + + cn vn ent .
Usando la matriz fundamental Y podemos escribir la solucin un problema de valor inicial
y0 = Ay, y(t0 ) = y0 como y(t) = Y (t)Y 1 (t0 )y0 .
4v1 v2 2v3 =0
v1 + v2 3v3 =0
4v1 + v2 + 2v3 =0
2v1 v2 2v3 =0
v1 v2 3v3 =0
4v1 + v2 =0
246 Sistemas de Ecuaciones Diferenciales Lineales
v1 v2 2v3 =0
v1 4v2 3v3 =0
4v1 + v2 3v3 =0
Ejercicios
1. Demostrar que todos los autovalores de la matriz de coeficientes de los siguientes
sistemas son
reales ydiferentes. Resolver el sistema.
5 3
a) y0 = y.
3 5
7.3 Sistemas Homogneos con Coeficientes Constantes 247
0 1 4 3
b) y = 5 y.
2 11
6 3
c) y0 = y.
1 2
6 4 8
d) y0 = 4 0 4 y.
8 4 6
4 1 4
e) y0 = 4 3 2 y.
1 1 1
3 1 1
f ) y = 3 5 1 y.
0
6 2 4
2. Demostrar que todos los autovalores de la matriz de coeficientes de los siguientes
sistemas son realesy diferentes. Resolver
el problema de valor inicial.
7 4 1
a) y0 = y, y(0) = .
6 7 3
7 2 1
b) y0 = y, y(0) = .
2 2 0
21 12 5
c) y0 = y, y(0) = .
24
15 3
2 2 3 1
0 1
d) y = 3 4 4 3 y, y(0) = 0 .
2 1 0 1
3 0 1 2
0
e) y = 11 2 7 y, y(0) = 1 .
1 0 3 3
2 5 1 8
f ) y0 = 4 1 1 y, y(0) = 10.
4 5 3 4
1 0 1
1 1
e) y(t) = 2 2 e + 2 1 e2t .
1 3 4t 5
1 1
y0 = Ay,
(A I)x = (A I)x = 0,
lo que implica que x es el autovector asociado al autovalor . De esta manera tenemos que
las funciones
y1 (t) = e x y y2 (t) = e x
son soluciones del sistema de ecuaciones. Usando la frmula de Euler podemos reescribir
la solucin generada por estas dos soluciones como
y01 + y1 5y2 =0
4y1 + y02 + 5y2 =0.
7.3 Sistemas Homogneos con Coeficientes Constantes 249
o, equivalentemente el sistema
(2 4i)v1 + 5v2 =0
4v1 (2 + 4i)v2 =0.
Note que la primera ecuacin es mltiplo de la segunda, (con mltiplo 2+4i 5 ) por lo
tanto fijando v1 = 5 obtenemos v2 = 2 + 4i y as el autovector es x = (5, 2 + 4i)T =
(5, 2)T + i(0, 4)T . Por lo tanto, del teorema anterior la solucin general de sistema es
3t 5 0
yg (t) =Ae cos(4t) sin(4t)
2 4
3t 5 0
+ Be sin(4t) cos(4t) .
2 4
Puede ocurrir que la matriz de coeficientes de un sistema de ecuaciones diferenciales
lineales con coeficientes constantes tenga autovalores tanto reales como complejos, en esta
situacin podemos hallar soluciones del sistema para cada caso y verificar que el conjunto
de soluciones que se obtiene para todos los casos es un conjunto fundamental de soluciones.
Ejercicios
1. Demostrar que la matriz de coeficientes de los siguientes sistemas tiene al menos un
autovalor complejo.
Resolver el sistema.
1 2
a) y0 = y.
5
5
0 1 1 2
b) y = 5 y.
4 5
5 6
c) y0 = y.
3 1
2 1 1
d) y0 = 0 1 1 y.
1 0 1
1 1 2
e) y0 = 4 0 1 y.
1 2 1
2. Demostrar que la matriz de coeficientes de los siguientes sistemas tiene al menos un
autovalor complejo.
Resolver el problema
de valor inicial.
4 2 1
a) y0 = 16 y, y(0) = .
5 2 1
0 4 6 2
b) y = y, y(0) = .
3 2 1
7 15 5
c) y0 = y, y(0) = .
3
1 1
5 2 1 4
d) y = 3 3 2 2 y, y(0) = 0.
0 1
1 3 2 6
252 Sistemas de Ecuaciones Diferenciales Lineales
4 4 0 1
e) y = 8 10 20 y, y(0) = 1.
0
2 3 2 0
4 4 4 16
0
f ) y = 10 3 15 y, y(0) = 14.
2 3 1 6
Para 2 = 5:
2 2 2 v1 0
(A + 5I)v2 = 2 2 2
v2 = 0 .
2 2 2 v3 0
De donde tenemos la ecuacin v1 + v2 + v3 = 0, podemos fijar dos valores, digamos v1 =
1, v2 = 1 y as v3 = 0, obtenemos el autovector (1, 1, 0)T . Otro vector linealmente
independiente se forma tomando v1 = 0, v2 = 1 y v3 = 1. Esto es v3 = (0, 1, 1)T . Las dos
soluciones linealmente independientes son
1 0
5t
y2 (t) = 1 e , y3 (t) = 1 e5t .
0 1
y01 4y1 + y2 =0
3y1 y02 + y2 y3 =0
y1 y03 + y3 =0.
7.3 Sistemas Homogneos con Coeficientes Constantes 255
2v1 v2 =0
3v1 v2 v3 =0
v1 v3 =0.
Ejercicios
1. Demostrar que todos los autovalores de la matriz de coeficientes de los siguientes
sistemas son
realescon multiplicidad mayor a 1. Resolver el sistema.
3 4
a) y0 = y.
1 7
3 1
b) y0 = 15 y.
1 1
10 9
c) y0 = y.
4 2
0 2 1
d) y0 = 4 6 1 y.
0 4 2
6 4 4
e) y0 = 2 1 1 y.
2 3 1
1 12 8
f ) y = 1 9 4 y.
0
1 6 1
2. Demostrar que todos los autovalores de la matriz de coeficientes de los siguientes
sistemas son
reales conmultiplicidad
mayor
a 1. Resolver el problema de valor inicial.
11 8 6
a) y0 = y, y(0) = .
2 3 2
0 3 4 3
b) y = y, y(0) = .
1 7 1
7 3 2
c) y0 = y, y(0) = .
3
1 0
1 1 0 6
d) y = 3 1 1 2 y, y(0) = 5 .
0 1
1 1 1 1
2 2 1 1
e) y = 2 2 1 y, y(0) = 3.
0
3 3 2 1
7.4 Sistemas no Homogneos con Coeficientes Constantes 257
3 3 4 6
0
f) y = 4 5 8 y, y(0) = 9 .
2 3 5 1
y = Y (t)c(t).
Para el sistema no homogneo con condicin inicial y(t0 ) = y0 , la solucin se puede escribir
como Z
t
y = Y (t)c(t) = Y (t) c + Y 1 (s)Q(s)ds
t0
con y(t0 ) = Y (t0 )c, lo que implica que c = Y 1 (t0 )y(t0 ), con lo cual tenemos que la solucin
es Z t
1 1
y(t) = Y (t) Y (t0 )y(t0 ) + Y (s)Q(s)ds .
t0
Ahora
Z Z t Z
1 (1 t)et (2t + 3)et 2e 2(1 t)
Y Q(t)dt = dt = dt
et 2et 0 2
2t t 2
= .
2t
Hallemos la solucin del sistema homogneo. Para esto calculemos los autovalores de la
matriz de coeficientes
3 3 1
3 2 2 = (2+ )(3 )+12(14+6 ) = ( +1)( 2 2 +2) = 0.
1 2
del cual, restando la segunda ecuacin a la primera, tenemos que 5v1 5v2 = 0, luego para
v1 = 1, hallamos v2 = 1 y v3 = 1. As un autovalor correspondiente es v1 = (1, 1, 1)T y
una primera solucin del sistema es
1
y1 (t) = 1 et .
1
Para 2 = 1 + i
31i 3 1 v1 0
(A + (1 + i)I)v2 = 3 2 1 i 2 v2 = 0 .
1 2 1 i v3 0
Tomemos el subsistema (
3v1 (3 + i)v2 + 2v3 = 0
v1 + 2v2 (1 + i)v3 = 0.
Sumamos la primera ecuacin con 3 veces la segunda y obtenemos (3i)v2 (1+3i)v3 = 0.
(1+3i)(3+i)
Hacemos v3 = 1 y as (3 i)v2 = 1 + 3i, de donde v2 = 1+3i
3i = (3i)(3+i) = i. Ahora usamos
la segunda ecuacin del subsistema anterior para hallar v1 . Esto es, v1 = 2v2 (1 + i)v3 =
i 1 y el autovalor es (1 + i, i, 1)T = (1, 0, 1)T + i(1, 1, 0)T . Las dos soluciones vienen
dadas por
1 1
y2 (t) =et 0 cos(t) 1 sin(t) ,
1 0
7.4 Sistemas no Homogneos con Coeficientes Constantes 261
1 1
y3 (t) =et 1 cos(t) + 0 sin(t) .
0 1
Ejercicios
1. Resolver los siguientes sistemas de ecuaciones no homogneos
a) y01 = 2y1 + y2 + 26 sin(t), y02 = 3y1 + 4y2 .
b) y01 = y1 + 8y2 + 9t, y02 = y1 + y2 + 3et .
3t
c) y01 = y1 + 2y2 , y02 = 3y1 + 4y2 + e2te +1 .
d) y01 = y1 + y2 + e2t , y02 = 2y1 + 3y2 .
e) y01 = y1 2y2 y3 , y02 = y1 + y2 + y3 + 12t, y03 = y1 y3 .
7.4 Sistemas no Homogneos con Coeficientes Constantes 263
Como hemos visto no siempre es posible hallar una expresin analtica para la solucin de
una ecuacin diferencial ordinaria o un sistema de ecuaciones diferenciales ordinarias, es
por esoto que se hace necesario disponer de otros mtodos para estudiar la solubilidad y
caracterizacin matemtica de la ecuacin o sistema de ecuaciones diferenciales. En algunas
ecuaciones se disponen de mtodos de anlisis numrico para encontrar una solucin
aproximada. Otro mtodo ampliamente utilizado en investigaciones y aplicaciones es el
mtodo cualitativo mediante el cual se busca establecer matemticamente las caractersticas
ms significativas de la solucin.
En este captulo estudiaremos los aspectos cualitativos de las soluciones de ecuaciones
diferenciales y sistemas de ecuaciones diferenciales autnomas.
El concepto de curva solucin sugiere el siguiente mtodo grfico para construir soluciones
aproximadas de una ecuacin diferencial dy/dx = f (x, y).
Definicin 8.2 Dada una ecuacin diferencial dy/dx = f (x, y), tracemos por cada uno
de los puntos (x, y) un segmento corto de recta que tenga la pendiente m = f (x, y). El
conjunto de todos estos segmentos de rectas se denomina campo de direcciones o campo de
pendientes asociado a la ecuacin diferencial dada.
dy
Figura 8.1: Campo de pendientes de la ecuacin dx = 10 14 y.
Ejemplo 8.1 En la Figura 8.2 se muestra la grfica del campo de pendientes de la ecuacin
dy/dx = x + y. En la Figura 8.3 aparece el campo de pendientes de misma ecuacin y una
solucin que pasa por el punto (1, 2), para la cual la recta tangente a la grafica de la solucin
en dicho punto posee pendiente 3.
8.1 Ecuaciones Diferenciales de Primer Orden 267
dy
Figura 8.2: Campo de pendientes de la ecuacin dx = x + y.
dy
Figura 8.3: Solucin de la ecuacin dx = x + y que pasa por el punto (1, 2).
Obs
dy
1. Si en la ecuacin diferencial dx = f (x, y), la funcin f no depende de y, entonces
las isoclinas son lineas verticales.
2. En general, para el caso en que f (x, y) dependa de x e y, las isoclinas pueden
ser cualquier tipo de curva. Por ejemplo, si f (x, y) = x2 2x y2 4y 3, las
isclinas de la ecuacin diferencial correspondientes a nmeros reales c > 0, son
hiprbolas con eje transverso horizontal y con centro en (1, 2).Tambin, de
los vrtices al centro hay c unidades y los focos estan ubicados a 2c unidades
desde el centro.
dy
Figura 8.4: Isoclinas de la ecuacin dx = y.
dy
Figura 8.5: Campo de pendientes de la ecuacin dx = y.
dy
Figura 8.6: Algunas soluciones de la ecuacin dx = y.
Ejercicios
1. Dada la ecuacin diferencial
dp
= p(p 2)(4 p)
dt
para la poblacin p (en miles) de cierta especie en el instante t.
a) Bosquejar el campo de direcciones usando el mtodo de isclinas.
b) Si la poblacin inicial es 5000 (es decir, p(0) = 5) Qu puede decir acerca de
la poblacin lmite lm p(t)?
t
270 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO
3. Trazar las isclinas y bosquejar varias curvas solucin, incluyendo la curva que
satisfaga las condiciones iniciales dadas.
dy
a) dx = ln |x| y(1) = 1.
dy
b) dx = 1x y(1) = 1.
dy
c) dx = x + 2y y(0) = 1.
dy
d) dx = x2 y(0) = 1.
8.1 Ecuaciones Diferenciales de Primer Orden 271
Trazar las isclinas y bosquejar varias curvas solucin, incluyendo la curva que
satisfase y(0) = 0. Qu ocurre con y(x) cuando x ?
b) 4.
c) 0.
d) 0.
3. a) Figura 8.10
c) Figura 8.11
272 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO
dh dg d dg
|x0 +k = |x0 (x k)|x0 +k = |x0 = f (g(x0 )) = f (h(x0 + k));
dx dx dx dx
dy
luego, h tambin es solucin de la ecuacin diferencial dx = f (y) en el intervalo J.
Supongamos que g satisface la condicin inicial g(x0 ) = y0 . Ahora, dada otra condicin
inicial, digamos y(u0 ) = v0 , consideremos k = u0 x0 y definamos h(x) = g(x k) + (v0
y0 ). El argumento anterior, demuestra que h es solucin del problema de valor inicial
dy
dx = f (y) sujeto a y(u0 ) = v0 . Luego, del Teorema Existencia y Unicidad, tenemos que h es
la nica solucin de este problema de valor inicial. Esto quiere decir que cualquier solucin
de un problema de valor incial para una ecuacin diferencial autnoma, determina las
soluciones para cualquier otro problema de valor inicial con la misma ecuacin diferencial.
dy
= (y 1)(y 2)(y 3). (8.3)
dx
Ejemplo 8.5 Para la ecuacin diferencial autnoma del Ejemplo 8.3, los equilibrios son
0
1 = 1, y2 = 2 y y3 = 3; en la Figura 8.13 se grfica y = (y 1)(y 2)(y 3)
y en el plano
yy0 . Las lneas continuas ilustran estas soluciones de equilibrio.
Note como todas las soluciones y(t) que estn lo suficientemente cerca al punto de equilibrio
2 = 2 tiende al mismo cuando x +. Todo lo contrario sucede con los otros dos puntos
y
de equilibrio y1 = 1 y y3 = 3 para los cuales las soluciones cercanas se alejan de ellos.
equilibrio inestables.
El siguiente teorema nos permite analizar de manera analtica la naturaleza asinttica de los
equilibrios de una ecuacin diferencial autnoma.
Si esbozamos el grfico de y0 = f (y) en el plano yy0 , el Teorema 8.1 dice que si la funcin f
corta el eje y en un equilibrio y de manera decreciente ( f 0 (y ) < 0) entonces este equilibrio
es asintticamente estable, mientras que si corta el eje y de manera creciente ( f 0 (y ) > 0),
y es inestable. Si el equilibrio es tal que f 0 (y ) = 0, el equilibrio se llama semi-estable.
8.1 Ecuaciones Diferenciales de Primer Orden 275
Ejemplo 8.6 Para la ecuacin diferencial (8.3), y0 = (y 1)(y 2)(y 3) tenemos que
f 0 (y) = 3y2 12y + 11 y
0
Obs Note que toda ecuacin autnoma y = f (y) es de variable separable y por tanto en
principio se puede resolver integrando, lo cual puede ser un proceso complicado. Sin
embargo, el anlisis cualitativo permite conocer el comportamiento de las soluciones
sin necesidad de resolver explcitamente la ecuacin diferencial.
Con esta informacin podemos podemos caracterizar los equilibrios de la siguiente manera:
Pozo: Todos los puntos de equilibrios asintticamente estables. (Tambin llamado sumidero
o atractor.)
Fuente: Todos los puntos de equilibrios inestables. (Tambin llamado repulsor.)
Nodo: Todos los puntos de equilibrio semi-estables.
Con esta informacin podemos estudiar el comportamiento de los puntos de equilibrios de
una ecuacin autnoma mediante la llamada lnea de fase.
dy
Definicin 8.6 La linea de fase asociada a la ecuacin diferencial autnoma dx = f (y),
es un diagrama que describe el comportamiento geomtrico de los puntos de equilibrio de
la misma y el cual se construye de la siguiente forma:
Primero debemos cerciorarnos de que la ecuacin diferencial en cuestin satisfaga las
condiciones del Teorema Existencia y Unicidad.
Hallar todas las soluciones de equilibrio de la ecuacin diferencial.
Trazar la recta real y sobre ella marcar todos los equilibrios.
276 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO
y y y
J 1I I 2J J 3I
Fuente Pozo Nodo
Ejemplo 8.7 Establecer los puntos de equilibrio para la ecuacin diferencial autnoma
dy
= y2 7y + 10, (8.4)
dx
determinar su tipo y graficar su lnea fase.
Solucin. En la Figura 8.16 se muestra la grfica de f (y) en el plano yy0 .
0
1 = 2 y y2 = 5. Adems f (2) =
Los puntos de equilibrio la ecuacin autnoma (8.4) son y
0
3 < 0 y f (5) = 3 > 0, luego y1 = 2 es un atractor mientras que y2 = 5 es una fuente. La
2 5
IJ JI
Pozo Fuente
Ejercicios
1. Hallar todos los equilibrios de las ecuaciones diferenciales dadas y determinar cuales
estables
son asintticamente
0
a) y = ry ln Ky .
ry(Ky)
b) y0 = K+ay .
y
c) y0 = ry 1 K , 0 < < 1.
d) y0 = y(re1y/K d).
2. a) El modelo de una poblacin est gobernado por la ecuacin diferencial
y0 = y(e3y 1).
278 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO
y0 = y(2y ey )
C0 = C( C).
a) y0 = y2 2y + 1.
b) y0 = cos(y).
c) y0 = (y 2)4 .
d) y0 = y2 (4 y2 ).
e) y0 = 10 + 3y y2 .
9. Utilizar la lnea fase para predecir el comportamiento asinttico (cuando t +)
de la solucin que satisfaga la condicin inicial dada.
a) y0 = y(y2 2y + 1), y(0) = 0,5.
b) y0 = y(y2 2y 8), y(0) = 10.
0
c) y = cos(y), y(0) = 0.
10. Sea g(x) una funcin tal que g(K) = 0 y g(x) > 0 para 0 < x < K y suponga que
0 < x0 < K.
a) Demostrar que la funcin x(t) definida de manera implcita por la relacin
Z x(t)
du
=t (8.5)
x0 ug(u)
es una solucin del problema de valor inicial x0 = xg(x), x(0) = x0 .
b) Demostrar que la integral dada en (8.5) es negativa si x(t) < x0 y positiva si
x(t) > x0 y deducir que la solucin del problema de valor inicial debe satisfacer
que x(t) > x0 para t > 0.
c) Demostrar que cuando t , la integral en (8.5) diverge. Deducir que g(x(t))
0 y as x(t) K.
11. Para la ecuacin diferencial
y0 = f (y)
con f (0) = 0, f (x) < 0 (0 < x < K0 ), f (K0 ) = 0, f (x) > 0 (K0 < x < K), f (K) = 0,
f (x) < 0 (x > K), demostrar que los equilibrios 0 y K son asintticamente estables y
que el equilibrio en K0 es inestable.
r < 0.
1 = 0 semi-estable, y2 = K asintticamente estable si r < 0, inestable r > 0.
c) y
3. 0.
8. a) Figura 8.19
1
JI
Nodo
d) Figura 8.20
2 0 2
JI JI IJ
Fuente Nodo Pozo
(
dy1
dt = f1 (y1 , y2 )
dy2 (8.6)
dt = f2 (y1 , y2 ).
Geomtricamente podemos interpretar la solucin {y1 (t), y2 (t)} de nos maneras: Primero
podemos graficar y1 (t), y2 (t) con respecto a t en el mismo plano como muestra la siguiente
figura:
y1 y2
y1 (t)
y2 (t)
Esta representacin nos permite conocer cmo varan las soluciones respecto a t. O, segundo,
podemos considerar a y1 (t), y2 (t) como ecuaciones paramtricas de una curva en el plano
y1 y2 , con t como parmetro a lo largo de la curva.
8.2 Sistemas de Ecuaciones Autnomas 281
y2
y1
b b
y1
Mientras que en el mismo plano son las funciones dadas en la Figura 8.24
dy2
dt = g(y1 , y2 )
t0 pertenece al dominio de y1 , y2 y (y10 , y20 )
y1 (t0 ) = y10 , y2 (t0 ) = y20
tiene una nica solucin y1 = y1 (t), y2 = y2 (t), definida en cierto intervalo real (a, b) que
contiene el punto t0 . Para cada punto en el plano fase existe una y slo una trayectoria
que contiene al punto.
Definicin 8.7 Punto de equilibrio. Un punto de equilibrio, o equilibrio, es una
solucin (y
1 , y2 ) del par de ecuaciones f 1 (y1 , y2 ) = 0, f 2 (y1 , y2 ) = 0. As un equilibrio es
Figura 8.25: En el plano fase el campo de vectores (3y z, y + z)T asociado al sistema
y0 = 3y z, z0 = y + z y varias curvas y(t), z(t) las cuales se alejan del origen. El campo
de vectores es tangente a las curvas solucin. Las rbitas se aproximan al infinito cuando
t .
284 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO
Las rbitas encajan de tal manera que los vectores tangentes coinciden con el campo de
vectores. Un retrato de fase es un bosquejo del plano fase y unas pocas trayectorias tpicas,
junto con sus puntos de equilibrio. El diagrama de fase puede, o no, tener al campo de
vectores.
Una forma para hallar el retrato de fase del sistema (8.6) es suponer que y1 y y2 son
soluciones del sistema, entonces para los valores de t para los cuales f1 (y1 (t), y2 (t)) 6= 0,
de la regla de la cadena tenemos
dy2
dy2 dt f2 (y1 (y), y2 (t))
= dy1
= .
dy1 f1 (y1 (t), y2 (t))
dt
dy2 f2 (y1 , y2 )
= .
dy1 f1 (y1 , y2 )
y00 = y y3 .
y01 =y
y0 =y1 y31 .
dy y1 y31
= .
dy1 y
y41
y2 = y21 + c.
2
Un retrato de fase que contiene algunas curvas viene dado en la Figura 8.26
8.2 Sistemas de Ecuaciones Autnomas 285
Ejemplo 8.11 Pndulo. Un peso con masa m est suspendido al final de una barra de
longitud L. Sea medido en el sentido positivo como indica la figura
Como la barra es rgida, la masa viaja a lo largo de una circunferencia de radio L. Ignorando
la resistencia del aire y la masa de la barra, aplicamos la segunda ley de Newton en la
direccin tangencial a la circunferencia
Figura 8.28: Diagrama del plano fase para la ecuacin del pndulo simple.
8.2 Sistemas de Ecuaciones Autnomas 287
Los Ejemplos 8.8 y 8.11 son casos especiales de los llamados sistemas Hamiltonianos.
Definicin 8.8 Un sistema bidimensional de la forma
h
y01 =
y2
h
y02 =
y1
Obs Note que una condicin necesaria para que el sistema (8.6) sea Hamiltoniano es que
f1 f2
= .
y1 y2
Solucin. Para hallar las funciones Hamiltanianas para este sistema consideremos
h
(y1 , y2 ) = 3y22 2y31 y2 .
y2
h
Integrando respecto a y2 tenemos que h(y1 , y2 ) = y32 y31 y22 + k(y1 ). Como y1 = (3y21 y22
2y1 ) obtenemos
y as,
k(y1 ) = y21 + c.
Finalmente, h(y1 , y2 ) = y22 y31 y22 + y21 + c son las funciones Hamiltonianas.
(y1 (t), y2 (t)) con (y1 (t0 ), y2 (t0 )) suficientemente cerca al equilibrio, permanece cerca al
equilibrio para todo t. Formalmente un equilibrio (y 1 , y2 ) es estable si para todo > 0, existe
un > 0, tal que si |y1 y1 (t0 )| < , |y2 y2 (t0 )| < , entonces la solucin y1 (t), y2 (t)
satisface que |y1 (t) y 1 | < , |y2 (t) y2 | < , para todo punto t.
soluciones con (y1 (t0 ), y2 (t0 )) suficientemente cerca al equilibrio tienden al equilibrio
cuando t . Esto es, si para algn > 0 se tiene que
|y1 (t0 ) y
1 | < , |y2 (t0 ) y2 | < entonces y1 (t) y1 , y2 (t) y2 , cuando t .
y01 = f1 (y1 , y2 )
y02 = f2 (y1 , y2 ), (8.10)
as el sistema
u0 = f1 (u + y
1 , v + y2 )
v0 = f2 (u + y
1 , v + y2 ),
f1 (u + y
1 , v + y2 ) = f 1 (y1 , y2 ) + f1 (y
1 , y2 )u + f1 (y
1 , y2 )v + h1
y1 y2
f2 (u + y
1 , v + y2 ) = f 2 (y1 , y2 ) + f2 (y
1 , y2 )u + f2 (y
1 , y2 )v + h2
y1 y2
h1 (u, v) h2 (u, v)
lm = lm = 0.
u0 2
u +v 2 u0 u2 + v2
v0 v0
8.2 Sistemas de Ecuaciones Autnomas 289
bado
y01 = f1 (y1 , y2 )
y02 = f2 (y1 , y2 )
para t 0.
290 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO
u0 = ( by ) bx v, v0 = cy u + ( + cx )v.
u0 = u, v0 = v
y en el equilibrio c, b tenemos
b c
u0 = v, v0 = u.
c b
Ejercicios
1. Dibujar el retrato del plano fase para cada ecuacin diferencial
a) x00 2xx0 = 0.
b) x00 + |x| = 0.
c) x00 + ex = 1
2. Demostrar que los sistemas son Hamiltonianos y hallar las ecuaciones para las rbitas.
Dibujar el plano fase.
a) y01 = y1 + 16y2 , y02 = 8y1 y2 .
1y2 +y2
b) y01 = y12 2 , y02 = 2y 1
y2 , y2 6= 0.
2
3. Para cada uno de los sistemas que sean Hamiltonianos, hallar la funcin Hamiltoniana.
Si el sistema no es Hamiltoniano decirlo.
8.2 Sistemas de Ecuaciones Autnomas 291
I
S0 (t) =N S S
N
I
I 0 (t) = S ( + )I
N
V 0 (t) =I V
excepto para las rbitas en el eje u, toda rbita tiene un mximo y un mnimo en
trminos de la variable u por lo que se vuelve al origen. El retrato de fase se muestra
en la Figura 8.31. Si > 0, el retrato de fase es el mismo pero con las flechas en
direccin opuesta.
8.34. si > 0, el retrato de fase es el mismo con la excepcin que las fechas estn al
contrario.
Estos seis casos se pueden clasificar como de cuatro tipos distintos dependiendo del com-
portamiento de las rbitas:
Nodo: En los casos I a III todas las rbitas se aproximan al origen cuando t + (o
cuando t dependiendo de los signos de y ) con una direccin lmite y el
origen se llama un nodo del sistema.
Silla: En el caso IV slo dos rbitas se aproximan al origen cuando t y las otras
rbitas se alejan del origen. En este caso el origen se llama un punto de silla.
Foco: En el caso VI toda rbita gira en forma espiral alrededor del origen en el sentido que
su argumento angular tiende a + o . El origen en este caso se llama foco, vrtice
o punto espiral.
Centro: En el caso V toda rbita es peridica y el origen se llama centro.
Ahora, podemos describir el comportamiento asinttico de las soluciones cerca del equilibrio
(0, 0) en cada uno de los 6 casos anteriores mediante los autovalores de las matrices de
coeficiente de estos sistemas. Para el caso general de un sistema 2 2 podemos utilizar esta
informacin para para describir el comportamiento de las rbitas de la matriz de coeficientes
de un sistema lineal, ya que las matrices semejantes tienen los mismos autovalores.
Comencemos por notar las siguientes caractersticas de los autovalores de una matriz A
2 2:
a b 1 0
2
c b 0 1 = (a + d) + (ad bc) = 0.
Obs Una descripcin ms simple es que el origen asintticamente estable si ambos auto-
valores tienen parte real negativa y es inestable si al menos un autovalor tiene parte
real positiva. si ambos autovalores tiene parte real cero, el origen es estable pero no
asintticamente estable. Como A es invertible, ad bc 6= 0 descarta la posibilidad que
= 0 sea un autovalor, por lo tanto, autovalores con parte real cero slo ocurre si los
autovalores son imaginarios puros como en el Caso V.
Tr(y
1 , y2 ) = f1 (y
1 , y2 ) + f2 (y
1 , y2 ) < 0,
y1 y2
|A(y
1 , y2 )| = f1 (y
1 , y2 ) f2 (y
1 , y2 ) f1 (y
1 , y2 ) f2 (y
1 , y2 ) > 0,
y1 y2 y2 y1
entonces el equilibrio (y
1 , y2 ) es asintticamente estable.
que las funciones f1 y f2 del sistema (8.12) son suficientemente suaves para poder ser
expandidas en polinomios de Taylor y omitir los trminos de orden superior.
Si la linealizacin tiene un nodo en un equilibrio, entonces el equilibrio del sistema no
lineal tambin es un nodo definido para significar que toda rbita tiende al equilibrio cuando
t con una direccin lmite. Si la linealizacin tiene tiene un punto foco en un
equilibrio, entonces el equilibrio del sistema no lineal tambin tiene un foco, definido para
significar que toda rbita tiende al equilibrio (t ) con su variable angular tendiendo
al infinito. En el caso cuando la linealizacin tiene un punto de silla en un equilibrio,
el sistema no lineal tambin tendr un punto de silla en el equilibrio. Un punto de silla
se define por la caracterstica que es una curva alrededor del equilibrio tal que las rbitas
que comienzan en esta curva tienden al equilibrio pero las rbitas que comienzan fuera de
esta curva no pueden estar cercas del equilibrio. Otra caracterstica es que hay dos rbitas
tendiendo al equilibrio cuando t + y hay dos rbitas alejndose desde el equilibrio o
tendiendo al equilibrio cuando t . Estas rbitas se llaman separatrices, las dos rbitas
que tienden al punto de silla son sepatatrices estables mientras que las dos rbitas que se
alejan del punto de silla son separatrices inestables, las otras rbitas parecen hiprbolas.
Un centro se define como un equlibrio para el cual existe una sucesin infinita de rbitas
peridicas alrededor del equilibrio con las rbitas acercndose a este equilibrio. Si la
linealizacin tiene un centro en el equilibrio, entonces el equilibrio del sistema no lineal
no necesariamente es un centro, tambin puede ser un foco asintticamente estable o un
foco inestable.
dy 18x
=
dx 8y
Como f (x, y) = 8y, las trayectorias del semiplano superior (y > 0) fluyen hacia
la izquierda y viceversa. El determinante de la matriz de coeficiente es 144 > 0 y
la traza es cero, por lo tanto el equilibrio es un centro el cual es estable pero no
asintticamente estable. Al resolver la ecuacin obtenemos 9x2 + 4y2 = c lo cual
indica que sus trayectorias son elipses que encierran el equilibrio.
2. Los puntos crticos son 2y = 0, 2x = 0 es decir x = 0, y = 0, por lo tanto (0, 0) el
nico equilibrio del sistema. El campo de direcciones para la ecuacin
dy 2x x
= =
dx 2y y
Como f (x, y) = 2y, las trayectorias del semiplano superior (y > 0) fluyen hacia la
derecha y viceversa. Al resolver la ecuacin obtenemos y2 x2 = k lo cual indica
que sus trayectorias son hiprbolas, que encierran al equilibrio. El determinante de la
matriz de coeficientes es 4 < 0 y as el origen es un punto de silla inestable pues los
autovalores de la matriz de coeficientes son 2.
3. Los equilibrios se hallan resolviendo el sistema.
(
2x + y = 0
5x 4y = 0
de donde x = 0, y = 0, son el conjunto de puntos crticos, por lo tanto (0, 0) es el
punto crtico. El campo de direcciones para la ecuacin
dy 5x 4y
= (8.14)
dx 2x + y
aparece en la Figura 8.38
Observe que las soluciones fluyen hacia la derecha para el semiplano 5x 4y > 0;
es decir, para todos los puntos debajo de la recta 5x 4y = 0. El determinante de la
matriz de coeficientes de este sistema es 3 > 0, la traza 6 < 0 y = 24 > 0, por lo
tanto el origen es un nodo asintticamente estable.
4. Los equilibrios se hallan resolviendo el sistema.
(
5x 3y 2 = 0
4x 3y 1 = 0
Observe que las soluciones fluyen hacia la derecha para el plano 4x 3y 1 > 0, es
decir para todos los puntos debajo de la recta 4x 3y 1 = 0.
Solucin. Los puntos crticos son las soluciones de y = 0 y 2(x2 a)y x = 0. El equilibrio
as es (x , y ) = (0, 0). La matriz de coeficientes del sistema en el equilibrio es
0 1 0 1
= .
1 2((x )2 1) 1 2
302 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO
Solucin. En el ejemplo 8.13 vimos que los equilibrios son (0, 0) cuya matriz de coeficien-
tes de la linealizacin
0
.
0
Como el determinante es negativo, (0, 0) es inestable. Para el equilibrio c , b , la matriz
de la linealizacin es !
b
0 c
c
b 0
con traza cero, es decir el equilibrio es un centro. Este caso no lo contempla el Teorema 8.5.
Sin embargo, resolviendo la ecuacin
dy y( + cx)
=
dx x( by)
obtenemos la solucin ln |x| ln |y| + cx + by = k y todas las rbitas son dadas im-
plcitamente por esta ecuacin. Hacemos el cambio de variable x = /c + u, y = /b + v
obtenemos
ln + v ln + v + c +u +b + v = k.
c b c b
8.2 Sistemas de Ecuaciones Autnomas 303
Note que
cu
ln + u = ln + ln 1 +
c c
y si k k0 es pequeo podemos usar la expresin ln(1 + x) x x2 /2 para aproximar esta
expresin por
cu c2 u2
ln + 2 .
c
De igual manera podemos aproximar ln | /b + v| por ln | /b| + bv/ b2 v2 / 2 . De esta
forma las rbitas se pueden aproximar por
c2 2 b2
ln cu + u ln bv + v2 + + cu + + bv = h,
c b
o, equivalentemente
c2 2 b2 2
u + v = h + ln + ln = h h0 ,
c b
lo cual representa una elipse si h > h0 con el equilibrio (x , y ) como su centro. Esto muestra
que para h h0 pequeo y positivo las rbitas son curvas cerradas alrededor del equilibrio,
como las soluciones viajan a travs de una curva cerrada, estn deben ser peridicas. As
este equilibrio no es asintticamente estable ni inestable. Este tipo de equilibrio es llamado
neutralmente estable.
Ejercicios
1. Determinar si los equilibrios de los sistemas dados son asintticamente estables o
inestables
a) y01 = y1 y2 , y02 = y1 + y2 2.
b) y01 = y2 , y02 = y1 + y2 1.
c) y01 = ey2 , y02 = ey1 .
2. Determinar el comportamiento de las rbitas de los siguientes sistemas
a) x0 = x( ax by), y0 = y( cx dy).
b) x0 = x( ax by), y0 = y( + cx).
c) x0 = y, y0 = x y3 .
3. La dinmica de un sistema biolgico de lapas y algas marinas es dada por el siguiente
sistema
ds
=s s2 sl
dt
dl l l2
=sl
dt s s
donde las densidades de las algas y las lapas son dadas por l 0 y s 0 respectiva-
mente.
304 Anlisis Cualitativo de Ecuaciones y Sistemas de EDO