../..

/22519778_294718504345109_1361201845594081250_

DEPARTAMENTO DE:

CIENCIAS EXACTAS

CARRERA DE:

MECATRÓNICA

Tarea

NOMBRE:

ˆTUGULINAGO LEMA WELLINGTON SANTIAGO

DOCENTE:

ROMAN VARGAS WILSON MARCELO

NRC:

4203

LATACUNGA

Noviembre 2020- Abril 2021

1
 Tarea 3. Modelado con EDO de primer orden
SANTIAGO TUGULINAGO
24 de junio de 2021

Índice
1. Resolver múltiplos de 3 a partir del ejercicio 7 del libro de Dennis Zill Cap
1.3 3
1.1. Ejericio 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2. Ejericio 12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3. Ejericio 15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4. Ejericio 18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.5. Ejericio 21 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6. Ejericio 24 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.7. Ejericio 27 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.8. Ejericio 30 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.9. Ejericio 33 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.10. Ejericio 36 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

2. Resolver múltiplos de 3 a partir del ejercicio 7 del libro de Dennis Zill Cap
3.1 11
2.1. Ejercicio 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2. Ejercicio 12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.3. Ejercicio 15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.4. Ejercicio 18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.5. Ejercicio 21 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.6. Ejercicio 24 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.7. Ejercicio 27 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.8. Ejercicio 30 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.9. Ejercicio 33 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.10. Ejercicio 36 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.11. Ejercicio 39 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.12. Ejercicio 42 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
2.13. Ejercicio 45 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.14. Ejercicio 48 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31

3. REFERENCIAS 33

2
1. Resolver múltiplos de 3 a partir del ejercicio 7 del
libro de Dennis Zill Cap 1.3
1.1. Ejericio 9
Suponga que un tanque grande de mezclado contiene inicialmente 300 galones de agua en
los que se disolvieron 50 libras de sal. Entra agua pura a una razón de 3 gal/min y cuando
la solución está bien revuelta, sale a la misma razón. Determine una ecuación diferencial que
exprese la cantidad A(t) de sal que hay en el tanque al tiempo t. ¿Cuánto vale A(0)?

Resolución:

Datos: Desarrollo
dA
V0 = 300 gal = Ce Fe − Cs Fs
dt
dA gl A Fe gl
A = 50 lb t=0 = (3 )(0) − ( + )(3 )
dt min 300 gl min
gal dA A
Fe = 3 min
Ce = 0 =−
dt 100
gal
Fe = 3 min
Ce =? Conclusión
Pregunta: Por lo tanto tenemos que la cantidad de
dA
=? sustancia cuandoA(0) es igual a 50 lb
dt

1.2. Ejericio 12
Generalice el modelo dado en la ecuación (8) de la página 23, suponiendo que el gran tanque
contiene inicialmente N0 número de galones de salmuera, r entra y r sale son las razones de
entrada y salida de la salmuera, respectivamente (medidas en galones por minuto), centra es la
concentración de sal en el fl ujo que entra, c(t) es la concentración de sal en el tanque ası́ como
en el fl ujo que sale al tiempo t (medida en libras de sal por galón), y A(t) es la cantidad de
sal en el tanque al tiempo t.

Resolución:

Datos: Desarrollo:
dA
V0 = No gal = Ce F e − Cs F s
dt
dA lb A(t)
A0 = A(t) = (rentra )(Centra ) min −rsalida
dt N0 + (rentrada − rsalida )t
gal
Fe = rentrada min Ce = Centrada
gal
Fs = rsalida min Cs = Csalida Conclusión
Pregunta: Por lo tanto tenemos que la ecuacion diferencial
dA
=? nos queda en función de las variables del problema
dt

3
1.3. Ejericio 15
Un circuito en serie tiene un resistor y un inductor como se muestra en la fi gura 1.3.13.
Determine una ecuación diferencial para la corriente i(t) si la resistencia es R, la inductancia
es L y el voltaje aplicado es E(t).

Figura 1: FIGURA:1.3.13. Circuito LR del problema 15

Resolución:

Datos. Pregunta
Corriente: i(t) i(t) =?
Resistencia: R
Inductancia: L
Voltaje aplicado: E(t)

Desarrollo
Aplicando Ley de Kirchof.

El voltaje aplicado E(t) en un circuito cerrado debe ser igual a la suma de los voltajes en el
circuito, en consecuencia.

sabemos que la corriente es:

dq
i=
dt
Caı́da de volatje en la resistencia.
 Aplicando
 ley de Ohm.
dq
V = Ri = R
dt
Caı́da de voltaje
  en
 2 inductor
el 
di dq
L =L
dt dt2
Obtenemos la ecuacion
  diferencial
di
E(t) = L + Ri
dt

Conclusión: Tenemos como resultado final una ecuacion diferencial de primer orden

1.4. Ejericio 18
Un barril cilı́ndrico de s pies de diámetro y w lb de peso, está fl otando en agua como se
muestra en la fi gura 1.3.16a. Después de un hundimiento inicial el barril presenta un movi-
miento oscilatorio, hacia arriba y hacia abajo, a lo largo de la vertical. Utilizando la fi gura
1.3.16b, defi na una ecuación diferencial para establecer el desplazamiento vertical y(t), si se
supone que el origen está en el eje vertical y en la superfi cie del agua cuando el barril está

4
en reposo. Use el Principio de Arquı́medes: la fuerza de fl otación o hacia arriba que ejerce el
agua sobre el barril es igual al peso del agua desplazada. Suponga que la dirección hacia abajo
es positiva, que la densidad de masa del agua es 62.4 lb/pies3 y que no hay resistencia entre el
barril y el agua.

