Capı́tulo 7. Sección 7.

1 [ clase 18 ] Cálculo Integral

Ecuaciones diferenciales y modelados Escuela de Matemáticas

13. Ecuaciones diferenciales y modelado mediante ecuaciones diferenciales [ Secciones 9.1-9.3 ]

Una ecuación diferencial es una ecuación que contiene una función desconocida y una o más de sus derivadas. Casi en todos los
casos una ecuación diferencial es la expresión matemática de una ley fı́sica. El orden de una ecuación diferencial es el orden de la
derivada de orden superior que se presente en la ecuación.
Ejemplo
Si y = y(x) entonces
y 00 + 3y 0 − 5y = 3x
es una ecuación diferencial de segundo orden.

Una función y = f (x) se llama solución de una ecuación diferencial si y y sus derivadas satisfacen la ecuación.
Ejemplo
y = xex es solución de la ecuación diferencial
y 0 − y = ex

Cuando se pide resolver una ecuación diferencial, se espera que encontremos todas las posibles soluciones de la ecuación.
Ejemplo
Resolver la ecuación diferencial
y 0 = x2
La solución general (todas las soluciones posibles) está dada por

x3
y= + C.
3

Cuando se trabaja con ecuaciones diferenciales, por lo general no estamos interesados en encontrar una familia de soluciones
(solución general ), sino que buscamos una solución que satisfaga algunas condiciones adicionales llamadas condiciones iniciales.
En otras palabras, en la mayorı́a de los casos se quiere resolver un problema de valor inicial (PVI) de la forma
 0
 y + p(x)y = q(x)

y(x0 ) = y0 condición inicial



Geométricamente, esto nos dice que toda la familia de curvas que es solución de la ecuación diferencial, escogemos aquella que
pasa por el punto (x0 , y0 ).
Ejemplo
Verifique que y = −t cos t − t es una solución del PVI
dy

 t
 = y + t2 sin t
dt


y(π) = 0

Solución: Observemos que y(π) = −π cos π − π = 0. Por otro lado,

dy
= −t(− sin t) + cos t(−1) − 1 = t sin t − cos t − 1
dt
Reemplazando en la ecuación diferencial nos queda
dy
t = t(t sin t − cos t − 1) = t2 sin t − t cos t − t = t2 sin t + y
dt

En general resolver una ecuación diferencial es una tarea difı́cil. Acá solo veremos un cierto tipo de ecuaciones las cuales
podremos resolver de manera explı́cita. Una ecuación separable es una ecuación diferencial de primer orden que se puede escribir
en la forma

1
Capı́tulo 7. Sección 7.1 [ clase 18 ] Cálculo Integral
Ecuaciones diferenciales y modelados Escuela de Matemáticas

dy
= g(x)f (y).
dx
Entonces si f (y) 6= 0 reescribimos lo anterior en forma diferencial como

h(y)dy = g(x)dx,
1
donde h(y) = . Luego, para resolver la ecuación basta con integrar a ambos lados de la igualdad en la forma
f (y)
Z Z
h(y)dy = g(x)dx.

Ejemplo
Resuelva las ecuaciones diferenciales
dy
(a) = xe−y
dx
(b) (y + sin y)y 0 = x + x3
Solución.

(a) En este caso podemos escribir

x2 x2
Z Z  
ey dy = xdx =⇒ ey dy = xdx =⇒ xdx =⇒ ey = + C =⇒ y = ln +C .
2 2

(b) Observemos que igual que en el caso anterior podemos escribir

(y + sin y)dy = (x + x3 )dx,

integrando en ambos lados

y2 x2 x4
Z Z
(y + sin y)dy = (x + x3 )dx =⇒ − cos y = + + C,
2 2 4
observar que aquı́ no podemos despejar la y entonces la solución se deja de manera implı́cita.

Aplicaciones
Problema 1.
Encuentre una ecuación de la curva que pasa por el punto (0, 1) y cuya pendiente en (x, y) es xy.
Solución.
dy
Tenemos que la pendiente en cada punto es = xy, luego despejamos y resolviendo la ecuación
dx
x2
Z Z
dy dy 2
= xdx =⇒ = xdx =⇒ ln |y| = + C =⇒ y = ±eC ex /2 ,
y y 2
2
ası́ y = kex /2
. Puesto que la curva para por el punto (0, 1), entonces 1 = ke0 , luego la curva está dada por
2
y = ex /2
.

Problema 2. Crecimiento y desintegración exponenciales

Uno de los modelos más simples es aquel en el que se supone que una población crece a una tasa proporcional al tamaño de la
población:
dP
= kP (t),
dt
donde P (t) es la población presente en un tiempo t. Este modelo se aplica cuanco la población se encuentra en unas condiciones
ideales.

2
Capı́tulo 7. Sección 7.1 [ clase 18 ] Cálculo Integral
Ecuaciones diferenciales y modelados Escuela de Matemáticas

En general, si y(t) es el valor de una cantidad y en un tiempo t y si la razón de cambio de y con respecto a t es proporcional a su
tamaño y(t) en cualquier tiempo, entonces
dy
= ky,
dt
donde k es una constante. Esta ecuación se llama ley de crecimiento natural si k > 0 o ley de desintegración natural si k < 0. En
cualquier caso, podemos resolverla como una ecuación diferencial separable y usarla para resolver algunos problemas:
Z
dy dy
= ky =⇒ = kdt =⇒ ln |y| = kt + C =⇒ y = ±eC ekt ,
dt y

