SEGUNDA PARTE
DR. ING. MANUEL KUROKAWA GUERREROS
SOLUCIÓN DE E.D. POR SERIE DE TAYLOR:
Serie de Taylor: Sea f(x) infinitamente diferenciable en 𝐼 𝜖 𝑅.
Sea: 𝑥0 ∈ 𝐼
𝑓 (𝑘) (𝑥 )
֜ 𝑓 𝑥 = σ∞
𝑘=0
0
(𝑥 − 𝑥0 )𝑘 𝑓 (𝑘) (𝑥 ) : Derivada k ésima
𝑘! 0
Desarrollando se obtiene:
𝒇′ (𝒙 ) 𝒇′′ (𝒙 ) 𝒇′′′ (𝒙 )
𝒇 𝒙 = 𝒇(𝒙𝟎 ) + 𝟎
(𝒙 − 𝒙𝟎 )𝟏 + 𝟎
(𝒙 − 𝒙𝟎 )𝟐 + 𝟎
(𝒙 − 𝒙𝟎 )𝟑 + ⋯
𝟏! 𝟐! 𝟑!
En particular:
𝑓(𝑘) (0)
Para 𝑥0 = 0 ֜ 𝒇 𝒙 = σ∞
𝑘=0 𝑥𝑘 [Serie de Maclaurin]
𝑘!
Desarrollando:
𝒇′ 𝟎
𝒙 𝒇′′ 𝟎
𝒙𝟐 𝒇′′′ (𝟎) 𝒙𝟑
𝒇 𝒙 =𝒇 𝟎 + + + +⋯
𝟏! 𝟐! 𝟑!
PROBLEMA 1 :
Resolver 𝑦´´ + 2 𝑥 − 1 𝑦 ′ − 𝑦 = 0 , 𝒚 𝟏 =𝒚´ 𝟏 =𝟏
SOLUCIÓN:
Aplicaremos el método de Taylor, donde nos dan las condiciones iniciales
alrededor del punto x = 1.
𝑓(1) 𝑘
Sea la solución: y x = σ∞
𝑘=0 𝑘! 𝑥−1 𝑘
… (1)
De la ecuación diferencial tenemos:
𝑦´´ = 𝑦 − 2 𝑥 − 1 𝑦´ … … … … (2)
Para x=1: 𝑦´´ 1 = 𝑦 1 − 0 → 𝑦´´ 1 = 1
Derivando la ecuación 2: 𝑦´´´ = 𝑦´ − 2𝑦´ − 2 𝑥 − 1 𝑦´´ … … … (3)
Para x=1 𝑦´´´ 1 = −𝑦´ 1 − 2 𝑥 − 1 𝑦´´ = −1
Derivando la ecuación 3: 𝑦 𝑖𝑣 = −𝑦´´ − 2𝑦´´ − 2 𝑥 − 1 𝑦´´´ … (4)
Para x=1 𝑦 𝑖𝑣 1 = −3𝑦´´ − 2 𝑥 − 1 𝑦´´´ = −3
Derivando la ecuación 4: 𝑦 𝑣 = −3𝑦´´´ − 2𝑦´´´ − 2 𝑥 − 1 𝑦 𝑖𝑣
Para x=1 𝑦 𝑣 1 = −5𝑦´´´ − 2 𝑥 − 1 𝑦 𝑖𝑣 = 5
La solución será
∞
𝑦 𝑘 (1) 𝑘
𝑦´´(1) 𝑦´´´(1)
y x = 𝑥−1 = 𝑦 1 + 𝑦´ 1 𝑥 − 1 + (𝑥 − 1)2 + (𝑥 − 1)3 …
𝑘! 2! 3!
𝑘=0
1 1 3 5
𝑦 𝑥 =1+ 𝑥−1 + (𝑥 − 1)2 − 𝑥−1 3 − 𝑥−1 4 + 𝑥−1 5…
2! 3! 4! 5!
Problema 2:
Resolver por serie de Taylor alrededor de 𝑥0 = 0 la ecuación:
𝑦 ′′ − 𝑒 𝑥 𝑦 = 0 , y(0) = 1 ; y’(0) = 1
Solución:
Para 𝑥0 = 0 , tenemos la serie de Maclaurin:
𝒙𝟐 𝒙𝟑
𝒚 𝒙 =𝒚 𝟎 + 𝒚′ 𝟎 𝒙+ 𝒚′′ 𝟎 + 𝒚′′′ 𝟎 + …
𝟐! 𝟑!
De los datos , y(0) = 1 ; y’(0) = 1 en la ecuación:
→ 𝑦 ′′ 0 − 𝑒 0 𝑦 0 = 0 → 𝑦 ′′ 0 = 1
Derivando la E.D.
𝑦 ′′′ − 𝑒 𝑥 𝑦 − 𝑒 𝑥 𝑦 ′ = 0
Analizando en 0:
𝑦 ′′′ (0) − 𝑒 0 𝑦(0) − 𝑒 0 𝑦 ′ (0) = 0 → 𝑦 ′′′ 0 = 2
Derivando la E.D. otra vez:
𝒚𝒊𝒗 − 𝒆𝒙 𝒚 − 𝒆𝒙 𝒚′ − 𝒆𝒙 𝒚′ − 𝒆𝒙 𝒚′′ = 𝟎
֜ 𝒚𝒊𝒗 − 𝒆𝒙 𝒚′′ − 𝟐𝒆𝒙 𝒚′ − 𝒆𝒙 𝒚 = 𝟎
Evaluando en 0:
𝟎 ′′ 𝟎 ′
𝒚𝒊𝒗 𝟎
(𝟎) − 𝒆 𝒚 𝟎 − 𝟐𝒆 𝒚 𝟎 − 𝒆 𝒚 𝟎 = 𝟎 → 𝒚𝒊𝒗
(𝟎) = 𝟒
En la serie:
𝒙𝟐 𝒙𝟑 𝒙𝟒
𝒚 𝒙 = 𝟏 + 𝟏𝒙 + 𝟏 + 𝟐 + 𝟒 …
𝟐! 𝟑! 𝟒!
𝒙𝟐 𝒙𝟑 𝒙𝟒
֜ 𝒚 𝒙 =𝟏+𝒙+ + + …
𝟐 𝟑 𝟔
Nota: Si no dieran las condiciones iniciales, entonces asumir:
′
𝑦(0) = 𝐶1 ; 𝑦(0) = 𝐶2 y obtener la solución general.
SOLUCIÓN DE E.D. ALREDEDOR DE PUNTOS
SINGULARES REGULARES
Sea la E.D:
𝑦´´ + 𝑃(𝑥)𝑦´ + 𝑄(𝑥)𝑦 = 0 ….…. 𝐼
Si se tiene un punto singular en x = x0. (𝑃(𝑥) y Q(𝑥) no son
funciones analíticas en x = x0 )
La ecuación 𝐼 se puede transformar en:
𝑝 𝑥 𝑞 𝑥
𝑦´´ + 𝑦´ + 𝑦 = 0 ….…. 𝐼𝐼
(𝑥 − x0) (𝑥 − x0)2
𝑝 𝑥 𝑞 𝑥
Siendo: 𝑃 𝑥 = y Q 𝑥 =
(𝑥−𝑥0 ) (𝑥−𝑥0 )2
Si de la ecuación 𝐼𝐼 : 𝑝 𝑥 y 𝑞(𝑥) son analíticas en x = x0 ;
entonces decimos que tenemos un punto singular regular para x = x0
(caso contrario se tiene un punto singular irregular).
Si 𝑥 = 𝑥0 es un punto singular regular de la ecuación diferencial
𝑎2 𝑥 𝑦 ′′ + 𝑎1 𝑥 𝑦 ′ + 𝑎0 𝑥 𝑦 = 0 existe al menos una solución
en series de potencias de la forma:
∞ ∞
𝑦 = (𝑥 − 𝑥0 )𝑟 𝑐𝑛 (𝑥 − 𝑥0 )𝑛 = 𝑐𝑛 (𝑥 − 𝑥0 )𝑛+𝑟
𝑛=0 𝑛=0
Donde el número r es una constante a determinar.
Ecuación indicial:
Sea 𝑥0 = 0 un punto singular regular ; reemplazando en la
ecuación 𝐼𝐼 tenemos:
𝑝 𝑥 𝑞 𝑥
𝑦´´ + 𝑦´ + 2 𝑦 = 0 ….…. 𝐼𝐼𝐼
𝑥 𝑥
Entonces la ecuación cuadrática:
𝑟 𝑟 − 1 + 𝑝(0) 𝑟 + 𝑞(0) = 0
se denomina ECUACIÓN INDICIAL de la ecuación diferencial y
nos dará las raíces 𝑟1 𝑦 𝑟2 con los cuales obtendremos la
solución de la ecuación (III) mediante la serie de Frobenius:
∞
𝑦 = 𝑎𝑘 𝑥 𝑘+𝑟
𝑘=0
DERIVADAS:
∞
𝑦′(𝑥) = (𝑘 + 𝑟)𝑎𝑘 𝑥 (𝑘+𝑟−1)
𝑘=0
∞
𝑦′′(𝑥) = (𝑘 + 𝑟)(𝑘 + 𝑟 − 1)𝑎𝑘 𝑥 (𝑘+𝑟−2)
𝑘=0
Nota: El método de Frobenius sólo se aplicará para cuando las
raíces de la ecuación indicial sean REALES.
CASOS:
Al utilizar el método de Frobenius para resolver la ecuación diferencial
de segundo orden se distinguen tres casos que corresponden a la
naturaleza de las raíces reales r1 y r2 :
CASO I :
𝑆𝑖: 𝑟1 ≠ 𝑟2 𝑦 𝑟1 − 𝑟2 ∉ 𝑍
Entonces:
∞ ∞
𝑘 𝑘
𝑦1(𝑥) = 𝑥 𝑟1
𝑎𝑘 𝑥 ; 𝑦2(𝑥) = 𝑥 𝑟2
𝑏𝑘 𝑥
𝑘=0 𝑘=0
CASO II :
𝑆𝑖: 𝑟1 = 𝑟2 = 𝑟
Entonces:
∞ ∞
𝑘 𝑘
𝑦1(𝑥) = 𝑥 𝑎𝑘 𝑥
𝑟
; 𝑟
𝑦2(𝑥) = 𝑥 𝑏𝑘 𝑥 + 𝑦1(𝑥) 𝑙𝑛𝑥
𝑘=0 𝑘=0
CASO III :
𝑆𝑖: 𝑟1 ≠ 𝑟2 𝑦 𝑟1 − 𝑟2 ∈ 𝑍
Entonces:
∞ ∞
𝑦1(𝑥) = 𝑥 𝑟1 𝑎𝑘 𝑥𝑘 ; 𝑦2(𝑥) = 𝑥 𝑟2 𝑏𝑘 𝑥 𝑘
+ 𝐶𝑦
1(𝑥) ln(𝑥)
𝑘=0 𝑘=0
Para los 3 casos:
S.G. : 𝑦(𝑥) = 𝐶1 𝑦1(𝑥) + 𝐶2 𝑦2(𝑥)
PROBLEMA :
Encontrar la solución general de la ecuación diferencial:
2𝑥 2 𝑦 ′′ − 𝑥𝑦 ′ + 1 + 𝑥 𝑦 = 0
Alrededor del punto singular 𝑥0 = 0
SOLUCIÓN:
′
𝑥𝑦 1+𝑥 𝑦
𝑦 ′′ − 2 + =0 (𝑥0 = 0 𝑒𝑠 𝑢𝑛 𝑝𝑢𝑛𝑡𝑜 𝑠𝑖𝑛𝑔𝑢𝑙𝑎𝑟 𝑟𝑒𝑔𝑢𝑙𝑎𝑟)
2𝑥 2𝑥 2
Luego:
∞
𝑌1 (𝑥) = 𝑐𝑛 𝑥 𝑛+𝑟
𝑛=0
∞
𝑌′1 (𝑥) = (𝑛 + 𝑟)𝑐𝑛 𝑥 𝑛+𝑟−1
𝑛=0
∞
𝑌′′1 (𝑥) = (𝑛 + 𝑟)(𝑛 + 𝑟 − 1)𝑐𝑛 𝑥 𝑛+𝑟−2
𝑛=0
Reemplazando en la ecuación diferencial dado:
∞ ∞ ∞
2𝑥 2 (𝑛 + 𝑟)(𝑛 + 𝑟 − 1)𝑐𝑛 𝑥 𝑛+𝑟−2 − 𝑥 (𝑛 + 𝑟)𝑐𝑛 𝑥 𝑛+𝑟−1 + (1 + 𝑥) 𝑐𝑛 𝑥 𝑛+𝑟 = 0
𝑛=0 𝑛=0 𝑛=0
∞ ∞ ∞ ∞
2(𝑛 + 𝑟)(𝑛 + 𝑟 − 1)𝑐𝑛 𝑥 𝑛+𝑟 − 𝑛 + 𝑟 𝑐𝑛 𝑥 𝑛+𝑟 − 𝑐𝑛 𝑥 𝑛+𝑟 + 𝑐𝑛 𝑥 𝑛+𝑟+1 = 0
𝑛=0 𝑛=0 𝑛=0 𝑛=0
∞ ∞
(𝑛 + 𝑟 − 1)(2𝑛 + 2𝑟 − 1)𝑐𝑛 𝑥 𝑛+𝑟 + 𝑐𝑛−1 𝑥 𝑛+𝑟 = 0
𝑛=0 𝑛=1
∞
(𝑟 − 1)(2𝑟 − 1)𝑐0 𝑥 𝑟 + [(𝑛 + 𝑟 − 1)(2𝑛 + 2𝑟 − 1)𝑐𝑛 + 𝑐𝑛−1 ]𝑥 𝑛+𝑟 = 0 . . . (𝐼)
𝑛=1
Aplicando el método de los coeficientes indeterminados:
𝑐𝑛−1
𝑟1 = 1 , 𝑟2 = 1/2 ; 𝑐𝑛 = − … … … (𝐼𝐼)
𝑛+𝑟−1 2𝑛+2𝑟−1
Para 𝑛 ≥ 1
PRIMERA SOLUCIÓN : Para r=1
𝑐𝑛−1
𝑐𝑛 = − Para 𝑛 ≥ 1
(𝑛)(2𝑛+1)
𝑐0
𝑛 = 1: 𝑐1 = −
1.3
𝑐1 𝑐0
𝑛 = 2: 𝑐2 =− =
2.5 1.2.3.5
𝑐2 𝑐0
𝑛 = 3: 𝑐3 =− =−
3.7 1.2. 32 . 5.7
𝑐3 𝑐0
𝑛 = 4: 𝑐4 =− =
4.9 1.2. 32 4.5.7.9
𝑐4 𝑐0
𝑛 = 5: 𝑐5 =− =−
5.11 1.2. 32 4. 52 . 7.9.11
Como:
∞ ∞
𝑌1 (𝑥) = 𝑐𝑛 𝑥 𝑛+𝑟 = 𝑥 𝑐𝑛 𝑥 𝑛
𝑛=0 𝑛=0
∞ ∞
֜ 𝑌1 (𝑥) = 𝑐𝑛 𝑥 𝑛+𝑟 = 𝑥 𝑐𝑛 𝑥 𝑛
𝑛=0 𝑛=0
֜ 𝑌1 (𝑥) = 𝑥(𝐶0 + 𝐶1 𝑥 + 𝐶2 𝑥 2 + 𝐶3 𝑥 3 + 𝐶4 𝑥 4 + 𝐶5 𝑥 5 + … … . . )
𝑐0 𝑐0 2
𝑐0 3
𝑐0
֜ 𝑌1 (𝑥) = 𝑥(𝐶0 − 𝑥+ 𝑥 − 2
𝑥 + 2
𝑥4 +
1.3 1.2.3.5 1.2. 3 . 5.7 1.2. 3 4.5.7.9
𝑐0
− 2 2 𝑥5 + …………..)
1.2.3 4.5 .7.9.11
1 1 1 1
֜ 𝑌1 (𝑥) = 𝑥𝐶0 (1 − 𝑥+ 𝑥2 − 𝑥 3
+ 𝑥 4
+
1.3 1.2.3.5 1.2. 32 . 5.7 1.2. 32 4.5.7.9
1
− 2 2 𝑥5 + …………..)
1.2.3 4.5 .7.9.11
SEGUNDA SOLUCIÓN :
1 1
𝑟1 −𝑟2 = 1 − = no es un entero. (CASO I)
2 2
∞
֜ 𝑌2 (𝑥) = 𝑥1/2 𝑏𝑛 𝑥 𝑛
𝑛=0
Reemplazando 𝑟2 en (II) se obtiene:
𝑏𝑛−1 𝑏𝑛−1
𝑏𝑛 = − =− , ∀𝑛 ≥ 1
1 𝑛 2𝑛 − 1
𝑛 − 2𝑛
2
Para:
𝑏0
𝑛 = 1: 𝑏1 = −
1
𝑏1 𝑏0
𝑛 = 2: 𝑏2 = − =
2.3 1.2.3
𝑏2 𝑏0
𝑛 = 3: 𝑏3 = − =−
3.5 1.2. 32 . 5
𝑏3 𝑏0
𝑛 = 4: 𝑏4 = − =
4.7 1.2. 32 . 4.5.7
Como:
∞
𝑌2 (𝑥) = 𝑥1/2 𝑏𝑛 𝑥 𝑛
𝑛=0
1
֜ 𝑌2 (𝑥) = 𝑥 2 (𝑏0 + 𝑏1 𝑥 + 𝑏2 𝑥 2 + 𝑏3 𝑥 3 + 𝑏4 𝑥 4 + 𝑏5 𝑥 5 + … … . . )
1 𝑏0 2 𝑏0 3
𝑏0
𝑌2 𝑥 = 𝑥2 𝑏0 − 𝑏0 𝑥 + 𝑥 − 2
𝑥 + 2
𝑥 4 + 𝑏5 𝑥 5 + … … . .
1.2.3 1.2. 3 . 5 1.2. 3 . 4.5.7
Luego, la solución general será:
𝑌 𝑥 = 𝐶0 𝑌1 𝑥 + 𝑏0 𝑌2 𝑥