APLICACIONES DE LAS ECUACIONES

DIFERENCIALES DE PRIMER ORDEN

INTRODUCCIÓN

Existe un gran número de aplicaciones de las ecuaciones diferenciales en diversas áreas

del conocimiento, ya sea en las ciencias exactas, e incluso en las ciencias sociales, en donde lo
determinístico pierde fuerza para dar paso a lo probabilístico; pero a pesar de ese aspecto, las
ecuaciones diferenciales en su más reciente etapa de evolución, han tratado de resolver problemas
de este tipo y con mucho éxito. El análisis cualitativo de las ecuaciones diferenciales ha sido
desde la mitad del siglo XIX hasta la fecha, una poderosa herramienta para tratar de comprender
como se comportan las variables en un proceso, ya sea determinístico o estocástico.

En la presente sección, se mostrarán algunas de las más comunes aplicaciones que se

imparten en el desarrollo de los primeros semestres de ingeniería, en cuanto a matemática
aplicada se refiere. Es claro que la mayoría de estas aplicaciones son modelos de baja resolución
(más adelante se explicará con detalle este término) y solo representarán una aproximación ó
idealización de un cierto proceso ya sea de la naturaleza física o humana.

Como fenómenos físicos se deben entender los procesos que ocurren en la naturaleza
física que rodea al ser humano, como los huracanes, las mareas, la caída de los cuerpos, la
fotosíntesis, el enfriamiento de los cuerpos, la difusividad, etc. Estos fenómenos pueden existir
sin la intervención humana; por ejemplo la tormenta en espiral de la mancha roja del planeta
Júpiter, o los feroces vientos de un valle en el planeta Marte. Por otro lado como fenómeno
humano se debe entender como el proceso que ocurre con la intervención directa o indirecta del
ser humano y la naturaleza física; es difícil concebir un fenómeno humano desligado totalmente
de la naturaleza física; hasta los proceso neuronales que ocurren en el cerebro humano son
afectados por el entorno físico. Otro ejemplo es el crecimiento demográfico que definitivamente
esta ligado con el entorno físico, como el clima, la alimentación, etc. En adelante se considerará
como fenómenos de la naturaleza, un fenómeno físico, un fenómeno humano o una combinación
de ambos.

Ya en el apartado sobre los modelos matemáticos se trató la importancia que han tenido
éstos a través de la historia del hombre. El estudio ahora se centra en la modelación de ciertos
fenómenos que ocurren en la naturaleza física y en el fenómeno humano.

¿CÓMO SE MODELA UN FENÓMENO DE LA NATURALEZA?

La descripción matemática de un sistema o fenómeno, llamado modelo matemático, se

forma con ciertos objetivos en mente; por ejemplo, tratar de comprender cómo varía la
temperatura de una aleta de enfriamiento, intentar fechar fósiles analizando la desintegración de
una sustancia radiactiva, ya sea en el fósil o en el estrato donde se encontraba.

La formulación de un modelo matemático de un fenómeno se inicia :

a) Mediante la identificación de las variables causantes del las variaciones que sufre el sistema.
Podremos elegir no incorporar todas las variables en el modelo desde el comienzo. En ese
paso especificamos el nivel de resolución del modelo.

Cabe señalar que si se elige erróneamente a las variables que afectan significativamente el
sistema, se logra un modelo de muy baja resolución y no serviría de nada obtenerlo.
b) Se establece un conjunto de hipótesis razonables acerca del sistema que se trata de describir.
Esas hipótesis también incluyen todas las leyes empíricas aplicables al sistema.

Y la obtención del modelo matemático final depende de cualquiera de las siguientes

circunstancias

c) Si la hipótesis del modelo se establece mediante ecuaciones matemáticas, éstas se resuelven

cuantitativa o cualitativamente, junto con las condiciones necesarias para tal efecto. Por
ejemplo una hipótesis puede ser una ecuación diferencial o un sistema de ecuaciones
diferenciales y las condiciones necesarias serían las condiciones iniciales o de Cauchy

d) Si la hipótesis del modelo es un algoritmo numérico, éste se aplica junto con las condiciones
necesarias para la obtención del modelo matemático, tal es el caso de los fenómenos de
propagación de ondas.

La presente obra, se circunscribe únicamente a la parte c), es decir, las hipótesis se establecen
como ecuaciones diferenciales hasta el orden n inclusive o sistemas de ecuaciones diferenciales .

El siguiente esquema (Fig. 8) muestra la formulación de un modelo considerando como hipótesis

una ecuación diferencial de primer orden:

Se formula la hipótesis del Se consideran las

Se toman en consideración modelo con la ecuación condiciones iniciales del
las variables que afectan diferencial proceso mediante
significativamente el proceso
Φ ( x , y , y' ) = 0 y (x0 ) = y 0

Reformulación de la
Hipótesis

Confrontación del modelo Se consideran datos

Aplicación de la metodología
con los hechos reales. Si es adicionales posteriores al
adecuada para la obtención
satisfactorio, el modelo es inicio del proceso,
del modelo de la forma
definitivo, si no es así, se generalmente dados por la
reformula la hipótesis ϕ( x , y ) = 0 experiencia.

Para algunos fines quizá baste contar con modelos de baja resolución; por ejemplo, en los
cursos básicos de física el lector habrá advertido que al modelar el movimiento de un cuerpo que
cae sobre la superficie de la Tierra, se hace caso omiso de la resistencia del aire y algunas otras
influencias colaterales como vientos cruzados etc. Pero si usted es un científico cuyo objeto es
predecir con exactitud la trayectoria de vuelo de un proyectil de largo alcance, deberá tener en
cuenta la resistencia del aire, esos vientos cruzados además de la curvatura de la Tierra, etc.
 Trayectorias Ortogonales
Se dice que dos curvas son ortogonales si se interceptan y en los puntos de corte
sus rectas tangentes son perpendiculares entre sí. Si todas las curvas de una familia de
curvas , son ortogonales a todas las curvas de otra familia
, entonces se dice que las familias son cada una, trayectorias
ortogonales de la otra.

Es una familia de curvas cuyas pendientes son perpendiculares entre sí, o de otra
manera son las curvas que se interceptan formando ángulos rectos. Si una familia de
curvas tiene la ecuación:

La ecuación diferencial de las trayectorias ortogonales a ella, es otra familia de

la forma:

Para obtener las trayectorias ortogonales de una ecuación diferencial, se toma:

Y como:
 Trayectorias ortogonales y sus aplicaciones

Suponga que nos dan una familia de curvas como en la Figura 3.16 (líneas
Gruesas). Podemos pensar en otra familia de curvas (líneas puxiteadas). Tal que cada
miembro de esta familia corte a cada miembro de la primera familia en ángulos rectos.
Por ejemplo, la curva AB se encuentra con varios Miembros de la familia punteada en
ángulos rectos en los puntos L, M, N, O, P. Decimos que las ‘familias son mutuamente
ortogonales, o que una familia Forma un conjunto trayectorias ortogonales de la otra
familia.

Figura 3.16

Considere la familia de todos los círculos con centro en el origen; Unos cuantos
círculos aparecen en, la Figura 3.17. Las trayectorias ortogonales Para esta familia de
círculos podrían ser miembros de la familia de las líneas rectas (líneas punteadas).
Similarmente las trayectorias ortogonales de la familia de líneas rectas que pasan por el
origen son los círculos con centro en el origen.
 Figura 3.17

Como una situación más complicada, considere la familia de elipses Figura 3.18
y la familia de curvas ortogonales a ellas. Las curvas de una familia Son las trayectorias
ortogonales de la otra familia. Las aplicaciones de Trayectorias ortogonales son
numerosas en física e ingeniería.

Figura 3.18
 Como una aplicación muy elemental, considere la Figura 3.19. Aquí NS
representa una barra Magnética, siendo N su polo norte, y S su polo sur. Si sus
limaduras de hierro se esparcen alrededor del magneto encontramos que ellas se
ordenan así mismas Como las curvas punteadas de la Figura 3.19. Estas curvas se
llaman Líneas de fuerza.* Las curvas perpendiculares a estas (líneas gruesas) se llaman
líneas equipotenciales, o curvas de igual potencial. Aquí, también los miembros de una
familia constituyen las trayectorias ortogonales de la otra familia.

Figura 3.19

Como otro ejemplo de física considere la Figura 3.20, la cual representa Un

mapa del clima tan familiar en muchos de nuestros periódicos diarios. Las curvas
representan isobaras, las cuales son curvas que conectan todas Las ciudades que
reportan la misma presión barométrica a la oficina metereológica. Las trayectorias
ortogonales de la familia de isobaras podrían indicar La dirección general del viento
desde áreas de alta a baja presión. En vez de Isobaras, la Figura 3.20 podría representar
curvas isotérmicas las cuales son Curvas que conectan puntos que tienen la misma
temperatura.

Figura 3.20
Ejemplos:
1.) Encontrar las trayectorias de familias de parábolas

Solución: de la ecuación de la familia dada se sigue que:

Sustituyendo

Obtenemos

Luego las trayectorias ortogonales deben cumplir

Resolviendo la ecuación diferencial encontramos

Así las trayectorias ortogonales de las parábolas con vértice en el origen y cuyo
eje es el eje y, son elipses con centro en el origen y eje mayor en el eje x.

2.) Obtener la familia de curvas ortogonales a

Derivamos:

Pendiente

Se cumple:

Aplicando doble c
 Integrando

Aplicando exponencial

Se eliminan los exponenciales con los logaritmos naturales y queda:

Ò Familia de curvas ortogonales.

3.) Halle las trayectorias ortogonales de la familia de parábolas

De obtenemos:

Y como despejado

Nos queda que:

Buscamos

Integramos

Ò
Ejemplo 3
a s
Halle la familia de trayectorias ortogonales a las exponenciales

y = Cex . á t ic
t em
Derivando obtenemos:

M
y 0 = Cex
a
d 0
y =y e
t o
n
La ecuación que satiface la familia ortogonal es:
e
t am y0 =
−1
y

a r
p
De la EDO
Resolvemos

y0 = − 1 a s
ic
át
y

t em
ydy = −dx
1 2 M a
2
de
y = −x + C

toábolas
La familia ortogonal son las par

en x = C − y 

am  1 2
2 .

r t
p a
D e
as
y
ic
5
át
em
4

t
2
Ma
de x
to
−8 −6 −4 −2 0 2 4 6 8

en −2

am
r t −4

p a −5

De Gráfica de las exponenciales y = Cex .

Gráfica de las parábolas 1 2
2 y = −x + C.
MA2115 Matemáticas IV (semi-presencial) Práctica 05
8 / 21
N
Ejemplo 4
a s
tic
Halle la familia de trayectorias ortogonales a la familia de c´ırculos

á
y 2 + x 2 = 2Cx.

e m
Derivando obtenemos:
a t
2yy + 2x = 2CM
0

y + x d=e(2yy + 2x)x
2 2 0

t o = y −x
2xyy 0 2 2

en y = y 2xy− x 0
2 2

a m
r t satiface la familia ortogonal es:
a
La ecuación que
p
D e y = 0 −1
=
2xy
y 2 −x 2
x −y 2 2
2xy

a s
ic
Resolvemos esta ecuaci´on homogenea:

2xy 0 2 yx át
em
y = 2 =
x − y2 1 − y 2

a
x
t y0 0
M
Hacemos el cambio: y = zx donde = z + xz .

ze
2z
d 1−z
+ xz = 0
2

to xz = z + z 3

en 1 − z 11 − z
0
2

a m 2

r t z +z
dz = dx
x 3

p a 
1
−
2z

dz =
dx

D e z z +1 x 2
2
ln(z) − ln(z + 1) = ln(Ax)
 a s
ln(z) − ln(z 2 + 1) = ln(Ax) á tic

z

t em
ln 2
z +1
M a
= ln(Ax)
z
de 2
= Ax
z +1

to
y
x

en
= Ax
y
2
x +1

am
y
=A
r t y2 + x2

pa 2By = y 2 + x 2

De
La familia ortogonal son las circunferencias
 
2 2
x + (y − B) = B .2
 
MA2115 Matemáticas IV (semi-presencial) Pr ´
11 / 21

actica 05
N

a s
y
tic
5
á
4
e m
2
M at
d e x
to
−8 −6 −4 −2 0 2 4 6 8

e n −2

am
r t −4

p a −5

De Gráfica de las circunferencias x 2 + y 2 = 2Cx.

Gráfica de las circunferencias x 2MA2115 2 = 2By .
+ yMatemáticas IV (semi-presencial) Práctica 05
12 / 21
N
Ejemplo 5
a s
i c
Halle la familia de trayectorias ortogonales a la familia de curvas.
t
y + x = Ce .
m á y

Derivando obtenemos:
a t e
y + 1 = Ce yM
0 y 0

0
d
y + 1 = (ye+ x )y 0

t o1 = (y +1x − 1)y 0

e n y= 0

m y +x −1
tasatiface la familia ortogonal es:
r
La ecuación que
a
e p 0
y = 1−y −x
D s
Resolvemos esta ecuacion
´ diferencial:
i c a
0
y = 1−y −x á t
Hacemos el cambio: u = x + y donde u = 1 + ye
m
a 0
t . 0 0

e M
u −1 = 1−
u = 2−u
u
0

1d
t o2−u
du = dx

enln(2− ln(2 − u) = x + A

a m − x − y ) = −x + B
r t 2 − x − y = Ce −x

a
p ortogonal son las curvas
La familia

D e y = 2 − x − Ce .

−x
 

a s
y
ic
5
át
em
4

a t
2
M
de
MA2115 Matemáticas IV (semi-presencial) Práctica 05
14 / 21
N
x
to
−8 −6 −4 −2 0 2 4 6 8

en −2

a m
r t −4

p a −5

De Gráfica de las curvas x + y = Cey .

Gráfica de las curvas y = 2 − x + Ce−x .
MA2115 Matemáticas IV (semi-presencial) Práctica 05
15 / 21
N
ENFRIAMIENTO DE LOS CUERPOS

Según la Ley empírica del enfriamiento de Newton establece que:

“ La velocidad o razón de cambio de la temperatura de un

cuerpo es proporcional a la diferencia de la temperatura de
dicho cuerpo y la temperatura del medio ambiente que lo
rodea”

matemáticamente esto puede representarse de la siguiente forma:

dT
= k (T − T0 ) (2.33)
dt

donde T es la temperatura del cuerpo, T0 es la temperatura del medio que rodea al cuerpo, y k es
una constante de proporcionalidad

Ejemplo 1.- Un cuerpo se saca de un reservorio donde la temperatura es de 37 ºC y se lleva al

exterior donde la temperatura del ambiente es de 15 ºC. La temperatura de cuerpo es medida
instantáneamente con un termopar y después de ½ minuto la temperatura del cuerpo es de 25 ºC.
¿Cuál es la temperatura del cuerpo después de 1 minuto?; ¿Cuánto tiempo transcurrirá
aproximadamente para que la lectura del termopar marque 15 ºC?. Elabore una gráfica de este
fenómeno.

Solución.- Si la temperatura del ambiente es de 15 ºC, entonces T0 = 15, luego, la hipótesis del
modelo es:
dT
= k (T − 15)
dt

esta es una ecuación diferencial de primer orden y de variables separables.

Separando las variables e integrando se tiene:
dT
∫T − 15
= kdt ∫
calculando las integrales se obtiene:
ln(T − 15) = kt + C1
aplicando leyes de logaritmos:
T − 15 = e kt +C1
si e C1 = C y despejando T se obtiene:

T (t ) = 15 + Ce kt

la cual representa el modelo general para el fenómeno citado, que determina la temperatura para
cualquier tiempo t. Pero este modelo está incompleto; se tienen que calcular los parámetros C y k.

Para la determinación de C son necesarios los siguientes datos: T (0 ) = 37 , que es la temperatura

inicial del cuerpo antes de ser sacado del reservorio. Aplicando dichos datos al modelo general
deducido se tiene:

T (0) = 37 = 15 + Ce k (0 )

de donde se deduce que C = 22. Sustituyendo este valor en el modelo general, se tiene:

T (t ) = 15 + 22e kt
Pero todavía falta determinar el valor de k, para obtener el modelo matemático definitivo. Para
calcular el valor de k se necesitan datos adicionales. Según el problema, éstos son: T ( 12 ) = 25 , lo
cual significa que después de medio minuto la temperatura del cuerpo es de 25 ºC. Aplicando
estas condiciones en la última relación se tiene:
T ( 12 ) = 25 = 15 + 22e 2
1
k

de donde se deduce que:

1
k 10
e2 =
22
despejando k se obtiene:
k = 2ln 10
22

lo cual arroja un valor aproximado de k = −1.577. Sustituyendo este valor en la última expresión
para T, se tiene:
T (t ) = 15 + 22e −1.577t

la cual representa el modelo matemático para la variación de la temperatura del cuerpo en

cualquier momento t. Conociendo el modelo definitivo para este caso en particular, es posible
ahora contestar las preguntas plateadas en el problema.

¿Cuál es la temperatura del cuerpo después de 1 minuto?

Si se sustituye t = 1 en el modelo definitivo, se tiene: T (1) = 15 + 22e −1.577 (1)

simplificando se tiene que: T(1) = 19.54, lo cual significa que pasado un minuto la temperatura
del cuerpo es de 19.54 ºC
t
(min)
T(t) Por otro lado ¿Cuánto tiempo
(ºC)
transcurrirá aproximadamente
para que la lectura del
termopar marque 15 ºC?. Para
contestar a la pregunta
supóngase que T(t) = 15. Pero
solo es posible calcular un
tiempo aproximado a través
de la gráfica.

La figura 9 muestra la gráfica

de la variación de la
temperatura en función del
tiempo y una tabla de valores
Fig. 9
de T, en donde es posible
visualizar el comportamiento
asintótico de T(t), es decir la temperatura del cuerpo tiende a 15 ºC conforme aumenta el tiempo.
Como puede observarse en la gráfica, aproximadamente para t = 3 , la temperatura alcanza un
valor muy próximo a 15. esto es, pasados 3 minutos la temperatura del cuerpo habrá alcanzado un
valor muy próximo a 15 ºC.

Ejemplo 2.- Unas pinzas metálicas que inicialmente están a 23 ºC, se dejan caer en
una cubeta con agua hirviendo( Fig. 10). Si se sabe que su temperatura aumentó 3
ºC en un segundo. ¿Cuánto tiempo tardarán las pinzas en alcanzar una temperatura
de 80 ºC?. ¿Cuánto tiempo tardarán en llegar a 90 ºC?

Solución.- Si se considera que la temperatura del agua hirviendo es de 100 ºC,

entonces T0 = 100, luego, la hipótesis del modelo es:

dT
= k (T − 100 ) Fig. 10
dt

esta es una ecuación diferencial de primer orden y de variables separables y es además, similar
al obtenido en el problema anterior. El procedimiento para solucionar esta ecuación es
obviamente parecido al aplicado en el problema 1.

Resolviendo la ecuación diferencial anterior se obtiene:

T ( t ) = 100 + Ce kt
la cual representa el modelo general para el fenómeno citado, que determina la temperatura para
cualquier tiempo t. Pero este modelo está incompleto; se tienen que calcular los parámetros C y k.

Para la determinación de C son necesarios los siguientes datos: Si se sabe que las pinzas están
inicialmente a 23 ºC, entonces, T ( 0 ) = 23 Aplicando dichos datos al modelo general deducido se
tiene:
T ( 0 ) = 23 = 100 + Ce k ( 0)

de donde se deduce que C = −77. Sustituyendo este valor en el modelo general, se tiene:

T ( t ) = 100 − 77e kt

Para calcular el valor de k se necesitan datos adicionales. Acorde con el problema, si se sabe que
la temperatura de las pinzas aumentó 3ºC en un segundo, entonces T (1) = 26 . Aplicando estas
condiciones en la última relación se tiene:
T (1) = 26 = 100 − 77e k (1)
de donde se deduce que:
−74
ek =
−77
despejando k se obtiene:
k = ln 77
74

lo cual arroja un valor aproximado de k = −0.0397. Sustituyendo este valor en la última expresión
para la temperatura T, se tiene:

T ( t ) = 100 − 77e −0.0397 t

la cual representa el modelo matemático para la variación de la temperatura de las pinzas en

cualquier momento t. Conociendo el modelo definitivo para este caso en particular, es posible
ahora contestar las preguntas plateadas en el problema.

¿Cuánto tiempo tardarán las pinzas en alcanzar una temperatura de 80 ºC?, es decir ¿Cuál es el
tiempo t para el cual la temperatura T sea igual a 80 ºC?. Matemáticamente esto se representa de
la siguiente manera:
80 = 100 − 77e−0.0397 t

ahora se despeja el valor de t, el cual está dado por:

1 20
t=− ln
0.0397 77

simplificando los cálculos se tiene t = 33.9, es decir que aproximadamente a los 34 segundos las
pinzas habrán alcanzado la temperatura de 80 ºC
Por otro lado, ¿Cuánto tiempo tardarán en llegar a 90 ºC? Matemáticamente esto se representa de
la siguiente manera:
90 = 100 − 77e−0.0397 t

ahora se despeja el valor de t, el cual está dado por:

1 10
t=− ln
0.0397 77

simplificando los cálculos se tiene t = 51.41, es decir que aproximadamente a los 51.4 segundos
las pinzas habrán alcanzado la temperatura de 90 ºC

La figura 11 muestra la gráfica de la variación de la temperatura en función del tiempo y una

tabla de valores de T, en donde es posible visualizar el comportamiento asintótico creciente de
T(t), es decir la temperatura del cuerpo tiende a los 100 ºC conforme aumenta el tiempo. También
puede observarse en la tabla de valores, que pasados los 110 segundos las temperaturas son muy
semejantes entre sí, esto es debido al comportamiento asintótico del modelo matemático que
representa al fenómeno.

t T(t)
(seg) (ºC)

Fig. 11

Es claro que la hipótesis matemática (2.33), sirve como base tanto para fenómenos de
enfriamiento como de calentamiento, como lo hace constar los dos ejemplos anteriores, claro está
que los modelos encontrados son ideales, puesto que no considera la variación simultánea de la
temperatura del medio ambiente que en la realidad física, son variables y no constantes como se
supuso desde el principio. Empero si se toma ésta como constante, los modelos se pueden
considerar como aceptables.
 ◦
on 100 C, se retira del fuego y se mantiene en
Ejemplo 3. Se calienta agua hasta su punto de ebullici´
un recinto que se encuentra a una temperatura constante de 60◦ C. Al final de 3 minutos la temperatura
del agua es de 90◦ C. Calcular la temperatura del agua al final de 6 minutos y el tiempo en la cual la
temperatura del agua ser´a de 75◦ C.

Figura 3.10: Agua a 100◦ C.

Solución: Sea y = f (t) la temperatura de un cuerpo en un tiempo t; la velocidad de variación de y con

respecto al tiempo t es proporcional a la diferencia entre su temperatura y la del medio ambiente
(Ley de Enfriamiento de Newton). Sea M la temperatura del medio ambiente. Luego la ecuación
diferencial que sirve de modelo a este problema es y  = −k(y − M ) o equivalentemente

y  + ky = kM, (3.36)

donde k es una constante positiva, esta ecuación es diferencial lineal de primer orden, luego su
solución viene dada por
 t
−kt −kt
f (t) = f (a)e +e kM ekt dt, (3.37)
a

donde f (a) es la temperatura del cuerpo en un tiempo t = a. Para nuestro problema se tiene
f (0) = 100, la cual es la condición inicial dada, además M = 60. Reemplazando estas condiciones
en (3.37) se tiene
 t
f (t) = 100e−kt + e−kt k60ekx dx
a
= 100e−kt + e−kt [60ekt − 60]
= 100e−kt − 60e−kt + 60
= 40e−kt + 60. (3.38)

Hallamos k utilizando la condición de f (3) = 90. En efecto, evaluando f definido en (3.38) en t = 3

se obtiene
f (3) = 40e−3k + 60 = 90,
 
lo anterior implica que k = 13 ln 43 . Entonces la solución a la ecuación diferencial (3.36) es
1 4
f (t) = 40e− 3 ln( 3 )t + 60.

Hallemos la temperatura al final de 6 minutos, calculamos

1 4
f (6) = 40e− 3 ln( 3 )(6) + 60
4
= 40e−2 ln( 3 ) + 60
16
= 40e− ln( 9 ) + 60
45
= + 60 = 82,5
2
Por lo tanto, después de 6 minutos la temperatura será de 83◦ C aproximadaente. El instante en el
cual la temperatura es de 75◦ C, lo podemos calcular ası́:

75 = f (t) = 40e−(1/3) ln(4/3)t + 60,

luego

75 = 40e−(t/3) ln(4/3) + 60
15 = 40e(t/3) ln(4/3) ,

la solución de la ecuación anterior es t = 10 minutos.

 CONCENTRACIÓN EN TANQUES AGITADOS

Cuando dos soluciones con diferentes concentraciones se ponen en contacto,

sucede el proceso de difusión o transferencia de masa y llega un momento en
el que la mezcla al final, tiende a una concentración uniforme. La difusión
puede ser molecular, es decir sin que
intervengan agentes externos y la difusión por
remolino o difusión turbulenta en donde existe la
intervención de un dispositivo externo tal como un agitador, por
mencionar un ejemplo. En este apartado se tratará el fenómeno del
mezclado de soluciones salinas o salmueras de diferente
concentraciones, en tanques agitados. Pero antes se deben mostrar
ciertos conceptos teóricos preliminares. El flujo volumétrico Q de una corriente de fluido que
V
viaja a través de una tubería está dado por Q = donde V es el volumen de fluido que fluye en
t
un tiempo t. La concentración C de una solución es la relación entre la masa del soluto y el
M
volumen del disolvente, esto es C = , donde en este caso M es la masa del soluto y V el
V
volumen del disolvente. Por otro lado, el flujo másico o flujo masivo m es la relación de la masa
M
que fluye en un intervalo de tiempo determinado, esto es m = , donde M es la masa del fluido
t
que fluye en el tiempo t. Es claro que existe una relación entre el flujo volumétrico y la
concentración de una solución. Esta relación está dada por el flujo másico, esto es:

⎛ V ⎞⎛ M ⎞
m = ⎜ ⎟⎜ ⎟ = QC (2.34)
⎝ t ⎠⎝ V ⎠
Q1
Supóngase que se tiene un tanque agitado (Fig. 12) con un volumen
C1
de agua V0 en donde se han disuelto M0 unidades de masa de sal.
Supóngase también que al tanque, entra una solución de una
salmuera de concentración C1 y flujo volumétrico Q1 y sale la V0
solución mezclada con una concentración C2 y con un flujo
volumétrico Q2. Si se aplica el principio de la conservación de la
materia 1, es decir si se hace un balance de materia al sistema, se
puede determinar la velocidad con la que varía la cantidad de sal
M(t) en cualquier momento t, en términos de los flujos másicos de C2
entrada y de salida. El balance matemáticamente puede ser Q2
expresado de la siguiente forma: Fig. 12

dM ⎛ flujo masico de ⎞ ⎛ flujo masico de ⎞

=⎜ ⎟−⎜ ⎟ = m 1 − m 2 (2.35)
dt ⎝ entrada ⎠ ⎝ salida ⎠
y en términos de las concentraciones y de los flujos volumétricos:

dM
= C1Q1 − C2Q2 (2.36)
dt
donde C2 es función de M.
M (t )
Esta dependencia está dada por C2 = , de tal manera que la ecuación (2.35) queda
V0
finalmente como:
dM M
= C1Q1 − Q2 (2.37)
dt V0

La ecuación (2.37) que es una ecuación diferencial lineal de primer orden, representa la hipótesis
del fenómeno de la mezcla de dos soluciones en un tanque agitado, que al resolverla, se obtiene el
modelo que determina la cantidad de sal M, que hay en el tanque en cualquier momento t. Esta
ecuación diferencial puede estar sujeta a la condición inicial dada por: M ( 0 ) = M 0

Es interesante comprender que si Q1 = Q2, entonces, el nivel del líquido en el tanque permanece
invariable, pero si Q1 > Q2, entonces el nivel del líquido aumenta; por otro lado si Q1 < Q2,
entonces el nivel del líquido en el tanque disminuye.

Ejemplo 3.- En un tanque que contiene 200 litros de agua, se

disuelven 100 gramos de sal y dicho tanque es agitado 5 lts/min
continuamente (Fig. 13). Al mismo tiempo, se introducen en 2 gr/lto
el tanque 5 lts/min de una solución cuya concentración es de
2 gr/lto; del tanque se bombea la solución mezclada al
mismo flujo volumétrico. Encuentre una expresión 200 lts
matemática que determine la cantidad de sal en el tanque en
cualquier momento t. ¿Cuál es la cantidad de sal en el tanque
pasado mucho tiempo?

Solución.- Se determinan los datos iniciales V0 = 200 y M(t)/200

M0 = 100. Luego se determinan los flujos másicos de entrada 5 lts/min
y de salida. El flujo másico de entrada esta dado por:
Fig. 13

⎛ 2 gr ⎞ ⎛ 5lto ⎞ gr
m 1 = ⎜ ⎟⎜ ⎟ = 10
⎝ lto ⎠ ⎝ min ⎠ min

El flujo másico de salida está dado por.

⎛ M ( t ) gr ⎞ ⎛ 5lto ⎞ M ( t ) gr
m 2 = ⎜ ⎟⎜ ⎟=
⎝ 200lto ⎠ ⎝ min ⎠ 40 min

Aplicando la fórmula (2.35) se tiene:

 dM M
= 10 −
dt 40

dM M
reacomodando los términos: + = 10
dt 40

esta es una ecuación diferencial lineal de primer orden, que representa la hipótesis del modelo
matemático que determinará la cantidad de sal que hay en el tanque en cualquier momento t.
Aplicando la metodología para resolver ecuaciones diferenciales de primer orden se tiene:

M ( t ) = e 40 ∫ ∫ 10e 40 ∫ dt + ce 40 ∫
− 1 dt 1 dt − 1 dt

simplificando se tiene:
M ( t ) = 400 + ce
− 40
1 t

esta integral general, representa el modelo general para el caso dado; solo falta determinar la
constante c, la cual se obtiene aplicando la condición de Cauchy M ( 0 ) = M 0 = 100 . Aplicando
( )
esta condición a la integral general se obtiene: M ( 0 ) = 100 = 400 + ce
− 40
1 0

de donde se deduce que c = −300. Finalmente el modelo queda determinado por la ecuación:

M ( t ) = 400 − 300e
− 40
1 t

el cual representa la cantidad de sal que hay en el tanque en cualquier momento t.

T M(t)
(min) (gr)

