Ecuaciones Diferenciales

Por
Eider Leandro Arcila Dager1

Primera Edición

2 de julio de 2020

1 Profesor investigador. Correos electrónicos: arcilaeider@[Link], earcila@[Link], [Link]@[Link],

© 2014-2020 Eider Leandro Arcila Dager

Cualquier persona que reciba una copia de este trabajo, tiene concebido el permiso para:

 Visualizar este trabajo en un dispositivo electrónico.

 Imprimir o fotocopiar este trabajo.
 Distribuir a terceros una copia de este trabajo, siempre y cuando no sufra
modificaciones en su totalidad, con 92 páginas en total.

12 de enero de 2020.
Sobre el autor
Eider Leandro Arcila Dager es caleño, nació el 31/5/1985, desde el 2012 es Matemático de la Universidad
del Valle en Cali-Colombia, en ese mismo año realizó estudios de iniciación científica en probabilidad y
combinatoria en el IMPA2 en Rio de Janeiro-Brasil y actualmente adelanta estudios de maestría en
enseñanza de la matemática en la UTP3 en Pereira-Colombia y es candidato a Doctor en Ciencias
Pedagógicas de la UCPEJV4 en La Habana-Cuba. En el 2012 fue Profesor Catedrático en la UNIVALLE5.
Desde el 2013 es Profesor Investigador UNIAJC6, desde el 2016 es profesor catedrático de la FUC7 y
desde el 2020 es profesor catedrático del FCECEP8. Además pertenece al grupo de investigación en
pedagogía GIP y al grupo de investigación y simulación en ciencias básicas aplicadas GISCBA. Su
principal área de investigación es la Didáctica de la Matemática, Ecuaciones Diferenciales y el Análisis
Matemático y su línea de investigación es el modelamiento matemático.

Sobre la primera edición

Esta guía fue elaborada a partir de las notas de clase de los cursos de Matemáticas III, Ecuaciones
Diferenciales y Matemáticas Especiales orientadas en UNIAJC y FUC.
Para seguir este libro no es necesario cualquier conocimiento previo en ecuaciones diferenciales. Los
pre-requisitos son matemática básica, cálculo diferencial e integral.

2 Instituto de Matemática Pura y Aplicada, Rio de Janeiro-Brasil

3 Universidad Tecnológica de Pereira, Pereira-Colombia
4 Universidad de Ciencias Pedagógicas Enrique José Varona, La Habana-Cuba
5 Universidad del Valle, Cali-Colombia
6 Institución Universitaria Antonio José Camacho, Cali-Colombia
7 Fundación Universitaria Católica “Lumen Gentium”, Cali-Colombia
8 Fundación Centro Colombiano de Estudios Profesionales, Cali-Colombia
Agradecimientos

Quiero agradecer inmensamente a todos los compañeros profesores del Departamento de Ciencias
Básicas de la Institución Universitaria Antonio José Camacho y de la Fundación Universitaria Católica en
especial a:

 Carlos Arturo Muñoz Vargas

 Víctor Manuel Uribe Villegas
 Emiliano Grueso Cárdenas
 Amarilis Báez Bravo
 Héctor García Rodríguez
 Sandra Esther Suárez Chávez
 Milton Fabián Castaño Muñoz
 Steev Romero Agredo
 Gil Ernesto Daza Pérez
 Alejandro Orejuela Benavides
 Mario German Domínguez Pombo
 Daira Patricia Rodríguez Morales
 Jackeline Valencia Giraldo

Además a mis profesores de posgrado

 Aurelio Quintana Valdéz
 Luis Fernando Plaza Gálvez
Para Kika, Mokita, Alejo, Wancha, Pisky, Karlotta y Simon
 Índice
1. Generalidades, soluciones y Problemas de valor inicial (PVI) ..................................................... 1
1.1. Ecuaciones diferenciales y sus soluciones............................................................................. 1
1.2. Problemas de valor inicial ......................................................................................................... 2
2. Ecuaciones diferenciales de primer orden ..................................................................................... 4
2.1. Variables separables ................................................................................................................. 5
2.1.1. Modelo de crecimiento poblacional y desintegración radiactiva .................................... 7
2.1.2. Ley del enfriamiento o calentamiento de Newton .......................................................... 11
2.2. Ecuaciones diferenciales lineales .......................................................................................... 13
2.2.1. Modelo de mezclas ............................................................................................................ 16
2.3. Soluciones por sustitución ..................................................................................................... 20
2.3.1. Ecuaciones homogéneas de primer orden ..................................................................... 20
2.3.2. Ecuación de Bernoulli ....................................................................................................... 24
2.4. Ecuaciones exactas ................................................................................................................. 26
3. Ecuaciones diferenciales de orden superior................................................................................ 31
3.1. Ecuaciones diferenciales lineales homogéneas ................................................................... 31
3.1.1. Sistema Masa-Resorte, movimiento libre no amortiguado ............................................ 34
3.1.2. Sistema Masa-Resorte, movimiento libre amortiguado ................................................. 37
3.2. Ecuaciones diferenciales lineales no homogéneas .............................................................. 40
3.2.1. Método de los coeficientes indeterminados ................................................................... 40
3.2.2. Sistemas Masa-Resorte: Movimiento forzado ................................................................ 43
4. La transformada de Laplace .......................................................................................................... 45
4.1. Solución de problemas de valor inicial .................................................................................. 45
4.2. Transformada de Laplace de funciones periódicas .............................................................. 50
4.3. Traslación en el eje 𝒕................................................................................................................ 52
4.4. Más sobre la transformada de Laplace .................................................................................. 55
EJERCICIOS DE ECUACIONES DIFERENCIALES .............................................................................. 60
1. Generalidades, soluciones y PVI................................................................................................... 61
2. Ecuaciones diferenciales de primer orden ................................................................................... 64
2.1. Variables separables ............................................................................................................... 64
2.1.1. Modelo de crecimiento poblacional y desintegración radiactiva .................................. 66
2.1.2. Ley de enfriamiento o calentamiento de Newton ........................................................... 69
2.2. Ecuaciones diferenciales lineales .......................................................................................... 71
 2.2.1. Modelo de mezclas ............................................................................................................ 74
2.3. Soluciones por sustitución ..................................................................................................... 75
2.4. Ecuaciones exactas ................................................................................................................. 75
3. Ecuaciones diferenciales de orden superior................................................................................ 76
3.1. Ecuaciones diferenciales lineales homogéneas ................................................................... 76
3.1.1. Sistema Masa-Resorte: Movimiento con/sin amortiguamiento ..................................... 78
3.2. Ecuaciones diferenciales no homogéneas ............................................................................ 79
3.2.1. Sistema Masa-Resorte: Movimiento forzado .................................................................. 80
4. La transformada de Laplace .......................................................................................................... 81
4.1. Solución de problemas de valor inicial .................................................................................. 81
4.2. Transformada de Laplace de funciones periódicas .............................................................. 83
4.3. Traslación eje 𝒕 ......................................................................................................................... 84
4.4. Más sobre la transformada de Laplace .................................................................................. 85
 CAPÍTULO 1
1. Generalidades, soluciones y Problemas
de valor inicial (PVI)
1.1. Ecuaciones diferenciales y sus soluciones

Una Ecuación Diferencial de primer orden es una igualdad que posee derivadas de orden uno o primera
derivada. Por ejemplo, las expresiones
𝑑𝑦
a) 𝑒 𝑥 𝑑𝑥 = 3𝑥
𝑑𝑠
b) = 𝑘(𝑠 − 𝑟)
𝑑𝑟
2 2
c) 𝑥 𝑦 𝑑𝑦 = (𝑦 + 2)𝑑𝑥
son ecuaciones diferenciales de primer orden. En la ecuación diferencial a), la variable dependiente es 𝑦
y la variable independiente es 𝑥, en la b), la variable dependiente es 𝑠 y la independiente es 𝑟. En la
ecuación diferencial c) no se puede establecer el tipo de variable.
La expresión 𝑥𝑦 ′′ − 7𝑥𝑦 ′ + 15𝑦 = 0 es una ecuación diferencial de segundo orden cuya variable
dependiente es 𝑦 y variable independiente 𝑥.
Soluciones de una Ecuación Diferencial: El objetivo de resolver una ecuación diferencial es encontrar
una función 𝒚(𝒙) (si 𝑦 es la variable dependiente y 𝑥 la independiente) llamada “solución general”, que
satisfaga la ecuación diferencial.
Ejemplo 1: Verificar que la función
1
𝑦(𝑥) =
1 + 𝑐𝑒 −𝑥
es una solución general de la ecuación diferencial de primer orden
𝑑𝑦
= 𝑦 − 𝑦2
𝑑𝑥
Solución: Verifiquemos que la función 𝑦(𝑥) satisface que su derivada es igual a la diferencia entre ella y
su cuadrado.
𝑑𝑦 𝑐𝑒 −𝑥
=
𝑑𝑥 (1 + 𝑐𝑒 −𝑥 )2
Por lo tanto

1
 2
𝑐𝑒 −𝑥 1 1
= − ( )
(1 + 𝑐𝑒 −𝑥 )2 ⏟
⏟ 1 + 𝑐𝑒 −𝑥 ⏟1 + 𝑐𝑒 −𝑥
𝑑𝑦 𝑦 𝑦2
𝑑𝑥

desarrollando el lado derecho

𝑐𝑒 −𝑥 1 1 1 + 𝑐𝑒 −𝑥 − 1 𝑐𝑒 −𝑥
= − = = ∎
(1 + 𝑐𝑒 −𝑥 )2 1 + 𝑐𝑒 −𝑥 (1 + 𝑐𝑒 −𝑥 )2 (1 + 𝑐𝑒 −𝑥 )2 (1 + 𝑐𝑒 −𝑥 )2

Ejercicio en la guía 1: Verificar que la función 𝑦 = 𝑐1 𝑒 2𝑥 + 𝑐2 𝑥𝑒 2𝑥 es una solución de la ecuación diferencial

𝑑2 𝑦 𝑑𝑦
− 4 𝑑𝑥 + 4𝑦 = 0
𝑑𝑥 2

1.2. Problemas de valor inicial

Diversas problemáticas que se modelan a través de ecuaciones diferenciales deben satisfacer

