Clase13 MM NVCT 2023
Clase13 MM NVCT 2023
Clase13 MM NVCT 2023
Método de Frobenius
∑[ ]
∞
2 k+λ
aλ(k + λ)(k + λ − 1) + aλ(k + λ) − n aλ x +
λ=0
∞
k+λ+2
∑
+ aλx =0
λ=0



Neil Vladimir Corzo Trejo 3
Limitaciones.
La ecuación indicial se obtiene, en este caso, para λ = 0 y del coe ciente de x k , que es la
potencia más baja que aparece en la ODE.
2
ao[k(k − 1) + k − n ] = 0
Como a0 ≠ 0 entonces:
2
[k − n ] = 0 2 ⟹ k=±n
k+1
Podemos también sacar el siguiente coe ciente, el que va con x , ya que éste también debe de
ser cero.
2
a1[(k + 1)k + (k + 1) − n ] = 0
1
Para k = ± n, la ecuación anterior es cero sólo si a1 = 0 (ó el caso especial k = ± n = − ).
2
2 k 2 k+1
ao[k(k − 1) + k − n ]x + a1[(k + 1)k + (k + 1) − n ]x +
∑[ ]
∞
2 k+λ
+ (k + λ)(k + λ − 1) + (k + λ) − n aλx +
λ=2
∞
k+λ+2
∑
+ aλx =0
Neil Vladimir Corzo Trejo λ=0 4
fi
fi
Limitaciones.
Una vez que ya sacamos los dos primeros términos, las dos series empiezan a tener potencias
iguales. Vamos a recorrer la primer serie usando λ = j + 2 que implica j = λ − 2 . En la segunda
serie sólo remplazaremos λ → j . También ya remplazaremos k = n , la raíz mayor de la ecuación
indicial. Primero vamos a simpli car:
∑[ ]
∞ ∞
2 k+λ k+λ+2
∑
(k + λ)(k + λ − 1) + (k + λ) − n aλx + aλx =0
λ=2 λ=0
∑[ ]
∞ ∞
2 n+λ n+λ+2
∑
(n + λ)(n + λ − 1) + (n + λ) − n aλx + aλx =0
λ=2 λ=0
∑[ ]
∞ ∞
2 2 n+λ n+λ+2
∑
(n + λ) − n aλx + aλx =0
λ=2 λ=0
∑[ ]
∞ ∞
n+λ n+λ+2
∑
(2n + λ)λ aλx + aλx =0
λ=2 λ=0
Neil Vladimir Corzo Trejo 5
fi
Limitaciones.
Una vez que ya sacamos los dos primeros términos, las dos series empiezan a tener potencias
iguales. Vamos a recorrer la primer serie usando λ = j + 2 que implica j = λ − 2 . En la segunda
serie sólo remplazaremos λ → j . También ya remplazaremos k = n , la raíz mayor de la ecuación
indicial.
∑[ ]
∞ ∞
n+λ n+λ+2
∑
(2n + λ)λ aλx + aλx =0
λ=2 λ=0
∑[ ]
∞ ∞
n+j+2 n+λ+2
∑
(2n + j + 2)( j + 2) aj+2x + aλx =0
j=0 λ=0
∑[ ]
∞ ∞
n+j+2 n+j+2
∑
(2n + j + 2)( j + 2) aj+2x + aj x =0
j=0 j=0
∑[ ]
∞
n+j+2
(2n + j + 2)( j + 2)aj+2 + aj x =0
j=0
Neil Vladimir Corzo Trejo 6
1 1
a2 = − a0 = − a0
Limitaciones. 2(2n + 2) 2 ⋅ 2(n + 1)
Entonces la relación de recurrencia es: (n + 1)! = (n + 1)n!
aj
(2n + j + 2)( j + 2)aj+2 + aj = 0 ⟹ aj+2 = −
(2n + j + 2)( j + 2)
Los primeros términos son:
1 a0n!
a2 = − a0 =− 2
2(2n + 2) 2 1!(n + 1)!
1 a0n!
a4 = − a2 = 4
4(2n + 4) 2 2!(n + 2)!
1 a0n!
a6 = − a4 =− 6
6(2n + 6) 2 3!(n + 3)!
En general: a0n!
p
a2p = (−1)
22p p!(n + p)!
Neil Vladimir Corzo Trejo 7
Limitaciones.
[ ]
2 4
La solución es:
n n!x n!x
y(x) = a0 x 1 − 2 + 4 −⋯
2 1!(n + 1)! 2 2!(n + 2)!
Y en forma de sumatoria tenemos: ∞ n+2j
j n!x
∑
y(x) = a0 (−1) 2j
j=0
2 j!(n + j)!
(2)
∞ n+2j
n j 1 x
∑
y(x) = a02 n! (−1)
j=0
j!(n + j)!
Esta solución es la función de Bessel Jn(x).
Si n no es un entero, se pueden generar dos soluciones distintas. La segunda solución se puede
expresar como J−n(x).
Si n es un entero, NO podemos generar una segunda solución con la otra raíz de la ecuación
indicial k = − n. Existe un problema porque la relación de recurrencia falla y diverge.
aj
aj+2 = − →∞ j + 2 = 2n
(−2n + j + 2)( j + 2)
Neil Vladimir Corzo Trejo 8
Limitaciones.
Es más, tenemos que para n entero:
n
J−n(x) = (−1) Jn(x)
Y no tenemos una segunda solución linealmente independiente.
Tenemos un ejemplo donde el método de Frobenius sólo nos dio una solución para el caso donde
las diferencia de las raíces de la ecuación indicial es un número entero.
∑[ ]
∞
k+n−2
an(k + n)(k + n − 1) − 6an x =0
n=0
La ecuación indicial se obtiene para n = 0:
a0 ≠ 0
a0k(k − 1) − 6a0 ⟹ 2
k(k − 1) − 6 = k − k − 6 = (k + 2)(k − 3) = 0
Además NO tenemos una relación de recurrencia porque el coe ciente que va con x k+n−2 no
involucra diferentes an’s.
3 1
La solución general es: y(x) = c x + c 1 2 2
x


