Aplicaciones de EDO Primer Orden 2023-2

Universidad Nacional Agraria La Molina
FACULTA D DE CIENCI AS
DEPARTAM EN TO ACADÉM ICO DE M ATEM ÁTICA
CURSO: CC2076 CÁLCULO PARA CIENCIAS

CICLO: 2023 - II
4.3.5 Aplicaciones de ecuaciones diferenciales de primer orden
CÁLCULO PARA CIENCIAS 2023 - II - UNALM

4.3.5 Aplicaciones de ecuaciones diferenciales de primer orden.
Trayectorias ortogonales (TO)
Se dice que dos curvas se cortan de manera ortogonal, si las rectas tangentes que
pasan por cada punto de intersección son perpendiculares entre sí.
Una trayectoria ortogonal de una

familia, es una curva que cruza con
cada una de las curvas de la familia
de manera ortogonal.
https://victoryepes.blogs.upv.es/tag/red-de-flujo/

Las curvas C₁: y = f(x), C₂: y = g(x) se cortan (C₁∩C₂ ) de manera
perpendicular en un punto 𝑥0 , 𝑦0 , si LT₁ es perpendicular a LT₂, donde LT₁ y LT₂
son respectivamente las rectas tangentes correspondientes a C₁ y C₂ en el punto
𝑥0 , 𝑦0 . En general, las pendientes m₁ y m₂ de LT₁ y LT₂ respectivamente
𝑑𝑦 1
satisfacen 𝑚1 . 𝑚2 = −1. Esto, es, si 𝑚1 = 𝑥0 se cumple 𝑚2 = − 𝑑𝑦 𝑥0 .
𝑑𝑥
𝑑𝑥
Ejemplo 1
Encuentre las trayectorias

ortogonales de la familia de curva:
𝑥 2 + 4𝑦 2 = 𝑐.

Derivado implícitamente en la ecuación anterior se tiene:
𝑥
𝑥 + 4𝑦𝑦 ′ = 0 → 𝑦 ′ = −
4𝑦
Para encontrar las trayectoria ortogonales (TO) a la familia de las curvas, hacemos el
cambio:
1
𝑦′ por − ′
𝑦
La ecuación diferencial que se tiene es:
1 𝑥
− ′=−
𝑦 4𝑦
Integrando obtenemos las TO.
𝑦 = 𝑐𝑥 4 , 𝑐 ∈ ℝ

Ejemplo 2
Encuentre las trayectorias ortogonales a la familia de curvas:
𝑦 2 + 2𝑐𝑦 + 𝑥 2 = 0,
donde c es un parámetro.
Solución.
Derivado implícitamente en la ecuación, tenemos:
𝑥 2 + 𝑦2 2𝑥𝑦
2𝑥 + 2𝑐𝑦 ′ + 2𝑦𝑦 ′ = 0, como 2𝑐 = − ′
→𝑦 = 2
𝑦 𝑥 − 𝑦2
Para hallar las TO a la familia de curvas, hacemos el cambio:
1 𝑥 2 − 𝑦2
𝑦′ ′
=− ′ ⇒ 𝑦 =−
𝑦 2𝑥𝑦

La ecuación diferencial obtenida es:
𝑥 2 − 𝑦 2 𝑑𝑥 + 2𝑥𝑦𝑑𝑦 = 0.
Esta ecuación se puede resolver como homogénea, Bernoulli o como exacta con F.I.
Como exacta con F.I.
Identificamos:
𝜕𝑃(𝑥,𝑦)
𝑃 𝑥, 𝑦 = 𝑥2 − 𝑦2 = −2𝑦
𝜕𝑦
⟹ 𝜕𝑃(𝑥,𝑦) 𝜕𝑄(𝑥,𝑦)
≠
𝜕𝑦 𝜕𝑥
𝑄 𝑥, 𝑦 = 2𝑥𝑦 𝜕𝑄(𝑥,𝑦)
= 2𝑦
𝜕𝑥
encontramos g(x):
−4𝑦 2
𝑔 𝑥 = =−
2𝑥𝑦 𝑥

El F.I. de la ecuación es:
−2
‫𝑥𝑑 𝑥 ׬‬
𝑢 𝑥 = 𝑒 = 𝑥 −2
Multiplicando por 𝑢 𝑥 se obtiene la ecuación diferencial exacta:
1 − 𝑦 2 𝑥 −2 𝑑𝑥 + 2𝑦𝑥 −1 𝑑𝑦 = 0,
𝑀(𝑥,𝑦) 𝑁(𝑥,𝑦)
hallando 𝐹 𝑥, 𝑦 ,
𝐹 𝑥, 𝑦 = න 1 − 𝑦 2 𝑥 −2 𝑑𝑥 + 𝜑 𝑦 = 𝑥 + 𝑦 2 𝑥 −1 + 𝜑 𝑦 … ∗
Derivando 𝐹 de la ecuación (*) respecto de “y” e igualamos a 𝑁:
𝜕𝐹(𝑥, 𝑦, 𝑐) −1
𝑑𝜑 𝑦
= 2𝑦𝑥 + = 2𝑦𝑥 −1
𝜕𝑦 𝑑𝑦
𝑑𝜑 𝑦 = 0, entonces, 𝜑 𝑦 = 𝑐,