condiciones iniciales sobre su solución o sus derivadas. A esta situación se le denomina problema de
valor inicial (PVI) y en general se describe de la siguiente forma
𝑑𝑛 𝑦
= 𝑓(𝑥, 𝑦, 𝑦 ′ , … , 𝑦 (𝑛−1) )
{ 𝑑𝑥 𝑛

𝑦(𝑥0 ) = 𝑦0 , 𝑦 ′ (𝑥0 ) = 𝑦1 , … , 𝑦 (𝑛−1) (𝑥0 ) = 𝑦𝑛−1

De acuerdo a la expresión anterior, un problema de valor inicial de primer orden sería
𝑑𝑦
= 𝑓(𝑥, 𝑦)
{𝑑𝑥
𝑦(𝑥0 ) = 𝑦0

Y uno de segundo orden:

𝑑2𝑦
{ = 𝑓(𝑥, 𝑦, 𝑦 ′ )
𝑑𝑥 2
𝑦(𝑥0 ) = 𝑦0 , 𝑦 ′ (𝑥0 ) = 𝑦1

2
Geométricamente un PVI de primer orden significa que la gráfica de la solución de la ecuación diferencial
𝑦(𝑥) pase por el punto (𝑥0 , 𝑦0 ) y 𝑥0 debe estar definido en un algún intervalo 𝐽.

Ejercicio en la guía 2: Describir geométricamente el significado de un PVI de segundo orden.

R/ Descripción por parte del estudiante

Ejemplo 2: Encuentre una solución del PVI de primer orden que consiste en la ecuación diferencial 𝑦 ′ =
𝑦 − 𝑦 2 y la condición inicial 𝑦(−1) = 2, donde 𝑦 = 1/(1 + 𝑐1 𝑒 −𝑥 ) es una familia uniparamétrica de
soluciones de la ecuación diferencial.
Solución: La solución 𝑦 debe pasar por el punto (−1,2), por lo tanto si 𝑦 = 1/(1 + 𝑐1 𝑒 −𝑥 ), entonces
1 1 1 1 1
2= → 2= → 1 + 𝑐1 𝑒 = → 𝑐1 𝑒 = − → 𝑐1 = −
1 + 𝑐1 𝑒 −(−1) 1 + 𝑐1 𝑒 2 2 2𝑒
Concluimos que la solución 𝑦 que satisface la ecuación diferencial y pasa por el punto (−1,2) es
1
𝑦= ∎
1
1 + (− 2𝑒) 𝑒 −𝑥

Ejemplo 3: Use el hecho de que 𝑦 = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 −𝑥 es una familia biparamétrica de soluciones de la

ecuación diferencial 𝑦 ′′ − 𝑦 = 0, para llegar a una solución del PVI formado por la ecuación diferencial y
las condiciones iniciales 𝑦(0) = 1 , 𝑦 ′ (0) = 2.
Solución: La solución 𝑦 debe pasar por el punto (0,1), por lo tanto si 𝑦 = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 −𝑥 , entonces
1 = 𝑐1 𝑒 0 + 𝑐2 𝑒 −0 → 1 = 𝑐1 + 𝑐2
Además la pendiente de la recta tangente en 𝑥 = 0 es 2. Calculando la primera derivada de 𝑦:
𝑦 ′ = 𝑐1 𝑒 𝑥 − 𝑐2 𝑒 −𝑥
Por lo tanto
2 = 𝑐1 𝑒 0 − 𝑐2 𝑒 −0 → 2 = 𝑐1 − 𝑐2
De las dos ecuaciones anteriores, obtenemos un sistema de ecuaciones lineales 2 × 2:
𝑐1 + 𝑐2 = 1
{
𝑐1 − 𝑐2 = 2
3
Sumando ambas ecuaciones, obtenemos 2𝑐1 = 3, donde 𝑐1 = 3/2. Reemplazando este valor de 𝑐1 en la
3
ecuación 1 obtenemos que 𝑐2 = 1 − = −1/2
2

Concluimos que la solución 𝑦 que satisface la ecuación diferencial y las condiciones iniciales es:
3 1
𝑦 = 𝑒 𝑥 − 𝑒 −𝑥 ∎
2 2
Ejercicio en la guía 3: Encuentre una solución del PVI de segundo orden que consiste en la ecuación
diferencial 𝑥 ′′ + 𝑥 = 0 y las condiciones iniciales 𝑥(𝜋/6) = 1/2 y 𝑥′(𝜋/6) = 0, donde 𝑥 = 𝑐1 cos 𝑡 + 𝑐2 sin 𝑡
es una familia biparamétrica de soluciones de la ecuación diferencial.

√𝟑 𝟏
R/ 𝒙(𝒕) = 𝐜𝐨𝐬 𝒕 + 𝟒 𝐬𝐢𝐧 𝒕
𝟒

CAPÍTULO 2
2. Ecuaciones diferenciales de primer
orden
Iniciamos el largo y hermoso camino de las ecuaciones diferenciales definiendo el concepto de ecuación
diferencial de primer orden y con ello diferentes métodos para su solución. Se plantean diferentes
aplicaciones que se modelan mediante ecuaciones diferenciales de primer orden como los modelos de
crecimiento y enfriamiento entre otros.

4
 2.1. Variables separables

Un primer método para encontrar soluciones de una ecuación diferencial se denomina el método de las
variables separables.
Una ecuación diferencial se denomina de variables separables si se puede expresar de la forma

𝑑𝑦
= 𝑓(𝑥) ∙ 𝑔(𝑦)
𝑑𝑥
𝑑𝑦
Ejemplo 3: Determinar si la ecuación diferencial 𝑒 𝑥 𝑦 𝑑𝑥 = 𝑒 −𝑦 + 𝑒 −2𝑥−𝑦 es de
variables separables.

Solución: Observemos que 𝑒 −2𝑥−𝑦 =𝑒 −2𝑥 𝑒 −𝑦 y usando factor común obtenemos:

𝑑𝑦
𝑒𝑥𝑦 = 𝑒 −𝑦 + 𝑒 −2𝑥 𝑒 −𝑦
𝑑𝑥
𝑑𝑦
𝑒𝑥𝑦 = 𝑒 −𝑦 (1 + 𝑒 −2𝑥 )
𝑑𝑥
𝑑𝑦
Despejando la expresión 𝑑𝑥

𝑑𝑦 𝑒 −𝑦 1 + 𝑒 −2𝑥
= ∙
𝑑𝑥 ⏟𝑦 ⏟ 𝑒𝑥
𝑔(𝑦) 𝑓(𝑥)

Después de verificar que una ecuación diferencial es de variables separables, se procede realizar una
“separación” de las variables (las expresiones que contienen la variable 𝑦 se ubican en el lado izquierdo
de la igualdad y las expresiones que contienen la variable 𝑥 se ubican en el lado derecho de la igualdad)9.
Ejemplo 4: (Separación de las variables). Separar las variables de la ecuación diferencial

𝑑𝑥 𝑥 2 𝑦 2
=
𝑑𝑦 𝑥 + 1

Solución: Las expresiones que contienen la variable 𝑥 se ubican en el lado izquierdo

𝑥 + 1 𝑑𝑥
= 𝑦2
𝑥 2 𝑑𝑦

Finalmente, trasladamos la expresión “𝑑𝑦” al lado derecho

𝑥+1
𝑑𝑥 = 𝑦 2 𝑑𝑦 ∎
𝑥2

9Si la ecuación diferencial posee otro tipo de variables (𝑧, 𝑤, 𝑡, 𝜃, … 𝑒𝑡𝑐), se realiza la misma separación usando estas
variables.
5
Para finalizar la resolución de una ecuación diferencial por el método de variables separables, después
de la separación, se procede a integrar con respecto a las variables indicadas en cada lado de la
ecuación.

Ejemplo 5: Resuelva la ecuación diferencial 𝑐𝑠𝑐𝑦 𝑑𝑥 + 𝑠𝑒𝑐 2 𝑥 𝑑𝑦 = 0.

Solución: Separando las variables
𝑑𝑦 𝑑𝑥
=−
𝑐𝑠𝑐𝑦 𝑠𝑒𝑐 2 𝑥
1 1
Integrando ambos lados y utilizando las relaciones trigonométricas básicas = 𝑠𝑒𝑛𝑦 y = 𝑐𝑜𝑠 2 𝑥,
𝑐𝑠𝑐𝑦 𝑠𝑒𝑐 2 𝑥
tenemos:

∫ 𝑠𝑒𝑛𝑦 𝑑𝑦 = − ∫ 𝑐𝑜𝑠 2 𝑥 𝑑𝑥

1
Usando la identidad del ángulo doble 𝑐𝑜𝑠 2 𝑥 = 2 (1 + 𝑐𝑜𝑠2𝑥)

1
∫ 𝑠𝑒𝑛𝑦 𝑑𝑦 = − ∫(1 + 𝑐𝑜𝑠2𝑥) 𝑑𝑥
2
1 1
−𝑐𝑜𝑠𝑦 = − (𝑥 + 𝑠𝑒𝑛2𝑥) + 𝐾
2 2
Despejando 𝑦, tenemos:
1 1
𝑦(𝑥) = 𝑐𝑜𝑠 −1 [ (𝑥 + 𝑠𝑒𝑛2𝑥) + 𝐾] ∎
2 2
Ejercicio en la guía 4: Verificar que la solución del ejemplo 5 satisface la ecuación diferencial

𝑑𝑦
Ejemplo 6: Resolver la ecuación diferencial 𝑥 2 𝑑𝑥 = 𝑦 − 𝑥𝑦.

Solución: Factorizando el lado derecho y separando variables, obtenemos:

6
 𝑑𝑦 1 − 𝑥
= 2 𝑑𝑥
𝑦 𝑥
Integrando ambos lados
𝑑𝑦 1−𝑥
∫ = ∫ 2 𝑑𝑥
𝑦 𝑥

1 𝑥 𝑥 −1 1
𝑙𝑛|𝑦| = ∫ 2 𝑑𝑥 − ∫ 2 𝑑𝑥 = − 𝑙𝑛|𝑥| + 𝐾 = − − 𝑙𝑛|𝑥| + 𝐾
𝑥 𝑥 −1 𝑥
Despejando la variable 𝑦, obtenemos
1
1 𝑒 −𝑥
𝑦(𝑥) = 𝑒 (−𝑥−𝑙𝑛|𝑥|+𝐾) =𝐾 ∎
𝑥
Ejercicio en la guía 5: Resolver la ecuación diferencial 2𝑦(𝑥 + 2)𝑑𝑦 = 𝑥 𝑑𝑥

R/ 𝒚𝟐 = 𝒙 + 𝟐 − 𝟐 𝐥𝐧|𝒙 + 𝟐| + 𝑪

Pérdida de una solución: Se debe tener precaución al separar las variables. Un ejemplo es la ecuación
𝑑𝑦 2+𝑘𝑒 4𝑥
diferencial = 𝑦 2 − 4 cuya solución general es 𝑦 = 1−𝑘𝑒 4𝑥. Sin embargo, 𝑦 = −2 es una solución
𝑑𝑥
constante y no se puede obtener de la solución general mediante la variación de la constante 𝑘. Esto
quiere decir que se ha perdido la solución durante los procesos de separación de variables.

2.1.1. Modelo de crecimiento poblacional y desintegración radiactiva

El modelo
𝑑𝑃
= 𝑘𝑃
{ 𝑑𝑡
𝑃(𝑡0 ) = 𝑃0

donde 𝑘 es una constante de proporcionalidad, sirve para resolver fenómenos de crecimiento o

decaimiento.
La ecuación diferencial es de variables separables por lo tanto

7
 𝑑𝑃 𝑑𝑃
= 𝑘𝑃 → = 𝑘 𝑑𝑡
𝑑𝑡 𝑃
integrando ambos lados
𝑑𝑃
∫ = ∫ 𝑘 𝑑𝑡
𝑃
ln 𝑃 = 𝑘𝑡 + 𝑐

𝑒 ln 𝑃 = 𝑒 𝑘𝑡+𝑐

𝑃(𝑡) = 𝑐𝑒 𝑘𝑡

Si utilizamos la condición inicial 𝑃(𝑡0 ) = 𝑃0 :

𝑃(𝑡0 ) = 𝑐𝑒 𝑘𝑡0 = 𝑃0

entonces despejando 𝑐, tenemos

𝑃0
𝑐= = 𝑃0 𝑒 −𝑘𝑡0
𝑒 𝑘𝑡0
La solución del problema de valor inicial es

𝑃(𝑡) = 𝑃0 𝑒 −𝑘𝑡0 𝑒 𝑘𝑡

𝑷(𝒕) = 𝑷𝟎 𝒆𝒌(𝒕−𝒕𝟎 )

Un caso particular es trabajar en 𝑡0 = 0, entonces la solución queda de la forma

𝑃(𝑡) = 𝑃0 𝑒 𝑘𝑡

Ejemplo 7: Se sabe que la población de una comunidad crece proporcionalmente al número de personas
presentes en el tiempo 𝑡. Si la población inicial 𝑃0 se duplicó en 6 años, ¿En cuánto tiempo se triplicará?

Solución: El problema de valor inicial a resolver es

𝑑𝑃
= 𝑘𝑃
{ 𝑑𝑡
𝑃(6) = 2𝑃0

Para asi, encontrar el valor de 𝑡 que satisfaga la condición 𝑃(𝑡) = 3𝑃0 .Para esto debemos encontrar la
función 𝑃(𝑡) que satisfaga la ecuación diferencial y la condición inicial.

La solución general de la ecuación diferencial es

𝑃(𝑡) = 𝑐𝑒 𝑘𝑡

como 𝑃(0) = 𝑃0 porque 𝑃0 es la población inicial, entonces

𝑃(0) = 𝑐𝑒 𝑘(0) = 𝑃0

donde 𝑐 = 𝑃0 . La solución queda

8
 𝑃(𝑡) = 𝑃0 𝑒 𝑘𝑡

con la condición 𝑃(6) = 2𝑃0 calculamos el valor de 𝑘

𝑃(6) = 𝑃0 𝑒 𝑘(6) = 2𝑃0

simplificando y despejando 𝑘
𝑙𝑛 2
𝑘= ≈ 0,1155
6
La solución queda

𝑃(𝑡) = 𝑃0 𝑒 0,1155𝑡

Debemos encontrar el valor de 𝑡 que cumple que 𝑃(𝑡) = 3𝑃0

𝑃(𝑡) = 𝑃0 𝑒 0,1155𝑡 = 3𝑃0

Simplificando y despejando 𝑡
𝑙𝑛 3
𝑡= ≈ 9,5117 𝑎ñ𝑜𝑠
0,1155
Ejercicio en la guía 6: Un cultivo de bacterias de población 𝑃 crece a un ritmo proporcional a 𝑃. Entre las 7
P.M. y las 8 P.M. la población se triplica. ¿A qué hora será ciento cincuenta veces mayor que la que había
a las 7 P.M.?

R/ 𝑨𝒑𝒓𝒐𝒙. 𝟏𝟏: 𝟑𝟓 𝒑. 𝒎

Ejemplo 8: Inicialmente había 200 miligramos de una sustancia radiactiva. Después de 5 horas la masa
disminuyó 4%. Si la razón de decaimiento, en cualquier momento, es proporcional a la cantidad de
sustancia presente al tiempo 𝑡, determine la cantidad que queda después de 30 horas.

Solución: El problema de valor inicial a resolver es

𝑑𝑃
= 𝑘𝑃
{ 𝑑𝑡
𝑃(0) = 200
Sabemos que su solución es
9
 𝑃(𝑡) = 200𝑒 𝑘𝑡

Como después de 5 horas la masa inicial de 200 mg disminuyó 4% entonces una nueva condición es

𝑃(5) = 192

Por lo tanto

𝑃(5) = 200𝑒 𝑘(5) = 192

Despejando el valor de la constante de proporcionalidad
192
𝑒 5𝑘 =
200
192
5𝑘 = ln ( )
200
192
ln (200)
𝑘= ≈ −0,008164398
5
La solución queda de la forma

𝑃(𝑡) = 200𝑒 −0,008164398𝑡

La cantidad que queda después de 30 horas es

𝑃(30) = 200𝑒 −0,008164398(30) ≈ 156.5515 mg

Observación: La vida media de una sustancia es el tiempo 𝑡𝑣𝑚 que transcurre para que la cantidad inicial de
𝑃0
sustancia se reduzca a la mitad, esto quiere decir, 𝑃(𝑡𝑣𝑚 ) = .
2

Ejercicio en la guía 7: Calcule la vida media de la sustancia del ejemplo 8.

R/ 𝟖𝟒. 𝟗 𝒉𝒐𝒓𝒂𝒔

10
 2.1.2. Ley del enfriamiento o calentamiento de Newton

La rapidez con la que cambia la temperatura de un cuerpo es proporcional a la diferencia entre la

temperatura del cuerpo y la del medio en que está rodeado, que se llama la temperatura ambiente. Si
𝑑𝑇
𝑇(𝑡) es la temperatura del cuerpo en el instante 𝑡, 𝑇𝑎 es la temperatura del medio que lo rodea y es la
𝑑𝑡
rapidez con que cambia la temperatura del cuerpo, entonces la ley de Newton de
Enfriamiento/Calentamiento se escribe matemáticamente
𝑑𝑇
= 𝑘(𝑇 − 𝑇𝑎 )
𝑑𝑡
donde 𝑘 es una constante de proporcionalidad negativa.

Para encontrar la solución 𝑇(𝑡) de la ecuación diferencial, usamos el método de las variables separables
𝑑𝑇
= 𝑘 𝑑𝑡
𝑇 − 𝑇𝑎
integrando ambos lados de la ecuación
𝑑𝑇
∫ = ∫ 𝑘 𝑑𝑡
𝑇 − 𝑇𝑎

ln(𝑇 − 𝑇𝑎 ) = 𝑘𝑡 + 𝑐

Usando propiedades exponenciales y despejando 𝑇

𝑇 − 𝑇𝑎 = 𝑒 𝑘𝑡+𝑐 = 𝑐𝑒 𝑘𝑡

𝑻(𝒕) = 𝑻𝒂 + 𝒄𝒆𝒌𝒕

Ejemplo 9: Un termómetro se cambia de una habitación donde la temperatura es de 72℉ al exterior,

donde la temperatura del aire es de 10℉. Después de medio minuto el termómetro indica 52℉. ¿Cuál es
la lectura del termómetro en 𝑡 = 2 min? ¿Cuánto tiempo le tomará al termómetro alcanzar los 16℉?

Solución: La temperatura ambiente es 𝑇𝑎 = 10℉, por lo tanto:

𝑇(𝑡) = 10 + 𝑐𝑒 𝑘𝑡

La temperatura inicial del termómetro es 72℉, esto quiere decir 𝑇(0) = 72 y después de medio minuto
1
indica 52℉, esto significa 𝑇 (2) = 52.

𝑇(0) = 10 + 𝑐𝑒 𝑘(0) = 72

Concluyendo que 𝑐 = 62, por lo tanto

𝑇(𝑡) = 10 + 62𝑒 𝑘𝑡

11
luego
1 1
𝑇 ( ) = 10 + 62𝑒 𝑘(2) = 52
2
despejando 𝑘
21
𝑘 = 2 𝑙𝑛 ( ) ≈ −0,78
31
Por lo tanto

𝑇(𝑡) = 10 + 62𝑒 −0,78𝑡

La temperatura del termómetro en 𝑡 = 2 min es

𝑇(2) = 10 + 62𝑒 −0,78(2) ≈ 23℉

Para responder la última pregunta, cuánto le toma a la temperatura alcanzar 16℉, planteamos la siguiente
ecuación

16 = 10 + 62𝑒 −0,78𝑡

despejando 𝑡, obtenemos
3
𝑙𝑛 (31)
𝑡= ≈3
−0,78
En aproximadamente 3 min le toma a la temperatura alcanzar 16℉. ∎

Ejercicio en la guía 8: Un cuerpo se calienta a 110℃ y se coloca en aire a 10℃. Tras una hora su
temperatura es de 60℃. ¿Cuánto tiempo tarda todavía en alcanzar los 30℃?

R/ 𝒕 = 𝟏. 𝟑𝟐 𝒉𝒐𝒓𝒂𝒔

12
 2.2. Ecuaciones diferenciales lineales
Una ecuación diferencial se denomina lineal si es de la forma

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝒅𝒚
𝑎𝑛 (𝑥) 𝑛
+ 𝑎 𝑛−1 (𝑥) 𝑛−1
+ ⋯ + 𝒂𝟏 (𝒙) + 𝒂𝟎 (𝒙)𝒚 = 𝒈(𝒙)
𝑑𝑥 𝑑𝑥 𝒅𝒙
Una ecuación diferencial lineal de primer orden es
𝑑𝑦
𝑎1 (𝑥) + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥)
𝑑𝑥
Para resolver la ecuación lineal de primer orden, debemos escribir la ecuación en su forma estándar, esto
es
𝑑𝑦 𝑎0 (𝑥) 𝑔(𝑥)
1 + 𝑦=
𝑑𝑥 𝑎1 (𝑥) 𝑎1 (𝑥)
𝑎 (𝑥) 𝑔(𝑥)
haciendo 𝑎0 (𝑥) = 𝑝(𝑥) y 𝑎 = 𝑞(𝑥), tenemos
1 1 (𝑥)

𝑑𝑦
+ 𝑝(𝑥)𝑦 = 𝑞(𝑥)
𝑑𝑥
El segundo paso es encontrar un factor integrante

𝐹. 𝐼 = 𝑒 ∫ 𝑝(𝑥)𝑑𝑥
Multiplicamos la ecuación por el factor integrante
𝑑𝑦
𝑒 ∫ 𝑝(𝑥)𝑑𝑥 ( + 𝑝(𝑥)𝑦) = 𝑒 ∫ 𝑝(𝑥)𝑑𝑥 𝑞(𝑥)
𝑑𝑥
Notemos que el lado izquierdo de la anterior expresión es la derivada con respecto a 𝑥 (variable
independiente) del factor integrante por 𝑦 (variable dependiente), esto quiere decir
𝑑𝑦 𝑑 ∫ 𝑝(𝑥)𝑑𝑥
𝑒 ∫ 𝑝(𝑥)𝑑𝑥 ( + 𝑝(𝑥)𝑦) = (𝑒 𝑦)
𝑑𝑥 𝑑𝑥
Por lo tanto
𝑑 ∫ 𝑝(𝑥)𝑑𝑥
(𝑒 𝑦) = 𝑒 ∫ 𝑝(𝑥)𝑑𝑥 𝑞(𝑥)
𝑑𝑥
Integrando ambos lados con respecto a 𝑥 (variable independiente)
𝑑 ∫ 𝑝(𝑥)𝑑𝑥
∫ (𝑒 𝑦)𝑑𝑥 = ∫ 𝑒 ∫ 𝑝(𝑥)𝑑𝑥 𝑞(𝑥) 𝑑𝑥
𝑑𝑥

𝑒 ∫ 𝑝(𝑥)𝑑𝑥 𝑦 = ∫ 𝑒 ∫ 𝑝(𝑥)𝑑𝑥 𝑞(𝑥) 𝑑𝑥 + 𝐶

13
despejando la variable dependiente 𝑦

∫ 𝒆∫ 𝒑(𝒙)𝒅𝒙 𝒒(𝒙) 𝒅𝒙 + 𝑪
𝒚=
𝒆∫ 𝒑(𝒙)𝒅𝒙
Esta última expresión es la solución de la ecuación diferencial lineal de primer orden. En el siguiente
ejemplo se hace una descripción paso a paso de la solución de una ecuación diferencial lineal de primer
orden.
𝑑𝑦
Ejemplo 10: Resolver la ecuación diferencial 𝑐𝑜𝑠𝑥 𝑑𝑥 + (𝑠𝑒𝑛𝑥)𝑦 = 1.

Solución: Llevamos la ecuación lineal a su forma estándar

𝑑𝑦 𝑠𝑒𝑛𝑥 1
+ 𝑦=
𝑑𝑥 𝑐𝑜𝑠𝑥 𝑐𝑜𝑠𝑥
𝑠𝑒𝑛𝑥
𝑠𝑒𝑛𝑥 1
Veamos que 𝑝(𝑥) = y 𝑞(𝑥) = 𝑐𝑜𝑠𝑥 , y por tanto el factor integrante es 𝐹. 𝐼 = 𝑒 ∫ 𝑝(𝑥)𝑑𝑥 = 𝑒 ∫𝑐𝑜𝑠𝑥 𝑑𝑥 =
𝑐𝑜𝑠𝑥
𝑒 −𝑙𝑛(𝑐𝑜𝑠𝑥) = (𝑐𝑜𝑠𝑥)−1 = 𝑠𝑒𝑐𝑥. Multiplicamos la ecuación en su forma estándar por el factor integrante:
𝑑𝑦 𝑠𝑒𝑛𝑥 1
𝑠𝑒𝑐𝑥 ( + 𝑦) = 𝑠𝑒𝑐𝑥 ( )
𝑑𝑥 𝑐𝑜𝑠𝑥 𝑐𝑜𝑠𝑥
El lado izquierdo de la anterior expresión es equivalente a:
𝑑𝑦 𝑠𝑒𝑛𝑥 𝑑
𝑠𝑒𝑐𝑥 ( + 𝑦) = (𝑦 ∙ 𝑠𝑒𝑐𝑥)
𝑑𝑥 𝑐𝑜𝑠𝑥 𝑑𝑥
Y por tanto
𝑑 1
(𝑦 ∙ 𝑠𝑒𝑐𝑥) = 𝑠𝑒𝑐𝑥 ( ) = 𝑠𝑒𝑐 2 𝑥
𝑑𝑥 𝑐𝑜𝑠𝑥
Integrando ambos lados con respecto a 𝑥:
𝑑
∫ (𝑦 ∙ 𝑠𝑒𝑐𝑥) 𝑑𝑥 = ∫ 𝑠𝑒𝑐 2 𝑥 𝑑𝑥
𝑑𝑥
𝑦 ∙ 𝑠𝑒𝑐𝑥 = 𝑡𝑎𝑛𝑥 + 𝐾
Finalmente
𝑡𝑎𝑛𝑥 𝐾
𝑦(𝑥) = + = 𝑠𝑒𝑛𝑥 ∙ 𝑠𝑒𝑐 2 𝑥 + 𝐾𝑐𝑜𝑠𝑥 ∎
𝑠𝑒𝑐𝑥 𝑠𝑒𝑐𝑥

La descripción hecha en el ejemplo 10 de la solución de la ecuación diferencial lineal de primer orden

𝑑𝑦
𝑐𝑜𝑠𝑥 𝑑𝑥 + (𝑠𝑒𝑛𝑥)𝑦 = 1. Pero si desde un principio determinamos las funciones 𝑝(𝑥) y 𝑞(𝑥), en el ejemplo

14
 sin 𝑥 1
anterior 𝑝(𝑥) = cos 𝑥 y 𝑞(𝑥) = cos 𝑥 y reemplazamos en la expresión
∫ 𝑒 ∫ 𝑝(𝑥)𝑑𝑥 𝑞(𝑥) 𝑑𝑥+𝐶
𝑦= obtenemos la misma solución del ejemplo 10.
𝑒 ∫ 𝑝(𝑥)𝑑𝑥
𝜋
Ejercicio en la guía 9: Encontrar la solución del ejemplo 10 que satisfaga la condición inicial 𝑦(𝜋) = 2

𝑑𝑦
Ejemplo 11: Resolver el problema de valor inicial + 2𝑥𝑦 = 𝑓(𝑥) ; 𝑦(0) = 2,
𝑑𝑥

𝑥, 0 ≤ 𝑥 < 1
donde 𝑓(𝑥) = { .
0, 𝑥 ≥ 1
Solución: Caso 1: 𝟎 ≤ 𝒙 < 𝟏.
𝑑𝑦
La ecuación se convierte en + 2𝑥𝑦 = 𝑥, vemos que es una ecuación diferencial separable y lineal.
𝑑𝑥
2
Calculando el factor integrante 𝐹. 𝐼 = 𝑒 ∫ 2𝑥 𝑑𝑥 = 𝑒 𝑥 y multiplicando la ecuación por este factor

2 𝑑𝑦 2
𝑒𝑥 ( + 2𝑥𝑦) = 𝑥𝑒 𝑥
𝑑𝑥
Escribiendo la forma equivalente del lado izquierdo e integrando

2 1 𝑥2 1 2
𝑒𝑥 𝑦 = 𝑒 +𝐾 ⟹ 𝑦(𝑥) = + 𝐾𝑒 −𝑥
2 2
1 2 3
Como 𝑦(0) = 2, entonces 𝑦(0) = 2 + 𝐾𝑒 −0 = 2, de aquí, 𝐾 = 2. La solución del caso 1 es

1 3 −𝑥 2
𝑦(𝑥) = + 𝑒
2 2
Caso 2: 𝒙 ≥ 𝟏
𝑑𝑦
La ecuación se convierte en + 2𝑥𝑦 = 0, vemos que es una ecuación diferencial separable y lineal.
𝑑𝑥
2
Calculando el factor integrante 𝐹. 𝐼 = 𝑒 ∫ 2𝑥 𝑑𝑥 = 𝑒 𝑥 y multiplicando la ecuación por este factor

2 𝑑𝑦
𝑒𝑥 ( + 2𝑥𝑦) = 0
𝑑𝑥
Escribiendo la forma equivalente del lado izquierdo e integrando
2 2
𝑒𝑥 𝑦 = 𝐾 ⟹ 𝑦(𝑥) = 𝐾𝑒 −𝑥
La solución de la unión de los dos casos quedaría:
15
 1 3 −𝑥 2
+ 𝑒 , 0≤𝑥<1
𝑦(𝑥) = {2 2
2
𝐾𝑒 −𝑥 , 𝑥≥1

La solución 𝑦(𝑥) debe ser una función continua en 𝑥 ≥ 0, por tanto, debemos encontrar el valor de 𝐾 para
que la función sea continua en 𝑥 = 1, para esto, los limites laterales alrededor de 𝑥 = 1 deben coincidir:
1 3 2 1 3
 𝑙𝑖𝑚− (2 + 2 𝑒 −𝑥 ) = 2 + 2 𝑒 −1
𝑥→1
2
 𝑙𝑖𝑚+ 𝐾𝑒 −𝑥 = 𝐾𝑒 −1
𝑥→1

1 3
Igualando los dos límites laterales, obtenemos que 𝐾 = 2 𝑒 + 2. Finalmente la solución del problema de
valor inicial es
1 3 −𝑥 2
+ 𝑒 , 0≤𝑥<1
𝑦(𝑥) = { 2 2 ∎
1 3 2
( 𝑒 + ) 𝑒 −𝑥 , 𝑥≥1
2 2

𝑑𝑦
Ejercicio en la guía 10: Resolver la ecuación diferencial 𝑥 𝑑𝑥 + 2𝑦 = 4 usando la separación de variables y
el método de las ecuaciones lineales.

2.2.1. Modelo de mezclas

La ecuación diferencial de primer orden que modela la cantidad de sal (en libras) en un tanque al mezclar
dos soluciones que contienen sal a distintas concentraciones es
𝑑𝐴
= (𝑇𝑎𝑠𝑎 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑠𝑡𝑎𝑛𝑐𝑖𝑎) − (𝑇𝑎𝑠𝑎 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑠𝑡𝑎𝑛𝑐𝑖𝑎)
𝑑𝑡

16
Donde cada razón es el producto de la velocidad en (gal/min) y la concentración10 en (lb/gal).

Ejemplo 12: Un tanque contiene 200 litros de un líquido en el que se han disuelto 30 g de sal. A la
salmuera entra 1 g de sal por litro a una razón (velocidad) de 4 L/min; la solución bien mezclada sale del
tanque con la misma razón. Encuentre la cantidad 𝐴(𝑡) de gramos de sal que hay en el tanque al tiempo
𝑡.

Solución: Inicialmente se han disuelto 30 g de sal, por lo tanto la condición inicial es 𝐴(0) = 30. El modelo
de mezclas es
𝑑𝐴
= (𝑇𝑎𝑠𝑎 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑠𝑡𝑎𝑛𝑐𝑖𝑎) − (𝑇𝑎𝑠𝑎 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑠𝑡𝑎𝑛𝑐𝑖𝑎) = 𝑅𝑒𝑛𝑡𝑟𝑎 − 𝑅𝑠𝑎𝑙𝑒
𝑑𝑡
Calculando las razones:
𝑅𝑒𝑛𝑡𝑟𝑎 = (𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑠𝑎𝑙)(𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎)
𝑅𝑒𝑛𝑡𝑟𝑎 = (1𝑔/𝐿)(4𝐿/𝑚𝑖𝑛) = 4 𝑔/𝑚𝑖𝑛

𝑅𝑠𝑎𝑙𝑒 = (𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑠𝑎𝑙)(𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎)

𝐴(𝑡) 𝐴(𝑡)
𝑅𝑠𝑎𝑙𝑒 = ( 𝑔/𝐿) (4𝐿/𝑚𝑖𝑛) = 𝑔/𝑚𝑖𝑛
200 50
Por lo tanto, el problema de valor inicial de mezclas es
𝑑𝐴 𝐴
= 4− ; 𝐴(0) = 30
𝑑𝑡 50
La ecuación diferencial resultante es lineal y su forma estándar es
𝑑𝐴 𝐴
+ =4
𝑑𝑡 50
1
Donde 𝑝(𝑡) = 50 y 𝑞(𝑡) = 4, la solución de la ecuación lineal es

∫ 𝑒 ∫ 𝑝(𝑡)𝑑𝑡 ∙ 𝑞(𝑡)𝑑𝑡 + 𝐶
𝐴(𝑡) =
𝑒 ∫ 𝑝(𝑡)𝑑𝑡
1 1 1
𝑑𝑡
∫ 𝑒 ∫50 ∙ 4 𝑑𝑡 + 𝐶 4 ∫ 𝑒 50𝑡 + 𝐶 200𝑒 50𝑡 + 𝐶 1
𝐴(𝑡) = 1 = 1 = 1 = 200 + 𝐶𝑒 −50𝑡
𝑒 ∫50𝑑𝑡 𝑒 50𝑡 𝑒 50𝑡
Tenemos que
1
𝐴(𝑡) = 200 + 𝐶𝑒 −50𝑡

𝐶𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑠𝑢𝑠𝑡𝑎𝑛𝑐𝑖𝑎
10 La fórmula de concentración de una sustancia es 𝐶 =
𝑉𝑜𝑙𝑢𝑚𝑒𝑛
17
Aplicando la condición inicial 𝐴(0) = 30:
1
𝐴(0) = 200 + 𝐶𝑒 −50(0) = 200 + 𝐶 = 30
Por lo tanto
𝐶 = 30 − 200 = −170
La cantidad de gramos de sal que hay en el tanque al tiempo 𝑡 es
1
𝐴(𝑡) = 200 − 170𝑒 −50𝑡 ∎

Ejercicio en la guía 11: Resuelva el ejercicio del ejemplo 11 suponiendo que al tanque entra agua pura.

Diferencia entre las velocidades de entrada y salida

 Cuando la solución bien mezclada sale con una razón mayor entonces el tanque desalojará la
𝑔𝑎𝑙
solución con una razón de (𝑟𝑒𝑛𝑡𝑟𝑎 − 𝑟𝑠𝑎𝑙𝑒 ) 𝑚𝑖𝑛 y después de 𝑡 minutos (𝑟𝑒𝑛𝑡𝑟𝑎 − 𝑟𝑠𝑎𝑙𝑒 )𝑡 galones se
desalojarán.
 Cuando la solución bien mezclada sale con una razón menor entonces se acumulará liquido en el
𝑔𝑎𝑙
tanque con una razón de (𝑟𝑒𝑛𝑡𝑟𝑎 − 𝑟𝑠𝑎𝑙𝑒 ) 𝑚𝑖𝑛 y después de 𝑡 minutos (𝑟𝑒𝑛𝑡𝑟𝑎 − 𝑟𝑠𝑎𝑙𝑒 )𝑡 galones se
acumularán.
Por lo tanto en el tanque habrá:

𝑉(𝑡) = 𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒𝑙 𝑡𝑎𝑛𝑞𝑢𝑒 + (𝑟𝑒𝑛𝑡𝑟𝑎 − 𝑟𝑠𝑎𝑙𝑒 )𝑡

Ejemplo 13: Un gran tanque de 500 galones está lleno de agua pura. Le entra salmuera que tiene 2 lb
de sal por galón a razón de 5 gal/min. La solución bien mezclada sale del tanque con una razón de 10
gal/min. Determine la cantidad 𝐴(𝑡) de libras de sal que hay en el tanque al tiempo 𝑡.

Solución: Inicialmente el tanque está lleno de agua pura, por lo tanto la condición inicial es 𝐴(0) = 0. La
razón de entrada de sal es
𝑅𝑒𝑛𝑡𝑟𝑎 = (𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑠𝑎𝑙)(𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎)
𝑅𝑒𝑛𝑡𝑟𝑎 = (2𝑙𝑏/𝑔𝑎𝑙)(5 𝑔𝑎𝑙/𝑚𝑖𝑛) = 10 𝑙𝑏/𝑚𝑖𝑛
18
Como la solución bien mezclada sale con una razón mayor entonces el tanque desalojará la solución con
𝑔𝑎𝑙
una razón de 𝑟𝑒𝑛𝑡𝑟𝑎 − 𝑟𝑠𝑎𝑙𝑒 = (5 − 10) = −5 𝑔𝑎𝑙/𝑚𝑖𝑛. Después de 𝑡 minutos (−5 𝑔𝑎𝑙/𝑚𝑖𝑛)∙(𝑡 𝑚𝑖𝑛)=
𝑚𝑖𝑛
−5𝑡 𝑔𝑎𝑙 se desalojaran, por lo que en el tanque habrá 500 − 5𝑡 galones de salmuera.
La razón de salida de sal es
𝑅𝑠𝑎𝑙𝑒 = (𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑠𝑎𝑙)(𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎)
𝐴(𝑡) 10𝐴(𝑡) 2𝐴(𝑡)
𝑅𝑠𝑎𝑙𝑒 = ( 𝑙𝑏/𝑔𝑎𝑙) (10 𝑔𝑎𝑙/𝑚𝑖𝑛) = 𝑙𝑏/𝑚𝑖𝑛 = 𝑙𝑏/𝑚𝑖𝑛
500 − 5𝑡 500 − 5𝑡 100 − 𝑡
Por lo tanto, el problema de valor inicial de mezclas es
𝑑𝐴 2𝐴
= 10 − ; 𝐴(0) = 0
𝑑𝑡 100 − 𝑡
La ecuación diferencial resultante es lineal y su forma estándar es
𝑑𝐴 2𝐴
+ = 10
𝑑𝑡 100 − 𝑡
2
Donde 𝑝(𝑡) = 100−𝑡 y 𝑞(𝑡) = 10, la solución de la ecuación lineal es

∫ 𝑒 ∫ 𝑝(𝑡)𝑑𝑡 ∙ 𝑞(𝑡)𝑑𝑡 + 𝐶
𝐴(𝑡) =
𝑒 ∫ 𝑝(𝑡)𝑑𝑡

2
𝑑𝑡
∫ 𝑒 ∫100−𝑡 ∙ 10 𝑑𝑡 + 𝐶 10 ∫ 𝑒 −2 ln(100−𝑡) 𝑑𝑡 + 𝐶 10 ∫ 𝑒 −2 ln(100−𝑡) 𝑑𝑡 + 𝐶 10 ∫(100 − 𝑡)−2 𝑑𝑡 + 𝐶
𝐴(𝑡) = 2 = = =
𝑒 −2 ln(100−𝑡) 𝑒 ln(100−𝑡)
−2
(100 − 𝑡)−2
𝑒 ∫100−𝑡𝑑𝑡
10(100 − 𝑡)−1 + 𝐶
= = 10(100 − 𝑡) + 𝐶(100 − 𝑡)2
(100 − 𝑡)−2
Tenemos que
𝐴(𝑡) = 1000 − 10𝑡 + 𝐶(100 − 𝑡)2
Aplicando la condición inicial 𝐴(0) = 0:
𝐴(0) = 1000 − 10(0) + 𝐶(100 − (0))2 = 1000 + 𝐶(10000) = 0
Por lo tanto
−1000 1
𝐶= =−
10000 10

La cantidad de gramos de sal que hay en el tanque al tiempo 𝑡 es

1
𝐴(𝑡) = 1000 − 10𝑡 − (100 − 𝑡)2 ∎
10

19
Ejercicio en la guía 12: Un gran tanque está parcialmente lleno con 100 galones de fluido en los que se
disolvieron 10 libras de sal. La salmuera tiene ½ libra de sal por galón que entra al tanque a razón de 6
gal/min. La solución bien mezclada sale del tanque a razón de 4 gal/min. Determine la cantidad de libras
de sal que hay en el tanque después de 30 minutos.

2.3. Soluciones por sustitución

En esta sección se trabajaran otros métodos para resolver ecuaciones diferenciales de primer orden
mediante diferentes tipos de sustituciones o transformaciones.

2.3.1. Ecuaciones homogéneas de primer orden

Antes de resolver las ecuaciones diferenciales, debemos verificar ¿cuándo una función es homogénea?.
Una función 𝑓(𝑥, 𝑦) es homogénea de grado 𝛼 si
𝑓(𝑡𝑥, 𝑡𝑦) = 𝑡 ∝ 𝑓(𝑥, 𝑦)

Ejemplo 14: Verificar que la función 𝑓(𝑥, 𝑦) = 𝑥 2 + 𝑦 2 es homogénea de grado ∝= 2.

Solución: Verifiquemos la igualdad
𝑓(𝑡𝑥, 𝑡𝑦) = (𝑡𝑥)2 + (𝑡𝑦)2 = 𝑡 2 𝑥 2 + 𝑡 2 𝑦 2 = 𝑡 2 (𝑥 2 + 𝑦 2 ) = 𝑡 ∝ (𝑥 2 + 𝑦 2 )
Por lo tanto, la función es homogénea de grado ∝= 2. ∎
Una ecuación diferencial de la forma
𝑀(𝑥, 𝑦)𝑑𝑥 + 𝑁(𝑥, 𝑦)𝑑𝑦 = 0

20
Se puede ver como una ecuación diferencial homogénea 11 si, tanto 𝑀(𝑥, 𝑦) como 𝑁(𝑥, 𝑦) son
homogéneas del mismo grado ∝.

Ejemplo 15: Verifiquemos que la siguiente ecuación no es homogénea, donde 𝑅, 𝐿 y 𝑉 son constantes.
(𝑅𝑖 − 𝑉)𝑑𝑡 + 𝐿𝑑𝑖 = 0
Solución:
Definimos 𝑀(𝑡, 𝑖) = 𝑅𝑖 − 𝑉
𝑀(𝑘𝑡, 𝑘𝑖) = 𝑅(𝑘𝑖) − 𝑉
Pero es imposible factorizar el valor de 𝑘.∎

Ejercicio en la guía 13: Determine si la ecuación diferencial (𝑦 2 + 𝑥𝑦)𝑑𝑥 + 𝑥 2 𝑑𝑦 = 0 es homogénea y

calcule su grado de homogeneidad.

Solución de una ecuación diferencial homogénea:

La siguiente ecuación diferencial
(𝑥 − 𝑦)𝑑𝑥 + 𝑥𝑑𝑦 = 0
Es un caso particular de una ecuación diferencial homogénea, ya que 𝑀(𝑥, 𝑦) = 𝑥 − 𝑦 y 𝑁(𝑥, 𝑦) = 𝑥 son
funciones homogéneas de grado ∝= 1.

11 Las ecuaciones homogéneas son también llamadas ecuaciones diferenciales reducibles a separables.
21
Ejercicio en la guía 14: Verificar la afirmación anterior.

Para empezar a encontrar su solución, utilizamos la siguiente transformación

𝑦 = 𝑢𝑥
𝑑𝑦 = 𝑥𝑑𝑢 + 𝑢𝑑𝑥

Por tanto la ecuación diferencial queda de la siguiente forma

[𝑥 − (𝑢𝑥)]𝑑𝑥 + 𝑥(𝑥𝑑𝑢 + 𝑢𝑑𝑥) = 0

Aplicando propiedad distributiva y sumando términos semejantes

𝑥𝑑𝑥 − 𝑢𝑥𝑑𝑥 + 𝑥 2 𝑑𝑢 + 𝑢𝑥𝑑𝑥 = 0

𝑥𝑑𝑥 + 𝑥 2 𝑑𝑢 = 0

Observemos que mediante el uso de la transformación la ecuación homogénea queda convertida en una
ecuación diferencial de variables separables.

Resultado importante: Toda ecuación diferencial homogénea de primer orden siempre puede ser reducida a una
ecuación diferencial de variables separables.

Siguiendo con el ejercicio, si separamos variables e integrando ambos lados

𝑥
− 𝑑𝑥 = 𝑑𝑢
𝑥2
− ln|𝑥| + 𝑐 = 𝑢
𝑦
De la transformación inicial sabemos que 𝑢 = 𝑥 , por tanto
𝑦
= − ln|𝑥| + 𝑐
𝑥
Y la solución explicita de la ecuación es

𝑦 = −𝑥 ln|𝑥| + 𝑐𝑥

22
Ejercicio en la guía 15: Verificar que la solución explícita anterior satisface la ecuación diferencial
(𝑥 − 𝑦)𝑑𝑥 + 𝑥𝑑𝑦 = 0

Ejemplo 16: Resolver la ecuación diferencial (𝑥 + √𝑥𝑦)𝑑𝑦 = 𝑦𝑑𝑥.

Solución: Observemos que 𝑀(𝑥, 𝑦) = −𝑦 y 𝑁(𝑥, 𝑦) = 𝑥 + √𝑥𝑦. Verifiquemos que 𝑀 y 𝑁 son funciones
homogéneas de grado ∝= 1.
 𝑀(𝑡𝑥, 𝑡𝑦) = −𝑡𝑦 = 𝑡(−𝑦), función homogénea de grado ∝= 1.
 𝑁(𝑡𝑥, 𝑡𝑦) = 𝑡𝑥 + √(𝑡𝑥)(𝑡𝑦) = 𝑡𝑥 + √𝑡 2 𝑥𝑦 = 𝑡𝑥 + 𝑡√𝑥𝑦 = 𝑡(𝑥 + √𝑥𝑦), función homogénea de grado
∝= 1.

Si utilizamos la transformación indicada, la ecuación diferencial queda

(𝑥 + √𝑥(𝑢𝑥)) (𝑥𝑑𝑢 + 𝑢𝑑𝑥) = 𝑢𝑥𝑑𝑥

(𝑥 + 𝑥√𝑢)(𝑥𝑑𝑢 + 𝑢𝑑𝑥) = 𝑢𝑥𝑑𝑥

𝑥 2 𝑑𝑢 + 𝑢𝑥𝑑𝑥 + 𝑥 2 √𝑢 𝑑𝑢 + 𝑥𝑢3/2 𝑑𝑥 − 𝑢𝑥𝑑𝑥 = 0

𝑢𝑥𝑑𝑥 + 𝑥𝑢3/2 𝑑𝑥 − 𝑢𝑥𝑑𝑥 = −𝑥 2 𝑑𝑢 − 𝑥 2 √𝑢 𝑑𝑢

Simplificando el lado izquierdo y Factorizando el lado derecho

𝑥𝑢3/2 𝑑𝑥 = −𝑥 2 (1 + √𝑢)𝑑𝑢

Separando variables

𝑥 1 + √𝑢
− 𝑑𝑥 = 𝑑𝑢
𝑥2 𝑢3/2
Simplificando ambos lados e integrando
1
− 𝑑𝑥 = (𝑢−3/2 + 𝑢−1 )𝑑𝑢
𝑥
− 𝑙𝑛|𝑥| + 𝑐 = −2√𝑢 + 𝑙𝑛|𝑢|
23
 𝑦
Como 𝑢 = 𝑥 , entonces la solución implícita es

𝑦 𝑦
− 𝑙𝑛|𝑥| + 𝑐 = −2√ + 𝑙𝑛 | | ∎
𝑥 𝑥

Ejercicio en la guía 16: Verificar la solución del ejemplo 15.

𝑑𝑦
Ejercicio en la guía 17: Resolver el problema de valor inicial: 𝑥𝑦 2 𝑑𝑥 = 𝑦 3 − 𝑥 3 ; 𝑦(1) = 2.

2.3.2. Ecuación de Bernoulli

Una ecuación diferencial de la forma

𝑑𝑦
+ 𝑔(𝑥)𝑦 = 𝑞(𝑥)𝑦 𝑛
𝑑𝑥
Es llamada una ecuación diferencial de Bernoulli12. Si utilizamos la sustitución

12Jacob Bernoulli (1654-1705) científico suizo, en 1690 se convirtió en la primera persona en desarrollar la técnica de la
separación de variables.
24
 𝑢 = 𝑦 1−𝑛

Entonces
𝑑𝑢 𝑑𝑦
= (1 − 𝑛)𝑦 −𝑛
𝑑𝑥 𝑑𝑥
𝑑𝑦
Despejando 𝑑𝑥 ,

𝑑𝑦 𝑦 𝑛 𝑑𝑢
=
𝑑𝑥 1 − 𝑛 𝑑𝑥
Reemplazando en la ecuación diferencial de Bernoulli
𝑦 𝑛 𝑑𝑢 1 𝑦
+ 𝑔(𝑥)𝑢1−𝑛 = 𝑞(𝑥)
1 − 𝑛 𝑑𝑥 𝑢
𝑛
Es fácil verificar que 𝑦 𝑛 = 𝑢1−𝑛 y por lo tanto
𝑑𝑢
+ (1 − 𝑛)𝑔(𝑥)𝑢 = (1 − 𝑛)𝑞(𝑥)
𝑑𝑥
La ecuación diferencial anterior es una ecuación lineal y se resuelve mediante el método estudiado en la
sección 2.2.
Resultado importante: Toda ecuación diferencial de Bernoulli se puede transformar en una ecuación diferencial lineal
de primer orden.
Para finalizar, la ecuación diferencial lineal depende de la variable 𝑢, por lo tanto se debe utilizar de nuevo
la sustitución 𝑢 = 𝑦1−𝑛 para retornar a la variable dependiente inicial 𝑦.

Ejemplo 17: Convertir la ecuación diferencial

𝑑𝑦
𝑡2 + 𝑦 2 = 𝑡𝑦
𝑑𝑡
en una ecuación diferencial lineal.
Solución: Observemos que la ecuación diferencial es de Bernoulli
𝑑𝑦 1 1
− 𝑦 = − 2 𝑦2
𝑑𝑡 𝑡 𝑡
Con 𝑛 = 2 y utilizando la sustitución 𝑢 = 𝑦 1−𝑛 = 𝑦 −1 , tenemos que
𝑑𝑦 𝑑𝑢
= −𝑦 2
𝑑𝑡 𝑑𝑡
Finalmente

25
 𝑑𝑢 1 1
+ 𝑢= 2
𝑑𝑡 𝑡 𝑡
La anterior ecuación es una ecuación diferencial lineal de primer orden.∎

Ejercicio en la guía 18: Resolver la ecuación lineal resultante del ejemplo 17 y posteriormente encontrar la
solución de la ecuación diferencial de Bernoulli.

2.4. Ecuaciones exactas

Una ecuación diferencial de la forma

𝑀(𝑥, 𝑦)𝑑𝑥 + 𝑁(𝑥, 𝑦)𝑑𝑦 = 0
Se denomina exacta si se satisface la siguiente igualdad:
𝜕𝑀 𝜕𝑁
=
𝜕𝑦 𝜕𝑥
Ejemplo 18: Verificar si la ecuación diferencial
(𝑠𝑖𝑛 𝑦 − 𝑦 𝑠𝑖𝑛 𝑥)𝑑𝑥 + (𝑐𝑜𝑠 𝑥 + 𝑥 𝑐𝑜𝑠 𝑦 − 𝑦)𝑑𝑦 = 0
es exacta.
Solución: Sea 𝑀(𝑥, 𝑦) = 𝑠𝑖𝑛 𝑦 − 𝑦 𝑠𝑖𝑛 𝑥 y 𝑁(𝑥, 𝑦) = 𝑐𝑜𝑠 𝑥 + 𝑥 𝑐𝑜𝑠 𝑦 − 𝑦
𝜕𝑀
Encontremos 𝜕𝑦

𝜕𝑀
= 𝑐𝑜𝑠 𝑦 − 𝑠𝑖𝑛 𝑥
𝜕𝑦

26
 𝜕𝑁
Ahora encontremos 𝜕𝑥

𝜕𝑁
= − 𝑠𝑖𝑛 𝑥 + 𝑐𝑜𝑠 𝑦 = 𝑐𝑜𝑠 𝑦 − 𝑠𝑖𝑛 𝑥
𝜕𝑥
𝜕𝑀 𝜕𝑁
Claramente se verifica que = ∎
𝜕𝑦 𝜕𝑥

𝑑𝑦
Ejercicio en la guía 19: Verificar que la ecuación diferencial (1 − 2𝑥 2 − 2𝑦) 𝑑𝑥 = 4𝑥 3 + 4𝑥𝑦 es exacta.

Ejemplo 19: Determinar el valor de 𝑘 de manera que la ecuación diferencial sea exacta

(6𝑥𝑦 3 + 𝑐𝑜𝑠𝑦)𝑑𝑥 + (2𝑘𝑥 2 𝑦 2 − 𝑥𝑠𝑒𝑛𝑦)𝑑𝑦 = 0

Solución: Sea 𝑀(𝑥, 𝑦) = 6𝑥𝑦 3 + 𝑐𝑜𝑠𝑦 y 𝑁(𝑥, 𝑦) = 2𝑘𝑥 2 𝑦 2 − 𝑥𝑠𝑒𝑛𝑦, basta igualar las derivadas parciales
y despejar 𝑘
18𝑥𝑦 2 − 𝑠𝑒𝑛𝑦 = ⏟
⏟ 4𝑘𝑥𝑦 2 − 𝑠𝑒𝑛𝑦
𝑀𝑦 𝑁𝑥

9
De aquí, 18 = 4𝑘 y por tanto 𝑘 = 2 . ∎

Solución de ecuaciones diferenciales exactas

Dada una ecuación diferencial de la forma
𝑀(𝑥, 𝑦)𝑑𝑥 + 𝑁(𝑥, 𝑦)𝑑𝑦 = 0
Existe una función 𝐹 tal que
𝜕𝐹
= 𝑀(𝑥, 𝑦)
𝜕𝑥
Podemos determinar la función 𝐹 integrando 𝑀 con respecto a 𝑥 mientras la variable 𝑦 se conserva
constante

27
 𝐹(𝑥, 𝑦) = ∫ 𝑀(𝑥, 𝑦)𝑑𝑥 + 𝑔(𝑦)

donde 𝑔(𝑦) es la constante de integración. Ahora derivando con respecto a 𝑦:

𝜕𝐹 𝜕
= ∫ 𝑀(𝑥, 𝑦)𝑑𝑥 + 𝑔′ (𝑦)
𝜕𝑦 𝜕𝑦
𝜕𝐹
Como 𝜕𝑦 = 𝑁(𝑥, 𝑦), entonces

𝜕
∫ 𝑀(𝑥, 𝑦)𝑑𝑥 + 𝑔′ (𝑦) = 𝑁(𝑥, 𝑦)
𝜕𝑦
De donde
𝜕
𝑔′ (𝑦) = 𝑁(𝑥, 𝑦) − ∫ 𝑀(𝑥, 𝑦)𝑑𝑥
𝜕𝑦
integrando 𝑔′ (𝑦) :
𝜕
𝑔(𝑦) = ∫ (𝑁(𝑥, 𝑦) − ∫ 𝑀(𝑥, 𝑦)𝑑𝑥) 𝑑𝑦
𝜕𝑦
reemplazando 𝑔(𝑦) en 𝐹(𝑥, 𝑦)
𝜕
𝐹(𝑥, 𝑦) = ∫ 𝑀(𝑥, 𝑦)𝑑𝑥 + ∫ (𝑁(𝑥, 𝑦) − ∫ 𝑀(𝑥, 𝑦)𝑑𝑥) 𝑑𝑦
𝜕𝑦
La solución implícita es 𝐹(𝑥, 𝑦) = 𝐾, entonces
𝜕
∫ 𝑀(𝑥, 𝑦)𝑑𝑥 + ∫ (𝑁(𝑥, 𝑦) − ∫ 𝑀(𝑥, 𝑦)𝑑𝑥) 𝑑𝑦 = 𝐾
𝜕𝑦

Ejemplo 20: Resolver el problema de valor inicial 𝑥𝑑𝑥 + (𝑥 2 𝑦 + 4𝑦)𝑑𝑦 = 0, 𝑦(4) = 0.

Solución: Sea 𝑀(𝑥, 𝑦) = 𝑥 y 𝑁(𝑥, 𝑦) = 𝑥 2 𝑦 + 4𝑦, entonces 𝑀𝑦 = 0 y 𝑁𝑥 = 2𝑥𝑦, por tanto la ecuación
2
diferencial no es exacta. Si utilizamos el factor integrante 𝜇(𝑦) = 𝑒 𝑦 y multiplicamos la ecuación inicial
2 2
𝑥𝑒 𝑦 𝑑𝑥 + 𝑒 𝑦 (𝑥 2 𝑦 + 4𝑦)𝑑𝑦 = 0
Se obtiene una ecuación diferencial exacta.

28
 2
Ejercicio en la guía 20: Verificar que el factor integrante 𝜇(𝑦) = 𝑒 𝑦 convierte la ecuación diferencial del
ejemplo 20 en exacta.

Siguiendo con el ejemplo 20:

𝜕𝐹 2 𝜕𝐹 2 𝜕𝐹
Sea 𝜕𝑥 = 𝑥𝑒 𝑦 y = 𝑒 𝑦 (𝑥 2 𝑦 + 4𝑦), integrando 𝜕𝑥 con respecto a 𝑥, tenemos
𝜕𝑦

𝜕𝐹 2 1 2
∫ 𝑑𝑥 = 𝐹(𝑥, 𝑦) = ∫ 𝑥𝑒 𝑦 𝑑𝑥 = 𝑥 2 𝑒 𝑦 + 𝑔(𝑦)
𝜕𝑥 2
Derivando 𝐹 con respecto a 𝑦
𝜕𝐹 2
= 𝑦𝑥 2 𝑒 𝑦 + 𝑔′(𝑦)
𝜕𝑦
𝜕𝐹
Sabiendo que 𝜕𝑦 = 𝑁(𝑥, 𝑦), tenemos
2 2
𝑦𝑥 2 𝑒 𝑦 + 𝑔′ (𝑦) = 𝑒⏟𝑦 (𝑥 2 𝑦 + 4𝑦)
𝑁(𝑥,𝑦)

2
Por lo tanto 𝑔′ (𝑦) = 4𝑦𝑒 𝑦 . Para encontrar 𝑔(𝑦)

2 2
𝑔(𝑦) = ∫ 𝑔′ (𝑦)𝑑𝑦 = 4𝑦𝑒 𝑦 𝑑𝑦 = 2𝑒 𝑦
La solución implícita es
1 2 2
𝐹(𝑥, 𝑦) = 𝑥 2 𝑒 𝑦 + 2𝑒 𝑦 = 𝐾
2
Si 𝑦(4) = 0, entonces
1 2 2
(4)2 𝑒 (0) + 2𝑒 (0) = 10 = 𝐾
𝐹(4,0) =
2
La solución implícita del problema de valor inicial es
1 2 𝑦2 2
𝑥 𝑒 + 2𝑒 𝑦 = 10 ∎
2

29
Ejercicio en la guía 21: Verificar que la solución implícita del ejemplo 19 satisface la ecuación diferencial.

2𝑥 𝑥2
Ejercicio en la guía 22: Resolver la ecuación diferencial 𝑦
𝑑𝑥 − 𝑦 2 𝑑𝑦 = 0

30
 CAPÍTULO 3
3. Ecuaciones diferenciales de orden
superior
Diversos fenómenos físicos se modelan mediante ecuaciones diferenciales de orden superior (𝑛 ≥ 2). En
este capítulo se mostrará la solución de ecuaciones diferenciales de orden superior y el fenómeno masa
resorte modelado mediante una ecuación diferencial de orden dos.

3.1. Ecuaciones diferenciales lineales homogéneas

Una ecuación diferencial lineal homogénea de orden 𝑛 es de la forma

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥) 𝑛
+ 𝑎 𝑛−1 (𝑥) 𝑛−1
+ ⋯ + 𝑎1 (𝑥) + 𝑎0 (𝑥)𝑦 = 𝟎
𝑑𝑥 𝑑𝑥 𝑑𝑥
La característica principal es su igualación a cero.

En este capítulo se trabajaran las ecuaciones lineales de orden superior cuyos coeficientes (A, B, C y D)
de las derivadas y la variable dependiente son constantes, esto es:
𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝐴 + 𝐵 + ⋯ + 𝐶 + 𝐷𝑦 = 𝟎
𝑑𝑥 𝑛 𝑑𝑥 𝑛−1 𝑑𝑥
La ecuación auxiliar de la ecuación diferencial se define de la forma
𝐴𝑚𝑛 + 𝐵𝑚𝑛−1 + ⋯ + 𝐶𝑚 + 𝐷 = 0

Finalmente, dada una ecuación diferencial homogénea con coeficientes constantes de orden superior, la
solución depende del número de soluciones de la ecuación auxiliar.
Tomemos un caso particular de una ecuación auxiliar de grado dos
 Si sus soluciones son números reales diferentes 𝑚1 ≠ 𝑚2 , entonces la solución de la ecuación
diferencial es
𝑦(𝑥) = 𝑐1 𝑒 𝑚1 𝑥 + 𝑐2 𝑒 𝑚2 𝑥
 Si sus soluciones son números reales iguales 𝑚1 = 𝑚2 , entonces la solución de la ecuación
diferencial es
𝑦(𝑥) = 𝑐1 𝑒 𝑚1 𝑥 + 𝑐2 𝑥𝑒 𝑚2 𝑥

31
  Si sus soluciones son números complejos 𝑚 =∝ ±𝛽𝑖, entonces la solución de la ecuación
diferencial es
𝑦(𝑥) = 𝑒 ∝𝑥 (𝑐1 𝑐𝑜𝑠𝛽𝑥 + 𝑐2 𝑠𝑒𝑛𝛽𝑥)

Ejemplo 21: Resuelva el problema con valores iniciales

𝑦 ′′ + 2𝑦 ′ + 2𝑦 = 0; 𝑦(0) = 2, 𝑦 ′ (0) = 1
Solución: La ecuación auxiliar asociada a la ecuación diferencial es
𝑚2 + 2𝑚 + 2 = 0
Utilizando la formula cuadrática para encontrar la solución

−2 ± √4 − 4(1)(2)
𝑚=
2(1)
Donde 𝑚 = −1 ± 𝑖, de aquí ∝= −1 y 𝛽 = 1. Reemplazando en la fórmula para solución compleja
𝑦(𝑥) = 𝑒 −𝑥 (𝑐1 𝑐𝑜𝑠𝑥 + 𝑐2 𝑠𝑒𝑛𝑥)
Utilizando la primera condición inicial

𝑦(0) = 𝑒 −(0) (𝑐1 𝑐𝑜𝑠(0) + 𝑐2 𝑠𝑒𝑛(0)) = 2

De aquí 𝑐1 = 2. Derivando la solución general
𝑦 ′ (𝑥) = −𝑒 −𝑥 (𝑐1 𝑐𝑜𝑠𝑥 + 𝑐2 𝑠𝑒𝑛𝑥) + 𝑒 −𝑥 (−𝑐1 𝑠𝑒𝑛𝑥 + 𝑐2 𝑐𝑜𝑠𝑥)
Utilizando la segunda condición inicial

𝑦 ′ (0) = −𝑒 −(0) (𝑐1 𝑐𝑜𝑠(0) + 𝑐2 𝑠𝑒𝑛(0)) + 𝑒 −(0) (−𝑐1 𝑠𝑒𝑛(0) + 𝑐2 𝑐𝑜𝑠(0)) = 1

De aquí, 𝑐2 − 𝑐1 = 1, si 𝑐1 = 2 entonces 𝑐2 = 3. Por lo tanto

𝑦(𝑥) = 𝑒 −𝑥 (2𝑐𝑜𝑠𝑥 + 3𝑠𝑒𝑛𝑥) ∎

Ejercicio en la guía 23: Resolver la ecuación diferencial 2𝑦 ′′ − 3𝑦 ′ + 4𝑦 = 0

𝟑
√𝟐𝟑 √𝟐𝟑
R/ 𝒚(𝒙) = 𝒆𝟒𝒙 (𝒄𝟏 𝐜𝐨𝐬 𝒙 + 𝒄𝟐 𝐬𝐢𝐧 𝒙) 32
𝟒 𝟒
Ejemplo 22: Encuentre la solución general de la ecuación diferencial de tercer orden
𝑑3𝑢 𝑑2 𝑢
+ − 2𝑢 = 0
𝑑𝑡 3 𝑑𝑡 2
Solución: La ecuación auxiliar asociada es
𝑚3 + 𝑚2 − 2 = 0
Es claro que 𝑚 = 1 es una solución de la ecuación. Para encontrar las otras soluciones, utilizamos división
sintética
1 1 0 −2 𝟏
1 2 2
1 2 2 0

Por lo tanto,
𝑚3 + 𝑚2 − 2 = (𝑚 − 1)(𝑚2 + 2𝑚 + 2) = 0
De aquí, 𝑚 − 1 = 0 ó 𝑚2 + 2𝑚 + 2 = 0. Utilizando la formula cuadrática en esta última ecuación

−2 ± √4 − 4(1)(2) −2 ± √−4 −2 √−4 2𝑖

𝑚= = = ± = −1 ± = −1 ± 𝑖
2(1) 2 2 2 2
De aquí, ∝= −1 y 𝛽 = 1.
Observemos que tenemos una solución real y dos soluciones complejas, por lo tanto la solución de la
ecuación diferencial del tercer orden es
𝑢(𝑡) = 𝑐1 𝑒 𝑡 + 𝑒 −𝑡 (𝑐2 𝑐𝑜𝑠𝑡 + 𝑐3 𝑠𝑒𝑛𝑡)∎

Ejercicio en la guía 24: Resolver la ecuación diferencial 𝑦 ′′ − 36𝑦 = 0

33
Ejercicio en la guía 25: Resolver la ecuación diferencial 𝑦 ′′ + 6𝑦 ′ + 9𝑦 = 0

R/ 𝒚(𝒙) = 𝒄𝟏 𝒆−𝟑𝒙 + 𝒄𝟐 𝒙𝒆−𝟑𝒙

3.1.1. Sistema Masa-Resorte, movimiento libre no amortiguado

El movimiento libre no amortiguado de un sistema masa-resorte es un caso particular de una ecuación

diferencial lineal homogénea de orden dos.
La ecuación diferencial del movimiento libre no amortiguado es
𝑑2𝑥
+ 𝜔2 𝑥 = 0
𝑑𝑡 2
𝑘
donde 𝜔2 = 𝑚 .
𝑊
𝑚 = 𝑔 : Masa del objeto sujeto al resorte. 𝑊 es el peso de la masa y 𝑔 es la gravedad13.

𝑘: Constante de resorte, que por la ley de Hooke14, el resorte mismo ejerce una fuerza restauradora 𝐹
opuesta a la dirección de elongación y proporcional a la cantidad de elongación 𝑠 y es expresada como
𝐹 = 𝑘𝑠.

Dos condiciones iniciales relacionadas con la ecuación son el desplazamiento inicial 𝑥(0) = 𝑥0 y la
velocidad inicial 𝑥 ′ (0) = 𝑥1 . Por ejemplo,
 si 𝑥0 > 0 y 𝑥1 < 0, la masa parte de un punto debajo de la posición de equilibrio con una velocidad
aplicada hacia arriba.
 Si 𝑥0 = 0, se dice que la masa se libera del reposo.
 si 𝑥0 < 0 y 𝑥1 = 0, la masa se libera desde el reposo de un punto |𝑥0 | unidades arriba de la
posición de equilibrio.

𝑙𝑖𝑏𝑟𝑎𝑠
13 𝑚= = 𝑠𝑙𝑢𝑔
𝑝𝑖𝑒/𝑠 2
14Robert Hooke (1635-1703), científico inglés, formuló la teoría del movimiento planetario como un problema de mecánica, y
mantuvo continuas disputas con su contemporáneo Isaac Newton respecto a la teoría de la luz y a la ley de la gravitación
universal.
34
Ejemplo 23: Una masa que pesa 24 libras, unida al extremo de un resorte, lo alarga 4 pulgadas15. Al
inicio, la masa se libera desde el reposo en un punto 2 pulgadas arriba de la posición de equilibrio.
Encuentre la ecuación del movimiento.
𝑊 24 3
Solución: Calculando la masa: 𝑚 = = 32 = 4 slug. La constante de resorte con la ley de Hooke y
𝑔
1
sabiendo que 4 pulgadas es equivalente a 3 pie,

1
𝐹 = 𝑘𝑠 → 2 = 𝑘 ( )
3
por lo tanto, 𝑘 = 6. Nuestra ecuación diferencial queda
𝑑2 𝑥
+ 8𝑥 = 0
𝑑𝑡 2
El desplazamiento inicial y la velocidad inicial son 𝑥(0) = −2, 𝑥 ′ (0) = 0.
Resolviendo la ecuación auxiliar de segundo grado
𝑚2 + 8 = 0

Sus soluciones son 𝑚 = ±√8𝑖 con 𝛼 = 0 y 𝛽 = √8, por lo tanto

𝑥(𝑡) = 𝑐1 𝑐𝑜𝑠√8𝑡 + 𝑐2 𝑠𝑒𝑛√8𝑡

Aplicando la condición inicial 𝑥(0) = −2 a 𝑥(𝑡):

𝑥(0) = 𝑐1 𝑐𝑜𝑠√8(0) + 𝑐2 𝑠𝑒𝑛√8(0) = −2

Obtenemos que 𝑐1 = −2

𝑥(𝑡) = −2𝑐𝑜𝑠√8𝑡 + 𝑐2 𝑠𝑒𝑛√8𝑡

Derivando 𝑥(𝑡):

𝑥 ′ (𝑡) = 2√8𝑠𝑒𝑛√8𝑡 + √8𝑐2 𝑐𝑜𝑠√8𝑡

Aplicando la condición inicial 𝑥′(0) = 0 a 𝑥′(𝑡):

𝑥 ′ (0) = 2√8𝑠𝑒𝑛√8(0) + √8𝑐2 𝑐𝑜𝑠√8(0) = 0

Obtenemos que 𝑐2 = 0, por lo tanto la ecuación del movimiento queda

𝑥(𝑡) = −2𝑐𝑜𝑠√8𝑡 ∎
Ejemplo 24: Una masa que pesa 64 libras alarga 0.32 pies un resorte. Al inicio la masa se libera desde
un punto que está 8 pulgadas arriba de la posición de equilibrio con una velocidad descendente de 5
pies/s.
a) Encuentre la ecuación del movimiento.
b) ¿Cuál es la amplitud y el periodo del movimiento?

15 1 pie ≈ 12 pulgadas
35
 c) ¿Cuántos ciclos completos habrá realizado la masa al final de 3𝜋 segundos?
Solución:
a) Para encontrar las ecuación del movimiento determinamos la masa y la constante del resorte:

𝑤 64
𝑤 = 𝑚𝑔 ⇒ 𝑚= = = 2 𝑠𝑙𝑢𝑔
𝑔 32
𝐹 64
𝐹 = 𝑘𝑥 ⇒ 𝑘= = = 200
𝑥 0.32
La ecuación diferencial del movimiento armónico simple es:
2𝑥 ′′ + 200𝑥 = 0 → 𝑥 ′′ + 100𝑥 = 0
2
Con condiciones iniciales 𝑥(0) = −8 𝑝𝑙𝑔 ≅ − 3 𝑝𝑖𝑒 𝑦 𝑥 ′ (0) = 5.

La ecuación de grado dos auxiliar es 𝑚2 + 100 = 0, y sus soluciones son 𝑚 = ±10𝑖. La solución general
es
𝑥(𝑡) = 𝑐1 𝑐𝑜𝑠 10𝑡 + 𝑐2 𝑠𝑖𝑛 10𝑡
reemplazando la condición inicial de posición
2
𝑥(0) = 𝑐1 𝑐𝑜𝑠 10(0) + 𝑐2 𝑠𝑖𝑛 10(0) = −
3
Deducimos que 𝑐1 = −2/3. Si derivamos la solución general y reemplazamos la segunda condición de
velocidad
𝑥 ′ (𝑡) = −10𝑐1 𝑠𝑖𝑛 10𝑡 + 10𝑐2 𝑐𝑜𝑠 10𝑡
𝑥 ′ (0) = −10𝑐1 𝑠𝑖𝑛 10(0) + 10𝑐2 𝑐𝑜𝑠 10(0) = 5
Deducimos que 𝑐2 = 1/2. La solución del problema de valor inicial es
2 1
𝑥(𝑡) = − 𝑐𝑜𝑠 10𝑡 + 𝑠𝑖𝑛 10𝑡 ∎
3 2
Forma alternativa de la solución: Escribamos la solución anterior de la forma alternativa 𝑥(𝑡) =
𝐴 𝑠𝑖𝑛(𝜔𝑡 + 𝜙). Para esto usamos la ecuación

𝐴 = √𝑐12 + 𝑐22

El ángulo de fase 𝜙 se define por

𝑐
𝑠𝑖𝑛 𝜙 = 𝐴1 𝑐1
𝑐2 } 𝑡𝑎𝑛 𝜙 = 𝑐
𝑐𝑜𝑠 𝜙 = 𝐴 2

Por lo tanto

36
 2 2 1 2 4 1 5
√
𝐴 = (− ) + ( ) = √ + =
3 2 9 4 6

4
𝑡𝑎𝑛 𝜙 = −4/3 → 𝜙 = 𝑡𝑎𝑛−1 (− ) = −0.9272 𝑟𝑎𝑑
3
El ángulo fase es válido ya que 𝜙 es un ángulo en el cuarto cuadrante donde 𝑠𝑖𝑛 𝜙 < 0 y 𝑐𝑜𝑠 𝜙 > 0. La
ecuación en su forma alternativa es
5
𝑥(𝑡) = 𝑠𝑖𝑛(10𝑡 − 0.9272) ∎
6
2𝜋 2𝜋 𝜋
b) La amplitud es 𝐴 = 5/6 y el periodo se define como 𝑇 = |𝜔| = = 5.
10
1 𝜔
c) La frecuencia del movimiento es 𝑓 = 𝑇 = 2𝜋 , y es el número de ciclos completados cada segundo.
10
Por lo tanto, si en un segundo realiza 𝑓 = 2𝜋 = 5/𝜋 ciclos, en 𝑡 = 3𝜋 segundos realiza

5
(3𝜋) = 15 𝐶𝑖𝑐𝑙𝑜𝑠 ∎
𝜋
Ejercicio en la guía 26: Una fuerza de 400 newtons alarga 2 metros un resorte. Una masa de 50 kilogramos
se une al extremo del resorte y se libera inicialmente desde la posición de equilibrio con una velocidad
ascendente de 10 m/s. encuentre la ecuación del movimiento.

3.1.2. Sistema Masa-Resorte, movimiento libre amortiguado

La ecuación diferencial del movimiento libre amortiguado es

𝑑2𝑥 𝑑𝑥
2
+ 2𝜆 + 𝜔2 𝑥 = 0
𝑑𝑡 𝑑𝑡
𝑘 𝜑
donde 𝜔2 = 𝑚 y 2𝜆 = 𝑚.

37
𝑘 es la constante de resorte, 𝑚 es la masa del objeto sujeto al resorte y 𝜑 es una constante de
amortiguamiento positiva.
La ecuación auxiliar es 𝑚2 + 2𝜆𝑥 + 𝜔2 = 0 y sus soluciones son

𝑚 = −𝜆 ± √𝜆2 − 𝜔 2
 Caso 1: 𝜆2 > ω2 , el sistema está sobreamortiguado porque el coeficiente de amortiguamiento 𝜑
es grande comparado con la constante de resorte 𝑘. Su solución es16

2 −𝜔2 𝑡 2 −𝜔2 𝑡
𝑥(𝑡) = 𝑒 −𝜆𝑡 (𝑐1 𝑒 √𝜆 + 𝑐2 𝑒 −√𝜆 )

 Caso 2: 𝜆2 = ω2 , el sistema está críticamente amortiguado porque cualquier disminución en la

fuerza de amortiguamiento daría un movimiento oscilatorio. Su solución es

𝑥(𝑡) = 𝑒 −𝜆𝑡 (𝑐1 + 𝑐2 𝑡)

 Caso 3: 𝜆2 < ω2 , el sistema está subamortiguado porque el coeficiente de amortiguamiento 𝜑 es

pequeño comparado con la constante de resorte 𝑘. Su solución es

𝑥(𝑡) = 𝑒 −𝜆𝑡 (𝑐1 𝑐𝑜𝑠√𝜔 2 − 𝜆2 𝑡 + 𝑐2 𝑠𝑒𝑛√𝜔 2 − 𝜆2 𝑡)

Ejemplo 25: Una masa de 2 kiligramos se fija a un resorte cuya constante es 16 N/m y luego el sistema
completo se sumerge en una sustancia que imparte una fuerza amortiguada igual a 20 veces la velocidad
instantánea. Encuentre la ecuación del movimiento si la masa se libera inicialmente desde un punto 2
metros debajo de la posición de equilibrio con una velocidad ascendente de 14 m/s.
Solución: La ecuación diferencial es
𝑑2 𝑥 𝑑𝑥
2 2 + 20 + 16𝑥 = 0
𝑑𝑡 𝑑𝑡
El desplazamiento inicial y la velocidad inicial son 𝑥(0) = 2, 𝑥 ′ (0) = −14.
Resolviendo la ecuación auxiliar de segundo grado
2𝑚2 + 20𝑚 + 16 = 0

Sus soluciones son 𝑚 = −5 ± √17, por lo tanto

𝑥(𝑡) = 𝑐1 𝑒 (−5+√17)𝑡 + 𝑐2 𝑒 (−5−√17)𝑡

Aplicando la condición inicial 𝑥(0) = 2 a 𝑥(𝑡):

𝑥(0) = 𝑐1 𝑒 (−5+√17)(0) + 𝑐2 𝑒 (−5−√17)(0) = 2

16 Esta ecuación representa un movimiento uniforme y no oscilatorio.

38
Obtenemos que 𝑐1 + 𝑐2 = 2
Derivando 𝑥(𝑡):

𝑥 ′ (𝑡) = 𝑐1 (−5 + √17)𝑒 (−5+√17)𝑡 + 𝑐2 (−5 − √17)𝑒 (−5−√17)𝑡

Aplicando la condición inicial 𝑥 ′ (0) = −14 a 𝑥′(𝑡):

𝑥 ′ (0) = 𝑐1 (−5 + √17)𝑒 (−5+√17)(0) + 𝑐2 (−5 − √17)𝑒 (−5−√17)(0) = −14

Obtenemos que 𝑐1 (−5 + √17) + 𝑐2 (−5 − √17) = −14. Resolviendo el sistema de ecuaciones:

𝑐1 + 𝑐2 = 2
{
𝑐1 (−5 + √17) + 𝑐2 (−5 − √17) = −14

2 2
Obtenemos que 𝑐1 = 1 − y 𝑐2 = + 1. La ecuación del movimiento es
√17 √17

2 2
𝑥(𝑡) = (1 − ) 𝑒 (−5+√17)𝑡 + ( + 1) 𝑒 (−5−√17)𝑡 ∎
√17 √17
Ejercicio en la guía 27: Una fuerza de 2 libras alarga 1 pie un resorte. Una masa que pesa 3.4 libras se une
al resorte y luego se sumerge el sistema en un medio que ofrece una fuerza de amortiguamiento igual a
0,6 veces la velocidad instantánea.

a) Encuentre la ecuación del movimiento si inicialmente se libera la masa desde el reposo en un punto
situado a 2 pies por encima de la posición de equilibrio.
b) Encuentre la primera vez en que la masa pasa a través de la posición de equilibrio en dirección hacia
arriba.

