Escuela de Ciencias Exactas e Ingenierı́a

Programas de Ingenierı́a
Modelación Matemática
Primer semestre 2020

Ecuaciones Separables

!Ejemplo
Según la ley de enfriamiento de Newton, la rapidez con que se enfrı́a (o calienta) un objeto es proporcional a la
diferencia entre su temperatura y la del medio que la rodea, que es la temperatura ambiente. Si T (t) representa la
temperatura del objeto en cualquier instante t y Tm la temperatura ambiente ¿Cuál es la ecuación diferencial que
modela un problema de enfriamiento o calentamiento?.

Se calienta agua a la temperatura del punto de ebullición de 100◦ C. El agua se retira, y se lleva a una habitación
la cual está a una temperatura constante de 50◦ C. Luego de 4 minutos la temperatura del agua es 80◦ C ¿Cuál es
la temperatura del agua en cualquier minuto t? ¿Cuál es la temperatura del agua después de 8 minutos? ¿Cuándo
la temperatura será de 60◦ C? ¿Cuál es la temperatura después de un largo tiempo?

• Variables
– Tiempo: t Variable Independiente.
– T (t): (Temperatura del objeto en un instante de tiempo t) variable dependiente

• Ecuación diferencial que modela un problema de enfriamiento (o calentamiento)

dT
– = k(T − Tm )
dt
• Ecuación diferencial que modela este problema
dT
– = k(T − 50)
dt
• Condiciones iniciales
– Temperatura inicial T (0) = 100.
– Después de 4 minutos : T (4) = 80

• Soporte teórico

Una ecuación diferencial de primer orden es separable si es de la forma

dy
= g(x)h(y)
dx
1
Si se multiplica la ecuación anterior por y éste se renombra p(y), se obtiene:
h(y)

dy
p(y) = g(x)
dx
Separando variables e integrando la igualdad se obtiene,

Z Z
p(y) · dy = g(x) · dx

1 Material recopilado, diseñado y/o editado

por Bibiana Patiño y Jhon Alex Alvarez.
 Y de aquı́,

H(y) = G(x) + C

1
Donde H(y) y G(x) son antiderivadas de p(y) = y g(x) respectivamente. Este último resultado expresa la
h(y)
solución (casi siempre de manera explı́cita) de la ecuación diferencial inicial.

!Solución
Ecuación diferencial
dT
= k(T − 50)
dt

Solución de la ecuación diferencial:

Z Z
dT
= kdt
T − 50

Por lo tanto:
ln(T − 50) = kt + D

Despejando T se obtiene:

T = Cekt + 50

Utilizamos las condiciones iniciales T (0) = 100 y T (4) = 80 para encontrar el valor de las constantes k y C.

Si T (0) = 100 entonces 100 = Cek(0) + 50 ası́ que C = 50

Si T (4) = 80 entonces 80 = 50ek(4) + 50 de donde k = −0, 12771.

De esta manera la temperatura en cualquier instante de tiempo t está dada por la función

T (t) = 50e−0,12771t + 50

¿Cuál es la temperatura del agua después de 8 minutos?

Respuesta: La temperatura después de 8 minutos es 68◦ C

¿Cuándo la temperatura será de 60◦ C?

Respuesta: La temperatura será de 60◦ C aproximadamente a los 12.6 minutos.

¿Cuál es la temperatura después de una largo tiempo?

Respuesta: Después de un largo tiempo la temperatura se mantiene en 50◦ C que es la temperatura ambiente.

En Geogebra con la ayuda del comando ResuelveEDO(Ecuación, Variable dependiente, Variable inde-
pendiente, Punto(s) de f ) en CAS se obtiene la solucion.

Material recopilado, diseñado y/o editado 2

por Bibiana Patiño y Jhon Alex Alvarez.
Comparación de la solución real con la solución obtenida con el campo
de direcciones y métodos numéricos
Para comparar la solución real con la arrojada con el campo de direccioens y con los métodos numéricos se tomará
la ecuación diferencial con el valor obtenido de k = −0, 12771. La ecuaión diferencial del ejemplo es
dT
= −0, 12771(T − 50)
dt
T (0) = 100

!Campo de Direcciones

3 Material recopilado, diseñado y/o editado

por Bibiana Patiño y Jhon Alex Alvarez.
Se observa que la solución real coincide con la tendencia del campo de direcciones.

!Método de Euler

!Método de Runge-Kutta orden 4

Material recopilado, diseñado y/o editado 4

por Bibiana Patiño y Jhon Alex Alvarez.
 Ejercicios
1. La cantidad de individuos de una comunidad se incrementa con una razón proporcional a la cantidad de
personas que tiene en cualquier instante de tiempo t. Si la población se duplico en tres años, ¿En cuánto
tiempo se triplicará y cuadruplicará?

