Tema5 EcuacionesDiferencialesSeparables Aplicaciones

5.
Ecuaciones diferenciales 202
5.6 Aplicaciones de las ecuaciones diferenciales

Las ecuaciones diferenciales, debido a que relacionan los valores de una función con los de su(s) derivada(s),
son una herramienta fundamental en el tratamiento matemático de cualquier fenómeno dinámico, es decir,
que involucre magnitudes que cambian con el tiempo (o con cualquier otra magnitud). Por ello, sus campos
de aplicación son numerosos en física, química, biología, economía, . . . Se presentan a continuación algunos
ejemplos.
Ejemplo 5.23
En 1990 se arrojaron a un lago 1000 ejemplares de cierta especie de peces, de la que
previamente no había ninguno. En 1997 se estimó que la cantidad de peces de esa especie que
había en el lago en aquel momento era de 3000. Suponiendo que la velocidad de crecimiento de
la población de peces es constante, calcular la cantidad de peces en los años 2000 y 2010.
Que la velocidad de crecimiento de la población sea constante significa que, si llamamos
p(t) ≡ número de peces en el instante t
se tiene que
p (t) = k (constante) (5.13)
El valor de esta constante, k, no lo conocemos, de momento, pero veremos cómo se puede deducir utilizando
adecuadamente el resto de la información de que disponemos.
La ecuación (5.13) se puede resolver (dejando la constante k como un parámetro) y se tiene
p(t) = kt + C, C ∈ R arbitraria (5.14)
Ahora tenemos dos constantes “desconocidas”: k y C. Pero también tenemos dos informaciones que utilizar:
sabemos que
1. p(0) = 1000 (inicialmente había 1000 peces)
2. p(7) = 3000 (7 años después había 3000 peces)

Sustituyendo estos valores en (5.14) se tiene:

1000 = p(0) = k · 0 + C = C ⇔ C = 1000
2000
3000 = p(7) = k · 7 + C = 7k + 1000 ⇔ 7k = 2000 ⇔ k =
7
Con esto ya se tiene la expresión exacta de la función que nos da el número de peces que hay en el lago en
cualquier instante t:
2000
p(t) = t + 1000
7
y, con ella, ya se puede calcular lo que nos piden:
2000 27000
p(10) = · 10 + 1000 = ≈ 3857
7 7
2000 47000
p(20) = · 20 + 1000 = ≈ 6714
7 7
Así pues, la solución es
En el año 2000 había 3857 peces.
En el año 2010 había 6714 peces.
Matemáticas Aplicadas a la Biología - Grado en Biología R. Echevarría - Dpto. EDAN - Univ. de Sevilla
5. Ecuaciones diferenciales 203
Ejemplo 5.24
Si el número de bacterias contenidas en 1 litro de leche se duplica en 4 horas y suponiendo
que la tasa de multiplicación es constante, calcular en cuánto tiempo se hará 25 veces mayor.
Sea y(t) el número de bacterias en el instante t.

Suponer que la tasa de multiplicación de la población de bacterias es constante consiste en suponer que
y (t) = k k = constante (5.15)
El valor de la constante k, que de momento es desconocido, se puede deducir a partir de la información adicional
que tenemos.
Comenzamos por resolver la ecuación diferencial (5.15):
y(t) = kt + C, C ∈ R arbitraria (5.16)
La información de que disponemos es

y(0) = y0 número inicial de bacterias
y(4) = 2y0 el número de bacterias se duplica en 4 horas
Sustituimos estos datos en (5.16)

y0 = y(0) = k · 0 + C ⇔ C = y0
y0
2y0 = y(4) = k · 4 + C = 4k + y0 ⇔ y0 = 4k ⇔ k =
4
En consecuencia la función que nos da el número de bacterias en cualquier instante t es
y0 y0
y(t) = t + y0 = (t + 4)
4 4
siendo y0 = número inicial de bacterias.
Lo que se desea saber es en qué instante, t, el número de bacterias será igual a 25 veces el número que había
inicialmente.
y0
25 y0 = y(t) = (t + 4) ⇔ 100 = t + 4 ⇔ t = 100 − 4 = 96
4
Así pues, la solución es 96 horas .
5.6.1 Dinámica de poblaciones: modelo de Malthus o exponencial

El comportamiento de una población de seres vivos cuyo número de individuos varía en el tiempo puede ser
matemáticamente modelada mediante ecuaciones diferenciales y constituye, de hecho, uno de los principales
campos de aplicación de las Matemáticas a la Biología.
Cuando una población no está sujeta a condicionantes externos (falta de alimentos, competencia por el espacio,
por los recursos, . . . ) su ritmo de crecimiento o decrecimiento es debido únicamente al equilibrio entre su tasa
de natalidad y su tasa de mortandad: la velocidad de crecimiento de la población (o de decrecimiento, si nacen
menos individuos de los que mueren) es proporcional al número de individuos que la componen.
Para expresar esto matemáticamente, denotemos
N = N (t) número de habitantes en el instante t.
Entonces, la velocidad de crecimiento de la población, N’(t), verifica la siguiente ecuación diferencial:
N = r N, (5.17)
donde r es una constante, que caracteriza la tasa de crecimiento de la población, y que usualmente se determina
experimentalmente.
Si r > 0 la población aumentará de tamaño, por ser la velocidad de crecimiento positiva, mientras que si r < 0
la población disminuirá de tamaño.
Si en el instante inicial t = 0, el número de individuos es N (0) = N0 , entonces N (t) es solución del siguiente
problema de valor inicial:

N = r N t ≥ 0
(5.18)
N (0) = N0 .
Esta ecuación se resuelve fácilmente, ya que es de variables sepa- 20