39
 3.2. Ecuaciones diferenciales lineales no homogéneas

3.2.1. Método de los coeficientes indeterminados

Dada la ecuación diferencial lineal de orden superior

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥) 𝑛
+ 𝑎 𝑛−1 (𝑥) 𝑛−1
+ ⋯ + 𝑎1 (𝑥) + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥)
𝑑𝑥 𝑑𝑥 𝑑𝑥
Se denomina no homogénea si 𝑔(𝑥) ≠ 0.
Para encontrar la solución general de una ecuación diferencial lineal no homogénea debemos encontrar
dos soluciones llamadas 𝑦ℎ (solución de la ecuación homogénea) y 𝑦𝑝 (solución particular) y así, la
solución es
𝑦(𝑥) = 𝑦ℎ + 𝑦𝑝

Para encontrar 𝑦ℎ usamos lo estudiado en la sección 2.1.

Para encontrar la solución particular 𝑦𝑝 , debemos observar que tipo de función es 𝑔(𝑥) y así suponer la
forma de 𝑦𝑝 .

La siguiente tabla contiene diferentes formas de 𝑦𝑝 dependiendo de la forma de 𝑔(𝑥).

𝑔(𝑥) Ejemplo 𝑦𝑝
𝑔(𝑥) = 𝑘, 𝑘 constante 𝑔(𝑥) = 4 𝐴
𝑔(𝑥) es lineal 𝑔(𝑥) = 2𝑥 − 1 𝐴𝑥 + 𝐵
𝑔(𝑥) es senoidal o cosenoidal 𝑔(𝑥) = 2𝑠𝑒𝑛3𝑥 𝐴𝑠𝑒𝑛3𝑥 + 𝐵𝑐𝑜𝑠3𝑥
𝑔(𝑥) es exponencial 𝑔(𝑥) = 2𝑒 5𝑥 𝐴𝑒 5𝑥