Fig. 14
 La figura 14 muestra el comportamiento asintótico de la función que determina la
cantidad de sal en el tanque en cualquier momento t.

Esta cantidad tiende a un valor de 400 gr y su explicación es muy simple. El flujo que entra en el
tanque desplaza a la cantidad de sal que había inicialmente y después de mucho tiempo solo esta
presente un flujo con una concentración de 2 gr/lto, que al multiplicarlo por 200 lts que es el
volumen constante del líquido en el tanque, resultan los 400 gr de sal presentes en el tanque. La
solución en ese instante tendría precisamente la concentración del flujo de entrada.

Ejemplo 4.- En un tanque que contiene 100 litros de agua, 4 lts/min

se disuelven 150 gramos de sal y dicho tanque es agitado
continuamente (Fig. 15). Al mismo tiempo, se introducen en 2 gr/lto

el tanque 4 lts/min de una solución cuya concentración es de

2 gr/lto; Del tanque se bombea la solución mezclada a una 100 lts
velocidad de 2 lts/min. Encuentre una expresión matemática
que determine la cantidad de sal en el tanque en cualquier
momento t. ¿Cuál es la cantidad de sal en el tanque pasado
mucho tiempo?
M(t)/100
Solución.- Se determinan los datos iniciales V0 = 100 y
2 lts/min
M0 = 150. Luego se determinan los flujos másicos de
entrada y de salida. El flujo másico de entrada esta dado Fig. 15
por:

⎛ 2 gr ⎞ ⎛ 4lto ⎞ gr
m 1 = ⎜ ⎟⎜ ⎟=8
⎝ lto ⎠ ⎝ min ⎠ min

El flujo másico de salida es diferente en este caso, puesto que los flujos volumétricos de entrada y
de salida no son los mismos. El flujo volumétrico neto en el tanque es de (4 − 2) lts/min, es decir
2 lts/min; en otras palabras el líquido dentro del tanque aumentará 2t litros cada minuto. El flujo
másico ahora está dado por:

⎛ M ( t ) gr ⎞ ⎛ 2lto ⎞ M ( t ) gr
m 2 = ⎜ ⎟⎜ ⎟=
⎝ 100lto + 2t ⎠ ⎝ min ⎠ 50 + t min

Aplicando la fórmula (2.35) se tiene:

dM M
= 8−
dt 50 + t
reacomodando los términos:
dM M
+ =8
dt 50 + t

que nuevamente es una ecuación diferencial lineal de primer orden, que representa la hipótesis
del modelo matemático que determinará la cantidad de sal que hay en el tanque en cualquier
momento t.
Aplicando la metodología para resolver ecuaciones diferenciales de primer orden se tiene:

dt dt dt
− ∫ 50+t ∫ 50+t − ∫ 50+t
M (t ) = e ∫ 8e dt + ce
ó
8 c
M (t ) = ∫ ( 50 + t ) dt +
50 + t 50 + t
simplificando:
c
M ( t ) = 200 + 4t +
50 + t

esta integral general, representa el modelo general para el caso dado; solo falta determinar la
constante c, la cual se obtiene aplicando la condición de Cauchy M ( 0 ) = M 0 = 150 . Aplicando
esta condición a la integral general se obtiene:
c
M ( 0 ) = 150 = 200 + 4 ( 0 ) +
50 + 0

de donde se deduce que c = −2500. Finalmente el modelo queda determinado por la ecuación:

2500
M ( t ) = 200 + 4t −
50 + t

el cual representa la cantidad de sal que hay en el tanque en cualquier momento t.

t M(t)
(min) (gr)

Fig. 16

En la figura 16 se muestran los valores de la cantidad de sal M(t) en cualquier momento t, dada
por la integral particular que se dedujo anteriormente así como su gráfica. Como se puede
observar en la tabla, los valores de M(t) crecen monótonamente sin una tendencia asintótica. La
gráfica de la función M(t) parece una línea recta y eso es debido al término − 2500
50 + t que tiende a
cero rápidamente cuando t crece monótonamente; así que para valores grandes de t, solo la
expresión 200 + 4t es significativa.

Los modelos matemáticos que se encontraron en estos dos ejemplos, es decir, las
expresiones matemáticas que determinan la cantidad de sal en cualquier instante en el tiempo, son
de baja resolución, ya que no se tomaron en cuenta factores como: la presión en las tuberías de
entrada y de salida, los factores de fricción en las orillas de los agujeros del tanque por donde
están conectadas las tuberías de entrada y de salida, los efectos gravitacionales, los coeficientes
de transferencia de masa de ambas soluciones, la viscosidad de ambas soluciones, la temperatura
de operación, etc. Empero a pesar de estas omisiones, los modelos se aproximan a una realidad
física salvo en el ejemplo 4.

La respuesta o salida del sistema, es decir el modelo matemático del ejemplo 4 fue
producto de las condiciones iniciales y obedece a esta circunstancia, pero contrasta con lo que
podría suceder en la realidad, ya que el tanque tiene una capacidad finita de almacenamiento y
ésta sería rebasada por el volumen creciente del líquido que entra. En el caso de que el tanque
esté abierto en la parte superior, el líquido se derramaría por esa parte, pero si el tanque solo tiene
los agujeros de entrada y salida, la bomba que transporta el líquido hacia el tanque se quemaría.

VACIADO DE UN TANQUE
En 1644 el matemático italiano Evangelista Torricelli (1608-1647) publicó
“De Motum Gravium”, artículo que contiene lo que se conoce como la Ley o
Teorema de Toriccelli la cual establece que:

“La velocidad con que el líquido sale o es drenado por un orificio

practicado en pared plana y delgada en un tanque es igual a la que
adquiriría un cuerpo cayendo libremente desde el nivel del líquido
al centro del orificio”

Lo anterior quiere decir que la velocidad de salida del líquido mencionado a una
profundidad h con relación al fondo del tanque, es v = 2 gh , donde g es la
aceleración debida a la gravedad. Supóngase que se
AL
tiene un tanque lleno de un líquido (Fig. 17) y éste se
deja drenar a través de un agujero por la acción de la
gravedad. Se desea encontrar la profundidad h del
h líquido en cualquier instante t.

A0
Si el área de la sección transversal del orificio es Ao y la velocidad de
salida del líquido es v = 2 gh , entonces el flujo volumétrico o el
Fig. 17 volumen del líquido por unidad de tiempo es A0 2 gh . Por lo tanto la
velocidad o razón de cambio del volumen del líquido que sale por el orificio esta dado por la
siguiente ecuación diferencial:
dV
= − A0 2 gh (2.38)
dt

donde el signo menos indica que el volumen del líquido esta disminuyendo con el tiempo. Si el
volumen del líquido en el tanque es V ( t ) = AL h , donde AL es el área del espejo del líquido (en
dV dh
este caso es constante) entonces = AL . Sustituyendo esta última expresión en (2.38) y
dt dt
reacomodando los términos se tiene:

dh A
=− 0 2 gh (2.39)
dt AL

Esta hipótesis correspondiente al fenómeno del vaciado del tanque, está restringida a un área AL
constante del espejo del líquido y al hecho de que no se consideran factores de fricción el las
orillas del orificio, los cuales pueden reducir el flujo volumétrico de salida.

La relación (2.39) es una ecuación diferencial de variables separables que puede estar sujeta a la
condición de que h ( 0 ) = h0 , donde h0 representa la altura inicial del líquido antes de comenzar el
drenado.

Ejemplo 5.- Un tanque con forma de cilindro circular recto, de radio interior r = 36 cm y con una
capacidad de 250 lts. contiene 200 lts. de alcohol etílico, el cual es drenado por un orificio de 2
cm de radio. Encuentre una expresión matemática que determine la altura del líquido drenado en
cualquier instante t.

Solución.- Con los datos proporcionados se puede calcular la altura inicial h0 del alcohol etílico
en el tanque antes del drenado. Si el volumen del líquido forma un cilindro circular recto, su
volumen está dado por V = π r 2 h0 ,
V
Luego, la altura inicial del líquido es h0 = 2 . Si 200 lts. = 200,000 cm3 , entonces,
πr
3
200000cm
h0 = , lo cual da un valor de h0 = 49.12 cm. Si el espejo del líquido forma un círculo
π ( 36cm )
2

perfecto entonces el área de éste esta dada por AL = π r 2 ; introduciendo los datos, el área del
espejo está dada por AL = π ( 36cm ) , lo cual arroja un valor de AL = 4,071.5 cm2. Por otro lado el
2

área del orificio, considerándolo también como un círculo perfecto, está dada por A0 = π r0 2 . Si el
radio del agujero r0 es de 2 cm, entonces A0 = π ( 2cm ) , lo cual arroja un valor de A0 = 12.56
2

cm2.
Teniendo ya todos estos datos, es posible aplicar la ecuación (2.39). Considerando g = 981
cm/seg2 se tiene:
dh 12.56
=− 2 ( 981) h
dt 4071.5
la cual se reduce a:
dh
= −0.1366 h
dt

esta es una ecuación diferencial de variables separables sujeta a la condición de que

h ( 0 ) = 49.12 y representa la hipótesis del modelo matemático que determinará la altura h del
alcohol etílico drenado en cualquier momento t. Separando las variables e integrando se tiene:

dh
∫ h
= −0.1366 ∫ dt

simplificando:
2 h = −0.1366t + C

esta es la integral general o modelo general para el caso dado. Ahora, si se aplica la condición
inicial h ( 0 ) = 49.12 , se tiene: 2 49.12 = −0.1366 ( 0 ) + C . De donde se deduce que C = 14.017 ,
sustituyendo este valor en la integral general anteriormente deducida, entonces, el modelo
matemático final que determina la profundidad del alcohol etílico en cualquier instante es:
2 h = 14.017 − 0.1366t
o en forma explícita:
h ( t ) = 14 (14.017 − 0.1366t )
2

t h(t)
(seg) (cm)

Fig. 18

Como una cuestión adicional, ¿Cuándo se vaciará totalmente el tanque?. Para responder esta
cuestión basta con trabajar la integral particular, sustituir h = 0 y despejar el valor de t deseado.
Esto es:
2 0 − 14.017
t=
−0.1366

el cual arroja un valor de t = 102.6 seg ó 1.71 min, prácticamente, en menos de 2 minutos el
tanque se habrá drenado totalmente.

La Fig. 18 muestra la gráfica y una tabla de valores de el modelo matemático anteriormente

deducido, que determina la altura o profundidad del alcohol etílico, cuando éste es drenado en un
tanque de almacenamiento. Como puede observarse en la tabla, los valores de h tienden a cero
conforme pasa el tiempo y analíticamente se comprobó que exactamente a los 102.6 seg. el
tanque se habrá vaciado en su totalidad.

Empero la expresión matemática del modelo si no es acotada, entonces podría arrojar valores
mayores a cero después de los 102.6 seg. transcurridos; esto definitivamente no tendría sentido
físico, ya que en esos momentos no habría líquido en el tanque. Esto deja patente la importancia
del acotamiento de una función que representa un fenómeno en particular.

como analíticamente se demostró, h llega a cero exactamente cuando t = 22. Es interesante

observar como los valores de h parecen decrecer rápidamente cuando t > 18 seg.
 CAIDA LIBRE DE LOS CUERPOS

En la dinámica de los cuerpos rígidos se estudia el movimiento de un

objeto que se encuentra dentro de un campo gravitacional, cuando es soltado
o lanzado desde una altura determinada. Generalmente surge la necesidad de
saber cuál es la velocidad de caída y cuanto tiempo tarda en caer el objeto al
suelo.

Supóngase que se desea saber cuál es la velocidad de caída de un cuerpo en el aire en función de
su altura referenciada a la posición inicial del cuerpo. Primeramente se consideran los siguientes
supuestos:

a) La causa del movimiento es el propio peso del cuerpo

b) Se desprecia toda resistencia producida por el aire u otros objetos
c) No se toma en cuenta la forma del cuerpo, es decir, se considera tal cuerpo como masa
puntual.
d) La velocidad inicial del cuerpo es cero, es decir el cuerpo cae libremente
 e) La altura es referenciada a la posición inicial del cuerpo y no al nivel del suelo.

Bajo estas condiciones, puede suponerse que las fuerzas involucradas son proporcionales a la
aceleración que sufre el cuerpo en su caída; luego entonces puede ser aplicada la segunda ley de
Newton del movimiento de los cuerpos:
∑ Fx = max (2.40)

Aquí, x representará la distancia vertical entre la posición inicial (donde comienza el

movimiento) y una posición final del mismo después de un tiempo t, tomando como referencia a
la misma posición inicial; La hipótesis del modelo en forma diferencial adquiere la forma:

W dv
∑F x = max = v =W
g dx
(2.41)

simplificando se obtiene:
dv
v =g (2.42)
dx

Donde g es la aceleración de la gravedad, W es el peso del cuerpo y v es la velocidad del mismo.

dv
Aquí se toma la aceleración como a = v .
dx

La ecuación (2.41), representa la hipótesis del modelo matemático para la caída libre de un
cuerpo bajo las condiciones antes descritas. Al resolver esta ecuación diferencial se encuentra el
modelo matemático que mide la velocidad de cuerpo en caída libre en función de la altura
referenciada a la posición inicial del mismo. Es interesante observar que la ecuación diferencial
(2.42) es independiente del peso del cuerpo y en consecuencia, la velocidad de caída de un
cuerpo es independiente del tamaño y forma del cuerpo. En otras palabras, la velocidad es la
misma para dos cuerpos totalmente diferentes en tamaño y forma arrojados desde la misma
altura. Claro está que esta conclusión está basada en los supuestos antes descritos.

A continuación se mostrará el caso en el que es importante considerar el factor de forma de los

cuerpos así como la resistencia del aire; es decir, se tomarán en cuenta los siguientes supuestos:

a) La causa del movimiento sigue siendo el peso del cuerpo

b) Se tomará en cuenta la resistencia del aire de la siguiente manera: Se utilizará la relación
1
empírica R = c x ρSv 2 como la fuerza de resistencia al movimiento; donde cx es un
2
coeficiente adimensional de resistencia que depende de la forma del cuerpo. Por ejemplo
para un paracaídas cx = 1.4, para una esfera cx = 0.5 etc., ρ es la densidad del medio
ambiente (aire en este caso) y S es la superficie de la proyección del cuerpo sobre el plano
perpendicular a la dirección del movimiento.
c) Nuevamente, la velocidad inicial del cuerpo es cero, es decir el cuerpo cae libremente.

Aplicando nuevamente la segunda ley de Newton al sistema se tiene:

 1
∑Fx = W − R = W − cx ρ Sv 2
2
(2.43)

La nueva hipótesis del modelo en su forma diferencial adquiere la forma:

W dv 1
v = W − cx ρ Sv 2 (2.44)
g dx 2

donde g es la aceleración de la gravedad, cx es el coeficiente de resistencia, ρ es la densidad del

aire, S es es la superficie de la proyección del cuerpo sobre el plano perpendicular a la dirección
del movimiento, W es el peso del cuerpo y v es la velocidad del mismo. Aquí también, se toma la
dv
aceleración como a = v .
dx

Ejemplo 7.- Una tuerca de acero que se encontraba en una de las

partes más altas de un evaporador (Fig. 22) de una refinería,
accidentalmente cae desde una altura de 30 m. Encuentre una
expresión matemática que determine la velocidad de la tuerca en
función de la distancia entre el punto de partida hasta cualquier
punto de la trayectoria de caída. ¿Cuál será la velocidad cuando la
tuerca se encuentra a 2 metros por encima del suelo?
30 m
Solución.- Aplicando directamente la relación (2.42) se tiene:

dv
v = 9.81
dx

donde v es la velocidad de caída, g es la aceleración de la gravedad Fig. 22

y x es la distancia medida desde el punto de partida. Esta es una
ecuación diferencial de variables separables sujeta a la condición de v ( 0 ) = 0 , ya que el punto
de partida o de referencia es el punto desde donde parte la tuerca. Separando las variables e
integrando se obtiene:
∫ vdv = ∫ gdx
calculando las integrales se tiene:
1
2 v 2 = 9.81x + C

esta integral general corresponde al modelo general para el caso dado. Aplicando la condición
v ( 0 ) = 0 se tiene:
( 0) = 9.81( 0 ) + C
1 2
2

de donde se deduce que C = 0. Sustituyendo este valor en la integral general anteriormente

deducida, entonces, el modelo matemático final que determina la velocidad de caída para
cualquier distancia x es:
 1
2 v 2 = 9.81x
o en forma explícita:
v ( x ) = 4.43 x

siempre y cuando la variable x se encuentre acotada dentro del intervalo de definición 0 ≤ x < 30 .
Este es el modelo que determina la velocidad de la tuerca en función de la distancia del punto de
partida hacia un punto cualquiera de la trayectoria de caída.

Ahora la cuestión ¿Cuál será la velocidad cuando la tuerca se encuentra a 2 metros por encima
del suelo? Debe contestarse tomando en cuenta que para x = 28 m la tuerca se encuentra a 2 m
por encima del suelo. Sustituyendo este valor en el modelo matemático se tiene v ( 28 ) = 4.43 28
el cual arroja un valor de v = 23.43 m/seg ; convirtiendo esta velocidad a kilómetros por hora, la
velocidad es de 84.37 km/h. Este resultado podría ser catastrófico para alguien que por descuido
se haya quitado el casco de seguridad.

x v(t)
(m) (m/seg)

Fig. 23

En la figura 23 se muestran la gráfica y una tabla de valores correspondiente al modelo

matemático que mide la velocidad de caída de la tuerca en función de la distancia referenciada al
punto de partida. Como puede observarse en la gráfica, los valores de la velocidad crecen
monótonamente hasta un valor aproximado de 24.26 m/seg, que representa la máxima velocidad
que adquiere la tuerca
Este es un modelo de muy baja resolución, puesto que en la realidad física los cuerpos que se
mueven a través del aire son sometidos a la resistencia que opone el mismo y, en altas
velocidades, dicha resistencia es considerable. La importancia del factor de forma es
verdaderamente relevante si se considera el diseño de los paracaídas o el diseño de aviones, más
no es así para el diseño de sondas espaciales.
 PROBLEMAS VARIADOS

Ejemplo 3.4.1. Una gota de aceite, cuya masa es de 0,2 g cae partiendo del reposo. Cuando su velocidad
es de 40 cm/s, la fuerza debida a la resistencia del aire es de 160 dinas. Suponiendo que dicha fuerza es
proporcional a la velocidad instantánea:
a) Hállense la velocidad y la distancia recorrida en función del tiempo.
b) Determinar la velocidad lı́mite.

Figura 3.12: Gota de aceite de masa 0, 2g.

Solución: La gota está afectada por dos fuerzas una debida a su peso w, la cual obliga a la gota a ir
hacia abajo y otra que se opone a esta y es la resistencia del aire r, la cual es proporcional a la
velocidad instantánea. Consideraremos positivo hacia abajo y negativo hacia arriba. Entonces según
la segunda ley de Newton tenemos F = w − kv, donde k es la constante de proporcionalidad y es
positiva, luego la ecuación diferencial correspondiente es mv  = mg − kv, siendo m la masa, v la
velocidad, g la gravedad, esta ecuación se puede transformar en v  + kvm = g. La cual es una ecuación
diferencial lineal de primer orden, su solución es
 t
−A(t) −A(t)
v(t) = v(0)e +e geA(x) dx, (3.45)
0
 donde  t
k k
A(t) = dx = t,
0 m m
luego
 t
kt kt kx
v(t) = v(0)e− m + e− m ge m dx
 0gm
gm 
− kt − kt
kt
= v(0)e m +e em −
m
k  k
kt gm  − kt
= v(0)e− m + 1−e m . (3.46)
k
Si v(0) = 0, como en nuestro caso entonces:
gm  kt

v(t) = 1 − e− m .
k
De acuerdo a la condición dada que la resistencia al movimiento es de 160 dinas cuando la velocidad
es de 40 centı́metros por segundo y dicha fuerza es proporcional a la velocidad instantánea, entonces
podemos hallar el valor de k, ya que 160 = 40k, luego k = 4. Reemplazando este valor en la solución
final tenemos que v(t) = 49[1 − e−20t ]. Para hallar la distancia recorrida en función del tiempo
podemos emplear la condición que v(t) = dsdt luego:

ds
= 49 − 49e−20t .
dt
Esta es una ecuación de variables separables, entonces
 s  t  t
ds = 49dt − 49e−20s ds
0 0 0
49
s(t) = 49t + [e−20t − 1].
20
La velocidad lı́mite se tiene cuando t se hace tan grande que podemos decir que t → ∞, entonces en
este caso e−20t tiende a cero, luego la velocidad lı́mite es 49 centı́metros por segundo.

Ejemplo 3.4.2. Un dı́a comenzó a nevar a una

tasa pesada y constante. Una máquina quitanieves
comenzó a trabajar a medio dı́a, avanzando a 2 millas
la primera hora y a una milla la segunda hora. ¿A qué
hora comenzó a nevar?
Figura 3.13: Dı́a nevado.

Solución: En primer lugar, los datos del problema están de acuerdo con la idea de que el quitanieves se
moverá más lentamente a medida que la capa de nieve aumenta. Sin pretender determinar cuáles
asunciones pueden ser buenas o malas, asumimos que el quitanieves quita la nieve a una tasa
constante de K metros cúbicos por hora. Luego t será el tiempo medido en horas a partir del medio
dı́a, y x denotará la profundidad en mm de nieve en un tiempo t. Además y denotará la distancia que
d
se mueve el quitanieves con respecto al tiempo. Luego dt es la velocidad de quitanieves y nuestra
dy
asunción será wx dt = K, donde w es la anchura que da el quitanieves. Para encontrar cómo x
 depende del tiempo, t0 será el número de horas antes del medio dı́a cuando comenzó a nevar, y S
será la tasa constante en mm en la cual la profundidad de la nieve se incrementa. Luego usando
t > −to
x = S(t + to ),
y obtenemos la ecuación diferencial:
dy k
= .
dt ws(t + to )
k
Esta ecuación diferencial tiene la forma dy
dt = f (t), y vemos que y = ws [ln(t + t0 ) + C] donde C es
una constante. Se puede sospechar o especular que el conocimiento de muchos pares de valores de y
y t deben permitirnos determinar t0 y ası́ obtener una solución para el problema. En y = 0 cuando
t = 0 encontramos que en la ecuación anterior C = − ln t0 . Entonces:
 
k t
y= ln 1 + .
ws t0

Luego usamos el hecho de que y = 2 cuando t = 1, y y = 3 cuando t = 2 para obtener:

 2  3
2 1
1+ = 1+ .
to to

Expandiendo los exponentes y simplificando la ecuación t2o + to − 1 = 0, donde to > 0, obtenemos

√
(−1 + 5)
to = = 0, 618 horas = 37 min.
2
Por lo tanto, la hora en que empezó a nevar fue a las 11:23 a.m.

Ejemplo 3.4.3. Una partı́cula se mueve en lı́nea recta hacia un punto fijo 0, bajo la acción de una
fuerza atractiva en 0 que varı́a directamente con su distancia a 0. Cuando t = 0 la partı́cula dista 4
cms de 0 y se mueve hacia 0 con una velocidad de 6 centı́metros por segundo, y una aceleración de 16
centı́metros por segundo cuadrado.

a) Hállense su posición y velocidad de partı́cula en función del tiempo.

b) Obténgase la amplitud, el periodo y la frecuencia del movimiento.

c) Determı́nese la velocidad y aceleración máxima.

Solución: Tomemos un punto fijo 0 como origen de un sistema de coordenadas, supongamos que la
partı́cula inicia su movimiento en un punto A, tomemos 0A como sentido positivo. Sea P la posición
de la partı́cula en un instante t. Como la magnitud de la fuerza de atracción F hacia 0 es proporcional
a su distancia a este punto tendremos F = kx, con k > 0. Supongamos x ≥ 0, luego la magnitud
de la fuerza es kx. Como la fuerza F está dirigida hacia 0 (en sentido negativo) tendremos que
F = −kx con x mayor o igual a cero. Si x < 0, luego la magnitud de la fuerza es −kx, como
la fuerza F está dirigida hacia la derecha (en sentido positivo) luego F = −kx con x menor que
2
cero; luego la ecuación diferencial que describe este fenómeno fı́sico es m ddt2x = −kx. Para nuestro
caso la velocidad y la aceleración está dirigida hacia O, cuando la partı́cula dista 4 cm o sea para
t = 0 tendremos x(0) = 4, v(0) = −6, a(0) = −16, cada una con sus respectivas unidades del
k k k
sistema CGS, ahora podemos calcular el cociente m ya que −16 = − m 4, luego m = 4. La ecuación
d2 x
diferencial original se transforma en dt2 = −4x. Para resolver esta ecuación de segundo orden la
 dx
transformamos en una de primer orden, mediante el siguiente cambio de variable, como v(t) = dt ,
de modo que
d2 x dv dv dx dv
= = · = v· .
dt dt dx dt dt
Luego la ecuación diferencial de segundo orden se convierte en:
dv
v· = −4x.
dx
La cual es una ecuación diferencial de variables separables, su solución está dada por
 v  x
v dv = − 4x dx.
−6 4

Entonces:
v 2 (t) 36
− = −2x2 + 32
2 2
v 2 (t) = 4(25 − x2 )

v(t) = ±2 25 − x2 .
dx
Como v(t) = dt , reemplazando en la última ecuación tenemos una ecuación de variables separables.
dx 
= 2 25 − x2 (3.47)
dt
dx 
= −2 25 − x2 . (3.48)
dt
La solución de la ecuación diferencial (3.47) es
 x  t
dx
√ = 2dt (3.49)
4 25 − x 0
x 4
arcsin − arcsin = 2t. (3.50)
5 5
Definimos u = arcsin 52 , luego sin u = 4
5 , y, cos u = 35 , luego la ecuación (3.50) se transforma en
x
arcsin = 2t + u.
5
Aplicando la función seno tenemos:

x(t) = 5 sin(u + 2t)

= 5 cos 2t sin u + sin 2t cos u.

Luego

x = 4 cos 2t + 3 sin 2t
v = −8 sin 2t + 6 cos 2t.

Podemos definir el signo, tanto para x como para v, usando la condición dada de que v(0) = −6, lo
cual implica que:

x(t) = 4 cos 2t − 3 sin 2t

v(t) = −8 sin 2t − 6 cos 2t
a(t) = −16 cos 2t + 12 sin 2t.
 Para calcular la amplitud, periodo y frecuencia del movimiento recurrimos a la fórmula x(t) =
5 sin(u − 2t). Luego v(t) = −10 cos(u − 2t) y a(t) = −20 sin(u − 2t). Es fácil ver que la amplitud
es 5, el periodo es π y la frecuencia es π1 , cada uno con sus respectivas unidades. La velocidad
máxima es de 10 centı́metros por segundo y la aceleración máxima es de 20 centı́metros por segundo
cuadrado.
Ejemplo 3.4.4. Para movimientos de gran rapidez en el aire, como el que realiza un paracaidista, que
está cayendo antes de que se abra el paracaı́das la resistencia del aire es cercano a una potencia de la
velocidad instantánea v(t). Determinar una ecuación diferencial para la velocidad v(t) de un cuerpo de
masa m que cae, si la resistencia del aire es proporcional al cuadrado de la velocidad instantánea (ver
Figura 3.14).

Figura 3.14: El paracaidista.

Solución: En la dinámica del movimiento actúan dos fuerzas f1 el peso y f2 la resistencia del aire. A
partir de las leyes de la fı́sica se establece que la fuerza total es

f = f1 + f2 .

El peso actúa sobre el cuerpo hacia abajo (dirección positiva). La resistencia (kv 2 ) con k > 0 actúa
para impedir el movimiento, va hacia arriba (dirección negativa).

f = mg − kv 2 .
Utilizando la segunda ley de Newton 
f = ma.
Igualando las sumatorias de fuerzas
mg − kv 2 = ma,
donde
dv
a= .
dt
Sustituyendo en la ecuación anterior se obtiene
dv
mg − mg − kv 2 = m .
dt
Ejemplo 3.4.5. Un paracaidista cuya masa es de 75kg se deja caer desde un helicóptero, que está a 4000
m de altura. La fuerza debida a la resistencia del aire es proporcional a la velocidad instantánea de caı́da,
siendo k = 15kg/s la constante de cuando el paracaı́das está cerrado y k = 105 kg/s cuando está abierto.
Si el paracaı́das se abre un minuto después de saltar. Calcular la velocidad v(t) y la correspondiente
posición s(t) en cada instante. ¿Cuánto tarda en llegar a la superficie?
Solución: Podemos observar que tenemos dos situaciones diferentes una (P1 ) antes de abrir el paracaı́das
y otra (P2 ) después de abrir el paracaı́das, al cabo de 60 s. Esto da lugar a dos problemas de valor
inicial diferentes:
dv
m = mg − 15v , v(0) = 0 (3.51)
dt
dv
m = mg − 105v , v(60) = v1 . (3.52)
dt
El PVI (3.51) corresponde a P1 y el PVI (3.52) corresponde a P2 . La solución del problema P1 , que
es válida para 0 ≤ t ≤ 60, viene dada por:
mg mg −15t −t
v(t) = − e m = 5g(1 − e 5 ).
15 15
Haciendo uso de esta expresión, obtenemos que v(60) ≈ 49m/s = v1 . Podemos ahora resolver el
problema P2 , obteniendo que
5g 7t
v(t) = (1 + 6e84− 5 ),
7

para t > 60. Si calculamos el desplazamiento s(t) = v(t)dt, tenemos que, para 0 ≤ t ≤ 60, y
teniendo en cuenta que s(0) = 0 (se mide a partir de la posición del helicóptero), se obtiene que
−t
s(t) = 5g(t + 5e 5 − 5),

para 0 ≤ t ≤ 60. Cuando t = 60, al abrir el paracaı́das, se tiene que el espacio recorrido es
5 84− 7t
s(t) = 5g(55 + 5e−12 ) ≈ 2695m. Para t > 60 la solución es s(t) = 5g
7 (t − 7 6e
5 ) + 2305. Por

último, cuando llega al suelo debe ser s(t) = 4000. Por tanto, resolvemos la ecuación en t,
 
5g 5 84− 7t
t − 6e 5 + 2305 = 4000
7 7
5 7t
t − 6e84− 5 = 242,
7
7t
y cuando t > 60 se verifica que t − 75 6e84− 5 ≈ t, luego el paracaidista tarda en llegar a la superficie
t ≈ 242 s

Ejemplo 3.4.6. Un cuerpo de 8 lb de peso cae partiendo del reposo desde una gran altura. Conforme cae,
actúa sobre él la resistencia del aire a la que supondremos (en libras) numéricamente igual a 2v, siendo
v la velocidad (en pies por segundo). Hallar la velocidad y la distancia recorrida al cabo de t segundos.