Neil Vladimir Corzo Trejo 13
fi
Ejemplo 76b
La segunda ecuación es la siguiente: 6
y′′ − 3 y = 0
x
Note que lo único diferente con la ecuación anterior es el término q(x).
Es fácil ver que x = 0 es un punto singular, además como :
2 2 −6
lim x q(x) = lim x 3 → − ∞
x→0 x→0 x
Por lo que x = 0 es un punto singular irregular.
∑[ ]
∞
2 k+n−2
an(k + n)(k + n − 1) + an(k + n) − c an x =0
n=0
No existe relación de recurrencia, y la ecuación indicial es:
2
a0k(k − 1) + a0k − c a0 = 0 ⟹ k=±c
a0 ≠ 0
La solución general es: c 1
y(x) = c1x + c2 c



Neil Vladimir Corzo Trejo 16
x
Ejemplo 76d 1 c 2
La cuarta ecuación es la siguiente: y′′ + 2 y′ − 2 y = 0
x x
Note que lo único diferente con la ecuación anterior es el término p(x).
Además, tenemos que: x
lim xp(x) = lim 2 → ∞ .
x→0 x→0 x
Por lo que x = 0 es un punto singular irregular.
Sustituyendo la solución propuesta en la ODE tenemos:
2
1 c
y′′ + 2 y′ − 2 y = 0 ⟹
x x
∞ ∞ 2 ∞
k+n−2 1 k+n−1 c k+n
∑ ∑ ∑
an(k + n)(k + n − 1)x + 2 an(k + n)x − 2 anx = 0
n=0
x n=0
x n=0
∞ ∞ ∞
k+n−2 k+n−3 2 k+n−2
∑ ∑ ∑
an(k + n)(k + n − 1)x + an(k + n)x −c anx =0
n=0 n=0 n=0