Ejemplo 26: Resolver la ecuación diferencial lineal no homogénea de orden dos

𝑦 ′′ − 4𝑦 = 4𝑥
Solución: La solución general es de la forma
𝑦(𝑥) = 𝑦ℎ + 𝑦𝑝

Donde 𝑦ℎ es la solución de la ecuación diferencial homogénea

𝑦 ′′ − 4𝑦 = 0
La ecuación de grado dos auxiliar es 𝑚2 − 4 = 0, cuya solución es 𝑚 = ±2. La solución de la ecuación
homogénea es
𝑦ℎ = 𝑐1 𝑒 2𝑥 + 𝑐2 𝑒 −2𝑥
40
Para encontrar 𝑦𝑝 , debemos observar el tipo de función 𝑔(𝑥), dicha función es lineal y por lo tanto,
debemos suponer una solución particular de la forma 𝑦𝑝 = 𝐴𝑥 + 𝐵. Derivando dos veces

𝑦𝑝′ = 𝐴 , 𝑦𝑝′′ = 0

y reemplazando en la ecuación diferencial inicial no homogénea

0 − 4(𝐴𝑥 + 𝐵) = 4𝑥
donde 𝐴 = −1 y 𝐵 = 0. Por lo tanto 𝑦𝑝 = −𝑥.

La solución general es
𝑦(𝑥) = 𝑦ℎ + 𝑦𝑝 = 𝑐1 𝑒 2𝑥 + 𝑐2 𝑒 −2𝑥 − 𝑥 ∎

Error en el método: Consideremos la ecuación diferencial lineal no homogénea

𝑦 ′′ − 𝑦 ′ = 1
Para determinar 𝑦ℎ utilizamos la ecuación auxiliar 𝑚2 − 𝑚 = 0, obteniendo las soluciones 𝑚 = 0 y 𝑚 = 1,
por lo tanto
𝑦ℎ (𝑥) = 𝑐1 + 𝑐2 𝑒 𝑥
Para determinar 𝑦𝑝 , observemos que 𝑔(𝑥) = 1 y es constante, por lo tanto la solución particular tiene la
forma 𝑦𝑝 = 𝐴 de acuerdo a la tabla. Reemplazando en la ecuación diferencial con con 𝑦′𝑝 (𝑥) = 0 y
𝑦′′𝑝 (𝑥) = 0

⏟
0 − ⏟
0 =1
𝑦′′𝑝 (𝑥) 𝑦 ′ 𝑝 (𝑥)

Claramente 0 = 1 es una inconsistencia y no podemos determinar el coeficiente 𝐴, esto es llamado error

en el método.
Para solucionarlo, la supuesta solución particular inicial por un factor 𝑥 (según cual sea la variable
independiente de la ecuación planteada), quedando una nueva solución particular de la forma
𝑦𝑝 (𝑥) = 𝐴𝑥

y retomamos el proceso.
Comentario: Si la supuesta solución particular es igual a algún miembro de la suma en la solución homogénea, traerá
consigo un error en el método.

Ejemplo 27: Resolver 2𝑧 ′′ + 𝑧 = 9𝑒 2𝑡

Solución: Iniciamos resolviendo la ecuación diferencial homogénea asociada 2𝑧 ′′ + 𝑧 = 0, usando la
1
ecuación auxiliar 2𝑚2 + 1 = 0. Despejando 𝑚 obtenemos 𝑚 = ±√2 𝑖 y por lo tanto

41
 1 1
𝑧ℎ (𝑡) = 𝑐1 𝑐𝑜𝑠 (√ ) 𝑡+𝑐2 𝑠𝑒𝑛 (√ ) 𝑡
2 2

Para encontrar 𝑧𝑝 , debemos observar la forma de 𝑔(𝑡) = 9𝑒 2𝑡 y comparando en la tabla, suponemos la

solución de la forma
𝑧𝑝 (𝑡) = 𝐴𝑒 2𝑡

si 𝑧′𝑝 (𝑡) = 2𝐴𝑒 2𝑡 y 𝑧′′𝑝 (𝑡) = 4𝐴𝑒 2𝑡 y reemplazando en la ecuación diferencial no homogénea

2(4𝐴𝑒 2𝑡 ) + 𝐴𝑒 2𝑡 = 9𝑒 2𝑡
9𝐴𝑒 2𝑡 = 9𝑒 2𝑡
De aquí 𝐴 = 1 y por lo tanto 𝑧𝑝 (𝑡) = 𝑒 2𝑡 . La solución general es

1 1
𝑧(𝑡) = 𝑐1 𝑐𝑜𝑠 (√ ) 𝑡+𝑐2 𝑠𝑒𝑛 (√ ) 𝑡 + 𝑒 2𝑡 ∎
2 2

Ejemplo 28: Resolver 2𝑥 ′′ + 𝑥 ′ = 3𝑡 2

Solución: Resolviendo la ecuación diferencial homogénea asociada 2𝑥 ′′ + 𝑥 ′ = 0, obtenemos la solución