Solución Elegimos el eje x positivo como vertical y hacia abajo, a lo largo de la trayectoria del cuerpo
B, el origen en el punto en que el cuerpo inicia su caı́da. Las fuerzas que actúan sobre el cuerpo
son:

1. F1 , su peso, 8lb, actuando hacia abajo y, por tanto positiva.

2. F2 , la resistencia del aire numéricamente igual a 2v, actuando hacia arriba y por lo tanto
negativa y con valor de −2v.

La segunda ley de Newton F = ma se transforma en m = dv dt = F 1 + F 2. Tomando g = 32 y

utilizando m = wg = 8, 32 = 14 , y puesto que el cuerpo se encontraba inicialmente en reposo se
obtiene
1 dv
= 8 − 2v (3.53)
4 dt
v(0) = 0. (3.54)
 Figura 3.15: Geometrı́a.

La ecuación (3.53) es separable. Separando variable tenemos

dv
= 4dt.
8 − 2v
Integrando hallamos que
1
− ln 8 − 2v = 4t + c0 .
2
Que se reduce a
8 − 2v = c1 e−8t .
Aplicando la condición de (3.54) encontramos que c1 = 8 por lo que la velocidad en el instante t
viene dada por
v = 4(1 − e−8t ). (3.55)
Para determinar ahora la distancia recorrida hasta el instante t escribimos la ecuación (3.55) en la
forma
dx
= v = 4(1 − e−8t )
dt
Observe que x(0) = 0. Integrando la ecuación anterior tenemos
 
1
x = 4 t + e−8t + c2
8

Puesto que x = 0 cuando t = 0 encontramos que c2 = 21 por lo que en definitiva, la distancia

recorrida viene dada por  
1 −8t 1
x=4 t+ e − . (3.56)
8 2
Interpretación de Resultados: La ecuación (3.54) demuestra que cuando t → ∞, la velocidad v se
aproxima a la velocidad lı́mite 4(pies/s). Observamos también que esta velocidad lı́mite se alcanza
aproximadamente en un tiempo muy corto. La ecuación (3.56) demuestra que cuando t → ∞,
x → ∞. ¿Implica esto que el cuerpo penetrará en la tierra y continuará eternamente? Por supuesto
que no, ya que cuando el cuerpo alcance la superficie de la tierra detendrá su movimiento. ¿Cómo
reconciliar este final evidente del movimiento con lo que enuncia la ecuación (3.56)? Es bien sencillo,
cuando el cuerpo alcanza la superficie terrestre la ecuación diferencial (3.53) y, por lo tanto, la
ecuación (3.56), dejan de ser aplicables.
 Ejemplo 3.4.9. Suponga que en principio un gran depósito de mezclado contiene 300 galones de agua
en la que se han disuelto 50 libras de sal. Otra solución de salmuera se bombea hacia el depósito a razón
de 3 galones por minuto, y cuando la solución está bien agitada, se bombea hacia afuera sólo 2 galones
por minuto. Si la concentración de la solución entrante es 2 libras por galón, determine una ecuación
diferencial para la cantidad de sal A(t) que se encuentra en el tanque en el instante t.

Figura 3.18: Gráfica del problema.

Solución: Sabemos que la cantidad de sal (A) que se encuentra en el tanque en un tiempo t, está dado
por la cantidad de sal que entra al recipiente y la cantidad de sal que sale de este; matemáticamente
esto es expresado como
dA
= Rentrada − Rsalida ,
dt
donde la cantidad de sal que entra (Rentrada ) está dada por

gal lib lib

Rentrada = 3 ∗2 =6 ,
min gal min

donde el primer término de esta expresión es la velocidad de entrada de la solución, y el segundo es

la concentración de este flujo. Y la cantidad de sal que sale es
 
A lib gal 2A lib
Rsalida = (2 )= .
300 + t gal min 300 + t min

En esta expresión el primer término es la concentración de sal en el flujo de salida, la concentración

de esta, está dada por la cantidad de soluto (sal) sobre la cantidad de solvente (agua), y está última
varı́a con el tiempo porque el flujo de salida es menor que el de entrada. Retomando tenemos que

dA
= Rentrada − Rsalida
dt
lb 2A lb
= 6 − .
min 300 + t min

Para resolver esta ecuación diferencial, entonces tenemos que:

dA 2A
=6− .
dt 300 + t
Llevando la ecuación a la forma estándar se obtiene

dA 2A
+ = 6,
dt 300 + t
entonces,
2
P (x) = , f (x) = 6.
300 + t
Factor Integrante   2
P (x)dx dt
e =e 200+t = e2ln|300+t| = (300 + t)2 .
Extrapolando con la solución estándar encontramos que

c 1
A(t) = + 6(300 + t)2 dt
(300 + t)2 (300 + t)2
c 2(300 + t)3
= + .
(300 + t)2 (300 + t)2

Simplificando
c
A(t) = + 2(300 + t).
(300 + t)2
3.4.6. Aplicaciones a volúmenes

Supongamos que el recipiente es el de la Figura 3.25, r

H
donde A(y) es el área de la sección recta del depósito
a una altura y y B es el área del orificio con bordes y
perfectamente pulimentados.

Figura 3.25: Volumen 2.

Si no existiera rozamiento no hay pérdida de energı́a potencial, toda se convertirı́a en energı́a cinética
la cual serı́a la única que aparecerı́a en este fenómeno fı́sico. Luego
1
mgy = mv 2 ,
2
1
donde y es la altura de caı́da, luego la velocidad de salida serı́a de (2gy) 2 metros por segundo. Como el
1
área de orificio es B entonces el producto B(2gy) 2 representa el número de metros cúbicos por segundo
de lı́quido que sale por el orificio. Sea c el coeficiente de rozamiento, es decir, el chorro no es completo por
1
el orifico, luego la velocidad de descarga es cB(2gy) 2 metros cúbicos por segundo. Sea V (y) el volumen
del lı́quido que está en el depósito a una altura y entonces:
 y
dv
V (y) = A(u)du, y, = A(y).
0 dy
Luego
dv 1
= −cB(2gy) 2 .
dt
Este resultado es negativo porque el volumen decrece con el tiempo. Usando la regla de la cadena tenemos:
dv dv dy dv
= . = A(y) ,
dt dy dt dt
entonces combinando la dos últimas ecuaciones tenemos la ecuación que describe nuestro suceso
dy 1
A(y) = −cB(2gy) 2 .
dt
La cual es una ecuación de variables separables de primer orden.
Ejemplo 3.4.18. Un embudo tiene la forma de un cono circular recto, cuyo vértice está hacia abajo y
está lleno de agua. Si la mitad del volumen se desocupa en un tiempo T , encontrar el tiempo requerido
para vaciarse completamente.
Solución: El gráfico del recipiente y sus elementos son: H la altura del cono, R el radio mayor, r el
radio de la sección recta A (y) a una altura y del vértice, como la altura depende del tiempo, sea
y = f (t) la función que determina la altura del lı́quido en el recipiente, de acuerdo a la condición
dada que la mitad del volumen sale en un tiempo T , podemos hallar la altura del lı́quido que queda
2 2
en función de H. Tenemos que Vtotal = πR3 H , y, 21 Vtotal = πR6 H . De acuerdo a la semejanza de
2h h3
triángulos, tenemos que: Rr = H h
. Luego, Rr2 H = H 3 como la razón entre los dos volúmenes es de
1:2, entonces
1 h3
= 3.
2 H
Luego
H 3 = 2h3 ,
de donde
H
f (T ) = h = √ .
2
Como
πR2 y 2
A(y) = πr2 = .
H2
dy 1
Luego la ecuación diferencial A(y) = dt = −CB(2gy) 2 se transforma en

πR2 y 2 dy 1

2
= −cB(2gy) 2 ,
H dt
cuya solución viene dada por
 y  t
πR2 3
√ y 2 dy = − cBdt.
H 2 2g H 0

Luego,
2πR2  5 5

√ y 2 − H2 = −cBt.
5H 2 2g
Para y = 0 tenemos el tiempo total, entonces
5
2πR2 H 2
√ = cBt total.
5H 2 2g

Luego, √
2πR2 H
ttotal = √ .
5cB 2g
El cual aparece en función de H, R, c y B, y lo importante es dejarlo en función únicamente de T ,
H
3 , luego podemos reemplazar esta condición en la solución
pero la condición dada de que f (T ) = √
2
a la ecuación diferencial y tenemos
1  
2πR2 H 2 1
√ √
6
− 1 = −cBT
5cB 2g 32
√ √
2πR2 H 6
32
√ =T√
6
.
5cB 2g 32 − 1
Luego,
ttotal = 2, 28T.
 h
Oricio
Circular

Figura 3.26: Tanque con agujero circular.

Ejemplo 3.4.19. Suponga que está saliendo agua de un tanque a través de un agujero circular de
área A0 que está en el fondo. Cuando el agua sale a través del agujero, la fricción y la contracción√ de
la corriente cerca del agujero reducen el volumen de agua que sale del tanque por segundo a cA0 2gh
donde c donde 0 < c < 1 es una constante empı́rica. Determine una ecuación diferencial para la altura h
del agua al tiempo t para el tanque cúbico, que se muestra en la figura. El radio del agujero es de 2 pulg
g = 32 pies
s2
, el área del águjero es Ah = π(2pulg)2 , el área de la base es Aw = 100, c es la constante que
√
depende de las propiedades fı́sicas del lı́quido. Partimos de: dv dt = −Ah 2gh . En primer lugar, podemos
obtener el volumen del tanque v = 100h. Como la altura y el volumen varı́an con respecto al tiempo,
tenemos
dv dh dh 1 dv
= 100 ⇔ = .
dt dt dt 100 dt
Pero
dv  dv √
= cAh 2gh ⇔ = cπ(2)2 2 ∗ 32 ∗ h.
dt dt
 2
dv 1 √
= cπ 64h
dt 6
 
dv 1 √
= cπ 8 h
dt 36
dv 2 √
= cπ h.
dt 9
dv dh
Sustituyendo dt en dt se obtiene
dh 1 2 √
= cπ h
dt 100 9
√
dh cπ h
= .
dt 450
Pero como el nivel del agua está disminuyendo concluimos
√
dh cπ h
=− .
dt 450
Ejemplo 3.8. Un termómetro se saca de una habitación, donde la temperatura
del aire es de 70°F, al exterior, donde la temperatura es de 10°F. Después de 1/2
minuto el termómetro marca 50°F. ¿Cuánto marca el termómetro cuando t = 1
minuto? ¿Cuánto demorará el termómetro en alcanzar los 15°F?

Solución. Es claro que Tm = 10°F, de modo que reemplazamos en la ecua-

ción (3.10):

dT
= k(T − 10) (3.11)
dt
T (0) = 70°F
1
 
T = 50°F.
2

Se debe resolver la ecuación (3.11) con condiciones iniciales, la cual es lineal

y tambien de variables separables. Empleando la técnica de resolución de
ecuaciones de variables separables obtenemos:

dT
=kdt
T − 10
w dT w
= kdt
T − 10
ln | T − 10 |=kt + c1
T − 10 =c2 ekt , c2 ∈ ℝ.

Entonces T = 10 + c2 e kt .
Cuando t = 0, T = 70°F, de modo que 70 = 10 + c2 , es decir, c2 = 60 y,
por lo tanto, T = 10 + 60e kt . Ahora podemos hallar k ya que T (1/2) = 50°F
1
50 =10 + 60e k 2
1 40
e2k =
60  
2
k =2 ln ≈ −0,81093.
3

En consecuencia,
T (t) = 10 + 60e −0,81093t .
100 · Algunas aplicaciones de edo de primer orden

Para t = 1 tenemos T (1) = 10 + 60e −0,81093·1 = 36,66 ≈ 37°F. El termó-

metro marca 37°F cuando t = 1 minuto. Ahora debemos hallar cuánto se
demora para alcanzar 15°F,

15 =10 + 60e −0,81093t

5
=e −0,81093t
60
ln 1/12


t= = 3,06427 ≈ 3,1 min.
−0,81093

Para alcanzar 15°F se demora 3,1 minutos. X



Ahora queremos formular un módelo matemático que describa el perfil

de temperatura dentro de un edificio durante 24 horas como función de la
temperatura exterior, el calor generado dentro del edificio y el calefactor o
aire acondicionado.
Sea T (t) la temperatura dentro del edificio en el instante t y veamos al
edificio como un único comportamiento. La razón de cambio en la tem-
peratura queda determinada por todos los factores que generan o disipan
calor.
Factores que afectan la temperatura dentro del edificio:
· Calor generado por las personas, las luces y las máquinas dentro del
edificio. Esto incrementa la temperatura y se denota por H (t), que
siempre es positiva.
· Calentamiento o enfriamiento proporcionado por calefacción o aire
acondicionado. Esta razón de incremento (o decremento) de tempera-
tura se denota por U (t), que es positiva para la calefacción y negativa
para el aire acondicionado. H (t) y U (t) están en términos de tempe-
ratura por unidad de tiempo.
· Efecto de la temperatura exterior M (t) sobre la temperatura dentro
del edificio.
El tercer factor lo podemos modelar mediante la ley de enfriamiento de
Newton. Esto es: la razón de cambio de la temperatura T (t) es proporcional
a la diferencia entre la temperatura exterior M (t) y la temperatura interior
T (t). Es decir,
dT
= K [M (t) − T (t)].
dt
La constante depende de las propiedades físicas del edificio, como cantidad de
puertas, ventanas y el tipo de aislamiento. Por lo tanto, el modelo matemático
que describe el perfil de temperatura dentro de un edificio considerando estos
factores que generan o disipan calor viene dado por la ecuación diferencial
dT
= K [M (t) − T (t)] + H (t) + U (t), (3.12)
dt
donde
dT
+ KT (t) = KM (t) + H (t) + U (t)
dt
es una ecuación diferencial lineal con

p(t) = K
Q (t) = KM (t) + H (t) + U (t).

El factor integrante es r
Kdt
𝜇(t) = e = e Kt .
Utilizando la fórmula (2.14) obtenemos:
hw i
T (t) = e −Kt e Kt [KM (t) + H (t) + U (t)]dt + c

y K1 es la constante de tiempo para el edificio (sin calefacción o aire acondi-

cionado).

Ejemplo 3.9. En una calurosa mañana de sábado, cuando las personas trabajan
dentro del edificio, el aire acondicionado mantiene la temperatura interior en 24°C. A
mediodía, el aire acondicionado se apaga y las personas se van a casa. La temperatura
exterior es constante e igual a 35°C durante el resto de la tarde. Si la constante de
tiempo del edificio es de 4 horas, ¿cuál será la temperatura dentro del edificio a las 2
pm? ¿Y a las 6 pm? ¿En qué momento llegará la temperatura interior del edificio a
27°C?

Solución. Se tiene M (t) = 35°C, T (0) = 240 °C, H (t) = 0, U (t) = 0,

1 1
K = 4, por lo tanto tenemos que K = 4 . Hallemos T (2). Reemplazamos en
la ecuación:
hw i
T (t) = e −Kt e Kt [KM (t) + H (t) + U (t)]dt + c

obtenemos,
w
1 1

