TALLER DE APLICACIÓN DE ECUACIONES DIFERENCIALES

1. Una taza de café se enfría según la ley de Newton para el enfriamiento (3).
Use los datos d la gráfica de temperaturas T(t) ilustrada en la figura.
Las constantes Tm T0 y K en un modelo del tipo de un problema de valor inicial
de primer orden

dT
=K ( T −T m )
dt
T ( 0 )=T 0

dT
∫ T −T =∫ kdt
m

ln ( T −T M ) =kt +C

T −T m=e Kt+C

T =T m +C e Kt (1)

Para T ( 0 )=T 0

T 0=T m +C e 0
T 0=T m +C

C=T 0−T m

Reemplazando C en la ecuación (1)

T =T m + ( T 0−T m ) e Kt

Para una T m=80

T =80+ ( T 0−80 ) e Kt

Para T ( 0 )=160

160=80+ ( T 0−80 ) e 0

160=80+T 0−80

T 0=160

Reemplazando T0 en la ecuación (1)

T =80+ 80 e Kt
Para T (100) =90

90=80+80 e 100 k

10=80 e100 k
10 100k
=e
80
1
ln =lne 100 k
8
1
ln =100 k
8
1
ln ( )
8
K=
100
K=−0,02079441542
Reemplazando K en la ecuación (1)

T =80+ 80 e−0,02079441542t
2. La población de una comunidad aumenta a una tasa que es proporcional al
número de personas presente en el tiempo t. Si una población inicial P 0 se ha
duplicado en 5 años ¿cuánto tardará en triplicarse?, ¿y en cuadruplicarse?

Formula general

p= p 0=∁ ℮kt ⟹ p=∁℮kt

p= p 0
t 0=0
t ( s ) =2 p0 ; se va aduplicar en 5 años

p0=∁ ℮ k (0 )
p0=∁ ℮0

p0=∁

Para t (5) = 2p0

2 p 0=c ℮k (5)

2 p 0= p0 ℮ 5k

2=℮5 k

ln 2=ln ℮5 k

ln 2=5 k

ln 2
k= =0,138
5

¿Cuánto tiempo tardará en triplicarse?

3 p0=∁ ℮(0,138) t
3 p0= p0 ℮(0,138 ) t
3=℮( 0,138) t

ln 3=ln ℮(0,138) t

ln 3=0,138 t

ln3
t= =7,9 ≃8 años .
0,138

¿Cuánto tiempo tardará en cuadriplicarse?

4 p 0=∁ ℮(0,138) t

4 p 0= p0 ℮(0,138) t

4=℮( 0,138) t

ln 4=ln ℮(0,138) t

ln 4=0,138 t

ln 4
t= =10,04 ≃ 10 años
0,138

3. un tanque contiene 200 litros de fluido en el cual se han disuelto 30

gramos de sal. La salmuera, que contiene un gramo de sal por litro se
bombea hacia el deposito a una velocidad de 4L/min: perfectamente
mezclada. La solución se bombea hacia fuera a la misma velocidad.
Encuentre el número A(t) de gramos de sal presentes en el tanque en el
tiempo t.

La velocidad con la que cambia la cantidad de sal es A´(t) en el tanque

de mezcla.
A(t) = cantidad de agua en el tanque en el instante t.
V(t) = Volumen de fluido en el tanque en el instante t.
t = tiempo en minutos

Formula
 dA
=razón de entrada de sal−razón de salida de sal
dt

De otra manera se puede ver como

dA
=R entrada−R salida
dt

Donde R es concentración por salida (velocidad)

Datos

 V ( t ) =200l
 La cantidad inicial de sal disuelta A0 =A 0=30 g
 Concentración de la solución entrante = 1g/L
 Velocidad de entrada = 4L/min
A
 Concentración de salida = ∙g/L
200
4L
 Velocidad de salida
min
Vamos hallar A(t) para eso definimos la velocidad de entrada y salida, entonces
1g 4 L 4 g
Rentrada = ∁℮ x V e = ∙ =
L min min

Ag 4 L A
R salida=∁ s x V s= ∙ = g /min
200 L min 50

dA
=R entrada−R salida
dt
Sustituyendo valores

dA g A g
=4 − ∙
dt min 50 min
dA 200−A
=
dt 50
200− A
dA= ∙ dt
50
50 dA=( 200−A ) dt
50
∙ dA=dt
200− A
 dA
50∫ = dt
200− A ∫
Hacemos un cambio de variable
u=200− A
d u=−dA
Sustituimos
du
−50∫ = dt
u ∫
−50 ln u=t+C
Reemplazamos a u
−50 ln (200− A)=t +C
−t
ln ( 200−A )= −C1
50
−t
−c 1
e ln (200− A )=e 50
−t
200−A=e 50 ∙ e−C 1 ; e−C =C 2
1

−t
− A=−200+C 2 e 50
−t
A=200−C 2 e 50
−t
50
A ( t )=200−C 2 e

Pero sabemos que la cantidad que ingresa de sal disuelta es A(0) = 30 gramos,
entonces.
−0
A ( 0 ) =200−C 2 e 50