𝑥ℎ (𝑡) = 𝑐1 + 𝑐2 𝑒 (−1/2)𝑡
Para encontrar 𝑥𝑝 , suponemos la solución de la forma

𝑥𝑝 (𝑡) = 𝐴𝑡 2 + 𝐵𝑡 + 𝐶

Pero observemos que la parte constante de 𝑥ℎ (𝑡) es equivalente al termino 𝐶 de 𝑥𝑝 (𝑡) y esto trae un error
en el método.
Si multiplicamos la 𝑥𝑝 (𝑡) inicial por un factor 𝑡, obtenemos una nueva 𝑥𝑝 (𝑡) de la forma

𝑥𝑝 (𝑡) = 𝑡(𝐴𝑡 2 + 𝐵𝑡 + 𝐶) = 𝐴𝑡 3 + 𝐵𝑡 2 + 𝐶𝑡

donde 𝑥′𝑝 (𝑡) = 3𝐴𝑡 2 + 2𝐵𝑡 + 𝐶 y 𝑥′′𝑝 (𝑡) = 6𝐴𝑡 + 2𝐵. Reemplazando en la ecuación diferencial inicial

2(6𝐴𝑡 + 2𝐵) + 3𝐴𝑡 2 + 2𝐵𝑡 + 𝐶 = 3𝑡 2

de aquí, 𝐴 = 1, 𝐵 = −6, 𝐶 = 24. Por lo tanto
𝑥𝑝 (𝑡) = 𝑡 3 − 6𝑡 2 + 24𝑡

La solución general de la ecuación no homogénea es

𝑥(𝑡) = 𝑐1 + 𝑐2 𝑒 (−1/2)𝑡 + 𝑡 3 − 6𝑡 2 + 24𝑡 ∎

42
Ejercicio en la guía 28: Resuelva la ecuación diferencial
𝑦 ′′′ − 6𝑦 ′′ = 3 − 𝑐𝑜𝑠𝑥

3.2.2. Sistemas Masa-Resorte: Movimiento forzado

Movimiento forzado con amortiguamiento: La ecuación diferencial del movimiento forzado con
amortiguamiento es
𝑑2 𝑥 𝑑𝑥
𝑚 2
= −𝑘𝑥 − 𝜑 + 𝑓(𝑡)
𝑑𝑡 𝑑𝑡
Dividiendo entre 𝑚, se obtiene
𝑑2𝑥 𝑑𝑥
2
+ 2𝜆 + 𝜔2 𝑥 = 𝐺(𝑡)
𝑑𝑡 𝑑𝑡
𝑓(𝑡) 17
donde 𝐺(𝑡) = .
𝑚
𝜑 𝑘
Además, 2𝜆 = 𝑚 y 𝜔2 = 𝑚. Si 𝜑 es igual a cero, tenemos un movimiento forzado sin amortiguamiento18.
8
Ejemplo 29: Una masa que pesa 16 libras alarga pie un resorte. Al inicio, la masa se libera desde el
3
reposo en un punto 2 pies abajo de la posición de equilibrio y el movimiento posterior ocurre en un medio
1
que ofrece una fuerza de amortiguamiento igual a de la velocidad instantánea. Encuentre la ecuación
2
del movimiento si se aplica a la masa una fuerza externa igual a 𝑓(𝑡) = 10 cos 2𝑡.

17 Cuando 𝐺 es una función periódica, la solución general de la ecuación diferencial es la suma de una función no periódica
𝑥ℎ (𝑡) y una función periódica 𝑥𝑝 (𝑡). Además lim 𝑥ℎ (𝑡) = 0. Por lo tanto, para valores de 𝑡 grandes, los desplazamientos de la
𝑡→∞
masa se aproximan mediante la solución 𝑥𝑝 (𝑡). 𝑥ℎ (𝑡) se denomina término transitorio o solución transitoria y la función
𝑥𝑝 (𝑡) se denomina término de estado estable o solución de estado estable.
18 Cuando se ejerce una fuerza periódica sin amortiguamiento, no hay término transitorio en la solución. También se observa

que una fuerza periódica con una frecuencia cercana o igual que la frecuencia de las vibraciones libres amortiguadas causa
un problema grave en un sistema mecánico oscilatorio.
43
 𝑊 16 1
Solución: Calculando la masa: 𝑚 = = 32 = 2 slug. La constante de resorte con la ley de Hooke,
𝑔

8
𝐹 = 𝑘𝑠 → 16 = 𝑘 ( )
3
por lo tanto, 𝑘 = 6. Nuestra ecuación diferencial queda
𝑑2 𝑥 𝑑𝑥
+ + 12𝑥 = 20 cos 2𝑡
𝑑𝑡 2 𝑑𝑡
El desplazamiento inicial y la velocidad inicial son 𝑥(0) = 2, 𝑥 ′ (0) = 0.
Resolviendo la ecuación auxiliar de segundo grado
𝑚2 + 𝑚 + 12 = 0
1 √47 1 √47
Sus soluciones son 𝑚 = − 2 ± 𝑖 con 𝛼 = − 2 y 𝛽 = , por lo tanto
2 2

√47 √47
𝑥ℎ (𝑡) = 𝑒 (−1/2)𝑡 (𝑐1 cos ( ) 𝑡 + 𝑐2 sin ( ) 𝑡)
2 2

Suponemos una solución 𝑥𝑝 de la forma

𝑥𝑝 (𝑡) = 𝐴 cos 2𝑡 + 𝐵 sin 2𝑡

Derivando dos veces y reemplazando en la ecuación diferencial

𝑥𝑝′ (𝑡) = −2𝐴 sin 2𝑡 + 2𝐵 cos 2𝑡

𝑥𝑝′′ (𝑡) = −4𝐴 cos 2𝑡 − 4𝐵 sin 2𝑡

−4𝐴 cos 2𝑡 − 4𝐵 sin 2𝑡 + (−2𝐴 sin 2𝑡 + 2𝐵 cos 2𝑡) + 12(𝐴 cos 2𝑡 + 𝐵 sin 2𝑡) = 20 cos 2𝑡
(8𝐴 + 2𝐵) cos 2𝑡 + (8𝐵 − 2𝐴) sin 2𝑡 = 20 cos 2𝑡
Resolviendo el sistema de ecuaciones lineales
8𝐴 + 2𝐵 = 20
{
−2𝐴 + 8𝐵 = 0
40 10
obtenemos, 𝐴 = 17 y 𝐵 = 17 . Por lo tanto la solución particular es

40 10
𝑥𝑝 (𝑡) = cos 2𝑡 + sin 2𝑡
17 17
La ecuación del movimiento queda:

√47 √47 40 10
𝑥(𝑡) = 𝑒 (−1/2)𝑡 (𝑐1 cos ( ) 𝑡 + 𝑐2 sin ( ) 𝑡) + cos 2𝑡 + sin 2𝑡 ∎
2 2 17 17

44
Ejercicio en la guía 29: Una carreta que pesa 128 libras está sujeta a un muro con un resorte de constante
𝑘 = 64 𝑙𝑖𝑏𝑟𝑎𝑠/𝑝𝑖𝑒. Se aparta la carreta 6 pulgadas del muro y se suelta sin velocidad inicial.
Simultáneamente, se le aplica una fuerza periódica externa 𝑓(𝑡) = 32𝑠𝑒𝑛4𝑡. Suponiendo que no hay
resistencia del aire, hallar la posición de la carreta en el instante 𝑡.

CAPÍTULO 4
4. La transformada de Laplace
La transformada de Laplace es un operador fundamental que permite resolver fenómenos naturales de
una manera más elemental con relación a otras formas estudiadas en los capítulos 1 y 2. Se estudiará
su definición y propiedades fundamentales para resolver problemas de valores iniciales.

4.1. Solución de problemas de valor inicial

Para utilizar el método usando la transformada de Laplace, utilizamos el símbolo ℒ para convertir
funciones que dependen de una variable 𝑡 a otra función que depende de la variable 𝑠 de Laplace. Para
realizar la transformación usamos la definición de la transformada de Laplace

∞
ℒ{𝑓(𝑡)} = ∫ 𝑒 −𝑠𝑡 𝑓(𝑡)𝑑𝑡
0

La integral es impropia y resalta que: “Si queremos calcular la transformada de Laplace de cualquier
función 𝑓 (de orden exponencial19), debemos a dicha 𝑓 multiplicarla por 𝑒 −𝑠𝑡 y después hacer la integral
desde cero hasta infinito”.

19 Una función es de orden exponencial 𝑐 si existen constantes 𝑐, 𝑀 > 0 y 𝑇 > 0 tales que |𝑓(𝑡)| ≤ 𝑀𝑒 𝑐𝑡 para toda 𝑡 > 𝑇.
45
Ejemplo 30: Calcular la transformada de Laplace de la función 𝑓(𝑡) = 𝑒 𝑎𝑡 .
Solución: Usando la definición de la transformada de Laplace
∞ 𝑏
1 𝑏 1 1
ℒ{𝑒 𝑎𝑡 }
=∫ 𝑒 𝑒 𝑑𝑡 = 𝑙𝑖𝑚 ∫ 𝑒 (−𝑠+𝑎)𝑡 𝑑𝑡 = lim
−𝑠𝑡 𝑎𝑡
𝑒 (−𝑠+𝑎)𝑡 |0 = 𝑙𝑖𝑚 ( 𝑒 (−𝑠+𝑎)𝑏 − )
0 𝑏→∞ 0 𝑏→∞ −𝑠 + 𝑎 𝑏→∞ −𝑠 + 𝑎 −𝑠 + 𝑎
1
Si 𝑠 > 𝑎, entonces −𝑠+𝑎 𝑒 (−𝑠+𝑎)𝑏 → 0 cuando 𝑏 → ∞. Por lo tanto

1
ℒ{𝑒 𝑎𝑡 } = ; 𝑠>𝑎 ∎
𝑠−𝑎

Ejercicio en la guía 30: Calcule la transformada de Laplace de la función 𝑓(𝑡) = cos 𝑡 usando la definición
de la transformada de Laplace.

𝒔
R/ 𝒔𝟐 +𝟏

Algunas fórmulas importantes

𝑛!
ℒ{𝑡 𝑛 } =
𝑠 𝑛+1

1
ℒ{𝑒 𝑎𝑡 } =
𝑠−𝑎
𝑠
ℒ{cos 𝑘𝑡} =
𝑠2 + 𝑘2

𝑘
ℒ{sin 𝑘𝑡} =
𝑠2 + 𝑘2

46
Propiedad de linealidad de la transformada de Laplace
 Linealidad: La transformada de Laplace es una transformación Lineal, esto quiere decir, que para
cualquier constante 𝛼 y 𝛽:

ℒ{𝛼𝑓(𝑡) ± 𝛽𝑔(𝑡)} = 𝛼ℒ{𝑓(𝑡)} ± 𝛽ℒ{𝑔(𝑡)} = 𝛼𝐹(𝑠) ± 𝛽𝐺(𝑠)

Transformada de la derivada:

Transformada de una derivada: Si 𝑓 (𝑛) (𝑡) es la 𝑛 −ésima derivada de una función, entonces

ℒ{𝑓 (𝑛) (𝑡) } = 𝑠 𝑛 𝐹(𝑠) − 𝑠 𝑛−1 𝑓(0) − 𝑠 𝑛−2 𝑓 ′ (0) − ⋯ − 𝑓 (𝑛−1) (0)

donde ℒ{𝑓(𝑡)} = 𝐹(𝑠).

A partir de la fórmula de la transformada de una derivada podemos concluir

𝑑𝑦
 ℒ { 𝑑𝑡 } = 𝑠𝑌(𝑠) − 𝑦(0)
𝑑2 𝑦
 ℒ { 𝑑𝑡 2 } = 𝑠 2 𝑌(𝑠) − 𝑠𝑦(0) − 𝑦′(0)

donde ℒ{𝑦(𝑡)} = 𝑌(𝑠).

Transformada inversa de Laplace:

Si ℒ{𝑦(𝑡)} = 𝑌(𝑠) entonces definimos la transformada inversa de

Laplace como:
ℒ −1 {𝑌(𝑠)} = 𝑦(𝑡)

Ejercicio en la guía 31: Usando las fórmulas importantes, calcular la transformada inversa de Laplace de
4 1
las siguientes funciones 𝐹(𝑠) = 2𝑠−1 y 𝐺(𝑠) = 𝑠2 +7

𝟏
𝟏
R/ 𝒇(𝒕) = 𝟐𝒆𝟐𝒕 y 𝒈(𝒕) = 𝐬𝐢𝐧 √𝟕 𝒕
√𝟕

47
Ejemplo 31: Encuentre 𝑌(𝑠) del siguiente problema de valores iniciales:
𝑡, 𝑡 < 2
𝑦 ′′ + 4𝑦′ = 𝑔(𝑡), donde 𝑔(𝑡) = { , sujeto a las condiciones 𝑦(0) = −1 ; 𝑦´(0) = 0.
5, 𝑡 > 2
Solución: Aplicando la transformada de Laplace en ambos lados de la ecuación diferencial
ℒ{𝑦′′} + ℒ{4𝑦} = ℒ{𝑔(𝑡)}
2 ∞
= 𝑠 2 𝑌(𝑠) − 𝑠𝑦(0) − 𝑦 ′ (0) + 4(𝑠𝑌(𝑠) − 𝑦(0)) = ∫ 𝑒 −𝑠𝑡 𝑡 𝑑𝑡 + ∫ 𝑒 −𝑠𝑡 5 𝑑𝑡
0 2

Resolviendo la primera integral por partes, la segunda integral impropia y reemplazando las condiciones
iniciales
2 1 1 5
= 𝑠 2 𝑌(𝑠) + 𝑠 + 4𝑠𝑌(𝑠) + 4 = − 𝑒 −2𝑠 − 2 𝑒 −2𝑠 + 2 + 𝑒 −2𝑠
𝑠 𝑠 𝑠 𝑠
Despejando 𝑌(𝑠)
3𝑠𝑒 −2𝑠 − 𝑒 −2𝑠 + 1 − 𝑠 3 − 4𝑠 2
𝑌(𝑠) = ∎
𝑠 3 (𝑠 + 4)

Ejemplo 32: Resolver la siguiente ecuación diferencial usando la transformada de Laplace

𝑦 ′ − 4𝑦 = 1 , 𝑦(0) = 0
Solución: Transformando ambos lados de la ecuación diferencial
ℒ{𝑦 ′ − 4𝑦} = ℒ{1}
Aplicando la linealidad de la transformada
ℒ{𝑦′} − 4ℒ{𝑦} = ℒ{1}
Utilizando la transformada de la derivada de primer orden y otras transformadas elementales
1
𝑠𝑌(𝑠) − 𝑦(0) − 4𝑌(𝑠) =
𝑠
Aplicando la condición inicial y despejando 𝑌(𝑠):
1
𝑌(𝑠) =
𝑠(𝑠 − 4)
Separando en fracciones parciales
1 −1 1
= +
𝑠(𝑠 − 4) 4𝑠 4(𝑠 − 4)
La solución del problema de valor inicial es
−1 1 −1 1 1 1
𝑦(𝑡) = ℒ −1 { + } = ℒ −1 { } + ℒ −1 { } = − + 𝑒 4𝑡 ∎
4𝑠 4(𝑠 − 4) 4𝑠 4(𝑠 − 4) 4 4
48
Ejercicio en la guía 32: Resolver la ecuación diferencial usando algún método del capítulo 3 y también la
transformada de Laplace
𝑦 ′′ − 2𝑦 ′ + 𝑦 = 𝑒 𝑥 ; 𝑦(0) = 1 , 𝑦 ′ (0) = 0

Ejercicio en la guía 33: Cómo resolver un problema de valor inicial usando la transformada de Laplace cuya
condición inicial no esté en 𝑡 = 0. Ejemplo:
𝑦 ′′ − 2𝑦 ′ + 𝑦 = 6𝑡 − 2 ; 𝑦(−1) = 3 , 𝑦 ′ (−1) = 7

49
 4.2. Transformada de Laplace de funciones periódicas

Una función 𝑓(𝑡) se denomina función periódica si cumple la igualdad

𝑓(𝑡 + 𝑇) = 𝑓(𝑡)
donde 𝑇 se denomina el periodo.
Si 𝑓(𝑡) es una función periódica, entonces la transformada de Laplace de dicha función se define
mediante la expresión

Si 𝑓 es una función periódica con periodo 𝑇, entonces

𝑇
1
ℒ{(𝑓(𝑡)} = −𝑠𝑇
∫ 𝑒 −𝑠𝑡 𝑓(𝑡)𝑑𝑡
1−𝑒 0

Como 𝑓(𝑡) = 𝑠𝑒𝑛𝑡 es una función periódica con periodo 𝑇 = 2𝜋, entonces
2𝜋
1 1
ℒ{𝑠𝑒𝑛𝑡} = −𝑠(2𝜋)
∫ 𝑒 −𝑠𝑡 𝑠𝑒𝑛𝑡 𝑑𝑡 = 2
1−𝑒 0 𝑠 +1

Ejercicio en la guía 34: Verificar el resultado anterior usando la transformada de una función periódica.

50
Ejemplo 33: Determine ℒ{𝑓(𝑡)}, donde 𝑓(𝑡) = 𝑒 𝑡 , 0 < 𝑡 < 1, y 𝑓(𝑡) es una función periódica con periodo
𝑇 = 1.
Solución: Usando el resultado para funciones periódicas
1
1
ℒ{𝑒 𝑡 } = ∫ 𝑒 −𝑠𝑡 𝑒 𝑡 𝑑𝑡
1 − 𝑒 −𝑠(1) 0

1 1 1
1 1 1
∫ 𝑒 −𝑠𝑡 𝑡
𝑒 𝑑𝑡 = ∫ 𝑒 −𝑠𝑡+𝑡
𝑑𝑡 = ∫ 𝑒 (−𝑠+1)𝑡 𝑑𝑡
1 − 𝑒 −𝑠(1) 0 1 − 𝑒 −𝑠 0 1 − 𝑒 −𝑠 0
Calculando la integral respectiva, obtenemos
1
1 1 (−𝑠+1)𝑡 1 1 (−𝑠+1) 1
( 𝑒 )| = ( 𝑒 − ) ∎
1 − 𝑒 −𝑠 −𝑠 + 1 0 1 − 𝑒 −𝑠 −𝑠 + 1 −𝑠 + 1

Ejercicio en la guía 35: Calcular la transformada de Laplace de la función y periodo dado

1 0≤𝑡<𝑎
𝑓(𝑡) = { ; 𝑇 = 2𝑎
0 𝑎 ≤ 𝑡 < 2𝑎

51
 4.3. Traslación en el eje 𝒕

Transformada de la función escalón unitario: La función escalón unitario o función Heaviside se

define como
0, 𝑡<𝑎
𝑈(𝑡 − 𝑎) = {
1, 𝑡≥𝑎
Cuando una función 𝑓 se multiplica por una función escalón unitario, lo que resulta es una anulación de
la función para los 𝑡 < 𝑎, esto quiere decir
0, 𝑡<𝑎
𝑓(𝑡) ∙ 𝑈(𝑡 − 𝑎) = {
𝑓(𝑡), 𝑡≥𝑎
0, 0 ≤ 𝑡 < 1
Podemos ver la función 𝑉(𝑡) = { de la forma 𝑉(𝑡) = 5𝑈(𝑡 − 1).
5, 𝑡 ≥ 1
Para encontrar la transformada de Laplace de funciones escalón unitario usamos el siguiente resultado

ℒ{𝑓(𝑡)𝑈(𝑡 − 𝑎)} = 𝑒 −𝑎𝑠 ℒ{𝑓(𝑡 + 𝑎)}

Si 𝑓(𝑡) = 5 entonces 𝑓(𝑡 + 𝑎) = 5, por lo tanto

−𝑠 −𝑠
5 5𝑒 −𝑠
ℒ{5𝑈(𝑡 − 1)} = 𝑒 ℒ{5} = 𝑒 ( )=
𝑠 𝑠
Por lo tanto el modelo de la corriente queda
5𝑒 −𝑠
𝑠𝐼(𝑠) + 𝐼(𝑠) = ℒ{𝑉(𝑡)} =
𝑠
Despejando 𝐼(𝑠)
5𝑒 −𝑠
𝐼(𝑠) =
𝑠(𝑠 + 1)
Para calcular la transformada inversa de Laplace para 𝐼(𝑠), utilizamos el siguiente resultado

ℒ −1 {𝑒 −𝑎𝑠 𝐹(𝑠)} = 𝑓(𝑡 − 𝑎)𝑈(𝑡 − 𝑎)

Por lo tanto
5𝑒 −𝑠
𝑖(𝑡) = ℒ −1 {𝐼(𝑠)} = ℒ −1 { }
𝑠(𝑠 + 1)

52
 5
Observemos que 𝐹(𝑠) = 𝑠(𝑠+1) , separando en fracciones parciales y la inversa de Laplace obtenemos
𝑓(𝑡) = 5 − 5𝑒 −𝑡 y por lo tanto 𝑓(𝑡 − 1) = 5 − 5𝑒 −(𝑡−1) . Finalmente tenemos
5𝑒 −𝑠
ℒ −1 { } = [5 − 5𝑒 −(𝑡−1) ]𝑈(𝑡 − 1)
𝑠(𝑠 + 1)
Por lo tanto la corriente en cualquier instante es

𝑖(𝑡) = [5 − 5𝑒 −(𝑡−1) ]𝑈(𝑡 − 1)

Otras formas de representar funciones a trozos con respecto a la función escalón unitario:

𝑔(𝑡), 0 ≤ 𝑡 < 𝑎
Si 𝑓(𝑡) = {
ℎ(𝑡), 𝑡 ≥ 𝑎
Entonces
𝑓(𝑡) = 𝑔(𝑡) − 𝑔(𝑡)𝑈(𝑡 − 𝑎) + ℎ(𝑡)𝑈(𝑡 − 𝑎)

0, 0≤𝑡<𝑎
Si 𝑓(𝑡) = {1, 𝑎≤𝑡<𝑏
0, 𝑡≥𝑏
Entonces
𝑓(𝑡) = ℎ(𝑡)[𝑈(𝑡 − 𝑎) − 𝑈(𝑡 − 𝑏)]

Ejemplo 34: Resolver el problema con valores iniciales

𝑦 ′′ + 𝑦 = 𝑈(𝑡 − 2) − 𝑈(𝑡 − 4) ; 𝑦(0) = 1, 𝑦 ′ (0) = 0
Solución: Usando la transformada de Laplace en ambos lados de la ecuación diferencial
ℒ{𝑦′′} + ℒ{𝑦} = ℒ{𝑈(𝑡 − 2)} − ℒ{𝑈(𝑡 − 4)}
Calculando las transformadas de Laplace
𝑒 −2𝑠 𝑒 −4𝑠
𝑠 2 𝑌(𝑠) − 𝑠𝑦(0) − 𝑦 ′ (0) + 𝑌(𝑠) = −
𝑠 𝑠
Utilizando las condiciones iniciales y despejando 𝑌(𝑠)
𝑒 −2𝑠 𝑒 −4𝑠 𝑠
𝑌(𝑠) = − +
𝑠(𝑠 2 + 1) 𝑠(𝑠 2 + 1) 𝑠 2 + 1

53
Calculando las transformadas inversas de Laplace
𝑒 −2𝑠
 ℒ −1 {𝑠(𝑠2 +1)} = 1 − 𝑠𝑒𝑛(𝑡 − 2)
𝑒 −4𝑠
 ℒ −1 {𝑠(𝑠2 +1)} = 1 − 𝑠𝑒𝑛(𝑡 − 4)
𝑠
 ℒ −1 {𝑠2 +1} = 𝑐𝑜𝑠𝑡

Por lo tanto
𝑦(𝑡) = 𝑠𝑒𝑛(𝑡 − 4) − 𝑠𝑒𝑛(𝑡 − 2) + 𝑐𝑜𝑠𝑡 ∎

2
(1+𝑒 −2𝑠 )
Ejercicio en la guía 36: Calcule la transformada inversa de Laplace de la función 𝐹(𝑠) = 𝑠+2

Ejemplo 35: Resolver el problema con valores iniciales

𝑦 ′′ + 5𝑦 ′ + 6𝑦 = 𝑔(𝑡)
0, 0 ≤ 𝑡 < 1
con 𝑔(𝑡) = { 𝑡, 1 ≤ 𝑡 < 5 , además, 𝑦(0) = 0, 𝑦 ′ (0) = 2.
1, 𝑡≥5
Solución: Si escribimos 𝑔(𝑡) en términos de la función escalón unitario
𝑔(𝑡) = 𝑡𝑈(𝑡 − 1) − 𝑡𝑈(𝑡 − 5) + 𝑈(𝑡 − 5) = 𝑡𝑈(𝑡 − 1) − (𝑡 − 1)𝑈(𝑡 − 5)
Por lo tanto
ℒ{𝑦′′} + 5ℒ{𝑦′} + 6ℒ{𝑦} = ℒ{𝑡𝑈(𝑡 − 1)} − ℒ{(𝑡 − 1)𝑈(𝑡 − 5)}

𝑠 2 𝑌(𝑠) − 𝑠𝑦(0) − 𝑦′(0) + 5(𝑠𝑌(𝑠) − 𝑦(0)) + 6𝑌(𝑠) = 𝑒 −𝑠 ℒ{𝑡 + 1} − 𝑒 −5𝑠 ℒ{𝑡}

2 −𝑠
1 1 𝑒 −5𝑠
𝑠 𝑌(𝑠) − 2 + 5𝑠𝑌(𝑠) + 6𝑌(𝑠) = 𝑒 ( 2+ )− 2
𝑠 𝑠 𝑠
despejando 𝑌(𝑠)

54
 𝑠+1 𝑒 −5𝑠
𝑌(𝑠) = 𝑒 −𝑠 [ ] −
𝑠 2 (𝑠 2 + 5𝑠 + 6) 𝑠 2 (𝑠 2 + 5𝑠 + 6)

Ejercicio en la guía 37: Encontrar 𝑦(𝑡) para finalizar el ejemplo 35.

4.4. Más sobre la transformada de Laplace

 Traslación en el eje 𝒔: Un resultado muy importante denominado la traslación en el eje 𝑠

ℒ{𝑒 𝑎𝑡 𝑓(𝑡)} = 𝐹(𝑠 − 𝑎)

donde 𝐹(𝑠 − 𝑎) = ℒ{𝑓(𝑡)}|𝑠→𝑠−𝑎

Ejemplo 36: Calcular ℒ{𝑒 −2𝑡 𝑐𝑜𝑠𝑡}

Solución:
𝑠 𝑠+2
ℒ{𝑒 −2𝑡 𝑐𝑜𝑠𝑡} = ℒ{𝑐𝑜𝑠𝑡}|𝑠→𝑠−(−2) = { }| = ∎
𝑠2 + 1 𝑠→𝑠+2 (𝑠 + 2)2 + 1
A partir de la traslación en el eje 𝑠, tenemos:

ℒ −1 {𝐹(𝑠 − 𝑎)} = 𝑒 𝑎𝑡 𝑓(𝑡)

 Derivada de una transformada: Se define la derivada de una transformada mediante la formula

𝑑𝑛
ℒ{𝑡 𝑛 𝑓(𝑡)} = (−1)𝑛 𝑛 {𝐹(𝑠)}
𝑑𝑠

55
 donde ℒ{𝑓(𝑡)} = 𝐹(𝑠)

Ejemplo 37: Calcular ℒ{𝑡 2 𝑒 3𝑡 }

Solución: Si observamos la derivada de una transformada, 𝑛 = 2 y por lo tanto
𝑑2 1
ℒ{𝑡 2 𝑒 3𝑡 } = (−1)2 2
{ }
𝑑𝑠 𝑠 − 3
1 𝑑2 1 2
Si 𝐹(𝑠) = 𝑠−3 entonces 𝑑𝑠2 {𝑠−3} = (𝑠−3)3 , por lo tanto

2 3𝑡 }
𝑑2 1 2
ℒ{𝑡 𝑒 = (−1)2 { }= ∎
2
𝑑𝑠 𝑠 − 3 (𝑠 − 3)3

Si en el ejercicio del ejemplo 37, utilizamos la traslación en el eje 𝑠 obtenemos el mismo resultado
2! 2
ℒ{𝑡 2 𝑒 3𝑡 } = ℒ{𝑡 2 }|𝑠→𝑠−3 = { 3 }| =
𝑠 𝑠→𝑠−3 (𝑠 − 3)3

Ejercicio en la guía 38: Calcular la transformada de Laplace de la función 𝑓(𝑡) = 𝑡 3 𝑒 −𝑡 usando la definición
de la traslación en el eje 𝑠 y la derivada de una transformada.

 Convolución: Sea 𝑓 y 𝑔 dos funciones continuas, se define un producto especial 𝑓 ∗ 𝑔 de la

forma
𝑡
𝑓 ∗ 𝑔 = ∫ 𝑓(𝜏)𝑔(𝑡 − 𝜏)𝑑𝜏
0

No es complicado verificar que 𝑓 ∗ 𝑔 = 𝑔 ∗ 𝑓. Además

56
 ℒ{𝑓 ∗ 𝑔} = ℒ{𝑓(𝑡)} ∙ ℒ{𝑔(𝑡)} = 𝐹(𝑠)𝐺(𝑠)
ℒ −1 {𝐹(𝑠)𝐺(𝑠)} = 𝑓 ∗ 𝑔

Ejemplo 38: Calcular ℒ{𝑒 2𝑡 ∗ 𝑠𝑒𝑛𝑡}

1 1
Solución: Sabemos que ℒ{𝑒 2𝑡 } = 𝑠−2 y ℒ{𝑠𝑒𝑛𝑡} = 𝑠2 +1 . Por lo tanto

1 1 1
ℒ{𝑒 2𝑡 ∗ 𝑠𝑒𝑛𝑡} = ℒ{𝑒 2𝑡 } ∙ ℒ{𝑠𝑒𝑛𝑡} = ∙ 2 = ∎
𝑠 − 2 𝑠 + 1 (𝑠 − 2)(𝑠 2 + 1)

1
Ejemplo 39: Calcular ℒ −1 {𝑠2 −𝑠}
1
Solución: Si separamos la fracción 𝑠2 −𝑠 en suma de fracciones parciales obtenemos

1 −1 1
= +
𝑠2 −𝑠 𝑠 𝑠−1
y por lo tanto
1 −1 1
ℒ −1 { 2
} = ℒ −1
{ } + ℒ −1
{ } = −1 + 𝑒 𝑡
𝑠 −𝑠 𝑠 𝑠−1
Si utilizamos la convolución:
1 1 1 1
ℒ −1 { } = ℒ −1
{ } = ℒ −1
{ ∙ } = 1 ∗ 𝑒 𝑡 = −1 + 𝑒 𝑡 ∎
𝑠2 − 𝑠 𝑠(𝑠 − 1) 𝑠 𝑠−1

Ejercicio en la guía 39: Mostrar que 1 ∗ 𝑒 𝑡 = −1 + 𝑒 𝑡

57
  Transformada de una integral: Utilizando el resultado de convolución

𝑡
𝑓 ∗ 𝑔 = ∫ 𝑓(𝜏)𝑔(𝑡 − 𝜏)𝑑𝜏
0
1
Si 𝑔(𝑡) = 1 entonces ℒ{1} = 𝑠 y por lo tanto
𝑡
1 𝐹(𝑠)
ℒ{𝑓 ∗ 1} = ℒ {∫ 𝑓(𝜏)𝑑𝜏} = 𝐹(𝑠) ∙ =
0 𝑠 𝑠

Además
𝑡
𝐹(𝑠)
ℒ −1 { } = ∫ 𝑓(𝜏)𝑑𝜏
𝑠 0

Ejemplo 40: Calcular

𝑡
ℒ {𝑡 ∫ 𝑠𝑒𝑛𝜏 𝑑𝜏}
0

Solución: Aplicando la derivada de una transformada con 𝑛 = 1

𝑡 𝑡
𝑑 𝑑 1 3𝑠 2 + 1
ℒ {𝑡 ∫ 𝑠𝑒𝑛𝜏 𝑑𝜏} = − (ℒ {∫ 𝑠𝑒𝑛𝜏 𝑑𝜏}) = − ( 2 )= 2 2
0 𝑑𝑠 0 𝑑𝑠 𝑠(𝑠 + 1) 𝑠 (𝑠 + 1)2

Por lo tanto
𝑡
3𝑠 2 + 1
ℒ {𝑡 ∫ 𝑠𝑒𝑛𝜏 𝑑𝜏} = 2 2 ∎
0 𝑠 (𝑠 + 1)2

Ejemplo 41: Resolver la ecuación integro-diferencial

𝑡
𝑡
𝑓(𝑡) = 𝑡𝑒 + ∫ 𝜏𝑓(𝑡 − 𝜏)𝑑𝜏
0

Solución: Debemos encontrar el valor de 𝑓(𝑡). Aplicando la transformada de Laplace en ambos lados
𝑡
𝑡
ℒ{𝑓(𝑡)} = ℒ {𝑡𝑒 + ∫ 𝜏𝑓(𝑡 − 𝜏)𝑑𝜏}
0
𝑡
𝐹(𝑠) = ℒ{𝑡𝑒 𝑡 } + ℒ {∫ 𝜏𝑓(𝑡 − 𝜏)𝑑𝜏}
0

Calculando las dos transformadas de Laplace indicadas

58
 𝑑 𝑑 1 1
 ℒ{𝑡𝑒 𝑡 } = − 𝑑𝑠 {ℒ{𝑒 𝑡 }} = − 𝑑𝑠 {𝑠−1} = (𝑠−1)2
𝑡 1
 ℒ {∫0 𝜏𝑓(𝑡 − 𝜏)𝑑𝜏} = ℒ{𝑡 ∗ 𝑓(𝑡)} = ℒ{𝑡} ∙ ℒ{𝑓(𝑡)} = 𝑠2 ∙ 𝐹(𝑠)

Por lo tanto
𝑡
1 1
𝐹(𝑠) = ℒ{𝑡𝑒 𝑡 } + ℒ {∫ 𝜏𝑓(𝑡 − 𝜏)𝑑𝜏} = + 2 ∙ 𝐹(𝑠)
0 (𝑠 − 1) 2 𝑠

despejando 𝐹(𝑠)
1 1
𝐹(𝑠) − ∙ 𝐹(𝑠) =
𝑠 2 (𝑠 − 1)2
𝑠2 𝑠2
𝐹(𝑠) = =
(𝑠 2 − 1)(𝑠 − 1)2 (𝑠 + 1)(𝑠 − 1)3
Separando en fracciones parciales
𝑠2 𝐴 𝐵 𝐶 𝐷
= + + +
(𝑠 + 1)(𝑠 − 1)3 𝑠 − 1 𝑠 + 1 (𝑠 − 1)2 (𝑠 − 1)3
Podemos encontrar que 𝐴 = 1/8, 𝐵 = −1/8, 𝐶 = 3/4, 𝐷 = 1/2. Reemplazando
𝑠2 1/8 1/8 3/4 1/2
𝐹(𝑠) = = − + +
(𝑠 + 1)(𝑠 − 1)3 𝑠 − 1 𝑠 + 1 (𝑠 − 1)2 (𝑠 − 1)3
Como 𝑓(𝑡) = ℒ −1 {𝐹(𝑠)}, tenemos
1 1 3 1
𝑓(𝑡) = 𝑒 𝑡 − 𝑒 −𝑡 + 𝑡𝑒 𝑡 + 𝑡 2 𝑒 𝑡 ∎
8 8 4 4

Ejercicio en la guía 40: Use la transformada de Laplace para resolver la ecuación integral o ecuación
integrodiferencial
𝑡
𝑓(𝑡) = 2𝑡 − 4 ∫ sin 𝜏 𝑓(𝑡 − 𝜏) 𝑑𝜏
0

59
EJERCICIOS DE ECUACIONES
DIFERENCIALES

60
1. Generalidades, soluciones y PVI

1) Resuelva la siguiente situación

a) Muestre que la familia de soluciones de la ecuación diferencial

𝑑𝑦 2
= (𝑥 − 3)(𝑦 + 1)3
𝑑𝑥

es
3
𝑥2
𝑦(𝑥) = −1 + ( − 𝑥 + 𝐾)
6

R/ Demostración

b) Muestre que 𝑦 = −1 satisface la ecuación diferencial del ítem (a)

R/ Demostración

c) Muestre que no hay una elección de la constante 𝐾 de modo que la solución del ítem (a) produzca
la solución 𝑦 = −1. Así, perdemos la solución 𝑦 = −1 al dividir entre (𝑦 + 1)2/3 .
R/ Demostración

2) Determine si la siguiente función

𝑐
𝑦 = 𝑐2 + ; 𝑐 𝑒𝑠 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝑥

es solución de la ecuación diferencial

𝑦 + 𝑥𝑦 ′ = 𝑥 4 (𝑦′)2
R/ Si

3) Resuelva el siguiente problema de análisis

a) Muestre que 𝑦 = 𝑒 −𝑥 es una solución de la ecuación lineal

𝑑𝑦
+𝑦 =0
𝑑𝑥
y 𝑦 = 𝑥 −1 es una solución de la ecuación no lineal
𝑑𝑦
+ 𝑦2 = 0
𝑑𝑥
R/ Demostración

61
 b) Muestre que para cualquier constante 𝐾, 𝐾𝑒 −𝑥 es una solución de la ecuación lineal, mientras que
𝐾𝑥 −1 es una solución de la ecuación no lineal sólo cuando 𝐾 = 0 ó 1.
R/ Demostración

c) Muestre que para cualquier ecuación lineal de la forma

𝑑𝑦
+ 𝑃(𝑥)𝑦 = 0
𝑑𝑥
Si 𝑔(𝑥) es una solución, entonces para cualquier constante 𝐾, la función 𝐾𝑔(𝑥) también es una
solución.
R/ Demostración

4) Muestre que 𝑦 = 𝑒 𝑥 − 𝑥 es una solución explícita de la ecuación diferencial