de donde y = Aekt , además y(0) = A, luego podemos decir que la solución del PVI
dy


 = ky
dt


y(0) = y0
es
y(t) = y0 ekt .

Ejemplo 1.
Un cultivo de bacterias contiene inicialmente 100 células y crece con una rapidez proporcional a su tamaño. Después de
una hora, la población ha aumentado a 420.
(a) Encuentre una expresión para el número de bacterias después de t horas.

(b) Encuentre el número de bacterias después de 3 horas.

(c) Encuentre la rapidez de crecimiento después de 3 horas.
(d) ¿Cuándo llegará a 10000 células la población?
Solución.
Sea y(t) el número de bacterias en un instante t. Puesto que el crecimiento de la población es proporcional a su tamaño,
entonces
dy
= ky,
dt
y ya vimos que la solución de dicha ecuación es y(t) = y0 ekt = 100ekt . Ahora, como y(1) = 420, tenemos que

100ek = 420 =⇒ k = ln(4,2).

(a) Ası́, el número de bacterias después de t horas es y(t) = 100eln(4,2)t = 100(4,2)t .

(b) Tenemos y(3) = 100(4,2)3 ≈ 7409. Después de 3 horas habrá aproximadamente 7409 bacterias.
dy
(c) La rapidez de crecimiento en cualquier instante t está dada por = ln(4,2)y(t), por tanto, después de 3 horas será
dt
dy
= ln(4,2)y(3) = 7409 ln(4,2) ≈ 10633 bacterias/hora.
dt

ln(100)
(d) Debemos hallar t tal que 100(4,2)t = 10000. Ası́, t = ≈ 3,2. La población llegará a 10000 en 3 horas y 12
ln(4,2)
minutos aproximadamente.

3
Capı́tulo 7. Sección 7.1 [ clase 18 ] Cálculo Integral
Ecuaciones diferenciales y modelados Escuela de Matemáticas

Ejemplo 2.
El torio 234 decae con una rapidez proporcional a la cantidad presente del mismo. Si 100 mg de este material se reducen
a 82.04 mg en una semana:
(a) Encuentre una expresión para la cantidad presente en cualquier instante t.

(b) Encuentre el intervalo de tiempo que debe transcurrir para que la masa decaiga hasta la mitad de su valor original
(vida media).
Solución.
Sea Q(t) la cantidad de torio (en mg) en el instante t (en dı́as). Puesto que el material decae proporcional a la cantidad
presente
dQ
= kQ,
dt
donode k < 0. Y como antes la solución de esta ecuación está dada por Q(t) = Q0 ekt = 100ekt . Además, Q(7) = 82,04
entonces
ln(0,8204)
100e7k = 82,04 =⇒ k = = −0,02828.
7
(a) La cantidad presente en cualquier instante t es Q(t) = 100e−0,02828t mg.
Q(0)
(b) La vida media es el tiempo t para que Q(t) = , en este caso
2
 
1 ln 2
100e−0,02828t = 50 =⇒ −0,02828t = ln =⇒ t = ≈ 24,5 dı́as.
2 0,02828

Problema 3. Ley de enfriamiento de Newton.

La ley de enfriamiento de Newton dice que la rapidez de enfriamiento de un objeto es proporcional a la diferencia de
temperatura entre el objeto y su entorno, siempre que esta diferencia no sea demasiado grande. (Esta ley también se aplica al
calentamiento). Ası́
dT
= k(T − Ts ),
dt
donde T (t) es la temperatura del objeto en el tiempo t, Ts es la temperatura del entorno, y k es una constante (de
proporcionalidad).

4
Capı́tulo 7. Sección 7.1 [ clase 18 ] Cálculo Integral
Ecuaciones diferenciales y modelados Escuela de Matemáticas

Ejemplo
En la investigación de un asesinato, la temperatura del cuerpo era de 32,5◦ C a la 1:30 p.m. y 30,3◦ C una hora des-
pués. La temperatura normal del cuerpo es 37◦ C y la temperatura del entorno era de 20◦ C. ¿A qué hora aconteció el crimen?

Solución.
Sean T (t) y Ts definidas como antes. En este caso T (0) = 32,5, T (1) = 30,30 y Ts = 20. Ası́

dT
= k(T − 20).
dt
De forma completamente análoga a lo hecho anteriormente podemos ver que la solución de esta ecuación es

T (t) = 20 + Cekt .

Ahora, T (0) = 32,5 =⇒ C = 12,5 y con T (1) = 30,3 tenemos

 
k k 10,3 10,3
20 + 12,5e = 30,3 =⇒ e = =⇒ k = ln = −0,19358.
12,5 12,5

Esto nos dice que la temperatura del cuerpo en cualquier instante t está dada por

T (t) = 20 + 12,5e−0,19358t .

Nos pregunta por t cuando T (t) = 37◦ C, ası́

 
17 17
20 + 12,5e−0,19358t = 37 =⇒ e−0,19358t = =⇒ −0,19358t = ln =⇒ t ≈ −1,6,
12,5 12,5

esto nos dice que el crimen tuvo lugar 1 hora y 36 minutos antes de la 1:30 p.m., es decir, a las 11:54 a.m.