Datos: Pregunta:
lb dy
Densidad: ρagua = 62,4 =?
pies3 dt
Diámetro del Barril: d = s
s
radio: r =
2
Peso:W lb de peso.
De acuerdo a la figura 1.3.16.b podemos establecer que el volumen del cilindro es:

 s 2
v=π y
2
De acuerdo al principio de Arquı́medes:la fuerza de empuje o hacia arriba que ejerce el agua
sobre el barril es igual al peso del agua desplazada. Por tanto la fuerza de empuje sobre un
objeto = peso del fluido desplazado por el objeto.

Fempuje = ρV aguaderramada
   
lb s 2
Fempuje = 62,4 π y
pies3 2
Fempuje = 15,6πs2 y N

Aplicando la Segunda Ley de Newtón. F = ma

d2 y
 
m = −mg
dt2
Recordar que: W = mg
 2 
dy
m = −W (1)
dt2
Aplicando la definición téorica del principio de Arquı́mides
Fempuje = w = 15,6πs2 y
Por tanto remplazando
 2  en la ecuación (1)
dy
m 2
= −15,6πs2 y
dt

5
Se puede apreciar que nuestra ecuación diferencial obtiene forma

w
Remplazando: m =
g
 2 
w dy
+ 15,6πs2 y = 0
g dt2
Donde nuestra
 2  ecuación diferencial es:
dy 15,6πs2 y
=− g
dt2 w
Conclusión: Tenemos como resultado final una ecuacion diferencial de segundo orden

1.5. Ejericio 21
De acuerdo con la ley de la gravitación universal de Newton, la aceleración de caı́da libre a
de un cuerpo, tal como el satélite que se muestra en la fi gura 1.3.18, que está cayendo desde una
gran distancia hacia la superficie no es la constante g. Más bien, la aceleración a es inversamente
proporcional al cuadrado de la distancia desde el centro de la Tierra a = k/r2 donde k es la
constante de proporcionalidad. Utilice el hecho de que en la superficie de la Tierra, r = R y
a = g, para determinar k. Si la dirección positiva se considera hacia arriba, utilice la segunda
ley de Newton y la ley de la gravitación universal para encontrar una ecuación diferencial para
la distancia r.
Resolución:

Datos
k
a= r2

r=R
a=g
k = constante
Desarrollo

Tenemos que:
k
a= 2
r
k
g= 2
r
Despejamos K

6
 k = gr2
Segunada ley de newton
F = ma
Reemplazamos a teniendo encuenta el signo negativo ya que la la masa se desliza hacia abajo
ma = −mg
d2 y
m = −mg
dt2
d2 y
= −g
dt2
Tenemos que g = a
d2 y
= −a
dt2
Remplazamos a
d2 y k
= −
dt2 r2
2
dy gR2
= −
dt2 r2
2 2
d y gR
+ 2 =0
dt2 r
Conclusión: Tenemos una ecuación no lineal de segundo orden

1.6. Ejericio 24
Falta de memoria Con los datos del problema anterior suponga que la razón con la cual
el material es olvidado es proporcional a la cantidad memorizada al tiempo t. Determine una
ecuación diferencial para A(t), cuando se considera la falta de memoria
Datos
M = Cantidad de tema a memorizar
A(t)= Cantidad de memoria en tiempo t
dA
= rapidez con la que se memoriza
dt
A=cantidad del tema que falta por memorizar
K= constante de proporcionalidad
dA
= K(M − A) : Teoria del apredizaje
dt
Desarrollo
dA
= K1 (M − A)−k 2 A
dt
Conclusión: Tenemos como resultado final una ecuacion diferencial de primer orden

1.7. Ejericio 27
Superficie reflectora Suponga que cuando la curva plana C que se muestra en la figura
1.3.21 se gira respecto al eje x genera una superficie de revolución, con la propiedad de que
todos los rayos de luz paralelos al eje x que inciden en la superficie son reflejados a un solo

7
punto O (el origen). Utilice el hecho de que el ángulo de incidencia es igual al ángulo de reflexión
para determinar una ecuación diferencial que describa la forma de la curva C. Esta curva C
es importante en aplicaciones como construcción de telescopios o antenas de satélites, faros
delanteros de automóviles y colectores solares. [Sugerencia: La inspección de la figura muestra
que podemos escribir φ = 2θ. ¿Por qué? Ahora utilice una identidad trigonométrica adecuada.]

Datos:
φ = 2θ
Desarrollo:

2θ + α = π (2)
y
tan(α) = (3)
−x
0
y = tan(θ) (4)
Reemplazamos (2) en (3)
y
tan(π − 2θ) =
−x
y
tan(π) − tan(2θ) =
−x
y
0 − tan(2θ) =
−x
y 2tan(θ)
=
x
1 − tan2 (θ)
x2tan(θ)
y= (5)
1 − tan2 (θ)
Reemplazamos (4) en (5)
2xy 0
y= 2
1 − y0
2
y − yy 0 = 2xy 0
2
yy 0 + 2xy 0 − y
Formula General del trinomio cuadrado perfecto
+
√
−b 2
− b − 4ac
y0 =
2a
Encontramos los valores de y 0 1

8
 p
−x + 4x2 − 4y 2
y01 =
2y
Encontramos los valores de y 0 2
p
0 −x − 4x2 − 4y 2
y2=
2y
Conclusión: Tenemos como resultado final una ecuacion diferencial de primer orden

1.8. Ejericio 30
Lea nuevamente el análisis que condujo a la ecuación (8). Si suponemos que inicialmente el
tanque conserva, digamos 50 libras de sal, es porque se le está agregando sal continuamente al
tanque para t > 0, A(t) será una función creciente. Analice cómo podrı́a determinar a partir de
la ED, sin realmente resolverla, el número de libras de sal en el tanque después de un periodo
largo.
Resolución:

Datos: Desarrollo
dA A
Ce = 2lb/gal =6−
dt 100
A
V0 = 300gal A(t) = 6 −
100
600 − A
Cs =? 0=
100
A0 = 50lb A = 600
Fe = 3gal/min
Fs = 3gal/min
t>0

Conclusión Si A(t) continúa aumentando sin alcanzar un máximo, entonces A0 (t) > 0 ,para
t > 0 y A no puede exceder de 600.En este caso, si A0 (t) se acerca a 0 cuando t aumenta hasta
el infinito, A(t) se acerca a 600 cuando t aumenta hasta el infinito

1.9. Ejericio 33
Cuerpo en caı́da En el problema 21 suponga que r = R + s donde s es la distancia desde la
superficie de la Tierra al cuerpo que cae. ¿Cómo es la ecuación diferencial que se obtuvo en el
problema 21 cuando s es muy pequeña en comparación con R? Sugerencia: Considere la serie
binomial para: (R + s)2 = R−2 (1 + s/R)−2
Resolución:
Datos:
d2 r −gR2
= (6)
dt2 r2
R = constante

Desarrollo

9
 r =R+S (7)
Encontramos la segunda derrivada
d2 r d2 s
= (8)
dt2 dt2
Remplazamos (6) en (8)
d2 s −gR2
= (9)
dt2 r2
Remplazamos (7) en (9)
d2 s −gR2
=
dt2 (R + s)2
d2 s
dt2
= −gR2 ∗ (R + s)2
Por la serie de Taylor se obtiene
2
ds
≈ −gR2 + [R−2 − 2sR−3 + .....]
dt2
d2 s
2
≈ −g + 2gs
R3
+ ...
dt
Conclusión

d2 s
Por lo tanto, como R mucho mayor que s, la ecuación diferencial se aproxima a = −g
dt2

1.10. Ejericio 36
Lea nuevamente esta sección y clasifique cada modelo matemático como lineal o no lineal
dp
dt
= kp Lineal

dA
dt
= kA Lineal

dT
dt
= k(T − Tm ) Lineal

dx
dt
= kx(n + 1 − x) No Lineal

dX
dt
= k(α − x)(β − x) No Lineal

dA A
dt
=6− 100
Lineal
√
dh
dt
= − AhAw2gh No Lineal
2
1 ddt2q + R dq
dt
+ 1c q = E(t) Lineal

d2 s
dt2
= −g Lineal

m dv
dt
= mg − kv Lineal
2
m ddt2s + k ds
dt
= mg Lineal

10
2. Resolver múltiplos de 3 a partir del ejercicio 7 del
libro de Dennis Zill Cap 3.1
2.1. Ejercicio 9
Cuando pasa un rayo vertical de luz por un medio transparente, la razón con que decrece su
intensidad I es proporcional a I(t), donde t representa el espesor, en pies, del medio. En agua
limpia de mar, la intensidad a 3pies debajo de la superficie es 25º/º de la intensidad inicial I0
del rayo incidente. ¿Cuál es la intensidad del rayo a 15 pies debajo de la superficie?

Datos:
t0 = 0 I = I0
t1 = 3pies I1 = 0,25
t2 = 15pies I2 =?
k= Constante
Pregunta
¿Cuál es la intensidad del rayo a 15 pies debajo de la superficie?
Desarrollo:
dI
= kI (10)
dt
Encontramos los valores de I(0) e I(3)
I(0) = I0
I(3) = 0,25I
Integramos (10)
R dI R
= kdt
I
ln(I) = kt + c
I = cekt
I(t) = cekt
Encontramos el valor de c
I0 = cek0
I0 = c
Encontramos el valor de k
I(3) = cekt
0,25I0 = I0 ek3
0,25 = ek3
3k = ln(0,25)
k = −0,462

11
 Encontramos el valor de I(15)
I(t) = I0 e−0,462t
I(15) = I0 e−0,462∗15
I(15)=0.00098 º/º I0

Conclisión: Porlo tanto tenemos que la intencidad del rayo de luz a 15 pies es del 0.00098 º/º
de la intencidad inicial I0

2.2. Ejercicio 12
El sudario de Turı́n muestra el negativo de la imagen del cuerpo de un hombre que parece
que fue crucifi cado, muchas personas creen que es el sudario del entierro de Jesús de Nazaret.
Vea la fi gura 3.1.9. En 1988 el Vaticano concedió permiso para fechar con carbono el sudario.
Tres laboratorios cientı́fi cos independientes analizaron el paño y concluyeron que el sudario
tenı́a una antigüedad de 660 años,* una antigüedad consistente con su aparición histórica.
Usando esta antigüedad determine qué porcentaje de la cantidad original de C-14 quedaba en
el paño en 1988.

Datos:
A0 = 100 t=0
1
A1 = A
2 0
t = 560
A2 =? t = 660
k= Constante
Pregunta
¿Qué porcentaje de la cantidad original de C-14 quedaba en el paño en 1988?
Desarrollo:
dA
= kA (11)
dt
Integramos (11)
R dI R
= kdt
A

12
 ln(A) = kt + c
A = cekt
I(t) = cekt
Encontramos el valor de c
A0 = cek0
A0 = c
Encontramos el valor de k
1
A0 = A0 ekt
2
1
= ek∗5600
2
5600k = ln(0,5)
k = −0,000124
Encontramos el valor de A(660)
A(t) = A0 e0,000124t
A(660) = 100e0,000124∗660
A(660) =92º/º

Conclusión: La muestra que se tomó del sudario tubo un 92 por ciento de carbono 14 en 1988
cuando fue estudiado.