𝑑𝑦
+ 𝑦 2 = 𝑒 2𝑥 + (1 − 2𝑥)𝑒 𝑥 + 𝑥 2 − 1
𝑑𝑥
R/ Demostración

5) Determine si la función dada es una solución de la ecuación diferencial correspondiente.

a) 𝑦 = sin 𝑥 + 𝑥 2 , 𝑦 ′′ + 𝑦 = 𝑥 2 + 2
R/ Si

b) 𝑥 = 3 cos 𝑡 − 5 sin 𝑡 , 𝑥 ′′ + 𝑥 = 0
R/ Si

c) 𝑥 = cos 2𝑡 , 𝑥 ′ + 𝑡𝑥 = sin 2𝑡
R/ No

d) 𝜃 = 2𝑒 3𝑡 − 𝑒 2𝑡 , 𝜃 ′′ − 𝜃𝜃 ′ + 3𝜃 = −2𝑒 2𝑡
R/ No

e) 𝑦 = 𝑒 2𝑥 − 3𝑒 −𝑥 , 𝑦 ′′ − 𝑦 ′ − 2𝑦 = 0
R/ Si

f) 𝑦 = 3 sin 2𝑥 + 𝑒 −𝑥 , 𝑦 ′′ + 4𝑦 = 5𝑒 −𝑥
R/ Si
2 𝑥 2 2
g) 𝑦 = 𝑒 −𝑥 ∫0 𝑒 𝑡 𝑑𝑡 + 𝑐𝑒 −𝑥 , 𝑦 ′ + 2𝑥𝑦 = 1
R/ Si

62
 h) 𝑦 = 𝑐1 𝑥 −1 + 𝑐2 𝑥 + 𝑐3 𝑥 ln 𝑥 + 4𝑥 2 , 𝑥 3 𝑦 ′′′ + 2𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ + 𝑦 = 12𝑥 2

R/ Si

𝑑𝑦 𝑥𝑦
i) 𝑥 2 + 𝑐𝑦 2 = 1 , = 𝑥 2 −1
𝑑𝑥

R/ Si
6) Determine si la relación dada es una solución implícita de la ecuación diferencial correspondiente.

𝑑𝑦 𝑥
a) 𝑥 2 + 𝑦 2 = 8 ; =𝑦
𝑑𝑥
R/ No

𝑑𝑦 2𝑥𝑦
b) 𝑦 − ln 𝑦 = 𝑥 2 + 1 , = 𝑦−1
𝑑𝑥
R/ Si

𝑑𝑦 𝑒 −𝑥𝑦 −𝑦
c) 𝑒 𝑥𝑦 + 𝑦 = 𝑥 − 1 , = 𝑒 −𝑥𝑦 +𝑥
𝑑𝑥
R/ Si

𝑑𝑦
d) 𝑥 2 − sin(𝑥 + 𝑦) = 1 , = 2𝑥 sec(𝑥 + 𝑦) − 1
𝑑𝑥
R/ Si

6𝑥𝑦´+(𝑦´)3 sin 𝑦−2(𝑦´)2

e) sin 𝑦 + 𝑥𝑦 − 𝑥 3 = 2 , 𝑦 ′′ = 3𝑥 2 −𝑦
R/ Si

2𝑥−1 𝑑𝑥
f) ln ( 𝑥−1 ) = 𝑡 , = (𝑥 − 1)(1 − 2𝑥)
𝑑𝑡
R/ Si

g) −2𝑥 2 𝑦 + 𝑦 2 = 1 , 2𝑥𝑦 𝑑𝑥 + (𝑥 2 − 𝑦)𝑑𝑦 = 0

R/ Si

7) Determine los valores de 𝑚 para los que la función 𝜇(𝑥) = 𝑒 𝑚𝑥 es una solución de la ecuación
diferencial

𝑦 ′′′ + 3𝑦 ′′ + 2𝑦 ′ = 0
R/ −𝟐, −𝟏, 𝟎

8) Determine los valores de 𝑚 para los que la función 𝜇(𝑥) = 𝑥 𝑚 es una solución de la ecuación
diferencial
𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ − 5𝑦 = 0
R/ 𝟏 + √𝟔 , 𝟏 − √𝟔

63
9) Encuentre una solución del PVI de primer orden que consiste en la ecuación diferencial 𝑦 ′ = 𝑦 − 𝑦 2
y la condición inicial 𝑦(−1) = 2, donde 𝑦 = 1/(1 + 𝑐𝑒 −𝑥 ) es una familia uniparamétrica de soluciones
de la ecuación diferencial.
𝟐
R/ 𝒚 = 𝟐−𝒆−(𝒙+𝟏)

2. Ecuaciones diferenciales de primer orden

2.1. Variables separables

1. Resuelva

𝑑𝑦
a) + 2𝑥𝑦 2 = 0
𝑑𝑥
𝟏
R/ 𝒚 = 𝒙𝟐 +𝒄

𝑑𝑁
b) + 𝑁 = 𝑁𝑡𝑒 𝑡+2
𝑑𝑡
𝒕+𝟐 −𝒆𝒕+𝟐 −𝒕
R/ 𝑵 = 𝒄𝒆𝒕𝒆

2. Resolver el siguiente problema de valor inicial.

𝑑𝑇
= 𝑘(𝑇 − 𝑇𝑚 ); 𝑇(0) = 𝑇0 , 𝑘, 𝑇𝑚 y 𝑇0 son constantes
𝑑𝑡

R/ 𝑻(𝒕) = 𝑻𝒎 + (𝑻𝟎 − 𝑻𝒎 )𝒆𝒌𝒕

3. Determine si la ecuación diferencial dada es separable

𝑑𝑦
a) = 2𝑦 4 + 𝑦 − 3
𝑑𝑥
R/ Si

𝑑𝑦 𝑦𝑒 𝑥+𝑦
b) =
𝑑𝑥 𝑥 2 +2
R/ Si

𝑑𝑠 𝑠+1
c) 𝑠 2 + 𝑑𝑡 = 𝑠𝑡
R/ No

d) (𝑥𝑦 2 + 3𝑦 2 )𝑑𝑦 − 2𝑥 𝑑𝑥 = 0
64
 R/ Si

4. Resuelva la ecuación

𝑑𝑦
a) = 𝑦(3 + sin 𝑥)
𝑑𝑥
R/ 𝒚 = 𝒄𝒆𝟑𝒙−𝐜𝐨𝐬 𝒙

𝑑𝑦 1
b) = 𝑥𝑦 3
𝑑𝑥
R/ 𝒚𝟒 = 𝟒 𝐥𝐧 𝒙 + 𝒄

𝑑𝑣 1−4𝑣 2
c) 𝑥 =
𝑑𝑥 3𝑣
𝟖
R/ 𝐥𝐧|𝟏 − 𝟒𝒗𝟐 | = − 𝟑 𝐥𝐧|𝒙| + 𝒄

𝑑𝑥
d) + 𝑥2 = 𝑥
𝑑𝑡
𝒄𝒆𝒕
R/ 𝒙 = 𝟏+𝒄𝒆𝒕

e) 𝑑𝑥 + 𝑒 3𝑥 𝑑𝑦 = 0
𝑑𝑃
f) = 𝑃 − 𝑃2
𝑑𝑡
𝑑𝑦 𝑥𝑦+2𝑦−𝑥−2
g) =
𝑑𝑥 𝑥𝑦−3𝑦+𝑥−3

𝑑𝑥
h) = 3𝑥𝑡 2
𝑑𝑡
𝟑
R/ 𝒙(𝒕) = 𝒄𝒆𝒕

i) 𝑥 2 𝑑𝑥 + 2𝑦 𝑑𝑦 = 0

5. Un modelo de poblaciones utilizado en las predicciones actuariales se basa en la ecuación de

Gompertz
𝑑𝑃
= 𝑃(𝑎 − 𝑏 ln 𝑃)
𝑑𝑡
donde 𝑎 y 𝑏 son constantes

a) Halle 𝑃(𝑡) en la ecuación de Gompertz.

−𝒃𝒕 +𝒂
R/ 𝑷(𝒕) = 𝒆𝒄𝒆 𝒃

b) Si 𝑃(0) = 𝑃0 > 0, dé una fórmula para 𝑃(𝑡) en términos de 𝑎, 𝑏, 𝑃0 y 𝑡.

𝒂 𝟏−𝒆−𝒃𝒕
−𝒃𝒕
R/ 𝑷(𝒕) = 𝑷𝟎 𝒆 (𝒆 )
𝒃

65
2.1.1. Modelo de crecimiento poblacional y desintegración
radiactiva

1) La población de una comunidad crece a una razón proporcional a la población en el tiempo 𝑡. La

población inicial 𝑃0 de 500 aumenta 20% en 10 años. ¿Cuál será la población pasados 30 años? ¿Qué
tan rápido está creciendo la población en 𝑡 = 35?
R/ 864 personas, 17.26 personas/año

2) Se sabe que la población de una comunidad crece con una razón proporcional al número de personas
presentes en el tiempo 𝑡. Si la población inicial 𝑃0 se duplicó en 5 años, ¿en cuánto tiempo se triplicará
y cuadruplicará?
R/ 7.92 años, 10 años

3) Suponga que se sabe que la población de una comunidad del problema 2 es de 10000 después de
tres años. ¿Cuál era la población inicial 𝑃0 ? ¿Cuál será la población en 10 años? ¿Qué tan rápido está
creciendo la población en 𝑡 = 10?
R/ 6598 personas, 26390 personas, 3658.45 personas/año

4) La población de un pueblo crece con una razón proporcional a la población en el tiempo 𝑡. La población
inicial de 500 aumenta 15% en 10 años. ¿Cuál será la población pasados 30 años? ¿Qué tan rápido
está creciendo la población en 𝑡 = 30?
R/ 760 personas, 10.63 personas/año

5) La población de bacterias en un cultivo crece a una razón proporcional a la cantidad de bacterias

presentes en el tiempo 𝑡. Después de dos horas se observa que hay 400 bacterias presentes. Después
de 9 horas hay 3000 bacterias presentes. ¿Cuál era la cantidad inicial de bacterias en el cultivo?
R/ 𝟐𝟐𝟒. 𝟗𝟑 𝒃𝒂𝒄𝒕𝒆𝒓𝒊𝒂𝒔

6) Un cultivo de bacterias de población 𝑃 crece a un ritmo proporcional a 𝑃. Entre las 6 P.M. y las 7 P.M.
la población se triplica. ¿A qué hora será cien veces mayor que la que había a las 6 P.M.?
R/ 𝑷𝒂𝒔𝒂𝒅𝒂𝒔 𝒍𝒂𝒔 𝟏𝟎 𝒑𝒎

7) La población de una ciudad minera crece a un ritmo proporcional a dicha población. En dos años la
población se ha doblado, y un año más tarde había 10000 habitantes. ¿Cuál era la población inicial?
R/ 3536 personas

