EDO Notas de Curso

6
5
4
3
2
1
10
15
8 0
x1 (t) =
>
>
>
<y 0 (t) =
1
>
x02 (t) =
>
>
: 0
y2 (t) =
x1 (t)( a1 + b1 y1 (t)),
y1 (t)(c1 d1 x1 (t) e1 y1 (t)),
x2 (t)( a2 + b2 y2 (t)),
y2 (t)(c2 d2 x2 (t) kx1 (t)y1 (y)y2 (t)),
20
25
30
Notas de curso
ECUACIONES DIFERENCIALES
ORDINARIAS
Tcnicas de Resolucin
Luz Marina Moya y Edixon Rojas
3
20
10
Luz. M. Moya
[email protected]
Pontificia Universidad Javeriana
Bogot-Colombia
Edixon M. Rojas
[email protected]
Universidad Nacional
Bogot-Colombia
c 2014 Luz Marina Moya & Edixon Rojas

Copyright
N OTAS DE CURSO
Licencia Creative Commons 3.0 Reconocimiento-No Comercial-Compartir Igual
ndice general
Nociones Fundamentales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1
Definiciones y Ejemplos
1.2
Solucin de una Ecuacin Diferencial.
1.2.1
Solucin General y Particular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.2.2
Problema de Valor Inicial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
Ecuaciones Diferenciales de Primer Orden . . . . . . . . . . . . . 21
2.1
Ecuaciones Diferenciales No Lineales
2.1.1
Ecuaciones Separables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.1.2
Ecuaciones Exactas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.1.3
Factor Integrante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.2
Ecuaciones Diferenciales Lineales
36
2.3
Soluciones por Sustituciones
39
2.3.1
Ecuaciones Diferenciales Homogneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.3.2
Ecuaciones con Coeficientes Lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
2.3.3
Ecuacin de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.3.4
Ecuacin de Ricatti . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.3.5
Ecuacin de Lagrange . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
2.3.6
Reduccin a Separacin de Variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
6
11
22
Aplicaciones EDO de Primer Orden . . . . . . . . . . . . . . . . . . . . . . 57
3.1
Crecimiento y Decrecimiento
3.1.1
3.1.2
Modelo Logstico de Crecimiento Poblacional . . . . . . . . . . . . . . . . . . . . . . 61

Determinacin de Edades por el Mtodo del Carbono 14 . . . . . . . . . . . . 65
3.2
Mezclas
67
3.3
Problemas de Temperatura
72
3.4
Mecnica de Newton
80
3.5
Trayectorias Ortogonales
85
Ecuaciones Diferenciales Lineales de Orden Superior . 91
4.1
Teora Preliminar
4.2
Ecuaciones lineales Homogneas con Coeficientes Constantes 100
4.3
Mtodo de Coeficientes Indeterminados
103
4.4
Reduccin de Orden.
108
4.5
Variacin de Parmetros
112
4.6
Ecuacin de Cauchy-Euler
117
4.7
Aplicaciones a la Fsica
124
4.7.1
4.7.2
Movimiento Armnico Simple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124

Movimiento Armnico Amortiguado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
57
92
1 Nociones Fundamentales
En este captulo el lector encontrar el lenguaje bsico de las ecuaciones diferenciales

ordinarias, a decir; trminos como solucin general y particular (implcitas y explcitas),
problemas de volr inicial, as como una primera clasificacin de este tipo de ecuaciones.
Antes de entrar en detalles tcnicos, precisaremos el significado de algunos trminos que,
en ocasiones, se suponen de conocimiento general.
Definicin 1.1. Un objeto es una unidad elemental sobre la que se pueden hacer observaciones 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.
Los problemas que son de inters en reas de las ciencias e ingeniera, son problemas dinmicos, que involucran variables independientes que cambian respecto de otras variables
dependientes y es precisamente el concepto de derivada, los que permiten describir cuantitativamente las relaciones entre las variables involucradas en la descripcin del problema 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
Nociones Fundamentales
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.
1.1
Definiciones y Ejemplos
Definicin 1.2 (Ecuacin diferencial). Se llama ecuacin diferencial a toda ecuacin que involucre una funcin desconocida y alguna de sus derivadas. Ms exactamente, una ecuacin
diferencial es una expresin de la forma

f
f
n1 f f
n2 f
nt f
F x1 , . . . , xt , f ,
, . . . , n1 ,
, . . . , n2 ,
, . . . , nt = 0,
(1.1)
x1
xt
xt
x1 x2
x2
que involucre una funcin desconocida f en las variables x1 , . . . , xt ; algunas de sus derivadas
y donde F denota una funcin de n = n1 + + nt variables.
A continuacin presentaremos algunos ejemplos sencillos de ecuaciones diferenciales que
aparecen en algunas reas de las ciencias.
Ejemplo 1.1. Un modelo sencillo para el crecimiento de poblaciones es la Ley de Malthus,
la cual afirma que la tasa de crecimiento de la poblacin de un pas aumenta en forma
proporcional a la poblacin total P(t) del mismo en cualquier momento t. En otras palabras,
mientras ms personas haya en el momento t, ms personas habr en el futuro. En trminos
matemticos
dP
= kP,
dt
(1.2)
donde k es la constante de propocionalidad.

Ejemplo 1.2. La segunda ley de Newton, la cual se encarga de cuantificar el concepto
de fuerza, afirma que la fuerza neta F aplicada sobre una partcula de masa m es igual al
cambio en su momentum lineal con respecto al tiempo. En el caso en que la masa del cuerpo
permanece constante, se deduce que la fuerza aplicada sobre el cuerpo es proporcional a la
aceleracin a(t) que adquiere dicho cuerpo, esto es,
F = ma.
Si x(t) denota la posicin de la partcula en el tiempo t, entonces se tiene la siguiente
ecuacin diferencial,

dx
d2x
F t, x,
m 2 = 0,
(1.3)
dt
dt
donde la fuerza F es una funcin que depende del tiempo t, de la posicin x(t) y de la
velocidad dx/dt.
1.1 Definiciones y Ejemplos
Ejemplo 1.3. El modelo de Hodgkin-Huxley para pulsos neuronales. En neurobiologa, la

comunicacin neuronal es un rea de activa investigacin. En este caso, se pretende estudiar
las seales elctricas, tambin conocidas como pulsos o disparos, realizados por clulas
nerviosas individuales o neuronas. El encargado de propagar las seales elctricas hacia
el exterior de una neurona es el axn, tambin conocido como cilindroeje o neurita y el
cual es un tubo cilndrico largo que se extiende desde cada neurona. Los pulsos elctricos
aparecen dado que la membrana de cada neurona tiene preferencias de permeabilidad, esto
es, es fcilmente permeable por ciertos iones qumicos para los cuales dicha permeabilidad
es afectada por las corrientes y potenciales presentes.
Los elementos ms importantes en este sistema son los iones de sodio (Na+ ), los de potasio
(K + ); y la diferencia de potencial V existente entre el interior de la membrana y el medio
circundante, debido a las altas concentraciones de iones de K + al interior del axn. El
modelo de Hodgkin-Huxley propone la siguiente ecuacin para el potencial V ,
C
dV
= Ie Ii ,
dt
(1.4)
donde C es la capacitancia de la membrana, Ie la corriente externa e Ii la corriente interna.

La corriente interna, a su vez, satisface una ecuacin del siguiente estilo
Ii = g0 (V V0 ) + g1 m3 h(V V1 ) + g2 n4 (V V2 ),
(1.5)
donde g0 , g1 , g2 > 0 y V0 ,V1 ,V2 son constantes, mientras que m, n, h son variables de bloqueo,
correspondientes a la activacion del sodio (Na+ ), activacin de potasio (K + ) e inactivacin de Na+ respectivamente. El modelo es normalizado, de forma tal que las variables
de bloqueo siempre tomen valores entre 0 y 1 y tal que puedan ser interpretadas como
probabilidades para que cierto tipo de canal est abierto. Al combinar las dos ecuaciones
(1.4) y (1.5) obtenemos
C
dV
= Ie (g0 (V V0 ) + g1 m3 h(V V1 ) + g2 n4 (V V2 )).
dt
(1.6)
Mientras Ie es tratado como un parmetro externo, la dinmica interna depende en gran

parte del signo de los tres tres trminos que aparecen en (1.5).
Clasificacin de las Ecuaciones Diferenciales
1. Ecuacin Diferencial Ordinaria: Una ecuacin diferencial es llamada ecuacin diferencial ordinaria, abreviando EDO, si es una ecuacin diferencial que slo contiene
derivadas ordinarias. Ms exactamente, si es una ecuacin diferencial de la forma

df
dn f
F x, f , , , n = 0,
dx
dx
2. Ecuacin Diferencial Parcial: Una ecuacin diferencial parcial es una ecuacin
diferencial de la forma (1.1) en la cual aparecen las derivadas parciales de una funcin
que depende de ms de una variable.
8
Ejemplo 1.4. La ecuacin diferencial
dx
+ 4x = cos(t),
dt
donde x es la variable dependiente y t la variable independiente, es una ecuacin diferencial

ordinaria.
Ejemplo 1.5. Las ecuaciones diferenciales (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.
Ejemplo 1.6. La ecuacin diferencial
x x
+
= 3,
y t
donde x y t son variables independientes y y es la variable dependiente, es un ejemplo de
ecuacin diferencial parcial.
Ejemplo 1.7. Otro ejemplo de ecuacin diferencial parcial viene dado por
2
2u
u
2 u
=
c
p qu,
2
2
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.
Ahora que estamos en capacidad de distinguir las ecuaciones diferenciales de acuerdo a su
tipo, es momento de refinar un poco ms nuestra clasificacin. Para tal fin, introduciremos a
continuacin el concepto de orden y linealidad de una ecuacin diferencial.
Definicin 1.3 (Orden). El orden de una ecuacin diferencial es el orden de la mayor
derivada que aparecen en la ecuacin.
Ejemplo 1.9. Consideremos la ecuacin diferencial ordinaria
F(x, y, y0 , y00 , , y(n) ) = 0,
(1.7)
donde F es una funcin de n + 2 variables. Si en la expresin explicita de la ecuacin (1.7)

aparece y(n) , entonces nos encontramos frente a una ecuacin diferencial ordinaria de orden
n.
1.1 Definiciones y Ejemplos
Ejemplo 1.10. A continuacin explicaremos un poco el ejemplo anterior e indicaremos en

cada caso a que funcin corresponde F.
d2y
dy
1. La ecuacin diferencial dx
2 + 2 dx + 3y = 0 es una ecuacin diferencial de orden 2 en
la cual F est dada por la funcin F(x, y, z, w) = 3y + 2z + w. Ntese que si reemplazamos z por la primera derivada de y y w por la segunda derivada de y e igualamos a
cero, obtenemos la ecuacin diferencial dada.
2
t
2. De la misma forma, la ecuacin diferencial ddt 2x + x dx
dt cos(t)x = e es una ecuacin
diferencial ordinaria de orden 2, en la cual la funcin F corresponde a F(t, x, u, v) =
w + xu cos(t)x et .
3
3. La ecuacin diferencial ddt 3x = sin(x) es una ecuacin diferencial ordinaria de orden 3,

en la cual la funcin F corresponde a F(t, x, y, z, w) = w sin(x).
Definicin 1.4 (Grado). El grado de una ecuacin diferencial es el exponente al cual est
elevada la derivada de mayor orden que aparece en ella. Si esta derivada est elevada a un
exponente no natural el grado de la ecuacin no est definido.
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 grado
2.
Definicin 1.5. Una ecuacin diferencial ordinaria es lineal, si es de la forma
an (x)
dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y + g(x) = 0,
n1
n
n1
dx
dx
dx
(1.8)
donde los coeficientes g(x) y ai (x) con i = 1, . . . , n son funciones que no dependen de la
funcin incgina y y de ninguna de sus derivadas.
Si la ecuacin diferencial no puede llevarse a la forma (1.8) decimos que la ecuacin es no
lineal.
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.
10
Ejemplo 1.13. La ecuacin diferencial (y00 )2 = sin(x) + yex es una ecuacin diferencial
ordinaria cuadrtica de segundo orden.
d2 y
dy
= sin(x), es un ejemplo de ecuacin
Ejemplo 1.14. La ecuacin diferencial e dx2 + dx
diferencial de segundo orden no lineal. Ntese que este es un ejemplo en el cual no tiene
sentido hablar del grado de la ecuacin.
Otros ejemplos de ecuaciones diferenciales no lineales son los siguientes:

Ejemplo 1.15. La ecuacin diferencial parcial
u
u 3u
+
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.
Ejemplo 1.16. La ecuacn de Van der Pol, es una ecuacin diferencial ordinaria cuadrtica
de segundo orden
d2x
dx
(1 x2 ) + x = 0,
2
dt
dt
donde es una constante, usada para describir el comportamiento de un oscilador no
conservativo con amortiguamiento no lineal.
dy
Ejemplo 1.17. Las ecuaciones diferenciales dx
= xy2 y
de ecuaciones diferenciales ordinarias no lineales.
d3y
dx3
+ y2 = 0, tambin son ejemplos
Ejercicios
1. Clasifique las siguientes ecuaciones diferenciales de acuerdo a su tipo, indique las
variables dependiente e independiente y proporcione el orden de la ecuacin. Si la
ecuacin diferencial es ordinaria, determine 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
cos(x) + 4xy dx
+ y2 ex = 0.
2
e)
f)
g)
h)
3 xy3 z3y = 0.
ln(x) + 6xy00 4xy = cos(x).
d3y
dy 4
x dx
3 ( dx ) + y = 0.
(sin( ))y00 (cos( ))y0 = 2.
1.2 Solucin de una Ecuacin Diferencial.
11
Solucin
a) Esta ecuacin es ordinaria porque no contiene derivadas parciales. La variable dependiente es y y la independiente x. La ecuacin es de 3 orden ya que se tiene y000 como
la mayor derivada, y como es la derivada de mayor orden de exponente 1 el grado es
1. La ecuacin es lineal.
b) Esta ecuacin es ordinaria porque no contiene derivadas parciales. La variable dependiente es y y la independiente x. La ecuacin es de 2 orden ya que se tiene y00 , y como
es la derivada de mayor orden de exponente 1 el grado es 1. La ecuacin es lineal.
c) Esta ecuacin es ordinaria porque no contiene derivadas parciales. La variable dependiente es y y la independiente x. La ecuacin es de 2 orden ya que se tiene y00 , y como
es la derivada de mayor orden de exponente 1 el grado es 1. La ecuacin es no lineal
debido al trmino y2 .
d) Esta ecuacin es ordinaria porque no contiene derivadas parciales. La variable dependiente es y y la independiente x. La ecuacin es de 2 orden ya que se tiene d 2 y/dx2 , y
como es la derivada de mayor orden de exponente 3 el grado es 3. La ecuacin es no
lineal.
e) Esta ecuacin es parcial de tercer orden porque contiene tercera derivada parcial. La
variable dependiente es y y la independiente x y z.
1.2
Solucin de una Ecuacin Diferencial.

Uno de los objetivos al estudiar una ecuacin diferencial es determinar si sta posee solucin
y, en caso de que tal de que exista, si es nica. A continuacin daremos la definicin precisa
de lo que entenderemos por solucin de una ecuacin diferencial ordinaria.
Definicin 1.6 (Solucin). Una solucin de la ecuacin diferencial ordinaria
F(x, y, y0 , y00 , , y(n) ) = 0,
(1.9)
es una funcin y = f (x) definida en un intervalo abierto I, tal que para todo x I las derivadas
f 0 (x), f 00 (x), . . . , f (n) (x) existen y al reemplazarlas en la ecuacin (1.9), la expresin se
convierte en una identidad.
Ejemplo 1.18. Consideremos la ecuacin diferencial
dP
= kP.
dt
(1.10)
Verificar que la funcin

P(t) = cekt ,
< t < ,
es solucin de la ecuacin, donde c es una constante arbitraria.
12
Solucin. Derivando la funcin dada obtenemos P0 (t) = ckekt = kP(t), que es expresin en
la ecuacin (1.10).

En el proceso de solucionar una ecuacin diferencial, en algn momento necesitamos de la
integral indefinida (o antiderivada) de alguna funcin. Por ejemplo, para la ecuacin
dy
= g(x)
dx
la solucin tiene la forma
y(x) =
R
g(x)dx + C
en algn intervalo I, donde g(x)dx es cualquier antiderivada de g(x).

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.
Ejemplo 1.19. Comprobar que la frmula y =

x + yy0 = 0,
p
(1 + x2 ) satisface la ecuacin diferencial
(1.11)
sin embargo no es una solucin de esta.

p
Solucin. Derivando y = (1 + x2 ) tenemos que y0 =
x
.
(1+x2 )
Sustituyendo en la
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.20. La ecuacin diferencial de primer orden

xy0 = 1
(1.12)
no tiene solucin en el intervalo (1, 1).

Solucin. Notemos que solo las funciones de la forma y = ln(|x|) + c con c R satisfacen
la identidad (1.12). Sin embargo, esta funcin es discontinua en x = 0 (por tanto no es
diferenciable en este punto) y por definicin de solucin esta debe satisfacer la ecuacin en
todo punto del intervalo (1, 1).
13
Definicin 1.7 (Solucin Explcita). Una solucin explcita de la ecuacin diferencial (1.9)
es una solucin y = f (x), en donde y est escrita explcitamente en trminos de x, esto es, y
se ha despejado en trminos de x.
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), obtenemos f 0 (x) = 1 ex . Luego
f 0 (x) + f (x) = 1 ex + x + ex = 1 + x.

Ejemplo 1.22. Verifique que para cualquier par de nmeros reales c1 y c2 , la funcin
(x) = c1 e2x + c2 e3x
definida en todo R, es una solucin explcita de la ecuacin
y00 5y0 + 6y = 0.
Solucin. Dado que,
0 (x) = 2c1 e2x + 3c2 e3x y 00 (x) = 4c1 e2x + 9c2 e3x ,
entonces
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
tan2 (x) = (x + tan(x) x) = tan2 (x).
Se tiene as la igualdad para cualquier x en cada uno de los intervalos en donde y est
definida.
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.
14
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.
2. La funcin (x) es diferenciable en el intervalo I y satisface
F(x, (x), 0 (x), . . . , (n) (x)) = 0,
para todo x I.
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 (x) definida por la relacin g(x, y) = 0
es la solucin.
Ejemplo 1.24. Sea C > 0 un nmero real. Verificar que la relacin

x2 + y2 = C
(1.13)
determina una solucin implcita de la ecuacin diferencial

dy
+x = 0
dx

en el intervalo ( C, C).
(1.14)
Solucin: Diferenciando implcitamente la ecuacin (1.13) obtenemos

2x + 2y
dy
= 0,
dx
de donde, dividiendo todo entre 2, obtenemos

la igualdad (1.14). Ahora, despejando
y de
2
la ecuacin (1.13) obtenemos y = C x . Consideremos la funcin (x) = C x2

15

2
definida
en el intervalo ( C, C). Puesto que d /dx = x/ C x para todo x
( C, C), sustituyendo en (1.14) obtenemos
p
x
C x2
+ x = 0.
C x2
Luego, la ecuacin (1.13) determina una solucin implcita de la ecuacin (1.14).

Ejemplo 1.25. Verificar que la relacin
x2 y + y2 C = 0,
(1.15)
determina una solucin implcita de la ecuacin diferencial

2xy + (x2 + 2y)
dy
= 0.
dx
(1.16)
Solucin: El Teorema de la Funcin Implcita garantiza la existencia de una funcin

diferenciable y = f (x) que satisface (1.15). Diferenciando implcitamente la ecuacin (1.15)
se obtiene
dy
dy
+ 2y = 0,
dx
dx
dy
2xy + (x2 + 2y) = 0,
dx
2xy + x2
que es igual a la ecuacin (1.16). Luego la relacin (1.15) determina una solucin implcita
de (1.16).
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.
16
1.2.1
Solucin General y Particular

Definicin 1.9 (Solucin General). La solucin general de una ecuacin diferencial es el
conjunto de todas las funciones que verifican la ecuacin diferencial.
Aqu diremos que la solucin general de una ecuacin diferencial consiste en una familia
n-paramtrica de funciones, parmetros que seran las constantes que aparecen al realizar
la operacin de integrar n-veces, siendo n el orden de la ecuacin. Cuando existe alguna
(inusual) solucin que no pertenence a dicha familia, entonces esta funcin, como es tpico,
recibir el nombre de solucin singular.
Definicin 1.10 (Solucin Particular). Se llama solucin particular de la ecuacin diferencial, a cualquier funcin que la satisfaga. Esto es, a cualquier elemento del conjunto solucin
general. Una solucin particular se puede obtener fijando valores a los parmetros de la
familia de funciones solucin de la ecuacin.
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 la expresin para y obtenemos y0 = 4kxe2x , es decir, y0 = 4xy.
Obs En el ejemplo anterior, la familia de soluciones dada es en realidad la solucin general

de la ecuacin diferencial en cuestin; pero para demostrar esto necesitamos de las
tcnicas bsicas para hallar soluciones de ecuaciones diferenciales ordinarias, las
cuales pospondremos hasta el prximo captulo.
Ejemplo 1.27. Demostrar que la funcin y = 3 3ex

ecuacin diferencial
dy
+ xy = 3x.
dx
Solucin: Derivando y obtenemos y0 = 3xex
2 /2
2 /2
es una solucin particular de la
. Luego,
x2
x2
x2
x2
dy
+ xy = 3xe 2 + x(3 3e 2 ) = 3xe 2 + 3x 3xe 2 = 3x.
dx
17
Ejemplo 1.28. Demostrar que la familia de funciones y = cx + c2 , definidas en toda la recta

real, hacen parte de la solucin general de la ecuacin diferencial
y = xy0 + (y0 )2
2
y que la funcin y1 = x4 es una solucin singular.

Solucin: Derivando la expresin para y obtenemos y0 = c, de donde se tiene la identidad
deseada, y as la familia de funciones y = cx + c2 forman parte de la solucin general. Por
otro lado, derivando y1 obtenemos y01 = 2x y sustituyendo en la ecuacin vemos que esta se
2
satisface. Sin embargo, no es posible obtener y1 = x4 como una solucin particular de la
2
ecuacin, no importa el valor de c en y = cx + c2 . As y1 = x4 es una solucin singular.
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
esta ecuacin. Sin embargo y = 0 tambin es una solucin la cual no se obtiene de
1
la familia de funciones de la forma y = (x+c)
2 y por lo tanto es una solucin singular.
2
c
Pero, por otro lado, cualquier funcin de la forma y = (cx+1)
2 tambin es solucin de la
ecuacin diferencial y la solucin y = 0 es parte de esta familia, por lo cual y = 0 es una
solucin tanto singular como no singular dependiendo de la familia de soluciones que
se considere. Esto lleva a cuestionarse la exactitud de estos trminos. Sin embargo, en
este texto nos abocaremos al estudio de algunas tcnicas de resolucin de ecuaciones
diferenciales que nos permitir establecer una nica familia de funciones, las cuales
llamaremos soluciones generales, y no consideraremos soluciones singulares o
familias de soluciones con otra configuracin en su frmula.
1.2.2
Problema de Valor Inicial

Definicin 1.11 (Problema de Valor Inicial). Una ecuacin diferencial ordinaria de orden n
con condiciones iniciales, es una ecuacin diferencial de la forma

dy
d n1 y d n y
F x, y, , . . . , n1 , n = 0,
dx
dx
dx
para la cual se busca una solucin definida en un intervalo I sujeta a n condiciones
y(x0 ) = y0 ,
dy
d n1 y
(x0 ) = y1 , . . . , n1 (x0 ) = yn1 ,
dx
dx
en un punto x0 I y donde y0 , y1 , . . . , yn1 R son constantes dadas.

El conjunto de datos que consta de la ecuacin diferencial y las n condiciones, recibe el
nombre de Problema de Valor Inicial (PVI) o tambin Problema de Cauchy.
18
Ejemplo 1.29. La ecuacin
4
dy
+ 4y ex = 0, sujeta a las condiciones y(0) = ,
dx
3
es un problema con condiciones iniciales.
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,
(xex + cx)
dy y
= 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,
por lo que c = 1. As, la solucin del problema con valor inicial es

y = xex x.
Figura 1.1: Algunas curvas solucin en la familia y = xex + cx.

19
Ejercicios
1. Compruebe 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)].
3x .
f ) y00 + y0 12y = 0; y = c1 e4x + c2 e

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 .

9 cos(2x)7 sin(2x)
1
2
000
00
2x
x
x
e2x .
j) y y e sin (x) = 0; y = c1 + c2 e + c3 e + 12 +
520
k)
l)
m)
n)
)
x3 y000 + 2x2 y00 xy0 + y = 12x2 ; y = c1 x + c2 x ln |x| + 4x2 , x > 0.

y2 1 (2y + xy)y0 = 0; y2 1 = (x + 2)2 .
dy
exy + eyx dx
= 0; e2y + e2x = 1.
xy0 = 2y; y = x2 .
y00 = 1 2 ; y = arcsin(x) + 1 x2 .
1x
o) y0 = 13 y; y = 2ex/3 .
dr
2r sin( ) = 0; r = a sec2 ( ).
p) cos( ) d
2. Demuestre que la ecuacin diferencial

