Capitulo 4
MODELOS
Y
PREDICCIÓN
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
Uno de los principales objetivos en la enseñanza del Cálculo para un alumno de bachillerato,
es que pueda manipular las operaciones que consideramos básicas en esta área de las
matemáticas: la derivación y la integración. Sin embargo también se ha descuidado la relación
de éstas con problemas de aplicación que muestren su utilidad; contemplado esto dentro del
Modelo del Colegio, representa una contradicción con los principios que lo han regido y este
punto es una de las razones principales de la existencia del capítulo 4.
En este sentido, poco se ha hecho en diversos textos sobre la asignatura que se refieran a un
enfoque más aplicado del Cálculo, a no ser en textos referidos más a áreas específicas como
Biología, Economía, Física y que se estudian para alumnos que se encuentran en nivel
licenciatura.
Uno de los objetivos que se persigue en el programa de estudios, es que el alumno pueda
introducirse a los temas que sean de su interés observando que, conociendo y manejando los
conceptos de derivación e integración, la información que puede obtener respecto a un
problema es además de importante y relevante para casos más apegados a la realidad.
Se presenta una clasificación de los problemas que se contemplarán en este capítulo en un
sentido cualitativo, dejando abierta la posibilidad de caracterizarlos por áreas de estudio
(ciencias biológicas, físicas, económicas). Comenzamos por presentar un caso sencillo de lo
que llamamos en los capítulos anteriores ecuación diferencial, donde es importante mencionar
que este concepto es un área basta e interesante de las matemáticas que requiere del
conocimiento de la derivada y de la integral y que al mismo tiempo aporta importantes
resultados en las matemáticas aplicadas.
Dado que la familiarización del alumno con los métodos de integración es muy temprana, las
limitaciones respecto a estudiar otro tipo de ecuaciones diferenciales es clara. En este capítulo
restringiremos nuestro objetivo a establecer una relación entre este modelo de ecuaciones y los
problemas de aplicación.
GRUPO 401C 95
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
dF
PROBLEMAS QUE INVOLUCRAN EL MODELO = kF
dt
Problema 4.1. Para cargar un capacitor a
través de una resistencia, se utiliza el circuito
mostrado en la figura, la relación matemática
entre el voltaje proporcionado por la batería y
el que se establece en la resistencia y el
capacitor es
dQ
CV − Q =
RC
dt
C indica la capacitancia del capacitor; R la resistencia del resistor, V el voltaje
proporcionado por la batería y t el tiempo durante el cual se esta cargando el capacitor.
dQ
Solución. La razón de cambio instantánea de la carga del capacitor respecto al tiempo
dt
dQ CV − Q
es: =
dt RC
dQ
Observar que la rapidez de variación de la carga depende de la misma carga Q , es por
dt
dQ
ello que el buscar graficar es prácticamente imposible, se requiere conocer la función
dt
explícita de cómo varía la carga respecto al tiempo Q ( t ) .
Para ello debemos utilizar un método matemático que permita despejar a la variable Q para
obtenerla en términos del tiempo t .
La expresión (1) nos permite calcular la diferencial de Q ( dQ ) de la siguiente manera:
dQ CV − Q CV − Q
dQ = dt que a su vez es igual a dt ; esto es dQ = dt , se observa que
dt RC RC
en la expresión la variable Q aparece en ambos lados de la expresión, la cual se puede
escribir como:
dQ
= RCdt
CV − Q
Observar que en esta expresión fue posible escribir del lado izquierdo lo que contiene a la
variable Q y del lado derecho lo que contiene la variable t . Es decir se hizo una separación de
variables.
Para obtener como varía la carga Q en función del tiempo debemos integrar con lo que se
denomina las condiciones iniciales. Para este problema suponemos que al tiempo t = 0 , la
carga en el capacitor es Qo = 0 y al tiempo t la carga es Q .
GRUPO 401C 96
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
Q t
dQ t
∫0 CV − Q = ∫0 RCdt ; de donde − ( In ( CV − Q ) − InCV ) =
RC
Despejando Q , se tiene=Q CV 1 − e ( −t
RC
)
Problema 4.2. La temperatura de un cuerpo T rodeado por aire a temperatura To varía de
modo que el ritmo de variación de su temperatura es proporcional a la diferencia de
temperaturas (Ley del enfriamiento de Newton). Encontrar la temperatura T en función del
tiempo t .
Solución. La expresión matemática que representa a este proceso es
dT
= k (T − T0 ).
dt
dT
Se tiene que dT
= dt k (T − T0 )dt , por lo tanto
=
dt
dT
= kdt ,
T − T0
Observar que en esta expresión fue posible escribir del lado izquierdo lo que contiene a la
variable T y del lado derecho lo que contiene la variable t . Es decir se hizo una separación de
variables.
Integrando se tiene ln T − T0 = kt + c, donde c es una constante
Por tanto T − T0 = Ae kt , A>0.
Como inicialmente T > T0 , entonces T − T0 siempre positivo y, por tanto, T − T0 = T − T0 .
Podemos eliminar así las barras verticales del valor absoluto, y se tiene
T − T0 = Ae kt .
Problema 4.3. Si la cantidad inicial de un fermento es de 1g y al cabo de 1 hora es de 1.2g,
¿cuánto fermento habrá al cabo de 5 horas de comenzar la fermentación si la velocidad de
incremento del fermento es proporcional a la cantidad existente?
Solución. La expresión matemática que expresa al problema es
dC
= kC ,
dt
dC
donde C (t ) representa la cantidad de fermento en el instante t. =
y dC = dt kC ,
dt
dC
Es una ecuación de variables separables de la forma = kt , integrando se tiene
C
GRUPO 401C 97
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
dC
∫ C
= ∫ kdt + α , α cte, de integración.
De donde InC= kt + α y C = Ae kt , donde A = eα > 0.
Para que la integral exista, es necesario que C > 0, o que C < 0, para todo valor de t. También
es solución la función C (t ) = 0, , para todo valor de t, como se comprueba por sustitución
directa en la ecuación diferencial. Las tres posibilidades se pueden recoger en una única
kt
expresión C (t ) = Ae ,
Si se admite que A pueda tomar valores reales cualesquiera. En este problema, el valor inicial
C (0) = 1 , lo que provoca que A = 1. Por tanto C (t ) = e kt .
Para obtener el valor de la constante k , se utiliza el dato de que C (1) = 1.2 . Sustituyendo en
k
la ecuación resulta 1,2 = e . por lo tanto k = In1.2 y la solución particular buscada es
C (t ) = et ln 1, 2 . La cantidad de fermento al cabo de 5 horas será
C (5) = e5 ln 1, 2 = 1,25 ≈ 2,49 g.
Problema 4.4. Cuatro insectos se encuentran cada uno en una esquina de una mesa cuadrada
de lado a. En el instante t = 0 empiezan a moverse, cada uno en la dirección del que tiene a
su derecha. Encontrar: a) La ecuación de la trayectoria seguida por los insectos hasta
encontrarse en el centro de la mesa.
Solución. En la mesa se pueden definir coordenadas polares con el polo en el centro de la
mesa y el eje polar pasando por uno de Ios vértices del cuadrado.
Por la simetría del problema es suficiente con obtener la trayectoria de uno de los insectos, por
ejemplo el que se encuentra inicialmente sobre el eje polar, puesto que la de los otros se
obtiene mediante rotación de la primera trayectoria desde centro de la mesa y ángulos
π 3π
de , π y radianes, respectivamente.
2 2
El ángulo Ψ que forma el radio vector con la dirección en que se mueve el insecto, es siempre
3π
de radianes, ya que, por simetría, el radio vector es siempre una diagonal del cuadrado
4
que forman los insectos en todo instante. Por tanto,
dθ 3π
r = tg ( ) = −1.
dr 4
GRUPO 401C 98
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
dθ dr
Pero dθ = dr , por lo tanto se tiene la separación de variables = −dθ , ahora
dr r
−θ
integrando se tiene r = Ae .
2 2 −θ
La condición inicial es r (0) = a . Sustituyendo queda finalmente r = a e .
2 2
Problema 4.5. Un esquiador se desliza por una pendiente alcanzando una velocidad
v0 en el instante t = 0, en el que la pista se vuelve horizontal. Desde dicho punto
hasta que finalmente se detiene hay 25 m. Sabiendo que la resistencia del aire es
proporcional a la velocidad del esquiador, que toda otra resistencia al movimiento del
esquiador es despreciable y que en el primer segundo, es decir en t = 1, recorre 10
m, se pide calcular la velocidad v0 .
Solución. La segunda ley de Newton nos permite establecer la relación entre la fuerza
y la velocidad
dv
m = kv, k es una constante; m es la masa y v la velocidad
dt
dF dv
Esta expresión es de la forma = kF , y=mdv m= dt mkvdt ,
dt dt
dv
= αdt , donde α = mk es constante, la expresión muestra nuevamente una
v
separación de variables
Integrando se obtiene la solución general v(t ) = Keαt , K una constante.
La condición inicial v(0) = v0 , al sustituirla en la solución general, nos proporciona la
solución particular
v = v0eαt .
Por otra parte, integrando la función velocidad respecto del tiempo, obtenemos el
espacio recorrido
eαt
x(t ) = ∫ v(t )dt = v0 + C , C es una constante
α
,
v0
y, como x(0) = 0, se sigue que C = − . Sustituyendo este valor resulta la ecuación
α
del movimiento
v0
x(t ) =
α
(e αt
−1 )
que, en particular, para t = 1 resulta
v0
10 =
α
(e αt
)
−1
GRUPO 401C 99
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
αt v0
Como e = , la ecuación del movimiento en función de la velocidad es
α
v0 v0
x (v ) = − 1.
α α
v0
Por tanto, como x = 25 cuando v = 0 , se tiene 25 = − , con lo que la ecuación
α
queda en la forma 10 = −25 e ( αt
)
− 1 , lo que permite despejar el valor de la constante
3
α = ln . Finalmente,
5
3
v0 = −25α = 25 ln = 12,77 m / s.
5
Problema 4.6. Un día comenzó a nevar por la mañana y siguió nevando de forma constante
durante todo el día. A mediodía, una máquina quitanieves comenzó a limpiar una carretera a
ritmo constante, en términos de cantidad de nieve retirada en cada hora. A las 2 de la tarde la
máquina había limpiado 2 km. de carretera y a las 4 de la tarde 1 km. más. ¿A qué hora
empezó a nevar?
RESOLUCIÓN. Consideremos que al tiempo t = 0 comienza a nevar. Al cabo de un tiempo t
la altura de la nieve es h = at ( a constante). Si la máquina quita un volumen V de nieve
constante por unidad de tiempo, la anchura de la pala es L, y designamos por x la longitud
recorrida por la máquina desde que empezó a trabajar, el volumen que quita
instantáneamente por unidad de tiempo es constante, entonces
dV dx
= hL = c, ( c constante)
dt dt
c dx
definiendo α= y utilizando dx = dt , se puede escribir , al tener una separación de
aL dt
variables se puede proceder a integrar.
2 t0 + 2
dt
Si se considera que al mediodía t = t0 , entonces ∫ dx = α
0
∫
t0
t
por lo que
t0 + 2
2 = α ln .
t0
3 t0 + 4
dt t0 + 4
y también ∫ dx = α
0
∫
t0
t
, entonces 3 = α ln
t0
.
Dividiendo ambas expresiones y tomando antilogaritmos, resulta
GRUPO 401C 100
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
t0 + 4
ln 2 3
t0 3 t0 + 4 t + 2
= ⇒ = 0 .
t +2 2 t0 t0
ln 0
t0
2
Al simplificar se llega a la ecuación t0 + 2t0 + 4 = 0. Como t0 > 0 la solución es t0 = 5 − 1.
En conclusión, comenzó a nevar 5 − 1 horas antes del mediodía.
Problema 4.7. Un depósito contiene 10 litros de salmuera con dos gramos de sal disueltos. En
un cierto instante comienza a entrar en el depósito salmuera que contiene I gramo de sal por
litro, a razón de 3 litros por minuto y al mismo tiempo, a salir del depósito la disolución,
supuesta homogénea en todo momento, a razón de 4 litros por minuto. Hallar la cantidad de sal
x(t ) que contiene el depósito en un instante cualquiera t.
RESOLUCIÓN. El volumen de agua en el depósito en el instante t es 10 − t litros, t en
minutos. Luego, en dicho instante, la concentración de sal en el depósito es de
x
g / I.
10 − t
x
Si salen 4 litros, saldrán instantáneamente 4 gramos de sal por minuto.
10 − t
dx
El aumento instantáneo de sal por unidad de tiempo en el depósito en el instante t es
dt
4x
Como entran en ese instante 3g/I y salen entonces:
10 − t
dx 4x
= 3− .
dt 10 − t
En esta expresión, no es posible aplicar el método que hemos venido utilizando en problemas
anteriores, ya que:
dx 4x 30 − 3t − 4 x
Cuando dx = 3−
dt = dt = dt , se tiene una expresión donde no es
dt 10 − t 10 − t
posible realizar una separación de variables como anteriormente se hizo. La solución utilizando
el procedimiento de solución de ecuaciones diferenciales es:
x(t ) = 10 − t + k (10 − t ) 4 .
La solución particular, correspondiente al valor inicial x(0) = 2, es
[ ]
x(t ) = (10 − t ) 1 − 8 ⋅ 10−4 (10 − t )3 , 0 ≤ t ≤ 10.
dF dF
A la expresión del tipo = kF o 0 donde k es una constante, se le conoce
− kF =
dt dt
GRUPO 401C 101
CALCULO DIFERENCIAL E INTEGRAL II
MODELOS Y PREDICCION
como ecuación diferencial homogénea de orden uno, y su solución la podemos obtener
utilizando el método de separación de variables, obteniéndose así una solución del tipo
F = Ae kt .
La justificación rigurosa de todos los métodos, más o menos elementales, de resolución de
ecuaciones diferenciales presentados en los problemas anteriores, no es parte del presente
texto, sin embargo en el siguiente problema se plantea esta justificación en el caso particular
que acabamos de exponer, el de las ecuaciones de variables separables.
Ejemplo 4.1. Justificar rigurosamente el método de separación de variables
dx
aplicando a la ecuación g ( x ) = f (t ),
dt
dx
Solución. Suponiendo que la ecuación diferencial g ( x) = f (t ), tiene solución
(aunque ésta dt
no se conozca), que denotaremos por x = ϕ (t ), tenemos que g (φ (t )φ´(t )) = f (t )
en un cierto intervalo I de t . Podemos efectuar la integración de ambos miembros en
dicho intervalo, obteniendo
∫ g (φ (t )φ=
´(t ))dt ∫ f (t )dt + c, c constante de integración.
El cambio de variable x = ϕ (t ) transforma el primer miembro en
∫ g ( x)dx = ∫ f (dt )dt + c
que define, en general implícitamente, la solución general de la ecuación diferencial
dada.
■ NOTA Las hipótesis anteriores no son suficientes para poder aplicar el teorema de
existencia y unicidad de soluciones. Aún así, podría ocurrir que la solución obtenida
fuese única, ya que las condiciones del teorema son suficientes, pero no necesarias.
El siguiente ejemplo muestra que con las hipótesis del Ejemplo 1.3 la unicidad no
está garantizada.
GRUPO 401C 102