8) Un moho crece a un ritmo proporcional a la cantidad presente. Inicialmente había 2 gramos. En dos
días ha pasado a haber 3 gramos.

a) Si 𝑥 = 𝑥(𝑡) es la masa de moho en el instante 𝑡, probar que 𝑥 = 2(3/2)𝑡/2

R/ Demostración

66
 b) Calcular la cantidad al cabo de 10 días.
R/ 15.19

9) La fisión nuclear produce neutrones en una pila atómica a un ritmo proporcional al número de
neutrones presentes en cada momento. Si hay 𝑛0 neutrones inicialmente y hay 𝑛1 y 𝑛2 ,
respectivamente, en los instantes 𝑡1 y 𝑡2 , demostrar que

𝑛1 𝑡2 𝑛2 𝑡1
( ) =( )
𝑛0 𝑛0

R/ Demostración

10) En 1980, el departamento de recursos naturales liberó 1000 ejemplares de una especie de pez en un
lago. En 1987, la población de estos peces en el lago se estimó en 3000. Use la ley de crecimiento
exponencial (Ley de Malthus) para estimar la población de estos peces en el lago en el año 2010.
R/ 𝟏𝟏𝟎𝟖𝟔𝟖 𝒑𝒆𝒄𝒆𝒔

11) En 1970, la población de lagartos en los terrenos del Centro Espacial Kennedy se estimó en 300. En
1980, la población había aumentado hasta un estimado de 1500. Use la ley de crecimiento
exponencial para estimar la población de lagartos en dichos terrenos, en el año 2010.
R/ 𝟏𝟖𝟕𝟓𝟎𝟎 𝒍𝒂𝒈𝒂𝒓𝒕𝒐𝒔

12) Cuando pasa un rayo vertical de luz por un medio transparente, la razón con que crece su intensidad
𝐼 es proporcional a 𝐼(𝑡), donde 𝑡 representa el espesor, en pies, del medio. En agua limpia de mar, la
intensidad a 3 pies debajo de la superficie es 25% de la intensidad inicial 𝐼0 del rayo incidente. ¿Cuál
es la intensidad del rayo a 15 pies debajo de la superficie?
R/

13) Cuando el interés es compuesto continuamente, la cantidad de dinero aumenta con una razón
proporcional a la cantidad presente 𝑆 al tiempo 𝑡, es decir, 𝑑𝑆⁄𝑑𝑡 = 𝑟𝑆, donde 𝑟 es la razón de interés
anual.

a) Calcule la cantidad reunida al final de 5 años cuando se depositan $5000 en una cuenta de
3
ahorro que rinde el 5 4 % de interés anual compuesto continuamente.
R/

b) ¿En cuántos años se habrá duplicado el capital inicial?

R/

67
14) El isótopo radiactivo del plomo Pb-209, decae con una razón proporcional a la cantidad presente al
tiempo 𝑡 y tiene una vida media de 4 horas. Si al principio había 2 gramos de plomo, ¿cuánto tiempo
debe trascurrir para que decaiga 80%?
R/ 𝟗. 𝟐𝟗 𝒉𝒐𝒓𝒂𝒔

15) Resuelva los siguientes problemas

𝑑𝐴
a) El problema de valores iniciales = 𝑘𝐴 , 𝐴(0) = 𝐴0 es el modelo de decaimiento de una
𝑑𝑡
− ln 2
sustancia radiactiva. Demuestre que, en general, la vida media 𝑇 de la sustancia es 𝑇 = 𝑘
R/ Demostración

b) Demuestre que la solución del problema con valores iniciales del ítem a) se puede escribir
como 𝐴(𝑡) = 𝐴0 2−𝑡/𝑇 .
R/ Demostración

c) Si una sustancia radiactiva tiene la vida media 𝑇 dada en el ítem a), ¿cuánto tiempo le tomará
1
a una cantidad 𝐴0 de sustancia decaer a 6 𝐴0 ?
R/

16) La población de una ciudad aumenta con un coeficiente de variación que es proporcional al número
de sus habitantes en cualquier instante 𝑡. Si la población de la ciudad era 30000 en 1960 y 35000 en
1970, ¿cuál será su población en 1990?
R/ 𝟒𝟎𝟖𝟑𝟑 𝒑𝒆𝒓𝒔𝒐𝒏𝒂𝒔

17) Si la mitad de cierta cantidad de radio se desintegra en 1600 años, ¿qué porcentaje de la cantidad
original quedará al cabo de 2400 años? ¿Y de 8000 años?
R/ 𝟑𝟓. 𝟑𝟓% , 𝟑. 𝟏𝟐𝟓%

18) Si la vida media (semivida) de una sustancia radiactiva es de 20 días, ¿cuánto tardará en
desintegrarse el 90% de ella?
R/ 𝟔𝟔. 𝟒𝟒 𝒅í𝒂𝒔

19) El uranio 238 se desintegra a un ritmo proporcional a la cantidad presente. Si hay 𝑥1 y 𝑥2 gramos en
los instantes 𝑡1 y 𝑡2 , probar que la vida media es

(𝑡2 − 𝑡1 ) log 2
𝑥
log (𝑥1 )
2
R/ Demostración

68
20) El radio se descompone a una velocidad proporcional a la cantidad presente. Si la mitad de la cantidad
original desaparece en 1600 años, hallar el porcentaje de pérdida en 100 años?
R/ 𝟒. 𝟐𝟒%

21) Si en un principio se tienen 50 g de una sustancia radiactiva y después de 3 días sólo restan 10 g,
¿qué porcentaje de la cantidad original quedará después de 4 días?
R/ 𝟏𝟏. 𝟕%

22) Si en un principio se tienen 300 g de una sustancia radiactiva y después de 5 años restan 200 g,
¿cuánto tiempo deberá transcurrir para que sólo queden 10 g?
R/ 𝟒𝟏. 𝟗𝟒 𝒂ñ𝒐𝒔

23) Con frecuencia, el fechado por carbono se usa para determinar la edad de un fósil. Por ejemplo,
en una cueva de Sudáfrica se halló un cráneo humano junto con los restos de una hoguera. Los
arqueólogos creen que la edad del cráneo sea igual a la edad de la hoguera. Se ha determinado que
sólo queda 2% de la cantidad original de carbono 14 en los restos de madera en la hoguera. Estime la
edad del cráneo, si la vida media del carbono 14 es de aproximadamente 5600 años.
R/ 𝟑𝟏𝟔𝟎𝟓. 𝟓𝟗 𝒂ñ𝒐𝒔

24) Suponga que el coeficiente de variación instantánea con que se desintegra un núcleo radiactivo es
proporcional al número de tales núcleos, presentes en una muestra dada. En una cierta muestra, el
10% del número original de núcleos radiactivos ha sufrido desintegración en un periodo de 200 años.

a) ¿Qué porcentaje de los núcleos radiactivos originales quedará al cabo de 1000 años?
R/ 𝟓𝟗. 𝟎𝟓%

b) ¿En cuántos años quedará solamente un cuarto del número original?

R/ 𝟐𝟔𝟑𝟏. 𝟓𝟐 𝒂ñ𝒐𝒔

2.1.2. Ley de enfriamiento o calentamiento de Newton

1) Un termómetro se lleva de una habitación hasta el ambiente exterior, donde la temperatura del aire
es 5℉. Después de 2 minutos, el termómetro indica 56℉ y después de 6 minutos indica 28℉. ¿Cuál
era la temperatura inicial de la habitación?
R/ 𝟖𝟎. 𝟗𝟒 ℉

2) Una pequeña barra de metal, cuya temperatura inicial era de 22℃, se deja caer en un gran tanque de
agua hirviendo. ¿Cuánto tiempo tardará la barra en alcanzar los 85℃, si se sabe que su temperatura
aumentó 3℃ en 1 segundo? ¿Cuánto tiempo tardará en alcanzar 98℃?
R/ 42.03 seg ; 93.41 seg

69
3) Un termómetro que indica 71℉ se coloca en un horno precalentado a una temperatura constante. A
través de una ventana de vidrio en la puerta del horno, un observador registra que el termómetro lee
120℉ después de 0,5 minutos y 145℉ después de 1 minuto. ¿Cuál es la temperatura del horno?
R/ 171.04 ℉

4) Un cuerpo de temperatura desconocida se coloca en un frigorífico que se mantiene a temperatura

constante de 0℉. Tras 15 minutos, el cuerpo está a 30℉, y después de 30 minutos está ya a 15℉.
¿Cuál era su temperatura inicial?
R/ 𝟔𝟎℉

5) Una olla de sopa, inicialmente hirviendo, se enfría en aire a 0℃, y a los 30 minutos está a una
temperatura de 20℃. ¿Cuánto se enfriará en los siguientes 30 minutos?
R/ 16 ℃

6) Por razones obvias, la sala de disección de un forense se mantiene fría a una temperatura constante
de 5℃ (41℉). Mientras se encontraba realizando la autopsia de la víctima de un asesinato, el propio
forense es asesinado, y el cuerpo de la víctima robado. A las 10 A.M. el ayudante del forense descubre
su cadáver a una temperatura de 23℃. A mediodía, su temperatura es de 18,5℃. Suponiendo que el
forense tenía en vida la temperatura normal de 37℃ (98,6℉), ¿a qué hora fue asesinado?
R/ 𝟔 𝒂. 𝒎

7) Dos grandes tanques 𝐴 y 𝐵 del mismo tamaño se llenan con fluidos diferentes. Los fluidos en los
tanques 𝐴 y 𝐵 se mantienen a 0℃ y a 100℃, respectivamente. Una pequeña barra de metal, cuya
temperatura inicial es de 100℃, se sumerge dentro del tanque 𝐴. Después de 1 minuto la temperatura
de la barra es de 90℃. Después de 2 minutos se saca la barra e inmediatamente se trasfiere al otro
tanque. Después de 1 minuto en el tanque 𝐵 la temperatura se eleva 10℃. ¿Cuánto tiempo medido
desde el comienzo de todo el proceso, le tomará a la barra alcanzar los 99.9℃?
R/ 𝟗 𝒎𝒊𝒏

8) Un contenedor de líquido caliente se coloca en un congelador que se mantiene a una temperatura

constante de 20℉. La temperatura inicial del líquido es 160℉. Después de 5 minutos, la temperatura
del líquido es 60℉. ¿Cuánto tiempo se necesitará para que su temperatura disminuya a 30℉?
R/

9) Un cadáver se encontró dentro de un cuarto cerrado en una casa donde la temperatura era constante
a 70℉. Al tiempo del descubrimiento la temperatura del corazón del cadáver se determinó de 85℉.
Una hora después una segunda medición mostró que la temperatura del corazón era de 80℉.
Suponga que el tiempo de la muerte corresponde a 𝑡 = 0 y que la temperatura del corazón en ese
momento era de 98.6℉. Determine ¿cuántas horas pasaron antes de que se encontrara el cadáver?
R/

70
10) Cuando un objeto se extrae del horno y se coloca en un entorno con una temperatura constante
de 80℉, la temperatura en el centro es 1500℉. Una hora después de extraerlo, la temperatura del
centro es 1120℉. Encontrar la temperatura del centro 5 horas después de extraer el objeto del horno.
R/ 𝟑𝟕𝟗. 𝟐𝟑℉

11) Según la ley de Newton, la velocidad a que se enfría una sustancia al aire libre es proporcional a
la diferencia entre la temperatura de la sustancia y la del aire. Si la temperatura del aire es 30℃ y la
sustancia se enfría de 100 a 70℃ en 15 minutos, ¿cuándo será 40℃ la temperatura de la sustancia?
R/

12) Un cuerpo se enfría de 60℃ a 50℃ en 15 min, encontrándose sumergido en aire que se mantiene
a 30℃. ¿Cuánto tiempo tardará este cuerpo en enfriarse de 100℃ a 80℃ en aire que se mantiene
a 50℃ ?
R/

2.2. Ecuaciones diferenciales lineales

1. La ecuación diferencial que describe la velocidad 𝑣 de una masa que cae sujeta a una resistencia del
aire proporcional a la velocidad instantánea es

𝑑𝑣
= 𝑚𝑔 − 𝑘𝑣 𝑚
𝑑𝑡
donde 𝑘 > 0 es una constante de proporcionalidad. La dirección positiva se toma hacia abajo.

a) Resuelva la ecuación sujeta a la condición inicial 𝑣(0) = 𝑣0 .

𝒌
𝒈𝒎 𝒈𝒎
R/ 𝒗(𝒕) = + (𝒗𝟎 − ) 𝒆− 𝒎 𝒕
𝒌 𝒌

b) Determinar la velocidad límite o terminal de la masa.

𝒈𝒎
R/ 𝒌

c) Si la distancia 𝑠, medida desde el punto en el que se suelta la masa se relaciona con la velocidad
𝑑𝑠
𝑣 por = 𝑣(𝑡), determine una expresión explicita para 𝑠(𝑡), si 𝑠(0) = 0.
𝑑𝑡
𝒌
𝒎 𝒈𝒎 𝒗𝟎 𝒎 𝒈𝒎𝟐
R/ 𝒔(𝒕) = (𝒈𝒕 + 𝒗𝟎 − )−( − ) 𝒆− 𝒎 𝒕
𝒌 𝒌 𝒌 𝒌𝟐

2. Cuando cae una gota de lluvia, ésta se evapora mientras conserva su forma esférica. Si se hacen
suposiciones adicionales de que la rapidez a la que se evapora la gota de lluvia es proporcional a su
área superficial y que se desprecia la resistencia del área, entonces un modelo para la velocidad 𝑣(𝑡)
de la gota de lluvia es

𝑑𝑣 3(𝑘/𝜌)
+ 𝑣=𝑔
𝑑𝑡 (𝑘/𝜌)𝑡 + 𝑟0
71
 Aquí 𝜌 es la densidad del agua, 𝑟0 es el radio de la gota de lluvia en 𝑡 = 0, 𝑘 < 0 es la constante de
proporcionalidad y la dirección hacia abajo se considera positiva. Determine 𝑣(𝑡) si la gota de lluvia
cae a partir del reposo.
𝒈𝝆 𝒌𝒕 𝒌𝒕 −𝟒
R/ 𝒗(𝒕) = ( + 𝒓𝟎 ) [𝟏 + 𝒓𝟒𝟎 ( 𝝆 + 𝒓𝟎 ) ]
𝟒𝒌 𝝆

3. La ecuación diferencial
𝑑𝑃
= (𝑘 cos 𝑡)𝑃
𝑑𝑡

Donde 𝑘 es una constante positiva, es un modelo matemático para una población 𝑃(𝑡) que
experimenta fluctuaciones anuales. Resuelva la ecuación diferencial sujeta a 𝑃(0) = 𝑃0 .
R/ 𝑷(𝒕) = 𝑷𝟎 𝒆𝒌 𝐬𝐢𝐧 𝒕

4. En un modelo del cambio de población de 𝑃(𝑡) de una comunidad, se supone que

𝑑𝑃 𝑑𝐵 𝑑𝐷
= −
𝑑𝑡 𝑑𝑡 𝑑𝑡
donde 𝑑𝐵/𝑑𝑡 y 𝑑𝐷/𝑑𝑡 son las tasas de natalidad y mortalidad, respectivamente.

𝑑𝐵 𝑑𝐷
a) Determine 𝑃(𝑡) si = 𝑘1 𝑃 y = 𝑘2 𝑃.
𝑑𝑡 𝑑𝑡
b) Analice los casos 𝑘1 > 𝑘2 , 𝑘1 = 𝑘2 , 𝑘1 < 𝑘2 .

5. Un modelo que describe la población de una pesquería en la que se cosecha con una razón constante
está dada por

𝑑𝑃
= 𝑘𝑃 − ℎ
𝑑𝑡
donde 𝑘 y ℎ son constantes positivas.

a) Resuelva la ecuación diferencial sujeta a 𝑃(0) = 𝑃0 .

𝒉
R/ 𝑷(𝒕) = 𝑷𝟎 𝒆𝒌𝒕 + 𝒌 (𝟏 − 𝒆𝒌𝒕 )

b) Describa el comportamiento de la población 𝑃(𝑡) conforme pasa el tiempo en los tres casos 𝑃0 >
ℎ/𝑘, 𝑃0 = ℎ/𝑘 y 0 < 𝑃0 < ℎ/𝑘.
c) Utilice los resultados del inciso b) para determinar si la población de peces desaparecerá en un
tiempo finito, es decir, si existe un tiempo 𝑡 > 0 tal que 𝑃(𝑡) = 0. Si la población desaparecerá,
entonces determine en qué tiempo 𝑡 pasará.

6. Un modelo matemático para la razón con la que se propaga una medicina en el torrente sanguíneo
está dado por

72
 𝑑𝑥
= 𝑟 − 𝑘𝑥
𝑑𝑡

donde 𝑟 y 𝑘 son constantes positivas. Sea 𝑥(𝑡) la función que describe la concentración de la medicina
en el torrente sanguíneo al tiempo 𝑡. Resuelva la ecuación diferencial sujeta a 𝑥(0) = 0.
𝒓 𝒓
R/ 𝒙(𝒕) = 𝒌 − 𝒌 𝒆−𝒌𝒕

7. Cuando se considera la falta de memoria, la razón de memorización de un tema está dada por
𝑑𝐴
= 𝑘1 (𝑀 − 𝐴) − 𝑘2 𝐴
𝑑𝑡

donde 𝑘1 > 0, 𝑘2 > 0, 𝐴(𝑡) es la cantidad memorizada al tiempo 𝑡, 𝑀 es la cantidad total a memorizarse
y 𝑀 − 𝐴 es la cantidad que falta por memorizar. Resuelva la ecuación diferencial sujeta a 𝐴(0) = 0.
𝒌𝟏 𝑴
R/ 𝑨(𝒕) = 𝒌 (𝟏 − 𝒆−(𝒌𝟐 +𝒌𝟏 )𝒕 )
𝟐 +𝒌𝟏

8. Resolver la ecuación diferencial (1 + 𝑥 2 )𝑦 ′ + 2𝑥𝑦 = 𝑓(𝑥) ; 𝑦(0) = 0, donde

𝑥, 0≤𝑥<1
𝑓(𝑥) = {
−𝑥, 𝑥≥1
Respuesta:

9. Resuelva el problema de valor inicial

3𝑦
a) 𝑦 ′ + + 2 = 3𝑥 , 𝑦(1) = 1
𝑥
𝟑𝒙𝟐 𝟏𝒙 𝟗
R/ 𝒚 = − + 𝟏𝟎 𝒙−𝟑
𝟓 𝟐

b) 𝑡 3 𝑥 ′ + 3𝑡 2 𝑥 = 𝑡 , 𝑥(2) = 0
𝟏 −𝟏
R/ 𝒙(𝒕) = 𝒕 − 𝟐𝒕−𝟑
𝟐

10. Con frecuencia, la secreción de hormonas en la sangre es una actividad periódica. Si una hormona
se secreta en un ciclo de 24 horas, entonces la razón de cambio del nivel de la hormona en la sangre
se puede representar mediante el problema con valor inicial

𝑑𝑥 𝜋𝑡
= 𝛼 − 𝛽 cos − 𝑘𝑥 , 𝑥(0) = 𝑥0
𝑑𝑡 12

donde 𝑥(𝑡) es la cantidad de la hormona en la sangre en el instante 𝑡, 𝛼 es la razón promedio de

secreción, 𝛽 es la cantidad de variación diaria en la secreción y 𝑘 es una constante positiva que refleja

73
 la razón con la que un cuerpo elimina la hormona de la sangre. Si 𝛼 = 𝛽 = 1, 𝑘 = 2 y 𝑥0 = 10, halle
𝑥(𝑡).
Respuesta:

11. Resuelva las siguientes ecuaciones diferenciales

𝑑𝑦
a) 𝑥 𝑑𝑥 + 4𝑦 = 𝑥 3 − 𝑥

b) 𝑥 2 𝑦 ′ + 𝑥(𝑥 + 2)𝑦 = 𝑒 𝑥

2.2.1. Modelo de mezclas

1) Un gran tanque de 500 galones está lleno de agua pura. Le entra salmuera que tiene 2 lb de sal por
galón a razón de 5 gal/min. La solución bien mezclada sale del tanque con la misma razón. Determine
la cantidad 𝐴(𝑡) de libras de sal que hay en el tanque al tiempo 𝑡.
Respuesta:

2) En el problema 1, ¿cuál es la concentración 𝑐(𝑡) de sal en el tanque al tiempo 𝑡? ¿Y al tiempo 𝑡 = 5

min? ¿Cuál es la concentración en el tanque después de un largo tiempo, es decir, conforme 𝑡 → ∞?
¿Para qué tiempo la concentración de sal en el tanque es igual a la mitad de este valor límite?
Respuesta:

3) Resuelva el problema 1 suponiendo que la solución sale con una razón de 10 gal/min. ¿Cuándo se
vacía el tanque?
Respuesta:

4) Una solución salina entra a una razón constante de 4 litros/minuto en un tanque de gran tamaño que
en un principio contenía 100 litros de agua pura. La solución dentro del tanque se mantiene bien
revuelta y sale del tanque a razón de 3 litros/minuto. Si la concentración de sal en la solución que
entra al tanque es de 0.2 kg/litro, determine la masa de sal (o cantidad de sal) en el tanque después
de 5 minutos. ¿Cuándo llegará la concentración de sal en el tanque a 0.1 kg/litro?
Respuesta:

5) Un tanque contiene 100 DecaL de salmuera obtenida disolviendo 80 kg de sal en agua. Se introduce
en el tanque agua pura a una velocidad de 4 DecaL/min y la mezcla, conservada homogénea mediante
agitación, sale a la misma velocidad, yendo a parar a un segundo tanque que contiene al principio 100
DL de agua pura. Agitando se mantiene homogénea la mezcla que sale de este segundo tanque a la
misma velocidad ya citada. Hallar la cantidad de sal en el segundo tanque al cabo de 1 hora.
4 𝑥
Ayuda: 𝑥 ′ = 4 (5 𝑒 −0.04𝑡 ) − 4 100 𝑝𝑎𝑟𝑎 𝑒𝑙 𝑠𝑒𝑔𝑢𝑛𝑑𝑜 𝑡𝑎𝑛𝑞𝑢𝑒.

Respuesta: 17,4 Kg

74
6) Un tanque contiene 1000 litros de salmuera obtenida disolviendo 60 kg de sal en agua. Se introduce
en el tanque, a una velocidad de 20 L/min, agua salada que contiene 10 kg de sal por litro, y la mezcla,
conservada homogénea mediante agitación, sale a una velocidad de 30 L/min. Hallar la cantidad de
sal en el tanque al cabo de 1 hora.
Respuesta:

7) Un gran tanque de 500 galones está lleno de agua pura. Le entra salmuera que tiene 2 lb de sal por
galón a razón de 5 gal/min. La solución bien mezclada sale del tanque con la misma razón.
Determine la cantidad 𝐴(𝑡) de libras de sal que hay en el tanque al tiempo 𝑡.
𝟏
Respuesta: 𝑨(𝒕) = 𝟏𝟎𝟎𝟎 (𝟏 − 𝒆−𝟏𝟎𝟎𝒕 )

2.3. Soluciones por sustitución

1) Resolver las siguientes ecuaciones diferenciales
𝑑𝑥 𝑥 2 +𝑡√𝑡 2 +𝑥2
a) =
𝑑𝑡 𝑡𝑥
𝑑𝑦 𝑦(𝑙𝑛𝑦−𝑙𝑛𝑥+1)
b) =
𝑑𝑥 𝑥

2) Una ecuación diferencial de la forma

𝑑𝑦
= 𝑃(𝑥)𝑦 2 + 𝑄(𝑥)𝑦 + 𝑅(𝑥)
𝑑𝑥

es una ecuación de Riccati20 generalizada. Si se conoce una solución de la ecuación, digamos 𝑢(𝑥),
1
muestre que la sustitución 𝑦 = 𝑢 + 𝑣 reduce la ecuación a una ecuación lineal en 𝑣.

3) Resolver
𝑑𝑦
a) = (𝑥 − 𝑦 + 5)2
𝑑𝑥
b) Sin utilizar el método de las ecuaciones exactas

(𝑥 + 𝑦 − 1)𝑑𝑥 + (𝑦 − 𝑥 − 5)𝑑𝑦 = 0