30=200−C 2 e0
−170=−C2

C 2=170

Por lo tanto
 −t
A ( t )=200−170 e 50

4. Resolvemos el punto anterior teniendo como condición que se bombea agua

pura.
Para lo cual partimos de:

−t
A ( t )=200−C e 50
−t
A ( t )−200=−C e 50
−t
50
200−A ( t )=C e

200− A ( t ) −t
=e 50
C
−t
200−A (t )
ln ( ¿)=lne 50 ¿
C
200−A (t) −t
ln ( ¿)= ¿
C 50

200−A (t )
50 ln ( ¿)=t ¿
C

5. Una fuerza electromotriz

E ( t )= 120 ,0 ≤ t ≤ 20
{0 ,t >20
Se aplica a un circuito LR en serie que tiene inductancia de 20 henrys y
resistencia de 2 ohms, Encuentre la, corriente i(t) si i(0) = 0

Datos
R=2 Ω
L=20 H
i(0)=0
V R +V L =E ( t )

Sabemos que: V R=Ri (t)

 di(t)
V L=L
dt
di(t)
Ri ( t )+ L =E(t )
dt
Ri (t) di(t) E (t)
+ =
L dt L

2 ( t ) di ( t ) E ( t )
+ =
20 dt 20
di ( t ) i ( t ) E ( t )
+ =
dt 10 20

Para 0 ≤ t ≤ 20 entonces E (20) = 120

di ( t ) i ( t )
+ =120
dt 10

Resolvemos la ecuación diferencial.

Factor integrante

1 t
∫ p ( t ) dt ∫ 10 dt
10
e =e =e

t
di ( t ) 10t t
10 10
e + e i ( t )=120 e
dt
Reescribimos nuestra ecuación diferencial
p ( t ) dt
d ∫ p ( t) dt e∫ E (t )
dt
( e i ( t ) )=
L

t
t
d 10 e 10 120
dt
(
e i(t) = )
20
 t t

∫ (e 10
)
i(t) dt=6∫ e dt 10

t t
e 10 i (t )=6 e 10 ⋅10+ c
t −t t
10 10 10
i ( t )=60 e ∙ e +C ∙ e
t
10
i ( t )=60+ C ⋅e
Para i(0)
0
10
0=60+C ∙ e
0=60+C

C=−60 ( 1 )

Reescribiendo la ecuación
−t
i ( t )=60−60 e 10

Para t ≥ 20 ; E=0

di(t) R E (t)
+ i ( t )=
dt L L

di(t) 2
+ i ( t )=0
dt 20

di(t) 1
+ i ( t )=0
dt 10

di(t) −1
= i(t)
dt 10
 di ( t ) −1
∫ dt
=
10
∫ dt
−1
ln ( i ( t ) ) = t+C
10

−1
t
ln (i ( t ) )
e =e 10
+e C

−1
t
10
i ( t )=C e
Para t >20 ; t=20

i ( 20 )=C e−2
−1
t
10
i ( 20 )=60−60 e

i ( 20 )=60−60 e−2

Igualamos i (20) entonces tenemos

i ( 20 )=i ( 20 )

C e−2=60−60 e−2

60−60 e−2
C=
e−2

C=60 ( 1−e−2) ∙ e 2

C=60 ( e 2−e−2 ∙ e2 )

C=60 ( e 2−1 )
Sustituyendo tenemos:

−t
i ( t )=60 ( e2−1 ) e 10
Por lo tanto
 −1

{
t
10
i(t) 60−60 e 0 ≤t ≤ 20
−t
2 10
60 ( e −1 ) e t>20

6. La siguiente ecuación es una ecuación de Bernoulli la cual será resuelta por

sustitución.
dy 1
x + y= 2
dx y

dy
x + y= y −2 3
dx

w= y 1−n= y 1−(−2)= y 1 +2= y 3

w= y 3
1
3
y=w
2
dy 1 3 d w
= w
dx 3 dx
Ahora reemplazamos en la ecuación inicial
2
1 −1
1
−2
3 dw (3 )
x w + w 3 =w
3 dx

1 −2
−2 3 3
1 3 dw w w
w + =
3 dx x x

1 2 −2 2
d w 3 w 3 w3 3 w 3 w3
+ =
dx x x

d w 3 w 3 w0
+ =
dx x x

d w 3w 3
+ =
dx x x
Ahora tengamos en cuenta el método de factor de integración

∫ 3x dx x
3

e∫ p x dx =e =e 3 ln x =e ln =x3
( )

dw
x3 +3 x 2 w=3 x 2
dx

dw
x3 +3 x 2 w=3 x 2
dx
dw 3
( x w ) =3 x 2
dx

dw ( x 3 w ) =3 x 2 dx

∫ ( x 3 w ) dw=∫ 3 x 2 dx
x 3 w=x 3 +C

W =1+ C x−3 ; pero W = y3

−3

y 3=¿1 +C x ¿

3
y= √ 1+ x−3

7. La ecuación diferencial
dy
= p ( x ) +Q ( x ) y + R ( x ) y 2
dx

Se conoce como ecuación de Riccati.

a) Una ecuación de Riccati se puede resolver mediante la sucesión de dos

sustituciones siempre y cuando conozcamos una solución particular y 1
de la ecuación: Demuestre que la sustitución y = y 1 + u reduce la
ecuación de Riccati a una ecuación de Bernoulli (4) con n = 2. La
 ecuación de Bernoulli puede reducirse entonces a una ecuación lineal
mediante la sustitución de w = u-1

dy
= p ( x ) +Q ( x ) y + R ( x ) y 2 ( 1 )
dx

Bernoulli y ´ + g ( x )=f ( x ) y 2 ; para n=2

Sea y= y1 +u (2)

dy d y 1 d u
= +
dx dx dx

y ´= y ´ 1 +u ´ ( 3 )

Ahora reemplazamos la ecuación (3) y ecuación (2) en la ecuación (1)

y ´ 1 +u ´= p ( x )+ Q ( x ) ( y 1+ u ) + R(x )( y 1+u)2

y ´ 1 +u ´= p ( x )+ Q ( x ) y 1+Q ( X ) u+ y 21 R ( x ) +2 y 1 uR ( x ) + R ( x ) u2

u ´ =( p ( x ) +Q ( x ) y 1 + R ( x ) y 21 ) − y ´ 1+ ( Q ( x ) +2 R ( x ) y 1 ) u+ R ( x ) u2

u ´ = y ´ 1− y ´ 1+ ( Q ( x ) +2 R ( x ) y 1 ) u+ R ( x ) u2

u ´ =( Q ( x ) +2 R ( x ) y 1 ) u+ R ( x ) u2

R ( x ) u2=u´ −( Q ( x )+ 2 R ( x ) y 1 ) u( 4)

Ahora comparando la ecuación (4) con la ecuación de Bernoulli podemos

reducir a una ecuación diferencial lineal de primer grado Efectuando una
sustitución pertinente:

w=u−1
−1
d w d (u )
=
dx dx
 du
w ´ =−u−2
dx

w ´ =−u−2 u ´
Entonces podemos decir que: u ´ =−u−2 w ´ invertido

u ´ =−w2 w ´ (5)

para w−1=u

Ahora vamos a sustituir la ecuación (5) en la ecuación (4) de lo cual obtenemos

dw
−w−2 − ( Q ( x ) +2 R ( x ) y 1 ) w−1 =R ( x ) w−2
dx
Dividimos toda la expresión por −w−2 de lo cuál obtenemos
dw
−( Q ( x )+ 2 R ( x ) y 1 ) w=R ( x ) (6)
dx
La ecuación (6) es una ecuación lineal de primer orden.
b) Encuentre una familia de soluciones de un parámetro para la ecuación
diferencial.