Neil Vladimir Corzo Trejo 17
Ejemplo 76d
La ecuación indicial se obtiene del primer término de la segunda serie:
∞ ∞ ∞
k+n−2 k+n−3 2 k+n−2
∑ ∑ ∑
an(k + n)(k + n − 1)x + an(k + n)x −c anx =0
n=0 n=0 n=0
a0(k) = 0 ⟹ k = 0
Ahora tenemos k = 0 y además tenemos que recorrer la segunda serie.
∞ ∞ ∞
k+n−2 k+n−3 2 k+n−2
∑ ∑ ∑
an(k + n)(k + n − 1)x + an(k + n)x −c anx =0
n=0 n=1 n=0
∞ ∞ ∞
k+n−2 k+n−2 2 k+n−2
∑ ∑ ∑
an(k + n)(k + n − 1)x + an+1(k + n + 1)x −c anx =0
n=0 n=0 n=0
∑[ ]
∞ 2
2 n−2 c − n(n − 1)
ann(n − 1) + an+1(n + 1) − c an x =0 an+1 = an
n=0
n+1
Neil Vladimir Corzo Trejo 18
Ejemplo 76d 2
c − n(n − 1)
La relación de recurrencia es: an+1 = an
n+1
El problema es que la serie resultante NO va a converger:
an+1 n(n − 1)
lim = lim →∞
n→∞ an n→∞ n + 1
[ ∫a ]
x
W(x) = W(a)exp − p(x1)dx1 .
La constante W(a) es simplemente una constantes, y como a las soluciones y1,2 se les puede
multiplicar por cualquier constante y seguirán siendo soluciones, entonces podemos tomar
W(a) = 1.
Ahora, el Wronskiano se puede reescribir como:
dx ( y1 )
2 d y2
W(x) = y1y′2 − y′1y2 = y1


Neil Vladimir Corzo Trejo











22
Segunda Solución
Combinando ambas ecuaciones tenemos:
[ ∫a ] dx ( y1 )
x
2 d y2
W(x) = exp − p(x1)dx1 . W(x) = y1y′2 − y′1y2 = y1 .
dx ( y1 ) y1 [ ∫ ]
x
d y2 1
= 2 exp − p(x1)dx1 .
∫
y2(x) = y1(x) dx2
[y1(x2)]2
Ésta es la misma ecuación que habíamos obtenido antes:
[ ∫ ]
1
∫
U = 2 exp − pdx y2 = uy1 = y1 Udx
y1


Neil Vladimir Corzo Trejo 23
Ejemplo 77
Encuentre la segunda solución para la Ecuación de Bessel, con n = 0.
Solución. La ecuación de Bessel es:
1
2 2 2
x y′′ + xy′ + (x − n )y = 0 ⟹ y′′ + y′ + y = 0
En este caso las funciones son: x
1
p(x) = q(x) = 1 n=0
x
La ecuación indicial queda como:
2 2 2
k =n → k =0
La primer solución, tomando a0 = 1, es:
2 4
x x 6
y1(x) = J0(x) = 1 − + − O(x )
4 64
La segunda solución puede ser encontrada como:
exp[ − ∫ x dx1]
x2 1
x
∫
y2(x) = J0(x) 1
dx2
J0 (x)
2






Neil Vladimir Corzo Trejo 24
Ejemplo 77
Sustituyendo J0(x) y p(x) tenemos:
x2
dx1 1
∫ x1
exp[ − ] = exp[−ln x2] =
x2
Por otra parte, utilizando la expansión binomial:
2 4 2 4
1 x2 x2 x2 5x2
[ ⋯]
−2
= 1 − + + =1+ + +⋯
[J0(x)] 2 4 64 2 32
La segunda solución es entonces: x 2 4
1 x2 5x2
∫ x2 [
y2(x) = J0(x) 1+ + + ⋯]dx2
2 32
x 3
1 x2 5x2
∫ x2
y2(x) = J0(x) [ + + + ⋯]dx2
2 32
2 4
x 5x
y2(x) = J0(x){ln x + + + ⋯}
2 128
Neil Vladimir Corzo Trejo 25
Ejemplo 77
La segunda solución es por lo tanto:
2 4
x 5x
y2(x) = J0(x){ln x + + + ⋯}
2 128
Podemos poner esta solución en una forma más estándar. Primero la podemos multiplicar por
una contante (2/π) y sumarle un múltiplo de la primera solución:
2
[ln 2 + γ]J0(x)
π
Donde γ es una constante conocida como la constante de Euler-Mascheroni.
p−1 ≠ 0 q−2 ≠ 0
Vamos a obtener los primeros términos de la primera solución de la ODE y, posteriormente,
usar la fórmula para encontrar la segunda solución. Vamos también a tomar el caso donde las
raíces de la ecuación indicial di eren por un número entero (si este no es el caso, el Método de
Frobenius nos da dos soluciones).