2.4. Ecuaciones exactas

1) Resolver las siguientes ecuaciones diferenciales
a) (2𝑦 2 + 2𝑦 + 4𝑥 2 )𝑑𝑥 + (2𝑥𝑦 + 𝑥)𝑑𝑦 = 0
b) (𝑥 2 𝑠𝑒𝑛𝑥 + 4𝑦)𝑑𝑥 + 𝑥𝑑𝑦 = 0

20Jacopo Francesco Riccati (1676-1754), matemático y filósofo italiano. Estudió detalladamente la hidrodinámica sobre la
base de la mecánica newtoniana.
75
 c) (2𝑦 2 + 3𝑥)𝑑𝑥 + 2𝑥𝑦𝑑𝑦 = 0
d) (6𝑥𝑦)𝑑𝑥 + (4𝑦 + 9𝑥 2 )𝑑𝑦 = 0
e) (𝑦 2 + 𝑥𝑦 3 )𝑑𝑥 + (5𝑦 2 − 𝑥𝑦 + 𝑦 3 𝑠𝑒𝑛𝑦)𝑑𝑦 = 0

2) Considere la ecuación

(5𝑥 2 𝑦 + 6𝑥 3 𝑦 2 + 4𝑥𝑦 2 )𝑑𝑥 + (2𝑥 3 + 3𝑥 4 𝑦 + 3𝑥 2 𝑦)𝑑𝑦 = 0

a) Muestre que la ecuación no es exacta

R/ Demostración

b) Multiplique la ecuación por 𝑥 𝑛 𝑦 𝑚 y determine valores de 𝑛 y 𝑚 que hagan exacta a la ecuación

resultante.
R/ 𝒏 = 𝟐 , 𝒎 = 𝟏

c) Use la solución de la ecuación exacta resultante para resolver la ecuación original.

3) Resuelva el problema de valor inicial

(𝑥 2 + 𝑦 2 − 5)𝑑𝑥 = (𝑦 + 𝑥𝑦)𝑑𝑦, 𝑦(0) = 1

3. Ecuaciones diferenciales de orden superior

3.1. Ecuaciones diferenciales lineales homogéneas

1) Resuelva

a) 𝑦 ′′ − 4𝑦 ′ + 5𝑦 = 0
R/ 𝒚 = 𝒆𝟐𝒙 (𝒄𝟏 𝐜𝐨𝐬 𝒙 + 𝒄𝟐 𝐬𝐢𝐧 𝒙)

b) 𝑥 ′′′ − 𝑥 ′′ − 4𝑥 = 0
√𝟕 √𝟕
R/ 𝒙 = 𝒆−𝟏/𝟐𝒕 (𝒄𝟏 𝐜𝐨𝐬 𝒕 + 𝒄𝟐 𝐬𝐢𝐧 𝒕) + 𝒄𝟑 𝒆𝟐𝒕
𝟐 𝟐

𝑑2 𝑦 𝜋 𝜋
c) + 𝑦 = 0, 𝑦 ( 3 ) = 0 , 𝑦′ ( 3 ) = 2
𝑑𝜃2
R/

2) Para determinar una solución general de la ecuación de Cauchy-Euler

76
 𝑑2 𝑦 𝑑𝑦
3𝑥 2 𝑑𝑥 2 + 11𝑥 𝑑𝑥 − 3𝑦 = 0 ; 𝑥 > 0

Complete los siguientes pasos:

a) Use la sustitución 𝑥 = 𝑒 𝑡 para obtener

𝑑2 𝑦 𝑑𝑦
3 2 +8 − 3𝑦 = 0
𝑑𝑡 𝑑𝑡
R/ Demostración

b) Muestre que la ecuación del ítem a) tiene la solución general

𝑦(𝑡) = 𝑐1 𝑒 𝑡/3 + 𝑐2 𝑒 −3𝑡

R/ Demostración

c) Exprese 𝑦 en términos de 𝑥 para obtener

𝑦(𝑥) = 𝑐1 𝑒 1/3 + 𝑐2 𝑥 −3

R/ Demostración

3) Resuelva el problema de valores iniciales

𝑦 ′′′ − 𝑦 ′ = 0; 𝑦(0) = 2; 𝑦 ′ (0) = 3; 𝑦 ′′ (0) = −1

4) Las raíces de una ecuación cúbica auxiliar son 𝑚1 = 4 y 𝑚2 = 𝑚3 = −6. ¿Cuál es la ecuación
diferencial lineal homogénea correspondiente? Analice: ¿su respuesta es única?

5) Dos raíces de una ecuación auxiliar cúbica con coeficientes reales son 𝑚1 = −1/2 y 𝑚2 = 3 + 𝑖. ¿Cuál
es la ecuación diferencial lineal homogénea correspondiente?

6) Encuentre la solución general de la ecuación diferencial

𝑦 (4) = 𝑦

7) Encuentre la solución general de la ecuación diferencial

𝑦 ′′′ = 𝑦

𝟏
√𝟑 √𝟑
R/ 𝒚(𝒙) = 𝒄𝟏 𝒆𝒙 + 𝒆−𝟐𝒙 (𝒄𝟐 𝐜𝐨𝐬 𝒙 + 𝒄𝟑 𝐬𝐢𝐧 𝒙)
𝟐 𝟐

8) Verifique que 𝑦 = 2𝑠𝑒𝑛3𝑡 + 𝑐𝑜𝑠3𝑡 es una solución del problema con valores iniciales

2𝑦 ′′ + 18𝑦 = 0; 𝑦(0) = 1; 𝑦 ′ (0) = 6

77
 R/ Demostración

9) Para ver el efecto de la variación del parámetro 𝑏 en el problema con valores iniciales

𝑦 ′′ + 𝑏𝑦 ′ + 4𝑦 = 0; 𝑦(0) = 1; 𝑦 ′ (0) = 0

resuelva el problema para 𝑏 = 5, 𝑏 = 4 y 𝑏 = 2 y bosqueje las soluciones.

10) Muestre que

𝑥(𝑡) = 𝑐1 𝑐𝑜𝑠𝑤𝑡 + 𝑐2 𝑠𝑒𝑛𝑤𝑡

es una solución de la ecuación diferencial

𝑥 ′′ + 𝑤 2 𝑥 = 0
Muestre además que una solución que satisfaga las condiciones iniciales 𝑥(𝑡0 ) = 𝑥0 y 𝑥 ′ (0) = 𝑥1
es
𝑥
𝑥(𝑡) = 𝑥0 𝑐𝑜𝑠[𝑤(𝑡 − 𝑡0 )] + 1 𝑠𝑒𝑛[𝑤(𝑡 − 𝑡0 )]
𝑤

R/ Demostración

3.1.1. Sistema Masa-Resorte: Movimiento con/sin

amortiguamiento

1) Una masa que pesa 4 libras se une a un resorte cuya constante es 14 lb/pie. ¿Cuál es el periodo del
movimiento armónico simple?

2) Una masa que pesa 36 libras alarga 0.18 pies un resorte. Al inicio la masa se libera desde un punto
que está 6 pulgadas arriba de la posición de equilibrio con una velocidad descendente de 4 pies/s.
d) Encuentre la ecuación del movimiento.
e) ¿Cuál es la amplitud y el periodo del movimiento?
f) ¿En qué momento la masa pasa por la posición de equilibrio con dirección hacia abajo por tercera
vez?
g) ¿En qué instantes la masa alcanza sus desplazamientos extremos en cualquier lado de la posición
de equilibrio?
h) ¿Cuál es la posición de la masa en 𝑡 = 4 s?
i) ¿Cuál es la velocidad instantánea en 𝑡 = 5 s?
j) ¿Cuál es la aceleración en 𝑡 = 3 s?
k) ¿Cuál es la velocidad instantánea en los instantes cuando la masa pasa por la posición de
equilibrio?
l) ¿En qué instantes la masa está 4 pulgadas debajo de la posición de equilibrio?
m) ¿En qué instantes la masa está 6 pulgadas debajo de la posición de equilibrio apuntando en
dirección hacia arriba?
78
3) Después de que una masa de 20 libras se sujeta a un resorte de 10 pies, éste llega a medir 12 pies.
Se retira la masa y se sustituye con una de 18 libras. Luego se coloca al sistema en un medio que
ofrece una fuerza de amortiguamiento igual a la velocidad instantánea.

a) Encuentre la ecuación del movimiento si la masa se libera inicialmente desde el reposo de un

punto situado 2 pies arriba de la posición de equilibrio.
b) Calcule los tiempos en los que la masa pasa por la posición de equilibrio con dirección hacia abajo.

4) Una masa que pesa 20 libras produce un alargamiento de 4 pies en un resorte. La masa se une a un
dispositivo amortiguador que ofrece una fuerza de amortiguamiento igual a 𝜇 veces la velocidad
instantánea. Determine los valores de 𝜇 por lo que el movimiento posterior sea
a) Sobre amortiguado
b) Críticamente amortiguado
c) Subamortiguado

5) Una masa que pesa 8 libras se une a un resorte. Cuando se pone en movimiento, el sistema masa-
resorte exhibe movimiento armónico simple. Determine la ecuación del movimiento si la constante del
resorte es 1 lb/pie y la masa se libera inicialmente desde un punto 6 pulgadas debajo de la posición
de equilibrio, con una velocidad descendente de 1.5 pies/seg.

3.2. Ecuaciones diferenciales no homogéneas

1) Resuelva las siguientes ecuaciones diferenciales
a) 𝑦 ′′ − 10𝑦 ′ + 25𝑦 = 30𝑥 + 3
𝟔 𝟑
Respuesta: 𝒚(𝒙) = 𝒄𝟏 𝒆𝟓𝒙 + 𝒄𝟐 𝒙𝒆𝟓𝒙 + 𝟓 𝒙 + 𝟓

b) 𝑦 ′′ + 𝑦 = 2𝑥𝑠𝑒𝑛𝑥
c) 𝑦 ′′ + 2𝑦 ′ + 𝑦 = 𝑠𝑒𝑛𝑥 + 3𝑐𝑜𝑠2𝑥
d) 𝑦 (4) − 𝑦 ′′ = 4𝑥 + 2𝑥𝑒 −𝑥
e) 𝑦 ′′ − 6𝑦 ′ + 10𝑦 = 𝑒 3𝑥 − 𝑥
f) 𝑦 ′′′ − 𝑦 ′′ − 4𝑦 ′ + 4𝑦 = 5 − 𝑒 𝑥 + 𝑒 2𝑥

2) Considere la ecuación diferencial 𝑎𝑦 ′′ + 𝑏𝑦 ′ + 𝑐𝑦 = 𝑒 𝑘𝑡 , donde 𝑎, 𝑏, 𝑐 y 𝑘 son constantes. La ecuación

auxiliar de la ecuación homogénea asociada es 𝑎𝑚2 + 𝑏𝑚 + 𝑐 = 0.
a) Si 𝑘 no es una raíz de la ecuación auxiliar, muestre que se puede encontrar una solución particular
de la forma 𝑦𝑝 = 𝐴𝑒 𝑘𝑡 , donde 𝐴 = 1/(𝑎𝑘 2 + 𝑏𝑘 + 𝑐).
R/ Demostración

79
 b) Si 𝑘 es una raíz de la ecuación auxiliar de multiplicidad uno, muestre que se puede encontrar una
solución particular de la forma 𝑦𝑝 = 𝐴𝑥𝑒 𝑘𝑡 , donde 𝐴 = 1/(2𝑎𝑘 + 𝑏). Explique cómo se sabe que
𝑘 ≠ −𝑏/2𝑎.
R/ Demostración

c) Si 𝑘 es una raíz de la ecuación auxiliar de multiplicidad dos, muestre que se puede encontrar una
solución particular de la forma 𝑦𝑝 = 𝐴𝑥 2 𝑒 𝑘𝑡 , donde 𝐴 = 1/(2𝑎).
R/ Demostración

3) Resuelva la siguiente ecuación diferencial

𝑦 ′′ − 5𝑦 ′ + 6𝑦 = 𝑥𝑒 𝑥

4) Determinar la forma de una solución particular a:

𝑦 ′′ + 2𝑦 ′ − 3𝑦 = 𝑓(𝑡)
Si 𝑓(𝑡) es igual a:
a) 2𝑡𝑒 𝑡 𝑠𝑒𝑛𝑡
b) 𝑡 2 𝑐𝑜𝑠3𝑡
c) 5𝑒 −𝑡
d) 3𝑡𝑒 𝑡
e) 𝑡2𝑒𝑡
f) 7 cos 2𝑡

3.2.1. Sistema Masa-Resorte: Movimiento forzado

1) Una masa de 1 slug, cuando se une a un resorte, causa en éste un alargamiento de 2 pies. En 𝑡 = 0,
una fuerza externa igual a 𝑓(𝑡) = 8 sin 4𝑡 es aplicada al sistema. Encuentre la ecuación del movimiento
si el medio ofrece una fuerza de amortiguamiento igual a 8 veces la velocidad instantánea.

2) Resuelva los siguientes problemas

a) Muestre que la solución del problema de valores iniciales

𝑑2𝑥
+ 𝜔2 𝑥 = 𝐹0 cos 𝜇𝑡 , 𝑥(0) = 0, 𝑥 ′ (0) = 0
𝑑𝑡 2
0 𝐹
es 𝑥(𝑡) = 𝜔2−𝜇 2
(cos 𝜇𝑡 − cos 𝜔𝑡)

R/ Demostración
0 𝐹
b) Evalúe lim 𝜔2 −𝜇 2
(cos 𝜇𝑡 − cos 𝜔𝑡)
𝜇→𝜔

80
 𝑭𝟎 𝒕 𝐬𝐢𝐧 𝝎𝒕
R/ 𝟐𝝎

c) Muestre que 𝑥(𝑡) del ítem a) se puede escribir de la forma

−2𝐹0 1 1
𝑥(𝑡) = 2 2
sin (𝜇 − 𝜔)𝑡 sin (𝜇 − 𝜔)𝑡
𝜔 −𝜇 2 2
R/ Demostración

4. La transformada de Laplace

4.1. Solución de problemas de valor inicial

1) Encontrar la transformada de Laplace de las siguientes funciones
𝑡, 0 ≤ 𝑡 < 1
a) 𝑉(𝑡) = {
1, 𝑡 ≥ 1

𝑠𝑒𝑛𝑡, 0 ≤ 𝑡 < 𝜋/2

b) 𝑉(𝑡) = {
0, 𝑡 ≥ 𝜋/2

2) Calcular las siguientes transformadas inversas

2𝑠−6
a) ℒ −1 {𝑠2 +9}
Respuesta: 𝟐(𝐜𝐨𝐬 𝟑𝒕 − 𝐬𝐢𝐧 𝟑𝒕)

1
b) ℒ −1 {5𝑠−2}
𝟏
Respuesta: 𝟓 𝒆𝟐/𝟓𝒕

1
c) ℒ −1 {𝑠3 +5𝑠}

𝑠
d) ℒ −1 {𝑠2 +2𝑠−3}

𝑠
e) ℒ −1 {(𝑠−2)(𝑠−3)(𝑠−5)}

81
 𝑠
f) ℒ −1 {(𝑠+2)(𝑠2 +4)}
𝟏
Respuesta: 𝟒 (𝐜𝐨𝐬 𝟐𝒕 + 𝐬𝐢𝐧 𝟐𝒕 − 𝒆−𝟐𝒕 )

1
g) ℒ −1 {𝑠4 −16}
𝟏 𝟏 𝟓 𝟏𝟏
Respuesta: 𝟑𝟐 𝒆𝟐𝒕 − 𝟑𝟐 𝒆−𝟐𝒕 + 𝟒𝟖 𝐜𝐨𝐬 𝟐𝒕 − 𝟗𝟔 𝐬𝐢𝐧 𝟐𝒕

𝑠−1
h) ℒ −1 {2𝑠2 +𝑠+6}

3) Resuelva las siguientes ecuaciones diferenciales

𝑑𝑦
a) − 𝑦 = 1; 𝑦(0) = 0
𝑑𝑡

Respuesta: 𝒚 = −𝟏 + 𝒆𝒕

𝑑𝑦
b) 2 𝑑𝑡 + 𝑦 = 0; 𝑦(0) = −3

Respuesta: 𝒚 = −𝟑𝒆−𝒕/𝟐

c) 𝑦 ′′ − 4𝑦 ′ = 6𝑒 3𝑡 − 3𝑒 −𝑡 , 𝑦(0) = 1, 𝑦 ′ (0) = −1
𝟏𝟓 𝟑 𝟏𝟏
Respuesta: 𝒚 = − 𝟐𝒆𝟑𝒕 − 𝟓 𝒆−𝒕 + 𝟏𝟎 𝒆𝟒𝒕
𝟔

d) 𝑦 ′′ + 𝑦 = 𝑡 2 + 2 , 𝑦(0) = 1 , 𝑦 ′ (0) = −1
Respuesta:

4) Resuelva la siguiente ecuación diferencial

𝑦 ′′ + 𝑦 = 𝑡 ; 𝑦(𝜋) = 0, 𝑦 ′ (𝜋) = 0

5) Resuelva la siguiente ecuación diferencial

𝑦 ′′ + 6𝑦 ′ + 5𝑦 = 12𝑒 𝑡 ; 𝑦(0) = −1 , 𝑦 ′ (0) = 7

Respuesta: 𝒚 = 𝒆𝒕 − 𝒆−𝟓𝒕 − 𝒆−𝒕

6) Encuentre ℒ{𝑓(𝑡)}

a) 𝑓(𝑡) = sin 2𝑡 cos 2𝑡

82
 𝟐
Respuesta: 𝒔𝟐 +𝟏𝟔

b) 𝑓(𝑡) = sin(3𝑡 + 2)
𝟑 𝒔
Respuesta: 𝐜𝐨𝐬 𝟐 (𝒔𝟐 +𝟗) + 𝐬𝐢𝐧 𝟐 (𝒔𝟐 +𝟗)

c) 𝑓(𝑡) = cos 2 𝑡
𝟏 𝒔
Respuesta: 𝟐𝒔 + 𝟐(𝒔𝟐 +𝟒)

7) Utilice las siguientes fórmulas:

𝑠−𝑎
ℒ −1 { } = 𝑒 𝑎𝑡 cos 𝑏𝑡
(𝑠 − 𝑎)2 + 𝑏 2

𝑏
ℒ −1 { } = 𝑒 𝑎𝑡 sin 𝑏𝑡
(𝑠 − 𝑎)2 + 𝑏 2

Para resolver los siguientes problemas de valor inicial

a) 𝑦 ′ + 𝑦 = 𝑒 −3𝑡 cos 2𝑡 , 𝑦(0) = 0

b) 𝑦 ′′ − 2𝑦 ′ + 5𝑦 = 0 , 𝑦(0) = 1 , 𝑦 ′ (0) = 3

8) Determine ℒ −1 {𝐹(𝑠)} si
5
𝑠 2 𝐹(𝑠) − 4𝐹(𝑠) =
𝑠+1
𝟓 𝟓 𝟓
Respuesta: 𝒇(𝒕) = − 𝟑 𝒆−𝒕 + 𝟏𝟐 𝒆𝟐𝒕 + 𝟒 𝒆−𝟐𝒕

9) Determine ℒ −1 {𝐹(𝑠)} si
𝑠2 + 4
𝑠 2 𝐹(𝑠) + 𝑠𝐹(𝑠) − 6𝐹(𝑠) =
𝑠2 + 𝑠

4.2. Transformada de Laplace de funciones periódicas

1) Encontrar la transformada de Laplace de las siguientes funciones periódicas con periodo 𝑇 dado.

1, 0 ≤ 𝑡 < 𝑎
a) 𝑓(𝑡) = { 𝑇 = 2𝑎
−1, 𝑎 ≤ 𝑡 < 2𝑎

𝟏−𝒆−𝒂𝒔
Respuesta: 𝒔(𝟏+𝒆−𝒂𝒔 )

83
 𝑡, 0 ≤ 𝑡 < 1
b) 𝑓(𝑡) = { 𝑇=2
2 − 𝑡, 1 ≤ 𝑡 < 2

𝑠𝑒𝑛𝑡, 0 ≤ 𝑡 < 𝜋
c) 𝑓(𝑡) = { 𝑇 = 2𝜋
0, 𝜋 ≤ 𝑡 < 2𝜋
𝑎
d) 𝑓(𝑡) = 𝑏 𝑡 0<𝑡≤𝑏 𝑇=𝑏

2) Encontrar la transformada de Laplace de la función “rectificación de onda completa de sin 𝑡”

4.3. Traslación eje 𝒕

1) Calcule

a) ℒ{(3𝑡 + 1)𝑈(𝑡 − 1)}

2
(1+𝑒 −2𝑠 )
b) ℒ −1 { }
𝑠+2

2) Escriba cada función en términos de la función escalón unitario y encuentre su transformada de

Laplace usando dos formas:

 Por medio de la definición de la transformada de Laplace

 Por medio de la función escalón unitario

1, 0 ≤ 𝑡 < 4
a) 𝑓(𝑡) = {0, 4 ≤ 𝑡 < 5
1, 𝑡≥5

𝑡, 0 ≤ 𝑡 < 2
b) 𝑓(𝑡) = {
0, 𝑡≥2

84
3) Use la transformada de Laplace para resolver el problema de valor de inicial

𝑡, 0 ≤ 𝑡 < 1
a) 𝑦 ′ + 2𝑦 = 𝑓(𝑡),𝑦(0) = 0, donde 𝑓(𝑡) = {
0, 𝑡≥1
′′ ′ ′ (0)
b) 𝑦 − 5𝑦 + 6𝑦 = 𝑈(𝑡 − 1), 𝑦(0) = 0, 𝑦 =1

4.4. Más sobre la transformada de Laplace

1) Calcule ℒ −1 {𝐹(𝑠)} si:

𝑠+2
a) 𝐹(𝑠) = ln (𝑠−5)
1
b) 𝐹(𝑠) = 𝑎𝑟𝑐𝑡𝑎𝑛 (𝑠 )

Sugerencia: Utilice la derivada de una transformada.

1
2) Sea 𝐹(𝑠) = 𝑠(𝑠2 +1) , calcule la transformada de Laplace inversa de 𝐹(𝑠)
a) Usando fracciones parciales
b) Usando el teorema de convolución

3) Determine la transformada de Laplace de las siguientes funciones

𝑡
a) 𝑓(𝑡) = ∫0 (𝑡 − 𝑣)𝑒 3𝑣 𝑑𝑣
𝑡
b) 𝑓(𝑡) = ∫0 𝑒 𝑣 𝑠𝑒𝑛(𝑡 − 𝑣) 𝑑𝑣

4) Resuelva las siguientes ecuaciones íntegro-diferencial para 𝑦(𝑡)

𝑡
a) 𝑦(𝑡) + 3 ∫0 𝑦(𝑣)𝑠𝑒𝑛(𝑡 − 𝑣)𝑑𝑣 = 𝑡
𝑡
b) 𝑦(𝑡) + ∫0 (𝑡 − 𝑣)2 𝑦(𝑣)𝑑𝑣 = 𝑡 3 + 3
𝑡
c) 𝑦 ′(𝑡) + 𝑦(𝑡) − ∫0 𝑦(𝑣)𝑠𝑒𝑛(𝑡 − 𝑣)𝑑𝑣 = −𝑠𝑒𝑛𝑡 ; 𝑦(0) = 1
𝑡
d) 𝑦 ′ (𝑡) − 2 ∫0 𝑒 𝑡−𝑣 𝑦(𝑣)𝑑𝑣 = 𝑡 ; 𝑦(0) = 2

5) Calcular

a) ℒ{𝑡 2 𝑐𝑜𝑠𝑡}

85
 b) ℒ{𝑡𝑒 2𝑡 𝑠𝑒𝑛6𝑡}

6) Resuelva el siguiente problema de valor inicial

𝑡𝑦 ′′ − 𝑦 ′ = 2𝑡 2 ; 𝑦(0) = 0

7) Dada la función 𝑓(𝑡) = 𝑡𝑒 2𝑡

Calcule su transformada de Laplace usando las dos fórmulas siguientes:

ℒ{𝑒 𝑎𝑡 𝑓(𝑡)} = 𝐹(𝑠 − 𝑎)

𝑛
𝑑𝑛
ℒ{𝑡 𝑓(𝑡)} = (−1)𝑛 𝐹
𝑑𝑠 𝑛

Verifique que ambos resultados son iguales

86