dy
+ |y| + 3 = 0
dx
no tiene soluciones.
Solucin
1. En los ejercicios c), d), g) e i), la funcin indicada no es solucin de la ecuacin
diferencial dada, en los dems casos la funcin indicada si satisface la ecuacin
diferencial.

dy
2. Toda funcin en valor absoluto es positiva o 0 y por tanto la identidad dx
+|y|+3 = 0
no puede ser satisfecha.
20
15
10
5
20
40
60
80
2 Ecuaciones Diferenciales de Primer Orden
Como vimos, una ecuacin de diferencial de primer orden es una relacin de la forma
y0 (t) = f (t, y).
Cuando esta relacin es lineal en las variables t, y y y0 ; es decir, que se pueda reescribir
como
a1 (t)y0 + a2 (t)y + g(t) = 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 soluciones de ciertos tipos de ecuaciones
de primer orden. De manera ms precisa, se estudiarn tcnicas de solucin de ecuaciones
de primer orden no lineales llamadas ecuaciones separables 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 tales como ecuaciones Homogneas, de Bernoulli,
de Ricatti, con coeficientes 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 < y < t0 + h contenido en
a < t < b existe una nica solucin y = (t) del problema de valor inicial
y0 = f (t, y),
y(t0 ) = y0 .
100
Ecuaciones Diferenciales de Primer Orden
22
2.1
Ecuaciones Diferenciales No Lineales

Comenzaremos el estudio de las ecuaciones diferenciales no lineales de primer orden
considerando la clase de ecuaciones llamada separables.
2.1.1
Ecuaciones Separables
La ecuacin diferencial de primer orden
dy
= r(x, y)
dx
se denomina Separable si la funcin r(x, y) se puede expresar de la forma r(x, y) = f (x)
g(y), donde f es una funcin en trminos de x y g es una funcin en y:
dy
= f (x) g(y)
dx
a veces por comodidad
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, es decir
Z
dy
p(y) dx =
dx
p(y)dy =
f (x)dx
f (x)dx.
Ejemplo 2.1. Resolver la ecuacin

dy
= x + 2.
dx
Solucin. Si integramos ambos lados de la ecuacin con respecto a la variable independiente
x, es decir
Z
Z
dy
dx = (x + 2)dx
dx
o,
dy =
(x + 2)dx.
2.1 Ecuaciones Diferenciales No Lineales
23
Evaluando las integrales se obtiene

y=
x2
+ 2x + c
2
donde c es una constante arbitraria.

La ltima expresin constituye una familia de curvas parablicas, cuyo grfico es la Figura
2.1
Figura 2.1: Familia de curvas de y =
x2
2
+ 2x + c.
Ejemplo 2.2. Resolver la ecuacin

dy sec2 (y)
=
.
dx
1 + x2
Solucin. Si reescribimos la ecuacin en la forma
1 dy
1
=
sec2 (y) dx 1 + x2
e integramos ambos lados de la ecuacin con respecto a la variable independiente x, es decir
Z
1 dy
dx =
sec2 (y) dx
Z
dy
=
sec2 (y)
1
dx
1 + x2
dx
.
1 + x2
24
Evaluando las integrales se obtiene
y sin(2y)
+
= arctan(x) + c
2
4
2y + sin(2y) =4 arctan(x) + c
donde c es una constante arbitraria.

Ejemplo 2.3. Resolver el problema de valor inicial
dy
= x3 (1 y),
dx
y(0) = 3.
Solucin. Al separar e integrar tenemos

dy
=x3 dx
1y
Z
Z
1 dy
dx = x3 dx
y 1 dx
x4
ln(y 1) = + c.
4
Despejamos a ambos lados y tenemos:
x4
y 1 = e 4 ec .
x4
Despejando y se tiene y(x) = 1 + e 4 k, donde por simplicidad en la presentacin hemos

definido k = ec .
Empleando la condicin inicial y(0) = 3 obtenemos k = 2, de este modo la solucin es:
x4
y(x) = 1 + 2e 4 .

Ejercicios
1. Encontrar La solucin general de las siguientes ecuaciones diferenciales, por variables
separables:
dy
a) dx
= xy .
b) (1 y2 )dx xydy = 0.
c) (y + 1)dx + (y 1)(1 + x2 )dy.
25
d) (2x + 1)y0 + y2 = 0.
e) y0 x3 = x3 y.
f ) xy2 dy
+ (x2 + 1)dx.
2
g) xy + 1 + x2 y0 = 0.
h) xy2 dx + xdy = 0.
dy
i) dx
= 1 + x + y + xy.
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.
b)
c)
d)
e)
dy
y2 1
dx = x2 1 ; y(2) = 2.
dy
x2 dx
= y xy; y(1) = 1.
dy
2 2
dx = 3x (y + 1); y(0) = 1.
dy
2y dx
= 2x ; y(5) = 2.
x 16

dy
tan(x) dx
= y; y 2 = 2 .
dy
2
2 2
dx = 2xy + 3x y ; y(1) = 1.
f)
g)
3. Hallar todas las funciones f (x) tales que su derivada es el cuadrado de la funcin.
Solucin
1.
a)
b)
c)
d)
y = cx.
c = x2 (1 y2 ).
tan1 (x) + c = y + 2 ln |y + 1|.
2
y = ln |2x+1|+c
.
4
e) y = cex /4 1.
f) y3 = 23 x2 3 ln |x| + c.
g) y = 1 2 .
h) y =
2.
1+x +c
1
x+c .
ln |1 + y| = x + x2 + c.
2
y = tan( x 1+cx
).
x
4 3
y = 3 x + 4x + 1.
y = x.
1/x
y = e ex .
y = tan( 4 + x3 ).
x2 16(x2 16+ x2 16)
e) y =
.
2
(i)
j)
a)
b)
c)
d)
x 16
f) y = 21 sin(x).
g) y = 1+x12 +x3 .
1
3. f (x) = x+c
, y 6= 0.
26
2.1.2
Ecuaciones Exactas
Una ecuacin diferencial de la forma
M(x, y)dx + N(x, y)dy = 0
(2.1)
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:
dF(x, y) = M(x, y)dx + N(x, y)dy.
En este caso se tendr que
F
= M(x, y),
x
F
= N(x, y).
y
Solucin de una ecuacin diferencial exacta
Si la ecuacin M(x, y)dx + N(x, y)dy = 0 es exacta entonces podemos encontrar F(x, y) tal
que:
dF = M(x, y)dx + N(x, y)dy = 0
es decir F(x, y) = k, con k R.
De esta forma una vez encontrada F, la solucin de la ecuacin es F(x, y) = k.
El siguiente resultado nos permite determinar si una ecuacin diferencial de la forma (2.1)
es exacta. Daremos su prueba pues en ella se describe el procedimiento para resolver este
tipo de ecuaciones.
Teorema 2.2 Criterio para exactitud de una ecuacin diferencial. Sean M(x, y)
y N(x, y) funciones continuas y con derivadas parciales de primer orden continuas en una
regin rectangular R definida por a < x < b, c < y < d. La ecuacin diferencial M(x, y)dx +
N(x, y)dy = 0 es exacta si y solamente si:
M N
=
.
y
x
Demostracin: La demostracin es en dos partes.
1. Supongamos que la ecuacin es exacta: entonces existir una funcin F(x, y) tal que:
F
= M(x, y),
x
F
= N(x, y).
y
Derivando en estas igualdades respecto de x y y respectivamente, se tiene:

M
2F
=
,
y x
y
27
2F
N
=
.
x y
x
Puesto que M/ y y N/ x son continuas, entonces 2 F/ x y y 2 F/ y x son

tambin lo son. Esto garantiza la igualdad
M N
=
.
y
x
(2.2)
2. Ahora mostramos que si M y N satisfacen la ecuacin (2.2) entonces la ecuacin

M(x, y)dx + N(x, y)dy = 0 es exacta. Para la demostracin se quiere determinar la
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
F(x, y) =
M(x, y)dx + g(y)
(2.3)
en donde la funcin arbitraria g(y) es la constante de integracin. Ahora debemos

determinar g(y). Para ello utilizaremos (2.3) y lo derivamos con respecto a y

Z
F
dg(y)
M(x, y)dx
=
+
.
y
y
dy
Como F/ y = N(x, y) y despejando g0 (y) se tiene
g (y) = N(x, y)
Z

M(x, y)dx
.
y
(2.4)
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)
R

Z
Z
( M(x, y)dx)
F(x, y) = M(x, y)dx +
N(x, y)
dy.
y
28

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)
siguiendo un procedimiento anlogo al expuesto anteriormente, pero intercambiando variables.

Para solucionar una ecuacin diferencial exacta se seguirn los pasos indicados en la
demostracin anterior, donde se ha establecido cmo encontrar la funcin F(x, y).
Obs Note que la ecuacin de variables separables
f (x)
dy
=
,
dx
p(y)
se puede escribir en forma equivalente como
p(y)dy f (x)dx = 0
la cual es exacta, como se ve fcilmente.

(ex y + xex y)dx + (xex + 2)dy = 0,
y(0) = 1.
(2.5)
Solucin. En este caso M(x, y) = ex y + xex y, y N(x, y) = xex + 2. Como

N
M
= 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)
Derivando parcialmente con respecto a y y sustituimos xex + 2 en vez de N(x, y), ya que
F/ y = N(x, y):
F
=xex + g0 (y)
y
xex + 2 =xex + g0 (y)
2 =g0 (y).
29
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
F(x, y) = xex y + 2y.
Es decir, la solucin de la ecuacin diferencial (2.5) est dada por
xex y + 2y = k.
Alternativamente, se puede integrar simultneamente las dos funciones
M(x, y) = ex y + xex y,
N(x, y) = xex + 2
y compararlas para determinar unvocamente F(x, y), salvo una constante. As

F(x, y) =
M(x, y)dx =
N(x, y)dy =
(ex y + xex y)dx
(ex x + 2)dy
=ex y + (xex ex )y = xyex ,

por otro lado,
F(x, y) =
=(ex x + 2)y = xyex + 2y.

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, en esta ecuacin, cuando x = 0, y = 1 obtenemos k = 2. De
esta manera la solucin es:
2
y= x
.
xe + 2

(y sin(xy) + x2 )dx + (x sin(xy) + y2 )dy = 0,
y(0) = 1.
Solucin. En este caso M(x, y) = y sin(xy) + x2 , y N(x, y) = x sin(xy) + y2 , calculamos

M
N
y
y
x
30
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 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
g(y) =
Por tanto,
y3
y dy = .
3
2
F(x, y) = cos(xy) +
x3 y3
+
3
3
es decir
x3 y3
+ .
3
3
Ahora usemos la condicin inicial y(0) = 1. Tenemos
k = cos(xy) +
k = cos(0(1)) +
k =1
03 (1)3
+
3
3
1
4
= .
3
3
As la solucin del problema de valor inicial es

4
x3 y3
= cos(xy) + + .
3
3
3
31
Ejercicios
1. Determine 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.
y
1
d) (y2 sin(x))dx +
x
x dy = 0.
e) (tan |y| 2)dx + x sec2 (y) 1y dy = 0.

3y2
2y
x2
1 dy = 0.
f ) 2x
dx
+
+
4
3
2
y
y
x
x
y
g) (ex sin(y) + tan(y))dx+ (ex cos(y)

+ x sec2 (y))dy = 0.

h) (cos(x) + ln(y))dx + xy + ey dy = 0.
i) (1 + yexy )dx + (2y + xexy )dy = 0.
2. Muestre que cualquier ecuacin de la forma
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. Determine la funcin ms general
M(x, y) de modo que la ecuacin sea exacta:
a) M(x, y)dx + sec2 (y) yx2 dy = 0.
b) M(x, y)dx + (x3 y3 )dy = 0.
c) M(x, y)dx + (e2x ey sin(x))dy = 0.
5. Determine 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
(Ax + By)dx + (Cx + Dy)dy = 0
sea exacta.
7. Muestre que la ecuacin y0 + y = 0 es exacta si la multiplicamos por ex y solucionela.
Solucin
1. a)
b)
c)
d)
e)
y = xk2 3x .
k = y2 x x2 .
x2 y2 ln |x| + sin(y) = c.
No es exacta.
c = x tan(y) 2x ln |y|.
32
3.
f)
g)
h)
i)
a)
b)
c)
k = xy + xy3 .
k = ex sin(y) + x tan(y).
k = ex sin(y) + x tan(y) + ey .
k = x + exy + y2 .
1 = x4 y2 2x3 y x2 3x.
0 = 4y sin(x) 2 cos2 (x) + tan(x) + 2y2 .
1 = y3 ex xex + ex + 34 y4 .
2
4.
5.
6.
7.
2.1.3
d) 21 = cos(x) + y cos(x) + y2 .
a) 1y .
b) 3x2 y.
c) 2e2x y.
a) x3 y2 32 x2 + 4xy.
2
b) cos(x) + x2 cos(y).
B = C, A y D cualquier nmero real.
y = cex .
Factor Integrante
Algunas ecuaciones diferenciales que no son exactas pueden ser transformadas en estas
mediante una funcin que multiplicada a la ecuacin original la transforma en una ecuacin
diferencial exacta. Estas funciones reciben el nombre de factores integrantes.
Mtodo para hallar factores integrantes especiales
Sea
M(x, y)dx + N(x, y)dy = 0
(2.7)
una ecuacin diferencial que no es exacta. Si multiplicamos esta ecuacin por una funcin
(x, y) adecuada que haga que la ecuacin resultante:
(x, y)M(x, y)dx + (x, y)N(x, y)dy = 0
sea exacta, se denominar a esta funcn factor de integracin para la ecuacin. Para determinar 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 esto se reduce a la ecuacin

N M
M
N
=
.
y
x
x
y
(2.8)
Por facilidad, consideraremos que depende de una sola variable. En este caso, suponemos
que depende de x, es decir = (x), de tal forma que la ecuacin (2.8) se reduce a la
ecuacin separable
=
x
donde
" M
M
y
Nx
33
#
Nx
N
depende solo de x. Entonces el factor integrante para (2.7) es:
!
Z M N
y
x
(x) = exp
dx .
N
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 el factor integrante para (2.7) es:
!
Z N M
x
y
(y) = exp
dy .
M
Ejemplo 2.6. Resolver la ecuacin diferencial
(2y2 + 2y + 4x2 )dx + (2xy + x)dy = 0.
(2.9)
Solucin. Ntese que la ecuacin diferencial no es de variables separables. Determinamos

si es exacta:
M
= 4y + 2
y
N
= 2y + 1,
x
como
M N
6=
y
x
entonces la ecuacin diferencial no es exacta. Calculamos

M
y
Nx
4y + 2 2y 1
2y + 1
1
=
= ,
2xy + x
x(2y + 1) x
34
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
!
Z
1
(x) = exp
dx = x.
x
Al multiplicar la ecuacin (2.9) por (x) = x, obtenemos la ecuacin exacta
(2xy2 + 2xy + 4x3 )dx + (2x2 y + x2 )dy = 0.
En este caso M(x, y) = 2xy2 + 2xy + 4x3 , y N(x, y) = 2x2 y + x2 . Calculamos
N
M
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).
Ahora derivamos parcialmente F con respecto a y y sustituimos 2x2 y + x2 en lugar de N(x, y)

ya que F/ y = N(x, y):
F
=2x2 y + x2 + g0 (y)
y
2x2 y + x2 =2x2 y + x2 + g0 (y)
0 =g0 (y).
As tenemos que
g(y) = c.
Por tanto,
F(x, y) = x2 y2 + x2 y + x4
35
es decir:
k = x2 y2 + x2 y + x4 .

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. Muestre 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 ad bc 6= 0, y sean m y n nmeros
reales arbitrarios. Demuestre que
(axm y + byn+1 )dx + (cxm+1 + dxyn )dy = 0
tiene un factor integrante de la forma (x, y) = x y .
5. Determine el factor integrante (x, y) = x y y resolver las siguientes ecuaciones
diferenciales:
1
a) ( yx2 xy
)dx y12 dy = 0.
b) (12 + 5xy)dx + (6x/y + 2x2 )dy = 0.
2
c) 2 xy dx + ( xy2 1)dy = 0.
Solucin
1. a)
b)
c)
d)
e)
(x) = x2/3 , y = 21 ex + cex .

2
(x) = x2 , y = x2 1 + cex .
(x) = x, ex (x 1)y = x + c.
(y) = y2 , c = x + x2 y1 .
(x) = xex , y = ex (x2c+x3 ) .
f ) (y) = 1y , x ln |y| + ex sin(y) = c.
36
3.
b)
c)
5.
2.2
1
, 1
xy 3y = c.
x2 y2
(xy) = x31y3 , 2x1
2 y2 3 ln |y| = c.
1
(xy) = xy(xy+1) , xy ln |y| = c.
(x, y) = xy2 , xy 31 x3 = c.
(y) = x7 y3 , 2x6 y3 x5 y2
a) (xy) =
a)
b)
c) (x, y) = x2 y2 ,
x2
y
= c.
+ y = c.
Ecuaciones Diferenciales Lineales

Una Ecuacin Diferencial Lineal de primer orden es una ecuacin de la forma:
a1 (x)
dy
+ a0 (x)y = g(x), a1 (x) 6= 0
dx
(2.10)
por lo tanto
dy a0 (x)
g(x)
+
y=
.
dx a1 (x)
a1 (x)
Si definimos
P(x) =
g(x)
a0 (x)
y Q(x) =
a1 (x)
a1 (x)
entonces la ecuacin diferencial lineal (2.10) toma la forma

dy
+ P(x)y = Q(x).
dx
(2.11)
Notemos que la ecuacin (2.11) no es de variables separables ni exacta. Ahora se trata de

encontrar un factor integrante (x). Asumimos que las funciones P y Q son continuas en
algn intervalo I. Luego multiplicamos por (x) la ecuacin (2.11) y tenemos
dy
(x) + (x)P(x)y = (x)Q(x)
dx

(x)P(x)y (x)Q(x) dx + (x) dy = 0.
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)
(2.12)
2.2 Ecuaciones Diferenciales Lineales
37
Al integrar ambos lados tenemos el factor integrante

R
(x) = e
P(x)dx
(2.13)
Como se acaba de ver
(x)P(x) = 0 (x)
con esto la ecuacin diferencial (2.12) se convierte en:

dy d(x)
y = (x)Q(x)
(x) +
dx dx
d (x)y
= (x)Q(x)
dx
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 (2.14) es la solucin general de la ecuacin lineal (2.11).
En resumen, la solucin de la ecuacin lineal (estndar)
dy
+ P(x)y = Q(x)
dx
siendo P y Q continuas en algn intervalo I R, est dada por

Z
1
(x)Q(x)dx + c .
y(x) =
(x)
Donde
(x) = e
P(x)dx
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).
Ejemplo 2.7. Resolver
x
dy
1
+ 2y = 3 ,
dx
x
x 6= 0.
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
38
Donde vemos que
2
1
y Q(x) = 4 ,
x
x
entonces podemos hallar el factor integrante sustituyendo en la formula (2.13)
P(x) =
R 2
x dx
(x) = e
= e2 ln(x) = x2
y reemplazando en la ecuacin (2.14) tenemos

"Z
#
1
1
x2 4 dx + c
y(x) = 2
x
x
"Z
#
1
1
y(x) = 2
dx + c
x
x2
y(x) =
c
1
+ 2.
3
x
x
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.
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) Muestre por sustitucin que y1 + cy0 es una solucin de y0 + p(x)y = q(x) para
todo nmero c.
b) Muestre 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) muestre que y0 es solucin de de y0 + p(x)y = 0.
b) Escriba toda solucin de y0 + p(x)y = q(x) en trminos de y1 y y2 .
2.3 Soluciones por Sustituciones
39
5. Verifique que las funciones dadas para 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 .
Solucin
1. a)
b)
c)
d)
e)
f)
2. a)
4.
5.
2.3
y(x) = 3e2x x.
y(x) = x + x42 .
C
.
y(x) = 3x + x1/3

y(x) = 21 x2 + 19
2 x.
y(x) = 1.
1
y(x) = 1 + e 2 x(x+2) .
2 +C
x(y) = (1/2)y
.
y4
y
1 2
b) x(y) = 2 y + c e .

y
1
arctan(y)
+
c
(1 + y2 ).
c) x(y) = 2(1+y
+
2)
2
b)
a)
b)
c)
2
y = y1 +y
2 + c(y1 y2 ).
y = 1 + cex .
y = x + c/x.
y = c + c(x2 x).
Soluciones por Sustituciones

Para resolver una ecuacin diferencial, primero identificamos qu tipo de ecuacin es. En
algunos casos es posible hacer una sustitucin para transformar la ecuacin en una de un
tipo conocido.
2.3.1
Ecuaciones Diferenciales Homogneas

Se puede cambiar la forma de una ecuacin diferencial de la siguiente manera. Sea y
cualquier solucin de la ecuacin
y0 = F(x, y)
(2.15)
se introduce una nueva funcin (desconocida) u tal que y(x) = g(u(x)), y u debe satisfacer
la nueva ecuacin
g0 (u)u0 = F(x, g(u)).
(2.16)
A la inversa, si u satisface (2.16), entonces satisface (2.15). Este proceso tiene sentido si
la ecuacin resultante es ms sencilla de resolver que la inicial. La sustitucin que se va a
utilizar es y = ux, y convierte la ecuacin en variables separables.
40
Una funcin F de dos variables F(x, y) es homognea de grado n si y slo si F(tx,ty) =

t n F(x, y) para todo nmero x, y,t. Si la funcin F(x, y) es de grado 0; esto es, F(tx,ty) =
F(x, y) para t = 1y , asumiendo y 6= 0, se tendr F(x, y) = F( xy , 1), es decir, la funcin puede
asumirse como una funcin que depende solamente de la variable xy .
La ecuacin diferencial
dy
= r(x, y)
dx
(2.17)
es Homognea si r(x, y) es una funcin que solo depende de xy . En este caso la ecuacin
(2.17) toma la forma
y
dy
= r 1, .
(2.18)
dx
x
Ejemplo 2.8. Determinar si la ecuacin
(y3 xy2 )dx + 2x2 ydy = 0
es homognea.
Solucin. Escribimos la ecuacin en la forma
dy
dx
dy
dx
dy
dx
dy
dx
xy2 y3
2x2 y
xy2
y3
= 2 2
2x y 2x y
y
y2
= 2
2x 2x
1 y 1 y 2
=
.
2x 2 x
=
Vemos que el lado derecho solo depende de
y
x
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
u+x
du
= r(1, u).
dx
La nueva ecuacin es de variables separables x y u, una vez resuelta reemplazamos u por

para dejar la solucin en trminos de las variables originales.
Ejemplo 2.9. Hallar la solucin general de la ecuacin diferencial
(y3 xy2 )dx + 2x2 ydy = 0.
y
x
41
Solucin. Como vimos en el ejemplo anterior esta ecuacin es homognea, ahora hacemos
dy
la sustitucin u = xy y dx
= u + x du
dx .
As:

dy 1 y 1 y 2
=
dx 2 x 2 x
1
du 1
u + x = u u2 .
dx 2
2
La ecuacin anterior es separable. Por tanto obtenemos,
Z
2
dx
du =
2
u+u
x

2 ln(u) ln(u + 1) = ln(x) + c.
Ahora reemplazando u =
y
x
tenemos

y+x
y
+ 2 ln
= ln(x) + c
2 ln
x
x
la solucin en forma implcita de la ecuacin.
Ejercicios
1. En las siguientes funciones, determine si la funcin dada es homognea. Si lo es,
indique su grado de homogeneidad.
4
a) f (x, y) = 3x2 y x3 + 2yx .
2x2 y+xy2
b) f (x, y) =
.
y
3
2 2
y
c) f (x, y) = 4xy +2x
.
xy
d) f (x, y) = (x + 3y 2)2 .
2. Resolver las siguientes ecuaciones diferenciales.
a) y0 = y+x
x .
2yx
0
b) y = y .
y
c) y0 = x+y
.
x2 +y2
xy .
x
.
x+y
2x2 +3y2 +y
y0 =
.
x
2
2
(y + xy + x )dx + (2x2 + 3xy + y2 )dy = 0.
d) y0 =
e) y0 =
f)
g)
h) (2x + y)dx + (3x + 2y)dy = 0; y(1) = 0.
3. Encontrar La solucin general de las siguientes ecuaciones diferenciales homogneas.
dy
x+y
a) dx
= xy
.
42
b)
dy
dx
dy
dx
y
y2x .
x+2y
y+2x .
c)
=
4. Demuestre que las curvas solucin de la ecuacin homognea
y(2x3 y3 )
dy
=
dx
x(y3 x3 )
son de la forma x3 + y3 = Cxy.

5. Muestre que si y0 = F(x, y) es una ecuacin homognea, y y = g(x) es una solucin,
entonces y = (1/a)g(ax) es una solucin para cualquier a. Interprete su resultado
geomtricamente.
6. Comprobar que las siguientes ecuaciones diferenciales son homogneas y hallar la
solucin general
a) 2xydxp+ (x2 + y2 )dy = 0.
b) (x + y2 xy)dy
ydx = 0.
y
0
c) xy y x sin x = 0.
d) (2x2 y +y3 )dx + (xy2 2x3 )dy =
0.
e) xy cos xy dx xy sin xy + cos xy dy = 0.
f ) 2yex/y dy + (y 2xx/y )dy = 0.
g) (xey/x y sin xy )dx + x sin xy dy = 0.
7. Hallar la solucin del problema de valor inicial dado
a) (x2 + y2 )dx = 2xydy; y(1) = 0.
b) (xey/x + y)dx = xdy; y(1) = 0.
c) (xy y2 )dx x2 dy = 0; y(1) = 1.
Solucin
1. a)
b)
c)
d)
2. a)
b)
c)
d)
Homognea de grado 3.
Homognea de grado 1/2.
Homognea de grado 2.
No homognea.
y = x(ln |x| + c).
x = (y x)[ln |y x| + c].
x = y(ln |y| + c).
y2= 2x2 (ln |x| + c).