La familia de las TO es:
𝑥 2 −𝑐𝑥 + 𝑦 2 = 𝟎, 𝒄∈ℝ
Ejemplo 3
Encuentre las trayectorias ortogonales a la familia de curvas:
(𝑥 2 + 𝑦 2 )2 = 𝑐(2𝑥 2 + 𝑦 2 ),
donde c es parámetro.
Solución.
Derivado implícitamente, tenemos:
(𝑥 2 + 𝑦 2 )2
2 𝑥2 + 𝑦2 2𝑥 + 2𝑦𝑦 ′ = 𝑐 4𝑥 + 2𝑦𝑦 ′ , como 𝑐 =
2𝑥 2 + 𝑦 2
Despejado 𝑦 ′ se obtiene

− 2𝑥 3
𝑦′ = 2 . . . (1)
3𝑥 𝑦 + 𝑦 3
Para hallar las TO a la familia de curvas, hacemos el cambio:
1 3𝑥 2𝑦 + 𝑦3
𝑦 ′ por − ′ ⇒ 𝑦′= ,
𝑦 2𝑥 3
de donde se obtiene la ecuación diferencial

3 1
𝑦′ − 𝑦 = 3 𝑦3 … (2)
2𝑥 2𝑥
En este caso la ecuación (2) es una ecuación de Bernoulli ( 𝑘 = 3 ) (también se

puede resolver como una ecuación diferencial homogénea y exacta con F.I.)
𝑑𝑧 𝑑𝑦
Para resolverla hacemos cambio de variable 𝑧 = 𝑦1−(3) = 𝑦 −2 ⟹ = −2𝑦 −3 ,
𝑑𝑥 𝑑𝑥

𝑑𝑧
reemplazamos 𝑧 y en (2), luego simplificando la ecuación se tiene:
𝑑𝑥
3 1
𝑧′ + 𝑧 = − (ecuación lineal) … (3)
𝑥 𝑥3
3
‫𝑥𝑑𝑥׬‬
factor integrante: 𝑢 𝑥 = 𝑒 = 𝑥3;
multiplicando por 𝑢(𝑥) a la ecuación (3) se tiene
𝑑𝑧
𝑥3 + 3𝑥 2 𝑧 = −1; entonces ‫ 𝑥 𝑑 ׬‬3 𝑧 = − ‫ 𝑥𝑑 ׬‬+ 𝑐,
𝑑𝑥
𝑑
𝑥3𝑧
𝑑𝑥
resolviendo y reemplazando 𝑧 se obtiene la solución general, que son las T.O. a la

familia de curvas dadas
𝑥 3 𝑦 −2 = −𝑥 + 𝑐, 𝑐 ∈ ℝ. ó 𝑥 3 + 𝑥𝑦 2 = 𝑐 𝑦 2 , 𝑐 ∈ ℝ.

Ley de enfriamiento de Newton
La Ley de enfriamiento de Newton afirma que la rapidez de enfriamiento de un

objeto es directamente proporcional a la diferencia de las temperaturas entre el
objeto y su alrededor (medio ambiente), es decir;
𝑑𝑇
= −𝑘 𝑇 − 𝑇𝑎 ; 𝑘 > 0.
𝑑𝑡
Donde:
T(t): temperatura del objeto en el tiempo t (en ºC, ºF, ...)

Ta, : temperatura del medio que rodea al objeto (en ºC, ºF, …)
k : constante de proporcionalidad.
t : tiempo (s, min, ...)
Para el tiempo inicial 𝑡0 = 0 , se tiene un problema de valor inicial:
𝑑𝑇
= −𝑘 𝑇 − 𝑇𝑎 , 𝑇 𝑡0 = 𝑇0
𝑑𝑡
La solución general es:
𝑇 = 𝑇𝑎 + 𝑐𝑒 −𝑘𝑡 ; 𝑘 > 0, 𝑐 ∈ ℝ,
y la solución particular es:
𝑇 = 𝑇𝑎 + (𝑇0 − 𝑇𝑎 )𝑒 −𝑘𝑡 ; 𝑘 > 0,

Ejemplo 1.
Cierta habitación tiene una temperatura constante de 60ºC. Si un objeto se coloca

dentro de la habitación, se observa que se enfría de 350ºC hasta 150ºC en 25
minutos. Calcule el tiempo que tiene que transcurrir para enfriar dicho objeto hasta
70ºC.
Solución:
Sean:
𝑇ℎ ∶ temperatura de la habitación (ºC)
𝑇 𝑡 ∶ temperatura del objeto en el instante “t” (t en minutos)
𝑇0 ∶ temperatura inicial del objeto (ºC).
Por la Ley de Enfriamiento de Newton, se tiene la ecuación diferencial:

𝑑𝑇
= −𝑘 𝑇 − 𝑇𝑎 ; 𝑘>0 … (1)
𝑑𝑡