1
T (t) = e − 4 t e 4 t 35dt + c
4
1
h 1 i
T (t) = e − 4 t e 4 t 35 + c
1
T (t) = 35 + ce− 4 t .
Como T (0) = 24°C, concluimos que
1
24 = 35 + ce− 4 0
c = −11.
1
Por lo tanto, T (t) = 35 − 11e − 4 t , donde reemplazando podemos hallar
1
T (2) = 35 − 11e − 4 2 = 28,32°C.
La temperatura a las 2 pm es de 28,32°C. La temperatura a las 6 pm es de
32,54°C. Para hallar en qué momento t0 la temperatura llega a 27°C reem-
1 1
plazamos en T (t0 ) = 35 − 11e − 4 t0 ; es decir, 27 = 35 − 11e − 4 t0 ; despejando
t0 obtenemos t0 = 1:16 p.m. X

Ejemplo 3.10. Se va a construir un almacén sin calefacción ni aire acondicionado.

Según la cantidad de aislamiento, la constante de tiempo para este edificio puede
variar de 1 a 5 horas. Para ilustrar el efecto del aislamiento sobre la temperatura
dentro del almacén, suponga que la temperatura exterior varía como una onda
sinusoidal, con un mínimo de 16°C a las 2 a.m. y un máximo de 32°C a las 2 p.m.
Suponiendo que el término exponencial (que implica la temperatura inicial T0 ) se
ha extinguido, ¿cuál es la temperatura mínima dentro del edificio si la constante
de tiempo es 1 hora? ¿Y si la constante de tiempo es 5 horas? ¿Cuál es la máxima
temperatura dentro del edificio si la constante de tiempo es 1 hora? ¿Y si es 5 horas?

Solución. Como M (t) varía como una onda sinusoidal en un periodo de

24 horas, con un mínimo 16°C a las 2 a.m. y un máximo de 32°C a las 2
p.m. t = 12. Es decir.
M (t) = M0 − B cos(wt),
donde B es una constante positiva, M0 es la temperatura exterior promedio:
32 + 16 2𝜋 𝜋 radianes
M0 = = 24°C y w = = .
2 24 12 hora
Además U (t) = 0 porque no hay calefacción ni aire acondicionado. Ade-
más, H (t) = H0 = 0, pues no hay razón de calentamiento dentro del edificio.
De la ecuación (3.12) tenemos:
dT
= K M0 − B cos(wt) − T (t) + H0
 
dt
dT
+ KT = K M0 − B cos(wt) + H0 . (3.13)
 
dt
Llamemos a la temperatura promedio interior B0 , en este caso B0 = M0 + HK0 .
Reemplazamos en (3.13):
dT
+ KT = K [B0 − B cos(wt)] ,
dt
donde KB0 representa el valor promedio diario de KB0 − KB cos(wt); es
decir,
24
1 w
KB0 = (KB0 − KB cos(wt)) dt.
24
0

La solución de la ecuación es
hw i
T (t) = e −Kt e Kt [KB0 − KB cos(wt)] dt + c
  
B K w sin(wt) + K
w cos(wt)
T (t) = B0 −  2 + ce−Kt .
K
1+ w

Para t = 0 y T0 = 16°C, reemplazamos en la última ecuación:

 2
K
B w
16 = B0 −   2 + c.
K
1+ w

Por lo tanto,
 2
K
B w
c = 16 − B0 +  2 .
K
1+ w

Es decir,
 2
K
B w
cos(wt) + K sin(wt)

B w
T (t) = B0 − + ­­16 − B0 +   2 ®® e .
© ª −Kt
2
1 + wt


« 1+ Kw ¬

Como H0 = 0 y B0 = M0 , cuando no hay razón de calentamiento adicional

dentro del edificio H0 = 0 , la temperatura promedio interior B0 es igual

a la temperatura promedio exterior M0 . Por lo tanto, B0 = 24°C. Como el

término exponencial se ha extinguido, entonces c = 0. Luego, reemplazamos
  
B Kw sin(wt) + K
w cos(wt)
T (t) = 24 −  2 .
K
1+ w

Para hallar el valor de B tomamos

M (t) = M0 + B cos(wt),
ahora reemplazamos los valores de M0 y w:
 𝜋 
M (t) = 24 − B cos t .
12
Para t = 0 y M (t) = 16 se tiene B = 8. Por lo tanto,
  
K K
8 w sin(wt) + w cos(wt)
T (t) = 24 −  2 .
1+ K w

Para hallar la temperatura mínima para t = 0, K1 = 1, se tiene K = 1 y

tomamos w = 2𝜋 1
24 que es aproximadamente 4 . Reemplazando en la última
ecuación obtenemos:
8(4) 2
T (0) = 24 − ≈ 16,5°C;
1 + (4) 2
1
para K = 5, se tiene K = 15 . La temperatura mínima es
 2
4
8 5
T (0) = 24 −   2 =≈ 20,8°C.
4
1+ 5

1
Para hallar la temperatura máxima t = 12, K = 1, se tiene K = 1

8(4) 2
T (12) = 24 + ≈ 31,5°C
1 + (4) 2
1
para K = 5, se tiene K = 15 ; la temperatura máxima es
 2
4
8 5
T (12) = 24 +   2 =≈ 27,1°C.
4
1+ 5

Obtenemos así las respuestas a todas las preguntas. X



Ejemplo 3.11. Durante el verano la temperatura dentro de una camioneta llega a

55°C, mientras que en el exterior es constante e igual a 35°C. Cuando la conductora
entra a la camioneta, enciende el aire acondicionado con el termostato en 16°C. Si
la constante de tiempo para la camioneta es K1 = 2 horas y para la camioneta con el
aire acondicionado es K11 = 13 hora, ¿en qué momento llegará la temperatura dentro
de la camioneta a los 27°C?
Solución. K11 es la constante de tiempo con calefacción y aire acondicio-
nado, donde K1 = K + Ku . En la camioneta se instala un termostato que
se utiliza para comparar la temperatura real dentro de la camioneta con
una temperatura deseada TD . Si la temperatura real es menor que la tempe-
ratura deseada, el calefactor comienza a funcionar y, en caso contrario, se
desconecta. Si la temperatura real es mayor que la temperatura deseada, el
aire acondicionado comienza a enfriar y en caso contrario se desconecta. La
cantidad de calentamiento es proporcional a la diferencia de temperatura

U (t) = Ku [TD − T (t)] ,

donde Ku es la constante de proporcionalidad (positiva). Así, M (t) = 35°C,

T (t) = 55°C, H (t) = 0, TD = 16°C, K1 = 2. Entonces K = 12 , K11 = 13 , luego
K1 = 3. Como U (t) = Ku [TD − T (t)], de la ecuación (3.12) se tiene:
dT
= K [M (t) − T (t)] + Ku [TD − T (t)]
dt
dT
+ T [K + Ku ] = KM (t) + KuTD . (3.14)
dt
Vemos que para este ejemplo, la cantidad p(t) = K es igual a K + Ku ; por lo
tanto,

K1 = K + Ku
1
3 = + Ku
2
5
Ku = .
2
Reemplazando en (3.14) obtenemos:
dT 1 5
+ 3T = 35 + 16
dt 2 2
dT 115
+ 3T = ,
dt 2
la cual es una ecuación diferencial lineal. Su factor integrante es 𝜇(t) = e 3t ,
por lo tanto la solución es
w
3t 115 3t 115 115
  
T =e −3t
e dt + c = e −3t
e dt + c = + ce−3t .
2 6 6
Como T (0) = 55, tenemos:
115
55 = +c
6
luego,
215
c=
6
y
115 215 −3t
T (t) = + e .
6 6
Para hallar en qué momento t0 la temperatura llegará a 27°C, reemplaza-
mos en la última ecuación y despejamos t0 .

115 215 −3t0

27 = + e ,
6 6
Obtenemos así que t0 = 0,5068 horas, o equivalentemente t0 = 30,4
minutos. X


Ejercicios
1. Un termómetro se lleva del interior de una habitación al exterior,
donde la temperatura del aire es 60°F. Después de un minuto, el
termómetro indica 58°F; cinco minutos después marca 32°F. ¿Cuál
era la temperatura interior?
2. Una barra metálica pequeña, cuya temperatura inicial es de 25°F, se
deja caer en un recipiente con agua hirviente cuya temperatura es de
212°F. ¿Cuánto tiempo tardará en alcanzar 100°F si se sabe que la
temperatura aumentó 3°F en un segundo? ¿Cuánto tiempo tardará en
llegar a 105°F?
3. Una taza de café caliente, inicialmente a 100°C, se enfría hasta 85°C
en 5 minutos al estar en un cuarto con temperatura de 21°C. Usar solo
la ley de enfriamiento de Newton y determinar el momento en que la
temperatura del café estará a unos 55°C.
4. Una cerveza fría, inicialmente a 33°F, se calienta hasta 45°F en 3
minutos, en un cuarto con temperatura 65°F. ¿Qué tan caliente estará
la cerveza si se deja ahí durante 10 minutos?
5. Una taza de té caliente está inicialmente a 190°F y se deja en un cuarto
que tiene una temperatura ambiente de 75°F. Suponga que a partir
del tiempo t = 0 se enfría a razón de 25°F por minuto.
a) Suponga que se aplica la ley de Newton sobre el enfriamiento.
Escribir un problema de valor inicial que modele la temperatura
del té caliente.
b) ¿Cuánto tiempo le toma al té caliente enfriarse a una temperatura
de 120°F?
 4. Mecánica de Newton
La mecánica es el estudio del movimiento de los objetos y el efecto de las
fuerzas que actúan sobre ellos. La mecánica de Newton, o clásica, trata del
movimiento de los objetos comunes; es decir, de los objetos que son grandes
en comparación con un átomo y lentos en comparación con la velocidad de
la luz.
Un modelo de la mecánica de Newton se puede basar en las leyes del
movimiento de Newton:
· En ausencia de fuerzas, un objeto (cuerpo) en reposo seguirá en reposo,
y un cuerpo moviéndose a una velocidad constante en línea recta, lo
continuará haciendo indefinidamente.
· Cuando un cuerpo es sujeto a una o más fuerzas externas, la razón de
cambio temporal del momento del cuerpo es igual a la suma vectorial
de las fuerzas externas que actúan sobre él.
· Cuando un cuerpo interactúa con otro, la fuerza del primer cuerpo
sobre el segundo es igual en magnitud, pero opuesta en dirección, a la
fuerza del segundo cuerpo sobre el primero.
La segunda ley de Newton, que solo se aplica a marcos de referencia
inerciales, nos permite formular las ecuaciones de movimiento para un
cuerpo. Podemos expresar la segunda ley de Newton como
dp
= F (t, x, v),
dt
donde F (t, x, v) es la fuerza resultante sobre el cuerpo en el instante t,
posición x, velocidad v y p(t) es el momentum del cuerpo en el instante t. El
momentum es el producto de la masa del cuerpo y su velocidad, es decir,
p(t) = mv(t).
De este modo, podemos expresar la segunda ley de Newton como
dv
m = ma = F (t, x, v),
dt
donde a = dv
dt es la aceleración del cuerpo en el instante t.
En esta sección nos centraremos en situaciones donde la fuerza F no
depende de x. Luego la ecuación de primer orden en v(t) es
dv
m = ma = F (t, v). (3.15)
dt
Ejemplo 3.12. Un objeto de masa 5 kg recibe una velocidad inicial hacia abajo de
50 m/s y luego se le permite caer bajo la influencia de la gravedad. Suponga que la
fuerza en newtons debida a la resistencia del aire es −10v, donde v es la velocidad
del objeto en m/s. Determinar la ecuación de movimiento del objeto. Si el objeto está
inicialmente a 500 m sobre el suelo, determinar el momento en que el objeto golpeará
el suelo (figura 3.3).

b
F2 = −bv(t)
t=0
x(t)
b

F1 = mg
b t

Figura 3.3: Caída libre.

Solución. Hay dos fuerzas actuando sobre el objeto: una fuerza constante
debida al empuje hacia abajo de la gravedad y una fuerza debida a la resis-
tencia del aire que es proporcional a la velocidad del objeto, que actuán en
forma opuesta al movimiento del objeto; por lo tanto, el movimiento del
objeto se realizará a lo largo de un eje vertical. En este eje, elegimos el origen
como el punto donde el objeto fue lanzado inicialmente y definimos x(t)
como la distancia que ha caido el objeto hasta el instante t.
Las fuerzas que actúan sobre el objeto a lo largo de este eje se pueden
expresar de la siguiente manera:
· La fuerza debida a la gravedad

F1 = mg ,

donde g es la aceleración de la gravedad en la superficie de la tierra.

· La fuerza producto de la resistencia del aire:

F2 = −bv(t).

Aquí b > 0 es la constante de proporcionalidad y el signo negativo

está presente dado que la resistencia del aire actúa en forma opuesta al
movimiento del objeto.
· La fuerza neta que actúa en el objeto es:

F = F1 + F2 = mg − bv(t). (3.16)
Ahora, utilizando la ecuación (3.15), y sustituyendo en (3.16) se obtiene:
dv
m = mg − bv, v(0) = v0 , (3.17)
dt
donde m = 5 kg, v0 = v(0) = 50 ms , b = 10, g = 9,81 m
s2
. Reemplazamos en
(3.17):
dv m
5kg = 5kg · 9,81 2 − 10v(t)N
dt s
dv 5kg m 10N
= · 9,81 2 − v(t)
dt 5kg s 5kg
dv m m
= 9,81 2 − 2 2 v(t)
dt s s
dv
+ 2v(t) = 9,81.
dt
Utilizando el método de resolución para ecuaciones de tipo variables separa-
bles o lineal, hallamos la solución de esta ecuación diferencial:

v(t) = 4,905 + ce−2t ,

como v(0) = 50 se obtiene, c = 45,095. Entonces,

v(t) = 4,905 + 45,095e−2t .

Para hallar la ecuación de movimiento, primero resolvemos la ecuación

diferencial:
dx
v(t) = = 4,905 + 45,095e −2t
dt
donde, integrando respecto a t concluimos:

x(t) = 4,905t − 22,55e−2t + k. (3.18)

Como hemos considerado que x = 0 cuando t = 0, reemplazando en (3.18)

obtenemos k = 22,55. Por lo tanto, la ecuación de movimiento es

x(t) = 4,905t − 22,55e −2t + 22,55.

Puesto que el objeto se liberó a 500 m sobre el suelo, podemos determinar

el momento en que el objeto toca el suelo haciendo x(t) = 500:

500 = 4,905t − 22,55e −2t + 22,55

despejando t obtenemos:

t − 4.59e −2t = 97,3. (3.19)

 La ecuación (3.19) no se puede resolver de manera explícita en términos
de t. Para resolverla tomamos el hecho que e −2t será muy pequeño para t
grande, es decir ignoramos el término e −2t y obtenemos como aproximación
t ≈ 97,3 segundos. X


Ejemplo 3.13. Un paracaidista cuya masa es de 75 kg se arroja de un helicóptero

que vuela a 2.000 m sobre el suelo y cae hacia este bajo la influencia de la gravedad.
Suponga que la fuerza debida a la resistencia del aire es proporcional a la velocidad
del paracaidista, con la constante de proporcionalidad b1 = 30 N.s/m si el paracaídas
está cerrado, y b2 = 90 N.s/m cuando está abierto. Si el paracaídas no se abre hasta
que la velocidad del paracaidista es de 20 m/s. ¿Después de cuántos segundos llegará
al suelo?

Solución. Solo nos interesa el momento en que el paracaidista toca el suelo,