2 1 51(2y/x)
y
y
e) ln1 x x 5 ln 5+1+(2y/x) = 2 ln |x| + c.

y 3+ 2x2 +3y2
f)
c = 0.
3
x
3.
g) 4xy + y2 + x2 = c.
x
h) 2(y+x)
ln |y + x| = 21 .

a) 21 ln(x2 + y2 ) + arctan y/x = c.
b) 23 ln |y| 13 ln | 3x + y| = c.
ln |xy|
c) ln |x+y|3
= c.
2

6.
7.
2.3.2
a) 3x2 y + y3= 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
a) y2 = x2 + x.
b) ln(|x|) + ey/x = 1.
c) x = ex/y1 .
43
y
= c.
x
Ecuaciones con Coeficientes Lineales

Hemos utilizado sustituciones de y que transforman la ecuacin diferencial para poder
resolverla, pero en algunos casos se necesita sustituir x y y en nuevas variables. Este es el
caso para las llamadas Ecuaciones Diferenciales con Coeficientes Lineales.
Consideremos la ecuacin
(a1 x + b1 y + c1 )dx + (a2 x + b2 y + c2 )dy = 0.
(2.19)
Para resolverla consideremos los siguientes casos:

dy
1. Si a1 b2 = a2 b1 , la ecuacin (2.19) se puede escribir de la forma dx
= g(ax + by + c)
que se resuelve por el mtodo Reduccin a Separacin de Variables.
2. Si c1 = c2 = 0, la ecuacin (2.19) se convierte en una ecuacin homognea.

3. Si a1 b2 6= a2 b1 la ecuacin se hace homognea si tomamos una translacin de ejes de
la forma x = u + h y y = v + k, donde h y k son constantes, y
dy dv
= .
dx du
Sustituyendo en la ecuacin (2.19), se tiene
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 )
=
.
du
a2 u + b2 v + (a2 h + b2 k + c2 )
(2.20)
Como a1 b2 6= a2 b1 , entonces se puede hallar h y k tales que

a1 h + b1 k + c1 = 0
a2 h + b2 k + c2 = 0
(2.21)
para h y k soluciones del sistema (2.21), se tiene que la ecuacin (2.20) es homognea.
44
Ejemplo 2.10. Hallar la solucin general de

(3x + 3y + 4)dx + (2x + 2y + 7)dy = 0.
Solucin. Como a1 b2 = 3 2 = 2 3 = a2 b1 . As,
(3x + 3y + 4)dx + (2x + 2y + 7)dy = 0
[3(x + y) + 4]dx + [2(x + y) + 7]dy = 0
dy
3(x + y) + 4
=
.
dx 2(x + y) 7
Sea u = x + y, por lo tanto, u0 = 1 + y0 y la ecuacin se convierte en
3u + 4
du
1 =
dx
2u 7
la cual es una ecuacin de variables separables
Z
dx =
2u + 7
du
3u
integrando:
x + c = 2u 9 ln |3 u|.
De esta forma, la solucin general de la ecuacin est dada por
x + c = 2(x + y) 9 ln |3 (x + y)|.

Ejemplo 2.11. Hallar una solucin de la ecuacin diferencial
(3x + y 1)dx + (x + y + 3)dy = 0.
(2.22)
Solucin. Como a1 b2 = 3 1 6= 1 1 = a2 b1 , usaremos la translacin x = u + h y y = v + k,

donde h y k satisfacen el sistema
3h + k 1 =0
h + k + 3 =0.
Al resolver este sistema en trminos de h y k tenemos que h = 1, k = 2, es decir x = u 1
y y = v 2. Como dy = dv y dx = du, sustituyendo en la ecuacin (2.22) se tiene
(3u + 3 + v 2 1)du + (u 1 + v 2 + 3)dv = 0
(3u + v)du + (u + v)dv = 0
dv 3u v
=
.
du
u+v

La ecuacin anterior es homognea, luego tomamos la sustitucin z = uv , donde

dz
z + u du
y al sustituir en la ltima ecuacin obtenemos
45
dv
du
dz 3u zu
=
du
u + zu
dz 3 z
z+u =
du 1 + z
dz 3 2z z2
u =
du
1+z
z+u
separando variables tenemos:

Z
du
z+1
dz =
2
3 2z z
u
Z
Z
z+1
du
dz =
2
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
ahora aplicando exponencial obtenemos:

z2 + 2z 3 = u2 ec1 = cu2 .
Al sustituir de nuevo z, u y v, vemos que
v 2
+2
v
3 = cu2
u
u
v + 2uv 3u2 = c
2
(y + 2)2 + 2(y + 2)(x + 1) 3(x + 1)2 = c.

La ltima ecuacin es la solucin implcita de (2.22).
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.
46
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.
Solucin
1.
2.
2.3.3
a)
b)
c)
d)
e)
a)
b)
c)
ln(4(y 1)2 + (x 2)2 ) 2 arctan

5
2 (x y) = c.
2y2
x2
= c.
ln(|5x + 10y 1|) +

x + 2y + ln(|x + y 2|) = c.
x + 2y + ln(|x + y 1|) = c.
x + 3y 3 ln(|x + 2y + 3|) = c.
x + 3y + 2 ln(|2 x y|) = 1.
(y + 7)2 (3x + y + 1) = 128.

ln((x 1)2 + (y + 3)2 ) + 2 arctan x1
y+3 = 2 ln(3).
Ecuacin de Bernoulli
dy
+ P(x)y = Q(x)yn , n R
(2.23)
dx
se denomina Ecuacin Diferencial de Bernoulli. Ntese en particular, que si n es 0 1, la
ecuacin es lineal o de variables separables y puede resolverse directamente. Si n no es 0
ni 1, la ecuacin es no lineal, en este caso, hacemos la sustitucin u = y1n , que reduce la
ecuacin de Bernoulli a una lineal. En efecto, derivando u con respecto a x
du
dy
= (1 n)yn ,
dx
dx
reemplazando en la ecuacin (2.23) tenemos
u = y1n ,
donde
dy
1 n du
=
y
dx 1 n dx
1 n du
y
+ P(x)y = Q(x)yn , n R
1 n dx
ahora multiplicamos la ecuacin (2.24) por (1 n)yn nos queda:
du
+ (1 n)P(x)y1n = (1 n)Q(x)
dx
sustituyendo u = y1n ,
du
+ (1 n)P(x)u = (1 n)Q(x)
dx
la cual es una ecuacin diferencial lineal.
(2.24)
47
Ejemplo 2.12. Hallar la solucin general de cada una de las siguientes ecuaciones diferenciales:
dy
y = ex y2 .
1. dx
1
dy
y = ex y 2 .
2. dx
Solucin.
1. En este caso n = 2, tomando la sustitucin obtenemos u = y1 , derivandy
dy
2 dy
2 du
do con respecto a x se tiene, du
dx = y dx y despejando dx tenemos dx = y dx
reemplazando en la ecuacin obtenemos:
y2
du
y = ex y2
dx
(2.25)
multiplicando la ecuacin (2.25) por (y2 ), nos da

du
+ y1 = ex
dx
(2.26)
en (2.26) sustituimos u = y1 :
du
+ u = ex ,
dx
(2.27)
como la ecuacin diferencial es lineal, el factor integrante para la ecuacin (2.27) es:
R
(x) = e
1dx
= ex
utilizamos la formula (2.14) nos resulta:

Z

x
x
x
u(x) =e
e (e )dx + c

2x
e
x
+c
u(x) =e
(2.28)
en (2.28) reemplazamos u = y1
ex + 2cex
2
2
y= x
.
e + 2cex
y1 =
dy
2. Para n = 12 , entonces u = y2 , por tanto dx
= 2u3 du
dx . Sustituyendo este diferencial
2
u = y en la ecuacin original tenemos
2u3
multiplicando por
1
2
3
u tenemos
du
u2 = ex u1 ,
dx
du 1
1
+ u = ex u2 ,
dx 2
2
48
una ecuacin tipo Bernoulli que requiere de una nueva sustitucin para solucionarla
como en el caso a), es decir, en este caso n = 2, tomando la sustitucin obtenemos
s = u1 , derivando con respecto a x se tiene,
ds
du
= u2
dx
dx
y despejando
du
dx
tenemos
ds
du
= u2
dx
dx
reemplazando en la ecuacin obtenemos:
u2
ds 1
1
+ u = ex u2 ,
dx 2
2
(2.29)
multiplicando la ecuacin (2.29) por (u2 ), nos da

1
ds 1 1
u = ex ,
dx 2
2
(2.30)
en (2.30) sustituimos s = u1 :
ds 1
1
s = ex ,
(2.31)
dx 2
2
como la ecuacin diferencial es lineal, el factor integrante para la ecuacin (2.31) es:
(x) = e
R 1
2 dx
Utilizamos la frmula (2.14) nos resulta:

Z

1
1
1 x
x
x
s(x) =e 2
e 2 e dx + c
2
1

1
s(x) =e 2 x e 2 x + c
= e 2 x .
(2.32)
en (2.32) reemplazamos s = u1
1
u1 = ex + ce 2 x
1
u=
1 .
ex + ce 2 x
Ahora sustituimos u = y 2 es decir:

1
y 2 =
1
1
ex + ce 2 x

2
1
y = ex + ce 2 x .
49
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.
Solucin
1. a) y ln(|x|) + y + cxy = 1.
b) y12 = ce2x + x + 21 .
2.
2.3.4
c)
d)
e)
a)
b)
c(1x)
3
y3/2 = 4(1+x+x
2 ) + 1+x+x2 .
1 2y(1 + ln(|x|)) = cxy.
x2
y = (x1)c+1
.
y = 1.
sec(x) = y2 (tan(x) + 1).
Ecuacin de Ricatti
Se denomina Ecuacin Diferencial de Ricatti a una ecuacin diferencial de primer orden de
la forma
dy
= P(x) + Q(x)y + R(x)y2 .
dx
(2.33)
Tal ecuacin se puede resolver con dos sustituciones consecutivas siempre que conozcamos
una solucin particular y1 de dicha ecuacin. Esto es
dy1
= P(x) + Q(x)y1 + R(x)y21 .
dx
(2.34)
Para resolverla tomamos la sustitucin y = y1 + u, la cual reduce la ecuacin a una de tipo

Bernoulli.
Primero derivamos con respecto a x, es decir:
y = y1 + u,
dy dy1 du
=
+ .
dx
dx dx
Reemplazando en la ecuacin (2.33) tenemos

dy1 du
+
= P(x) + Q(x)(y1 + u) + R(x)(y1 + u)2 .
dx dx
(2.35)
50
Desarrollando el cuadrado y agrupando (2.35) nos queda:

dy1 du
+
= P(x) + Q(x)y1 + Q(x)u + R(x)y21 + R(x)2y1 u + R(x)u2
dx dx
sustituyendo (2.34) y simplificando tenemos:

du
= Q(x) + R(x)2y1 u + R(x)u2
dx

du
Q(x) + R(x)2y1 u = R(x)u2 .
dx
La cual es una ecuacin de Bernoulli con n = 2.
Ejemplo 2.13. Hallar una familia mono paramtrica de soluciones de la ecuacin diferencial
dy
4 1
= 2 y + y2
dx
x
x
en donde y1 =
2
x
es una solucin conocida de la ecuacin.
Solucin. Tomando la sustitucin y = y1 + u = 2x + u, derivando con respecto a x se tiene,

dy
2
du
dx = x2 + dx y reemplazando en la ecuacin obtenemos:

2
4 1 2
2
2 du
= 2
+u +
+u .
2+
x
dx
x
x x
x
Desarrollando el cuadrado y simplificando tenemos:
2 du
4
2 u 4 4u
= 2 2 + 2 + + u2
+
2
x
dx
x
x
x x
x
du 3u
=
+ u2
dx x
du
3
u = u2 .
dx
x
(2.36)
Como la ecuacin (2.36) es de Bernoulli con n = 2, tomando la sustitucin obtenemos

ds
du
s = u1 , derivando con respecto a x se tiene, dx
= u2 du
dx y despejando dx tenemos
du
ds
= u2 .
dx
dx
Reemplazando en la ecuacin (2.36) obtenemos:

3
2 ds
u
u = u2 ,
dx
x
(2.37)
51
multiplicando la ecuacin (2.37) por (u2 ), nos da

ds
3 1
+
u = 1.
dx
x
En (2.38) sustituimos s = u1 :

ds
3
+
s = 1.
dx
x
(2.38)
(2.39)
Como la ecuacin diferencial es lineal, el factor integrante para la ecuacin (2.39) es:
R 3
x dx
(x) = e
= e3 ln |x| = x3 .
Utilizando la frmula (2.14) nos resulta:

Z

3
3
s =x
x (1)dx + c
4

x
3
s =x
+c
4
4
x + c
s=
4x3
(2.40)
en (2.40) reemplazamos s = u1
x4 + c
4x3
4x3
u= 4
.
x + c
u1 =
(2.41)
Sustituyendo u = y 2x en (2.41) y despejando y, tenemos la solucin:

4x3
2
y=u=
+ .
4
x + c x
Ejercicios
1. 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.
52
Solucin
1.
a) u0 +
b)
c)
d)
2.3.5
2
x

+ 2x u = 1x ;
2 x2
y = x2 + 2xx2 e .
u0 2u tan(x) = sin(x);
u0 3x u = 1;
u0 1x u = x12 ;
y=
e +c
3 cos2 (x)
1
.
+
cos(x)
ccos3 (x)
1
x
y = + cx32x .
y = x + cx2x
2 1 .
Ecuacin de Lagrange
y = g(y0 )x + f (y0 ),
se denomina ecuacin de Lagrange. Para resolverla tomamos la sustitucin y0 = p para
obtener la solucin general de la ecuacin en forma paramtrica mediante la cual se reduce
a:
y = g(p)x + f (p).
(2.42)
Derivando en (2.42) con respecto a x, obtenemos:

y0 = g0 (p)p0 x + g(p) + f 0 (p)p0 .
(2.43)
Reemplazando y0 = p en la ecuacin (2.43) tenemos

p = g0 (p)p0 x + g(p) + f 0 (p)p0
p g(p) = g0 (p)p0 x + f 0 (p)p0
p g(p) = p0 (g0 (p)x + f 0 (p))
dp 0
(g (p)x + f 0 (p)),
p g(p) =
dx
despejando obtenemos:
dp
p g(p)
= 0
.
dx g (p)x + f 0 (p)
Como y0 = p y
dp
dx
son funciones inversas, intercambiando estas en (2.44) nos da
dx g0 (p)x + f 0 (p)
=
dp
p g(p)
0
dx
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)
(2.44)
53
donde (p) es el factor integrante, y p es el parmetro. As la expresin para y es

y = g(p)x(p, c) + f (p).
(2.45)
Las ecuaciones paramtricas de la solucin general de la ecuacin de Lagrange (2.45) son

x(p) = h(p, c) y y = g(p)x(p, c) + f (p) donde c es la constante de la familia de curvas.
Ejemplo 2.14. Hallar las ecuaciones paramtricas de la solucin de la ecuacin diferencial
y = x(y0 + 3) 2(y0 )2 .
Solucin. Tomando la sustitucin y0 = p y derivando con respecto a x se tiene
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
intercambiando y despejando se tiene
dx
x 4
= + p
dp
3 3
dx x 4
+ = p.
dp 3 3
(2.46)
Como la ecuacin (2.46) es lineal, obtenemos P(p) = 13 , Q(p) = 34 p, con factor integrante
R 1
3dp
= e 3 p . La solucin de esta ecuacin viene dada por la frmula (2.14), es
(p) = e
decir;
13 p
x(p) = e
Z
1
3p

4
pd p + c
3
1
x(p) = 4(p 3) + ce 3 p .
Las ecuaciones paramtricas de la solucin de la ecuacin diferencial son

1
x(p) = 4(p 3) + ce 3 p

13 p
y = 4(p 3) + ce
(p + 3) 2p2 .
54
Ejercicios
1. Hallar las ecuaciones paramtricas de la soluciones de las ecuaciones diferenciales
siguientes
a) y = x(y0 )2 + ln(y0 ).
0 2
0 3
b) y = x(y
)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 .
Solucin
a) x(p) =
1.
1p ln(p)+c
;
(p1)2
p3 + 23 p2 +c
;
(p1)2
y=
1p ln(p)+c
(p1)2
p3 + 3 p2 +c
p2 + ln(p).
2
y = p2 (p1)
+ p3 .
2
c) x(p) = (p+2 p+2 ln( p1)+c)(p1/2 +1);

1) + c)(p1/2 + 1) + 21 p2 23 p3/2 .
b) x(p) =
d)
2.3.6
y=
p(p+2 p+2 ln( p

2
1 2 cos (p)+4 sin(p)4c+1
+
y = (cos(p)+ p) 4
sin(p)+1
2
sin(p)4c+1
x(p) = 14 2 cos (p)+4
;
sin(p)+1
1
1
2 sin(p) + 2 p.
Reduccin a Separacin de Variables

Estudiaremos ahora una clase ecuaciones que se pueden resolver como ecuaciones de variables separables haciendo una sustitucin conveniente. Consideremos la ecuacin diferencial
dy
= g(ax + by + c).
dx
Aplicando la sustitucin u = ax + by + c sta se reduce a una ecuacin con variables
separables.
Ejemplo 2.15. Resolver
dy
= sin(x + y).
dx
Solucin. Si hacemos u = x + y, entonces
se transforma en
du
dx
dy
= 1 + dx
, por lo que la ecuacin diferencial
du
du
1 = sin(u), o sea
= 1 + sin(u)
dx
dx
separamos variables, e integrando obtenemos:
Z
1
du = dx
1 + sin(u)
tan(u) sec(u) =x + c
55
ahora como u = x + y sustituimos en la ltima ecuacin y nos queda como resultado:

tan(x + y) sec(x + y) = x + c.

Ejercicios
1. Hallar las soluciones de las ecuaciones diferenciales siguientes
a) y0 = (x y)2 .
b) y0 = sin(x y + 1).
dy
= (x + y 1)2 .
c) dx
dy
d) dx
= tan2(x + y).
e) y0 = 2 + y 2x + 3.
dy
= 1xy
f ) dx
x+y .
g)
dy
dx
= 1 + eyx+5 .
Solucin
2x ce2x x1
1. a) y = cxe ce
.
2x 1

b) y = x 1 + 2 arctan cx2
cx .
c) y = x 1 tan(c x).
tan(x+y)
d) 1+tan
2 (x+y) + 2y = x + c.
e) x 2 y 2x + 3 c = 0.
f ) |x y| = 2x + c.
g) y = 5 ln(c + ex ).
3 Aplicaciones EDO de Primer Orden
El Propsito de este captulo es demostrar la utilidad de las ecuaciones diferenciales

ordinarias de primer orden. Para ello mostraremos algunas situaciones tpicas que se modelan
mediante este tipo de ecuaciones.
3.1
Crecimiento y Decrecimiento
Ahora nos concentraremos en aplicaciones modeladas por EDO de primer orden llamadas
ecuaciones autnomas; esto es, ecuaciones de la forma
dy
= f (y).
dt
La ley de Malthus viene descrita por una ecuacin de este tipo. Esta ley se expresa en
smbolos matemticos como
dx
= kx,
dt
x(t0 ) = x0
(3.1)
en donde k es una constante. Esta aparece en muchas aplicaciones de distintas reas del
conocimiento que involucran crecimiento o decrecimiento de una cantidad modelada en
el tiempo. Por ejemplo, en biologa a menudo se observa que la rapidez con que ciertas
bacterias se multiplican es proporcional al nmero de bacterias presentes en cierto instante.
Para cortos intervalos de tiempo, la magnitud de una poblacin de animales pequeos,
como roedores, puede predecirse con bastante exactitud mediante la solucin de (3.1). En
fsica, un problema de valor inicial como (3.1) proporciona un modelo para aproximar la
cantidad que queda de una sustancia radioactiva que se desintegra, esta ecuacin en forma
simplificada sirve como modelo para determinar la temperatura de un cuerpo que se enfra.
Aplicaciones EDO de Primer Orden
58
En qumica, el mismo modelo diferencial sirve para determinar la cantidad de sustancia que
queda durante una reaccin.
La constante de proporcionalidad k puede ser negativa o positiva y se puede determinar si
se establece otro valor de la variable x en un tiempo t1 > t0 .
Ejemplo 3.1. Se sabe que la poblacin de cierta comunidad aumenta en un instante 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 fenmeno es el siguiente:
dP
= kP,
dt
(3.2)
sujeta a P(0) = P0 , esto produce el modelo matemtico llamado ley exponencial o de

Malthus para el crecimiento de poblaciones. Se pretende hallar el valor de t para el cual
P(t) = 3P0 .
La ecuacin diferencial en (3.2) es de variables separables, por tanto:
dP
= kdt
P
integrando tenemos, ln |P| = kt + c de donde despejando P obtenemos P(t) = Cekt .
Reemplazando la condicin inicial P(0) = P0 se tiene C = P0 por tanto:
P(t) = P0 ekt .
(3.3)
Ahora como P(5) = 2P0 remplazando en (3.3) se tiene:

2P0 =P0 e5k
2 =e5k
ln(2) =k5
ln(2)
k=
5
P(t) =P0 e
ln(2)
5 t
Para determinar el valor de t para el cual la poblacin se ha triplicado, despejamos t de

3P0 = P0 e
ln(2)
5 t
3.1 Crecimiento y Decrecimiento
59
luego cancelando y tomando la funcin logaritmo en ambos miembros de la ecuacin

tenemos:
ln(2)
t
5
ln(3)
t =5
7,9.
ln(2)
ln(3) =
As, se demora en triplicarse t 7,9 aos.
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

si el inters es compuesto es semestralmente,

P =A 1 +
2 100
4

si el inters es compuesto es trimestralmente,

P =A 1 +
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/r nr
r m in
P = lm A 1 +
= lm A 1 +
= Aenr .
m
m
m
m
Finalmente, reemplazando n por t obtenemos que
P = Aert
es el capital al final de un tiempo t si la cantidad A es compuesta instantneamente o
continuamente a un inters del r por ciento anual. La ecuacin anterior es la solucin de la
ecuacin diferencial
dP
= rP.
dt
60
Ejemplo 3.2. Cuanto tiempo toma duplicar $50 si el inters es compuesto continuamente
al 6 % anual?.
Solucin. En este caso r = 6/100 = 0, 06. La condicin inicial es P(0) = 50 y queremos
hallar t tal que P(t ) = 100. Tenemos
dP
= 0, 06P
dt
que es una ecuacin de variable separable (tambin lineal) y as
Z
dP
=
P
0, 06dt,
ln(P) = 0, 06t + k
por lo tanto P(t) = Ke0,06t , usando la condicin inicial obtenemos 50 = P(0) = Ke0,060 =
K. Sustituimos esto en la solucin de la ecuacin para t ; esto es: 100 = P(t ) = 50e0,06t .
Finalmente, ln(2) = 0, 06t lo cual implica t 11, 55 aos.

Ejercicios
1. 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.
2. 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 crecimeinto de poblaciones y estime la poblacin de
lagartos en dichos terrenos en el ao 2010.
3. 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?.
4. A qu tasa de inters se duplicarn $300 en 8 aos?.
5. Cual ser el rdito de $1.000 al 4,5 % de inters luego de 20 aos?.
6. 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
retira continuamente desde la fecha del deposito una cantidad de 3.600/365 al da?.
Solucin
1. Aproximadamente 27.416 peces.
2. 4.828,7.
3. 9.885,6.
4. Aproximadamente 8,7 %.
5. 568,31.
6. A =45.513..

3.1.1
61
Modelo Logstico de Crecimiento Poblacional

El modelo anterior (ley de Malthus) debido a su simplicidad, es solo aplicable a pequeas poblaciones. En aquellas que son grandes, se tendrn algunas restricciones: recursos,
espacio, alimento, etc., por lo tanto esto debera ser modelado. El modelo logstico tiene
en cuenta esta situacin: Supongamos que se ha estimado que cierto espacio puede tolerar
mximo una poblacin de tamao N 1 (este se puede encontrar por experimentacin, simulacin, etc.) y que si la poblacin es pequea (con respecto a N), la tasa de crecimiento de
sta crecer proporcionalmente a su tamao. Si la poblacin es demasiado grande (respecto
a N), la tasa de crecimiento es negativa. El modelo viene dado en la siguiente forma:

P
dP
= kP 1
.
dt
N
Ntese que el factor 1 NP es positivo cuando P < N y negativo cuando P > N. La ecuacin
dada es de variables separables. Veamos cmo se resuelve:

dP
P
= kP 1
dt
N
dP
= kdt
P 1 NP
Z
Z
dP
= kdt = kt +C1 .
P 1 NP
La integral del lado izquierdo se hace por fracciones parciales

Z
Z
dP
1
1/N
=
+
dP
P 1 P/N
P 1 NP

P

P

= ln |P| ln 1 = ln