Datos:
𝑇0 = 350 ∘𝐶, 𝑇ℎ = 𝑇𝑎 = 60 ∘𝐶, 𝑇 25 = 150 ∘𝐶
Resolviendo la ecuación diferencial (1) y reemplazando datos, tenemos:
𝑇 𝑡 = 𝑇𝑎 + 𝑇0 − 𝑇𝑎 𝑒 −𝑘𝑡 = 60 + 350 − 60 𝑒 −𝑘𝑡 = 60 + 290𝑒 −𝑘𝑡 … (2)
Para t = 25 minutos se tiene una temperatura de 𝑇 25 = 150 ºC, entonces, se tiene:

1
9 25
𝑒 −𝑘 =
29
Entonces, en la ecuación (2), para una temperatura 𝑇 𝑡 = 70°C, se tiene:
𝑡 1
𝑙𝑛
9 25
𝑡 = 25 29 = 71.946 minutos
70 = 60 + 290 ⟹ 9
29 𝑙𝑛
29

Ejemplo 2.
A las 15 h 00 min pm un cuerpo que tiene una temperatura de 70℃ es depositado en

un lugar donde la temperatura se mantiene a 40℃ . Después de 3 min, la temperatura
del cuerpo ha disminuido a 60℃. Determine:
a) ¿Cuál es la temperatura del cuerpo después de 5 minutos?
b) ¿Cuánto tiempo pasará para que el cuerpo tenga 50℃?
Solución:
a) Si 𝑇 𝑡 es la temperatura del cuerpo en ℃ después de t minutos, entonces la
ecuación diferencial
que modela a 𝑇 𝑡 es
𝑑𝑇(𝑡)
= −𝑘(𝑇 𝑡 − 𝑇𝑎 ),
𝑑𝑡
Datos:

𝑇 5 ≈ 55.262 ℃
𝑡
2 3
b) 𝑇 𝑡 = 40 + 30 = 50 → 𝑡 ≈ 8.1258 𝑚𝑖𝑛
3

Ley de crecimiento y decrecimiento exponencial
(Dinámica de poblaciones: modelo de Malthus)
Si consideramos que la rapidez de crecimiento de una población de individuos es

proporcional a la cantidad de individuos presentes en un instante determinado,
podremos modelar dicho comportamiento de la población a través del problema de
valor inicial.
𝑑𝑃 𝑡
= 𝑘𝑃 𝑡
𝑑𝑡 ,
donde 𝑃 𝑡0 = 𝑃0
P(t): población de individuos en el tiempo t.
k : constante de proporcionalidad.
t : tiempo transcurrido (s, min, días, años, ...), y t₀ el tiempo inicial.
P₀ : Población inicial de individuos.

La solución general es dada por:
𝑃 𝑡 = 𝑐𝑒 𝑘𝑡 ; 𝑐∈ℝ
• Si 𝑘 > 0, el modelo representa una función de crecimiento.
• Si 𝑘 < 0, el modelo representa una función de decrecimiento.

Ejemplo 1. En un estudio socioeconómico anual (Ministerio de la Producción)
vinculados a la cadena de valor de la Trucha Arco Iris en el Perú (2011 − 2016),
siendo la región de Puno el mayor productor. La producción de la Trucha en el año
2013 fue 29 mil toneladas (TM) y en el año 2015 fue de 34 mil TM. Si la
producción crece exponencialmente en cualquier momento 𝑡.
a) Determine una expresión matemática para la cantidad de producción de
Truchas Arco Iris en el Perú.
b) Cuántas TM de Truchas Arco Iris que se obtuvo en el año 2016.
Solución.
𝑑𝑃 𝑡
Modelo de crecimiento exponencial: = 𝑘𝑃 𝑡 , 𝑘 > 0.
𝑑𝑡
Solución general:
𝑃 𝑡 = 𝑐𝑒 𝑘𝑡 ; 𝑐∈ℝ
donde: 𝑃 𝑡 es la producción de Truchas Arco Iris en miles TM en el tiempo 𝑡.

a) De la solución general:𝑃 𝑡 = 𝑃0 𝑒 𝑘𝑡 , 𝑡 ≥ 0; cuando la población inicial es 𝑃0 .
Año Tiempo ⟹ Producción
2011 𝑡 = 0 ⟹ 𝑃 0 = 𝑃0
2013 𝑡 = 2 ⟹ 𝑃 2 = 29
2015 𝑡 = 4 ⟹ 𝑃 4 = 34
2016 𝑡 = 5 ⟹ 𝑃 5 =?
Con los datos se tiene: 29 = 𝑃0 𝑒 2𝑘 𝑦 34 = 𝑃0 𝑒 4𝑘
1
34 2 (29)2
calculamos, 𝑒𝑘 = ; entonces, 𝑃0 = = 24.74
29 34
𝑡
(29)2 34 2
Por lo tanto, 𝑃 𝑡 = =, 𝑡 ≥ 0
34 29

b) Para la producción del año 2016 se tiene, 𝑡 = 5, entonces;
5
(29)2 34 2
𝑃 5 = toneladas
34 29
En el año 2016 la producción de Truchas Arco Iris en el Perú, fue de 36.82
mil TM.
Ejemplo 2
Sea P(t) la población de una ciudad en el instante t. Si se sabe que la población

crece con una razón proporcional a la población presente, además la población
inicial es de 500 habitantes y aumenta en 20% durante 8 años, determine:
a) P(t) de forma explícita.
b) ¿En cuántos años habrá una población de 1000 habitantes?.