no el lugar. Así, únicamente consideramos la componente vertical de su
descenso. Para esto, necesitamos usar dos ecuaciones, una para describir el
movimiento antes de abrir el paracaídas y la otra para aplicar después de
abrirlo.
· Antes de abrirse: m = 75 kg, v0 = v(0) = 0, b = b1 = 30 N.s/m,
m
g = 9,81 seg 2 y x(0) = 0. Reemplazando en (3.15) tenemos:

dv
75 = 75 · 9,81 − 30v(t)
dt
dv m
= 9,81 2 − 0,4v(t)
dt seg
dv
+ 0,4v(t) = 9,81.
dt
Utilizando el método de resolución de ecuaciones de variables separa-
bles, o el hecho de que esta ecuación es lineal, resolvemos y hallamos
que
v1 (t) = 24,53 + ce−0,4t .
Como v(0) = 0, se obtiene que c = −24,53. Así obtenemos que
v1 (t) = 24,53 − 24,53e −0,4t y

x1 (t) = 24,53t + 61,31e −0,4t − 61,31.

Dado que el paracaídas se abre cuando se alcance la velocidad de

20 m/s, es decir la velocidad final, se puede hallar el tiempo que
recorrió:
20 = v(t1 ) = 24,53 − 24,53e −0,4t1 .
 Despejando de esta ecuación obtenemos el tiempo t1 :
t1 = 4,22s.
Luego, el paracaidista ha caído
x1 (4,22) = 24,53(4,22) + 61,31e −0.4(4,22) − 61,31 = 53,54
(en estos y otros cálculos, redondearemos nuestras cifras a dos cifras
decimales). Ahora, al abrirse el paracaídas, el paracaidista está a 2.000-
53,54=1.946,46 metros sobre el suelo.

· Consideramos el caso después de abrirse: m = 75 kg, v0 = v(0) = 20 ms ,

b2 = 90N · seg/m, g = 9,81 m
s2
y x2 (0) = 0; por lo tanto,

x2 (t) = 8,12t + 9,85 − 9,85e −1,2t .

Para determinar el momento en que el paracaidista llega al suelo, hacemos
x2 (t2 ) =1.946,46 m, la altura a la que estaba el paracaidista al abrir su
paracaídas. Esto da como resultado
1.946,46 = x2 (t2 ) = 8,12t2 + 9,85 − 9,85e −1,2t2
238,5 = t2 − 1,21e−1,2t2 . (3.20)
La ecuación (3.20) no se puede resolver de manera explícita en términos de
t. Para resolverla tomamos el hecho de que e −1,2t será muy pequeño para t
grande; esto significa que ignoramos el término e−1,2t y obtenemos como
aproximación t2 ≈ 238,5 segundos. Por lo tanto, golpeará el suelo 238,5 s
después de arrojarse desde el helicóptero. X


Ejercicios
1. Un objeto de 400 libras se libera desde el reposo a 500 pies sobre
el suelo y se le permite caer bajo la influencia de la gravedad. Supo-
niendo que la fuerza en libras debida a la resistencia del aire es −10v,
donde v es la velocidad del objeto en pies/s, determinar la ecuación de
movimiento del objeto. ¿En qué momento tocará el objeto el suelo?
2. Un objeto de masa 8 kg recibe una velocidad inicial hacia abajo de
20 m/s y luego se le permite caer bajo la influencia de la gravedad.
Suponga que la fuerza en newtons debida a la resistencia del aire
es −16v, donde v es la velocidad del objeto en m/s. Determinar la
ecuación de movimiento del objeto. Si el objeto está inicialmente a 100
m sobre el suelo, determinar el momento en que el objeto golpeará el
suelo.
3. Un objeto con masa de 2 kg se lanza desde el reposo de una plataforma
a 30 m sobre el agua y se deja caer bajo la influencia de la gravedad.
Después de que el objeto golpea el agua, comienza a hundirse, con la
gravedad jalándolo hacia abajo y una fuerza de flotación empujándolo
hacia arriba. Suponga que la fuerza de gravedad es constante, que
no hay cambios en el momento del objeto al golpear el agua, que la
fuerza de flotación es 1/2 del peso (peso = mg) y que la fuerza debida
a la resistencia del aire o del agua es proporcional a la velocidad del
objeto, con constante de proporcionalidad b1 = 10N.seg/m en el aire
y b2 = 30N.s/m en el agua. Determinar la ecuación de movimiento
del objeto. ¿Cuál es la velocidad del objeto al minuto después de ser
arrojado?
4. Un paracaidista cuya masa es de 100 kg se arroja de un helicóptero que
vuela a 3.000 m sobre el suelo y cae bajo la influencia de la gravedad.
Suponga que la fuerza debida a la resistencia del aire es proporcional a
la velocidad del paracaidista, con la constante de proporcionalidad b3 =
30 N.s/m cuando el paracaídas está cerrado y b4 = 50 N.s/m cuando
se abre. Si el paracaídas no se abre hasta 30 segundos después de que el
paracaidista sale del helicóptero, ¿después de cuántos segundos llegará
él al suelo? Si el paracaídas no se abre hasta 1 minuto después de que el
paracaidista sale del helicóptero, ¿después de cuántos segundos llegará
él al suelo?
5. Un objeto con masa de 100 kg se lanza desde el reposo de una lancha
hacia el agua y se deja hundir. Aunque la gravedad jala el objeto hacia
abajo, una fuerza de flotación de 1/40 veces el peso del objeto lo
empuja hacia arriba (peso = mg). Si suponemos que la resistencia del
agua ejerce sobre el objeto una fuerza proporcional a la velocidad del
objeto, con constante de proporcionalidad 10 N.s/m. Determinar la
ecuación de movimiento del objeto. ¿Después de cuantos segundos
ocurrirá que la velocidad del objeto es igual a 70 m/s?
6. Un objeto de masa m se libera desde el reposo y cae bajo la influencia
de la gravedad. Si la magnitud de la fuerza debida a la resistencia
del aire es bv n , donde b y n son constantes positivas, determinar la
velocidad límite del objeto (suponiendo que este límite existe). (Su-
gerencia: justifique que la existencia de una velocidad límite (finita)
implica dvdt → 0 cuando t → +∞).
7. Cuando la velocidad v de un objeto es muy grande, la magnitud de la
fuerza debida a la resistencia del aire es proporcional a v 2 y la fuerza
actúa en dirección opuesta al movimiento del objeto. Un proyectil con
masa de 3 kg se lanza desde el suelo hacia arriba y con una velocidad
 inicial de 500 m/s. Si la magnitud de la fuerza debida a la resistencia
del aire es 0,1v 2 , ¿en qué momento alcanzará el proyectil su máxima
altura sobre el suelo?. ¿Cuál es esa máxima altura?
8. Si un objeto cae desde un aeroplano, la velocidad hacia abajo después
de x seg es aproximadamente y 0 = gk (1 − e −kx ) en donde g = 9,8 m
s2
y
k = 0,2s . y(x) es la distancia que ha caído, por lo tanto y(0) = 0. Si
−1

este objeto cae desde una altura de 5.000 m, ¿cuánto le tomará llegar
a la tierra?

5. Trayectorias ortogonales
Dos curvas se cortan en ángulo recto si las respectivas rectas tangentes,
en los puntos de intersección, son rectas perpendiculares. En esta sección
utilizaremos ecuaciones diferenciales de primer orden para establecer cuando
dos familias de curvas son ortogonales; es decir, cuando cada curva de una
familia corta a las curvas pertenecientes a la otra familia en un ángulo recto.

Definición 3.1. Cuando todas las curvas de una familia F (x, y, c1 ) = 0

cortan ortogonalmente a todas las curvas de otra familia G (x, y, c2 ) = 0, se
dice que las familias son ortogonales. Esto significa que cada curva de la
familia F (x, y, c1 ) = 0 corta en ángulo recto a toda curva de la otra familia
G (x, y, c2 ) = 0.

Ejemplo 3.14. En la figura 3.4 se ve que la familia de rectas que pasan por
el origen y = c1 x y la familia de círculos concéntricos con centro en el origen
x 2 + y 2 = c2 son trayectorias ortogonales.
 Figura 3.4: Círculos concéntricos y las rectas y = c1 x.

Para encontrar las trayectorias ortogonales de una familia de curvas F (x, y, c1 ) =

dy
0, se halla de esta ecuación (por diferenciación implícita) la derivada dx , que usual-
mente se puede expresar de la forma siguente:

dy
= r (x, y)
dx
 
/𝜕x
r (x, y) = − 𝜕F
𝜕F /𝜕y , en todo punto (x, y) de F (x, y, c 1 ) = 0, donde 𝜕F
𝜕y ≠ 0 .
La ecuación de la familia ortogonal es dada, en forma diferencial, por la ecuación
diferencial
dy 1
=− .
dx r (x, y)
· Esta expresión está inspirada en el hecho que dos rectas son perpendiculares si
el producto de sus pendientes es −1.
Al eliminar la constante c1 entre el par de ecuaciones

F (x, y, c1 ) = 0
dy 1
=−
dx r (x, y)

y resolver la ecuación diferencial resultante (en caso de ser posible), obtenemos una
representación algebraica de la familia de curvas que es ortogonal a la familia
F (x, y, c1 ) = 0. X

Ejemplo 3.15. Hallar las trayectorias ortogonales de la familia de curvas y = e kx .

Solución. La derivada de y = e kx es
dy
= kekx .
dx
ln(y)
Despejamos k de y = e kx , ln(y) = kx, k = x , x ≠ 0 y reemplazamos:
dy ln(y) ln(y)
= e
dx x
dy y ln(y)
= .
dx x
La ecuación diferencial de la familia ortogonal es
dy x
=− .
dx y ln(y)
Esta última ecuación se resuelve por el método de resolución de ecuaciones
de variables separables:
dy
y ln(y) = −x
dx
w dy w
y ln(y) dx = −xdx
dx
y 2 ln(y) y 2 x2
− = − + c2
2 4 2
2y ln(y) − y + 2x 2 = 4c2
2 2

2y 2 ln(y) − y 2 + 2x 2 = c3 ,

donde por conveniencia se reemplazó 4c2 por c3 . X



Ejemplo 3.16. Hallar las trayectorias ortogonales a las curvas

y = cx n , n ∈ ℝ. (3.21)

Solución. La derivada de y = cx n es y 0 = cnx n−1 . Despejamos c de la

ecuación (3.21) de modo que tenemos c = yx −n , x, y ≠ 0. Reemplazando
concluimos que
dy
= yx −n nx n−1 = nx −1 y.
dx
dy x
La ecuación diferencial de la familia de curvas es dx = − ny . Esta última
ecuación es de variables separables, de donde tenemos
dy
ny = −x
dx
integramos la igualdad para obtener la solución:
dyw w
n y
dx = − xdx
dx
y 2 x2
n = − + c.
2 2
ny 2 x2
Así, las trayectorias ortogonales vienen dadas por la fórmula 2 + 2 = c.

Figura 3.5: Las curvas y = cx 2 para c = ±1, ±2 y sus trayectorias ortogonales.

Note que las trayectorias ortogonales incluyen, circunferencias, elipses e

hipérbolas. X


Ejercicios
1. Encontrar las trayectorias ortogonales a la familia:
a) x 2 + 2xy − y 2 = k.
b) x 2 + (y − k) 2 = k 2 .
c) y = −x − 1 + c1 e x .
d) y 2 = 4px.
e) xy = kx − 1.
1
f ) y = x+c 1
.
k
g) y = x , para x > 0, c > 0.
 h) xy = c.
i) e x cos(y) = k.
2
j) sin(y) = ke x .
k) x 2 + xy + y 2 = k.
2. Determinar el valor de a para que las familias de curvas y 3 = C1 x y
x 2 + ay 2 = C2 sean ortogonales.
3. Hallar las trayectorias ortogonales de cada una de las siguientes familias
de curvas
a) Una familia de líneas rectas que cortan el origen.
b) Una familia de círculos con radio variable, centradas en el eje x y
pasando a través del origen.
c) Una familia de elipses con centros en el origen y vértices en
(±1, 0).

Respuesta de ejercicios seleccionados

1. a) x 2 − 2xy − y 2 = c.
b) x 2 + y 2 = cx.
d) 2x 2 + y 2 = c.
e) x 3 + 3y = c.
i) e y sin(y) = c.
j) x = cos2 (y).
3. a) x 2 + y 2 = c.
b) x 2 + y 2 = cy.
2 2
c) x 2 = ce x +y .
2. Calefacción y enfriamiento de un edicio La calefacción y enfriamiento de
un edicio puede ser modelado por una ecuación diferencial

T 0 = k1 (T − Tα ) + k2 (T − Tω ) + P
donde T = T (t) es la temperatura del edicio en cualquier tiempo t, t, Tα la
temperatura fuera del edicio y Tω la temperatura deseada en el edicio, y P la
tasa de incremento de T debido a las máquinas y a las personas dentro del edicio,
k1 y k2 son constantes negativas. Solucione esta EDO, asumiendo que P = const,
Tw = const, y Ta varía sobre las 24 horas, digamos, Ta = A − Ccos(2π/24)t.
Discuta sobre el efecto que cada término de la ecuación tiene sobre la solución.

Solución
dT
= k1 (T − Ta ) + k2 (T − Tω ) + P
dt
dT
= k1 T − k1 Tα + k2 T − k2 Tω + P
dt
dT
= (k1 + k2 )T − k1 Tα − k2 Tω + P
dt
dT
− (k1 + k2 )T = −k1 (A − Ccos(2π/24)t) − k2 Tω + P (5)
dt
π


Utilizando factor integrante, aquí p(t) = (k1 + k2 ), r = k3 − k1 C cos t donde
12
k3 = k1 A + k2 Tω + P siendo k1 y k2 son constantes negativas, de esta forma
Z Z
h = p(t)dt = (k1 + k2 )dt = (k1 + k2 )t, (6)

así

eh = e(k1 +k2 )t . (7)

Por tanto la solución deseada es Solución Sea y la cantidad de aire

Z 
−h h
T (t) = e e r(t)dt + c , (8)

π



siendo r(t) = k1 A − C cos t + k2 Tω + P , en consecuencia
12
Z h   π  i 
−(k1 +k2 )t (k1 +k2 )t
T (t) = e e k1 A − C cos t + k2 Tω + P dt
12
π π



−144(k1 + k2 )2 k1 C cos 12 t − 12π(k1 + k2 )k1 C sin 12 t
=
144((k1 + k2 )2 + π 2 )
b
+ + C1 e−(k1 +k2 )t .
a
Mas ejercicios de modelación Matemática

1. Considere un tanque que es usado en ciertos experimentos de hidrodinámica,

Después de uno de los experimentos el tanque contiene 200L de una solución
de colorante con una concentración de 1g/L. Para preparar para el siguiente
experimento, el tanque es enjuagado con agua fresca uyendo a razón de 2 L/min,
la solución bien mezclada sale del tanque a igual velocidad. Encuentre el tiempo
en el cual la concentración de colorante en el tanque alcanza el 1 % de su cantidad
original.

Solución: Sea
q= cantidad de colorante en el tanque

t=tiempo en minutos

Debido al principio de entrada-salida se obtiene

dq
= ujo de entrada − ujo de salida ,
dt
donde

ujo de entrada = velocidad de entrada × concentración de entrada =0

ujo de salida = velocidad de salida × concentración de salida
cantidad
= velocidad de salida ×
volumen
q
= .
100
 De esta forma se tiene el siguiente modelo