P
1 N
N

NP

= ln
N P
de donde
o equivalentemente,

NP
= kt +C1
ln
N P
NP
= ekt+C1 = ekt eC1 = C2 ekt .
N P
Debemos considerar dos situaciones:
1 Usualmente
recibe el nombre de soporte.
62
(1.) Si N P > 0, entonces
NP
= C2 ekt
N P
NP = NC2 ekt PC2 ekt

P N +C2 ekt = NC2 ekt
P=
NC2 ekt
=
N +C2 ekt
P=
N
NC3 ekt + 1
N kt
C2 e
+1
donde C3 = C12 .
(2.) Si N P < 0, entonces
NP
= C2 ekt
N P
NP = NC2 ekt + PC2 ekt

P N C2 ekt = NC2 ekt
P=
N
NC2 ekt
= N kt
kt
N C2 e
1
C e
2
N
P=
1 NC3 ekt
donde, otra vez, hemos hecho C3 = C12 .

Si la poblacin inicial es P(0) = P0 < N entonces la poblacin crecer (pues
suponiendo k > 0) hasta el soporte:
lm P = lm
t 1 NC3 ekt
dP
dt
es positivo
= N.
En cualquiera de las dos situaciones la poblacin tiende a estar en equilibrio. Si P(0) =

P0 = 0, la poblacin se mantendr en cero, naturalmente. Finalmente, si P(0) = P0 =
N entonces la poblacin se mantendr constante en P(t) = N. En resumidas cuentas la
poblacin viene dada mediante la frmula
, si 0 < P0 < N
NC3 ekt +1
N
, si P0 > N
P(t) = 1NC3 ekt
0,
si P0 = 0
N,
si P0 = N.
En la Figura 3.1 se presentan dos situaciones hipotticas en que el soporte es 3 (asntota
horizontal). La primera, la poblacin inicial es mayor que el soporte y a lo largo del tiempo
63
decrecer. Y en la segunda ecuacin, la poblacin inicial es menor que el soporte (pero

positiva) y la poblacin crecer hasta el soporte.
Figura 3.1: Asntota horizontal modelo logstico.
Ejemplo 3.3. En cierta ciudad hay 4.800 personas susceptibles de contraer una cierta enfermedad 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 enfermedad).
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
NC3 ekt + 1
4.800
=
.
4.800C3 ekt + 1
P(t) =
64
La condicin inicial es P(0) = 300, as
4.800
= 300
4.800C3 + 1
4.800 =1.440.000C3 + 300.
P(0) =
De donde C3 = 0,003125 y actualizando la ecuacin tenemos

P(t) =
4.800
.
15ekt + 1
Se sabe tambin que P(10) =1.200, lo que permite encontrar el parmetro k.

P(10) =
4
4.800
= 1.200
15ek10 + 1
15ek10 + 1
=1
1
5
1
1
k = ln
0,16094.
10
5
ek10 =
Por tanto la expresin para la poblacin de enfermos est dada por

P(t) =
4.800
15e0,16094t
+1
1. Para calcular el nmero aproximado de enfermos al cabo de 20 das, hacemos

P(20) =
4.800
15e0,16094(20) + 1
3.000.
2. Para encontrar el tiempo (en das) en que la mitad de la poblacin se ha infectado

resolvemos
P(t) =2.400
4.800
4. 800
=
= 2. 400
2
15e0,16094t + 1
2
=1
0,16094t
15e
+1
15e0,16094t + 1 =2
1
e0,16094t =
15

1
t = 0,16094 ln
18,8 das.
15
65
Ejercicios
1. Las observaciones sobre el crecimiento de tumores en animales indican que el tamao
y(t) del tumor al tiempo t puede describirse mediante la ecuacin diferencial
y
dy
= ky ln
,
dt
a
donde k y a son constantes positivas. Esta ecuacin diferencial en ocasiones se conoce
como ley de crecimiento de Gompertz. Resuelva esta ecuacin diferencial.
2. Un modelo de poblaciones utilizado en las predicciones actuariales se basa en la
ecuacin de Gompertz
dP
= P(a b ln(P)),
dt
donde a y b son constantes.
a) Halle P(t) en la ecuacin de Gompertz.
b) Si P(0) = P0 > 0, d una frmula para P(t) en trminos de a, b, P0 y t.
c) Describa el comportamiento de P(t) cuando t +. (Sugerencia: Considere
los casos para b > 0 y b < 0).
3. 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?.
Solucin
kt
1. y(t) = aee c .
2.
a) p(t) = e
b) p(t) = e
3. 468,4.
3.1.2
aebtbc
b
aaebt
b
0
p ebt .
Determinacin de Edades por el Mtodo del Carbono 14

La teora de la datacin por carbono se basa en que el isotopo Carbono 14 se produce en
la atmsfera por la accin de la radiacin csmica sobre el nitrgeno. El cociente de la
cantidad de 14C y la cantidad de carbono ordinario presentes en la atmsfera es constante
y, en consecuencia, la proporcin de isotopo presente en todos los organismos vivos es la
misma que en la atmsfera. Cuando un organismo muere, la absorcin de 14C cesa. As,
comparando la proporcin de 14C que hay en un fsil con la proporcin constante encontrada
en la atmsfera, es posible obtener una estimacin razonable de su edad. El mtodo usa
adems el hecho que la vida media del 14C radioactivo es de aproximadamente 5600 aos.
Ejemplo 3.4. El isotopo radioactivo de plomo, 209 Pb, se desintegra, en un instante cualquiera, con una rapidez proporcional a la cantidad presente en dicho instante, y tiene una
vida media de 3,3 horas. Si inicialmente hay 1 gramo de plomo, cunto tiempo transcurrir
para que se desintegre el 90 % de dicho elemento?.
66
Solucin. Sea A(t) la cantidad de plomo que queda en un instante cualquiera. La ecuacin
diferencial con condicin inicial es:
dA
= kA,
dt
A(0) = 1 gramo.
Como en el Ejemplo 3.1, se tiene que A(t) = Cekt donde notamos por la condicin inicial que
C = 1 y as A(t) = ekt . Como el isotopo 209 Pb tiene una vida media de 3,3 horas, tenemos
que A(3,3h) = 12 gramos, entonces para evaluar k se debe resolver 12 = e3,3k . Despejando k
resulta

ln 1/2
k=
0,210045.
3,3
De esta forma, A(t) = e0,210045t .
Como se desintegra el 90 % de la cantidad inicial, entonces queda 10 % de la sustancia. Para

ln(0,1)
hallar t se debe resolver 0,1 1gramos = e0,210045t . Despejando t resulta t = 0,210045
10,96 11 horas. Es decir, para que se desintegre pasan 11 horas.

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.
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)
= 0,00012378.
k =
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.
3.2 Mezclas
67
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), 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 %?.
5. Muchos creen que el sudario de Turn, que muestra un negativo de la imagen de
un cuerpo de un hombre crucificado, es la mortaja de Jess de Nazaret. En 1988,
el Vaticano otorg autorizacin para que se fechara el carbono del manto. Tres
laboratorios cientficos independientes, que analizaron la tela, llegaron a la conclusin
que tiene unos 600 aos, edad que coincide con su aparicin histrica. Con esta edad
determine qu porcentaje de la cantidad original de 14C queda en la tela en 1988.
Solucin
1. a) M(t) = 150e0,18932t .
b) 1,17 aos.
2. 1.512,9 aos.
3. a) 3,12 %
b) 31 aos 219 das.
4. 14,3 horas.
5. Queda el 92,8 % de 14C.
3.2
Mezclas
Algunos problemas que involucran mezclas de fluidos dan lugar a plantear una ecuacin
diferencial lineal de primer orden que modele la cantidad de mezcla o la concentracin de
una sustancia en un fluido.
68
En los problemas se requiere determinar la cantidad A(t) de una sustancia que hay en un
tanque en cada instante t. La razn de cambio de la cantidad de la sustancia presente en el
tanque es igual a la velocidad de entrada menos la velocidad de salida
dA
= (Rapidez de entrada) (Rapidez de salida)
dt
dV
= R1 R2
dt
(3.4)
donde
Rapidez de entrada = (flujo de entrada) (concentracin).
Suponiendo que el flujo de entrada es igual al de salida, la concentracin se halla dividiendo
la cantidad A(t) por el volumen de la mezcla que hay en el instante t. As
Rapidez de salida = (flujo de salida) (concentracin).
Flujo de entrada
A(t) = Cantidad en el tiempo t
dA
dt
=(Rapidez entrada)x(Rapidez salida)
Flujo de salida
Ejemplo 3.6. La corriente sangunea lleva un medicamento hacia el interior de un rgano a

razn de 3 cm3 /seg y sale de l a la misma velocidad. Se sabe que el volumen del rgano
es de 125 cm3 y la concentracin del medicamento que entra es de 0, 2 g/cm3 . Cul es
la concentracin del medicamento en el rgano en el instante t si inicialmente no haba
vestigios alguno del medicamento?. Cundo la concentracin ser de 0, 1 g/cm3 ?.
Solucin. Primero debemos determinar la rapidez con que un medicamento entra al interior
de un rgano, sabemos que el medicamento entra a una razn constante de 3 cm3 /seg. La
concentracin del medicamento que entra es de 0, 2 g/cm3 , concluimos as que la rapidez
de entrada en el rgano es:
R1 = 3
cm3
g
g
0, 2 3 = 0, 6
.
seg
cm
seg
3.2 Mezclas
69
Ahora debemos determinar la rapidez de salida. Como el flujo de entrada es igual al de

salida, la concentracin es A(t)/(el volumen del rgano). Por lo tanto la rapidez con que el
medicamento sale del rgano es:
R2 =
A(t) g
cm3 3A(t) g
3
=
.
125 cm3
seg
125 seg
En un principio no hay medicamento en el rgano, de modo que A(0) = 0 g, al sustituir R1

y R2 en (3.4) tenemos:
dA
3A(t)
= 0, 6
dt
125
A(0) = 0
dA 3A
+
= 0, 6.
dt 125
(3.5)
3
La ecuacin (3.5) es lineal, donde P(t) = 125
y Q(t) = 0,6, entonces podemos hallar el
factor integrante sustituyendo en la formula (2.13)
R 3
125 dt
(t) = e
= e 125 t
y reemplazando en la frmula (2.14) tenemos

Z

3
3
125
t
t
A(t) = e
0,6e 125 dt + c
h
i
3
3
125
t
t
125
=e
25e + c
3
= 25 + ce 125 t .
Al usar la condicin inicial A(0) = 0 para hallar la constante, obtenemos

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) =
.
5
(3.6)
70
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
En t0 29 seg la concentracin ser de 0, 1 g/cm3 .
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 %, determine 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
dt
(3.7)
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
20
tanque tiene cido ntrico al 20 %, la concentracin es 100
, concluimos que la rapidez de
entrada en el tanque es:
R1 = 6
6 litros
litros 20
=
.
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
=
,
dt
5 100 t
dV
4V
6
+
= .
dt
100 t 5
V (0) = 1litro
(3.8)
3.2 Mezclas
71
4v
La ecuacin (3.8) es lineal donde P(t) = 100t
y Q(t) = 65 , entonces podemos hallar el
factor integrante sustituyendo en la frmula (2.13),
R
(t) = e
4
100t dt
= e4 ln(100t) = (100 t)4
y reemplazando en la frmula (2.14) tenemos

"Z
#
1
6
(100 t)4 dt + c
V (t) =
(100 t)4
5

1
2(100 t)3
=
+c
(100 t)4
5
= 0,4(100 t) + c(100 t)4 .
Al usar la condicin inicial V (0) = 1 para hallar la constante, obtenemos

V (0) = 0,4(100 0) + c(100 0)4 = 1
40 + 105 c = 1
c = 3,9 104 .
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.
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 extremadamente 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. Calcule la
cantidad de libras de sal que hay en el tanque a los 30 minutos.
72
4. Un depsito contiene 50 litros de una solucin compuesta por 90 % de agua y 10 %

de alcohol. Se vierte 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 ?.
3.3
Problemas de Temperatura
La ley del enfriamiento de Newton dice que en un cuerpo que se enfra, la rapidez con que
la temperatura T (t) cambia es proporcional a la diferencia entre la temperatura del cuerpo y
la temperatura constante Tm del medio que lo rodea. Esto es,
dT
= k(T Tm )
dt
(3.9)
en donde k es una constante de proporcionalidad.

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)
dt
T (0) = 700 F

1
T
= 500 F.
2
(3.10)
Se debe resolver la ecuacin (3.10) con condiciones iniciales, la cual es lineal y tambien de
3.3 Problemas de Temperatura
73
variables separables. Empleando la tcnica de variables separables obtenemos:

dT
=kdt
T 10
Z
Z
dT
= kdt
T 10
ln | T 10 |=kt + c1
T 10 =c2 ekt .
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
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.
Ahora queremos formular un mdelo matemtico que describa el perfil de temperatura

dentro de un edificio durante 24 horas como funcin de la temperatura exterior, el calor
generado dentro del edificio y el calefactor o aire acondicionado.
Sea T (t) la temperatura dentro del edificio en el instante t y veamos al edificio como un
nico comportamiento. La razn de cambio en la temperatura queda determinada por todos
los factores que generan o disipan calor.
Factores que afectan la temperatura dentro del edificio:
Calor generado por las personas, las luces y las mquinas dentro del edificio. Esto
incrementa la temperatura y se denota por H(t). Siempre es positiva.
74
Calentamiento o enfriamiento proporcionado por calefaccin o aire acondicionado.

Esta razn de incremento (o decremento) de temperatura se denota por U(t). Que es
positiva para calefaccin y negativa para aire acondicionado. Se tiene que H(t) y U(t)
estn en trminos de temperatura por unidad de tiempo.
Efecto de la temperatura exterior M(t) sobre la temperatura dentro del edificio.
El tercer factor lo podemos modelar mediante la ley de enfriamiento de Newton. Esto es: la
razn de cambio de la temperatura T (t) es proporcional a la diferencia entre la temperatura
exterior M(t) y la temperatura interior T (t). Es decir:
dT
= K[M(t) T (t)]
dt
la constante depende de las propiedades fsicas del edificio, como cantidad de puertas,
ventanas y el tipo de aislamiento. Por lo tanto,
dT
= K[M(t) T (t)] + H(t) +U(t)
dt
(3.11)
donde
dT
+ KT (t) = KM(t) + H(t) +U(t)
dt
es una ecuacin diferencial lineal con
p(t) = K
Q(t) = KM(t) + H(t) +U(t).
El factor integrante es
(t) = e
Kdt
= eKt .
Utilizando la frmula (2.14) nos resulta

Z

Kt
Kt
T (t) = e
e [KM(t) + H(t) +U(t)]dt + c
y
1
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?.
Solucin. Se tiene M(t) = 350 C, T (0) = 240 C, H(t) = 0, U(t) = 0,
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
1
K
= 4 por lo tanto

obtenemos:
75

1
e 35dt + c
T (t) = e
4
h 1
i
1
T (t) = e 4 t e 4 t 35 + c
14 t
Z
1
4t
T (t) = 35 + ce 4 t .
Como T (0) = 240 C, tenemos:

1
24 = 35 + ce 4 t
c = 11.
1
Por lo tanto, T (t) = 35 11e 4 t de donde reemplazando podemos hallar T (2) = 35

1
11e 4 2 = 28,320 C.
La temperatura a las 2 : 00 p.m. es de 28,320 C. La temperatura a las 6 : 00 p.m. es de
32,540 C. Para hallar en qu momento t0 la temperatura llega a 270 C reemplazamos en
1
1
T (t0 ) = 35 11e 4 t0 ; es decir, 27 = 35 11e 4 t0 , despejando t0 obtenemos t0 = 1 : 16
p.m..

Ejemplo 3.10. Se va a construir un almacn sin calefaccin ni aire acondicionado. Segn la
cantidad de aislamiento, la constante de tiempo para este edificio puede variar de 1 a 5 horas.
Para ilustrar el efecto del aislamiento sobre la temperatura dentro del almacn, suponga que
la temperatura exterior vara como una onda sinusoidal, con un mnimo de 160 C a las 2 : 00
a.m. y un mximo de 320 C a las 2 : 00 p.m.. Suponiendo que el trmino exponencial (que
implica la temperatura inicial T0 ) se ha extinguido, cul es la temperatura mnima dentro
del edificio, si la constante de tiempo es 1 hora?. Y si la constante de tiempo es 5 horas?.
Cul es la mxima temperatura dentro del edificio si la constante de tiempo es 1 hora?. Y
si es 5 horas?.
Solucin. Como M(t) vara como una onda sinusoidal en un periodo de 24 horas, con un
mnimo 160 C a las 2 : 00 a.m., y un mximo de 320 C a las 2 : 00 p.m., t = 12. Es decir
M(t) = M0 B cos(wt)
donde B es una constante positiva, M0 es la temperatura exterior promedio,
32 + 16
2
radianes
= 240 C y w =
=
.
2
24 12 hora
Adems U(t) = 0 porque no hay calefaccin ni aire acondicionado. H(t) = H0 = 0 no hay
razn de calentamiento dentro del edificio. De la ecuacin (3.11) tenemos

dT
= K M0 B cos(wt) T (t) + H0
dt

dT
+ KT = K M0 B cos(wt) + H0
(3.12)
dt
M(t) =
76
hacemos B0 = M0 + HK0 , reemplazando en (3.12) tenemos que:

dT
+ KT = K [B0 B cos(wt)]
dt
donde KB0 representa el valor promedio diario de KB0 KB cos(wt), es decir:
1
24
KB0 =
Z 24
0
(KB0 KB cos(wt)) dt
la solucin de la ecuacin es:

Z
Kt
Kt
T (t) = e
e [KB0 KB cos(wt)] dt + c

B cos(wt) + Kw sin(wt)
T (t) = B0
+ ceKt

w 2
1+ t
para t = 0 y T0 = 160 C reemplazando en la ltima ecuacin se tiene:

16 = B0
B
1+
por lo tanto
c = 16 B0 +
es decir

w 2
t
B
1+
+c

w 2
t
!

B cos(wt) + Kw sin(wt)
B
Kt
T (t) = B0
+ 16 B0 +

e
w 2
w 2
1+ t
1+ t
como H0 = 0, B0= M0 es decir, 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 cos(wt) + wk sin(wt)
T (t) = 24
2
1 + wt
para hallar el valor de B tomamos
M(t) = M0 B cos(wt)

M(t) = 24 B cos
t
12
77
para t = 0 y M(t) = 16 se tiene B = 8. Por lo tanto:

8 cos(wt) + Kw sin(wt)
.
T (t) = 24
2
1 + wt
Para hallar la temperatura mnima para t = 0, K1 = 1, se tiene K = 1 y tomamos w =

es aproximadamente 14 , reemplazando en la ltima ecuacin obtenemos:
8
T (0) = 24
para
1
K

5 2
4
= 20,80 C.
1+
8
1+
Para hallar la temperatura mxima t = 12,

T (12) = 24 +
1
K
= 16,40 C
que
= 5, se tiene K = 15 . La temperatura mnima es:

T (0) = 24
para

1 2
4
2
24
= 5, se tiene K =
1
5
1
K
= 1, se tiene K = 1
8
1+

1 2
4
= 31,50 C

5 2
4
= 27,10 C.
la temperatura mxima es:

8
T (12) = 24 +
1+
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?.
Solucin. K11 es la constante de tiempo con calefaccin y aire acondicionado donde K1 = K +
Ku. En la camioneta se instala un termostato que se utiliza para comparar la temperatura real
dentro de la camioneta con una temperatura deseada TD . Si la temperatura real es menor que
la temperatura deseada, el calefactor comienza a funcionar y en caso contrario se desconecta.
Si la temperatura real es mayor que la temperatura deseada, el aire acondicionado comienza
a enfriar y en caso contrario se desconecta. La cantidad de calentamiento es proporcional a
la diferencia de temperatura; es decir;
U(t) = Ku [TD T (t)]
78
donde Ku es la constante de proporcionalidad (positiva). As se tiene M(t) = 350 C, T (t) =

550 C, H(t) = 0, TD = 160 C, K1 = 2. Entonces K = 12 , K11 = 31 , luego K1 = 3. Como U(t) =
Ku [TD T (t)] de la ecuacin (3.11) se tiene
dT
= K [M(t) T (t)] + Ku [TD T (t)]
dt
dT
+ T [K + Ku ] = KM(t) + Ku TD
dt
vemos que para este ejemplo, la cantidad p(t) = K es igual a K + Ku , por lo tanto
(3.13)
K1 = K + Ku
1
3 = + Ku
2
5
Ku =
2
reemplazando en (3.13) tenemos:
dT
1
5
+ 3T = 35 + 16
dt
2
2
115
dT
+ 3T =
dt
2
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
+ ce3t .
T=
6
Como T (0) = 55, tenemos
115
+c
6
215
c=
6
115 215 3t
T (t) =
+
e ,
6
6
para hallar en qu momento t0 la temperatura llegara a 270 C reemplazamos en la ltima
ecuacin y despejamos t0 .
55 =
115 215 3t0

+
e
6
6
t0 = 0,5068 horas
t0 = 30,4 minutos.
27 =
79
Ejercicios
1. Un termmetro se lleva del interior de una habitacin al exterior, donde la temperatura
del aire es 60 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 determine 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.
(a.) Suponga que se aplica la ley de Newton sobre el enfriamiento. Escriba un
problema de valor inicial que modele la temperatura del t caliente.
(b.) Cunto tiempo le toma al t caliente enfriarse a una temperatura de 1200 F?.
6. En una fra maana de Bogot con una temperatura ambiente de 100C, se encontr
un cadver a las 10 : 00 a.m.. El detective sac un termmetro y midi la temperatura
del cuerpo: 330C. Luego sali a revisar la escena del crimen. Al regresar, a las 11 : 00
a.m., hall que la temperatura del cuerpo era de 320C. En qu momento ocurri el
asesinato?. (Sugerencia: la temperatura normal del cuerpo es de 370C.)
7. En una fresca maana de sbado, mientras las personas trabajan en el interior, el
calefactor mantiene la temperatura interior del edificio en 210 C. A medioda, el
aparato se apaga y los empleados se van a casa. La temperatura exterior es constante
e igual a 120 C durante el resto de la tarde. Si la constante de tiempo para el edificio es
de 3 horas, en qu momento llegar la temperatura del edificio a 160 C?. Si algunas
ventanas se dejan abiertas y la constante de tiempo se reduce a 2 horas, en qu
momento llegar la temperatura interior a 160 C?.
8. Una cochera sin calefaccin ni aire acondicionado tiene una constante de tiempo de 2
horas. Si la temperatura exterior vara como una onda sinusoidal con un mnimo de
500 F a las 2 : 00 a.m. y un mximo de 800 F a las 2 : 00 p.m., determine los instantes en
que el edificio alcanza su temperatura mxima y mnima, suponiendo que el trmino
exponencial se extingue.
9. Un lunes temprano por la maana, la temperatura en la sala de lectura ha descendido
hasta 400 F, igual a la temperatura exterior. A las 7 : 00 a.m., el conserje enciende el
calefactor con el termostato puesto en 700 F. La constante de tiempo para el edificio
80
es 1/K = 2 horas y la constante de tiempo para el edificio junto con su sistema de

calentamiento es 1/K1 = 1/2 hora. Suponiendo que la temperatura exterior permanece
constante, cul ser la temperatura dentro de la sala de lectura a las 8 : 00 a.m.?. En
qu momento llegar la temperatura dentro de la sala a 650 F?.
Solucin
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.
3.4
Mecnica de Newton
La mecnica es el estudio del movimiento de los objetos y el efecto de las fuerzas que
actan sobre tales objetos. la mecnica de Newton o clsica, trata del movimiento de los
objetos comunes; es decir, de los objetos que son grandes en comparacin con un tomo y
lentos en comparacin con la velocidad de la luz.
Un modelo de la mecnica de Newton se puede basar en las leyes del movimiento de
Newton:
(1.) En ausencia de fuerzas, un objeto (cuerpo) en descanso seguir en descanso, y un
cuerpo movindose a una velocidad constante en lnea recta, lo continuar haciendo
indefinidamente.
(2.) Cuando un cuerpo se sujeta a una o ms fuerzas externas, la razn de cambio temporal
del momento del cuerpo es igual a la suma vectorial de las fuerzas externas que actan
sobre l.
(3.) Cuando un cuerpo interacta con otro, la fuerza del primer cuerpo sobre el segundo
es igual en magnitud, pero opuesta en direccin, a la fuerza del segundo cuerpo sobre
el primero.
La segunda ley de Newton que solo se aplica a marcos de referencia inerciales, nos permite
formular las ecuaciones de movimiento para un cuerpo. Podemos expresar la segunda ley
de Newton como
dp
= F(t, x, v),
dt
donde F(t, x, v) es la fuerza resultante sobre el cuerpo en el instante t, posicin x, velocidad
v y p(t) es el momentum del cuerpo en el instante t. El momentum es el producto de la masa
3.4 Mecnica de Newton
81
del cuerpo y su velocidad, es decir,

p(t) = mv(t)
de modo que podemos expresar la segunda ley de Newton como
m
donde a =
dv
dt
dv
= ma = F(t, x, v)
dt
es la aceleracin del cuerpo en el instante t.
En esta seccin nos centraremos en situaciones donde la fuerza F no depende de x. Luego