rables (ver la Sección 5.3): 18
16
1
dN =
14
r dt
N 12
r=0.06
10
ln |N | = rt + C 8 r=0.04
6 r=0
N = C er t 4
2 r=−0.1
e, imponiendo la condición inicial N (0) = N0 , se obtiene 0
0 5 10 15 20
N = N0 e r t ,
Figura 5.2: Representación gráfica de la
cuya gráfica, para algunos valores de r, se representa en la Figu- función N = 5 ert , solución de (5.18)
ra 5.2. con N0 = 5, para varios valores de r.
Obsérvese que cuanto mas grande sea r, mas rápido es el crecimiento de la población, y que cuando r < 0 la
población decrece. Para r = 0 el tamaño de la población permanece constante.
Este modelo de crecimiento de poblaciones recibe su nombre de Thomas Malthus (1766-1843), un clérigo y
economista británico considerado el padre de la demografía. Basándose en este modelo, él dedujo que el creci-
miento (exponencial) del número de seres humanos sobre la Tierra conduciría a épocas de grandes hambrunas,
ya que la cantidad disponible de alimentos no aumentaría en la misma proporción que la población humana.
Este modelo de crecimiento de poblaciones es, como resulta obvio, excesivamente simple para reflejar situaciones
tan complejas como la de la población humana sobre la tierra. Sin embargo, resulta útil para modelizar mate-
máticamente algunos experimentos controlados en laboratorio con determinadas especies de microorganismos,
en sus etapas iniciales de desarrollo. Por ejemplo, si se inicia el cultivo de una pequeña colonia de bacterias
sobre un sustrato rico en nutrientes, entonces las bacterias pueden crecer y reproducirse sin restricciones, al
menos durante un cierto periodo de tiempo. (Un modelo más elaborado de dinámica de poblaciones, en el que
se imponen restricciones al crecimiento de la población, teniendo en cuenta otros aspectos vitales, se expone en
la Sección 5.6.6).
Ejemplo 5.25
(Cultivo de bacterias en laboratorio) Se sabe que la tasa de crecimiento de una determinada
población de bacterias es directamente proporcional al número de bacterias existentes. Se
realiza un cultivo en laboratorio, introduciendo 2.5 millones de bacterias en un recipiente. Se ob-
serva que la población se duplica cada 3 horas. Calcular la población existente al cabo de 11 horas.
Denotemos por P (t) al número de bacterias (en millones) que forman la población en el instante de tiempo t.
Se comienza a medir el tiempo (t = 0) en el instante en que se inicia el cultivo en el laboratorio.
Según se indica en el enunciado, la tasa de crecimiento de la población (velocidad a la que crece), P (t), es
directamente proporcional al número de bacterias de la población, es decir a P (t), lo que significa que es de la
forma kP (t) para alguna constante k que, de momento, no conocemos.
Esto significa que la población considerada sigue la ley (de Malthus):
P = kP ecuación diferencial cuyas soluciones son P (t) = C ekt

Para determinar las dos constantes C y k hay que utilizar las dos informaciones dadas:

P (0) = 2.5 (millones de bacterias)
P (3) = 2 × 2.5 = 5 (millones de bacterias)
De la primera de ellas se tiene

2.5 = P (0) = C ⇔ P (t) = 2.5 ekt
y de la segunda
5 ln(2)
5 = P (3) = 2.5 e3k ⇔ e3k = =2 ⇔ k= ≈ 0.231.
2.5 3
Luego, finalmente, la ley seguida por la población de bacterias es
P (t) = 2.5 e0.231 t .
El conocimiento de esta función nos permite conocer el número de bacterias que habrá en el cultivo en
cualquier instante (siempre y cuando, naturalmente, el modelo siga siendo válido). Por ejemplo, para saber
cuántas bacterias habrá 11 horas después de iniciar el experimento, bastará calcular
P (11) = 2.5 e0.231×11 ≈ 31.75.
Al cabo de 11 horas habrá aproximadamente 31.75 millones de bacterias
Ejemplo 5.26
(Población mundial). La población mundial en el año 1985 era de aproximadamente 4830
millones de personas y, en aquel momento, crecía a un ritmo de un 1.73 % por año. Suponiendo
que el crecimiento de la población se rigiera por el modelo exponencial, calcular el valor
estimado de la población mundial en el año 2010.
La ley de Malthus (o de crecimiento exponencial) dice que el número de individuos de la población en el

instante t, P (t), verifica la ecuación diferencial:
P (t) = kP (t), cuya solución general es P (t) = C ekt
En esta expresión hay dos constantes que no se conocen (de momento): k y C. Para determinar su valor
utilizaremos el resto de la información:
1. P (1985) = 4830 millones.
2. La población crece un 1.73 % cada año, de donde, por ejemplo, en el año 1986, la población se habría
incrementado en un 1.73 % de 4830 millones, es decir

1.73 1.73
P (1986) = 4830 + 4830 = 1 + 4830 = 4913 millones.
100 100
De ambos datos se tiene:

4830 = P (1985) = C e1985 k 4830 C e1985 k e1985 k
=⇒ = = = e1985 k · e−1986 k = e−k ,
4913 = P (1986) = C e1986 k 4913 C e1986 k e1986 k
y de aquí
4830 4830
ln = −k ⇔ k = − ln ≈ 0.0170
4913 4913
Ahora, una vez conocido el valor de k, se tiene:
4830
4830 = P (1985) = C e0.0170×1985 = C e33.7450 ⇔ C = ≈ 1.0683 × 10−11
e33.7450
Así, gracias a la información proporcionada se tienen ya los valores de las constantes C y k y por tanto la
expresión de P (t):
P (t) = 1.0683 × 10−11 e0.0170 t
Utilizando esta expresión se deduciría que el número de seres humanos en la tierra en el año 2010 sería:
P (2010) = 1.0683 × 10−11 e0.0170×2010 ≈ 7388 millones de personas
(la población real en el año 2010 era de 6972 millones de personas).