Solución:
Datos:
Si:
𝑡 = 0 → 𝑃 0 = 𝑃0 = 500 habitantes … 1
𝑡: 8 años → 𝑃8 = 1.2𝑃0 = 600 habitantes … (2)
a) La solución general de la ecuación diferencial es:

𝑃 𝑡 = 𝑐𝑒 𝑘𝑡 ; 𝑐 ∈ ℝ, 𝑡 en años.
Por la condición inicial, ecuación (1) se tiene:
𝑐 = 𝑃0 = 500,
por la condición de la ecuación (2) tenemos
600 = 500𝑒 8𝑘 … (4)
Entonces,

1/8
𝑘
6
𝑒 = ,
5
Por tanto, ecuación para expresar a la población es:
𝑡/8
6
𝑃 𝑡 = 500 habitantes
5
b) 𝑆𝑖 𝑃(𝑡) = 1000 habitantes, se tiene:
𝑡/8
6
1000 = 500 ,
5
entonces, la población será de 1000 habitantes a los 30.4 años.

Modelo logístico
Este modelo se utiliza para estudiar el crecimiento poblacional. Puede aplicarse

en la producción biológica, pesquera o pecuaria, etc.
La siguiente ecuación diferencial, denominada ecuación logística o modelo

logístico:
𝑝′ 𝑡 = 𝑟𝑝 𝑡 − 𝑚𝑝2 (𝑡)
Donde:
𝑝 𝑡 : representa el número de individuos de la población existentes en el
instante 𝑡;
𝑟𝑝 𝑡 : expresa el crecimiento natural de la población, debido a la reproducción
(la población crece de forma proporcional al número de individuos
presentes);
𝑚𝑝2 (𝑡): representa las limitaciones de los recursos (alimentos. espacio, etc.)
son limitados;
𝑟
𝑚 = , donde 𝑟 es tasa de crecimiento intrínseca y 𝑘 es la capacidad de carga
𝑘
o capacidad sustentable.
Si en el instante inicial 𝑡 = 0, el número de individuos es 𝑝 0 = 𝑃0 , donde

0 < 𝑃0 < 𝑘, entonces la población 𝑝 𝑡 es la solución del siguiente problema
de valor inicial:
𝑑𝑝 𝑡 𝑟 2
= 𝑟𝑝 𝑡 − 𝑝 𝑡
𝑑𝑡 𝑘 … 1
𝑝 0 = 𝑃0
Resolviendo la ecuación diferencial (1) de variables separables
𝑑𝑝 𝑡
න = න𝑟 𝑑𝑡 = 𝑟𝑡 + 𝑐 … 2
1
𝑝 𝑡 1− 𝑝 𝑡
𝑘

Descomponiendo en fracciones parciales la siguiente fracción:
1
1 1 𝑘 1 1
= + = +
1 𝑝 𝑡 1 𝑝 𝑡 𝑘−𝑝 𝑡
𝑝 𝑡 1− 𝑝 𝑡 1− 𝑝 𝑡
𝑘 𝑘
Resolviendo la ecuación (2), la solución general es:

𝑙𝑛 𝑝 𝑡 − 𝑙𝑛 𝑘 − 𝑝 𝑡 = 𝑟𝑡 + 𝑐, 𝑐 ∈ ℝ
Como, 0 < 𝑝 𝑡 < 𝑘 y 𝑟 > 0, simplificando y despejando 𝑝 𝑡 se tiene

la solución general de la ecuación diferencial.
𝐶𝑘𝑒 𝑟𝑡
𝑝 𝑡 = , 𝐶∈ℝ … 3
1 + 𝐶𝑒 𝑟𝑡

Por la condición inicial para 𝑡 = 0, 𝑝 0 = 𝑃0
𝐶𝑘 𝑃0
𝑝 0 = = 𝑃0 ⇒ 𝐶 =
1+𝐶 𝑘 − 𝑃0
Reemplazando la constante en la solución general obtenemos la solución del

problema de valor inicial (solución particular). Esta ecuación representa una
población 𝑝 𝑡 bajo ciertas condiciones y restricciones.
𝑘𝑃0
𝑝 𝑡 =
𝑘 − 𝑃0 𝑒 −𝑟𝑡 + 𝑃0
Ejemplo. Díaz (2010) en un estudio de poblaciones de erizo de mar durante 12

meses de experimento, en el primer mes del experimento (abril) tenía una
población de 213 erizos con una tasa de crecimiento del 21% si el área de
experimento tiene una capacidad sustentable de 840 erizos. Utilizando el modelo
logístico ¿Con cuánto de población de erizos de mar comenzó el estudio?

Solución.
Datos del problema:
tasa de crecimiento 𝑟 = 0.21
capacidad sustentable 𝑘 = 840
población en el primer mes 𝑡 = 1, 𝑝 1 = 213
población inicial 𝑃 0 = 𝑃0 = 𝑥
reemplazando
𝑑𝑝 𝑡 0.21 2
= 0.21𝑝 𝑡 − 𝑝 (𝑡)
𝑑𝑡 840
… (1)
𝑝 0 = 𝑃0