la ecuacin de primer orden en v(t) es:
m
dv
= ma = F(t, v).
dt
(3.14)
Ejemplo 3.12. Un objeto de masa 5 kg recibe una velocidad inicial hacia abajo de 50 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 10v, donde v es la velocidad del objeto en m/seg.
Determine la ecuacin de movimiento del objeto. Si el objeto est inicialmente a 500 m
sobre el suelo, determine el momento en que el objeto golpear el suelo. (Figura 3.2)
F2 = bv(t)
t=0
x(t)
b
F1 = mg
b
Figura 3.2: Cada libre

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:
82
La fuerza debida a la gravedad es
F1 = mg
donde g es la aceleracin debida a la gravedad en la superficie de la tierra.
La fuerza debida a la resistencia del aire es
F2 = bv(t)
donde b > 0 es la constante de proporcionalidad y el signo negativo est presente
debido a que la resistencia del aire acta en forma opuesta al movimiento del objeto.
La fuerza neta que acta en el objeto es:
F = F1 + F2 = mg bv(t).
(3.15)
Ahora utilizando la ecuacin (3.14) y sustituyendo en (3.15) se obtiene:

m
dv
= mg bv,
dt
v(0) = v0
(3.16)
m
m
donde m = 5 kg, v0 = v(0) = 50 seg
, b = 10, g = 9,81 seg
2 . Reemplazando en (3.16) obtenemos:
m
dv
= 5kg 9,81 2 10v(t)N
dt
seg
5kg
m
10N
=
9,81 2
v(t)
5kg
seg
5kg
m
m
= 9,81 2 2 2 v(t)
seg
seg
5kg
dv
dt
dv
dt
dv
+ 2v(t) = 9,81
dt
utilizando variables separables o lineal obtenemos:

v(t) = 4,905 + ce2t ,
como v(0) = 50 se obtiene, c = 45,095. Entonces
v(t) = 4,905 + 45,095e2t .
Para hallar la ecuacin de movimiento primero resolvemos la ecuacin diferencial
v(t) =
dx
= 4,905 + 45,095e2t
dt
donde
x(t) = 4,905t 22,55e2t + k.
(3.17)
3.4 Mecnica de Newton
83
Como hemos considerado que x = 0 cuando t = 0, reemplazando en (3.17) obtenemos

k = 22,55. Por lo tanto, la ecuacin de movimiento es:
x(t) = 4,905t 22,55e2t + 22,55.
Como 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
x(t) = 4,905t 22,55e2t + 22,55
o
t 4,59e2t = 97,3.
(3.18)
La ecuacin (3.18) no se puede resolver de manera explcita en trminos de t. Para resolverla

tomamos el hecho que e2t ser muy pequeo para t, es decir ignoramos el trmino e2t y
obtenemos como aproximacin t 97,3 segundos.

Ejemplo 3.13. Un paracaidista cuya masa es de 75 kg se arroja de un helicptero que vuela
a 2.000 m sobre el suelo y cae hacia ste 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 b1 = 30 N.seg/m cuando se abre, si el paracadas est
cerrado, y b2 = 90 N.seg/m cuando est abierto. Si el paracadas no se abre hasta que la
velocidad del paracaidista es de 20 m/seg. Despus de cuntos segundos llegar al suelo?.
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
x(0) = 0.
Reemplazando en (3.14) tenemos
dv
= 75 9,81 30v(t)
dt
dv
m
= 9,81 2 0,4v(t)
dt
seg
dv
+ 0,4v(t) = 9,81.
dt
75
Utilizando variables separables, o el hecho de que esta ecuacin es lineal, resolvemos y

obtenemos:
v1 (t) = 24,53 + ce0,4t
como v(0) = 0 se obtiene, c = 24,53. As v1 (t) = 24,53 24,53e0,4t y
x1 (t) = 24,53t1 + 61,31e0,4t1 61,31.
84
Como el paracadas se abre cuando se alcance la velocidad de 20 m/seg, es decir la velocidad

final. Se puede hallar el tiempo que recorri. Por lo tanto
20 = 24,53 24,53e0,4t1
t1 = 4,22seg
y ha cado
x1 (4,22) = 24,53(4,22) + 61,31e0,4(4,22) 61,31 = 53,54
(en estos y otros clculos, redondearemos nuestras cifras a dos cifras decimales). Ahora, al
abrirse el paracadas, el paracaidista est a 2.000-53,54=1.946,46 metros sobre el suelo.
m
Consideramos el caso despus de abrirse: m = 75 kg, v0 = v(0) = 20 seg
, b2 = 90N
m
seg/m, g = 9,81 seg2 y x2 (0) = 0, por lo tanto
x2 (t) = 8,12t2 + 9,85 9,85e1,2t2 .
Para determinar el momento en que el paracaidista llega al suelo, hacemos x2 (t) =1.946,46
m, la altura a la que estaba el paracaidista al abrir su paracadas. Esto da como resultado
1.946,46 = 8,12t2 + 9,85 9,85e1,2t2
238,5 = t2 1,21e1,2t2 .
(3.19)
La ecuacin (3.19) no se puede resolver de manera explcita en trminos de t. Para resolverla

tomamos el hecho de que e1,2t ser muy pequeo para t, es decir, ignoramos el trmino
e1,2t y obtenemos como aproximacin t2 238,5 segundos. Por lo tanto, golpear el suelo
238,5 seg despus de arrojarse desde el helicptero.

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,
determine 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.
Determine la ecuacin de movimiento del objeto. Si el objeto est inicialmente a 100
m sobre el suelo, determine 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
3.5 Trayectorias Ortogonales
4.
5.
6.
7.
8.
3.5
85
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. Determine la ecuacin de
movimiento del objeto. Cul es la velocidad del objeto al minuto despus de ser
arrojado?.
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?.
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. Determine la ecuacin de
movimiento del objeto. Despus de cuantos segundos ocurrir que la velocidad del
objeto es igual a 70 m/seg?.
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, determine la velocidad lmite del objeto (suponiendo que este
lmite existe). (Sugerencia: Justifique que la existencia de una velocidad lmite (finita)
implica dv
dt 0 cuando t +).
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
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?.
Si un objeto cae desde un aeroplano la velocidad hacia abajo despus de x seg es
m
1 . y(x) es la
aproximadamente y0 = gk (1 ekx ) en donde g = 9,8 seg
2 y k = 0, 2seg
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?.
Trayectorias Ortogonales
Dos curvas se cortan en ngulo recto si las respectivas rectas tangentes, en los puntos de
interseccin, son rectas perpendiculares.
Definicin 3.1. Cuando todas las curvas de una familia F(x, y, c1 ) = 0 cortan ortogonalmente a todas las curvas de otra familia G(x, y, c2 ) = 0, se dice que las familias son ortogonales.
86
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.
Figura 3.3: Crculos concntricos y las rectas y = c1 x.

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

F/ x
F
r(x, y) = F/ y , en todo punto (x, y) de F(x, y, c1 ) = 0 donde y 6= 0 .
La ecuacin de la familia ortogonal es dada, en forma diferencial, por la ecuacin 2
dy
1
=
.
dx
r(x, y)
Al eliminar la constante c1 entre el par de ecuaciones
F(x, y, c1 ) = 0
dy
1
=
dx
r(x, y)
2 Esta
expresin est inspirada en el hecho que dos rectas son perpendiculares si el producto de sus
pendientes es 1.
87
y resolver la ecuacin diferencial resultante (si es que es posible), obtenemos una representacin algebraica de la familia ortogonal a la familia F(x, y, c1 ) = 0.
Ejemplo 3.15. Hallar las trayectorias ortogonales de la familia de curvas y = ekx .
Solucin. La derivada de y = ekx es:
dy
= kekx .
dx
Despejamos k de y = ekx , ln(y) = kx, k =
ln(y)
x ,
x 6= 0 y reemplazando se obtiene:
dy ln(y) ln(y)
=
e
dx
x
dy y ln(y)
=
.
dx
x
La ecuacin diferencial de la familia ortogonal es
dy
x
=
dx
y ln(y)
esta ltima ecuacin se resuelve por variables separables:
y ln(y)dy = xdx
Z
y ln(y)dy =
xdx
x2
y2 ln(y) y2
= + c2
2
4
2
2
2
2
2y ln(y) y + 2x = 4c2
2y2 ln(y) y2 + 2x2 = c3
en donde por conveniencia se reemplaz 4c2 por c3 .
Ejemplo 3.16. Hallar las trayectorias ortogonales a las curvas

y = cxn ,
n R.
(3.20)
Solucin. La derivada de y = cxn es y = cnxn1 . Despejamos c de la ecuacin (3.20) y

tenemos c = yxn , x, y 6= 0. Reemplazando tenemos
dy
= xn nyxn = nx1 y.
dx
88
La ecuacin diferencial de la famila de curvas es

por variable separables, de donde tenemos:
dy
dx
n
n
x
= ny
. Esta ltima ecuacin se resuelve
nydy = xdx
ydy =
y2
x2
= + c.
2
2
xdx
La trayectorias ortogonales vienen dada por la frmula
ny2
2
+ x2 = c.
Figura 3.4: Las curvas y = cx2 para c = 1, 2 y sus trayectorias ortogonales

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 .
89
2. Determine 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).
Solucin
1. a)
b)
d)
e)
i)
j)
3. a)
b)
c)
x2 2xy y2 = c.
x2 + y2 = cx.
2x2 + y2 = c.
x3 + 3y = c.
ey sin(y) = c.
x = cos2 (y).
x2 + y2 = c.
x2 + y2 = cy.
2
2
x2 = cex +y .
4 Ecuaciones Diferenciales Lineales de Orden Superior
En este captulo presentamos algunos mtodos de solucin de ecuaciones diferenciales

lineales de orden mayor o igual a dos.
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
Si g(x) = 0 se dice que es homognea.
Si g(x) 6= 0 se dice que es no homognea.
Si las funciones ai (x) para i = 1, 2, , n son constantes; es decir, ai (x) ai R para
0 i n, entonces se dice que la ecuacin lineal es de coeficientes constantes.
Estudiaremos los mtodos ms conocidos y de mayor utilidad para hallar soluciones de
ecuaciones diferenciales como las antes mencionadas. El caso ms sencillo de ecuaciones
diferenciales lineales de orden superior, es el de ecuaciones diferenciales homogneas con
coeficientes constantes. El mtodo de solucin de este tipo de ecuaciones se reduce a la
factorizacin de polinonios de orden n con los coeficientes dados en la ecuacin diferencial.
El caso de ecuaciones diferenciales lineales con coeficientes constantes no homogneas se
resolver primero mediante el llamado mtodo de los coeficientes indeterminados el cual es
aplicable slo para ciertos trminos no homogneos, y luego usando el mtodo de reduccin
de orden y variacin de parmetros los cuales, a diferencia del mtodo de coeficientes
indeterminados, pueden ser usados para hallar soluciones de ecuaciones diferenciales
lineales no homogneas con coeficientes no constantes. A diferencia de las ecuaciones
homogneas con coeficientes constantes, las soluciones de otras clases de ecuaciones
dependern de al menos una solucin conocida de dicha ecuacin (o una solucin de la
ecuacin homognea asociada), la construccin de la solucin general en este caso se
Ecuaciones Diferenciales Lineales de Orden Superior
92
justificar mediante el llamado Principio de Superposicin.
4.1
Teora Preliminar
Antes de mostrar las tcnicas de solucin de ecuaciones lineales de orden superior,
daremos los resultados tericos que garantizan que dichas tcnicas son aplicables, y que las
eventuales soluciones son nicas.
Teorema 4.1 Existencia y unicidad de soluciones para ecuaciones lineales de
orden n. Dada la ecuacin diferencial lineal de orden n:
an (x)
d n1 y
dy
dny
+
a
(x)
+ + a1 (x) + a0 (x)y = g(x).
n1
n
n1
dx
dx
dx
Supongamos que g(x), ai (x), con i = 0, 1, 2, . . . , n, son funciones continuas en un intervalo

(a, b), que contiene al punto x0 y an (x) 6= 0 en (a, b). Entonces, para cualquier eleccin de
los valores iniciales y0 , y1 , . . . , yn1 , existe una nica solucin y(x) en todo el intervalo (a, b)
del problema con valor inicial
y(x0 ) = y0 , y0 (x0 ) = y1 , . . . , y(n1) (x0 ) = yn1 .
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).
Ejemplo 4.1. Considere la ecuacin diferencial

y00 y = 0.
(a.) Muestre que todas las funciones de la forma y(x) = ex son soluciones de dicha
ecuacin.
(b.) Muestre que todas las funciones de la forma y(x) = ex son soluciones de dicha
ecuacin.
(c.) Muestre que todas las funciones de la forma y(x) = 1 ex + 2 ex son soluciones de
dicha ecuacin.
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 prueba (a.). Para y(x) = ex , tenemos que y00 (x) = ex y as esta
funcin tambin es una solucin de la ecuacin, de donde tenemos (b.). Para mostrar (c.)
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.
4.1 Teora Preliminar
93
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 tambin es solucin de la
ecuacin. Lo anterior nos lleva a a preguntar cuando dos (o ms) soluciones de una ecuacin
diferencial son realmente diferentes, para ello usaremos las siguientes nociones.
Definicin 4.1 (Dependencia Lineal). Un conjunto de funciones { f1 , f2 , . . . , fn } es linealmente dependiente en un intervalo I si existen constantes c1 , c2 , . . . , cn , no todas iguales a 0,
tales que
c1 f1 + c2 f2 + + cn fn 0,
es decir
c1 f1 (x) + c2 f2 (x) + + cn fn (x) = 0
para todo x en el intervalo I. Por otro lado, si todas las constantes son iguales a cero, se dice
que { f1 , f2 , . . . , fn } es un conjunto linealmente independiente.
Definicin 4.2 (Wronskiano). Dadas las funciones f1 (x), f2 (x), . . . , fn (x) se define su
Wronskiano como el determinante:

f1 (x)
f2 (x)
fn (x)
0
f (x)
f20 (x)
fn0 (x)
1
W [ f1 , f2 , , fn ](x) =
..
..
..
.
...

.
.
.
(n1)

(n1)
(n1)
f

(x)
f
(x)
f
(x)
n
1
2
El Wronskiano es una herramienta muy til para determinar si un conjunto de funciones es
linealmente independiente o linealmente dependiente.
Si un conjunto de funciones f1 , . . . , fn continuas con derivadas de orden

n continuas es linealmente independiente en un intervalo (a, b), entonces el Wronskiano es
idnticamente cero en en (a, b).
Proposicin 4.2
La proposicin anterior nos da condiciones necesarias para la dependencia lineal de un

conjunto de funciones, pero esta condicin no es suficiente. Si el Wronskiano de un conjunto
de funciones es idnticamente cero en (a, b) no podemos concluir que el conjunto es
linealmente dependiente.
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|
= x3 |x|0 + x2 |x| 2x2 |x|.
W [ f1 , f2 ](x) =
0
2x x|x| + |x|
Cuando x > 0, |x| = x y |x|0 = x0 = 1. As, el Wronskiano queda x3 + x3 2x3 = 0. Cuando

x = 0, |x| = x = 0 y de nuevo el Wronskiano es 0.
94
Para x < 0, |x| = x y |x|0 = x0 = 1. El Wronskiano en este caso es x3 (1)x3 +2x3 = 0.

Por lo tanto, el Wronskiano es idnticamente cero en un intervalo que contiene al punto x = 0.
Sin embargo, estas funciones son linealmente independientes en I. En efecto, consideremos
la combinacin lineal
c1 x2 + c2 |x| = 0.
Si x = 1 entonces c1 = c2 y si x = 1 se tiene que c1 = c2 , luego la nica solucin es
c1 = c2 = 0. Es decir, {x2 , x|x|} es un conjunto linealmente independiente en I.

Condiciones necesarias y suficientes para que un conjunto de funciones { f1 , . . . , fn } sea
linealmente independiente en I viene dado mediante el llamado Grammiano del conjunto
{ f1 , . . . , fn }. Ver pgina 98.
Para un conjunto de funciones que son soluciones de una ecuacin diferencial lineal de
orden n homognea, el Wronskiano si aporta condiciones necesarias y suficientes para
independencia lineal de estas funciones:
Un conjunto de soluciones {y1 , y2 , . . . , yn } de la ecuacin diferencial
lineal de orden n homognea
Proposicin 4.3
an (x)
dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y = 0
n1
n
n1
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,
a2 (x)y00 + a1 (x)y0 + a0 (x)y = 0
el Wronskiano se puede calcular mediante la llamada frmula de Abel:
W [y1 , y2 ](x) = Ce
R a1 (x)
a2 (x)
dx
C constante.
Ms aun, tenemos lo siguiente.

Si cada funcin y1 , y2 , . . . , yn es una solucin de la ecuacin diferencial
lineal de orden n homognea
Proposicin 4.4
an (x)
dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y = 0
n1
n
n1
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.
95
Definicin 4.3 (Conjunto fundamental de soluciones). Las soluciones y1 , y2 , . . . , yn de la

ecuacin diferencial
an (x)
dny
d n1 y
dy
+ a0 (x)y = 0,
+
a
(x)
+
+
a
(x)
n1
1
dxn
dxn1
dx
forman un conjunto fundamental de soluciones en I si

W [y1 , y2 , , yn ](x) 6= 0.
El hecho de tener un conjunto fundamental de soluciones es importante por lo siguiente:
Proposicin 4.5
cin diferencial
Si {y1 , y2 , . . . , yn } en un conjunto fundamental de soluciones de la ecuadny

d n1 y
dy
an (x) n + an1 (x) n1 + + a1 (x) + a0 (x)y = 0,
dx
dx
dx
y si y es cualquier otra solucin de esta ecuacin, entonces existen constantes c1 , c2 , . . . , cn

tales que
y = c1 y1 + c2 y2 + + cn yn .
Es decir, un conjunto fundamental de soluciones es realmente una base para el conjunto de
todas las soluciones de una ecuacin lineal homognea de orden n. De ah que
y = c1 y1 + c2 y2 + + cn yn
reciba el nombre de solucin general de tal ecuacin.
Corolario 4.6
Sean y1 , y2 , . . . , yn soluciones en I de
an (x)
dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y = 0,
n1
n
n1
dx
dx
dx
entonces las siguientes afirmaciones son equivalentes:

(i.) {y1 , y2 , . . . , yn } es un conjunto fundamental de soluciones en I.
(ii.) y1 , y2 , . . . , yn son linealmente independientes en I.
(iii.) El Wronskiano W [y1 , y2 , . . . , yn ](x) 6= 0 en I.
Ejemplo 4.3. Determine si y1 (x) = ex cos(2x),

dependientes en el intervalo (0, 1).
y2 (x) = ex sin(2x) son linealmente
Solucin. Para determinar si son linealmente independientes hallamos el Wronskiano de y1

y y2

ex cos(2x)
ex sin(2x)

= 2e2x .
W [y1 , y2 ](x) = x
x
x
x
e cos(2x) 2e sin(2x) e sin(2x) + 2e cos(2x)
Donde 2e2x es diferente de cero en todo el conjunto de los nmeros reales, por lo tanto
{y1 , y2 } es un conjunto linealmente independiente.
96
Ahora consideraremos el caso de diferentes soluciones de una ecuacin diferencial

con distintos 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
d n1 y
dy
dny
+
a
(x)
+ + a1 (x) + a0 (x)y = g1 (x)
n1
n
n1
dx
dx
dx
y sea y2 una solucin de la ecuacin diferencial
an (x)
an (x)
dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y = g2 (x),
n1
n
n1
dx
dx
dx
entonces para cualquier par de constantes c1 y c2 , la funcin c1 y1 + c2 y2 es una solucin de

la ecuacin diferencial
an (x)
dny
d n1 y
dy
+ a0 (x)y = c1 g1 (x) + c2 g2 (x).
+
a
(x)
+
+
a
(x)
n1
1
dxn
dxn1
dx
Como consecuencia de lo anterior tenemos el siguiente resultado.

Corolario 4.8
La solucin general de la ecuacin diferencial

an (x)
dny
d n1 y
dy
+ a0 (x)y = g(x)
+
a
(x)
+
+
a
(x)
n1
1
dxn
dxn1
dx
se puede escribir como

y(x) = c1 y1 (x) + + cn yn (x) +Y (x)
donde c1 y1 + + cn yn es la solucin general de la ecuacin homognea correspondiente
an (x)
d n1 y
dy
dny
+
a
(x)
+ + a1 (x) + a0 (x)y = 0
n1
n
n1
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:
(1.) Encontramos la solucin general yh de la ecuacin homognea
an (x)
d n1 y
dy
dny
+
a
(x)
+ + a1 (x) + a0 (x)y = 0 y despus
n1
n
n1
dx
dx
dx
97
(2.) encontramos una solucin particular y p de la ecuacin no homognea

dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y = g(x).
n1
n
n1
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 .
an (x)
Ejemplo 4.4. Dada la ecuacin diferencial y00 2y0 + 5y = 0

(a.) Determine si y1 (x) = ex cos(2x) y y2 (x) = ex sin(2x) son soluciones linealmente independientes de la ecuacin dada.
(b.) Determine una solucin general de la ecuacin diferencial.
(c.) Determine la solucin que satisface las condiciones iniciales y(0) = 2,
Solucin:
y0 (0) = 0.
(a.) Hallamos el Wronskiano de y1 y y2

x cos(2x)
x sin(2x)

e
e
= 2e2x .
W [y1 , y2 ](x) = x
x
x
x
e cos(2x) 2e sin(2x) e sin(2x) + 2e cos(2x)
Como 2e2x es diferente de cero para todos los reales, entonces y1 , y2 son linealmente
independientes.
(b.) 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.
Es decir, y1 y y2 son soluciones de la ecuacin diferencial, las cuales son linealmente
independientes en todos los reales (ninguna es mltiplo escalar de la otra en R). Por
lo tanto, por el Corolario 4.6, {ex cos(2x), ex sin(2x)} es un conjunto fundamental de
soluciones en R, de modo que una solucin general es
y(x) = c1 ex cos(2x) + c2 ex sin(2x).
(c.) 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 .
De donde, c1 = 2 y c2 = 1. Por lo tanto la solucin que satisface las condiciones
iniciales es:
y(x) = 2ex cos(2x) 1ex sin(2x).
98
Ejercicios
1. En los problemas a) a d), cada familia de funciones es la solucin general de la

ecuacin diferencial en el intervalo indicado. Determine 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 (, ), demuestre 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 contradice el Teorema 4.1.
3. Determine a dos miembros de la familia de soluciones en el ejercicio anterior que
satisfagan las condiciones iniciales y(0) = 0, y0 (0) = 0.
4. En los ejercicios a) y b) determine un intervalo con centro x = 0, para el cual el
problema de valor inicial tenga una solucin nica.
a) (x 2)y00 + 3y = x, y(0) = 0, y0 (0) = 1.
b) y00 + (tan(x))y = ex , y(0) = 1, y0 (0) = 0.
5. Use el Teorema 4.1 para analizar la existencia y unicidad de una solucin de la
ecuacin diferencial que satisfaga las condiciones iniciales y(1) = y0 , y0 (1) = y1 ,
donde y0 y y1 son constantes reales.
a) (1 + x2 )y00 + xy0 y = tan(x).
y0
+ y = ln(x).
b) ex y00 x3
6. Determine si se aplica el Teorema 4.1. En caso afirmativo, analice las conclusiones
que pueden extraerse. En caso negativo, explique por qu.
a) y00 + yy0 = x2 1; y(0) = 1, y0 (0) = 1.
b) (1 t)x00 + tx0 2 = s sin(t); x(0) = 1, x0 (0) = 1.
7. Sabiendo que y(x) = e2x es una solucin de (2x+1)y00 4(x+1)y0 +4y = 0, encontrar
la otra solucin linealmente independiente. Escribir la solucin general.
8. Grammiano: Un conjunto de funciones { f1 , . . . , fn } cada una continua en un intervalo I = [a, b] es linealmente dependiente en I si, y slo si,
Rb 2

Rb
Rb

f
(x)dx
f
(x)
f
(x)dx
.
.
.
f
(x)
f
(x)dx
n
1
2
1
a
a
a
1
R

Rb 2
Rb
b f (x) f (x)dx
. . . a f2 (x) fn (x)dx
1
a 2
a f 2 (x)dx
G[ f1 , . . . , fn ] =
= 0,
..
..
...

.
.
R

Rb 2
b f (x) f (x)dx R b f (x) f (x)dx . . .
f (x)dx
a
en I. Este determinante se llama el Grammiano del conjunto { f1 , . . . , fn }.

a) Demuestre que ex y e2x son linealmente independientes en [1, 1].
b) Demuestre que x y |x| son linealmente independientes en [1, 1].
9. Determine si las funciones y1 y y2 son linealmente independientes en el intervalo
(0, 1). Calcule adems el Wronskiano W [y1 , y2 ](x):
99
a) y1 (x) = ex cos(2x), y2 (x) = ex sin(2x).