Observación: este ejercicio también se puede hacer (y, de hecho, los cálculos son más fáciles) situando el
origen, t = 0, de la variable independiente en el año 1985, de modo que el año 1986 correspondería a t = 1 y
el año 2010 correspondería a t = 25. Entonces tendríamos la información P (0) = 4830 y P (1) = 4913 y lo que
se desea es calcular P (25).
5.6.2 Ley de enfriamiento de Newton

En determinadas condiciones, la velocidad a la que cambia la temperatura de un objeto es proporcional a la
diferencia entre la temperatura del ambiente que lo rodea y su propia temperatura. Si se denota por T (t)
la temperatura del objeto en el instante t, la ley anterior se expresa matemáticamente mediante la siguiente
ecuación diferencial ordinaria:
T (t) = k(M − T (t)), (5.19)
donde M es la temperatura del medio (que se supone constante)

y k es la constante de proporcionalidad, propia del objeto.
50 T0=55
Si en el instante inicial, t = 0, la temperatura toma el valor T0 ,
entonces la temperatura del objeto en cualquier instante posterior 40
T (t), viene dada por la solución del problema de valor inicial: 30

T0=30

T = k(M − T ), 20
(5.20)
T (0) = T0 . T0=15
10
Esta ecuación es de variables separables y su solución general es
0
T (t) = M + Ce−kt ,
0 50 100 150 200 250
C ∈ R arbitraria.
Figura 5.3: Representación gráfica de
La solución particular que verifica T (0) = T0 es
la solución de (5.20), para M = 25,
T (t) = M + (T0 − M )e−kt . k = 0.02 y varios valores del dato inicial
T0 .
En la Figura 5.3 están representadas las soluciones del problema (5.20) para diversos valores del dato inicial T0 .
Obsérvese que, como es obvio intuitivamente, la temperatura del objeto varía más rápidamente cuanto mayor
es la diferencia entre la temperatura inicial del objeto y la temperatura del medio.
Por otro lado, sea cual sea su temperatura inicial, la temperatura del objeto tiende, cuando pasa el tiempo, a
igualarse con la temperatura del medio: todas las soluciones tienen una asíntota horizontal en T = M .
Ejemplo 5.27
(Ley de enfriamiento de Newton) Un recipiente con agua hirviendo (100◦ C) se retira del
fuego en el instante t = 0 y se deja enfriar en una habitación grande que se encuentra a una
temperatura constante de 20◦ C. Sabiendo que pasados 5 minutos la temperatura del agua se ha
enfriado hasta 80◦ C:
a) Determinar la constante de proporcionalidad k.
b) Determinar el tiempo que tardará el agua del recipiente en descender hasta una temperatura
de 30◦ C.
a) Sea y = y(t) la temperatura del agua (en grados Celsius) en el instante de tiempo t (medido en minutos).
Según la ley de enfriamiento de Newton, la temperatura del objeto sigue la ley
y = k(20 − y),
donde k es una constante propia del objeto.

Comenzamos observando que esta ecuación tiene la solución trivial y = 20 (constante).
La ecuación es de variables separables y se integra fácilmente:

1
dy = k dt ⇔ − ln |20 − y| = kt + C ⇔ y = 20 − C e−kt .
20 − y
La solución trivial y = 0 está contenida en esta familia para el valor de C = 0.
En la expresión de y hay 2 constantes que determinar: k y C. Para determinarlas disponemos de 2 datos:
y(0) = 100 e y(5) = 80
(1) De 100 = y(0) = 20 − C se tiene que C = −80

(2) De 80 = y(5) = 20 + 80 e−5k se tiene
3
80 − 20 60 3 1 3
= = = e−5k ⇔ −5k = ln ⇔ k = − ln ⇔ k ≈ 0.0575
80 80 4 4 5 4
En consecuencia, la función que da la temperatura del agua es:
y(t) = 20 + 80 e−0.0575 t
b) Se trata ahora de averiguar para qué valor de t alcanza y(t) (descendiendo) el valor 30◦ C. Es decir, para
qué valor de t se tiene
30 = 20 + 80 e−0.0575 t
Operando en esta ecuación se tiene
30 − 20 10 1 1 −1 1
= = = e−0.0575 t ⇔ −0.0575 t = ln ⇔ t= ln ≈ 36.1642
80 80 8 8 0.0575 8
Es decir, aproximadamente 36 minutos .
Ejemplo 5.28
Un cadáver es encontrado en una nave industrial que está a una temperatura constante de
20◦ C. En el momento de ser encontrado, la temperatura del cadáver es de 35◦ C. Al cabo de
una hora su temperatura ha descendido a 34◦ C. Suponiendo que en el momento de la muerte
la temperatura del cuerpo era de 37◦ C, y que se cumple la Ley de Enfriamiento de Newton,
calcular a qué hora se produjo la muerte.
Denotamos por T = T (t) la temperatura del cadáver en el instante t, comenzando a contar el tiempo en el
momento del crimen. Puesto que sigue la ley de Newton y en el momento inicial (t = 0) era de 37◦ C, la función
T (t) es la solución del siguiente problema de valor inicial:

T = k(M − T ) = k(20 − T )
(P)
T (0) = 37
La solución general de la anterior ecuación es (véase el Ejemplo 5.27) T (t) = 20 − C e−kt . La solución trivial
T = 20 está incluída para C = 0.
Lo que queremos saber es el tiempo pasado desde el momento de la muerte hasta que se encontró el cadaver.
Si situamos el momento de la muerte en el instante t = 0, y denotamos por t̃ al instante (desconocido de
momento) en que se encontró el cadaver, la información que tenemos es la siguiente:
⎧
⎨ T (0) = 37
T (t̃) = 35
⎩
T (t̃ + 1) = 34
Con estos 3 datos debemos ser capaces de encontrar los valores de k, de C y de t̃.
37 = T (0) = 20 − C ⇔ C = 20 − 37 =⇒ C = −17
35 − 20 15
35 = T (t̃) = 20 + 17 e−kt̃ ⇔ e−kt̃ = =
17 17
15 −k 34 − 20 14
34 = T (t̃ + 1) = 20 + 17 e−k(t̃+1) = 20 + 17 e−kt̃ e−k = 20 + 17 e = 20 + 15 e−k ⇔ e−k = =
17 15 15
De la última igualdad se tiene que

