Taller 2° Seguimiento 2021 II.
Taller 2° Seguimiento 2021 II.
Taller 2° Seguimiento 2021 II.
FACULTAD DE INGENIERÍA
ECUACIONES DIFERENCIALES
PROF: DEUD SOTO PALOMINO
TALLER SEGUNDO SEGUIMIENTO
R//
𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ + 2𝑦 = 0 ; 𝑦1 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥)
𝑒 − ∫ 𝑃(𝑥)𝑑𝑥
𝑦2 = 𝑦1 ∫ 𝑑𝑥
(𝑦1 )2
(𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ + 2𝑦 = 0)
𝑥2
1 2𝑦
𝑦 ′′ − 𝑦 ′ + 2 = 0
𝑥 𝑥
1
𝑒 − ∫ −𝑥𝑑𝑥
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ 2 𝑑𝑥
( 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥))
𝑒 𝑙𝑛𝑥
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ 𝑑𝑥
𝑥 2 𝑠𝑒𝑛 2 (𝑙𝑛𝑥)
𝑥
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ 𝑑𝑥
𝑥 2 𝑠𝑒𝑛 2 (𝑙𝑛𝑥)
𝑥
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ 𝑑𝑥
𝑥 2 𝑠𝑒𝑛 2 𝑢
𝑢 = 𝑙𝑛𝑥
1
𝑑𝑢 = 𝑑𝑥
𝑥
𝑑𝑥 = 𝑥 𝑑𝑢
𝑥
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ 𝑥 𝑑𝑢
𝑥2 𝑠𝑒𝑛 2 𝑢
𝑥2
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ 𝑑𝑢
𝑥 2 𝑠𝑒𝑛 2 𝑢
1
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ 𝑑𝑢
𝑠𝑒𝑛2 𝑢
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) ∫ csc 2 𝑢 𝑑𝑢
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) (− cot 𝑢)
cos 𝑢
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) (− )
𝑠𝑒𝑛 𝑢
cos ln 𝑥
𝑦2 = 𝑥 𝑠𝑒𝑛(𝑙𝑛𝑥) (− )
𝑠𝑒𝑛 ln 𝑥
𝑦2 =−𝑥 cos ln 𝑥
𝑦𝐺 = 𝑐1 𝑦1 + 𝑐2 𝑦2
𝑦𝐺 = 𝑐1 𝑥 cos ln 𝑥 + 𝑐2 − 𝑥 cos ln 𝑥
R//
𝐷2 𝑦 + 𝑦 = 𝑥𝑒 −𝑥 + 3 𝑠𝑒𝑛 𝑥
Resolviendo por el método de coeficiente indeterminados
𝐷2 𝑦 + 𝑦 = 0
Ecuación característica
𝑓(𝑥) = 𝑥𝑒 −𝑥 + 3 𝑠𝑒𝑛 𝑥
𝑥𝑒 −𝑥 + 3 𝑠𝑒𝑛 𝑥 → (𝑎𝑥 + 𝑏)𝑒 −𝑥 + 𝑑𝑐𝑜𝑠𝑥 + 𝑒𝑠𝑒𝑛𝑥
Verificamos que el 𝑦𝑝 este correcto
Igualamos a la ecuación:
−2𝑎 𝑒 −𝑥 + (𝑎𝑥 + 𝑏)𝑒 −𝑥 − 2𝑑𝑠𝑒𝑛𝑥 + 2𝑒𝑐𝑜𝑠𝑥 − 𝑥𝑒𝑠𝑒𝑛𝑥 − 𝑥𝑑𝑐𝑜𝑠𝑥 + (𝑎𝑥 + 𝑏)𝑒 −𝑥
+ 𝑑𝑥𝑐𝑜𝑠𝑥 + 𝑒𝑥𝑠𝑒𝑛𝑥 = 𝑥𝑒 −𝑥 + 3 𝑠𝑒𝑛 𝑥
−2𝑎 𝑒 −𝑥 + 2𝑎𝑥𝑒 −𝑥 + 2𝑏𝑒 −𝑥 − 2𝑑𝑠𝑒𝑛𝑥 + 2𝑒𝑐𝑜𝑠𝑥 = 𝑥𝑒 −𝑥 + 3 𝑠𝑒𝑛 𝑥
2𝑎𝑥𝑒 −𝑥 + (2𝑏 − 2𝑎)𝑒 −𝑥 − 2𝑑𝑠𝑒𝑛𝑥 + 2𝑒𝑐𝑜𝑠𝑥 = 𝑥𝑒 −𝑥 + 3 𝑠𝑒𝑛 𝑥
Tenemos el sistema de ecuaciones:
1
2𝑎 = 1 → 𝑎 =
2
1 1
2𝑏 − 2𝑎 = 0 → 𝑏 + 2𝑏 − 2 = 0 → 2𝑏 − 1 = 0 → 2𝑏 = 1 → 𝑏 =
2 2
3
−2𝑑 = 3 →= −
2
2𝑒 = 0 → 𝑒 = 0
𝑦𝑝 = (𝑎𝑥 + 𝑏)𝑒 −𝑥 + 𝑥(𝑑𝑐𝑜𝑠𝑥 + 𝑒𝑠𝑒𝑛𝑥)
1 1 3
𝑦𝑝 = ( 𝑥 + )𝑒 −𝑥 + 𝑥(− 𝑐𝑜𝑠𝑥 + 0𝑠𝑒𝑛𝑥)
2 2 2
𝑒 −𝑥 3
𝑦𝑝 = (𝑥 + 1) − 𝑥𝑐𝑜𝑠𝑥
2 2
1
𝑦𝑝 = ((𝑥 + 1)𝑒 −𝑥 − 3𝑥𝑐𝑜𝑠𝑥) Solución particular
2
𝑦𝐺 = 𝑦ℎ + 𝑦𝑝
1
𝑦𝐺 = 𝑐1 𝑐𝑜𝑠 𝑥 + 𝑐2 𝑠𝑒𝑛 𝑥 + ((𝑥 + 1)𝑒−𝑥 − 3𝑥𝑐𝑜𝑠𝑥) Solución General
2
R//
𝑦 = 𝑦ℎ + 𝑦𝑝
(𝐷2 + 3𝐷 + 2)𝑦 = 0
𝐷2 + 3𝐷 + 2 = 0
−3 ± √32 − 4 ∗ 1 ∗ 2
𝐷1,2 =
2∗1
−3 ± 1
𝐷1,2 =
2∗1
−3 + 1 −3 + 1
𝐷1 = 𝐷2 =
2∗1 2∗1
𝐷 = −1 𝐷 = −2
𝑦 = 𝑐1 𝑒 −𝑥 + 𝑐2 𝑒 −2𝑥
𝑔(𝑥) = 𝑠𝑒𝑛(𝑒 𝑥 ) −→ 𝑦𝑝 = 𝑢1 𝑦1 ∗ 𝑢2 𝑦2
𝑦2 𝑔(𝑥)
𝑢1 = ∫ − 𝑑𝑥
𝑊(𝑦1 , 𝑦2 )
𝑦1 𝑔(𝑥)
𝑢2 = ∫ − 𝑑𝑥
𝑊(𝑦1 , 𝑦2 )
𝐷𝑜𝑛𝑑𝑒 𝑒𝑙 𝑤𝑟𝑜𝑛𝑠𝑘𝑖𝑎𝑛𝑜 𝑊(𝑦1 , 𝑦2 ) = 𝑦1 𝑦2 − 𝑦1 𝑦2
𝑦1 = 𝑒 −𝑥
𝑦2 = 𝑒 −2𝑥
𝑦´1 = 𝑒 −𝑥
𝑦´2 = −2𝑒 −2𝑥
𝑊(𝑦1 , 𝑦2 ) = 𝑦1 𝑦´2 − 𝑦´1 𝑦2
𝑊(𝑦1 , 𝑦2 ) = −𝑒−3𝑥
𝑦2 𝑔(𝑥)
𝑢1 = ∫ − 𝑑𝑥
𝑊(𝑦1 , 𝑦2 )
𝑒 −2𝑥 𝑠𝑒𝑛 (𝑒 𝑥 )
𝑢1 = ∫ − 𝑑𝑥
−𝑒 −3𝑥
𝑢1 = − cos(𝑒 𝑥 )
𝑦1 𝑔(𝑥)
𝑢2 = ∫ − 𝑑𝑥
𝑊(𝑦1 , 𝑦2 )
𝑒 −𝑥 𝑠𝑒𝑛 (𝑒 𝑥 )
𝑢2 = ∫ − 𝑑𝑥
−𝑒 −3𝑥
𝑒 −𝑥 𝑠𝑒𝑛 (𝑒 𝑥 )
∫− 𝑑𝑥 = 𝑒 𝑥 cos 𝑒 𝑥 − sen 𝑒 𝑥 + 𝐶
−𝑒 −3𝑥
𝑢2 = 𝑒 𝑥 cos 𝑒 𝑥 − sen 𝑒 𝑥
𝑦𝑝= 𝑢1 𝑦1 + 𝑢2 𝑦2
𝑦𝑝 = (cos(𝑒𝑥 )) 𝑒−𝑥 + (𝑒𝑥 cos(𝑒𝑥 ) − sen(𝑒𝑥 )) 𝑒−2𝑥
R//
𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ + 2𝑦 = cos(ln 𝑥 3 )
aplicando la sustitución
𝑧 = 𝑙𝑛𝑥 ⟶ 𝑒 𝑧 = 𝑥
ln 𝑥 3 = 3 ln 𝑥 → 3 ln 𝑥 = 3𝑧
𝑥 2 𝑦 ′′ = Dz(Dz − 1)y
𝑥𝑦 ′ = 𝐷𝑍 𝑦
Dz 2 𝑦 − 𝐷𝑍 y + 𝐷𝑍 𝑦 + 9𝑦 = cos(3𝑧)
Dz 2 𝑦 + 9𝑦 = cos(3𝑧)
Buscamos el 𝑦ℎ
Dz 2 𝑦 + 9𝑦 = 0
𝑚2 + 9 = 0
𝑚2 = −9
√𝑚2 = √−9
𝑚 = ±3𝑖
𝛼=0 𝛽=3
𝑦ℎ = 𝑐1 𝑠𝑒𝑛3𝑧 + 𝑐2 𝑐𝑜𝑠3𝑧
Ahora buscamos el 𝑦𝑝
𝑓(𝑥) → 𝑦𝑝
𝑦𝑝 = 𝑧(𝑎𝑠𝑒𝑛(3𝑧) + 𝑏𝑐𝑜𝑠(3𝑧))
Sustituimos en la ecuación:
Dz 2 𝑦 + 9𝑦 = cos(3𝑧)
1
6𝑎𝑐𝑜𝑠(3𝑧) = cos (3𝑧) ⟶ 6𝑎 = 1 → 𝑎 =
6
−6𝑏𝑠𝑒𝑛(3𝑧) = 0 ⟶ 𝑏 = 0
𝑦𝑝 = 𝑧(𝑎𝑠𝑒𝑛(3𝑧) + 𝑏𝑐𝑜𝑠(3𝑧))
1
𝑦𝑝 = 𝑧( 𝑠𝑒𝑛(3𝑧) + (0)𝑐𝑜𝑠(3𝑧))
6
1
𝑦𝑝 = 𝑧 𝑠𝑒𝑛(3𝑧)
6
La solución general es:
𝑦𝐺 = 𝑦ℎ + 𝑦𝑝
1
𝑦𝐺 = 𝑐1 𝑠𝑒𝑛3𝑧 + 𝑐2 𝑐𝑜𝑠3𝑧 + 𝑧 𝑠𝑒𝑛(3𝑧)
6
Volviendo a la sustitución original:
1
𝑦𝐺 = 𝑐1 𝑠𝑒𝑛 ln 𝑥3 + 𝑐2 𝑐𝑜𝑠 ln 𝑥3 + ln 𝑥 𝑠𝑒𝑛(ln 𝑥3 )
6
5. Muestre que un peso 𝑊, dada una velocidad inicial 𝑣0 , se desliza una distancia 𝑠 hacia abajo
por un plano inclinado sin fricción de inclinación 𝛼 en el tiempo
R//
y
y
x
α
α
α x
w
w
wy
Para estudiar este caso se utiliza la componente en y del peso ya que horizontalmente no se está
ejerciendo ninguna fuerza contraria.
𝑦 = 𝑊𝑠𝑒𝑛𝑎 → 𝑦 = 𝑚𝑔 𝑠𝑒𝑛 𝑎
Aplicando la segunda ley de newton
𝑑𝑣
Σ𝐹 = 𝑎𝑚 → Σ𝐹 = 𝑚
𝑑𝑡
𝑑𝑣 𝑑𝑣
𝑚 = 𝑚𝑔 𝑠𝑒𝑛 𝑎 → = 𝑔 𝑠𝑒𝑛 𝑎 → 𝑑𝑣 = (𝑔 𝑠𝑒𝑛 𝑎)𝑑𝑡
𝑑𝑡 𝑑𝑡
∫ 𝑑𝑣 = ∫ 𝑔 𝑠𝑒𝑛 𝑎 𝑑𝑡 → ∫ 𝑑𝑣 = 𝑔 𝑠𝑒𝑛 𝑎 ∫ 𝑑𝑡
𝑣 = 𝑔𝑡 𝑠𝑒𝑛 𝑎 + 𝑐
𝑣0 = 𝑔(0) 𝑠𝑒𝑛 𝑎 + 𝑐
𝑣0 = 𝑐
𝑣 = 𝑔𝑡 𝑠𝑒𝑛 𝑎 + 𝑣0
∫ 𝑑𝑠 = 𝑔 𝑠𝑒𝑛 𝑎 ∫ 𝑡 𝑑𝑡 + 𝑣0 ∫ 𝑑𝑡
𝑡2
𝑠 = 𝑔 𝑠𝑒𝑛 𝑎 + 𝑣0 𝑡 + 𝑐
2
02
0 = 𝑔 𝑠𝑒𝑛 𝑎 + 𝑣0 (0) + 𝑐
2
𝑐=0
Entonces:
𝑡2
𝑠 = 𝑔 𝑠𝑒𝑛 𝑎 + 𝑣0 𝑡
2
Despejando el tiempo 𝑡:
2𝑠 = 𝑔 𝑠𝑒𝑛 𝑎𝑡 2 + 2𝑣0 𝑡
0 = 𝑔 𝑠𝑒𝑛 𝑎𝑡 2 + 2𝑣0 𝑡 − 2𝑠
√𝑣0 2 + 2𝑔 𝑠𝑒𝑛 𝑎𝑠 − 𝑣0
𝑡=
𝑔 𝑠𝑒𝑛 𝑎
6. En cada punto 𝑃(𝑥, 𝑦) de una curva, el ángulo formado por la tangente y la ordenada es
bisecado por la recta que une al punto con el origen. Encuentre la ecuación de la curva sabiendo
que pasa por el punto (1,2).
R//
Hay que resolver los ángulos de los triángulos para poder determinar la ecuación de la recta
tangente y de esta forma al aplicar la ecuación diferencial obtener la familia de curvas. Esto quiere
𝑑𝑦
decir que la ecuación diferencial será la tangente del ángulo de la recta tangente 𝑑𝑥 = tan 𝛽
𝜃 + 𝛼 + 90 = 180
𝜃 + 𝛼 = 90
𝛼 = 90 − 𝜃
los ángulos del triángulo 1 tenemos:
𝛽 + 2𝛼 + 90 = 180
𝛽 + 2𝛼 = 90
2𝛼 = 90 − 𝛽
Ahora 𝛼 = 90 − 𝜃
2(90 − 𝜃) = 90 − 𝛽
180 − 2𝜃 = 90 − 𝛽
90 − 2𝜃 = −𝛽
2𝜃 − 90 = 𝛽
Con esta expresión podemos despejar los catetos del triángulo 1 en función del ángulo 𝛽:
tan 𝛽 = tan 2𝜃 − 90
𝜋
tan 2𝜃 − 90 → tan 2𝜃 −
2
𝜋
𝜋 sen 2𝜃 −
tan 2𝜃 − = 2
2 cos 2𝜃 − 𝜋
2
𝜋 𝜋 𝜋
sen 2𝜃 − 2 sen 2𝜃 cos 2 − sen 2 cos 2𝜃 sen 2𝜃 (0) − (1) cos 2𝜃 − cos 2𝜃
𝜋= 𝜋 𝜋 = cos 2𝜃 (0) + sen 2𝜃(1) = sen 2𝜃
cos 2𝜃 − cos 2𝜃 cos + sen 2𝜃 sen
2 2 2
− cos 2𝜃
tan 𝛽 =
sen 2𝜃
−(2 cos2 𝜃 − 1)
tan 𝛽 =
2 sen 𝜃 cos 𝜃
1 − 2 cos2 𝜃
tan 𝛽 =
2 sen 𝜃 cos 𝜃
Recordando la imagen del triángulo 1, vemos a que catetos se puede igualar las identidades
trigonométricas de seno y coseno
2
𝑥
1− 2( )
√𝑥 2 + 𝑦 2
tan 𝛽 =
𝑦 𝑥
2( )( )
√𝑥 2
+𝑦 2 √𝑥 + 𝑦 2
2
2𝑥 2 𝑥 2 + 𝑦 2 − 2𝑥 2
1−
𝑥2 + 𝑦 2 𝑥2 + 𝑦 2 (𝑥 2 + 𝑦 2 )(𝑥 2 + 𝑦 2 − 2𝑥 2 ) 𝑥 2 + 𝑦 2 − 2𝑥 2 𝑦 2 − 𝑥 2
= = = =
2𝑦𝑥 2𝑦𝑥 (2𝑦𝑥)(𝑥 2 + 𝑦 2 ) 2𝑦𝑥 2𝑦𝑥
𝑥2 + 𝑦 2 𝑥2 + 𝑦 2
𝑑𝑦 𝑦 2 − 𝑥 2
=
𝑑𝑥 2𝑦𝑥
𝑦 2 𝑥2
𝑑𝑦 𝑥 2 − 𝑥 2
= 𝑦𝑥
𝑑𝑥 2 2
𝑥
𝑦2
𝑑𝑦 𝑥 2 − 1
= 𝑦
𝑑𝑥 2𝑥
𝑦 𝑑𝑦 𝑑𝑧
𝑆𝑒𝑎 𝑧 = → 𝑧𝑥 = 𝑦 → =𝑧+𝑥
𝑥 𝑑𝑥 𝑑𝑥
𝑑𝑧 𝑧 2 − 1
𝑧+𝑥 =
𝑑𝑥 2𝑧
𝑑𝑧 𝑧 2 − 1
𝑥 = −𝑧
𝑑𝑥 2𝑧
𝑑𝑧 𝑧 2 − 1 − 2𝑧 2
𝑥 =
𝑑𝑥 2𝑧
𝑑𝑧 −1 − 𝑧 2
𝑥 =
𝑑𝑥 2𝑧
2𝑧 𝑑𝑥
∫ 𝑑𝑧 = ∫
−1 − 𝑧 2 𝑥
7. Determine la ecuación de la curva que pasa por el punto (0,1) la cual es ortogonal a cada
miembro de la familia 𝑥 2 + 𝑦 2 = 𝑐𝑒 𝑥 .
R//
La ecuación de la curva que es ortogonal a la familia de curvas y pasa por el punto 0,1 es la función
que tiene como pendiente la recta normal a la recta tangente a la familia de curvas
𝑥 2 + 𝑦 2 = 𝑐𝑒 𝑥
𝑥2 + 𝑦 2
=𝑐
𝑒𝑥
Derivando implícitamente:
(2𝑥 + 2𝑦𝑦´)𝑒 𝑥 − (𝑥 2 + 𝑦 2 )𝑒 𝑥
=0
𝑒 2𝑥
𝑒 𝑥 (2𝑥 + 2𝑦𝑦´ − 𝑥 2 − 𝑦 2 )
=0
𝑒 2𝑥
(2𝑥 + 2𝑦𝑦´ − 𝑥 2 − 𝑦 2 )
=0
𝑒𝑥
2𝑥 + 2𝑦𝑦´ − 𝑥 2 − 𝑦 2 = 0
2𝑦𝑦´ = 𝑥 2 − 2𝑥 + 𝑦 2
𝑥 2 − 2𝑥 + 𝑦 2
𝑦´ =
2𝑦
Con esto ya se tiene la pendiente de la función que es ortogonal a la familia de curvas, entonces
𝑑𝑥 2𝑦
=− 2
𝑑𝑦 𝑥 − 2𝑥 + 𝑦 2
𝑑𝑦 2𝑦
=− 2
𝑑𝑥 𝑥 − 2𝑥 + 𝑦 2
𝑦
𝑧= 𝑧𝑥 = 𝑦 𝑑𝑥 = 𝑦𝑑𝑧 + 𝑧𝑑𝑦
𝑥
(𝑥 2 − 2𝑥 + 𝑦 2 )𝑑𝑦
=0
−2𝑦𝑑𝑥
( 𝑧 2 𝑦 2 − 2𝑧𝑦 + 𝑦 2 )𝑑𝑦 = −2𝑦(𝑦𝑑𝑧 + 𝑧 𝑑𝑦)
𝑧 2 𝑦 2 𝑑𝑦 − 2𝑧𝑦𝑑𝑦 + 𝑦 2 𝑑𝑦 = −2𝑦𝑧𝑑𝑦 − 2𝑦 2 𝑑𝑧
𝑧 2 𝑦 2 𝑑𝑦 + 𝑦 2 𝑑𝑦 = −2𝑦 2 𝑑𝑧
(𝑧 2 + 1)𝑦 2 𝑑𝑦 = −2𝑦 2 𝑑𝑧
𝑦 2 𝑑𝑦 −2
2 = 2 𝑑𝑧
𝑦 𝑧 +1
−2
𝑑𝑦 = 𝑑𝑧
𝑧2+1
−2
∫ 𝑑𝑦 = ∫ 𝑑𝑧
𝑧2+1
𝑑𝑧
∫ 𝑑𝑦 = −2 ∫
𝑧2+1
𝑦 = −2 𝑡𝑎𝑛−1 𝑧 + 𝑐1
𝑥
𝑦 = −2 𝑡𝑎𝑛−1 ( ) + 𝑐1
𝑦
0
1 = −2 𝑡𝑎𝑛−1 ( ) + 𝑐1
1
1 = 𝑐1
Entonces la curva que es ortogonal a la familia de curvas 𝑥 2 + 𝑦 2 = 𝑐𝑒 𝑥 y que pasa por el punto
0,1 es:
𝑥
𝑦 = −2 𝑡𝑎𝑛−1 ( ) + 1
𝑦
𝐶𝐸0 𝑡
𝑞= (cos 𝜔𝑡 + 𝜔𝑅𝐶 𝑠𝑒𝑛 𝜔𝑡 − 𝑒 −𝑅𝐶 )
𝑅 𝐶 2𝜔2
2 +1
R//
𝑉 = 𝐸0 cos 𝜔𝑡
tenemos que el voltaje es 𝑉 = 𝐸0 cos 𝜔𝑡 y que es la suma de resistencia más el capacitor:
𝑑𝑞 𝑞
𝑅 + = 𝐸0 cos 𝜔𝑡
𝑑𝑡 𝐶
𝑑𝑞 𝑞
𝑅( + ) = 𝐸0 cos 𝜔𝑡
𝑑𝑡 𝐶𝑅
𝑑𝑞 𝑞 𝐸0 cos 𝜔𝑡
+ =
𝑑𝑡 𝐶𝑅 𝑅
1 𝐸0 cos 𝜔𝑡
La E.D. de Kirchhoff con 𝑃(𝑡) = 𝐶𝑅 y 𝑄(𝑡) = . Entonces el factor integrante es:
𝑅
1 1 𝑡
∫𝐶𝑅𝑑𝑡 ∫ 𝑑𝑡
𝑒 ∫ 𝑃(𝑡)𝑑𝑡 = 𝑒 =𝑒 𝐶𝑅 =𝑒 𝐶𝑅 . Multiplicamos la E.D. dada por el factor integrante
Tenemos que:
𝑞𝑢 = ∫ 𝑄(𝑡)𝑢𝑑𝑡
Reemplazando:
𝑡 𝐸0 cos 𝜔𝑡 𝑡
𝑞𝑒 𝐶𝑅 = ∫ 𝑒 𝐶𝑅 𝑑𝑡
𝑅
𝐸0 cos 𝜔𝑡 𝑡 𝐸0 𝑡
∫ 𝑒 𝐶𝑅 𝑑𝑡 → ∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡
𝑅 𝑅
𝑡
∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡
𝑢 = cos 𝜔𝑡 𝑑𝑢 = −𝜔𝑠𝑒𝑛𝜔𝑡
𝑡
𝑡 𝑒 𝐶𝑅
𝑑𝑣 = 𝑒 𝐶𝑅 𝑣=
𝐶𝑅
𝑡
𝑡 cos 𝜔𝑡 𝑒 𝐶𝑅 𝜔 𝑡
∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡 = + ∫ 𝑒 𝐶𝑅 𝑠𝑒𝑛𝜔𝑡
𝐶𝑅 𝐶𝑅
𝑢 = 𝑠𝑒𝑛𝜔𝑡 𝑑𝑢 = 𝜔𝑐𝑜𝑠𝜔𝑡
𝑡 𝑡
𝑑𝑣 = 𝑒 𝐶𝑅 𝑣 = 𝐶𝑅𝑒 𝐶𝑅
𝑡 𝑡 𝑡 𝑡
∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡 = 𝐶𝑅cos 𝜔𝑡 𝑒 𝐶𝑅 + 𝐶𝑅𝜔 (𝐶𝑅𝑠𝑒𝑛𝜔𝑡𝑒 𝐶𝑅 − 𝐶𝑅𝜔 ∫ 𝑒 𝐶𝑅 𝑐𝑜𝑠𝜔𝑡𝑑𝑡)
𝑡 𝑡 𝑡 𝑡
∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡 = 𝐶𝑅cos 𝜔𝑡 𝑒 𝐶𝑅 + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔𝑡𝑒 𝐶𝑅 − (𝐶𝑅)2 𝜔2 ∫ 𝑒 𝐶𝑅 𝑐𝑜𝑠𝜔𝑡𝑑𝑡
𝑡 𝑡 𝑡 𝑡
∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡 + (𝐶𝑅)2 𝜔2 ∫ 𝑒 𝐶𝑅 𝑐𝑜𝑠𝜔𝑡𝑑𝑡 = 𝐶𝑅cos 𝜔𝑡 𝑒 𝐶𝑅 + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔𝑡𝑒 𝐶𝑅
𝑡 𝑡 𝑡
∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡 (1 + (𝐶𝑅)2 𝜔2 ) = 𝐶𝑅cos 𝜔𝑡 𝑒 𝐶𝑅 + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔𝑡𝑒 𝐶𝑅
𝑡
𝑡 𝑒 𝐶𝑅 (𝐶𝑅cos 𝜔𝑡 + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔𝑡)
∫ cos 𝜔𝑡 𝑒 𝐶𝑅 𝑑𝑡 = + 𝑐1
(1 + (𝐶𝑅)2 𝜔 2 )
Tenemos que:
𝑡 𝐸0 𝐸0 cos 𝜔𝑡 𝑡
𝑞𝑒 𝐶𝑅 = (∫ 𝑒 𝐶𝑅 𝑑𝑡)
𝑅 𝑅
𝑡
𝑡 𝐸0 𝑒 𝐶𝑅 (𝐶𝑅cos 𝜔𝑡 + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔𝑡)
𝑞𝑒 𝐶𝑅 = ( ) + 𝑐1
𝑅 (1 + (𝐶𝑅)2 𝜔 2 )
𝑡
𝑡 𝐸0 𝑒 𝐶𝑅
𝑞𝑒 𝐶𝑅 = ( 2 2 2 (𝐶𝑅cos 𝜔𝑡 + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔𝑡)) + 𝑐1
𝑅 𝑅 𝐶 𝜔 +1
Si 𝑞 = 0 en 𝑡 = 0
𝑡
𝑡 𝐸0 𝑒 𝐶𝑅
𝑞𝑒 𝐶𝑅 = ( 2 2 2 (𝐶𝑅cos 𝜔𝑡 + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔𝑡)) + 𝑐1
𝑅 𝑅 𝐶 𝜔 +1
(0)
(0) 𝐸0 𝑒 𝐶𝑅
(0)𝑒 𝐶𝑅 = ( 2 2 2 (𝐶𝑅cos 𝜔(0) + (𝐶𝑅)2 𝜔𝑠𝑒𝑛𝜔(0))) + 𝑐1
𝑅 𝑅 𝐶 𝜔 +1
𝐸0 𝐶𝑅
0= + 𝑐1
𝑅 𝑅 𝐶 𝜔2 + 1
2 2
𝐸0 𝐶
𝑐1 = −
𝑅 𝐶2𝜔2
2 +1
𝑡
𝑡 𝑒 𝐶𝑅 𝐸 𝐶
𝑞𝑒 𝐶𝑅 = 𝐸0 ( 2 2 2 (𝐶cos 𝜔𝑡 + 𝐶 2 𝑅𝜔𝑠𝑒𝑛𝜔𝑡)) − 2 2 0 2
𝑅 𝐶 𝜔 +1 𝑅 𝐶 𝜔 +1
−𝑡
𝐸0 𝐶 𝐸 𝐶𝑒 𝐶𝑅
𝑞=( 2 2 2 (cos 𝜔𝑡 + 𝐶𝑅𝜔𝑠𝑒𝑛𝜔𝑡)) − 2 02 2
𝑅 𝐶 𝜔 +1 𝑅 𝐶 𝜔 +1
𝐸0 𝐶 −𝑡
𝑞= (cos 𝜔𝑡 + 𝐶𝑅𝜔𝑠𝑒𝑛𝜔𝑡 − 𝑒 𝐶𝑅 )
𝑅 𝐶 2𝜔2
2 +1
9. Un cuerpo se calienta a 110 °C y se expone al aire libre a una temperatura de 10 °C. Si al cabo
de una hora su temperatura es de 60 °C. ¿Cuánto tiempo adicional debe transcurrir para que
se enfríe a 30 °C?
R// Datos:
𝑡 = 0 →→ 𝑇 = 110°
𝑇𝑚 = 10°
𝑑𝑇
= −𝑘(𝑇 − 𝑇𝑚)
𝑑𝑡
𝑑𝑇
∫ = −𝑘 ∫ 𝑑𝑡
(𝑇 − 𝑇𝑚)
𝑢 = 𝑇 − 𝑇𝑚
𝑑𝑢 = 𝑑𝑡
𝑑𝑢
∫ = −𝑘𝑡
𝑢
𝑒 ln (𝑇−𝑇𝑚) = 𝑒 (−𝑘𝑡+𝐶)
𝑇 − 𝑇𝑚 = 𝑒 −𝑘𝑡 ∗ 𝑒 𝑐 →→ 𝐶
𝑇 − 𝑇𝑚 = 𝐶𝑒 −𝑘𝑡
𝑇 = 𝐶𝑒 −𝑘𝑡 + 𝑇𝑚
Reemplazamos:
110 = 𝐶𝑒 0 + 10
100 = 𝐶
𝑇 = 100 𝑒 −𝑘𝑡 + 𝑇𝑚
𝑡 = 1 ℎ −→ 𝑇 = 60
60 = 100𝑒 −𝑘 + 10
50 = 100 𝑒 −𝑘
50 1
= 𝑒 −𝑘 − −→ = 𝑒 −𝑘
100 2
1
ln ( ) = 𝑙𝑛 𝑒 −𝑘
2
ln(1) − ln(2) = −k
− ln(2) = k −→ k = ln (2)
𝑇 = 100 𝑒 − ln(2)𝑡 + 𝑇𝑚
𝑇 = 30
Reemplazamos:
30 = 100 𝑒 − ln(2)𝑡 + 10
20 = 100 𝑒 − ln(2)𝑡
20
= 𝑒 − ln(2)𝑡
100
1
= 𝑒 − ln(2)𝑡
5
1
ln ( ) = 𝑙𝑛 𝑒 − ln(2)𝑡
5
− ln(5) = − ln(2) 𝑡
ln(5)
𝑡=
ln(2)
𝑡 = 2,3 ℎ
10. Si en un análisis de una botella de leche se encuentran 500 organismos (bacterias), un día
después de haber sido embotelladas y al segundo día se encuentran 8000 organismos. ¿Cuál
es el número de organismos en el momento de embotellar la leche?
R//
𝑡1 = 1 𝑑𝑖𝑎 𝑝1 = 500
𝑡2 = 2 𝑑𝑖𝑎𝑠 𝑝2 = 8000
𝑑𝑝
= 𝑎𝑝
𝑑𝑡
𝑑𝑝
= 𝑎𝑑𝑡
𝑝
𝑑𝑝
∫ = ∫ 𝑎𝑑𝑡
𝑝
𝑑𝑝
∫ = 𝑎 ∫ 𝑑𝑡
𝑝
ln 𝑝 = 𝑎𝑡 + 𝑐
𝑒 ln 𝑝 = 𝑒 𝑎𝑡+𝑐
𝑝 = 𝑒 𝑎𝑡 𝑒 𝑐
𝑝 = 𝑒 𝑎𝑡 𝑘 𝑘 ∴ 𝑒𝑐
500 = 𝑒 𝑎 𝑘
500
=𝑘
𝑒𝑎
Ahora, tenemos que 𝑡2 = 2 𝑑𝑖𝑎 𝑝2 = 8000
8000
=𝑘
𝑒 2𝑎
8000 500
= 𝑎
𝑒 2𝑎 𝑒
500𝑒 2𝑎
8000 =
𝑒𝑎
8000 = 500𝑒 −𝑎 𝑒 2𝑎
8000
= 𝑒𝑎
500
8000
ln | | = ln 𝑒 𝑎
500
8000
ln | |=𝑎
500
500
=𝑘
8000
500
125
=𝑘
4
Tenemos la población en función del tiempo
8000
125𝑒 ln| 500 |𝑡
𝑝(𝑡) =
4
Igualando 𝑡 = 0 tenemos la población de bacterias en el momento de embotellar la leche
8000
125𝑒 ln| 500 |(0)
𝑝(0) =
4
125
𝑝(0) =
4