b) y1 (x) = xe2x , y2 (x) = e2x .
c) y1 (x) = x2 cos(ln x), y2 (x) = x2 sin(ln(x)).
d) y1 (x) = 0, y2 (x) = ex .
10. (a.) Verifique si las funciones y1 y y2 son linealmente independientes de la ecuacin
diferencial dada. (b.) Determine una solucin general de la ecuacin dada y (c.)
determine la solucin que satisface las condiciones iniciales.
a) x2 y00 2y = 0; y1 (x) = x2 , y2 (x) = x1 ; y(1) = 2; y0 (1) = 7.
b) y00 5y0 + 6y = 0; y1 (x) = e2x , y2 (x) = e3x ; y(0) = 1; y0 (0) = 4.
c) ty00 (t +2)y0 +2y = 0; y1 (t) = et , y2 (t) = t 2 +2t +2; y(1) = 0; y0 (1) =
1.
11. Sean y1 (x) = x3 y y2 (x) = |x3 |. Son y1 y y2 linealmente independientes en los
siguientes intervalos?
a) [0, ).
b) (, 0).
c) (, ).
d) Calcule el Wronskiano W [y1 , y2 ](x) en el intervalo (, ).
12. Encontrar las regiones donde el Teorema de existencia y unicidad garantiza la unicidad
para problemas de Cauchy:
y00 + 3xy0 + x3 y = ex ;
y(0) = 0,
y0 (0) = 0.
13. Dado que y1 (x) = cos(x) es una solucin de

y00 y0 + y = sin(x),
y y2 (x) = e2x /3 es una solucin de
y00 y0 + y = e2x ,
determine soluciones de las siguientes ecuaciones diferenciales:
a) y00 y0 + y = 5 sin(x).
b) y00 y0 + y = sin(x) 3e2x .
c) y00 y0 + y = 4 sin(x) + 18e2x .
14. dado que 1+x, 1+2x y 1+3x2 son soluciones de la ecuacin diferencial y00 + p(x)y0 +
q(x)y = g(x), determine la solucin de esta ecuacin que satisface y(1) = 1, y0 (1) = 0.
Solucin
1. b)
c)
d)
4. a)
9. b)
c)
y(x) = 33 e4x + 52 ex .
y(x) = 3x 4x ln(x).
y(x) = 1 + cos(x) 2 sin(x).
(, 2) o (2, ).
Linealmente independiente.
Linealmente independiente.
100
10.
11.
13.
4.2
d)
a)
c)
b)
d)
b)
c)
Linealmente dependiente.
y(x) = 3x2 + x1 .
y(x) = 5et1 t 2 2t 2.
No son linealmente independientes.
0.
y(x) = cos(x) e2x .
y(x) = 4 cos(x) + 6e2x .
Ecuaciones lineales Homogneas con Coeficientes Constantes

En esta seccin mostraremos la tcnica para resolver ecuaciones diferenciales ordinarias de
la forma:
an yn + an1 yn1 + + a1 y0 + a0 y = 0, con a0 , a1 , . . . , an R.
(4.1)
Sabemos que la solucin de una ecuacin diferencial de primer orden de la forma y0 + ay = 0

es: y(x) = ceax , entonces podemos generalizar este resultado para la ecuacin (4.1). Para
resolver esta ecuacin es necesario determinar n soluciones linealmente independientes.
Comenzaremos suponiendo que la ecuacin tiene una solucin de la forma y(x) = emx . El
problema se reduce a encontrar los valores del exponente m que hace que las funciones
y(x) = emx sean soluciones de la ecuacin. As,
y(x) = emx , y0 (x) = memx , y00 (x) = m2 emx , . . . , yn (x) = mn emx .
Sustituyendo en la ecuacin tenemos:
an mn emx + an1 mn1 emx + + a1 memx + a0 emx = 0,
(an mn + an1 mn1 + + a1 memx + a0 )emx = 0,

an mn + an1 mn1 + + a1 m + a0 = 0.
Esta ltima ecuacin es llamada ecuacin auxiliar (asociada a la ecuacin diferencial). Los
exponentes buscados son las races del polinomio caracterstico
p(m) = an mn + an1 mn1 + + a1 m + a0
del cual podemos tener varios casos:
El polinomio caracterstico tiene n races reales simples y distintas: m1 , m2 , . . . , mn .
Por lo tanto, tenemos n soluciones fundamentales de la ecuacin:
y1 (x) = em1 x , y2 (x) = em2 x , . . . , yn (x) = emn x .
Puede verificarse que el Wronskiano es no nulo para todo x R, y por la Definicin
4.3 se puede concluir que es un conjunto fundamental de soluciones y la solucin
general de la ecuacin es:
y(x) = C1 em1 x +C2 em2 x + +Cn emn x .
4.2 Ecuaciones lineales Homogneas con Coeficientes Constantes
101
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 otra
utilizamos el mtodo de reduccin de orden, con la frmula,
R e R p(x)dx
2
y2 (x) = y1 (x) [y (x)]2 dx en este caso p(x) = ba y m = b
2a ya que b 4ac = 0
1
sustituyendo en la ecuacin obtenemos:
mx
y2 (x) = e
R b
a dx
b
e 2a x
2 dx
resolviendo la integral obtenemos y2 (x) = xemx estas dos soluciones son linealmente
independientes y la solucin general de la ecuacin es:
y(x) = C1 emx +C2 xemx
y para la ecuacin (4.1), si la multiplicidad de la raz es n se obtienen n soluciones
linealmente independientes. As, la solucin general viene dada por:
y(x) = C1 emx +C2 xemx + +Cn xn1 emx .
El polinomio tiene un par de races complejas conjugadas :
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:
e(a+ib)x = eax eibx = eax [cos(bx) + i sin(bx)].
De esta forma podemos escribir las soluciones de la siguiente manera:
C1 ex [cos(bx) + i sin(bx)] +C2 ex [cos(bx) + i sin(bx)]
= C1 ex [cos(bx) + i sin(bx)] +C2 ex [cos(bx) i sin(bx)]
= ex [(C1 +C2 ) cos(bx) + (iC1 iC2 ) sin(bx)]
= ex [K1 cos(bx) + K2 sin(bx)].
Si la multiplicidad de la raz compleja es r la solucin general es:
y =ex [K1 cos(bx) + K2 sin(bx) + x(K3 cos(bx) + K4 sin(bx))
+ + xr1 (K2r1 cos(bx) + K2r sin(bx))].
Ejemplo 4.5. Determine la solucin general de la ecuacin diferencial:
(a.) y000 6y00 y0 + 6y = 0.
102
(b.) 4y00 + 20y0 + 25y = 0.

(c.) y000 y00 + y0 + 3y = 0.
Solucin: (a.) 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,
factorizando la ecuacin anterior obtenemos (m 6)(m 1)(m + 1) = 0. Por lo
tanto las races de la ecuacin auxiliar son 6,1 y -1, de modo que tres soluciones
fundamentales de la ecuacin son e6x , ex y ex , y como sabemos estas funciones
forman un conjunto fundamental de soluciones, por lo tanto solucin general es:
y(x) = c1 e6x + c2 ex + c3 ex .
(b.) 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.
La ecuacin anterior se reescribe como (2m + 5)(2m + 5) = 0. De donde se tiene
que 52 es una raz con multiplicidad dos, de modo que se tiene dos soluciones
5
5
fundamentales las cuales son e 2 x y xe 2 x . Luego la solucin general es:
5
y(x) = c1 e 2 x + c2 xe 2 x .
(c.) La ecuacin auxiliar en este caso es
m3 m2 + m + 3 = 0
la cual se factoriza como (m + 1)(m2 2m + 3) = 0. Por lo tanto las races de la
ecuacin auxiliar son:
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 Mtodo de Coeficientes Indeterminados
103
Ejercicios
1. Resolver las siguientes ecuaciones diferenciales homogneas:

a) 3y00 + 11y0 7y = 0.
b) y000 + y00 6y0 + y = 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. Examine 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
y00 + by0 + 4y = 0; y(0) = 1, y0 (0) = 0,
resuelva el problema para b = 5,4 y 2 y bosqueje las soluciones.
Solucin
1. b) y(x) = c1 e(1,8919)x+
c2
e(0,1724)x + c3 e(3,0644)x
.

x
x
5
5
d) y(x) = c1 e 2 cos 2 x + c2 e 2 sin 2 x .

g)
i)
k)
m)
n)
o)
4.3
y(x) = c1 ex + c2 xex + c3 x2 ex .
y(x) = 15 e3x + 15 e2x + ex .
y(x) = 12 sin(4x) + 2 cos(4x).
y(x) = c1 ex + c2 xex + c3 ex + c4 xex .
y(x) = 0.
y(x) = ex cos(2x) + (e 21) ex sin(2x).
Mtodo de Coeficientes Indeterminados

En esta seccin y las siguientes nos dedicaremos a hallar una solucin particular de una
ecuacin diferencial lineal de orden superior no homognea. Primero estudiaremos el
mtodo de los coeficientes indeterminados que nos permite hallar de una forma directa
104
una solucin particular de una ecuacin diferencial lineal no homognea con coeficientes
constantes. La tcnica consiste en asumir una forma particular de dicha solucin y p (x) la
cual depender directamente de la forma del trmino no homogneo; por ejemplo, si el
trmino no homogneo es g(x), la forma de la correspondiente solucin particular podra ser
y p (x) = g(x). As, al sustituir g(x) en la ecuacin obtendremos un sistema de ecuaciones
algebraico que nos permitir determinar .
A modo de ejemplo, consideremos la ecuacin
y00 + 2y0 + y = ex ,
dado que g(x) = ex , entonces la solucin particular 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 ,
de donde concluimos que el coeficiente es = 14 , y por tanto la solucin particular de la
ecuacin es
1
y p (x) = ex .
4
Debido a la naturaleza del mtodo, este es til en el caso de coeficientes constantes y
trminos no homogneos particulares. As que en este caso solo consideraremos trminos
no homogneos g(x) de los de la forma dada en la siguiente tabla.
g(x)
xn + + a
1.
2.
3.
4.
5.
pn (x) = an
1 x + a0
x
ae
a cos( x) + b sin( x)
pn (x)ex
pn (x) cos( x) + qm (x) sin( x), donde, qm (x) = bm xm + + b1 x + b0 .
6.
7.
aex cos( x) + bex sin( x)

ex [pn (x) cos( x) + qm (x) sin( x)]
y p (x)
s
s
n
x Pn (x) = x (An x + + A1 x + A0 )
xs Aex
xs (A cos( x) + B sin( x))
xs (An xn + + A1 x + A0 )ex
xs {PN (x) cos( x) + QN (x) sin( x)}, donde QN (x) = BN xN + + B1 x + B0 y N =
max(m, n).
xs (Aex cos( x) + Bex sin( x))
xs {ex (PN (x) cos( x) + QN (x) sin( x))},
y N = max(m, n).
Cuadro 4.1: Coeficiente no homogneo para el mtodo de los coeficientes indeterminados y

la correspondiente forma de la solucin.
Donde el entero no negativo s es el menor entero tal que ningn trmino de la solucin y p
sea solucin de la ecuacin homognea.
Ejemplo 4.6. Determinar una solucin general de la ecuacin:
105
(a.) y00 4y0 + 4y = xe2x .

(b.) y00 + 5y0 + 6y = sen(x) cos(2x).
Solucin: (a.) Primero resolvemos la ecuacin homognea y00 4y0 + 4y = 0 y obtenemos
como solucin
yh (x) = c1 e2x + c2 xe2x .
Observe que g(x) = xe2x es del tipo 4, es decir y p (x) = xs (Ax + B)e2x donde s = 2 ya
que el el trmino Ax2x es parte de solucin de la homognea por lo tanto,
y p (x) = x2 (Ax + B)e2x .
Para determinar A y B derivamos dos veces y remplazamos en la ecuacin diferencial:
(4Ax3 + (12A + 4B)x2 + (6A + 8B)x + 2B 4(2Ax3 + (3A + 2B)x2 + 2Bx)
+ 4(Ax3 + Bx2 ))e2x = xe2x .
Simplificando tenemos:
(4A 8A + 4A)x3 + (12A + 4B 12A 8B + 4B)x2
+ (6A 8B + 8B)x + 2B = x
de donde concluimos que:
A=
1
y B = 0.
6
Por lo tanto y p (x) = 16 x3 e2x , es decir la solucin general es

1
y = c1 e2x + c2 xe2x + x3 e2x .
6
(b.) 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
y00 + 5y0 + 6y = sin(x)
y00 + 5y0 + 6y = cos(2x).
Para resolver la ecuacin (4.2), g1 (x) = sin(x) es decir
y p1 (x) = xs [A cos(x) + B sin(x)]
(4.2)
(4.3)
106
donde s = 0 ya que y p1 no es una solucin de la homognea, por lo tanto y p (x) =

A cos(x) + B sin(x), para determinar A y B derivamos dos veces y remplazamos en la
ecuacin diferencial:
A cos(x) B sin(x) 5A sin(x) + 5B cos(x) + 6A cos(x) + 6B sin(x) = sin(x).
5A cos(x) + 5B cos(x) 5A sin(x) + 5B sin(x) = sin(x)
y concluimos que:
A=
1
1
yB= .
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
y p2 (x) = xs [C cos(2x) + D sin(2x)]

donde s = 0 ya que el y p2 no es una solucin de la ecuacin homognea, por lo tanto
y p2 (x) = C cos(2x) + D sin(2x),
para determinar C y D derivamos dos veces y reemplazamos en la ecuacin diferencial.
4C cos(2x) 4D sin(2x) 10C sin(2x) + 10D cos(2x) + 6C cos(2x) + 6D sin(2x)
= cos(2x).
2C cos(2x) + 2D sin(2x) 10C sin(2x) + 10D cos(2x) = cos(2x)
y concluimos que
C=
5
1
y D= .
52
52
1
5
Por lo tanto y p2 (x) = 52
cos(2x) 52
sin(2x), es decir la solucin general es
y(x) = C1 e3x +C2 e2x
1
1
1
5
cos(x) + sin(x) cos(2x) sin(2x).
10
10
52
52
Ejercicios
1. Determine 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.
7
c) y00 y0 2y = cos(x) sin(2x); y(0) = 20
, y0 (0) = 15 .
107
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 escribir de la forma an Dn y + an1 Dn1 y + + a1 Dy + a0 y = g(x) en donde Dk =
dk y
, k = 1, 2, , n. Cuando nos convenga, representaremos tambin la ecuacin de
dxk
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:
(D 1)(D + 1)[y] = y00 y = e3x .
(4.4)
a) Muestre que al aplicar D 3 a ambos lados de la ecuacin (4.4) tenemos

(D 3)(D 1)(D + 1)[y] = y00 y = 0
(4.5)
(como (D 3)[e3x ] = 0, decimos que D 3 anula a e3x ).

b) Si y p satisface la ecuacin diferencial (4.4) tambin satisface la ecuacin (4.5).
Verifique que
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).
c) Use el hecho de que c2 ex + c3 ex es una solucin de la ecuacin homognea
(D 1)(D + 1)[y] = 0, y concluya que existe una solucin particular de (4.4) de
la forma y p = c1 e3x .
d) Use el procedimiento anterior para mostrar que existe una solucin particular
(D 1)(D + 1)[y] = ex , de la forma y p = cxe3x .
108
4.4
Reduccin de Orden.
Este mtodo se usa para hallar una segunda solucin de una ecuacin diferencial lineal
no homognea de orden superior cuando se conoce una solucin y1 (x) = g(x) de sta.
Tambin, conocida una solucin de la ecuacin homognea correspondiente, podemos usar
este mtodo para hallar la solucin general de la ecuacin no homognea.
Dada la ecuacin diferencial de orden n:
dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y = f (x)
n1
n
n1
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
an (x)
d2y
dy
(4.7)
+ a1 (x) + a0 (x)y = 0
2
dx
dx
y una solucin no nula y1 (x) = g(x). Tenemos y2 (x) = g(x)u(x) la segunda solucin entonces:
a2 (x)
y0 =g0 u + gu0
y00 =gu00 + 2g0 u0 + g00 u.
Sustituyendo en la ecuacin (4.7):
a2 (x)[gu00 + 2g0 u0 + g00 u] + a1 (x)[g0 u + gu0 ] + a0 (x)[gu] = 0
[a2 (x)g00 + a1 (x)g0 + a0 (x)g]u + a2 (x)gu00 + 2a2 (x)g0 u0 + a1 (x)gu0 = 0.
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:
a2 (x)gw0 + 2a2 (x)g0 w + a1 (x)gw = 0.
Luego, al separar las variables e integrar, tenemos

w0
2a2 (x)g0 a1 (x)
=
+
w
a2 (x)g
a2 (x)
Z
Z
Z
0
dw
2g
a1 (x)
= dx
dx
w
g
a2 (x)
Z
a1 (x)
2
ln |w| = ln(g )
dx
a2 (x)
2
w =g e
R a1 (x)
a2 (x)
dx
4.4 Reduccin de Orden.
109
Reemplazando u0 en lugar de w en la ltima ecuacin e integrando obtenemos

u(x) =
R a1 (x)
a2 (x)
dx
[g(x)]2
dx.
Esta frmula nos da la segunda solucin linealmente independiente la cual es

y2 (x) = g(x)
R a1 (x)
dx
a2 (x)
[g(x)]2
dx.
Ejemplo 4.7. Dada la ecuacin diferencial tx00 (t + 1)x0 + x = 0, t > 0 y una solucin no
trivial x1 (t) = et . Determine una segunda solucin linealmente independiente.
Solucin: En este caso a2 (t) = t y a1 (t) = (t + 1) de modo que

Z
Z
1
(t + 1)
dt =
1+
dt = t + ln(t).
t
t
Donde hemos considerado que la constante de integracin se anula. La frmula para
reduccin de orden implica que
u(t) =
Z
t
x2 (t) =e
et+ln(t)
dt =
(et )2
et+ln(t)
dt = (t + 1)et
2t
e
et+ln(t)
dt = et [et (t + 1)] = (t + 1).
e2t
De manera similar, el cambio y = uy1 nos permite hallar la solucin general de una ecuacin
no homognea
an (x)
d n1 y
dy
dny
+
a
(x)
+ + a1 (x) + a0 (x)y = f (x)
n1
n
n1
dx
dx
dx
si conocemos una solucin y1 de la correspondiente ecuacin homognea.

an (x)
dny
d n1 y
dy
+
a
(x)
+ + a1 (x) + a0 (x)y = 0.
n1
n
n1
dx
dx
dx
Ejemplo 4.8. Hallar la solucin general de la ecuacin

xy00 (2x + 1)y0 + (x + 1)y = x2
dada la funcin y1 (x) = ex una solucin de la correspondiente ecuacin homognea
xy00 (2x + 1)y0 + (x + 1)y = 0.
110
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 .
As, y(x) = u(x)ex es una solucin de la ecuacin no homognea si
(xu00 u0 )ex = x2 .
Si hacemos z = u0 , tenemos que la ecuacin anterior es una ecuacin de primer orden:
(xz0 z)ex = x2 .
Esta ecuacin se escribe como
1
z0 z = xex ,
x
x
la cual tiene como solucin z(x) = xe +C1 x. Regresamos el cambio y tenemos u0 (x) =
xex +C1 x, integrando esto obtenemos u(x) = (x + 1)ex + C21 x2 +C2 .
Finalmente,
C1
y(x) = u(x)ex = x + 1 + x2 ex +C2 ex .
2

Ejemplo 4.9. Hallar la solucin del problema de valor inicial
x2 y00 + xy0 y = x2 + 1,
y(1) = 2, y0 (0) = 3,
dado y1 (x) = x una solucin de la ecuacin homognea correspondiente.

Solucin: Sea y(x) = u(x)y1 (x) = u(x)x. Tenemos que y0 = u0 x + u y y00 = u00 x + 2u0 . Reemplazamos en la ecuacin
x2 y00 + xy0 y =x2 (u00 x + 2u0 ) + x(u0 x + u) ux
=x3 u00 + 3x2 u0 .
Luego, y(x) = u(x)x es una solucin de la ecuacin no homognea si x3 u00 + 3x2 u0 = x2 + 1.

Haciendo el cambio de variable z = u0 tenemos que la ecuacin anterior se escribe como la
ecuacin de primer orden:
3
z0 + z = x1 + x3 .
x
Resolvemos esta ecuacin y obtenemos z(x) = 13 + x2 + c1 x3 . Regresando el cambio
u0 (x) = z(x) e integrando obtenemos
1
1 1
u(x) = x c1 x2 + c2 .
3
x 2
4.4 Reduccin de Orden.
111
Como y = ux, entonces tenemos que

1
1
y(x) = x2 1 c1 x1 + c2 x.
3
2
Resolvamos el problema de valor inicial. y0 (x) = 32 x + 12 c1 x2 + c2 , ahora usemos las
condiciones iniciales:
2 1
2 = y(1) = c1 + c2
3 2
2
1
3 = y0 (1) = + c1 + c2 .
3 2
La solucin de este sistema es c1 = 19/3 y c2 = 1/2, as la solucin del problema de
valor inicial es
1
19
1
y(x) = x2 1 + x1 x.
3
6
2

Ejercicios
1. Se da una ecuacin diferencial y una solucin no trivial f . Determine 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 .
g) y00 + 2y0 + 2y = ex sec(x); f (x) = ex cos(x).
h) x2 y00 + xy0 4y = 6x 4; f (x) = x2 .
i) (1 2x)y00 + 2y0 + (2x 3)y = (1 4x + 4x2 )ex ; f (x) = ex .
j) 4x2 y00 + 4x(x + 1)y0 + (2x + 3)y = 4x5/2 e2x ; f (x) = x1/2 .
k) x2 y00 5xy0 + 8y = 4x2 ; f (x) = x2 .
2
l) y00 + 4xy0 + (4x2 + 2)y = 8ex(x+2) ; f (x) = ex .
2. En fsica matemtica, muchos problemas con simetra esfrica implican el estudio de
la ecuacin de Legendre,
(1 x2 )y00 2xy0 + ( + 1)y = 0;
1 < x < 1,
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.
112
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
transforma la ecuacin diferencial
y00 + p(x)y0 + q(x)y = 0
en una ecuacin de la forma
u00 + f (x)u = 0.
La ltima ecuacin se llama forma normal de una ecuacin lineal de segundo orden
homognea.
5. En mecnica cuntica, el estudio de la ecuacin de Schdinger para el tomo de
hidrgeno conduce a considerar la ecuacin de Laguerre: xy00 + (1 x)y0 + y = 0,
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 Laguerre para = 2, f (x) = x2 4x + 2.
4.5
Variacin de Parmetros
Ahora vamos a estudiar un mtodo, llamado variacin de parmetros, que nos permite hallar
una solucin particular de una ecuacin diferencial lineal no homognea cuyos coeficientes
pueden ser funciones y el trmino no homogneo no tiene una forma especial.
Este mtodo nos proporciona una solucin particular y p de la ecuacin
y(n) + P1 (x)yn1 + + Pn (x)y = f (x),
para cualquier funcin continua f (x).
Dada la solucin general de la ecuacin homognea
y(n) + P1 (x)yn1 + + Pn (x)y = 0
digamos
yh = c1 y1 + c2 y2 + + cn yn .
El mtodo de variacin de parmetros propone una solucin particular de la forma

y p (x) = y1 (x)u1 (x) + y2 (x)u2 (x) + + yn (x)un (x),
donde las funciones u1 , u2 , . . . , un vienen dadas por la frmula
uk (x) :=
f (x)Wk (x)
dx,
W [y1 , . . . , yk ](x)
k = 1, . . . , n
con Wk (x) = (1)nkW [y1 , . . . , yk1 , yk+1 , . . . , yn ](x), k = 1, . . . , n.
4.5 Variacin de Parmetros
113
Ejemplo 4.10. Determinar una solucin general de la ecuacin:

(a) y00 + y = tan2 (x).
(b) y00 + 4y0 + 4y = e2x ln(x).
Solucin: (a) Primero resolvemos la ecuacin homognea y00 + y = 0 y obtenemos como
solucin
yh (x) = c1 cos(x) + c2 sin(x).
Ahora,
y p (x) = u1 (x) cos(x) + u2 (x) sin(x)
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)
= cos2 (x) + sin2 (x) = 1.
W [y1 , y2 ](x) =
sin(x) cos(x)
Por lo tanto,
tan2 (x) sin(x)

tan2 (x) cos(x)
dx, u2 (x) =
dx
1
1
(cos2 (x) + 1)
u1 (x) =
, u2 (x) = ln | sec(x) + tan(x)| sin(x)
cos(x)
u1 (x) =
es decir
(cos2 (x) + 1)
cos(x) + (ln | sec(x) + tan(x)| sin(x)) sin(x)
cos(x)
y p (x) = sin(x)(ln | sec(x) + tan(x)|) 2
y p (x) =
as, la solucin general de la ecuacin es:

y(x) = c1 cos(x) + c2 sin(x) + (ln | sec(x) + tan(x)|) sin(x) 2.
(b) Resolvemos la ecuacin homognea y00 + 4y0 + 4y = 0 y obtenemos como solucin
yh (x) = c1 e2x + c2 xe2x .
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

2x
e

xe
= e4x 2xe4x + 2xe4x = e4x .
W [y1 , y2 ](x) =
2x
2x
2x
2e
e 2xe
114
Por lo tanto,
Z
u1 (x) =
e2x ln(x)xe2x
e2x ln(x)e2x
dx,
u
(x)
=
dx
2
4x
e4x
Ze
Z
u1 (x) =
u1 (x) =
x ln(x)dx,
u2 (x) =
ln(x)dx
2x2 ln(x) + x2
, u2 (x) = x ln |x| x
4
es decir:
2x2 ln(x) + x2 2x
e + (x ln |x| x)xe2x
4
e2x x2 (2 ln |x| 3)
y p (x) =
4
y p (x) =
luego, la solucin general de la ecuacin es:

y(x) = c1 e2x + c2 xe2x +
e2x x2 (2 ln |x| 3)
.
4

y000 2y00 y0 + 2y = e3x ;
y(0) = 1, y0 (0) = 1, y00 (0) = 0.
Solucin: Primero resolvamos la ecuacin homognea

y000 2y00 y0 + 2y = 0
cuya solucin es yh (x) = c1 e2x + c2 ex + c3 ex . As, y1 (x) = e2x , y2 (x) = ex y y3 (x) = ex .
Calculemos Wk (x), k = 1, 2, 3.
x

e
ex
31

W1 (x) = (1) W [y2 , y3 ](x) = x x = 2
e
e
2x

x
e
e
W [y1 , y3 ](x) = 2x x = e3x
2e
e
32
W2 (x) = (1)
2x

e
ex

W [y1 , y2 ](x) = 2x
= 3ex .
2e
ex
33
W3 (x) = (1)
El Wronskiano de {y1 , y2 , y3 } es
2x

x
x
e
e
e

W [y1 , y2 , y3 ](x) = 2e2x ex ex = 6e2x .
4e2x ex ex
4.5 Variacin de Parmetros
115
Ahora calculamos
Z
2e3x
1
dx = ex ,
2x
6e
3
Z 3x 3x
e e
1
u2 (x) =
dx = e4x ,
2x
6e
24
Z
x
3x
3e e
1
u3 (x) =
dx = e2x .
2x
6e
4
u1 (x) =
De donde la solucin general de la ecuacin no homognea es

1
1
1
yg (x) =c1 e2x + c2 ex + c3 ex + e3x + e3x e3x
3
24
4
1
=c1 e2x + c2 ex + c3 ex + e3x .
8
Finalmente, calculemos la solucin del problema de valor inicial. Usando las condiciones
iniciales tenemos el sistema:
1 = yg (0) =c1 + c2 + c3 +
1
8
3
8
9
0 = y00g (0) =4c1 + c2 + c3 + .
8
1 = y0g (0) =2c1 c2 + c3 +
3
cuya solucin es c1 = 32 , c2 = 13
24 y c3 = 4 . As, la solucin del problema de valor inicial
es
2
13
3
1
y(x) = e2x + ex + ex + e3x .
3
24
4
8
Ejemplo 4.12. Hallar una solucin particular y p de la ecuacin

(x 1)y00 xy0 + y = (x 1)2 ,
dado que y1 (x) = x y y2 (x) = ex son soluciones de la correspondiente ecuacin homognea
(x 1)y00 xy0 + y = 0.
Solucin: La solucin particular y p tiene la forma
y p (x) = u1 (x)x + u2 (x)ex .
Para calcular u1 y u2 primero escribamos la ecuacin en forma la cannica
y00
x 0
1
y+
y = x 1,
x1
x1
116
ahora tenemos que

x ex
= (x 1)ex ,
W [y1 , y2 ](x) =
1 ex
W1 (x) = (1)21W [y2 ](x) = ex , W2 (x) = (1)22W [y1 ](x) = x, de donde tenemos
Z
(x 1)ex
dx = x,
ex (x 1)
Z
(x 1)x
dx = ex (x + 1).
u2 (x) =
ex (x 1)
u1 (x) =
As, la solucin particular de la ecuacin es

y p (x) = x2 x 1.

Ejercicios
1. Resuelva 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 .
000
00
0
h) y y + y y = ex sin(x).
i) y(4) + 2y00 + y = sin(t).
2. Hallar la solucin del problema de valor inicial dado.
a) y000 + y0 = sec(x); y(0) = 2, y0 (0) = 1, y00 (0) = 2.
b) y(4) + 2y00 + y = sin(x); y(0) = 2, y0 (0) = 0, y00 (0) = 1, y000 (0) = 1.
c) y000 y0 = csc(x); y(/2) = 2, y0 (/2) = 1, y00 (/2) = 1.
3. Hallar una solucin particular de la ecuacin, dadas las soluciones de la correspondiente ecuacin homognea.
a) 4x2 y00 + (4x 8x2 )y0 + (4x2 4x 1)y = 4x1/2 ex , x > 0; y1 (x) = x1/2 ex , y2 (x) =
x1/2 ex .
b) x2 y00 4xy0 + 6y = x5/2
, x > 0; y1 (x) = x2 , y2 (x)
= x3 .
c) 4xy00 + 2y0 + y = sin( x), x > 0; y1 (x) = cos( x), y2 (x) = sin( x).
d) sin(x)y00 + (2 sin(x) cos(x))y0 + (sin(x) cos(x))y = ex ; y1 (x) = ex ,
y2 (x) = ex cos(x).
e) x2 y00 x(x + 4)y0 + 2(x + 3)y = x4 ex ; y1 (x) = x2 , y2 (x) = x2 ex .
f ) (x 1)y00 xy0 + y = 2(x 1)2 ex ; y1 (x) = x, y2 (x) = ex .
g) 4x2 y00 4x(x + 1)y0 + (2x + 3)y = x5/2 ex , x > 0; y1 (x) = x, y2 (x) = xex .
4.6 Ecuacin de Cauchy-Euler
117
4. Use el mtodo de variacin de parmetros para demostrar que

y(x) = c1 cos(x) + c2 sin(x) +
Z b
a
f (s) sin(x s)ds
es una solucin general de la ecuacin diferencial y00 + y = f (x), donde f (x) es

una funcin continua en (, ). (Sugerencia: Use la identidad trigonomtrica
sin(x s) = sin(x) cos(s) sin(s) cos(x) ).
5. Determinar a para que y(x) = ax3 sea solucin de la ecuacin xy00 y0 = 3x2 y usar
el resultado para escribir una solucin general.
4.6
Ecuacin de Cauchy-Euler
Una ecuacin diferencial lineal de la forma:
an xn
n1 y
dny
dy
n1 d
+
a
x
+ + a1 x + a0 y = g(x),
n1
n
n1
dx
dx
dx
con an , an1 , . . . , a0 constantes es llamada una ecuacin de Cauchy-Euler. La caracterstica

de este tipo de ecuacin es que el grado de xk , k = 0, 1, . . . , n coincide con el orden k de la
dk y
derivada dx
k.
Primero resolvamos la ecuacin Cauchy-Euler homognea:
an xn
n1 y
dny
dy
n1 d
+
a
x
+ + a1 x + a0 y = 0.
n1
n
n1
dx
dx
dx
(4.8)
Para ello supongamos que una solucin es de la forma y(x) = xm . El problema se reduce a
encontrar los valores del exponente m que hace que las funciones y(x) = xm sean solucin
de la ecuacin. As,
y(x) = xm implica que y0 (x) = mxm1 , y00 (x) = m(m 1)xm2 , . . .
sustituyendo en la ecuacin tenemos
an xn xmn m(m 1)(m 2) [m (n 1)]+ + a1 xmxm1 + a0 xm = 0
an xm m(m 1)(m 2) [m (n 1)]+ + a1 mxm + a0 xm = 0
xm [an m(m 1)(m 2) [m (n 1)]+ + a1 m + a0 ] = 0.
Luego y(x) = xm es una solucin de la ecuacin diferencial siempre que m sea una solucin
de la ecuacin (caracterstica):
an m(m 1)(m 2) [m (n 1)] + + a1 m + a0 = 0.
De lo cual podemos tener los siguientes casos:
118
Caso I. El polinomio caracterstico tiene n races reales simples y distintas: m1 , m2 , . . . , mn .

En este caso tenemos n soluciones de la ecuacin que sern de la forma:
y1 (x) = xm1 , y2 (x) = xm2 , . . . , yn (x) = xmn .
Puede verificarse que el Wronskiano es no nulo para todo x R , por lo tanto por la
Definicin 4.3 estas funciones forman un conjunto fundamental de soluciones y la
solucin general de la ecuacin es:
y(x) = C1 xm1 +C2 xm2 + +Cn xmn .
Caso II. Las races del polinomio son reales pero algunas con multiplicidad mayor que uno.
Supongamos, para ejemplificar, que tenemos la ecuacin de segundo orden ax2 y00 +
bxy0 + cy = 0, sustituyendo tenemos:
ax2 m(m 1)xm2 + bxmxm1 + cxm =0
am(m 1)xm + bmxm + cxm =0
xm [am(m 1) + bm + c] =0
am2 + m(b a) + c =0
con dos races iguales, digamos, m1 . Entonces y1 (x) = xm1 es una solucin es una
solucin de la ecuacin diferencial. Para encontrar la otra solucin utilizamos el
mtodo de reduccin de orden (ver Pgina 109)
y2 (x) = y1 (x)
b
en este caso p(x) = ax
y m1 =
ecuacin obtenemos:
y2 (x) =x
m1
y2 (x) =x
m1
Z
Z
e
x
(ba)
2a
e p(x)dx
dx,
[y1 (x)]2
ya que (b a)2 4ac = 0 sustituyendo en la
R b
ax dx
(ba)
2a
e a ln(x)
b
x a +1
2 dx
dx
resolviendo la integral obtenemos que y2 (x) = ln(x)xm1 . Estas dos soluciones son
linealmente independientes y la solucin general de la ecuacin es:
y(x) = C1 xm1 +C2 xm1 ln(x),
si para la ecuacin (4.8) la multiplicidad de la raz en la ecuacin caracterstica es n,
se extiende estas soluciones obtenindose n soluciones linealmente independientes de
la forma ln2 (x)xm1 , ln3 (x)xm1 , . . . , lnn1 (x)xm1 . De esta forma, la solucin general es:
y(x) = C1 xm1 +C2 xm1 ln(x) +C3 xm1 ln2 (x) + +Cn lnn1 (x)xm1 .
119
Caso III. El polinomio tiene un par de races complejas conjugadas :

m1 =a + ib
m2 =a ib.
Como y(x) = xm es una solucin de la ecuacin diferencial y m1 = a + ib, m2 = a ib
son soluciones de la ecuacin caracterstica se tiene que, y1 (x) = xa+ib , y2 (x) = xaib ,
son soluciones particulares de la ecuacin diferencial.
Puede tomarse la expresin de estas exponenciales segn la frmula de Euler:
xa+ib = xa xib = xa eib ln(x) = xa [cos(b ln(x)) + i sin(b ln(x))].
De esta forma podemos escribir las soluciones de la siguiente manera:
C1 xa [cos(b ln(x)) + i sin(b ln(x))] +C2 xa [cos(b ln(x)) + i sin(b ln(x))] =
= C1 xa [cos(b ln(x)) + i sin(b ln(x))] +C2 xa [cos(b ln(x)) i sin(b ln(x))]
= xa [(C1 +C2 ) cos(b ln(x)) + (iC1 iC2 ) sin(b ln(x))]
= xa [K1 cos(b ln(x)) + K2 sin(b ln(x))].
Si la multiplicidad de la raz compleja es r, la solucin general viene dada por:
y(x) = xa [{K1 cos(b ln(x)) + K2 sin(b ln(x))} + ln(x){K3 cos(b ln(x))+
K4 sin(b ln(x))} + + (ln(x))r1 {K2r1 cos(b ln(x)) + K2r sin(b ln(x))}].
Ejemplo 4.13. Determine la solucin general de la ecuacin diferencial:
x3 y000 6y = 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.
Factorizando obtenemos (m 3)(m2 + 2) = 0. Por lo tanto, las races de la ecuacin auxiliar

son:
m = 3, m = i 2.
Una raz real y una imaginaria pura, luego a = 0, b = 2, y entonces una solucin general
de la ecuacin es
y(x) = c1 x3 + c2 cos( 2 ln(x)) + c3 sin( 2 ln(x)).

120
Ahora consideremos la ecuacin de Cauchy-Euler no homognea

an xn
n1 y
dy
dny
n1 d
+
a
x
+ + a1 x + a0 y = g(x),
n1
n
n1
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.
Ejemplo 4.14. Encontrar la solucin general de la ecuacin
x2 y00 + 3xy0 + y = x1 .
Solucin: Se trata de una ecuacin de Cauchy-Euler. Haciendo x = et , suponiendo x > 0,
tenemos
t = ln(x)
dt 1
=
dx x
dy 1
dy dy dt
=
=
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 =
dx dt
d 2 y d 2 y dy
x2 2 = 2 .
dx
dt
dt
(4.9)
(4.10)
Sustituyendo (4.9) y (4.10) en la ecuacin x2 y00 + 3xy0 + y = x1 obtenemos:

d 2 y dy
dy
+ 3 + y(t) = et
2
dt
dt
dt
2
d y
dy
+ 2 + y(t) = et ,
2
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
y obtenemos como solucin
yh (t) = c1 et + c2tet .
Usaremos variacin de parmetros para hallar una solucin particular de la ecuacin no

homognea:
y p (t) = u1 (t)et + u2 (t)tet
121
con y1 (t) = et y y2 (t) = tet . Hallamos el Wronskiano,

t
t

e
te
= e2t te2t + te2t = e2t
W [y1 , y2 ](t) = t t
t
e
e te
por lo tanto
u1 (t) =
et tet
et et
dt,
u
(t)
=
dt
2
2t
e2tZ
Z e
t, u2 (t) =
u1 (t) =
u1 (t) =
dt
t 2
, u2 (t) = t
2
es decir:
t 2 t 2 t
e +t e
2
t 2 et
y p (t) =
.
2
As, la solucin general de la ecuacin con coeficientes constantes es:
y p (t) =
t 2 et
2
t
sustituimos, x = e y se obtiene la solucin general de la ecuacin de Cauchy-Euler:
y(t) = c1 et + c2tet +
y(x) = c1 x1 + c2 ln(x)x1 +
ln2 (x)x1
.
2
Ejemplo 4.15. Hallar la solucin del problema de valor inicial

ln(x)
, y(1) = 0, y0 (1) = 1, y00 (1) = 1.
x
Solucin: Esta ecuacin es del tipo Cauchy-Euler ya que esta se puede reescribir como
x2 y000 xy00 + y0 =
x3 y000 x2 y00 + xy0 = ln(x).

Haciendo el cambio de variable x = et tenemos por la regla de la cadena
dy
dy dt
= dt =
dx dx
dy
dt
= et
d t
dt (e )
2
et ddt 2y
dy
dt
2

dy0
et dy
d 2 y dt
dy
dt
2t d y
=
=
=e
dx2 dx
et
dt 2 dt
dt

3

2t d 3 y d 2 y 2e2t d 2 y dy
3
e
dt
d2y
dy
d y
dt 3
dt 2
dt 2
3t d y
=
=e
3 2 +2
.
dx3
et
dt 3
dt
dt
(4.11)
122
Reemplazando esto en la ecuacin (4.11) tenemos

d3y
d2y
dy d 2 y dy dy
3
+
2
+ +
=t
dt 3
dt 2
dt dt 2 dt dt
d3y
d2y
dy
4 2 + 4 =t,
3
dt
dt
dt
(4.12)
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
m3 4m2 + 4m = m(m 2)2 = 0,
as, el conjunto fundamental de soluciones es {y1 (t) = 1, y2 (t) = e2t , y3 (t) = te2t }. Hallaremos una solucin particular de la ecuacin no homognea usando variacin de parmetros.
Calculemos el Wronskiano

1 e2t

te2t

2t
2t

W [y1 , y2 , y3 ](t) = 0 2e (2t + 1)e = 4e4t .
0 4e2t (4t + 4)e2t
Ahora hallemos ui , i = 1, 2, 3 tales que y p = u1 y1 + u2 y2 + u3 y3 . Para ello calculemos

2t

2t
e

te
31
= e4t
W1 (t) =(1) W [y2 , y3 ](t) = 2t 2t
2t
2e e + 2te

1
te2t
32

W2 (t) =(1) W [y1 , y3 ](t) =
= e2t (2t + 1)
0 e2t + 2te2t

1 e2t
33

= 2e2t ,
W3 (t) =(1) W [y1 , y2 ](t) =
0 2e2t
con esto tenemos que
Z
1
te4t
u1 (t) =
dt = t 2 ,
4t
4e
8
Z
1 2t
te2t (2t + 1)
u2 (t) =
dt
=
e (3 + 6t + 4t 2 ),
4e4t
16
Z
2te2t
1
dt = e2t (2t + 1).
u3 (t) =
4t
4e
8
As, la solucin general de la ecuacin (4.12) es
1
1
1
y(t) =c1 + c2 e2t + c3te2t + t 2 + e2t (3 + 6t + 4t 2 )e2t e2t (2t + 1)te2t
8
16
8
1
1
=c1 + c2 e2t + c3te2t + t 2 + t.
8
4
123
Regresamos el cambio de variable x = et y obtenemos t = ln(x), luego la solucin general

de la ecuacin (4.11) es
1
1
y(x) = c1 + c2 x2 + c3 x2 ln(x) + ln2 (x) + ln(x).
8
4
Ahora resolvamos el problema de valor inicial.
1 ln(x) 1
+
4 x
4x

1
1 1 ln(x)
00
2,
y (x) =2c2 + c3 (3 + 2 ln(x)) +
2
4
x
4x
y0 (x) =2c2 x + c3 (x + 2x ln(x)) +
usando las condiciones iniciales obtenemos el sistema de ecuaciones

0 = y(1) =c1 + c2
1 = y0 (1) =2c2 + c3 +
1 = y00 (1) =2c2 + 3c3
1
4
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 CauchyEuler 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. Determine una solucin general de la ecuacin diferencial dada
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.
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
3
000
2
00
k) x y x y + 2xy0 2y = x3 .
l) x3 y000 + 5x2 y00 + 2xy0 2y = x4 .
124
m) x3 y000 4x2 y2 + 8xy0 9y = 4 ln(x).

n) x3 y000 + x2 y0 6xy0 + 6y = 30x.
3. Determine una solucin general de la ecuacin diferencial dada
a) x2 y00 + 73xy0 + y = x1 .
b) 21 y00 + 2y = tan(2x) 21 ex .
c) (1 + x2 )y00 + xy0 y + 1 = 0 sabiendo que y1 (x) = x es una solucin particular.
d) x3 y000 4x2 y00 + 8xy0 + 8y = 4 ln(x).
4. La ecuacin de Bessel de orden un medio,
5
1
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). Determine una solucin general de la ecuacin no homognea.
4.7
Aplicaciones a la Fsica
En esta seccin estamos estudiaremos aplicaciones modeladas mediante ecuaciones diferenciales lineales de orden 2 con coeficientes constantes de la forma
d2x
dx
+ 2 x = f (t),
+
2r
dt 2
dt
con f una funcin continua en un intervalo I, r y constantes positivas.
4.7.1
Movimiento Armnico Simple

Muchos objetos tienen un movimiento vibratorio natural oscilando hacia adelante y atrs
sobre un punto de equilibrio fijo. Una partcula que oscila de esta manera en un medio con
resistencia, o factor de amortiguamiento, negativo se dice que ejecuta un movimiento sin
amortiguamiento o movimiento armnico simple. Ejemplos de objetos con movimiento
como el descrito son los resortes helicoidales y los pndulos.
As, formalmente, una partcula se dice que ejecuta un movimiento armnico simple si la
ecuacin del movimiento satisface una ecuacin diferencial de la forma
d2x
+ 2x = 0
2
dt
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
x(t) = c1 cos(t) + c2 sin(t).
Tomamos A 0, haciendo c1 = A sin( ) y c2 = A cos( ), sustituyendo en la ecuacin
anterior obtenemos:
c1 cos(t) + c2 sin(t) =A sin( ) cos(t) + A cos( ) sin()
=A sin(t + ).
4.7 Aplicaciones a la Fsica
125
q
Al hallar A y en trminos de c1 y c2 , tenemos A = c21 + c22 , tan( ) =
solucin de la ecuacin del movimiento armnico simple queda como:
c1
c2 ;
es decir, la
x(t) = A sin( + t).

La constante A es llamada la amplitud del movimiento y es el ngulo de fase del moviw
miento oscilatorio. El movimiento es peridico con periodo 2
w y frecuencia natural 2 .
El centro del segmento de recta en el cual la partcula se mueve de ida y vuelta se llama
posicin de equilibrio. El periodo se mide en unidades de tiempo, la frecuencia natural
ciclos/segundo y la frecuencia angular radianes/segundo.
Obs Note que si reemplazamos por en el procedimiento anterior tenemos que la
solucin de la ecuacin diferencial que modela el movimiento de una partcula que
ejecuta un movimiento armnico simple queda como
x(t) = A cos( + t).
As, la solucin de la ecuacin del movimiento armnico simple puede ser escrita en
cualquiera de las siguientes formas:
x(t) =c1 cos(t) + c2 sin(t),
x(t) =A sin( + t),
x(t) =A cos( + t).
Figura 4.1: Movimiento armnico simple de vibraciones libres no amortiguado.

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:
126
(a.)
(b.)
(c.)
(d.)
(e.)
La ecuacin del movimiento de la partcula.

La amplitud del movimiento.
El ngulo fase.
El periodo del movimiento.
La frecuencia del movimiento en ciclos por segundo.
Solucin: (a.) Como la partcula ejecuta un movimiento armnico simple, su ecuacin

de movimiento viene dada por
x(t) = A cos( + t),
v(t) =
dx
= A sin( + t).
dt
La frecuencia es 4 rad/seg, es decir = 4 y la ecuacin de movimiento queda como

x(t) = A cos( + 4t),
v(t) =
dx
= A sin( + 4t).
dt
(4.13)
Las condiciones iniciales son x(0) = 0 y v(0) = 4. Reemplazando estos valores en

(4.13) obtenemos
0 = A cos( ),
(b.)
(c.)
(d.)
(e.)
4 = 4A sin( ).
como la amplitud A 6= 0 tenemos que 0 = A cos( ) si = 2 . Evaluamos este

valor en la ecuacin 4 = 4A sin( ) y concluimos que A = 1. As, la ecuacin de
movimiento es:

x(t) = cos + 4t
o, x(t) = cos
+ 4t .
2
2
La amplitud del movimiento es 1.
El ngulo fase es 2 .
El periodo del movimiento es 2 radianes.
La frecuencia del movimiento en ciclos por segundo es 2 .
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, muestre que la partcula ejecuta un movimiento armnico simple. Describa el
movimiento.
Solucin: Por hiptesis tenemos que
F = kx
donde k > 0 es la constante de proporcionalidad. El signo negativo es necesario pues cuando
la partcula en el punto P1 (ver Figura 4.2), x es positiva y F acta en una direccin negativa;
cuando la partcula est en P2 , x es negativa y F acta en una direccin positiva. F y x
siempre tienen signos opuestos.
d2 x
m dt2 = kx,
4.7 Aplicaciones aFla= Fsica
d2 x
dt2
k
= m
x.
127
+
P2
P1
Figura 4.2: Descripcin del movimiento.

Luego, de la ley de Newton tenemos
d2x
= kx,
dt 2
d2x
k
= x.
2
dt
m
F =m
Como k y m son constantes positivas entonces podemos escribir 2 = k/m, as la ecuacin

diferencial queda como
d2x
+ 2x = 0
dt 2
(4.14)
la cual es la ecuacin del movimiento armnico simple. Ahora describamos el movimiento.