14 14
−k = ln =⇒ k = − ln ≈ 0.0690
15 15
15
Una vez conocido el valor de k, de la igualdad e−kt̃ = se puede despejar t̃ tomando logaritmos en
17
ambos miembros:

15 15 1 15
e−kt̃ = ⇔ −k t̃ = ln ⇔ t̃ = − ln ⇔ t̃ ≈ 1.8141 horas ≈ 1 hora 49 minutos
17 17 k 17
Así pues, el cadáver fué encontrado 1 hora y 49 minutos después de su muerte.
5.6.3 Dinámica de crecimiento de un individuo: modelo de Bertalanffy.

En los años 50 del siglo XX, el biólogo austriaco L. von Bertalanffy (1901-1972) desarrolló un modelo matemático
para la talla de un individuo en función de su edad, que se utiliza con frecuencia para predecir el tamaño de los
peces.
Sea L(t) la longitud del individuo en la edad t y sea A la talla máxima de la especie, es decir la talla máxima
alcanzable por un pez adulto.
La ley de crecimiento de este modelo dice que la velocidad de crecimiento es proporcional a la diferencia entre
la longitud actual y la longitud máxima:
Figura 5.4: Modelo de Bertalanffy.
L (t) = k(A − L(t)),

siendo k > 0, la constante de proporcionalidad, propia de cada especie.
Si en el instante inicial, t = 0, la longitud del individuo es 0 < L0 < A , entonces la función L(t), talla en el
instante t, será solución del siguiente problema de valor inicial:

L = k(A − L)
(5.21)
L(0) = L0 .
Como la diferencia entre la longitud actual y la longitud máxima alcanzable disminuye con el tiempo, la velocidad
de crecimiento disminuye también con el tiempo, lo que implica que los ejemplares de menor edad crecen a mayor
velocidad que los de mayor edad. En este modelo, la velocidad de crecimiento es siempre positiva. Esto significa
que los peces crecen durante toda su vida, que es lo que ocurre en la realidad.
La ecuación diferencial de (5.21) se puede integrar fácilmente, ya que es de variables separables:

dL
= k dt ⇐⇒ − ln |A − L| = kt + C ⇐⇒ A − L = Ce−kt .
A−L
Por tanto, la solución general de la ecuación es
L = A + Ce−kt , C ∈ R, arbitraria.
Imponiendo la condición inicial, L(0) = L0 , se tiene
L0 = L(0) = A + Ce0 = A + C ⇐⇒ C = L0 − A,
luego la solución del problema (5.21) es
L(t) = A + (L0 − A)e−kt .
En la Figura 5.5 está representada la solución del

problema (5.21) para A = 50, k = 0.5 y L0 = 0. 50
Obsérvese que la recta horizontal L = A es una

40
asíntota horizontal de la solución, es decir,
30
lı́m L(t) = A,
t→+∞
20
lo que expresa matemáticamente el hecho de que
la talla de los peces tiende, cuando pasa el tiem- 10
po, a aproximarse al valor A, pero sin nunca al- L0=0
canzarlo. 0
0 2 4 6 8 10 12
Por ello se puede decir que A es la longitud
asintótica de la especie. Figura 5.5: Representación gráfica de la solución
de (5.21), para A = 50, k = 0.5 y L0 = 0 .
Ejemplo 5.29
(Modelo de Bertalanffy ) Sea L(t) la longitud (en centímetros) de un pez en el tiempo t, medido
en meses. Se supone que el pez crece de acuerdo con la siguiente ley (de von Bertalanffy):

L = k(34 − L)
L(0) = 2.
1) Sabiendo que a la edad de 4 meses, el pez mide 10 centímetros, determinar la constante de

crecimiento k.
2) Calcular la longitud del pez a los 10 meses.
3) Calcular lı́m L(t) y dar una interpretación del resultado en el marco de la dinámica del
t→∞
crecimiento del pez.
La solución del problema de valor inicial se calcula fácilmente por ser la ecuación de variables separables:

1
L = k(34 − L) ⇔ dL = k dt ⇔ − ln |34 − L| = kt + C
34 − L
de donde se tiene L = 34 − Ce−kt e, imponiendo la condición inicial L(0) = 2, se encuentra el valor de la

constante C = 32.
Luego la longitud del pez viene dada por
L(t) = 34 − 32 e−kt .
Para determinar el valor de k es necesario utilizar más información: L(4) = 10. Entonces,

−4k −4k 24 3 1 3
10 = L(4) = 34 − 32 e ⇔ e = = ⇔ k = − ln = 0.0719.
32 4 4 4
Una vez conocido el valor de k se puede calcular la longitud del

pez en cualquier instante t > 0:
L(10) = 34 − 32 e−10k ≈ 18.4 cm.
Por último, es obvio que

34
−kt −kt
lı́m L(t) = lı́m 34 − 32 e = 34 − 32 lı́m e = 34,
t→+∞ t→+∞ t→+∞
18.4
lo cual significa que la curva que representa la longitud del pez
tiene una asíntota horizontal en L = 34. El pez sigue creciendo,
pero cada vez a menor velocidad, y su longitud tiende a acercarse 2
al valor 34, aunque sin nunca llegar a alcanzarlo. 0 10 100
5.6.4 Problemas de mezclas