dy −4 1
= − y + y 2 (1)
dx x 2 x

donde y1 = 2/x sea una solución conocida de la ecuación.

Partimos de la ecuación diferencial de Riccati

y ´= p ( x ) +Q ( x ) y+ R ( x ) y 2
Entonces:
−4
p ( x) =
x2
−1
Q ( x) =
x
R ( x )=1
Suponemos que:
y= y1 +u ( 2 )

Remplazamos la solución conocida

 2
y= +u ´
x

−2
y ´= +u ´ (3 )
x2

Reemplazamos las ecuaciones (3) y (2) en la ecuación (1)

2
u 2
−2 −4 x 2
x 2
+u ´ =
x 2
−
x ( )
+ +u
x

−2 2+ux 4 u 2
u´= − 2 + 2 4 +u
x2 x x x

−2 2 ux 4 u
u´= − 2 − 2 + 2 + 4 + u2
x 2
x x x x

−4 4
u´= +u +u2
x x
u 2
u ´ −3 =u ; ( 4 ) ecuación de Bernoulli para n=2
x

Sea w=u−1 ⇒ w ´=−u−2 u ´

Invertido u ´ =−w−2 w´ (5)

Para w-1 = u

Reemplazamos la ecuación (5) en la ecuación (4)

3
−w−2 w ´ − w−1=w−2
x
 3
w ´ − w−1 w 2=w−2 w2
x
3
w ´ + w=−1 ecuación diferencial de primer grado
x
3
p(x)=
x
Q ( x ) =−1
Ahora solucionamos por factor de integración

∫ P ( x ) dx ∫ 3x dx 3

e =e =e 3 ln x =e ln x =x 3

3
x 3 w ´− ∙ x 3 w=−x 3
x

x 3 w ´−3 x 2 w=−x3

∫ ( x 3 w ) ´ dw=−∫ x 3 dx
3 −x 4
x w= +C
4

−1 4 3
w= x ∙ x +C x−3
4
Pero w= u-2
−1
u−1 = x+C x 3
4
−1
−1
u= ( 4
x +C x−3 )
2
Pero u= y−
x
−1
2 −1
y− =
x 4 (
x +C x−3 )
 −1
2 −1
y= +
x (
4
x+ C x−3 )