2.3. Ejercicio 15
Una pequeña barra de metal, cuya temperatura inicial era de 20° C, se deja caer en un gran
tanque de agua hirviendo. ¿Cuánto tiempo tardará la barra en alcanzar los 90° C si se sabe que
su temperatura aumentó 2° en 1 segundo? ¿Cuánto tiempo tardará en alcanzar los 98° C?

Datos:
T0 = 20C T m = 100C
T1 = 90C t =?
T3 = 98C t =?
I = 2C t = 1s
K=constante
Preguntas
¿Cuánto tiempo tardará la barra en alcanzar los 90° C si se sabe que su temperatura aumentó
2° en 1 segundo?
¿Cuánto tiempo tardará en alcanzar los 98° C?
Desarrollo
dT
= k(T − T m)
dt
dT
= k(T − 100C) (12)
dt

13
 R dT R
Integramos (12) = kdt
(T − 100C)

ln(T − 100C) = kt + c
T − 100C = cekt T − 100C = cekt
T = cekt + 100C
T (t) = cekt + 100C
Encontramos el valor de c
T0 = cek0 + 100C
20 = cek0 + 100C
−80 = c
T (t) = −80ekt + 100C
Encontramos el valor de k
T1 = −80ek∗1 + 100C
22C = −80ek∗1 + 100C
−78C = −80ek∗1
78
= ek v
80
k = −0,0253
T (t) = −80e−0,0253∗t + 100C
Encontramos el tiempo de T2
T2 = −80e−0,0253∗t + 100C
90 = −80e−0,0253∗t + 100Cv
−10 = −80e−0,023∗5t
−10
= e−0,0253∗t
−80
t = 82,19s
Encontramos el tiempo de T3
T3 = −80e−0,0253∗t + 100C
98 = −80e−0,0253∗t + 100C
−2 = −80e−0,0253∗t
−2
= e−0,025∗t
−80
t = 146s

Conclusión: Para que el objeto llegue a los 90 grados, incrementando su temperatura 2 grados
cada segundo necesita 82.19 segundos y por consiguiente lo mismo pasa para que el objeto
llegue alos 98 grados, igual incrementando 2 grados cada segundo necesita 146 segundos

14
2.4. Ejercicio 18
Al tiempo t = 0 un tubo de ensayo sellado que contiene una sustancia quı́mica está inmerso
en un baño lı́quido. La temperatura inicial de la sustancia quı́mica en el tubo de ensayo es de
80° F. El baño lı́quido tiene una temperatura controlada (medida en grados Fahrenheit) dada
por Tm (t) = 100 − 4e−0,1t , t ≥ 0 donde t se mide en minutos.

a). Suponga que k = −0,1 en la ecuación (2). Antes de resolver el PVI, describa con palabras
cómo espera que sea la temperatura T (t) de la sustancia quı́mica a corto plazo. Y a largo plazo.

b).Resuelva el problema con valores iniciales. Use un programa de graficación para trazar la
gráfica de T (t) en diferentes intervalos de tiempo. ¿Las gráficas concuerdan con sus predicciones
del inciso a?
Datos:
T0 =80ºF t=0
Tm (t) = 100 − 4e−0,1t t≥0
k = −0,1
Desarrollo
Literal a)
Tm (t) = 100 − 4e−0,1t
Tm (0) = 100 − 4e−0,1∗0
Tm = 60ºF
Tenemos que la temperatura ambiente es de 60ºF y la del quimico es de 80ºF, a corto plazo la
temperatura del quimico disminuira ya que se adapta ala temperatura ambiente, mientras que
a larrgo plazo la temperatura debera incrementarse hasta llegar a los 100ºF, el cual es el punto
de ebullición.
Literal b)
dT
= k(T − 100 + 4e−kt )
dt
dT
= −0,1(T − 100 − 4e−0,1t )
dt
dT
= −0,1T + 10 + 4e−0,1t )
dt
dT
+ 0,1T = 10 + 4e−0,1t ) Edo lineal
dt
dT
+0,1T = 10 − 4e−0,1t
dt P (t) Q(t)

Buscamos el Factor Integrante

R
P (x)dt
FI = e
R
0,1dt
FI = e
F I = e0,1t
Resolvemos la integral
R
F I ∗ T = F IQ(t)

15
 e0,1t (10 − 4e−0,1t )
R
e0,1t T =
R
e0,1t T = 10e0,1t − 4
e0,1t T = 100e0,1t − 4t + c
T (t) = 100 − 4te−0,1t + ce−0,1t
Buscamos el valor de C
T (0) = 100 − 4 ∗ 0 ∗ e−0,1∗0 + ce−0,1∗0
80 = 100 + c
c = −20
T (t) = 100 − 4te−0,1t − 20e−0,1t
T (t) = 100 − e−0,1t (4t + 20)

Figura 2: Aumento de temperatura

Conclusión: Tenemos que la temperatura del quimico siempre llegara a uan temperatura final
de 100ºF

2.5. Ejercicio 21
Un tanque contiene 200 litros de un lı́quido en el que se han disuelto 30 g de sal. Salmuera
que tiene 1 g de sal por litro entra al tanque con una razón de 4 L/min; la solución bien
mezclada sale del tanque con la misma razón. Encuentre la cantidad A(t) de gramos de sal que
hay en el tanque al tiempo t.
Resolución:
Datos:
V0 = 200l
A0 = 30g
Ce = 1 Lg L
Fe = 4 min
L
Cs =? Fs = 4 min
Pregunta:
dA
=?
dt
Desarrollo
dA
= Ce F e − Cs F s
dt
dA L

g
L

A


= 4m 1L − 4m 200
dt
16
 dA
= 4 mg − 50 A



dt
Separamos las variables
dA
  = dt