Resolviendo la ecuación diferencial (1) por variables separables se tiene
1 1
න + 𝑑𝑝 𝑡 = න 0.21 𝑑𝑡 + 𝑐
𝑝 𝑡 840 − 𝑝 𝑡
𝑝 𝑡 𝑝 𝑡
𝑙𝑛 = 0.21𝑡 + 𝑐 ⇒ = 𝐶𝑒 0.21𝑡
840 − 𝑝(𝑡) 840 − 𝑝(𝑡)
Despejando 𝑝 𝑡 se tiene:
840𝐶𝑒 0.21𝑡
𝑝 𝑡 = 2
1 + 𝐶𝑒 0.21𝑡

𝑃0
Por la condición inicial; 𝑡 = 0, 𝑝 0 = 𝑃0 , obtenemos 𝐶 = ,
840−𝑃0
remplazando 𝐶 en la ecuación (2) se tiene la solución particular
840𝑃0
𝑝 𝑡 = … 3
840 − 𝑃0 𝑒 −0.21𝑡 + 𝑃0
Del dato:
𝑡 = 1 → 𝑃 1 = 213 erizos
Reemplazando en la ecuación 3 se tiene,
840𝑃0
213 = ⇒ 𝑃0 = 181
840 − 𝑃0 𝑒 −0.21(1) + 𝑃0
El experimento comenzó con una población 𝑃0 = 181 erizos.

Desintegración Radiactiva
Se puede utilizar el mismo modelo de decrecimiento exponencial para modelar el

comportamiento de la desintegración radiactiva, debido a que la velocidad de
descomposición o transmutación de cierto material radiactivo M es proporcional a la
cantidad del elemento presente.
𝑑
𝑀 𝑡 = 𝑘𝑀 𝑡 , 𝑀 𝑡0 = 𝑀0
𝑑𝑡
donde:
M(t) : Cantidad del elemento en un tiempo determinado t.
t : tiempo transcurrido (min, horas, días, años, ...), y t₀ el tiempo inicial.
k : constante de proporcionalidad de la desintegración.
M₀ : cantidad inicial del elemento en el tiempo t₀.

La solución general es: 𝑀 𝑡 = 𝑐𝑒 𝑘𝑡 ; 𝑐 ∈ ℝ,
donde 𝑐 es la cantidad inicial de material radiactivo (𝑐 = M0).
Ahora, en física se llama semivida o vida media de una sustancia radiactiva, al

tiempo necesario para que se desintegre o transmute la mitad del elemento
radioactivo de una cantidad inicial M₀ en el tiempo t₀.
Ejemplo 1.
Una sustancia radiactiva se desintegra a una razón directamente proporcional a su

cantidad presente. Si inicialmente se tiene M₀ mg de esta sustancia radiactiva y
después de 20 minutos queda el 80% y al cabo de una hora sólo queda el 10mg.
Determine la cantidad de esta sustancia radiactiva transcurridos 100 minutos.

Solución:
Datos:
Tasa de crecimiento: Solución general:
𝑑𝑀 𝑡 ⟹ 𝑀 𝑡 = 𝑐𝑒 𝑘𝑡 ; 𝑐∈ℝ
= 𝑘𝑀 𝑡
𝑑𝑡
Si:
𝑡=0 → 𝑀 0 = 𝑀0 𝑚𝑔 … 1
𝑡 = 20 𝑚𝑖𝑛. → 𝑀(20) = 0.8𝑀0 mg … (2)
𝑡 = 60 𝑚𝑖𝑛 → 𝑀(60) = 10 mg … (3)

𝑡 = 100 𝑚𝑖𝑛 → 𝑀(100) = 𝑥 mg … (4)

Por la condición inicial, ecuación (1) se tiene:
𝑐 = 𝑀0
Por la condición de la ecuación (2)
0.8𝑀0 = 𝑀0 𝑒 20𝑘 … (5)
1/20
Entonces: 4
𝑒𝑘 = ,
5
4 𝑡/20
luego M(𝑡) = 𝑀0
5
Por la condición de la ecuación (3)
60/20
4
10 = 𝑀0
5

de donde 𝑀0 = 19.53 𝑚𝑔
Por tanto la cantidad de sustancia radioactiva en el tiempo 𝑡 es:
𝑡/20
4
𝑀 𝑡 = 19.53 mg,
5
Por lo tanto, cuando 𝑡 = 100 minutos tenemos: 6.40 mg de sustancia.

Ejemplo 2.
La velocidad con que se desintegra cierto material radiactivo es proporcional a la

cantidad presente. Se sabe que inicialmente hay 40 miligramos de material y la
mitad de material ha experimentado la desintegración en un periodo de 14 horas.
Determine:
a) Una expresión matemática para la masa del material restante en un tiempo 𝑡.
b) ¿Cuál es la cantidad de material radiactivo después de 10 horas?
c) ¿Después de cuántas horas se perdió el 20 % del material radiactivo?
Solución:
a) Sea 𝑀 𝑡 la cantidad de material radiactivo presente en tiempo t. De acuerdo al
enunciado, se trata del siguiente problema de valor inicial:
𝑑𝑀(𝑡)
= 𝑘𝑀 𝑡 , 𝑀 𝑡0 = 𝑀0
𝑑𝑡
La solución de la ecuación es:
𝑀 𝑡 = 𝑀0 𝑒 𝑘𝑡
donde, 𝑀0 = 40 𝑚𝑔

𝑡
1 14
𝑀 𝑡 = 40 ,𝑡 ≥ 0
2
b) 𝑀 10 =24.38 mg
c) Se pierden 8 miligramos
𝑡 = 4.506
Se pierde el 20 % de la masa original después de 4.506 horas

Ejemplo 3.
La razón con que se desintegra cierto material radiactivo es proporcional a la cantidad
presente. Se sabe que inicialmente hay 60 miligramos de material y después de 5 horas
se observa que el material ha perdido el 25% de material radiactivo. Determine:
a) Una expresión matemática para la cantidad de material radiactivo en un tiempo 𝑡.
b) ¿Cuál es la cantidad de material después de 15 horas?
c) ¿Cuál es la vida media de este material radiactivo?
Solución:
a) Sea 𝑀 𝑡 la cantidad de material radiactivo en el tiempo 𝑡. De acuerdo al
enunciado, se trata del siguiente problema de valor inicial:
𝑑𝑀(𝑡)
= 𝑘𝑀 𝑡 , 𝑀 𝑡0 = 𝑀0
𝑑𝑡
La solución de la ecuación es:
𝑀 𝑡 = 𝑀0 𝑒 𝑘𝑡

donde, 𝑀0 = 60 𝑚𝑔
𝑀 𝑡 = 60𝑒 𝑘𝑡 … (1)
Como se pierde 25%, entonces queda el 75% de la cantidad inicial, entones:
3
(60) = 60𝑒 5𝑘
4
1
3 5
= 𝑒𝑘
4
al reemplazar en (1), se tiene
𝑡
3 5
𝑀 𝑡 = 60 ,𝑡 ≥ 0
4

b) Después de 15 horas la cantidad de material es,
15
3 5
𝑀 15 = 60 = 25.3125 mg
4
Por tanto la cantidad de material después de 15 horas es de 25.3125 mg
aproximadamente
c) Para hallar el tiempo de vida media, planteamos

𝑡
3 5
30 = 60
4
1
ln
𝑡=5 2
3
𝑙𝑛
4
𝑡 = 12.05
La vida media es de 12.05 horas.

Problemas de Mezclas
Consideremos un tanque que contiene L litros de un líquido donde se encuentra 𝑥0

gramos de una sustancia C (sal, azúcar u otras disuelta en agua). Si la mezcla sale
del tanque a una razón de a litros por minuto, de modo que el volumen V en el
tanque permanece constante. En este caso lo que se quiere es determinar la cantidad
de gramos 𝑥 de C que hay en el tanque en un tiempo t cualquiera.
Donde 𝑥 = 𝑓 𝑡 , 𝑡 ≥ 0, y la cantidad inicial de sustancia en el tanque es 𝑓 0 = 𝑥0.

Para encontrar la solución a este problema hacemos un balance de la cantidad de sustancia
en el tanque en un instante t, de la siguiente manera:
La cantidad de sustancia C (g/l) que entra al tanque la podemos expresar en gramos/min,

como:
𝑔𝑟𝑎𝑚𝑜𝑠 𝑙𝑖𝑡𝑟𝑜𝑠 𝑔𝑟𝑎𝑚𝑜𝑠 𝑔𝑟𝑎𝑚𝑜𝑠

Cantidad de 𝐶 que entra al tanque = 𝑎 𝑏 = 𝑎𝑏
𝑙𝑖𝑡𝑟𝑜 𝑚𝑖𝑛𝑢𝑡𝑜 𝑙𝑖𝑡𝑟𝑜 𝑚𝑖𝑛𝑢𝑡𝑜

Similarmente para la salida, la cantidad de C (g/l) que sale del tanque se puede dar en
gramos/min, de la siguiente manera:
𝑔𝑟𝑎𝑚𝑜𝑠 𝑙𝑖𝑡𝑟𝑜𝑠 𝑥 𝑔𝑟𝑎𝑚𝑜𝑠 𝑎𝑥 𝑔𝑟𝑎𝑚𝑜𝑠

Cantidad de 𝐶 que sale del tanque = 𝑎 = .
𝑙𝑖𝑡𝑟𝑜 𝑚𝑖𝑛𝑢𝑡𝑜 𝑉 𝑙𝑖𝑡𝑟𝑜 𝑉 𝑚𝑖𝑛𝑢𝑡𝑜
Entonces, la rapidez de variación de C (gramos/min) es:
𝑑𝑥 𝑎
= 𝑎𝑏 − 𝑥
𝑑𝑡 𝑉
Donde observamos que es una ecuación diferencial de variables separables.

Ejemplo 1
Un tanque mezclador inicialmente contiene 20 kg de sal disueltos en 300 litros de agua.

1
Se sabe que al tanque ingresa agua que contiene 𝑘𝑔 de sal por litro, a razón de 4
5
litros/min. Si la solución está saliendo del tanque con la misma rapidez que ingresa,
encontrar:
a) Una expresión que relacione la cantidad de sal que hay en el tanque con el tiempo 𝑡.
b) Después de cuanto tiempo, desde el inicio, el tanque contiene 40 kg de sal?
Solución.
Sea 𝑥 𝑡 : cantidad de sal que hay en el tanque en el instante 𝑡.

𝑥0 : cantidad inicial de sal que hay en el tanque.
Al inicio, para 𝑡 = 0,
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑉 = 300 𝐿
ቐ
𝑠𝑎𝑙 𝑥0 = 20 𝑘𝑔

𝑚𝑒𝑧𝑐𝑙𝑎 𝑎 = 4 𝐿/𝑚𝑖𝑛
e𝑛𝑡𝑟𝑎𝑑𝑎 = ൞
1
𝑠𝑎𝑙 𝑏= 𝑘𝑔/𝐿
5
𝑚𝑒𝑧𝑐𝑙𝑎 𝑎 = 4 𝐿/𝑚𝑖𝑛
𝑠𝑎𝑙𝑖𝑑𝑎 = 𝑎
𝑠𝑎𝑙 𝑥(𝑡) 𝑘𝑔/𝑚𝑖𝑛
𝑉
Haciendo del balance de sal, planteamos la EDO
𝑑𝑥(𝑡) 𝑎
= 𝑎𝑏 − 𝑥(𝑡)
𝑑𝑡 𝑉
𝑑𝑥(𝑡) 4 4
= − 𝑥(𝑡),
𝑑𝑡 5 300

𝑑𝑥 1 𝑑𝑥 1
= 60 − 𝑥 𝑡 ⇒ = 𝑑𝑡,
𝑑𝑡 75 60 − 𝑥 𝑡 75
Resolvemos la ecuación diferencial se tiene
1
− 𝑙𝑛 60 − 𝑥(𝑡) = 𝑡 + 𝑘 ⇒ 𝑥(𝑡) − 60 = 𝐶𝑒 −𝑡/75 ,
75
𝑥(𝑡) = 60 + 𝐶𝑒 −𝑡/75 ⇒ 𝑥(0) = 60 + 𝐶 ⇒ 20 = 60 + 𝐶 ⇒ 𝐶 = −40.
Entonces,
a) 𝑥(𝑡) = 60 − 40𝑒 −𝑡/75 , 𝑡 ≥ 0
b) para 40 kg, se tiene: 40 = 60 − 40𝑒 −𝑡/75 ⇒ 𝑡 = 75 𝑙𝑛 2 = 51.986 minutos

el tanque contiene 40 kg de sal.

Ejemplo 2.
En un tanque hay 1000 litros de agua contaminada con un compuesto M. Si al inicio

en el tanque 1500 gramos de M. Cuando 𝑡 = 0, ingresa al tanque agua pura a una
razón de 15L/min y la mezcla, bien agitada, sale del tanque a la misma razón.
Determine el tiempo en el que la concentración del contaminante en el tanque sea de
0.02 g/L.
Solución.
Sea 𝑀 𝑡 : cantidad de 𝑀que hay en el tanque en el instante 𝑡.
𝑑𝑴(𝑡)
: rapidez de variación de 𝑀 respecto al tiempo
𝑑𝑡

Datos:
𝑉 : Volumen del tanque
𝑉 = 1000 litros
𝑀0 : cantidad inicial (𝑡 = 0), de 𝑀 que hay en el tanque.
𝑀(0) = 𝑀0 = 1500 g.

𝑚𝑒𝑧𝑐𝑙𝑎 𝑎 = 15 𝐿Τ𝑚 𝑖𝑛
e𝑛𝑡𝑟𝑎𝑑𝑎 = ቐ
compuesto 𝑀 𝑏 = 0 𝑔Τ𝑚 𝑖𝑛
𝑚𝑒𝑧𝑐𝑙𝑎: 𝑎 = 15 𝐿Τ𝑚 𝑖𝑛
𝑠𝑎𝑙𝑖𝑑𝑎 = 𝑎
compuesto 𝑀: 𝑀 𝑡 𝑔Τ𝑚 𝑖𝑛
𝑉
Haciendo del balance de 𝑀, planteamos la EDO
𝑑𝑀(𝑡) 𝑎 𝑑𝑀 15
= 𝑎𝑏 − 𝑀 𝑡 ⇒ =0− 𝑀(𝑡),
𝑑𝑡 𝑉 𝑑𝑡 1000
Luego se tiene
𝑑𝑀 3
=− 𝑑𝑡 (1)
𝑀 𝑡 200

Integrando en la ecuación (1):
3𝑡
3 −
𝑙𝑛 𝑀(𝑡) = − 𝑡 + 𝑘 ⇒ 𝑀 𝑡 = 𝑐𝑒 200 ,
200
Para 𝑡 = 0
3𝑡
−
𝑀 0 = 𝑐 = 1500 ⇒ 𝑀 𝑡 = 1500𝑒 200 (2)
Con una concentración de 0.02 g/L, en el tanque se tendrá, 𝑀(𝑡) =20g,

reemplazando en la ecuación (2), se tiene:
200
20 = 1500𝑒 −3𝑡/200 ⇒ 𝑡 = 𝑙𝑛 75 = 287.83 minutos,
3
Por tanto, el tiempo para obtener una concentración de 0.02 g/L en el tanque
de aprox. 4 horas 48 min.

Ejemplo 3.
Un tanque contiene 300 litros de una solución, preparada con 200 kg de sal,
disueltos en el agua. Se bombea agua pura hacia el tanque a razón de 4 l/s y la
mezcla homogénea, producto de la agitación, se extrae del tanque a la misma
razón. Determine:
a) El número de kilogramos 𝑥 𝑡 de sal que hay en el tanque en un instante 𝑡.
b) ¿Qué cantidad de sal aproximadamente habrá en el tanque después de 5

minutos de ingresar agua pura?
c) ¿Cuánto tiempo pasará para que la cantidad de sal inicial en el tanque se

reduzca al 2 %?

Ejemplo 4.
Un tanque contiene inicialmente 800 𝑙𝑖𝑡𝑟𝑜𝑠 de salmuera con una concentración

de sal de 200 gramos por litro. Al tanque ingresa una salmuera con una
concentración de sal de 100 gramos por litro a razón de 8 𝑙𝑖𝑡𝑟𝑜𝑠/𝑚𝑖𝑛𝑢𝑡𝑜. La
mezcla agitada y homogenizada sale del tanque a razón de 5 𝑙𝑖𝑡𝑟𝑜𝑠/𝑚𝑖𝑛𝑢𝑡𝑜.
Asumiendo que la capacidad del tanque es mayor a 2000 litros, calcule:
a) Halle la cantidad de sal en el tanque, después de 𝑡 𝑚𝑖𝑛𝑢𝑡𝑜𝑠 de iniciado el
proceso.
b) Calcule la cantidad de sal dentro del tanque después de 2 horas de iniciado el
proceso.

Gracias

Aplicaciones de EDO Primer Orden 2023-2

Cargado por

Copyright:

Formatos disponibles

Aplicaciones de EDO Primer Orden 2023-2

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Aplicaciones de EDO Primer Orden 2023-2

Cargado por

Copyright:

Formatos disponibles

Universidad Nacional Agraria La Molina

CURSO: CC2076 CÁLCULO PARA CIENCIAS

4.3.5 Aplicaciones de ecuaciones diferenciales de primer orden

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Trayectorias ortogonales (TO)

Una trayectoria ortogonal de una

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Encuentre las trayectorias

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Integrando obtenemos las TO.

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Encuentre las trayectorias ortogonales a la familia de curvas:

Derivado implícitamente en la ecuación, tenemos:

Para hallar las TO a la familia de curvas, hacemos el cambio:

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Como exacta con F.I.

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Multiplicando por 𝑢 𝑥 se obtiene la ecuación diferencial exacta:

Derivando 𝐹 de la ecuación (*) respecto de “y” e igualamos a 𝑁:

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Derivado implícitamente, tenemos:

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Para hallar las TO a la familia de curvas, hacemos el cambio:

de donde se obtiene la ecuación diferencial

En este caso la ecuación (2) es una ecuación de Bernoulli ( 𝑘 = 3 ) (también se

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

multiplicando por 𝑢(𝑥) a la ecuación (3) se tiene

resolviendo y reemplazando 𝑧 se obtiene la solución general, que son las T.O. a la

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

La Ley de enfriamiento de Newton afirma que la rapidez de enfriamiento de un

T(t): temperatura del objeto en el tiempo t (en ºC, ºF, ...)

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

t : tiempo (s, min, ...)

Para el tiempo inicial 𝑡0 = 0 , se tiene un problema de valor inicial:

La solución general es:

y la solución particular es:

𝑇 = 𝑇𝑎 + (𝑇0 − 𝑇𝑎 )𝑒 −𝑘𝑡 ; 𝑘 > 0,

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Cierta habitación tiene una temperatura constante de 60ºC. Si un objeto se coloca

𝑇ℎ ∶ temperatura de la habitación (ºC)

𝑇 𝑡 ∶ temperatura del objeto en el instante “t” (t en minutos)

𝑇0 ∶ temperatura inicial del objeto (ºC).

Por la Ley de Enfriamiento de Newton, se tiene la ecuación diferencial:

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Resolviendo la ecuación diferencial (1) y reemplazando datos, tenemos:

𝑇 𝑡 = 𝑇𝑎 + 𝑇0 − 𝑇𝑎 𝑒 −𝑘𝑡 = 60 + 350 − 60 𝑒 −𝑘𝑡 = 60 + 290𝑒 −𝑘𝑡 … (2)

Para t = 25 minutos se tiene una temperatura de 𝑇 25 = 150 ºC, entonces, se tiene:

Entonces, en la ecuación (2), para una temperatura 𝑇 𝑡 = 70°C, se tiene:

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

A las 15 h 00 min pm un cuerpo que tiene una temperatura de 70℃ es depositado en

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

Si consideramos que la rapidez de crecimiento de una población de individuos es

P(t): población de individuos en el tiempo t.

t : tiempo transcurrido (s, min, días, años, ...), y t₀ el tiempo inicial.

P₀ : Población inicial de individuos.

CÁLCULO PARA CIENCIAS 2023 - II - UNALM

• Si 𝑘 > 0, el modelo representa una función de crecimiento.