En esta sección se estudian ciertas ecuaciones diferenciales que aparecen en problemas en los que se mezclan
dos fluidos.
Más concretamente, se considera un recipiente que contiene una cantidad de V litros de cierto fluido, en el que
se encuentra disuelta una cantidad, x0 , de cierta sustancia. En el recipiente entra constantemente fluido con
una concentración de ce gramos por litro y a una velocidad de ve litros por minuto. Se supone que los fluidos
en el recipiente se mezclan de forma instantánea y que la mezcla sale del recipiente a una velocidad de vs litros
por minuto.
Lo que se desea es determinar una función que indique la cantidad de sustancia que hay en el interior del
recipiente en cada instante, t.
Llamemos v(t) a la cantidad de fluido (litros) presente en el recipiente en el instante t, y x(t) a la cantidad de
sustancia disuelta (gramos) en el instante t, de forma que la concentración de sustancia disuelta en el instante
t es x(t)/v(t) gramos por litro.
La variación de la magnitud x(t) por unidad de tiempo es x (t) y viene dada por la diferencia entre la cantidad
de sustancia que entra (por unidad de tiempo) y la cantidad de sustancia que sale (por unidad de tiempo):
Variación de x(t) Cantidad de sustancia que Cantidad de sustancia que

x (t) = = −
por unidad de tiempo entra por unidad de tiempo sale por unidad de tiempo
Puesto que entran ve litros por minuto, que contienen una concentración ce gramos de sustancia por litro, se
tiene que entran ce · ve gramos por minuto de sustancia.
La concentración de sustancia en el fluido que sale es la del fluido en el interior del recipiente, es decir x(t)/v(t)
gramos por litro. Puesto que salen vs litros por minuto, se tiene que salen x(t)vs /v(t) gramos por minuto de la
sustancia disuelta.
Así pues, la variación de la concentración, x (t), verifica:
x(t)
x (t) = ce ve − vs
v(t)
La expresión de v(t), cantidad de fluido en el recipiente en el instante t, deberá ser determinada en cada caso,
ya que depende de la cantidad inicial y de las velocidades de entrada y salida del fluido en el recipiente. Si, por
ejemplo, la velocidad de entrada de fluido es igual a la velocidad de salida, entonces el volumen en el interior
del recipiente permanecerá constante.
Ejemplo 5.30
(Problema de mezclas) Un depósito contiene 100 litros de una disolución salina cuya concen-
tración es de 2.5 gramos de sal por litro. Una disolución conteniendo 2 gramos de sal por litro
entra en el depósito a razón de 5 litros por minuto y la mezcla (que se supone uniforme de
forma instantánea) sale del depósito a la misma velocidad. Encontrar la cantidad de sal que
hay en cada instante en el depósito.
Puesto que la velocidad a la que entra el líquido en el depósito es la misma a la que sale, en el depósito siempre
hay la misma cantidad de líquido: 100 litros.
Sea x(t) la cantidad de sal en el depósito en el instante t.
La variación por unidad de tiempo de la cantidad de sal en el depósito es:
x (t) = cantidad que entra por unidad de tiempo − cantidad que sale por unidad de tiempo
En el depósito entran 5l. por minuto de una disolución con 2gr. por litro, luego entran 10gr. de sal por minuto.
Puesto que la cantidad de sal en el depósito es x(t) y la cantidad de líquido que hay es 100l., la concentración
de la disolución en el depósito es de x(t)/100 gramos por litro. Esta disolución sale a una velocidad de 5 litros
por minuto, por lo tanto la sal sale a una velocidad de 5x(t)/100 gramos por minuto.
Así pues, se tiene:
5x
x = 10 −
100
Esta ecuación es de variables separables:
5x 1000 − 5x 1 1
x = 10 − = ⇔ x =
100 100 1000 − 5x 100

1 1 1 1
⇔ dx = dt ⇔ − ln |1000 − 5x| = t+C
1000 − 5x 100 5 100
5 1
⇔ ln |1000 − 5x| = − t + C = − t + C = −0.05t + C ⇔ 1000 − 5x = C e−0.05t
100 20
1000 − C e−0.05t
⇔ 5x = 1000 − C e−0.05t ⇔ x = = 200 − Ce−0.05t
5
Así pues, la solución general de la ecuación diferencial es
x = 200 − Ce−0.05t
Puesto que, inicialmente, la concentración de sal en el depósito era de 2.5 gramos por litro, la cantidad de sal
inicial era de
x(0) = 2.5 × 100 = 250
Sustituyendo esta condición inicial en la expresión de la solución general se tiene
250 = x(0) = 200 − C ⇔ C = −50
Luego la función que nos da la cantidad de sal en cualquier instante t es:
x(t) = 200 + 50e−0.05t
Ejemplo 5.31
(Problema de mezclas) La corriente sanguínea lleva un medicamento hacia el interior de un
órgano a razón de 3 cm3 /sg y sale de él a la misma velocidad. El órgano tiene un volumen de
125 cm3 . Si la concentración del medicamento en la sangre que entra en el órgano es de 0.2
gr/cm3 , se pide:
1) ¿Cuál es la concentración del medicamento en el órgano en cada instante si inicialmente no
había vestigio alguno del medicamento?
2) ¿Cuándo la concentración del medicamento en el órgano será de 0.1 gr/cm3 ?
La cantidad de medicamento que entra en el órgano por segundo es:
0.2 × 3 = 0.6 gramos
Si denotamos por x(t) la cantidad de medicamento presente en el órgano en el instante t se tendrá, puesto que
la sangre abandona el órgano a la misma velocidad a la que entra (3 cm3 /sg), que la cantidad de medicamento
que abandona el órgano por segundo será de
x(t) 3
3 = x(t)
125 125
En consecuencia, puesto que la variación por unidad de tiempo (i.e., por segundo) de la cantidad de medica-
mento viene dada por:
x (t) = cantidad que entra por segundo − cantidad que sale por segundo
se tiene
3 75 − 3x
x = 0.6 − x=
125 125
Esta ecuación es de variables separables:

1 1 1 t 3t
dx = dt ⇔ − ln |75−3x| = +C ⇔ ln |75−3x| = − +C ⇔ 75−3x = C e−3t/125
75 − 3x 125 3 125 125
Despejando aquí x se obtiene la solución general de la ecuación:
x = 25 − C e−3t/125
Puesto que, inicialmente, no había ninguna cantidad de medicamento en el órgano, la condición inicial para
x(t) es x(0) = 0, lo que conduce, sustituyendo, a:
0 = x(0) = 25 − C ⇔ C = 25
En consecuencia la función que nos da la cantidad de medicamento en el órgano en cada instante es
x(t) = 25(1 − e−3t/125 )
La concentración es la cantidad de medicamento dividido por el volumen del órgano, es decir

25 1
x(t)/125 = (1 − e−3t/125 ) = (1 − e−3t/125 )
125 5
Por lo tanto, la contestación a la primera pregunta es que
1
La concentración en el instante t es (1 − e−3t/125 )
5
Para contestar a la segunda pregunta hay que calcular para qué valor de t se verifica
1 3t
0.1 = (1 − e−3t/125 ) ⇔ 0.5 − 1 = −0.5 = −e−3t/125 ⇔ e−3t/125 = 0.5 ⇔ − = ln 0.5
5 125
125
⇔ t=− ln 0.5 ≈ 28.88 segundos
3
5.6.5 Dinámica de epidemias
Ejemplo 5.32
(Dinámica de epidemias) Un modelo simple de propagación de epidemias se obtiene cuando
se supone que la rapidez de contagio entre la población es directamente proporcional al número
de individuos contagiados multiplicado por el número de individuos no contagiados. Hallar la
solución general de esta ecuación.
Denotamos por I(t) el número de infectados por la epidemia en el instante t y por P (constante) el número
total de habitantes de la población, de forma que P − I(t) es el número de individuos no infectados. El modelo
establece que la velocidad de contagio I (t) es proporcional al número de infectados I(t) multiplicado por el
de no infectados P − I(t). En consecuencia se tiene
I = k I (P − I) (5.22)
donde k es la constante de proporcionalidad.

Esta ecuación es de variables separables y tiene las soluciones triviales I = 0 e I = P . Para calcular las demás:

1 dI 1
I = k I (P − I) ⇔ =k ⇔ dI = k dt = kt + C
I(P − I) dt I(P − I)
Para calcular la integral del primer miembro, que es racional, hay que escribir el integrando como una suma
de fracciones simples:
1 A B A = 1/P
= + ⇔
I(P − I) I P −I B = 1/P
En consecuencia, se tiene:

1 1/P 1/P 1 1 I
dI = + dI = (ln I − ln(P − I)) = ln = kt + C
I(P − I) I P −I P P P −I
I
⇔ ln = P (kt + C) = kP t + P C = kP t + C
P −I
I
⇔ = ekP t+C = ekP t eC = C ekP t
P −I
Operamos a continuación para despejar I en esta igualdad:
I = C ekP t (P − I) = CP ekP t − CekP t I
CP ekP t
⇔ I + CekP t I = I (1 + CekP t ) = CP ekP t ⇔ I =
1 + CekP t
Con esto ya tenemos la expresión de la solución general de la ecuación (5.22), que es mejor escribir dividiendo
numerador y denominador por CekP t :
P
I=
1 + Ce−kP t
La solución trivial I = P está contenida en esta familia para C = 0. Sin embargo la solución I = 0 no lo está:
para ningún valor que demos a la constante arbitraria C obtendremos la función I = 0.
En consecuencia, el conjunto de todas las soluciones de la ecuación es:
P
I= , ∀C ∈ R y además I = 0.
1 + Ce−kP t
Ejemplo 5.33
(Gripe aviar ) En una granja de 40.000 aves hay un pollo contagiado con la gripe aviar. Si
suponemos que la rapidez de contagio es directamente proporcional al número de aves conta-
giadas multiplicado por el número de no contagiadas, siendo la constante de proporcionalidad
k = 4 × 10−5 (midiendo el tiempo en días), determinar en cuánto tiempo un 75 % de los pollos de
la granja quedarían infectados.
Denotando por I(t) el número de pollos contagiados y por P el número total de pollos de la granja (población
total) se tiene la siguiente ecuación diferencial
I = k I (P − I)
donde k es la constante de proporcionalidad.

En este caso, P = 40000 y k = 4 × 10−5 = 0.00004 (de donde kP = 16 × 104 × 10−5 = 1.6).
Nos dicen, además, que inicialmente hay un pollo infectado, es decir, que se tiene I(0) = 1. En consecuencia, el
problema que hay que resolver para obtener la expresión de la función que representa el número de individuos
infectados en cualquier instante t es:
I = k I (P − I)
(5.23)
I(0) = 1
La solución general de esta ecuación diferencial es (véase Ejemplo 5.32):
P
I= (y además I = 0).
1 + Ce−kP t
Buscamos ahora la solución que verifica la condición inicial, I(0) = 1.
P 40000
1 = I(0) = ⇔ C = P −1 = 39999 =⇒ la solución del problema (5.23) es I(t) =
1+C 1 + 39999 e−1.6t
Buscamos ahora el valor del tiempo t para el cual I(t ) = 0.75 P = 30000. Para este t se tendrá
40000 40000 4
30000 = I(t ) = −1.6t
⇔ 1 + 39999 e−1.6t = =
1 + 39999 e 30000 3

1 4 1 1 1 1
⇔ e−1.6t = −1 = ⇔ −1.6t = ln ⇔ t = − ln
39999 3 119997 119997 1.6 119997
de donde se deduce que
El tiempo que tarda en estar contagiados el 75 % de los pollos es t ≈ 7.3
Ejemplo 5.34
Se sabe que la velocidad de propagación de una epidemia es proporcional al número de per-
sonas infectadas multiplicado por el número de personas no infectadas. Si denotamos por I(t)
el número de personas infectadas en el tiempo t, medido en días, y por P la población total, la
dinámica de la infección viene dada por
I = k I(P − I),
donde k > 0 es el coeficiente de proporcionalidad. En una población de 10000 habitantes se