A
(4) −
50
Resolvemos la integral
R dA R
  = dt
A
(4) −
50
  
A
−50ln (4) − =t+c
50
  
A
ln (4) − = − 50t + c
50
 
A t
(4) − = ce −50
50
 
A t
− = ce −50 − 4
50
t
A = 200 + ce −50
Encontramos el valor de c
0
A(0) = 200 + ce −50
30 = 200 + c
c = −170
t
A = 200 − 170e −50
Conclusión: obtenemos la ecucion que nos permitira calcular la concentracion con la que sale
nuestra salmuera en un tiempo t

2.6. Ejercicio 24
En el problema 23, ¿cuál es la concentración t) de sal en el tanque al tiempo t? ¿Y al
tiempo t = 5 min? ¿Cuál es la concentración en el tanque después de un largo tiempo, es decir,
conforme t → ∞? ¿Para qué tiempo la concentración de sal en el tanque es igual a la mitad de
este valor lı́mite?

Datos:
Del Ejercicio 23 tenemos:
A
C(t) =
500
t
A(t) = 1000 − 1000e− 100
Pregunta:
C(t) =?
C(5) =?

17
 C(∞) =?
Desarrollo
t
A(t) = 1000 − 1000e− 100
 t

− 100
A(t) = 500 2 − 2e
A(t) t
= 2 − 2e− 100
500
A
Remplazamos por C(t)
500
t
C(t) = 2 − 2e− 100
Encontramos el valor de C(5) =?
t
C(5) = 2 − 2e− 100
5
C(5) = 2 − 2e− 100
lb
C(5) = 0,0975 gal
Encontramos el valor de C(∞)
 t

C(t) = lim 2 − 2e− 100
t→∞
 t

− 100
C(t) = lim (2) − lim 2e
t→∞ t→∞

C(t) = 2 − 0
C(t) = 2
La mitad de C(∞) =1 por lo tanto
t
−
1 = 2 − 2e 100
t
−
−1
−2
= e 100
t
ln(0,5) = −
100
t = 69,31min
lb
Conclusión: Por lo tanto para un tiempo de 5 minutos la concentración es de 0,0975 gal y el
tiempo que se demora la concentración de sal en el tanque es de 69,31min

2.7. Ejercicio 27
Un gran tanque está parcialmente lleno con 100 galones de fluido en los que se disolvieron
1
10 libras de sal. La salmuera tiene de sal por galón que entra al tanque a razón de 6 gal/min.
2
La solución bien mezclada sale del tanque a razón de 4 gal/min. Determine la cantidad de
libras de sal que hay en el tanque después de 30 minutos.
Datos:
Volumen = 100 galones
A0 = 10 libras; t0 =0 min
1
centra = lb/min
2
18
fentra = 6 gal/min
fsale = 4 gal/min

Pregunta:

Determine la cantidad de libras de sal que hay en el tanque después de 30 minutos.

Desarrollo:

dA
= rentra − rsale
dt
dA
= centra fentra − csale fsale
dt
dA 1
= (6) − csale (4)
dt 2
dA A A
=3− (4) (Reemplazo, c = )
dt 100 + (6 − 4)t V
dA 4A
=3−
dt 100 + 2t
dA 4A
=3−
dt 2(50 + t)
dA 2A
=3−
dt 50 + t
dA 2A
+ =3
dt 50 + t
0
EDO lineal de primer orden de la forma y (x) + p(x)y = q(x) en función de t

2
p(t) =
50 + t
q(t) = 3

Determinar el Factor integrante (F.I):

R
p(t)dx
u(t) = e
R 2
dt
u(t) = e 50 + t Integración por sustitución
2ln(50+t)
u(t) = e
u(t) = (50 + t)2 Factor integrante

19
Establecer la ecuación:
Z
u(t)A = (u(t)q(t))dt
Z
(50 + t) A = ((50 + t)2 )3dt
2

Z
(50 + t) A = 3 (2500 + 100t + t2 )dt
2

t3
 
2 2
(50 + t) A = 3 2500t + 50t +
3
2 2 3
(50 + t) A = 7500t + 150t + t + C
7500t + 150t2 + t3 + C
A= Solución general
(50 + t)2

Reemplazo las condiciones iniciales:

A0 = 10t0 =0
7500(0) + 150(0)2 + (0)3 + C
10 =
(50 + 0)2
C
10 =
2500
C = 25000

Por lo tanto,

7500t + 150t2 + t3 + 25000

A= Solución particular
(50 + t)2

Calcúlo la cantidad de libras de sal que hay en el tanque despúes de 30 min:

t = 30
7500(30) + 150(30)2 + (30)3 + 25000
A=
(50 + 30)2
A = 64,375

Conclusión:

La cantidad de libras de sal que habrá en el tanque después de 30 minutos es 64.375 lb.

2.8. Ejercicio 30
Resuelva la ecuación (7) suponiendo que E(t) = E0 sen(ωt) y que i(0) = i0 .

Datos:
di
L + Ri = E(t) (7)
dt

20
E(t) = E0 Senωt

i(0) = i0

Desarrollo

Dada la ecuación (7) reemplazamos el valor de E(t) en la ecuación:

di
L + Ri = E(t)
dt
di
L + Ri = E0 Senωt
dt
E0 R E0 Lω
i= 2 2 2
Senωt − 2 2 Cosωt + ce−Rt/L
L ω +R L ω + R2
Donde:

i(0) = i0

En conclusión obtenemos:
E0 Lω
c = i0 +
L ω 2 + R2
2

2.9. Ejercicio 33
Se aplica una fuerza electromotriz:
(
12, 0 ≤ t ≤ 20
E(t) =
0, t > 20
a un circuito en serie LR en el que la inductancia es de 20 henrys y la resistencia es de 2 ohm.
Determine la corriente i(t), si i(0)=0.
Datos:
L=20 h
R=20 ohm
v= ? i(t), si i(0)=0

Desarrollo:
Ley de kirchhoff→ i(t)
di
L+ + Ri = E(t)
dt
Remplazomos los valores que tenemos en ecuación diferencial lineal
di
20 + 2i = 120
dt
i = 60

21
Aplicamos factor integral
t t
 
d 
e 10  = 6e 10
dt

t
−
i = 60 + c1 e 10
evaluamos la condición inicial i(0) = 0 → c1 = −60
t
−
i = 60 − 60e 10

tomamos encuenta la segunda condición t > 20 y aplicamos la Ecuación diferencial:

di
20 + 2i = 0
dt
aplicar factor integrante
t
 
d 
e 10  = 0
dt

t
−
i = c2 e 10
tomamos t=20
t t
− −
60 − 60e 10 = ce 10
c2 e−2 = 60 − 60e−2
c2 = 60(e2 − 1)
Respuesta:  t

 −
 60 − 60e 10 , 0 ≤ t ≤ 20
i(t) = t
−


60(e2 − 1)e 10 , t > 20


2.10. Ejercicio 36
¿Qué tan alto? (Sin resistencia del aire) Suponga que una pequeña bala de cañón que pesa
16 libras se dispara verticalmente hacia arriba, como se muestra en la figura 3.1.10, con una
pies
velocidad inicial de v0 = 300 . La respuesta a la pregunta “¿Qué tanto sube la bala de
s
cañón?”, depende de si se considera la resistencia del aire.
a) Suponga que se desprecia la resistencia del aire. Si la dirección es positiva hacia arriba,
d2 s
entonces un modelo para la bala del cañón está dado 2 (ecuación (12) de la sección 1.3).
dt
ds dv
Puesto que = v(t) la última ecuación diferencial es la misma que la ecuación = −g donde
dt dt
pies
se toma g = 32 2 . Encuentre la velocidad v(t) de la bala de cañón al tiempo t.
s

22
Datos:

w = 16lb
pies
vo = 300
s

Desarrollo:
Formula inicial

d2 s
dt2
= −mg
Formula simplificada

d2 s
= −g
dt2

(1)Conociendoque :
d2 s dv
2
= =a
dt dt

d2 s dv
Reemplazamos: dt2
= dt
en(1)

dv
= −g
dt

Separamos las variables e integramos

Z Z
dv = −g dt
Z Z
dv = −g dt
v(t) = −gt

23
 ds ds
(2)Reemplazamos : porv = −gt + c
dt dt
Reemplzamosnuestrodatov(0) = 300
pies
s
300 = −g(0) + c
c = 300

Dato 1
g = 32
Reemplazamos nuestros datos en la ecuación (2)

v(t) = −32t + 300 (3)

b) Utilice el resultado que se obtuvo en el inciso a) para determinar la altura s(t) de la bala de
cañón medida desde el nivel del suelo. Determine la altura máxima que alcanza la bala.
Utilizamos nuestra ecuación (3) v(t) = −32t + 300

ds
= −32t + 300
dt
Sabemos que nuestra altura máxima se da cuando la velocidad es igual a 0.

∴ 0 = −32t + 300
Ası́ que calculamos nuestro tiempo.

32t = 300
300
t=
32
t = 9,375s

Ahora integramos nuestra ecuación de la velocidad para obtener una ecuación de posición.

ds
= −32t + 300
Z Zdt
ds = (−32t + 300)dt
32 2
s=− t + 300t
2

Reemplazamos nuestro tiempo de 9,375s.

s = −16(9,375)2 + 300(9,375)
s = 104,25f t

24
2.11. Ejercicio 39
Evaporación de una gota de lluvia Cuando cae una gota de lluvia, ésta se evapora
mientras conserva su forma esférica. Si se hacen suposiciones adicionales de que la rapidez a
la que se evapora la gota de lluvia es proporcional a su área superficial y que se desprecia la
resistencia del aire, entonces un modelo para la velocidad v(t) de la gota de lluvia es
k
3( )
dv ρ
+ v=g
dt k
( )t + r0
ρ
Aquı́ ρ es la densidad del agua, r0 es el radio de la gota de lluvia en t = 0, k < 0 es la constante
de proporcionalidad y la dirección hacia abajo se considera positiva.

a) Determine v(t) si la gota de lluvia cae a partir del reposo.

b) Vuelva a leer el problema 34 de los ejercicios 1.3 y demuestre que el radio de la gota de
lluvia en el tiempo t es r(t) = (k/ρ)t + r0 .

c) Si r0 = 0,01 pies y r = 0,007 pies, 10 segundos después de que la gota cae desde una
nube, determine el tiempo en el que la gota de lluvia se ha evaporado por completo.

Datos
k
3( )
dv ρ
Ecuación Deferencial lineal: + v=g
dt k
( )t + r0
ρ
v0 = 0 t0 = 0

k<0 Constante de proporcionalidad.

r0 = 0,01 pies t0 = 0 seg.

r = 0,007 pies t=10 seg.

Desarrollo
Literal(a):
k
3( )
dv ρ
+ v=g
dt k
( )t + r0
ρ
Sacamos el Factor Integrante:

3(k/ρ)
R
dt
u(t) = e (k/ρ)t + r0

R (k/ρ)
3 dt
u(t) = e (k/ρ)t + r0

25
u(t) = e3ln[(k/ρ)t+r0 ]

u(t) = [(k/ρ)t + r0 ]3
Empezamos a integral la EDO:
R
[(k/ρ)t + r0 ]3 v(t) = [(k/ρ)t + r0 ]3 gdt
R
[(k/ρ)t + r0 ]3 v(t) = g [(k/ρ)t + r0 ]3 dt
ρR k
[(k/ρ)t + r0 ]3 v(t) = g [(k/ρ)t + r0 ]3 dt
k ρ
ρ
[(k/ρ)t + r0 ]3 v(t) = g [(k/ρ)t + r0 ]4 + c
4k
v0 = 0 t0 = 0
ρ
[(k/ρ)(0) + r0 ]3 (0) = g [(k/ρ)(0) + r0 ]4 + c
4k
ρ
0=g (r0 )4 + c
4k
ρ
c = −g (r0 )4
4k
ρ ρ
[(k/ρ)t + r0 ]3 v(t) = g [(k/ρ)t + r0 ]4 − g (r0 )4
4k 4k
Despues de resolverla obtenemos:

gρ gρ(r0 )4
v(t) = [(k/ρ)t + r0 ] − g
4k 4k[(k/ρ)t + r0 ]3
Literal(b):
m 4
ρ= v = πr3
v 3
m=ρv

dm d(ρv)
=
dt dt
4
m = ρ( πr3 )
3
dm
= 4kπr2
dt
4
d(ρ( πr2 ))
3 = 4kπr2
dt
1
ρdr
3 =k
dt
3k
dr = dt
ρ
k
dr = dt
ρ

26
R R k
dr = dt
ρ
k
r(t) = t + c
ρ
Si t = 0 entonces tenemos que:
k
r0 = (0) + c
ρ
c = r0

Reemplazando el valor de C nos queda:

k
r(t) = t + r0
ρ
Literal(c):

r0 = 0,01 pies t0 = 0 seg.

r = 0,007 pies t = 10 seg.

k
r(t) = t + r0
ρ
k
r(10) = (10) + 0,01
ρ
k
0,007 = (10) + 0,01
ρ
k
(10) = −0,003
ρ
k
= −0,0003
ρ
Se Evapora por completo cuando r(t) = 0

0 = −0,0003t + 0,01

0,0003t = 0,01

t = 33,33 seg

2.12. Ejercicio 42
Modelo de cosecha constante Un modelo que describe la población de una pesquerı́a en
la que se cosecha con una razón constante está dada por
dp
= kp − h
dt
donde k y h son constantes positivas.

27
 a) Resuelva la ED sujeta a P (0) = P0 .

b) Describa el comportamiento de la población P(t) conforme pasa el tiempo en los tres

casos P0 > h/k, P0 = h/k y 0 < P0 < h/k.

c) Utilice los resultados del inciso b) para determinar si la población de peces desaparecerá
en un tiempo finito, es decir, si existe un tiempo T > 0 tal que P(T)=0. Si la población
desaparecerá, entonces determine en qué tiempo T
Datos
Ecuación Deferencial lineal.
dp
= kp − h
dt
k y h Constantes positivos.

P (0) = P0 t=0

P0 > h/k, P0 = h/k y 0 < P0 < h/k

T >0 P(T)=0

Desarrollo
a)
dp
= kp − h
dt
dp
− kp = −h
dt
Sacamos el Factor Integrante.
R
−kdt
u(t) = e
u(t) = e−kt
Integramos la E.D Lineal.
e−kt P = e−kt (−h)dt
R

e−kt P = −h e−kt dt
R

hR
e−kt P = −ke−kt dt
k
h
e−kt P = e−kt + c
k
P (0) = P0 t=0
h −k(0)
e−k(0) P0 = e +c
k
h
P0 = + c
k
h
c = P0 −
k
h h
e−kt P (t) = e−kt + P0 −
k k
P (t) nos queda:

28
 h h
P (t) = + (P0 − )ekt
k k
b)
P0 > h/k Esto será positivo ya que reemplazando en la ecuación obtenida, como P0 es mayor
que h/k, saldrá positivo.
Para: P0 = h/k
h h h
P (t) = + ( − )ekt
k k k
h
P (t) =
k
Para 0 < P0 < h/k En este caso la función se irá disminuyendo a medida que el tiempo t va
aumentando.

c)
Como sabemos que función disminuye. Se puede decir que existe un tiempo t > 0 tal que
exista P(t)=0
h h
0 = + (P0 − )ekt
k k
h h
− = (P0 − )ekt
k k
h
−
kt
e = k
h
P0 −
k
h
−
ekt = k
kP0 − h
k
hk
ekt = − 2
k P0 − kh
hk
ekt =
kh − k 2 P0
hk
ekt =
k(h − kP0 )
h
ekt =
h − kP0
h
kt = ln( )
h − kP0
1 h
t = ln( )
k h − kP0

2.13. Ejercicio 45
Marcapasos de corazón. En la figura se muestra un marcapasos de corazón, que consiste
en un interruptor, una baterı́a, un capacitor y el corazón como resistor. Cuando el interruptor
S está en P , el capacitor se carga; cuando S está en Q, el capacitor se descarga, enviando
estı́mulos electricos al corazón.En el problema 47 de los ejercicios 2.3 vimos que durante este
tiempo en que se están aplicando estı́mulos eléctricos al corazón, el voltaje E a través del

29
corazón satisface la ED lineal
dE 1
=− E.
dt RC
a) Suponga que en el intervalo de tiempo de duración t1 , 0 < t < t1 , el interruptor S esta en la
posición P como se muestra en la figura y el capacitor se está cargando. Cuando el interruptor
se mueve a la posición Q al tiempo t1 el capacitor se descarga, enviando un impulso al corazón
durante el intervalo de tiempo de duración t2 : t1 ≤ t < t1 + t2 . Por lo que el intervalo inicial
de carga descarga 0t < t1 + t2 el voltaje en el corazón se modela realmente por la ecuación
diferencial definida por tramos.

