ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN
DEFINICIÓN: Una ecuación diferencial de segundo orden puede escribirse de forma implícita como
𝐹(𝑥, 𝑦, 𝑦 ′ , 𝑦′′) = 0
Las ecuaciones diferenciales de segundo orden que son muy importantes en la física son las ecuaciones lineales.
DEFINICIÓN: Una ecuación de segundo orden lineal es de la forma:
𝑦 ′′ + 𝑝 𝑥 𝑦 ′ + 𝑞 𝑥 𝑦 = 𝑟 𝑥
Cuando la función r es idénticamente nula diremos que la ecuación lineal es homogénea y si esto no ocurre, diremos
que es no homogénea.
EJEMPLOS:
La ecuación diferencial 𝑦 ′′ − cos 𝑥 𝑦′ + 𝑥𝑦 = 0 es una ecuación lineal homogénea.
La ecuación 𝑦 ′′ + 𝑙𝑛𝑥 𝑦′ = 2𝑒 𝑥 es una ecuación lineal no homogénea.
Ecuaciones lineales de segundo orden homogéneas y con coeficientes constantes
DEFINICIÓN: Una ecuación diferencial lineal de segundo orden homogénea y con coeficientes constantes tiene la
forma
𝑎𝑦′′ + 𝑏𝑦′ + 𝑐𝑦 = 0 (∗)
donde 𝑎 ≠ 0, 𝑏 y 𝑐 son constantes.
El polinomio
𝑝 𝑥 = 𝑎𝑥 2 + 𝑏𝑥 + 𝑐
asociado a la ecuación homogénea (*) es llamado el polinomio característico de la ecuación.
Claramente, las raíces 𝑟1 y 𝑟2 del polinomio característico 𝑝 pueden ser reales distintas, reales iguales o
complejas. En el caso que las raíces de 𝑝 sean complejas, una debe ser el conjugado de
la otra, pues recordemos que las raíces complejas vienen de a pares.
El siguiente resultado nos indica como encontrar todas las soluciones de esta ecuación diferencial teniendo en cuenta
las naturaleza de las raíces de su polinomio característico.
TEOREMA 1:
Si 𝑟1 y 𝑟2 son reales y distintas entonces cualquier solución de (2) es de la forma
𝑦(𝑥) = 𝐶1 𝑒 𝑟1 𝑥 + 𝐶2 𝑒 𝑟2𝑥
Si 𝑟1 y 𝑟2 son reales y 𝑟1 = 𝑟2 entonces cualquier solución de (2) es de la forma
𝑦(𝑥) = 𝐶1 + 𝐶2 𝑥 𝑒 𝑟1 𝑥
Si 𝑟1 = 𝛼 + 𝛽𝑖 y 𝑟2 = 𝛼 − 𝛽𝑖 entonces cualquier solución de (2) es de la forma
𝑦(𝑥) = 𝐶1 cos 𝛽𝑥 + 𝐶2 sin 𝛽𝑥 𝑒 𝛼𝑥
Al igual que con las ecuaciones de primer orden, en ocasiones buscamos resolver un problema con valores iniciales, es
decir buscamos una función que sea solución de la ecuación diferencial dada pero además queremos que cumpla
alguna condición extra.
Cuando la ecuación es de segundo orden las condiciones deben ser dos.
EJEMPLO: Encontrar la solución del siguiente problema de valores iniciales:
𝑦′′ + 2𝑦′ + 𝑦 = 0
𝑦(0) = 5
𝑦 ′ 0 = −3
El polinomio característico asociado es
𝑝(𝑥) = 𝑥 2 + 2𝑥 + 1
y sus raíces son 𝑟1 = −1 = 𝑟2 .
Así, la solución general de la ecuación es de la forma:
𝑦(𝑥) = 𝐶1 + 𝐶2 𝑥 𝑒 −𝑥 (∗)
Ahora usamos las condiciones iniciales para encontrar los valores de 𝐶1 y 𝐶2 .
Reemplazando la primera condición en (*) obtenemos:
5 = 𝑦 0 = 𝐶1 + 𝐶2 0 𝑒 −0 = 𝐶1 ⇒ 𝐶1 = 5
Si derivamos (*) y usamos la segunda condición obtenemos:
𝑑
−3 = 𝑦 ′ 0 = 𝐶1 + 𝐶2 𝑥 𝑒 −𝑥 = 𝐶2 𝑒 −𝑥 − 𝐶1 + 𝐶2 𝑥 𝑒 −𝑥 = 𝐶2 𝑒 −0 − 𝐶1 + 𝐶2 0 𝑒 −0 = 𝐶2 − 𝐶1 ⇒ 𝐶2 = 2
𝑑𝑥 𝑥=0
𝑥=0
Finalmente, la solución del problema de valores iniciales es:
𝑦(𝑥) = 5 + 2𝑥 𝑒 𝑥
Ecuaciones diferenciales lineales de segundo orden a coeficientes constantes no
homogéneas
Las ecuaciones diferenciales lineales de segundo orden a coeficientes constantes no homogéneas son de la forma:
𝑦 ′′ + 𝑏 𝑦 ′ + 𝑐 𝑦 = 𝑟 𝑥 (∗)
Si conocemos alguna solución particular 𝑦𝑝 de la ecuación dada, el siguiente resultado nos indica como encontrar
todas las soluciones de ésta ecuación.
TEOREMA 2: Si 𝑦𝑝 una solución particular de la ecuación (*) entonces cualquier solución de la misma es de la forma:
𝑦(𝑥) = 𝑦ℎ (𝑥) + 𝑦𝑝 (𝑥)
donde 𝑦ℎ es la solución general de la ecuación homogénea asociada,
𝑦 ′′ + 𝑏 𝑦 ′ + 𝑐 𝑦 = 0
dada por el TEOREMA 1.
A continuación vamos a estudiar métodos que nos permitan encontrar una solución particular 𝑦𝑝 para una ecuación
de la forma (*).
Método de los coeficientes indeterminados
La siguiente tabla nos da una elección inicial de la solución particular 𝑦𝑝 correspondiente a cada tipo de función 𝑟.
𝒓(𝒙) 𝒚𝒑 (𝒙)
𝐶𝑒 𝑎𝑥 𝐾𝑒 𝑎𝑥
𝐶𝑥 𝑛 𝑎0 + 𝑎1 𝑥 + ⋯ + 𝑎𝑛 𝑥 𝑛
𝐶 cos 𝑎𝑥 ó 𝐶 sen 𝑎𝑥 𝐾1 cos 𝑎𝑥 + 𝐾2 𝑠𝑒𝑛 𝑎𝑥
𝐶 𝑒 𝑎𝑥 cos 𝑏𝑥 ó 𝐶 𝑒 𝑎𝑥 sen 𝑏𝑥 𝑒 𝑎𝑥 𝐾1 cos 𝑏𝑥 + 𝐾2 𝑠𝑒𝑛 𝑏𝑥
Los coeficientes en 𝑦𝑝 serán determinados sustituyendo 𝑦𝑝 en la ecuación no homogénea. Podrá ocurrir que la
solución particular sugerida por la tabla resulte ser solución de la ecuación homogénea asociada, en este caso se
procede de la siguiente manera:
Si el polinomio característico de la ecuación homogénea asociada tiene dos raíces distintas se debe multiplicar por
𝑥 la solución particular 𝑦𝑝 y elegir este producto como solución particular de la no homogénea.
Si el polinomio característico de la ecuación homogénea asociada tiene raíces iguales se debe multiplicar por 𝑥 2 la
solución particular 𝑦𝑝 y elegir este producto como solución particular de la no homogénea.
Además, si 𝑟 es suma de funciones que figuran en la tabla entonces la solución particular 𝑦𝑝 será la suma de las
soluciones particulares correspondientes sugeridas en la tabla.
EJEMPLO: Encontrar la solución general de la ecuación 𝑦 ′′ − 4𝑦′ + 4𝑦 = 𝑒 2𝑥 .
Primero buscamos la solución general de la ecuación lineal homogénea asociada
𝑦 ′′ − 4𝑦′ + 4𝑦 = 0
Como el polinomio característico asociado 𝑝 𝑥 = 𝑥 2 − 4𝑥 + 4 tiene raíces 𝑟1 = 2 = 𝑟2 entonces la solución general
es:
𝑦ℎ (𝑥) = 𝐶1 + 𝐶2 𝑥 𝑒 2𝑥
Segundo, buscamos una solución particular de la ecuación no homogénea. Como
𝑟 (𝑥) = 𝑒 2𝑥 ,
la tabla propone
𝑦𝑝 (𝑥) = 𝑘𝑒 2𝑥
como solución particular. Pero observamos que esta función 𝑦𝑝 es solución de la ecuación homogénea asociada, por lo
tanto multiplicamos por 𝑥 2 dado que las raíces son iguales. Así,
𝑦𝑝 (𝑥) = 𝑘𝑥 2 𝑒 2𝑥
Reemplazamos,
𝑦𝑝 𝑥 = 𝑘𝑥 2 𝑒 2𝑥 ; 𝑦 ′ 𝑝 = 𝑘 2𝑥𝑒 2𝑥 + 2𝑥 2 𝑒 2𝑥 ; 𝑦 ′′ 𝑝 (𝑥) = 𝑘(2𝑒 2𝑥 + 8𝑥𝑒 2𝑥 + 4𝑥 2 𝑒 2𝑥 )
en la ecuación no homogénea para determinar 𝑘.
Así,
𝑦𝑝′′ − 4𝑦′𝑝 + 4𝑦 = 𝑒 2𝑥
𝑘 2𝑒 2𝑥 + 8𝑥𝑒 2𝑥 + 4𝑥 2 𝑒 2𝑥 − 4𝑘 2𝑥𝑒 2𝑥 + 2𝑥 2 𝑒 2𝑥 + 4𝑘𝑥 2 𝑒 2𝑥 = 𝑒 2𝑥
𝑘 2 + 8𝑥 + 4𝑥 2 − 8𝑥 − 8𝑥 2 + 4𝑥 2 𝑒 2𝑥 = 𝑒 2𝑥
𝑘2𝑒 2𝑥 = 𝑒 2𝑥
1
𝑘=2
1
y por lo tanto 𝑦𝑝 𝑥 = 2 𝑥 2 𝑒 2𝑥 .
Finalmente, la solución general de la ecuación no homogénea dada es:
1 2 2𝑥
𝑦 𝑥 = 𝐶1 + 𝐶2 𝑥 𝑒 2𝑥 + 𝑥 𝑒
2
Método de series de potencias
Muchas ecuaciones diferenciales de segundo orden
𝑦 ′′ + 𝑝 𝑥 𝑦 ′ + 𝑞(𝑥) 𝑦 = 𝑟 𝑥
tienen al menos una solución que admite desarrollo en serie de potencias
∞
𝑦 𝑥 = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + ⋯ = 𝑎𝑖 𝑥 𝑖
𝑖=0
Este método consiste en encontrar las soluciones, determinando los coeficientes del desarrollo en serie.
Este método puede utilizarse cuando 𝑝, 𝑞 y 𝑟 son polinomios o admiten desarrollo en serie de potencias.
Recordemos que, si 𝑦 𝑥 = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + ⋯ entonces
∞
𝑦 ′ 𝑥 = 𝑎1 + 𝑎2 2𝑥 + 𝑎3 3𝑥 2 + ⋯ = 𝑎𝑖 𝑖 𝑥 𝑖−1
𝑖=1
∞
𝑦 ′′ 𝑥 = 𝑎2 2 + 𝑎3 3 2 𝑥 + 𝑎4 4 3 𝑥 2 + ⋯ = 𝑎𝑖 𝑖(𝑖 − 1) 𝑥 𝑖−2
𝑖=2
El método consiste en reemplazar 𝑦 y sus derivadas en la ecuación no homogénea y encontrar los coeficientes aj del
desarrollo en serie.
EJEMPLO: Encontrar las soluciones de la ecuación 𝑦 ′′ + 𝑥 𝑦 = 0.
∞ 𝑖
Suponemos que la ecuación admite una solución 𝑦 𝑥 = 𝑖=0 𝑎𝑖 𝑥 .
Derivamos, ∞ ∞
𝑦′ 𝑥 = 𝑎𝑖 𝑖𝑥 𝑖−1 𝑦′′(𝑥) = 𝑎𝑖 𝑖(𝑖 − 1)𝑥 𝑖−2
𝑖=1 𝑖=2
Reemplazamos en la ecuación dada,
∞ ∞
𝑎𝑖 𝑖 𝑖 − 1 𝑥 𝑖−2 + 𝑥 𝑎𝑖 𝑥 𝑖 = 0
𝑖=2 𝑖=0
Llevamos el lado derecho de la igualdad a una única serie de potencias,
∞ ∞
𝑎𝑖 𝑖 𝑖 − 1 𝑥 𝑖−2 + 𝑎𝑖 𝑥 𝑖+1 = 0
𝑖=2 𝑖=0
∞ ∞
𝑎𝑗+2 (𝑗 + 2) 𝑗 + 1 𝑥 𝑗 + 𝑎𝑗−1 𝑥 𝑗 = 0
𝑗=0 𝑗=1
∞ ∞
𝑎2 2 + 𝑎𝑗+2 (𝑗 + 2) 𝑗 + 1 𝑥 𝑗 + 𝑎𝑗−1 𝑥 𝑗 = 0
𝑗=1 𝑗=1
∞
2𝑎2 + 𝑎𝑗+2 𝑗 + 2 𝑗 + 1 + 𝑎𝑗−1 𝑥 𝑗 = 0
𝑗=1
De aquí, obtenemos las siguientes condiciones para los coeficientes:
2𝑎2 = 0
𝑎𝑗+2 𝑗 + 2 𝑗 + 1 + 𝑎𝑗−1 = 0, ∀𝑗 ≥ 1
Las cuales son equivalentes a:
𝑎2 = 0
𝑎𝑘−3
𝑎𝑘 = − , ∀𝑘 ≥ 3
𝑘(𝑘 − 1)
Esto deja 𝑎0 y 𝑎1 como variables libres y a partir de estos quedan completamente determinados todos los coeficientes:
𝑎0 𝑎1 𝑎0 𝑎1
𝑎2 = 0, 𝑎3 = − , 𝑎4 = − , 𝑎5 = 0, 𝑎6 = , 𝑎7 = ,⋯
3∙2 4∙3 6∙5∙3∙2 7∙6∙4∙3
De manera más general,
𝑗
𝑎0 𝑗
𝑎1
𝑎3𝑗 = −1 𝑗 , 𝑎3𝑗+1 = −1 𝑗 , 𝑎3𝑗+2 = 0 ∀𝑗 ≥ 0
𝑖=1 3𝑖(3𝑖
− 1) + 1)3𝑖𝑖=1(3𝑖
Finalmente, las soluciones que admiten desarrollo en serie son de la forma
∞ ∞
𝑗 𝑗
−1 −1
𝑦 𝑥 = 𝑎0 𝑗 𝑥 3𝑗 + 𝑎1 𝑗 𝑥 3𝑗+1
𝑗=0 𝑖=1 3𝑖(3𝑖 − 1) 𝑗=0 𝑖=1(3𝑖 + 1)3𝑖
con 𝑎0 y 𝑎1 constantes cualesquiera.
Si a nuestro problema le sumamos las condiciones iniciales 𝑦 0 = 1, 𝑦 ′ 0 = 0 entonces vamos a determinar valores
particulares para las constantes 𝑎0 y 𝑎1 ,
∞ ∞
𝑗 𝑗
−1 −1
1 = 𝑦 0 = 𝑎0 𝑗
03𝑗 + 𝑎1 𝑗
03𝑗+1 = 𝑎0
𝑗=0 𝑖=1 3𝑖(3𝑖 − 1) 𝑗=0 𝑖=1(3𝑖 + 1)3𝑖
∞ ∞
−1 𝑗 −1 𝑗
0 = 𝑦′ 0 = 𝑎0 𝑗
3𝑗 03𝑗−1 + 𝑎1 𝑗
(3𝑗 + 1) 03𝑗 = 𝑎1
𝑗=1 𝑖=1 3𝑖(3𝑖 − 1) 𝑗=0 𝑖=1(3𝑖 + 1)3𝑖
Por lo tanto, la solución particular que resuelve la ecuación
diferencial dada y las condiciones dadas es:
∞
𝑗
−1
𝑦 𝑥 = 𝑗
𝑥 3𝑗
𝑗=0 𝑖=1 3𝑖(3𝑖 − 1)
𝑥3 𝑥6𝑥9
=1− + − +⋯
3∙2 6∙5∙3∙2 9∙8∙6∙5∙3∙2
En la grafica se muestran algunas funciones que son las
primeras sumas parciales de la serie que define la solución.
Se observa que éstas funciones convergen o se aproximan a
una única función para valores de 𝑥 < 1. Esto nos lleva a
concluir que el problema de valores iniciales admite solución
sobre el intervalo 𝐼 = (−∞, 1).
Método de Frobenius
Este método se usa para resolver ecuaciones de la forma
𝑏(𝑥) ′ 𝑐(𝑥)
𝑦 ′′ +
y + 2 𝑦=0 (∗)
𝑥 𝑥
en la cual 𝑏 y 𝑐 son polinomios o admiten desarrollo en serie de potencias. Cualquier ecuación de la forma (*) tiene
al menos una solución que puede ser expresada como
∞
𝑦 𝑥 = 𝑥𝑟 𝑎𝑖 𝑥 𝑖
𝑖=0
donde 𝑟 debe ser determinado y 𝑎 ≠ 0.
Bajo las hipótesis enunciadas para 𝑏 y 𝑐 podemos escribir:
𝑏 𝑥 = 𝑏0 + 𝐵 𝑥 𝑐 𝑥 = 𝑐0 + 𝐶(𝑥)
∞ 𝑖 ∞ 𝑖
donde 𝑏0 y 𝑐0 son constantes y, 𝐵 𝑥 = 𝑖=1 𝑏𝑖 𝑥 y𝐶 𝑥 = 𝑖=1 𝑐𝑖 𝑥 .
Luego, la ecuación (*) es equivalente a:
𝑥 2 𝑦 ′′ + 𝑥 𝑏 𝑥 𝑦 ′ + 𝑐 𝑥 𝑦 = 0 (∗)
Si derivamos 𝑦 obtenemos
∞ ∞ ∞ ∞ ∞
𝑦′ 𝑥 = 𝑟𝑥 𝑟−1 𝑎𝑖 𝑥 𝑖 + 𝑥 𝑟 𝑎𝑖 𝑖 𝑥 𝑖−1 = 𝑥 𝑟−1 𝑟𝑎𝑖 𝑥 𝑖 + 𝑥 𝑟−1 𝑎𝑖 𝑖 𝑥 𝑖 = 𝑥 𝑟−1 𝑟 + 𝑖 𝑎𝑖 𝑥 𝑖
𝑖=0 𝑖=1 𝑖=0 𝑖=1 𝑖=0
Ahora derivamos 𝑦′ y obtenemos,
∞ ∞
𝑦′′ 𝑥 = (𝑟 − 1)𝑥 𝑟−2 𝑟 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑥 𝑟−1 𝑟 + 𝑖 𝑎𝑖 𝑖 𝑥 𝑖−1
𝑖=0 𝑖=1
∞ ∞
= 𝑥 𝑟−2 (𝑟 − 1)(𝑟 + 𝑖)𝑎𝑖 𝑥 𝑖 + 𝑥 𝑟−2 (𝑟 + 𝑖)𝑎𝑖 𝑖 𝑥 𝑖
𝑖=0 𝑖=1
∞
= 𝑥 𝑟−2 (𝑟 − 1 + 𝑖)(𝑟 + 𝑖)𝑎𝑖 𝑥 𝑖
𝑖=0
Reemplazamos los desarrollos en serie que definen a 𝑦, 𝑦′ y 𝑦′′ en la ecuación (*),
∞ ∞ ∞
𝑥 2 𝑥 𝑟−2 (𝑟 − 1 + 𝑖)(𝑟 + 𝑖)𝑎𝑖 𝑥 𝑖 + 𝑥 𝑏 𝑥 𝑥 𝑟−1 𝑟 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑐 𝑥 𝑥 𝑟 𝑎𝑖 𝑥 𝑖 = 0
𝑖=0 𝑖=0 𝑖=0
∞ ∞ ∞
𝑥𝑟 (𝑟 − 1 + 𝑖)(𝑟 + 𝑖)𝑎𝑖 𝑥 𝑖 + 𝑏 𝑥 𝑥 𝑟 𝑟 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑐 𝑥 𝑥 𝑟 𝑎𝑖 𝑥 𝑖 = 0
𝑖=0 𝑖=0 𝑖=0
∞ ∞ ∞
𝑥𝑟 (𝑟 − 1 + 𝑖)(𝑟 + 𝑖)𝑎𝑖 𝑥 𝑖 + 𝑏0 + 𝐵(𝑥) 𝑥 𝑟 𝑟 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑐0 + 𝐶(𝑥) 𝑥 𝑟 𝑎𝑖 𝑥 𝑖 = 0
𝑖=0 𝑖=0 𝑖=0
∞ ∞ ∞
𝑥𝑟 𝑟 − 1 + 𝑖 𝑟 + 𝑖 + 𝑏0 𝑟 + 𝑖 + 𝑐0 𝑎𝑖 𝑥 𝑖 + 𝑥 𝑟 𝐵 𝑥 𝑟 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝐶(𝑥) 𝑎𝑖 𝑥 𝑖 = 0
𝑖=0 𝑖=0 𝑖=0
Notar que 𝐵 𝑥 ∞ 𝑖
𝑖=0 𝑟 + 𝑖 𝑎𝑖 𝑥 + 𝐶(𝑥)
∞
𝑖=0 𝑎𝑖 𝑥
𝑖
es una serie de la forma ∞
𝑖=1 𝑑𝑖 𝑥
𝑖
donde los coeficientes 𝑑𝑖 dependen
de 𝐵 y 𝐶. Así, tenemos que
∞ ∞
𝑥𝑟 𝑟 − 1 + 𝑖 𝑟 + 𝑖 + 𝑏0 𝑟 + 𝑖 + 𝑐0 𝑎𝑖 𝑥 𝑖 + 𝑥 𝑟 𝑑𝑖 𝑥 𝑖 = 0
𝑖=0 𝑖=1
∞ ∞
𝑥 𝑟 𝑟 𝑟 − 1 + 𝑏0 𝑟 + 𝑐0 + 𝑥 𝑟 𝑟 − 1 + 𝑖 𝑟 + 𝑖 + 𝑏0 𝑟 + 𝑖 + 𝑐0 𝑎𝑖 𝑥 𝑖 + 𝑥 𝑟 𝑑𝑖 𝑥 𝑖 = 0
𝑖=1 𝑖=1
∞
𝑥 𝑟 𝑟 𝑟 − 1 + 𝑏0 𝑟 + 𝑐0 + 𝑥 𝑟 𝑟 − 1 + 𝑖 𝑟 + 𝑖 + 𝑏0 𝑟 + 𝑖 + 𝑐0 𝑎𝑖 + 𝑑𝑖 𝑥 𝑖 = 0
𝑖=1
Esta igualdad se cumple si y solo si cada coeficiente de la serie del lado derecho se anula. Lo que nos lleva a la primera
igualdad
𝑟 𝑟 − 1 + 𝑏0 𝑟 + 𝑐0 = 0
conocida como ecuación indicial, sus raíces nos darán los valores posibles para 𝑟.
Una vez que encontramos 𝑟 procedemos de la misma manera que en el método de series de potencias, es decir,
encontramos los coeficientes aj que satisfacen
𝑟 − 1 + 𝑖 𝑟 + 𝑖 + 𝑏0 𝑟 + 𝑖 + 𝑐0 𝑎𝑖 + 𝑑𝑖 = 0, ∀𝑖 ≥ 1
EJEMPLO: ECUACIÓN DE BESSEL
Aplicaremos el método de Frobenius para resolver la ecuación de Bessel
2 ′′ ′ 2 2
1
𝑥 𝑦 + 𝑥𝑦 + 𝑥 − 𝑛 𝑦 = 0, 𝑛=±
En este caso, 2
1
𝑏 𝑥 = 1 ∧ 𝑐 𝑥 = 𝑥 2 − 𝑛2 ⇒ 𝑏0 = 1 ∧ 𝑐0 = −𝑛2 = − 4
y la ecuación indical es:
1
𝑟 𝑟−1 +𝑟−=0
4
1 1 1
cuyas raíces son 𝑟1 = 2 y 𝑟2 = − 2. Tomamos 𝑟 = 2 entonces
∞ ∞ ∞
1 1
−2 1 3
−2 1 1
𝑦 𝑥 = 𝑥2 𝑎𝑖 𝑥 𝑖 , 𝑦′ 𝑥 = 𝑥 + 𝑖 𝑎𝑖 𝑥 𝑖 , 𝑦 ′′ 𝑥 = 𝑥 − +𝑖 + 𝑖 𝑎𝑖 𝑥 𝑖
2 2 2
𝑖=0 𝑖=0 𝑖=0
reemplazamos en la ecuación de Bessel y obtenemos
∞ ∞ ∞
2
3
−2 1 1 1
−2 1 1
𝑥 𝑥 − +𝑖 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑥𝑥 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑥 2 − 2
𝑛 𝑥2 𝑎𝑖 𝑥 𝑖 = 0
2 2 2
𝑖=0 𝑖=0 𝑖=0
∞ ∞ ∞ ∞
1 1 1 𝑖
1 1 𝑖 2
1 1
𝑥2 − +𝑖 + 𝑖 𝑎𝑖 𝑥 + 𝑥 2 + 𝑖 𝑎𝑖 𝑥 − 𝑛 𝑥 2 𝑎𝑖 𝑥 + 𝑥 2 𝑥 2
𝑖 𝑎𝑖 𝑥 𝑖 = 0
2 2 2
𝑖=0 𝑖=0 𝑖=0 𝑖=0
∞ ∞
1 1 1
− +𝑖 +𝑖 + + 𝑖 − 𝑛2 𝑎𝑖 𝑥 𝑖 + 𝑎𝑖 𝑥 𝑖+2 = 0
2 2 2
𝑖=0 𝑖=0
∞ ∞
2
1 1 1
𝑖 − + + 𝑖 − 𝑎𝑖 𝑥 𝑖 + 𝑎𝑖−2 𝑥 𝑖 = 0
4 2 4
𝑖=0 𝑖=2
∞ ∞
1 1 1 1 1
− + − 𝑎0 + 𝑖2 − + 𝑖 + 𝑎𝑖 𝑥 𝑖 + 𝑎𝑖−2 𝑥 𝑖 = 0
4 2 4 4 4
𝑖=1 𝑖=2
∞ ∞
1 1
1− + 1 + 𝑎1 𝑥 + 𝑖 2 + 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑎𝑖−2 𝑥 𝑖 = 0
4 4
𝑖=2 𝑖=2
∞
2𝑎1 𝑥 + 𝑖 𝑖 + 1 𝑎𝑖 + 𝑎𝑖−2 𝑥 𝑖 = 0
𝑖=2
Para que la última igualdad se cumpla deben ser todos los coeficientes nulos, por lo tanto 𝑎0 es libre y los otros
coeficientes deben satisfacer:
𝑎1 = 0, 𝑖 + 1 𝑖 𝑎𝑖 + 𝑎𝑖−2 = 0 ∀𝑖 ≥ 2
Es decir,
𝑎𝑖−2
𝑎1 = 0, 𝑎𝑖 = − ∀𝑖 ≥ 2
𝑖+1 𝑖
𝑎0 𝑎0 𝑎2 𝑎0 𝑎0 𝑎4 𝑎0
⇒ 𝑎1 = 0, 𝑎2 = − = − , 𝑎3 = 0, 𝑎4 = − = = , 𝑎5 = 0, 𝑎6 = − = − ,⋯
3∙2 3! 5 ∙ 4 5 ∙ 4 ∙ 3 ∙ 2 5! 7∙6 7!
De manera más general,
𝑖
𝑎0
𝑎2𝑖+1 = 0 ∧ 𝑎2𝑖 = −1 ∀𝑖 ≥ 0
2𝑖 + 1!
Por lo tanto, la solución buscada es:
∞ ∞
1
𝑖
𝑎0 1
−2 𝑥 2𝑖+1 1
−2
𝑦 𝑥 = 𝑥2 −1 𝑥 2𝑖 = 𝑎0 𝑥 −1 𝑖
= 𝑎0 𝑥 𝑠𝑒𝑛(𝑥)
2𝑖 + 1! 2𝑖 + 1!
𝑖=0 𝑖=0
2 1
Si tomamos 𝑎0 = , la solución es conocida como la función de Bessel de orden 𝑛 = 2 y se denota por 𝐽1 , así
𝜋 2
2 1
𝐽1 𝑥 = 𝑥 2 𝑠𝑒𝑛(𝑥)
2 𝜋
1
De manera análoga se encuentra la solución asociada a 𝑟 = − 2.
EJEMPLO: ECUACIÓN DE LAGUERRE
Vamos a resolver la ecuación de Laguerre
𝑥𝑦 ′′ + 1 − 𝑥 𝑦 ′ + 𝑛𝑦 = 0
donde 𝑛 es un entero positivo o nulo, usando el método de Frobenius.
Observemos que la ecuación de Laguerre es equivalente a:
𝑥 2 𝑦 ′′ + 𝑥 1 − 𝑥 𝑦 ′ + 𝑛𝑥 𝑦 = 0
En este caso,
𝑏 𝑥 = 1 − 𝑥 ∧ 𝑐 𝑥 = 𝑛𝑥 ⇒ 𝑏0 = 1 ∧ 𝑐0 = 0
y la ecuación indical es:
𝑟 𝑟−1 +𝑟 =0
y tiene como única raíz a 𝑟 = 0 con multiplicidad 2. Tomamos 𝑟 = 0 entonces
∞ ∞ ∞
𝑦 𝑥 = 𝑎𝑖 𝑥 𝑖 , 𝑦′ 𝑥 = (𝑖 + 1) 𝑎𝑖+1 𝑥 𝑖 , 𝑦 ′′ 𝑥 = (𝑖 + 2)(𝑖 + 1)𝑎𝑖+2 𝑥 𝑖
𝑖=0 𝑖=0 𝑖=0
reemplazamos en la ecuación de Laguerre y obtenemos
∞ ∞ ∞
𝑥2 (𝑖 + 2)(𝑖 + 1)𝑎𝑖+2 𝑥 𝑖 + 𝑥 1 − 𝑥 (𝑖 + 1) 𝑎𝑖+1 𝑥 𝑖 + 𝑛𝑥 𝑎𝑖 𝑥 𝑖 = 0
𝑖=0 𝑖=0 𝑖=0
∞ ∞ ∞ ∞
(𝑖 + 2)(𝑖 + 1)𝑎𝑖+2 𝑥 𝑖+2 + (𝑖 + 1) 𝑎𝑖+1 𝑥 𝑖+1 − 𝑥 (𝑖 + 1) 𝑎𝑖+1 𝑥 𝑖+1 + 𝑛𝑥 𝑎𝑖 𝑥 𝑖 = 0
𝑖=0 𝑖=0 𝑖=0 𝑖=0
∞ ∞ ∞ ∞
𝑖(𝑖 − 1)𝑎𝑖 𝑥 𝑖 + 𝑖𝑎𝑖 𝑥 𝑖 − 𝑥 𝑖 𝑎𝑖 𝑥 𝑖 + 𝑛𝑥 𝑎𝑖 𝑥 𝑖 = 0
𝑖=2 𝑖=1 𝑖=1 𝑖=0
∞ ∞ ∞
𝑖(𝑖 − 1)𝑎𝑖 𝑥 𝑖 + 𝑖𝑎𝑖 𝑥 𝑖 + (𝑛 − 𝑖) 𝑎𝑖 𝑥 𝑖+1 + 𝑛𝑥𝑎0 = 0
𝑖=2 𝑖=1 𝑖=1
∞ ∞ ∞
𝑖(𝑖 − 1)𝑎𝑖 𝑥 𝑖 + 𝑎1 𝑥 + 𝑖𝑎𝑖 𝑥 𝑖 + (𝑛 − 𝑖 + 1) 𝑎𝑖−1 𝑥 𝑖 + 𝑛𝑥𝑎0 = 0
𝑖=2 𝑖=2 𝑖=2
∞ ∞ ∞
𝑖(𝑖 − 1)𝑎𝑖 𝑥 𝑖 + 𝑖𝑎𝑖 𝑥 𝑖 + (𝑛 − 𝑖 + 1) 𝑎𝑖−1 𝑥 𝑖 + 𝑥(𝑛𝑎0 + 𝑎1 ) = 0
𝑖=2 𝑖=2 𝑖=2
∞
𝑖 𝑖 − 1 𝑎𝑖 + 𝑖 𝑎𝑖 + (𝑛 − 𝑖 + 1) 𝑎𝑖−1 𝑥 𝑖 + 𝑥(𝑛𝑎0 + 𝑎1 ) = 0
𝑖=2
∞
𝑖 2 𝑎𝑖 + (𝑛 − 𝑖 + 1) 𝑎𝑖−1 𝑥 𝑖 + 𝑥(𝑛𝑎0 + 𝑎1 ) = 0
𝑖=2
Para que la igualdad se cumpla es necesario que todos los coeficientes de la serie se anulen. Por lo tanto, debe ser:
𝑛𝑎0 + 𝑎1 = 0, 𝑖 2 𝑎𝑖 + (𝑛 − 𝑖 + 1) 𝑎𝑖−1 = 0 ∀𝑖 ≥ 2
Es decir,
𝑎1 = −𝑛𝑎0
𝑛 − 𝑖 + 1 𝑎𝑖−1
𝑎𝑖 = − , 𝑖≥2
𝑖2
𝑛−1 𝑛−1 𝑛 𝑛−2 𝑛−2 𝑛−1 𝑛
𝑎1 = −𝑛𝑎0 , 𝑎2 = − 𝑎 = 𝑎 , 𝑎 = − 𝑎 = − 𝑎0 , ⋯
22 1 22 0 3
32 2 32 2 2
1 𝑛
2 ∙3⋯ 𝑛 − 2 𝑛 − 1 𝑛
𝑎𝑛 = − 2 𝑎𝑛−1 = −1 𝑎0 , 𝑎𝑛+1 = 0, 𝑎𝑛+2 = 0,
𝑛 𝑛 2 𝑛 − 1 2 ⋯ 32 22
De manera más general,
𝑖
(𝑛 − 𝑖 + 1)(𝑛 − 𝑖 + 2) ⋯ 𝑛 − 2 𝑛 − 1 𝑛
𝑎𝑖 = −1 𝑎0 ∀1 ≤ 𝑖 ≤ 𝑛 ∧ 𝑎𝑖 = 0 ∀𝑖 ≥ 𝑛 + 1
𝑖 2 𝑖 − 1 2 ⋯ 32 2 2
La solución buscada es el polinomio de grado 𝑛 definido por
𝑛 𝑛
𝑛!/ 𝑛 − 𝑖 ! 𝑖 𝑛 1 𝑖
𝑦 𝑥 = 𝑎0 + 𝑎0 −1 𝑖 𝑥 = 𝑎0 −1 𝑖 𝑥
𝑖! 𝑖! 𝑖 𝑖!
𝑖=1 𝑖=1
conocido como los polinomios de Laguerre.