2. Un cultivo tiene inicialmente una cantidad P0 de microbios. Para t = 1 hora, el número de microbios medio
es 32 P0 . Si la rapidez de multiplicación es proporcional al número de microbios presentes, determine el tiempo
necesario para que el número de microbios se triplique.

3. En un parque de exhibición animal se ha observado que la cantidad de animales aumenta proporcionalmente

al número actual estos mismos. Si después de 5 años su número se ha duplicado y después de siete años el
número de animales es 576, hallar el número de animales con que se contaba el dı́a de la inauguración del
parque.

4. La población de una comunidad de jaguares en cierto lugar de la selva crece con una tasa proporcional a la
población en cualquier momento. Su población inicial es de 400 y aumenta el 15% en 10 años. ¿Cuál es la
población pasados 40 años?

5. El xe − 125, isótopo radioactivo del xenón, se desintegra con una razón proporcional a la cantidad presente
en cualquier momento y tiene un periodo medio de vida de 17 horas. Si al principio habı́a 1 gramo de xenón,
¿Cuánto tiempo debe transcurrir para que se desintegre el 80%?

6. Cuando t = 0, habı́a 100 miligramos de una sustancia radiactiva. Al cabo de 5 horas, esa cantidad disminuyó
el 2%. Si la razón de desintegración, en cualquier momento, es proporcional a la cantidad de la sustancia
presente, calcule la cantidad que queda después de 5 horas. Calcule además el periodo medio de vida.

7. Cuando pasa un rayo vertical de luz por una sustancia transparente, la razón con la que decrece su cantidad I
es proporcional a I(t), donde t representa el espesor, en metros, del medio. En agua de mar clara, la intensidad
a 5 metros bajo la superficie, es el 25% de la cantidad inicial I0 del rayo incidente ¿Cuál es la intensidad del
rayo a 20 metros bajo la superficie?

8. Si P (t) es la cantidad de dinero en una cuenta de ahorros que paga una tasa de interés anual de r% compuesto
continuamente, entonces la variación de la cantidad de dinero con respecto al tiempo t (en años) es igual al
r% de la cantidad P (t).
dP r
= P
dt 100
Suponga que el interés es de 5% anual, si la cantidad inicial de dinero es de $1000 y que no hay retiros.

(a) ¿Cuánto dinero habrá en la cuenta después de 4 años?

(b) ¿En qué momento tendrá la cuenta $5500?

9. Al apagar un motor su temperatura es de 98◦ C y el medio en que se encuentra se conserva a 21◦ C. Si después
de 15 minutos el motor se ha enfriado a 90◦ C, encuentre:

(a) La temperatura del motor como función del tiempo.

(b) El instante en el cual la temperatura es 30◦ C.

10. Un cuerpo a una temperatura desconocida se coloca en una habitación que se mantiene a una temperatura
constante a 30◦ F. Si después de 10 minutos la temperatura del cuerpo es 0◦ F y después de 20 minutos la
temperatura del cuerpo es de 15◦ F, ¿Cuál es la temperatura inicial del cuerpo?

11. Si una barra metálica, cuya temperatura inicial es de 20◦ C, se deja caer en un recipiente de agua hirviente
¿Cuánto tiempo tardará en alcanzar 90◦ C, si se sabe que su temperatura aumentó 2◦ C en un segundo?
¿Cuánto tiempo tardará en llegar a 98◦ C?

12. Un dı́a lunes del año 2017 a las 7:00 AM, un detective de la ciudad de Santa Marta encuentra el cuerpo de una
persona en un parque público que conserva una temperatura de 26o C. En ese momento la temperatura del
cuerpo es de 28o C y pasada hora y media la temperatura es de 27.5o C. Considere la temperatura del cuerpo
en el momento de la muerte de 37o C. ¿Cuál fué la hora de la muerte ?

5 Material recopilado, diseñado y/o editado

por Bibiana Patiño y Jhon Alex Alvarez.
 13. A las 5 P.M. se encuentra una persona muerta en un cuarto cuya temperatura es de 20o C. Se mide la
temperatura y resulta ser de 25o C y despues de 20 minutos su temperatura es 22o C. Suponga que la tasa de
enfriamiento del cuerpo es proporcional a la diferencia entre la temperatura del cuerpo y la temperatura del
ambiente. Estime el momento de la muerte de esa persona (hora en la cual la temperatura era de 37.5o C)

14. La razón con que se disemina una medicina en el torrente sanguı́neo se describe con la ecuación diferencial

dx
= r − kx
dt

r y k son constantes positivas. La función x(t) describe la concentración del fármaco en sangre en el momento
t. Determine el valor lı́mite de x(t) cuando t → ∞. ¿En cuánto tiempo la concentración es la mitad del valor
lı́mite? Suponga que x(0)=0.

15. La ecuación diferencial dP

dt = (k cos t)P , en que k es una constante positiva, se usa con frecuencia para modelar
una población que sufre de fluctuaciones estacionales anuales. Determine P (t) y grafique la solución. Suponga
que P (0) = 100000