Neil Vladimir Corzo Trejo 27
fi
Segunda Solución para una Serie
Proponiendo la solución como:
∞ ∞
k+n k+n−1
∑ ∑
y(x) = anx , y′ = an(k + n)x ,
n=0 ∞ n=0 a0 ≠ 0.
k+n−2
∑
y′′ = an(k + n)(k + n − 1)x .
n=0
Y sustituyendo en la ecuación diferencial tenemos: y′′ + p(x)y′ + q(x)y = 0
∞ ∞ ∞ ∞ ∞
k+n−2 i k+n−1 j k+n
∑ ∑ ∑ ∑ ∑
an(k + n)(k + n − 1)x + pi x an(k + n)x + qj x anx =0
n=0 i=−1 n=0 j=−2 n=0
k−2
La ecuación indicial se obtiene con la potencia x y aparece en todos los términos. Ésta es:
2 (k − α)(k − α + m) = 0
k + (p−1 − 1)k + q−2 = 0 ⟹
Las raíces son k = α y k = α − m. Tomaremos la raíz mayor para la primera solución.
Además, tenemos que:
2
(k − α)(k − α + m) = k + (m − 2α)k + α(α − m)
⟹
p−1 − 1 = m − 2α
La primera solución queda entonces como:
∞
α n
∑
y1(x) = x anx ,
exp[ − ∫ p(x1)dx1]
x x2
∫
n=0
Para encontrar la segunda solución tenemos que usar: y2(x) = y1(x) dx2
[y1(x2)]2
Neil Vladimir Corzo Trejo 29
fi
Segunda Solución para una Serie
exp[ − ∫ p(x1)dx1]
Entonces: x x2
∫
y2(x) = y1(x) dx2
[y1(x2)]2
[ 1]
x exp − ∫ ∑ ∞x2 i
p x dx
∫
i=−1 i 1
y2(x) = y1(x) ∞ dx2
x22α[ ∑n=0 anx2n)]2
[ ∫ ∑ ] ( ∑k+1 )
x2 ∞ ∞
i −p−1 pk k+1
exp − pi x1dx1 = x2 exp − x2
i=−1 k=0
Neil Vladimir Corzo Trejo 30
Segunda Solución para una Serie
Entonces:
[ ∫ ∑ ] ( )
x2 ∞ ∞
i −p−1 pk k+1
∑k+1
exp − pi x1dx1 = x2 exp − x2
i=−1 k=0
[ ∫ ∑ ] [ ( ∑k+1 ) 2! ( ∑ k + 1 ) ]
x2 ∞ ∞ ∞ 2
i −p−1 pk k+1 1 pk k+1
exp − pi x1dx1 = x2 1+ − x2 + − + x2 +⋯
i=−1 k=0 k=0
[ ] (∑ )
∞ −1 ∞ −2 ∞
2α n 2 −2α n −2α λ
∑ ∑
x2 [ anx2 )] = x2 anx2 = x2 bλx2
n=0 n=0 λ=0
(∑ )
x ∞
∫
(−p−1−2α) n
y2(x) = y1(x) x2 cnx2 dx2
n=0
p−1 − 1 = m − 2α ⟹ (−p−1−2α)
x2 =x −m−1
Sustituyendo:
x
∫ (cox2
y2(x) = y1(x) −m−1
+ −m
c1x2 + m+1
c2x2 +⋯+ −1
cm x2 + ⋯)dx2
Entonces la segunda solución va a tener dos propiedades:
−m
1. La serie empieza en x , donde m es la diferencia de las raíces de la ecuación indicial.
2. Existe un término logarítmico que siempre aparecerá (salvo si cm = 0).
Arfken
Sección 9.6: 4, 9, 11, 13, 15, 18, 19, 20, 21, 26
Tarea Extra
Sección 9.6: 3, 6, 19, 23.