La solucin de la ecuacin (4.14) es
x(t) = A sin(t + ).
Derivando la ecuacin anterior tenemos
dx
= v(t) = A cos(t + ),
dt
donde v es la velocidad de la partcula. Como | sin( )| 1, entonces |x(t)| A. As la
partcula nunca pasa de los puntos A y A; por lo tanto, estos puntos representan los
desplazamientos mximos de la partcula desde el origen. Cuando |x(t)| = A, se tiene que
| sin(t + )| = 1 lo cual implica que cos(t + ) = 0, por lo tanto cuando |x(t)| = |A| la
velocidad v es cero. Esto prueba que la velocidad de la partcula en los puntos finales A es
cero.
cuando x(t) = 0; es decir, cuando la partcula est en el origen, se tiene que sin(t + ) = 0,
de donde se sigue que | cos(t + )| = 1. Esto significa que la partcula alcanza su velocidad
mxima vmax = |A| en el origen. Tambin se verifica que para x(t) entre 0 y |A| la velocidad
de la partcula est entre 0 y su valor mximo |A|, esta velocidad crece cuando la partcula
parte de A, donde la velocidad es cero, hasta el origen donde alcanza su mximo, mientras
que decrece cuando la partcula cruza el origen hasta que sea de nuevo cero en la posicin
A. De esta forma, la partcula oscila de ida y vuelta sin parar en un ciclo A, desde A
hasta A.
128
El movimiento de una masa sujeta a un resorte, que es un ejemplo de las vibraciones que
ocurren en sistemas mecnicos, se describe de la siguiente manera:
Consideremos un resorte helicoidal flexible con una longitud natural ` suspendido de un
soporte fijo. A l se le sujeta un objeto de masa m y luego entra en equilibrio. Debido al
estiramiento producto del objeto sujeto, se crea una fuerza o tensin F en el resorte la cual
trata de regresar el resorte a su longitud natural. La ley de Hooke dice que el resorte mismo
ejerce una fuerza F de restitucin proporcional a la distancia y que el resorte ha sido estirado
opuesta a la direccin del alargamiento que sucede con la masa m. As,
F = ky
donde k > 0 es la constante de proporcionalidad, llamada constante del resorte (o rigidez).
Como el resorte entra en equilibrio, la fuerza de restitucin debe ser igual a la fuerza de
peso (haca abajo) ejercida por la masa m. Esto es,
ky = mg.
Sea y = 0 la posicin de equilibrio del resorte con el objeto de masa m sujeto a l. Si el
resorte se estira una distancia adicional y, entonces las siguientes fuerzas actuarn en el
resorte
1.) Una fuerza de tensin del resorte (haca arriba) la cual por la ley de Hooke es k(` + y).
2.) Una fuerza debido a la fuerza de peso de la masa m sujeta al resorte la cual es igual a
mg.
longitud del resorte
y = 0 posicion de equilibrio
y
m
Figura 4.3: Ley de Hooke

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
m
d2y
=mg k(` + y)
dt 2
= ky.
129
o, equivalentemente,
d2y k
+ y = 0.
dt 2 m
(4.15)
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
t + c2 sin
t .
y(t) =c1 cos
m
m
r

k
=A sin
t + .
m
q
c21 + c22 , ngulo fase dado por tan( ) =
k/m
miento es peridico con periodo 2 y frecuencia natural 2 .
con amplitud de movimiento A =
(4.16)
(4.17)
c1
c2 .
El movi-
k/m
Ejemplo 4.18. Una masa de 3 Kg est unida a un resorte con rigidez k = 48 Nn . La masa se
desplaza 12 mts a la izquierda del punto de equilibrio y recibe una velocidad de 2 mts/seg
hacia la derecha. La fuerza de amortiguamiento es despreciable. Determine la ecuacin
de movimiento de la masa junto con su amplitud, periodo y frecuencia. Cunto tiempo
despus de su liberacin pasa la masa 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
,
w=
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
Como A = c21 + c22 , tenemos que A =
21 + 12 , luego la amplitud es A = 22 mts.
130
Para hallar el ngulo fase tomamos tan( ) = cc21 , reemplazando tenemos tan( ) = 1/2
1/2 = 1,
de donde = arctan(1) = 4 radianes que es el ngulo fase. Periodo
frecuencia es
w
4
2 ciclos
=
=
.
2
2
seg
2
4
seg, la
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
sin(4t + )
y(t) = cos(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
, donde n es
4
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

Ejemplo 4.19. Un cuerpo sujeto a un resorte helicoidal ejecuta un movimiento armnico
simple. La frecuencia natural del movimiento es 2 ciclos por segundo y su amplitud es de 1
mts. Hallar la velocidad del cuerpo cuando este pasa por el punto y = 12 mts.
Solucin: Como el cuerpo ejecuta un movimiento armnico simple, la ecuacin del movimiento es
!
r
k
t + .
y(t) = A cos
m
La amplitud
q es 1 mts; esto es, A = 1. La frecuencia es 2 ciclos por segundo de donde tenemos
q
k
que 2 = 2m , o mk = 4. Sustituyendo estos valores en la ecuacin del movimiento
tenemos
dy
y(t) = cos( + 4t), v(t) =
= 4 sin( + 4t).
dt
Cuando y = 12 obtenemos de la primera ecuacin
1
= cos( + 4t),
2
Luego, cuando y = 21 , sin( + 4t) =
3
2 .
+ 4t = arc cos(1/2).
Finalmente, obtenemos que
3
v(t) =
,
2
131
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 equilibrio.
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
por la posicin de equilibrio con una velocidad de 8 mts/seg.
(a.) Hallar su ecuacin de movimiento.
(b.) Hallar el periodo, la frecuencia y amplitud del movimiento.
4. Una partcula ejecuta un movimiento armnico simple con periodo a seg. Su velocidad en t = 0 cuando pasa por el punto x = x1 es v1 . Hallar su ecuacin de
movimiento.
5. Una partcula ejecuta un movimiento armnico simple, su frecuencia es 3 ciclos por
segundo. En t = 1 seg la partcula est a 3 mts de la posicin de equilibrio a una
velocidad de 6 mts/seg. Hallar la ecuacin del movimiento.
6. Un objeto que pesa 8 Kg que se mueve en linea recta es atrado haca el origen por
una fuerza F que es proporcional a la distancia del objeto al origen. Si esta fuerza es
de 6 Kg a una distancia de a una distancia de -2 mts, hallar la frecuencia natural del
sistema.
7. Un objeto de masa m que se mueve en linea recta es repelido desde el origen por
una fuerza F. Si la fuerza es proporcional a la distancia del objeto al origen, hallar
la posicin del objeto en funcin del tiempo. (Note que el objeto no ejecuta un
movimiento armnico simple.)
8. Un cuerpo de 12 Kg se engancha a un resorte helicoidal y lo estira 6 cm, luego
este regresa a su posicin inicial, se estira 4 cm y se suelta. Hallar la ecuacin del
movimiento, tambin el periodo, frecuencia y amplitud del movimiento. Con qu
velocidad el cuerpo cruza la posicin de equilibrio?.
9. Un cuerpo de 10 Kg estira un resorte 3 cm, luego de regresar al reposo se estira 6 cm
y comienza a oscilar.
132
(a.) Hallar la posicin y velocidad al final de 1 seg.

(b.) Cundo y con qu velocidad el cuerpo pasa por primera vez por la posicin de
equilibrio?.
(c.) Cual es la velocidad del cuerpo cuando est a 3 cm del equilibrio?.
(d.) Cundo su velocidad ser de 2 mts/seg y en qu posicin estar en ese instante?.
(e.) Cual es el periodo, frecuencia y amplitud del movimiento?.
Solucin
1. (a.)
(b.)
(c.)
(d.)
(e.)
2. (a.)
(b.)
(c.)
3. (a.)
(b.)

x(t) = 4 2 sin 4 + t .
4 2.
1 rad/seg.
4 .
4 segundos.
x(t) = 5 cos (8t); v(t) = 40 sin(8t).
8 rad/seg.
1
16 seg, 40 mts/seg.

x(t) = 6 sin 43 t .
3
seg, 34 rad/seg o 32 ciclos por segundo, 6 mts.
2q

4. x(t) = x12 + a2 v21 /4 sin + 2ta , donde = 2 x1 2 2 .
x1 +a v1 /4

p
5. x(t) = 9 + 1/ 2 sin ( + 6t), donde = arcsin 32
.
9 +1
o 3 ciclos
6. 12 rad/seg
por segundo.
q
q
k
7. x(t) = c1 e m + c2 e m , donde k es una constante de proporcionalidad.

8. y(t) = 13 cos(8t), 4 seg, 8 rad/seg o 4 ciclos por segundo, 13 mts, 83 mts/seg.
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.
4.7.2
Movimiento Armnico Amortiguado

Se dice que una partcula ejecuta un movimiento armnico amortiguado si su ecuacin de
movimiento satisface una ecuacin diferencial de la forma
d2x
dx
dx
d2x
m 2 + 2mr + m 2 x = 0,
+
2r
+ 2x = 0
(4.18)
2
dt
dt
dt
dt
donde el coeficiente 2mr > 0 es llamado el coeficiente de resistencia del sistema. Como
antes es la frecuencia (no amortiguada) del sistema y m la masa de la partcula.
2
2
La
ecuacin auxiliar de (4.18) es z + 2rz + = 0 la cual tiene como races z1,2 = r
r2 2 , estas muestran tres casos posibles, en las que las races son diferentes, iguales y
complejas.
133
Caso I: r2 2 > 0. Se tiene races reales distintas con 6= r y r > .
r2 2 t
x(t) = c1 ez1t + c2 ez2t = c1 er+
r2 2 t
+ c2 er
(4.19)
Como ambos exponentes son cantidades negativas, podemos escribir (4.19) como
x(t) = c1 et + c2 et ,
, < 0.
Si c1 , c2 6= 0 tal que c1 y c2 con igual signo entonces dado que et > 0 para todo t, no existe
t tal que x(t ) = 0 (es decir, el grfico de la funcin x no cruza el origen, ver Figura 4.4).
Figura 4.4: c1 c2 > 0.

Por otro lado, si c1 , c2 6= 0 adems c1 y c2 tienen diferente signo entonces existe t tal
que x(t ) = 0 (ver Figura 4.5). El punto t ser:
0 = x(t ) = c1 et + c2 et
c2
e( )t =
c1

c2
1
ln .
t =
c1
Figura 4.5: c1 c2 < 0.
134
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 ,
dt
, < 0.
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:
x(t) =c1 ert + c2tert
x0 (t) = c1 rert + c2 ert c2 rtert .
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.
Caso
III: r2 2 < 0. En este caso las races son imaginarias conjugadas z1,2 = r
2 r2 i. La ecuacin resultante del movimiento es:

p
p

rt
x(t) = e
c1 cos 2 r2t + c2 sin 2 r2t ,
la cual podemos escribir como
p
Aert sin( 2 r2t + ),
con A = c21 + c22 . Debido al trmino sin( 2 r2t + ) en la solucin, el movimiento es

oscilatorio. La amplitud (amortiguada) del movimiento es Aert y como r > 0 este factor
decrece cuando t crece y tiende a cero cuando t . As en el tiempo la partcula vibra
con oscilacin cada vez menor alrededor de la posicin de equilibrio. Las funciones que
aparecen en la solucin no son peridicas ya que sus valores no se repiten, sin embargo dado
que el movimiento es oscilatorio, decimos que la funcin es amortiguada peridicamente y
definimos su periodo (amortiguado) como el tiempo que toma la partcula, comenzando desde la posicin de equilibrio, en hacer una oscilacin completa. As, su periodo amortiguado
viene dado por
2
T=
.
2 + r2
135
Lafrecuencia (amortiguada) del movimiento es 2 r2 radianes por unidad de tiempo,

2
2
o 2r ciclos por unidad de tiempo. El factor ert se llama factor de amortiguamiento.
Como este factor decrece con el tiempo, el movimiento eventualmente acaba. Cuando t = 1r
el factor de amortiguamiento es 1e . El tiempo que toma el factor de amortiguamiento en
alcanzar este valor 1/e es llamado constante de tiempo. As la constante de tiempo es
= 1/r.
Decimos que el movimiento est subamortiguado porque la fuerza de amortiguamiento r es
pequea comparada con la fuerza de restitucin .
Figura 4.6: Movimiento oscilatorio cuya amplitud decrece con el tiempo.

El desplazamiento de una masa sujeta a un resorte como se presenta en la parte inicial
del movimiento armnico simple se amortigua por fuerzas externas como por fuerzas
internas. Por lo menos habr la fuerza de amortiguacin debida al medio que lo rodea. Se
ha determinado que las fuerzas de amortiguamiento son mltiplos de dx
dt , es decir esta
fuerza es proporcional a la velocidad instantnea de la masa. Llamamos b a la constante de
proporcionalidad. La segunda ley de Newton dice que:
d2x
dx
m 2 + b + kx = 0.
dt
dt
Dividiendo por m, 2 =
b
m
y haciendo w2 =
k
m
tenemos que
d2x
dx
+ 2 + w2 x = 0,
2
dt
dt
Ejemplo 4.20. El movimiento de un sistema masa-resorte con amortiguamiento est descrito por
y00 (t) + by0 (t) + 64y(t) = 0, y(0) = 1, y0 (0) = 0.
136
Determine la ecuacin de movimiento y bosqueje la grfica para b = 0, 10, 16 y 20.
Solucin. La ecuacin auxiliar est dada por:

m2 + bm + 64 = 0
b b2 256 b
b2 256
m=
=
.
2
2
2
(a.) b = 0. En este caso las raices de la ecuacin son complejas (puras) es decir:
16i
= 8i
2
m2 = 8i
m1 =
La ecuacin de movimiento (subamortiguado) est dada por

y(t) = c1 cos(8t) + c2 sin(8t).
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.
Por lo tanto, la ecuacin de movimiento de la masa es:

y(t) = cos(8t) = cos(8t + )
de donde A =
q
c21 + c22 = 1 es el factor de amortiguamiento.
Para hallar el ngulo , hacemos tan( ) = cc12 el cual es indeterminado, ademas es

un ngulo en el primer cuadrante, as = 2 . El grfico que describe el movimiento
de la masa se da en la Figura 4.7:
137

Figura 4.7: y(t) = cos 8t + 2 .
(b.) b = 10. Las races de la ecuacin son complejas conjugadas.
m1 = 5 + 39i
m2 = 5 39i.
La ecuacin del movimiento (subamortiguado) est dada por
y(t) =e5t (c1 cos( 39t) + c2 sin( 39t))
y0 (t) =e5t (c1 (5 cos( 39t) 39 sin( 39t))
+ c2 ( 39 cos( 39t) 39 sin( 39t))).

Usando las condiciones iniciales tenemos:
1
=y(0) = c1
2
2 =y0 (0) = 5c1 + 39c2 ,

de donde se tiene c1 = 21 y c2 =
este caso es:
y(t) =e5t
39
78 .
Por lo tanto, la ecuacin del movimiento en
1
39
cos( 39t)
sin( 39t)
2
78
78 5t
=
e (sin( 39t) + ).
39
Para hallar el ngulo , hacemos tan( ) = cc12 y tenemos que = 0, 15 rad. El

grfico que describe el movimiento de la masa se da en la Figura 4.8:
138
Figura 4.8: y(t) =
78 5t
(sin( 39t) + ).
39 e
(c.) b = 16. En este caso las races son reales y se repiten m1 = m2 = 8. La ecuacin del
movimiento (crticamente amortiguado) es
y(t) =c1 e8t + c2te8t

y0 (t) = 8c1 e8t + c2 e8t 8c2tet .
Usando las condiciones iniciales tenemos
1
=y(0) = c1
2
2 =y0 (0) = 8c1 + c2 .
As, c1 = 12 y c2 = 2. La ecuacin finalmente queda como
1
y(t) = e8t 2te8t .
2
El grfico que describe el movimiento de la masa se da en la Figura 4.9:
139
Figura 4.9: y(t) = 12 e8t 2te8t .
(d.) b = 20. Las races son reales y diferentes m1 = 4, m2 = 16. La ecuacin del
movimiento (sobreamortiguado) es
y(t) =c1 e4t + c2 e61t

y0 (t) = 4c1 e4t 16c2 e16t .
Resolvemos el problema de valor inicial usando las condiciones iniciales:
1
=y(0) = c1 + c2
2
2 =y0 (0) = 4c1 16c2 .
Resolviendo el sistema tenemos c1 = 12 y c2 = 0. La ecuacin del movimiento queda
de la forma:
1
y(t) = e4t .
2
El grfico que describe el movimiento de la masa se da en la Figura 4.10:
140
Figura 4.10: y(t) = 12 e4t .

Ejemplo 4.21. Un resorte helicoidal se estira 32 pulgadas al ser colgado de l un peso
de 2 libras y luego entra en equilibrio. Se le da ahora un tirn de 1 pie y se libera. Si el
resorte se sumerge en un medio cuyo coeficiente de resistencia es 1/2, hallar la ecuacin de
movimiento del objeto. Asuma que la fuerza de resistencia es proporcional a la velocidad.
Solucin. La segunda ley de Newton, debido al factor de resistencia, en este caso queda
como
1 dy
d2y
.
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,
3
g = 32 pies/seg2
la masa del objeto por lo tanto es

m=
1
ky
= .
g
16
De donde tenemos la ecuacin diferencial

1 d 2 y 1 dy 3
+
+ y = 0,
16 dt 2 2 dt 4
d2y
dy
+
8
+ 12y = 0
dt 2
dt
cuya solucin es
y(t) = c1 e2t + c2 e6t ,
y0 (t) = 2c1 e2t 6c2 e6t .
141
Las condiciones iniciales son y(0) = 1, y0 (0) = 0. Usando las condiciones iniciales nos
queda el sistema
1 =c1 + c2
0 = 2c2 6c2
cuya solucin es c1 =
3
2
y c2 = 21 . Luego, ecuacin del movimiento del objeto es

1
3
y(t) = e2t e6t .
2
2
Ejemplo 4.22. Si el coeficiente de resistencia en el ejemplo anterior es

(a.) La ecuacin del movimiento del sistema.
(b.) El factor de amortiguamiento.
(c.) La amplitud amortiguada del movimiento.
(d.) El periodo amortiguado del movimiento.
(e.) La frecuencia amortiguada del movimiento.
(f.) La constante del tiempo.
Solucin.
3
8
en vez de 12 , hallar:
(a.) La ecuacin diferencial del movimiento en este caso es

1 d 2 y 3 dy 3
+
+ y = 0,
16 dt 2 8 dt 4
La solucin es
d2y
dy
+
6
+ 12y = 0.
dt 2
dt
y(t) = Ae3t sin( 3t + )
y la velocidad es
y0 (t) = 3Ae3t sin( 3t + ) + 3Ae3t cos( 3t + ).

Las condiciones iniciales son y(0) = 1, y0 (0) = 0. De donde tenemos el sistema
1 =A sin( )
0 = 3A sin( ) + 3A cos( ).
Luego, A =
1
sin( )
0 = 3 + 3 arctan( ),
de donde tenemos
3
arctan( ) = = 3,
3
7
1
o
, sin( ) = .
6
6
2
1
Si asumimos sin( ) = 2 , tenemos que A = 2 y la ecuacin del movimiento queda

3t
.
y(t) = 2e sin
3t +
6
=
142
(b.) El factor de amortiguamiento es e3t .

(c.) La amplitud amortiguada del movimiento es 2e3t .
2
seg.
(d.) El periodo amortiguado del movimiento es
3
(e.) La frecuencia amortiguada del movimiento 3 rad/seg
(f.) La constante del tiempo = 13 seg.
3
2
ciclos por segundo.

Ejercicios
1. Una partcula se mueve en linea recta segn la ley
d2x
dx
+ 2r + x = 0,
2
dt
dt
2.
3.
4.
5.
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?.
Una partcula ejecuta un movimiento armnico amortiguado. en 10 segundo el factor
de amortiguamiento decrece en un 80 %. Su periodo de amortiguamiento es 2 seg.
Hallar la ecuacin diferencial del movimiento.
Un peso de 16 libra estira un resorte helicoidal 3/5 pies. El coeficiente de resistencia
del resorte es 8. Despus que el resorte regresa al origen, este es estirado 3 pulgadas
adicionales y luego se libera. Halle la ecuacin del movimiento.
En el ejercicio anterior cambie el coeficiente de resistencia a 10, luego halle la
ecuacin del movimiento.
La frecuencia natural de un resorte es 1 ciclo por segundo. Despus de que el resorte
se sumerge en un medio resistente su frecuencia se reduce a 23 de ciclo por segundo.
a) Cul es el factor de amortiguamiento?.
b) Cul es la ecuacin diferencial del movimiento?.
Solucin
1.
a) subamortiguado si 0 < r < 1, crticamente amortiguado si r = 1 y sobreamortiguado si r > 1.
143
b)
x(t) =(c1 + c2t)et .
1
x(t) =Ae 2 t sin( 3/2t + ).
1
x(t) =Ae 2 t sin( 3/2t + ).
x(t) =(c1 + c2t)et .
3)t
x(t) =c1 e(2+
2.
3.
4.
5.
3)t
+ c2 e(2
c) 35 .
d) No.
d2y
+ 0,322 dy
= 0.
dt + 9.896y
dt 2

1
8t
y(t) = e
4 + 2t .
1 4t
1 16t
y(t) = 3 e 12
e
.
4,68t
a) e
.
00
b) y + 9, 37y0 + 39, 5y = 0.
ndice alfabtico
Abel
frmula de, 94
anulador
mtodo, 107
aplicaciones
ecuacin de primer orden
crecimiento de poblaciones, 57
ley de enfriamiento, 72
mtodo del carbono 14, 65
mezclas, 67
trayectorias ortogonales, 85
ecuacin de segundo orden
movimiento armnico amortiguado,
132
movimiento armnico simple, 124
Bessel
ecuacin de, 124
carbono 14, 66
Cauchy
problema de, 17
coeficientes indeterminados, 91, 103
conjunto de soluciones
fundamental, 95
dependencia lineal, 93
derivadas parciales, 26
ecuacin
auxiliar, 100
ecuacin diferencial, 6
autnoma, 57
de Cauchy-Euler, 117
de primer orden, 21
con coeficientes lineales, 43
de Bernoulli, 46
de Lagrange, 52
de Ricatti, 49
exacta, 26
homognea, 39
lineal, 36
separable, 22
exacta
criterio para, 26
solucin, 26
homognea
criterio, 40
lineal de orden superior, 9, 91
con coeficientes constantes, 91, 100,
104
forma normal de, 112
homognea, 91, 100
146
NDICE ALFABTICO
no homognea, 91, 103
solucin general, 95
ordinaria, 7
parcial, 7
factor integrante, 32
mtodos para hallar, 32
familia de curvas, 85
frecuencia natural de movimiento oscilatorio, 125
grado de una ecuacin diferencial, 9
Grammiano, 94, 98
Hook
ley de, 128
independencia lineal, 94
isotopo radioactivo, 65, 67
Laguerre
ecuacin de, 112
Legendre
ecuacin de, 111
Malthus
ley de, 6, 57, 61
Newton
ley de enfriamiento, 72, 74
mecnica de, 80
segunda ley, 6
operador diferencial, 107
orden de una ecuacin diferencial, 8
poblacin, 60
polinomio caracterstico, 100, 118
reduccin de orden, 91, 108
Schdinger
ecuacin de, 112
solucin, 11
explcita, 13
general, 16
implcita, 14
particular, 16
por sustituciones, 39
singular, 16
sudario de Turin, 67
trayectorias ortogonales, 85
variables separables
reduccin a, 54
variacin de parmetros, 91, 112
Wronskiano, 93, 100

EDO Notas de Curso

Cargado por

Copyright:

Formatos disponibles

EDO Notas de Curso

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

EDO Notas de Curso

Cargado por

Copyright:

Formatos disponibles

6

Luz Marina Moya y Edixon Rojas

c 2014 Luz Marina Moya & Edixon Rojas

Solucin de una Ecuacin Diferencial.

Solucin General y Particular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

Problema de Valor Inicial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

Ecuaciones Diferenciales de Primer Orden . . . . . . . . . . . . . 21

Ecuaciones Diferenciales No Lineales

Ecuaciones Diferenciales Lineales

Soluciones por Sustituciones

Ecuaciones Diferenciales Homogneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39

Ecuaciones con Coeficientes Lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43

Reduccin a Separacin de Variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54

Aplicaciones EDO de Primer Orden . . . . . . . . . . . . . . . . . . . . . . 57

Modelo Logstico de Crecimiento Poblacional . . . . . . . . . . . . . . . . . . . . . . 61

Ecuaciones Diferenciales Lineales de Orden Superior . 91

Ecuaciones lineales Homogneas con Coeficientes Constantes 100

Mtodo de Coeficientes Indeterminados

Movimiento Armnico Simple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124

En este captulo el lector encontrar el lenguaje bsico de las ecuaciones diferenciales

embargo ya nos hemos encontrado con el problema de resolver ecuaciones diferenciales en

donde k es la constante de propocionalidad.

1.1 Definiciones y Ejemplos

Ejemplo 1.3. El modelo de Hodgkin-Huxley para pulsos neuronales. En neurobiologa, la

donde C es la capacitancia de la membrana, Ie la corriente externa e Ii la corriente interna.

Mientras Ie es tratado como un parmetro externo, la dinmica interna depende en gran

donde x es la variable dependiente y t la variable independiente, es una ecuacin diferencial

donde x y y son las conductividades trmicas del solido en las direcciones x y y, T es la

donde F es una funcin de n + 2 variables. Si en la expresin explicita de la ecuacin (1.7)

1.1 Definiciones y Ejemplos

Ejemplo 1.10. A continuacin explicaremos un poco el ejemplo anterior e indicaremos en

3. La ecuacin diferencial ddt 3x = sin(x) es una ecuacin diferencial ordinaria de orden 3,

Otros ejemplos de ecuaciones diferenciales no lineales son los siguientes:

+ y2 = 0, tambin son ejemplos

1.2 Solucin de una Ecuacin Diferencial.

Solucin de una Ecuacin Diferencial.

Verificar que la funcin

es solucin de la ecuacin, donde c es una constante arbitraria.

en algn intervalo I, donde g(x)dx es cualquier antiderivada de g(x).

Ejemplo 1.19. Comprobar que la frmula y =

sin embargo no es una solucin de esta.

Ejemplo 1.20. La ecuacin diferencial de primer orden

no tiene solucin en el intervalo (1, 1).

1.2 Solucin de una Ecuacin Diferencial.

Ejemplo 1.24. Sea C > 0 un nmero real. Verificar que la relacin

determina una solucin implcita de la ecuacin diferencial

Solucin: Diferenciando implcitamente la ecuacin (1.13) obtenemos

de donde, dividiendo todo entre 2, obtenemos

1.2 Solucin de una Ecuacin Diferencial.

determina una solucin implcita de la ecuacin diferencial

Solucin: El Teorema de la Funcin Implcita garantiza la existencia de una funcin

Solucin General y Particular

Solucin: Derivando la expresin para y obtenemos y0 = 4kxe2x , es decir, y0 = 4xy.

Obs En el ejemplo anterior, la familia de soluciones dada es en realidad la solucin general

Ejemplo 1.27. Demostrar que la funcin y = 3 3ex

es una solucin particular de la

1.2 Solucin de una Ecuacin Diferencial.

Ejemplo 1.28. Demostrar que la familia de funciones y = cx + c2 , definidas en toda la recta

y que la funcin y1 = x4 es una solucin singular.

Problema de Valor Inicial

g) (ex sin(y) + tan(y))dx+ (ex cos(y)