dq q
=− , (14)
dt 100
utilizando que la concentración en el tiempo t=0 es de 1 se concluye que q(0) =
200. Solucionando este problema de valores iniciales por el método de separación
de variables se concluye
t
q(t) = 200e− 100 . (15)

Se desea encontrar el tiempo en el cual se tiene una concentración de 1% de su

cantidad original es decir una concentración de 0,01g/L, como

cantidad q
concentración = = ,
volumen 200
se inere que
q
0,01 = (16)
200
Combinando (15) y (16) se tiene

t
0,01 = e− 100 ,

de lo cual t = 460min.

2. Un tanque con una capacidad de 500 galones originalmente contiene 200 galones
de agua con 100 libras de sal en solución. El agua que contiene 1 libra de sal por
galón entra a una velocidad de 3 gal / min, y la mezcla se deja salir del tanque
a una velocidad de 2 gal /min. Encuentre la cantidad de sal en el tanque en
cualquier momento antes del instante en que la solución comienza a desbordarse.
Encuentre la concentración (en libras por galón) de sal en el tanque cuando está
en el punto de desbordamiento. Comparar esta concentración con la concentración
teórica límite si el tanque tenía una capacidad innita.

Solución Observe que del ejercicio se tienen los siguientes datos

Capacidad del tanque: 500gal.
Velocidad de entrada: 3gal/min.
Concentración: 1lb de sal.

Velocidad de salida: 2gal/min.

S(t): Cantidad de sal en el instante t.
 S(0): 100lb de sal.

Además de esto el cambio de volumen esta dado por la función V (t) = 200 + t
Utilizando nuevamente el principio de entrada-salida, se obtiene que la EDO que
modela la situación es la siguiente

dS S
=3−2 con condición inicial S(0) = 100, (17)
dt 200 + t
reorganizando se tiene que

dS S
+2 = 3,
dt 200 + t
lo cual se resuelve por el método de factor integrante donde
Z Z
2 2
p(t) = , r = 3 h(t) = p(t)dt = dt = 2ln(200 + t) ⇒
200 + t 200 + t
1
eh = e2ln(200+t) = (200 + t)2 y e−h = ,
(200 + t)2
y reemplazando en la fórmula general de la solución
Z 
−h h
S(t) = e e rdt + C ,

se obtiene Z 
1 2
S(t) = 3(200 + t) dt + C
(200 + t)2
así
C
S(t) = (200 + t) +
(200 + t)2
usando la condición inicial se obtiene que

C
S(0) = (200) + = 100 ⇒ C = −4 × 106
(200)2
por tanto se tiene que la cantidad de sal en el tanque en cualquier momento antes
del instante en que la solución comienza a desbordarse es

4 × 106
S(t) = (200 + t) − . (18)
(200 + t)2
 Para saber la concentración (en libras por galón) de sal en el tanque cuando está
en el punto de desbordamiento, se resuelve la ecuación V (t) = 500 ( para conocer
el tiempo que toma en llenarse), de la cual t = 300 y reemplazando en (18) se
obtiene que la cantidad de sal en el momento de desbordamiento es

S(300) = 484,

y por tanto la concentración es

cantidad 484 121

c= = = ≈ 0,97lb/gal.
volumen 500 125

En el caso en el que el tanque tuviese una capacidad innita, queremos conocer

que sucede con la concentración cuando t → ∞, esto es

S(t)
lı́m c(t) = lı́m = 1. (19)
t→∞ t−→∞ 200 + t
por lo tanto la concentración sería de 1 lb/gal.

3. Un tanque contiene 100 galones de agua y 50 onzas de sal. El agua contiene una
1 1
concentración de sal de (1+ sin t) oz/gal y entra al tanque a razón de 2 gal/min
4 2
y la mezcla del tanque sale a la misma razón.

(a) Encuentre la cantidad de sal en el tanque en cada tiempo t.

(b) Graque la solución para un largo periodo de tiempo, lo suciente para que
veas el comportamiento nal de la gráca.

(c) El comportamiento a lo largo del tiempo de la solución es una oscilación

sobre un cierto nivel constante. ¾Cuál es este nivel? ¾Cuál es la amplitud de
la oscilación?

Solución: Sea y(t) la cantidad de sal en el tanque en el instante t. El modelo es

y0 = sal que entra − sal que sale

1 1
Como el agua contiene (1 + sin t) oz/gal y entran 2 galones cada minuto, la
4 2
1 1 1
sal que entra es de (1 + sin t) oz/min. La sal que sale será 2 gal/min × y.
2 2 100
Por tanto, el modelo es
 
0 1 1 1
y = 1 + sin t − y
2 2 50
y la condición inicial es y(0) = 50.

Reorganizando la ecuación diferencial tenemos que

 
1
0 1 1
y + y= 1 + sin t
50 2 2
1 1
R
dt
El factor integrante es µ=e 50 = e 50 t . Multiplicando la EDO por µ, tenemos
que  
1 1t 0 1
(µy) = e 50 1 + sin t
2 2
Integrando a ambos lados, nos queda que

25et/50 (− sin t + 50 cos t − 5002)

µy = − +C
5002
Luego,
25(− sin t + 50 cos t − 5002) C
y=− + t/50
5002 e
Aplicando la condición inicial y(0) = 50, tenemos que

y(0) = 24,75 + C = 50

Luego, C = 25,25. Es decir,

25(− sin t + 50 cos t − 5002) 25,25

y=− + t/50
5002 e
Usando Geogebra, vemos que la gráca es similar a una exponencial decreciente:
 Sin embargo, haciendo un zoom, vemos las oscilaciones alrededor del nivel 25, con
amplitud de
25
5002
= 5 × 10−3 .

4. Suponga que una suma S0 es invertida a una tasa de rendimiento anual r com-
puesta continuamente.

a) Encuentre el tiempo T, en función de r, requerido para que la suma original

se duplica

Solución Sean
S = cantidad de dinero.
S0 = cantidad de dinero invertida inicialmente.

t=tiempo medido en años.

El modelo para la inversión de capital viene dado por

dS
= rS, S(0) = S0 , (20)
dt
solucionando este problema de valores iniciales se llega a

S(t) = S0 ert .

Se desea conocer el tiempo T en el cual S(T ) = 2S0 es decir

2S0 = S0 erT ,
 despejando se obtiene
1
T = ln 2. (21)
r
b) Determine T si r = 7 %.
Solución En el caso que r = 7 % = 0,07, reemplazando en (21) se llega a
1
T = ln 2.
7

c) Encuentre la tasa de rendimiento anual si la inversión inicial se duplica en 8

años.

Solución En este caso T = 8, despejando r de la ecuación (21) se tiene

1 1
r= ln 2 = ln 2.
T 8
5. Un joven sin capital inicial invierte k dólares anuales a una tasa anual r. Suponga
que las inversiones se realizan continuamente y que el interés también se compone
de manera continua.

a) Determine la suma S(t) acumulada en cualquier momento t.

Solución Sean
S = cantidad de dinero.
S0 = cantidad de dinero invertida inicialmente.

t=tiempo medido en años.

La ecuación diferencial que modela la situación es la siguiente

dS
= k + rS
dt
reorganizando y solucionando por factor integrante se tiene que

dS
− rS = k
dt
Z
donde p(t) = −r, h(t) = −rdt = −rt,
Z 
−rt −rt k
S(t) = e ke dt + c = − + Cert
r
 teniendo en cuenta que en t = 0, S(0) = 0 por que está sin capital, entonces
C = kr , de ahí que la suma S(t) en cualquier momento t es

k rt
S(t) = (e − 1). (22)
r
b) Si r = 7,5 %, determine k de modo que un millón de dólares esté disponible
para la jubilación en 40 años.

Solución Del inciso 5a . se tiene que

Sr
k=
ert −1
donde r = 0,075, t = 40, S = $1000000
(1000000)(0,075)
k= = 3929, 68.
e(0,075)(40) − 1

c) Si k = $2000 anuales, determine la tasa de interés r que se debe obtener

para contar con un millón de dólares disponibles en 40 años.

Solución En esta caso S = $1000000, k = 2000, t = 40 de (22) se tiene que

r = 0,0977 = 9,77 %.

6. Un cierto estudiante pide prestado $8000 para comprar un carro. El prestamista

carga intereses a una tasa anual de 10 %. Suponiendo que el interés se carga conti-
nuamente y que el deudor realiza pagos continuamente a una tasa anual constante
k, determina la tasa de pago k que se requiere para pagar el préstamo en 3 años.
También determine cuántos intereses se pagan durante el período de 3 años.

Solución:
Sea y(t) lo que paga el deudor al cabo de t años.

y 0 = ry + k

donde r es la tasa de interés y k los pagos. Es decir,

y 0 = 0, 1y + k

o
y 0 − 0,1y = k.
 El factor integrando es µ = e−0,1t . Luego,

y = Ce0,1t − 10k

con C una constante arbitraria.

Usando la condición inicial y(0) = 0 obtenemos que

y(0) = C − 10k = 0

Es decir, C = 10k . Luego,

y = 10ke0,1t − 10k
Después de 3 años, tenemos y(3) = 8000, es decir,

y(3) = 10ke0,3 − 10k = 8000

Despejando, tenemos k = 2286,64. Es decir, se requieren $2286 anual para pagar

el préstamo en 3 años.

La cantidad de intereses al cabo de tres años es 0,1y(3), es decir,

0,1(10 × 2286,64e0,3 − 10 × 2286, 64) = 800

Luego, la cantidad aproximada de intereses que paga en tres años es de $800.

7. Una bola de masa 0.15 kg es lanzada hacia con una velocidad inicial de 20m/s
desde el techo de un edicio de 30m de altura. Despreciando la resistencia del
aire.

a) Encuentre la máxima altura sobre el suelo que alcanza la pelota.

Solución Sean
y = altura con respecto al tiempo.
v = velocidad con respecto al tiempo
g = gravedad.
t = tiempo en segundos
Utilizando tercera ley de Newton se obtiene F = ma, de lo cual se obtiene
el modelo
dv
= −g,
dt
así v = −gt + C . Como v(0) = 20 se llega a C = 20, así

v(t) = 20 − 9,8t.
 Se integra la anterior ecuación para obtener y, así

y(t) = 20t − 4, 9t2 + K,

sabemos que y(0) = 30, de donde k = 30, por tanto la función para la altura
es
y(t) = 20t − 4, 9t2 + 30. (23)

Para saber la máxima altura hacemos v(t) = 0, de esta forma tenemos que
t = 2, 04. La altura en este instante de tiempo es

y(2, 04) = 50, 4m

b) Asumiendo que la bola no golpea el edicio en su caída, encuentre el tiempo

que golpea el suelo.

Solución Para saber el tiempo en el cual la bola golpea el suelo se iguala la

ecuación (23) a cero, así

−4,9t2 + 20t + 30 = 0,

solucionando esta ecuación t ≈ 5, 25.

c) Realice la gráca de la velocidad y posición con respecto al tiempo.

Solución

8. Supongamos que las condiciones son como en el problema 7, excepto que hay una
fuerza debida a la resistencia del aire de magnitud |v|/30 dirigida opuesta a la
velocidad, donde la velocidad v se mide en m / s.

a) Encuentre la altura máxima sobre el suelo que la bola alcanza.

Solución Según el ejercicio 7, las condiciones dadas son las siguientes

 Masa: 0.15 kg

v(0) = 20m/s
h=30m,

además de esto hay una fuerza debida a la resistencia del aire de magnitud
v/30, por lo tanto la ED que modela la situación es la siguiente

dv v
m = −mg − ,
dt 30
separando variables queda

dv dt
=− ,
v + 44,1 4,5
solucionando
t
ln(v + 44,1) = + C,
4,5
t
v(t) = Ce− 4,5 − 44,1,
utilizando la condición inicial se tiene que C = 64,1 por lo tanto la velocidad
en el instante t es
t
v(t) = 64,1e− 4,5 − 44,1. (24)

Al integrar (24) se observa que la función para la altura es

t
x(t) = −288,45e 4,5 − 44,1t + K, (25)

como x(0) = 30 se tiene que K = 318,45, de lo cual

t
x(t) = 288,45e 4,5 − 44,1t + 318,45. (26)

Ahora para determinar la máxima altura que alcanza la bola se debe tener
en cuenta que en este punto v(t) = 0, solucionando esta ecuación se llega
a que t = 1,68sg , reemplazando este valor en (26) se obtiene que la altura
máxima es x(1,68) = 45,78m.
b) Encuentre el tiempo que la bola golpea el suelo.

Solución El tiempo en que la bola golpea al piso x(ts ) = 0 por lo que se

tiene
ts
288,45e 4,5 − 44,1ts + 318,45 = 0
lo cual da ts = 5,13 sg
 c) Trazar los grácos de velocidad y posición en función del tiempo.

Solución
9. Suponga que se dan las condiciones del problema 7 excepto que existe una fuerza
2
debida a la resistencia del aire de magnitud v /1325 en dirección opuesta a la
velocidad, donde la velocidad v se mide en m/s.

(a) Encuentre la altura máxima sobre el suelo que alcanza la bola

(b) Encuentre el tiempo en el que la bola golpea el suelo

(c) Dibuja las grácas de velocidad y posición versus tiempo. Compara estas
grácas con las de los ejercicios 7 y 8.

Solución: Sea y(t) la altura que alcanza la bola en el instante t. El peso de la

dv
bola es w(y) = −mg La aceleración que lleva la bola es dt
, luego

dv v2
m = −mg −
dt 1325
 v2
Es decir, v 0 = −g − 1325m
. Separando variables,

1
v2
v 0 = −1
g+ 1325m

Reemplazando m = 0,15 kg y g = 9,8 m/s2 ,

1
v2
v 0 = −1
9,8 + 198,75

Reescribiendo,
1
198,75 v 0 = −1
1947,45 + v 2
Integrando,  
v
4,5 arctan = −t + C
44,13
siendo C una constante arbitraria.

Despejando v, tenemos que

 
−t
v = 44,13 tan +C
4,5
Usando la condición inicial v(0) = 20, tenemos

v(0) = 44,13 tan(C) = 20

Despejando C,
C = arctan(0,45) = 0,42
Por tanto,  
−t
v = 44,13 tan + 0,42
4,5
 
−t
Como v = y0, tenemos
0
que y = 44,13 tan
4,5
+ 0,42 . Integrando tenemos que

  
−t
y = 198,585 ln cos + 0,42 +D
4,5
siendo D una constante arbitraria. Usando la condición inicial y(0) = 30, tenemos

y(0) = −18,55 + D = 30
Entonces D = 48,55. Así,

  
−t
y = 198,585 ln cos + 0,42 + 48,55
4,5

La altura es máxima cuando v = 0:

 
−t
0 = tan + 0,42
4,5

Tenemos que t = 1,9 segundos. La altura en este instante de tiempo es

y(1,9) = 48,55m

Es decir, la altura máxima es aproximadamente 48,55m.

La bola golpea el suelo cuando y = 0, es decir,

  
−t
0 = 198,585 ln cos + 0,42 + 48,55
4,5

Despejando el coseno, tenemos que

 
−t
cos + 0,42 = 0,78
4,5

Aplicando cos−1 , tomando el ángulo negativo, nos queda

−t
+ 0,42 = −0,67
4,5
Es decir, t = 4,9 segundos. Tenemos entonces que la bola golpea el suelo a los 4,9
segundos.

Las grácas de velocidad vs tiempo y distancia vs tiempo son, respectivamen-

te