detecta una enfermedad que afecta inicialmente a 50 personas. Al cabo de tres días, se observa
que son 250 las personas afectadas. Averiguar el número de enfermos que habrá pasados 12 días.
La ecuación I = k I(P − I) es de variables separables y su solución es (véase el ejercicio 5.32):
P
I(t) = (y además I = 0)
C e−kP t + 1
donde P = 10000. Para determinar las constantes C y k disponemos de la siguiente información:
I(0) = 50 e I(3) = 250.
En primer lugar,
P P
50 = I(0) = ⇔ C= − 1 = 199.
C +1 50
En segundo lugar,

P P 1 P
250 = I(3) = ⇔ 199 e−3kP + 1 = ⇔ e−3kP = −1
199 e−3kP + 1 250 199 250
de donde, tomando logaritmos en ambos miembros de la igualdad se tiene

1 P 1 1 P 0.5432
−3kP = ln −1 ⇔ k=− ln −1 = .
199 250 3P 199 250 P
En consecuencia, el número de infectados en cualquier instante t > 0 viene dado por
P 104
I(t) = −0.5432t = −0.5432t
199 · e +1 199 · e +1
y se tiene
104
I(12) = −0.5432×12 ≈ 7730
199 · e +1
Pasados 12 días habrá 7730 enfermos.
Ejemplo 5.35
(Dinámica de epidemias) En un campus universitario que tiene 1000 estudiantes hay un
único estudiante portador del virus de la gripe. Sea y(t) el número de estudiantes contagiados
en el día t. Si la velocidad con la que el virus se propaga es proporcional al producto entre el
número de estudiantes contagiados y el número de estudiantes no contagiados, se pide:
1) Determinar el número de personas enfermas en el día t si se sabe que pasados 4 días hay 50
enfermos.
2) Calcular cuándo habrá 500 estudiantes enfermos.
3) Si los estudiantes enfermos no se tratan con medicamentos, ¿qué número de enfermos habrá
cuando pase mucho tiempo? ¿Llegará a desaparecer la enfermedad?
Por lo que se indica, la función y(t) = número de estudiantes contagiados en el día t es solución de la ecuación
diferencial
y = ky(P − y)
donde P = 1000 es el número de individuos de la población. La solución general de esta ecuación es (véase
ejercicio 5.32):
P
y= (y además y = 0)
C e−kP t + 1
en cuya expresión hay dos constantes desconocidas (de momento) : k y C. Para determinarlas debemos usar
el resto de la información:
De y(0) = 1 se tiene
1000
1 = y(0) = ⇔ C + 1 = P = 1000 ⇔ C = 999
C +1
Por otra parte, de y(4) = 50 se tiene:
P
50 = y(4) = ⇔ 50 Ce−4kP + 1 = 50 C e−4kP + 50 = P
C e−4kP +1
Despejando de aquí e−4kP y tomando luego logaritmos en ambos miembros:

P − 50 P − 50
e−4kP = ⇔ ln e−4kP = −4kP = ln
50C 50C

1 P − 50 1 950
⇔ −kP = ln = ln ≈ −0.9906
4 50C 4 49950
Así pues,
1000
El número de personas enfermas el día t es y(t) =
999 e−0.9906 t + 1
Para saber cuándo habrá 500 estudiantes enfermos tenemos que calcular para qué valor de t se tiene

1000 −0.9906 t −0.9906 t 1 1
= 500 ⇔ 2 = 999 e + 1 ⇔ e = ⇔ −0.9906 t = ln
999 e−0.9906 t + 1 999 999

1 1
⇔ t=− ln ≈ 6.9723
0.9906 999
Por último, puesto que
P
lı́m y(t) = lı́m =P
t→+∞ t→+∞ C e−kP t
+1
resulta obvio que esta ley conduce a que, a la larga, la población entera resulta infectada.
5.6.6 Dinámica de poblaciones: ecuación logística

En la Sección 5.6.1, se ha considerado un modelo simple de la dinámica de poblaciones, en el que se supone
que no hay limitaciones de alimentos y, por tanto la población puede crecer de manera exponencial. El modelo
que se presenta ahora es un poco más complicado. En él se tiene en cuenta la existencia de circunstancias que
limitan el crecimiento exponencial de la población.
En determinadas condiciones, el crecimiento de algunas poblaciones se rige por la siguiente ley, denominada
logística:
p (t) = r p(t) − m p2 (t). (5.24)
En esta ecuación p(t) representa el número de individuos de la población existentes en el instante t. El primer
término de la derecha de esta ecuación (r p(t)) expresa matemáticamente el crecimiento natural de la población,
debido a la reproducción: la población crece de forma proporcional al número de individuos de la misma. El
segundo término (−m p2 (t)) intenta expresar el hecho de que, si los recursos (alimentos) son limitados, entonces
los individuos de la población “compiten” por ellos, impidiendo un crecimiento ilimitado. Este término hace
disminuir la velocidad a la que crece la población, razón por la que lleva signo menos.
Si en el instante inicial t = 0, el número de individuos es p(0) = p0 , entonces p = p(t) es solución del siguiente
problema de valor inicial:

p = r p − m p2 ,
(5.25)
p(0) = p0 .
La ecuación (5.24) es de variables separables, luego:

dp 1
= p(r − mp) ⇔ dp = dt.
dt p(r − mp)
Para calcular la integral de la izquierda hay que escribir el integrando como suma de fracciones simples:

1 A B A = 1/r
= + ⇔ 1 = A(r − mp) + Bp ⇐⇒
p(r − mp) p r − mp B = m/r
de donde, A = 1/r y B = m/r. Por lo tanto:

1 1/r m/r 1 1 m
dp = + dp = + dp = dt.
p(r − mp) p r − mp r p r − mp
Integrando, se obtiene
1
(ln |p| − ln |r − mp|) = t + C, con C ∈ R arbitraria
r
o, lo que es lo mismo,

p