dE 0, 0 ≤ t < t1
= 1
dt − RC E, t1 ≤ t < t1 + t2

Al moverse S entre P y Q, los intervalos de carga y descarga de duraciones t1 y t2 se repiten

indefinidamente. Suponga que t1 = 4s, t2 = 2s, E0 = 12v, E(0) = 0, E(4) = 12, E(6) = 0,
E(10) = 10, E(12) = 0, etc. Determine E(t) para 0 ≤ t ≤ 24.
b) Suponga para ilustrar que R = C = 1. Utilice un programa de graficación para trazar la
gráfica de la solución del PVI del inciso a) para 0 ≤ t ≤ 24.

Datos:
Cuando el interruptor S está en P , el capacitor se carga.
Cuando S está en Q, el capacitor se descarga.

Pregunta:
Determine E(t) para 0 ≤ t ≤ 24.
Utilice un programa de graficación para trazar la gráfica de la solución del PVI del inciso
a) para 0 ≤ t ≤ 24.
Desarrollo:

a) Escribimos la ecuación diferencial en el intervalo de (0,24):

1

dE − RC E , t ∈ [4, 6i ∪ [10, 12i ∪ [16, 18i ∪ [22, 24i
=
dt 0 , t ∈ [0, 4i ∪ [6, 10i ∪ [12, 16i ∪ [18, 22i
(1)
dE
dt
=0

30
Con la solución: E = c
y algunas condiciones iniciales:

E(0) = E(6) = E(12) = E(18) = 0

Por lo tanto, E = 0 en:

[0, 4i ∪ [6, 10i ∪ [12, 16i ∪ [18, 22i

En [4, 6i ∪ [10, 12i ∪ [16, 18i ∪ [22, 24i tenemos:

dE 1
=− E
dt RC
1 1
dE = − dt
Z E Z RC
1 1
dE = − dt
E RC
1
ln E = − t+c
RC
1
E = de− RC t

Reemplazo los valores iniciales:

4 4
[4, 6] E(4) = 12 =⇒ 12 = de− RC =⇒ d = 12e RC
10 10
[10, 12] E(10) = 12 =⇒ 12 = de− RC =⇒ d = 12e RC
16 16
[16, 18] E(16) = 12 =⇒ 12 = de− RC =⇒ d = 12e RC
22 22
[22, 24] E(22) = 12 =⇒ 12 = de− RC =⇒ d = 12e RC

Por lo tanto, la ecuación diferencial en el intervalo []0, 24 es:

 4−t

 12e RC , t ∈ [4, 6i
 10−t
 12e πC , t ∈ [10, 12i


16−t
E= 12e RC , t ∈ [16, 18i
 22−t
 12e RC , t ∈ [22, 24i



0 , otros


b) La gráfica de la solución del PVI del inciso a) para 0 ≤ t ≤ 24 es:

2.14. Ejercicio 48
Qué sube . . . a) Es bien conocido que el modelo que desprecia la resistencia del aire,
inciso a) del problema 36, predice que el tiempo ta que tarda la bala de cañón en alcanzar su
altura máxima es el mismo tiempo td que tarda la bala de cañón en llegar al suelo. Además la
magnitud de la velocidad de impacto vi es igual a la velocidad inicial v0 de la bala de cañón.
Compruebe ambos resultados.

31
 Figura 3: gráfica de la solución del PVI del inciso a) para 0 ≤ t ≤ 24

b) Después, utilizando el modelo del problema 37 que considera la resistencia del aire, compare
el valor de ta con td y el valor de la magnitud de vi con v0 . Aquı́ puede ser útil un programa
para determinar raı́ces con un SAC (o una calculadora graficadora).

Datos

Del ejercicio 36 el tiempo de ascenso es: ta = 9,375

v = 300

s(t) = −16t2 + 300t

v(t) = −32t + 300

Del ejercicio 37 el tiempo de ascenso es: ta = 9,162

s(t) = 1340000 − 6400t − 1340000e−0,005t

v(t) = −6400 + 6700e−0,005t

Desarrollo
Literal (a)

Usaremos el tiempo de ascenso del ejercicio 36 ta = 9,375, para encontrar el tiempo de descenso
con el que la bala de cañon llega al suelo, para eso resolveremos primero la siguiente ecuación:
S(t) = −16t2 + 300t = 0, de donde obtenemos un tiempo total de vuelo de: ttv = 18,75 s

Partiendo de aquı́ tendremos que el tiempo de descenso será:

td = ttv − ta

td = 18,75 − 9,375

td = 9,375

Ahora bien la velocidad de impacto de la bala de cañon será:

vi = v

32
vi = −300, por ir en descenso

Literal (b)

Usaremos el tiempo de ascenso del ejercicio 37, el cual tiene resistencia del aire ta = 9,375, para
encontrar el tiempo de descenso con el que la bala de cañon llega al suelo, para eso resolveremos
primero la siguiente ecuación:
s(t) = 1340000 − 6400t − 1340000e−0,005t = 0

Resolviendo el sistema con la ayuda de un programa se obtiene que el tiempo es:

ttv = 0,1846 s

Partiendo de aquı́ tendremos que el tiempo de descenso será:

td = ttv − ta

td = 0,1846 − 9,162

td = −8,9774 s

Ahora bien la velocidad de impacto de la bala de cañon será:

vi = v

vi = −291,00226, por ir en descenso

En conclusión: Los tiempos td y la velocidad vi , tanto del literal (a) y (b), no varian mucho

3. REFERENCIAS
Referencias
[1] D.G.Zill;Differential Equations with Boundary-Vary Problems CAP 1.3 y CAP 3.1.
BROOKS/COLE. 2009.

33