16. Con frecuencia, un cambio radical en la solución de una ecuación diferencial corresponde a un cambio muy
pequeño en la condición inicial o en la ecuación misma. Compare las soluciones de los siguientes problemas
de valor inicial. Haga además un análisis de la solución a partir del campo de direcciones y por medio del
método de Runge-Kutta de orden 4.
dy
(a) = (y − 1)2 y(0) = 1
dx
dy
(b) = (y − 1)2 y(0) = 1.01
dx
dy
(c) = (y − 1)2 + 0.01 y(0) = 1
dx
dy
(d) = (y − 1)2 − 0.01 y(0) = 1
dx
17. El análisis de las ecuaciones diferenciales no lineales a veces empieza omitiendo términos no lineales de la
ecuación o reemplazándolos con términos lineales. La ecuación diferencial que resulta, se llama linealización
de la primera ecuación; por ejemplo la ecuación diferencial no lineal
 
dP P
= rP 1 −
dt K

Donde r y K son constantes positivas, se usa con frecuencia como modelo de una población creciente pero
acotada. Se puede razonar que cuando P se acerca a cero, el término no lineal P 2 es insignificante. Entonces,
una linealización de la anterior ecuación es
dP
= rP
dt
Suponga que r = 0.02 y que K = 300. Compárese la solución de las dos ecuaciones, con el mismo valor inicial
P (0) por ejemplo: P (0) = 0.5, P (0) = 2 y P (0) = 200. Escriba sus observaciones.

18. La tasa de variación de una población de bacterias y(t) varı́a conjuntamente (k = 1) con el número de bacterias
y y 1 − t. Si la población inicial es de 10.

(a) Plantee la ecuación diferencial.

(b) Resuelve y muestra la gráfica de la solución a la ecuación diferencial.
(c) Basado en la gráfica, qué comportamiento presenta la población de bacterias.
(d) ¿Cuál es la población máxima de bacterias que se presenta?
dy
19. Considere la ecuación de crecimiento logı́stico = y(10 − y) junto con la condición inicial y(0) = 1.
dx
(a) Obtenga la solución exacta a este problema de valor inicial.

Material recopilado, diseñado y/o editado 6

por Bibiana Patiño y Jhon Alex Alvarez.
 (b) Solucionar numéricamente el problema de valor inicial anterior por medio de los métodos de Euler y
Runge-Kutta de orden 4. Utilice un tamaño de paso h = 0.1 y obtenga la solucion en el intervalo [0, 2].
(c) Realice un gráfico comparando la solución exacta con las dos soluciones numéricas obtenidas. ¿Cual
solución numérica se aproxima mejor a la solución exacta?, ¿Qué puede decir del comportamiento de las
soluciones cuando t → ∞?
20. En una población de tamaño L se ha esparcido una epidemia por causa de un individuo enfermo que llegó de
un viaje en el extranjero. Virólogos expertos estiman que en las próximas semanas la tasa de propagación de
la epidemia es conjuntamente proporcional al número de personas infectadas y al número de personas sanas
que quedan. Escriba una ecuación diferencial junto a la condición que describen la situación anterior.
21. Para la situación anterior, si la población infectada N es mucho menor que la poblacion total L, la ecuación
diferencial que describe el crecimiento de la epidemia se puede aproximar por:
dN
= KN
dt
donde K > 0. Basado en el enunciado anterior responda verdadero o falso a las siguientes afirmaciones
(respuesta sin justificacion no será tenida en cuenta)
(a) En la ecuación anterior la incógnita es N (t).
(b) La solución general a la ecuación anterior es de la forma N (t) = No e−kt
22. Un zoológico planea llevar un león marino a otra ciudad. El animal irá cubierto durante el viaje con una manta
mojada. En cualquier tiempo t, la manta perderá humedad debido a la evaporación a una razón proporcional
a la cantidad y(t) de agua presente en la manta. Inicialmente, la manta contendrá 40 litros de agua de mar.

(a) Escriba una ecuación diferencial junto a la condición que describen la situación anterior.
(b) Resuelva y grafique la solución a la ecuación diferencial.
23. Cuando se combinan dos productos quı́micos, la razón a la que se forma un nuevo compuesto está regido por
la ecuación diferencial
dy
= k(y 2 − by + c)
dt
Suponga que k = 0.025, b = 130 y que c = 4000. Aquı́ y(t) denota el número de gramos del nuevo compuesto
formado en el tiempo t.
(a) Calcule los puntos de equilibrio.
(b) Traze el campo de direcciones donde se observen los puntos de equilibrios y clasifı́quelos.
(c) Describa el comportamiento de y cuando t tiende a infinito.
(d) Resuelva la ecuación diferencial suponiendo la condición y(0) = 70. Grafı́quela sobre el campo de
direcciones.

7 Material recopilado, diseñado y/o editado

por Bibiana Patiño y Jhon Alex Alvarez.