ln = rt + C, con C ∈ R arbitraria.
r − mp
Tomando ahora exponenciales en ambos miembros de esta igualdad se tiene:
p Cr ert
= C ert ⇐⇒ p = Cr ert − Cm ert p ⇐⇒ p= .
r − mp 1 + Cm ert
Y de aquí, dividiendo numerador y denominador por Cert y renombrando la constante arbitraria C, se tiene,
finalmente, la expresión siguiente para la solución general de la ecuación logística:
r
p= .
m + C e−rt
Por tanto, la solución general de (5.24) es:

r
p(t) = , C ∈ R arbitraria. (5.26)
m + C e−rt
Esta ecuación tiene, además, las soluciones constantes p = β, para los valores de β que anulen el segundo
miembro de la ecuación diferencial, en este caso:

β=0
β(r − mβ) = 0 ⇐⇒
β = r/m,
La solución constante p = r/m está incluida en la expresión de la solución general, para el valor de C = 0. En
cambio, la solución constante p = 0 no se obtiene de la expresión de la solución general para ningún valor de la
constante arbitraria C: la ecuación logística tiene todas las soluciones dadas por (5.26) y, además, la solución
constante p = 0.
220
p0=20
200
La solución particular que verifica la condición p0=200
180
inicial p(0) = p0 se obtiene para el valor de la p =120
0
r − mp0 160
constante arbitraria C = y es: 140
p0 120
r p0
p(t) = . 100
mp0 + (r − mp0 ) e−rt 80
Su comportamiento cualitativo puede observarse 60

40
en la Figura 5.6 para varios valores de la condi-
20
ción inicial p0 .
0
0 50 100 150 200 250
Obsérvese que, sea cual sea el número de indi-
viduos de la población inicial, esta tiende, con Figura 5.6: Gráfica de la solución del problema
el tiempo, a estabilizarse en el valor constante (5.25) con r = 0.05 y m = 0.0003125, para varios
r
P = (asíntota horizontal de p(t)). valores de p0 .
m

Tema5 EcuacionesDiferencialesSeparables Aplicaciones

Cargado por

Copyright:

Formatos disponibles

Tema5 EcuacionesDiferencialesSeparables Aplicaciones

Cargado por

Información del documento

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

Tema5 EcuacionesDiferencialesSeparables Aplicaciones

Cargado por

Copyright:

Formatos disponibles

5.

Ecuaciones diferenciales 202

5.6 Aplicaciones de las ecuaciones diferenciales

Que la velocidad de crecimiento de la población sea constante signiﬁca que, si llamamos

p(t) ≡ número de peces en el instante t

p(t) = kt + C, C ∈ R arbitraria (5.14)

2. p(7) = 3000 (7 años después había 3000 peces)

En el año 2000 había 3857 peces.

En el año 2010 había 6714 peces.

Sea y(t) el número de bacterias en el instante t.

y (t) = k k = constante (5.15)

y(t) = kt + C, C ∈ R arbitraria (5.16)

La información de que disponemos es

Sustituimos estos datos en (5.16)

5.6.1 Dinámica de poblaciones: modelo de Malthus o exponencial

Para expresar esto matemáticamente, denotemos

N = N (t) número de habitantes en el instante t.

Entonces, la velocidad de crecimiento de la población, N’(t), veriﬁca la siguiente ecuación diferencial:

Esta ecuación se resuelve fácilmente, ya que es de variables sepa- 20

P = kP ecuación diferencial cuyas soluciones son P (t) = C ekt

De la primera de ellas se tiene

P (t) = 2.5 e0.231 t .

P (11) = 2.5 e0.231×11 ≈ 31.75.

Al cabo de 11 horas habrá aproximadamente 31.75 millones de bacterias

La ley de Malthus (o de crecimiento exponencial) dice que el número de individuos de la población en el

P (t) = kP (t), cuya solución general es P (t) = C ekt

De ambos datos se tiene:

P (2010) = 1.0683 × 10−11 e0.0170×2010 ≈ 7388 millones de personas

(la población real en el año 2010 era de 6972 millones de personas).

5.6.2 Ley de enfriamiento de Newton

donde M es la temperatura del medio (que se supone constante)

T (t), viene dada por la solución del problema de valor inicial: 30

donde k es una constante propia del objeto.

y(0) = 100 e y(5) = 80

(1) De 100 = y(0) = 20 − C se tiene que C = −80

En consecuencia, la función que da la temperatura del agua es:

Es decir, aproximadamente 36 minutos .

Así pues, el cadáver fué encontrado 1 hora y 49 minutos después de su muerte.

5.6.3 Dinámica de crecimiento de un individuo: modelo de Bertalanﬀy.

Figura 5.4: Modelo de Bertalanﬀy.

L (t) = k(A − L(t)),

La ecuación diferencial de (5.21) se puede integrar fácilmente, ya que es de variables separables:

Imponiendo la condición inicial, L(0) = L0 , se tiene

luego la solución del problema (5.21) es

L(t) = A + (L0 − A)e−kt .

En la Figura 5.5 está representada la solución del

Obsérvese que la recta horizontal L = A es una

1) Sabiendo que a la edad de 4 meses, el pez mide 10 centímetros, determinar la constante de

de donde se tiene L = 34 − Ce−kt e, imponiendo la condición inicial L(0) = 2, se encuentra el valor de la

Una vez conocido el valor de k se puede calcular la longitud del

L(10) = 34 − 32 e−10k ≈ 18.4 cm.

Por último, es obvio que

al valor 34, aunque sin nunca llegar a alcanzarlo. 0 10 100

5.6.4 Problemas de mezclas

Variación de x(t) Cantidad de sustancia que Cantidad de sustancia que

250 = x(0) = 200 − C ⇔ C = −50

Luego la función que nos da la cantidad de sal en cualquier instante t es:

x(t) = 200 + 50e−0.05t

La cantidad de medicamento que entra en el órgano por segundo es: