Métodos Matemáticos para la Física

Sonia Montecinos Geisse

8. Ecuaciones diferenciales. Muchos problemas físicos se plantean en la forma de una ecuación

diferencial, generalmente de varias variables.
Ejemplos:
 Dinámica
𝑑2 𝑥⃗
𝐹⃗ (𝑥⃗, 𝑡) = 𝑚
𝑑𝑡 2
 Ecuación de ondas
1 𝜕2
∇2 𝜓(𝑥⃗, 𝑡) − 𝜓(𝑥⃗, 𝑡) = 0
𝑐 2 𝜕𝑡 2
𝜕2 𝜕2 𝜕2
Donde ∇2 = 𝜕𝑥 2 + 𝜕𝑦2 + 𝜕𝑧2 es el Laplaciano.

 Mecánica Cuántica:
ℏ2 2 𝜕
− ∇ 𝜓(𝑥⃗, 𝑡) + 𝑉(𝑥⃗, 𝑡)𝜓(𝑥⃗, 𝑡) = 𝑖ℏ 𝜓(𝑥⃗, 𝑡)
2 𝜕𝑡
Ecuaciones diferenciales que dependen de una variable, se llaman ecuaciones diferenciales ordinarias,
y si dependen de más de una variable, ecuaciones diferenciales parciales.
8.1. Ecuaciones diferenciales parciales
En algunos casos, se pueden transformar en ecuaciones diferenciales ordinarias, usando el método
de
Separación de variables
Lo ilustraremos con un ejemplo. Supónganos que tenemos la ecuación
∇2 𝜓(𝑥⃗) + 𝑘 2 𝜓(𝑥⃗) = 0 (8.1)
Esta ecuación aparece en muchos problemas físicos.
i. Si 𝑘 2 = 0: ecuación de Laplace
ii. Si 𝑘 2 > 0: ecuación de Helmholtz (ecuación de onda)
iii. Si 𝑘 2 < 0: problemas de difusión
iv. 𝑘 2 = 𝐸 > 0; 𝐸 la energía del sistema: ecuación de Schrödinger de partícula libre
independiente del tiempo (Mecánica Cuántica).

Solución: Resolveremos la ecuación (8.1). El sistema de coordenadas que se usa para resolver el
problema, depende de las simetrías del problema, la cual está definida en las condiciones de borde
(CB).

1
 8.1.1. Coordenadas cartesianas
La ecuación (8.1) se escribe

𝜕 2 𝜓(𝑥⃗) 𝜕 2 𝜓(𝑥⃗) 𝜕 2 𝜓(𝑥⃗)

+ + + 𝑘 2 𝜓(𝑥⃗) = 0 (8.2)
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2
Postulamos una solución del tipo:

𝜓(𝑥⃗) = 𝑋(𝑥)𝑌(𝑦)𝑍(𝑧); 𝑥⃗ = (𝑥, 𝑦, 𝑧) (8.3)

Es directo ver que
∇2 𝜓(𝑥⃗) = 𝑋 ′′ (𝑥)𝑌(𝑦)𝑍(𝑧) + 𝑋(𝑥)𝑌 ′′ (𝑦)𝑍(𝑧) + 𝑋(𝑥)𝑌(𝑦)𝑍 ′′ (𝑧)

Reemplazamos en (8.2) y dividimos por 𝜓(𝑥⃗):

𝑋 ′′ 𝑌 ′′ 𝑍 ′′
+ + = −𝑘 2 (8.4)
𝑋 𝑌 𝑍
Como cada uno de los tres términos de la izquierda en (8.4) depende de variables distintas: La
ecuación se satisface sólo si cada uno de los términos es constante.
Llamamos a las constantes −𝑙 2 , −𝑚2 y −𝑛2, respectivamente. De esta forma hemos transformado
el problema de ecuaciones diferenciales parciales en un conjunto de tres ecuaciones diferenciales
ordinarias

𝑋 ′′ + 𝑙 2 𝑋 = 0; 𝑌′′ + 𝑚2 𝑌 = 0; 𝑍 ′′ + 𝑛2 𝑍 = 0 (8.5)
𝑙 + 𝑛 + 𝑚 = 𝑘2
2 2 2

Si 𝑙 2 , 𝑚2 y 𝑛2 son positivos (si 1, 𝑚, 𝑛 son complejos, los cuadrados pueden ser negativos o, en
general, números complejos), las soluciones son:

𝑋(𝑥) = 𝐴 sin 𝑙𝑥 + 𝐴1 cos 𝑙𝑥 ; 𝑌(𝑦) = 𝐵 sin 𝑚𝑦 + 𝐵1 cos 𝑚𝑦 ; 𝑍(𝑧) = 𝐶 sin 𝑛𝑧 + 𝐶1 cos 𝑛𝑧

El valor de las constantes y de 𝑙, 𝑚, 𝑛 depende de las CB.

En el caso particular de que la CB sea que la función se anula en los bordes de un cubo de lados a,
b, c, 𝑙 2 , 𝑚2 y 𝑛2 tienen que ser positivos.

Por ejemplo:
𝑋(𝑥) = 𝐴 sin 𝑙𝑥 + 𝐴1 cos 𝑙𝑥
Aplicamos las CB:

𝑋(𝑥 = 0) = 0 → 𝐴1 = 0
𝑋(𝑥 = 𝑎) = 0 = 𝐴 sin 𝑙𝑎
La única solución no nula es eligiendo la constante 𝑙 tal que se datisface:
𝑙𝑎 = 𝑛1 𝜋; 𝑛1 = 1,2, … ∞
Los 𝑛1 negativos no conducen a soluciones l.i. de la anterior.
Procedemos de manera similar en los otros ejes. Puede mostrarse que la solución general de la
ecuación (8.1) que satisfacen las CB es:

2
 𝑛1 𝜋 𝑛2 𝜋 𝑛3 𝜋
𝑋(𝑥) = 𝐴 sin 𝑥 ; 𝑌(𝑦) = 𝐵 sin 𝑦; 𝑍(𝑧) = 𝐶 sin 𝑧; 𝑛𝑖
𝑎 𝑏 𝑐
= 1,2, … ∞ (𝑖 = 1,2,3) (8.6)
La solución general se escribe

𝑛1 𝜋 𝑛2 𝜋 𝑛3 𝜋
𝜓(𝑥⃗) = 𝐴 sin ( 𝑥) sin ( 𝑦) sin ( 𝑧) ; 𝑛𝑖 = 1,2, … ∞ (8.7)
𝑎 𝑏 𝑐
Donde hemos renombrado la constante
𝐴𝐵𝐶 → 𝐴

8.1.2. Coordenadas esféricas

La ecuación (8.1) se escribe:

1 𝜕 2 𝜕𝜓 1 𝜕 𝜕𝜓 1 𝜕2𝜓
(𝑟 ) + (sin 𝜃 ) + + 𝑘 2 𝜓 = 0 (8.8)
𝑟 2 𝜕𝑟 𝜕𝑟 𝑟 2 sin 𝜃 𝜕𝜃 𝜕𝜃 𝑟 2 sin2 𝜃 𝜕𝜑2
Las coordenadas esféricas 𝑟, 𝜃, 𝜑 están definidas por las siguientes ecuaciones:
𝑥 = 𝑟 sin 𝜃 cos 𝜑 ; 𝑦 = 𝑟 sin 𝜃 sin 𝜑 ; 𝑧 = 𝑟 cos 𝜃 (8.9)

Para encontrar una solución de (8.8) en coordenadas esféricas procedemos de manera similar a lo
que hicimos en Coordenadas cartesianas. Postulamos:
𝜓(𝑟, 𝜃, 𝜑) = 𝑅(𝑟)𝛩(𝜃)𝛷(𝜑) (8.10)
2
Reemplazamos en (8.8) y multiplicamos por 𝑟 :

𝛩𝛷 𝑑 2 𝑑𝑅 𝑅𝛷 𝑑 𝑑𝛩 𝛩𝑅 𝑑2 𝛷
(𝑟 ) + (sin 𝜃 ) + + 𝑘 2 𝑅𝛩𝛷 = 0
𝑟 2 𝑑𝑟 𝑑𝑟 𝑟 2 sin 𝜃 𝑑𝜃 𝑑𝜃 𝑟 2 sin2 𝜃 𝑑𝜑2
Por simplicidad, hemos omitido el argumento de las funciones.

Después de un poco de ejercicio algebraico, se obtiene:

1 𝑑 2 𝑑𝑅 2 2
1 𝑑 𝑑𝛩 1 1 𝑑2 𝛷
{ (𝑟 ) + 𝑟 𝑘 +
} { (sin 𝜃 ) + }=0 (8.11)
𝑅 𝑑𝑟 𝑑𝑟 𝛩sin 𝜃 𝑑𝜃 𝑑𝜃 sin2 𝜃 𝛷 𝑑𝜑2

En el lado izquierdo de la ecuación tenemos dos términos que dependen de distintas variable. Para
que la ecuación se satisfaga debe tenerse:

1 𝑑 2 𝑑𝑅
(𝑟 ) + 𝑟 2 𝑘 2 = 𝑞; (8.12)
𝑅 𝑑𝑟 𝑑𝑟
1 1 𝑑 𝑑𝛩 1 1 𝑑2 𝛷
(sin 𝜃 ) + 2 = −𝑞 (8.13)
sin 𝜃 𝛩 𝑑𝜃 𝑑𝜃 sin 𝜃 𝛷 𝑑𝜑2
Con 𝑞 una constante.
Procedemos de manera similar en (8.13). Se encuentra
1 𝑑2 𝛷
= −𝑚2 ; ⇔ 𝛷′′ + 𝑚2 𝛷 = 0 (8.14)
𝛷 𝑑𝜑 2
Con 𝑚2 una constante.

Si 𝑚2 < 0; 𝑚 = 𝑖𝜅, 𝜅 real, la solución de (8.14) es de la forma:

𝛷 = 𝐴𝑒 𝜅𝜑 + 𝐵𝑒 −𝜅𝜑 ⇔ 𝛷 = 𝐴1 cosh 𝜅𝜑 + B1 sinh 𝜅𝜑 (8.15)

3
Ambas soluciones son equivalentes.
Si 𝑚2 > 0; (𝑚 real) la solución de (8.6) es de la forma:
𝛷 = 𝐴𝑒 𝑖𝑚𝜑 + 𝐵𝑒 −𝑖𝑚𝜑 ⇔ 𝛷 = 𝐴1 cos 𝑚𝜑 + B1 sin 𝑚𝜑 (8.16)

Debido a la definición del ángulo azimutal 𝜑 debe tenerse que 𝛷(𝜑 + 2𝜋) = 𝛷(𝜑). De esta
condición se deduce que 𝑚2 > 0. La condición de periodicidad se satisface 𝑚 = 0, ±1, ±2, ±3, ….
La solución general se escribe:

𝛷(𝜑) = 𝛷𝑚 (𝜑) = 𝐴𝑒 𝑖𝑚𝜑 ; = 0, ±1, ±2, ±3, … … (8.17)

Notemos que (8.4) + (8.17) indican que hay dos soluciones l.i.

Reemplazamos La ecuación (8.14) en (8.13). Se encuentra

1 𝑑 𝑑𝛩(𝜃) 𝑚2
(sin 𝜃 ) − 2 𝛩(𝜃) + 𝑞𝛩(𝜃) = 0 (8.18)
sin 𝜃 𝑑𝜃 𝑑𝜃 sin 𝜃

Propuesto. Plantear la ecuación (8.1) en coordenadas cilíndricas. Lugo usar separación de variables,
y plantear la solución de las funciones individuales.

Las ecuaciones (8.8) y (8.18) son ecuaciones diferenciales ordinarias. Sin embargo, la solución no
es directa. A continuación, veremos un método para resolver ecuaciones diferenciales ordinarias
(lineales) de segundo orden.

8.2. Ecuaciones diferenciales ordinarias

Consideremos una ecuación diferencial del tipo

𝑑2 𝑦 𝑑𝑦
+ 𝑃(𝑥) + 𝑄(𝑥)𝑦 = 0 (8.19)
𝑑𝑥 2 𝑑𝑥
Con CB
𝑦(𝑥 = 0) = 𝑦𝑜 ; 𝑦 ′ (𝑥 = 0) = 𝑦1 ; (8.19𝑎)

La ecuación (8.19) es una ecuación diferencial lineal, homogénea.

8.2.1. Puntos singulares y regulares

Se dice que 𝑥 = 𝑥𝑜 es un punto ordinario o regular de la ecuación diferencial (8.19) si 𝑃(𝑥) y 𝑄(𝑥)
son funciones regulares en 𝑥𝑜 .
Se dice que 𝑥 = 𝑥𝑜 es un punto singular si no es ordinario.
Un punto singular 𝑥𝑜 es singular regular, si es un punto singular y además las funciones
(𝑥 − 𝑥𝑜 )𝑃(𝑥); (𝑥 − 𝑥𝑜 )2 𝑄(𝑥)
son regulares en 𝑥 = 𝑥𝑜 .

4
Si 𝑥 = 𝑥𝑜 es un punto ordinario ó singular regular, en algunos casos puede encontrarse una solución
del tipo

𝑦 = ∑ 𝐴𝑛 (𝑥 − 𝑥𝑜 )𝑛 (8.20)
𝑛≥0

Es fácil ver que (hacerlo)

𝑦 ′ = ∑ 𝑛𝐴𝑛 (𝑥 − 𝑥𝑜 )𝑛−1 ; 𝑦 ′′ = ∑ 𝑛(𝑛 − 1)𝐴𝑛 (𝑥 − 𝑥𝑜 )𝑛−2 (8.21)

𝑛≥0 𝑛≥0

Notemos que hay algunos términos de la serie en (8.21) que se anulan. Notemos que hay dos
constantes arbitrarias que dependen de las CB. De (8.19a), encontramos
𝑦(𝑥 = 𝑥𝑜 ) = 𝐴𝑜 = 𝑦𝑜 ; 𝑦 ′ (𝑥 = 𝑥𝑜 ) = 𝐴1 = 𝑦1 ; (8.22)
Notar que (8.20) es una expansión en serie de Taylor de la solución en torno a 𝑥 = 𝑥𝑜 .
Ejemplo. Encontrar una solución de la ecuación de Airy.
𝑦 ′′ − 𝑥𝑦 = 0 (8.23)
Esta ecuación aparece en Mecánica Cuántica asociado a una partícula en presencia de una fuerza
constante (potencial lineal).
Postularemos una solución tipo (8.20) en torno a 𝑥 = 0. Mostrar primero que 𝑥 = 0 es un punto
ordinario ó singular-regular de (8.23).
Si comparamos (8.12) con (8.10) s directo ver que 𝑃(𝑥) = 0, y 𝑄(𝑥) = −𝑥. Ambas funciones son
regulares en 𝑥 = 0. Postulamos una solución tipo (8.20), con 𝑥𝑜 = 0:

𝑦 = ∑ 𝐴𝑛 𝑥 𝑛
𝑛≥0

𝑦 ′ = ∑ 𝑛𝐴𝑛 𝑥 𝑛−1 ; 𝑦 ′′ = ∑ 𝑛(𝑛 − 1)𝐴𝑛 𝑥 𝑛−2

𝑛≥0 𝑛≥0

Derivamos y reemplazamos. Usamos (8.21) para 𝑦 ′′ . Encontramos:

∑ 𝑛(𝑛 − 1)𝐴𝑛 𝑥 𝑛−2 − ∑ 𝐴𝑛 𝑥 𝑛+1 = 0

𝑛≥02 𝑛≥0

Ordenamos los índices. Vemos que la serie de la izquierda parte desde 𝑛 = 2

∑(𝑛 + 2)(𝑛 + 1)𝐴𝑛+2 𝑥 𝑛 − ∑ 𝐴𝑛−1 𝑥 𝑛 = 0

𝑛≥0 𝑛≥1

O, equivalentemente,

2𝐴2 + ∑((𝑛 + 1)(𝑛 + 2)𝐴𝑛+2 − 𝐴𝑛−1 )𝑥 𝑛 = 0

𝑛≥1

5
Como {𝑥 𝑛 }∞
𝑛=0 es un conjunto l.i., se obtiene:

1
𝐴2 = 0; 𝐴𝑛+2 = 𝐴 ; 𝑛≥1
(𝑛 + 1)(𝑛 + 2) 𝑛−1
El término de la derecha, se llama Relación de Recurrencia (RR).
En este caso, relaciona el valor de una constante de la serie con una constante anterior (3 ordenes
menos). Desarrollamos:
1 1 1
𝑛 = 1: 𝐴3 = 𝐴 ; 𝑛 = 2: 𝐴4 = 𝐴 ; 𝑛 = 3: 𝐴5 = 𝐴 =0
2∙3 0 3∙4 1 4∙5 2
Recordemos que 𝐴0 y 𝐴1 son constantes arbitrarias, cuyo valor depende de las CB (8.22).
Los siguientes términos
1 1 1 1
𝑛 = 4: 𝐴6 = 𝐴3 = 𝐴 𝑛 = 5: 𝐴7 = 𝐴4 = 𝐴 ; 𝐴8 = 0
5∙6 2∙3∙5∙6 0 6∙7 3∙4∙6∙7 1
Etcétera. En este ejemplo no es fácil encontrar un término general para los coeficientes de la serie.
Propuesto. Encontrar los 8 primeros términos de la solución de la ecuación de Airy para las siguientes
CB
a. 𝑦(0) = 1; 𝑦 ′ (0) = 0
b. 𝑦(0) = 0; 𝑦 ′ (0) = 1
Intente en esos casos una expresión general para los coeficientes. Por ejemplos, que coeficientes se
anulan?
Ejercicio. . Identifique las singularidades, en el caso de que existan. Verificar si el punto 𝑥 = 0
satisface las condiciones necesarias para encontrar una solución tipo (8.21) en torno a este punto.
Ejemplos de ecuaciones diferenciales de segundo orden
i. Ecuación de Bessel: 𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ + (𝑥 2 − 𝑛2 )𝑦 = 0
ii. Ecuación de Legendre (1 − 𝑥 2 )𝑦 ′′ − 2𝑥𝑦 ′ + 𝑙(𝑙 + 1)𝑦 = 0
iii. Ecuación de Laguerre 𝑥𝑦 ′′ + (1 − 𝑥)𝑦 ′ + 𝑎𝑦 = 0
iv. Oscilador armónico 𝑦 ′′ + 𝜔2 𝑦 = 0
v. Ecuación hypergeométrica (o de Gauss) 𝑥(𝑥 − 1)𝑦 ′′ + [(1 + 𝑎 + 𝑏)𝑥 − 𝑐]𝑦 ′ + 𝑎𝑏𝑦 = 0
(EH)
La ecuación hypergeométrica (EH) es una ecuación general, y muchas de las funciones que aparecen
en Mecánica Cuántica y en Física en general, son casos particulares de esta ecuación, con una elección
inteligente de los parámetros y variables.
Propuesto. Mostrar que
1−𝑥
𝑥→ ; 𝑎 = −1; 𝑏 = 𝑙 + 1; 𝑐 = 1
2
Transforma la EH en la ecuación de Legendre.
8.3. Método de Frobenius

6
Consideremos una ecuación diferencial del tipo (8.19). Sea 𝑥𝑜 un punto regular o singular regular de
la ecuación Diferencial.
El Método de Frobenius (MF) consiste en postular una solución del tipo

𝑦(𝑥) = (𝑥 − 𝑥𝑜 )𝑠 ∑ 𝐴𝑛 (𝑥 − 𝑥𝑜 )𝑛 = ∑ 𝐴𝑛 (𝑥 − 𝑥𝑜 )𝑛+𝑠 (8.24)

𝑛≥0 𝑛≥0

Derivamos

𝑦 ′ (𝑥) = ∑ 𝐴𝑛 (𝑛 + 𝑠)(𝑥 − 𝑥𝑜 )𝑛+𝑠−1 ; (8.25𝑎)

𝑛≥0

𝑦 ′′ (𝑥) = ∑ 𝐴𝑛 (𝑛 + 𝑠)(𝑛 + 𝑠 − 1)(𝑥 − 𝑥𝑜 )𝑛+𝑠−2 (8.25𝑏)

𝑛≥0

En lo que sigue consideraremos (sin pérdida de generalidad) que 𝑥𝑜 = 0.

Es importante recordar que existen dos soluciones l.i. de (8.19), 𝑦1 , 𝑦2 . La solución general se
escribe:
𝑦 = 𝐶1 𝑦1 + 𝐶2 𝑦2
Donde las constantes 𝐶1 y 𝐶2 dependen de las CB.
En el ejemplo que vimos antes (ecuación de Airy)
𝑦 ′′ − 𝑥𝑦 = 0
la solución es un caso particular de (8.25), con 𝑠 = 0.
Puede mostrarse que, si 𝑃(𝑥) y 𝑄(𝑥) son funciones analíticas en 𝑥 = 0 (o sea el punto 𝑥 = 0 es un
punto ordinario), entonces 𝑠 = 0.
Veremos a continuación algunos ejemplos
1.- Oscilador armónico.
𝑦 ′′ + 𝜔2 𝑦 = 0 (8.26)
La solución de esta ecuación es conocida, 𝑦1 = sin 𝜔𝑥; 𝑦2 = cos 𝜔𝑥, pero la buscaremos ahora con
el MF.
Postulamos una solución tipo (8.24) en torno a 𝑥𝑜 = 0 que, claramente, es un punto regular de (8.26)

𝑦(𝑥) = ∑ 𝐴𝑛 𝑥 𝑛+𝑠 (8.27)

𝑛≥0

Derivamos

𝑦 ′ (𝑥) = ∑(𝑛 + 𝑠)𝐴𝑛 𝑥 𝑛+𝑠−1 ; 𝑦 ′′ (𝑥) = ∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑥 𝑛+𝑠−2 (8.28)
𝑛≥0 𝑛≥0

Reemplazamos en (8.26) y factorizamos por 𝑥 𝑠 :

7
 𝑥 𝑠 (∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑥 𝑛−2 + 𝜔2 ∑ 𝐴𝑛 𝑥 𝑛 ) = 0
𝑛≥0 𝑛≥0

Hacemos en la serie de la izquierda el cambio de índice: 𝑛 → 𝑛 + 2.

𝑥 𝑠 ( ∑ (𝑛 + 2 + 𝑠)(𝑛 + 𝑠 + 1)𝐴𝑛+2 𝑥 𝑛 + 𝜔2 ∑ 𝐴𝑛 𝑥 𝑛 ) = 0
𝑛≥−2 𝑛≥0

Desarrollamos:

𝑥 𝑠 (𝑠(𝑠 − 1)𝐴𝑜 𝑥 −2 + (𝑠 + 1)𝑠𝐴1 𝑥 −1 + ∑((𝑛 + 2 + 𝑠)(𝑛 + 𝑠 + 1)𝐴𝑛+2 + 𝜔2 𝐴𝑛 ) 𝑥 𝑛 ) 0 = 0

𝑛≥0

Por la independencia de las funciones 𝑥 𝑛 , que incluye 𝑛 < 0 (ver serie de Laurent), encontramos las
ecuaciones:
𝑠(𝑠 − 1)𝐴𝑜 = 0; (𝑠 + 1)𝑠𝐴1 = 0. (8.28)
(𝑛 + 2 + 𝑠)(𝑛 + 𝑠 + 1)𝐴𝑛+2 + 𝜔2 𝐴𝑛 = 0; 𝑛 ≥ 0 (8.29)
La ecuación (8.28) se denomina Ecuación Indicial (EI). En este caso son dos ecuaciones que deben
satisfacerse simultáneamente.
De la ecuación (8.29) obtenemos la RR. En este caso
𝜔2
𝐴𝑛+2 = − 𝐴
(𝑛 + 2 + 𝑠)(𝑛 + 1 + 𝑠) 𝑛
Una solución posible de la EI es 𝑠 = 0, con 𝐴𝑜 y 𝐴1 constantes arbitrarias.
Reemplazamos en (8.29) y se encuentra:

𝜔2
𝐴𝑛+2 = − 𝐴 ; 𝑛≥0 (8.30)
(𝑛 + 2)(𝑛 + 1) 𝑛
La ecuación (8.30) conduce a dos soluciones l. i. Una serie con términos pares (los valores de 𝐴2𝑛 ),
y otra con impares (los valores de 𝐴2𝑛+1 ).
Primera solución 𝑦1 : 𝐴𝑜 arbitrario. De (8.30)
𝜔2 𝜔2 𝜔4
𝐴2 = − 𝐴 ; 𝐴4 = − 𝐴2 = 𝐴
2∙1 𝑜 4∙3 4! 𝑜
Usando inducción (hacerlo) puede mostrarse que
𝜔2𝑛
𝐴2𝑛 = (−1)𝑛 𝐴 ;
(2𝑛)! 𝑜
Demostrarlo. De donde,

8
 (−1)𝑛
𝑦1 (𝑥) = ∑ 𝐴2𝑛 𝑥 2𝑛 = 𝐴𝑜 ∑ (𝜔𝑥)2𝑛 = cos 𝜔𝑥
(2𝑛)!
𝑛≥0 𝑛≥0

(Hemos impuesto 𝐴𝑜 = 1). Esto no pierde generalidad.

Segunda solución. 𝑦2 : 𝐴1 arbitrario. Se encuentra:

𝜔2 𝜔2 𝜔4
𝐴3 = − 𝐴 ; 𝐴5 = − 𝐴3 = 𝐴
3∙2 1 5∙4 5! 1
𝜔2𝑛
𝐴2𝑛+1 = (−1)𝑛 𝐴
(2𝑛 + 1)! 1
𝐴1
Demostrarlo. Definiendo 𝐵 = ; se encuentra
𝜔

𝜔2𝑛+1 (−1)𝑛
𝐴2𝑛+1 = (−1)𝑛 𝐵; 𝑦2 (𝑥) = ∑ (𝜔𝑥)2𝑛+1 = sin 𝜔𝑥
(2𝑛 + 1)! (2𝑛 + 1)!
𝑛≥0

La solución general es:

𝑦(𝑥) = 𝐴𝑦1 (𝑥) + 𝐵𝑦2 (𝑥) = 𝐴 cos 𝜔𝑥 + 𝐵 sin 𝜔𝑥
Con la elección que hemos hecho en la EI, encontramos dos soluciones l.i. de (8.28) por lo tanto el
problema está solucionado. Si se encuentra una tercera solución, ésta necesariamente es
linealmente dependiente (l.d.), es decir, una combinación lineal de las otras dos.
En efecto, existen otras soluciones de la EI:
a. 𝑠 = 1; 𝐴1 = 0.
b. 𝑠 = −1; 𝐴𝑜 = 0.
Mostrar que a. conduce a la solución 𝑦2 y b. a la solución 𝑦1 .
2. Use el método de Frobenius para encontrar la solución de:
𝑥𝑦 ′′ + (1 − 𝑥)𝑦 ′ − 𝑦 = 0 (8.31)
En torno a 𝑥 = 0.
Solución. Primero hay que mostrar que 𝑥 = 0 es un punto regular, o singular regular de la ecuación
diferencial.
Identificamos 𝑃(𝑥) y 𝑄(𝑥).
1−𝑥 1
𝑃(𝑥) = ; 𝑄(𝑥) = −
𝑥 𝑥
Por lo tanto, 𝑥 = 0 es un punto singular regular.
Postulamos una solución del tipo (8.27). Usamos (8.28) y reemplazamos en (8.31)

𝑥 𝑠 (∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑥 𝑛−1 + ∑(𝑛 + 𝑠)𝐴𝑛 𝑥 𝑛−1 − ∑(𝑛 + 𝑠)𝐴𝑛 𝑥 𝑛 − ∑ 𝐴𝑛 𝑥 𝑛 ) = 0

𝑛≥0 𝑛≥0 𝑛≥0 𝑛≥0

9
Juntamos las dos primeras sumas ((𝑛 + 𝑠)(𝑛 + 𝑠 − 1) + (𝑛 + 𝑠) = (𝑛 + 𝑠)2 ).

Las dos primeras series se escriben:

∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑥 𝑛−1 + ∑(𝑛 + 𝑠)𝐴𝑛 𝑥 𝑛−1 = ∑(𝑛 + 𝑠)2 𝐴𝑛 𝑥 𝑛−1
𝑛≥0 𝑛≥0 𝑛≥0

∑(𝑛 + 𝑠)𝐴𝑛 𝑥 𝑛 + ∑ 𝐴𝑛 𝑥 𝑛 = ∑(𝑛 + 𝑠 + 1)𝐴𝑛 𝑥 𝑛

𝑛≥0 𝑛≥0 𝑛≥0

Juntamos términos y hacemos un cambio de índice. Encontramos:

𝑥 𝑠 ( ∑ (𝑛 + 1 + 𝑠)2 𝐴𝑛+1 𝑥 𝑛 − ∑(𝑛 + 𝑠 + 1)𝐴𝑛 𝑥 𝑛 ) = 0

𝑛≥−1 𝑛≥0

Juntamos términos de igual potencia:

𝑥 𝑠 (𝑠 2 𝐴𝑜 𝑥 −1 + ∑[(𝑛 + 1 + 𝑠)2 𝐴𝑛+1 − (𝑛 + 𝑠 + 1)𝐴𝑛 ]𝑥 𝑛 ) = 0

𝑛≥0

Igualamos a cero todos los términos que acompañan a 𝑥 𝑠+𝑛 .

Ecuación indicial:

𝑠 2 𝐴𝑜 = 0 (8.32𝑎)
Relación de Recurrencia:
1
𝐴𝑛+1 = 𝐴 ; 𝑛 ≥ 0; 𝑛 + 𝑠 + 1 ≠ 0 (8.32𝑏)
(𝑛 + 1 + 𝑠) 𝑛
Importante: La relación de recurrencia (RR) sirve para encontrar los términos de orden mayor de la
serie en término de los anteriores. Por eso, el término de orden mayor se pone a la izquierda
(recomendación).
La única solución de la EI (8.32a) es: 𝑠 = 0. Mostrar que la segunda solución (𝐴𝑜 = 0) conduce a la
solución nula.
Con esto la RR se escribe:
1
𝐴𝑛+1 = 𝐴 ; 𝑛≥0 (8.32𝑐)
𝑛+1 𝑛
Resolvemos
1 1
𝑛 = 0: 𝐴1 = 𝐴𝑜 ; 𝑛 = 1: 𝐴2 = 𝐴1 = 𝐴0
2 2
Término general (demostrarlo)
1
𝐴𝑛 = 𝐴 (8.32𝑑)
𝑛! 𝑜

10
Demostración. Usamos inducción. Ya está demostrado para 𝑛 = 0,1. Asumimos (8.32𝑐) es válido
para 𝑛. Debemos mostrar que es válida para 𝑛 + 1.
De (8.32𝑐) y (8.32d)
1 1 1 1
𝐴𝑛+1 = 𝐴𝑛 = 𝐴𝑜 = 𝐴
𝑛+1 (𝑛 + 1) 𝑛! (𝑛 + 1)! 𝑜
Queda demostrado.
Hemos mostrado que la (una) solución de (8.31) es:
1 𝑛
𝑦1 (𝑥) = ∑ 𝑥 = 𝑒𝑥
𝑛!
𝑛≥0

Queda propuesto mostrar que 𝑒 𝑥 es solución de (8.31).

En este caso, el método de Frobenius nos permite encontrar sólo una solución.
8.4. Una segunda solución
Recordemos que para que dos funciones 𝑦1 (𝑥) e 𝑦2 (𝑥) sean l.i. debe satisfacerse que la única solución
de
𝜆1 𝑦1 (𝑥) + 𝜆2 𝑦2 (𝑥) = 0
Es
𝜆1 = 𝜆2 = 0
La ecuación diferencial que queremos resolver es de la forma (ecuación (8.19)):
𝑦 ′′ + 𝑃(𝑥)𝑦 ′ + 𝑄(𝑥)𝑦 = 0
Si 𝑦1 (𝑥) e 𝑦2 (𝑥) son l.i. entonces se satisface que Wronskiano 𝑊 no se anula
𝑦1 𝑦2 ′ ′
𝑊 = |𝑦 ′ 𝑦2 ′ | = 𝑦1 𝑦2 − 𝑦1 𝑦2 ≠ 0 (8.33)
1

Ejemplo. Sea 𝑦1 (𝑥) = sin 𝑥 e 𝑦2 (𝑥) = cos 𝑥. Mostrar que: 𝑊 = −1, por lo tanto son l. i.
Derivando (8.33) se encuentra:
𝑊 ′ = 𝑦1 𝑦2 ′′ − 𝑦1 ′′ 𝑦2
Usando la ecuación diferencial (8.19) se encuentra
𝑊 ′ = −𝑃(𝑥)[𝑦1 𝑦2 ′ − 𝑦1 ′ 𝑦2 ] = −𝑃(𝑥)𝑊
O, equivalentemente
𝑊 ′ = −𝑃(𝑥)𝑊 (8.34)

En el caso particular 𝑃(𝑥) = 0 (por ejemplo oscilador armónico) la ecuación (8.34) toma la forma
simple 𝑊 ′ = 0, es decir

11
 𝑊 = 𝑦1 𝑦2 ′ − 𝑦1 ′ 𝑦2 = 𝑐𝑡𝑒 = 1 (8.35)
Como la ecuación diferencial que estamos resolviendo es homogénea, sin pérdida de generalidad
puede elegirse 𝑊 = 1. ¿Porqué?
Si 𝑃(𝑥) ≠ 0 la solución de la ecuación (8.34) (que también es válida cuando 𝑃(𝑥) = 0) se obtiene
integrando entre 𝑎 y 𝑥𝑜 . Se encuentra:
𝑥 𝑥
𝑑𝑊
∫ = − ∫ 𝑃(𝑞)𝑑𝑞
𝑎 𝑊 𝑎

Cuya solución es:

𝑥
ln 𝑊|𝑎𝑥 = − ∫ 𝑃(𝑞)𝑑𝑞
𝑎
𝑥
𝑊(𝑥) = 𝑊(𝑎) exp (− ∫ 𝑃(𝑞)𝑑𝑞) (8.36)
𝑎

Por otro lado, puede mostrarse que:

𝑑 𝑦2
𝑊(𝑥) = 𝑦1 𝑦2 ′ − 𝑦1 ′ 𝑦2 = 𝑦1 2 ( )
𝑑𝑥 𝑦1
Reemplazando en (8.35) y haciendo 𝑊(𝑎) = 1 (es importante convencerse que esto se puede hacer
sin pérdida de generalidad)
𝑥
𝑑 𝑦2
𝑦1 2 ( ) = exp (− ∫ 𝑃(𝑞)𝑑𝑞)
𝑑𝑥 𝑦1 𝑎
𝑥
𝑑 𝑦2 exp(− ∫𝑎 𝑃(𝑞)𝑑𝑞 )
( )=
𝑑𝑥 𝑦1 𝑦1 2 (𝑥)
Integramos entre b y x
𝑥 exp(− ∫𝑦 𝑃(𝑞)𝑑𝑞 )
𝑎
𝑦2 (𝑥) = 𝑦1 (𝑥) ∫ 𝑑𝑦 (8.37)
𝑏 𝑦1 2 (𝑦)

Notar que (8.36) nos sirve para encontrar una segunda solución, ¡aunque ya se conozca!
Por ejemplo, consideremos el oscilador armónico (ecuación (14)). En este caso se satisface 𝑃(𝑥) = 0
Sea
𝑦1 (𝑥) = cos 𝑥
una solución. Construiremos la segunda solución usando el método del Wronskiano. De (8.23),
tenemos que el numerador dentro de la integral es constate (= 1). Tenemos que
𝑥 𝑥
𝑑𝑦
𝑦2 (𝑥) = cos 𝑥 ∫ = cos 𝑥 ∫ 𝑠𝑒𝑐 2 (𝑦)𝑑𝑦 = cos 𝑥 (tan 𝑥 − tan 𝑏)
𝑏 𝑐𝑜𝑠 2 (𝑦) 𝑏

Por lo tanto

12
 𝑦2 (𝑥) = sin 𝑥 + 𝑐𝑡𝑒 cos 𝑥 : . 𝑦1 (𝑥) → 𝑦2 (𝑥) = sin 𝑥
Notar que el segundo término es de la forma: 𝑦2 (𝑥) = 𝑐𝑡𝑒 ∗ 𝑦1 (𝑥) y puede omitirse pues la solución
general es:
𝑦(𝑥) = 𝐴𝑦1 (𝑥) + 𝐵𝑦2 (𝑥)
Notar que en (8.23) no está explícitamente Q(x). Sin embargo, está implícitamente en la solución
𝑦1 (𝑥).

Retomemos ahora la ecuación (8.31).

𝑥𝑦 ′′ + (1 − 𝑥)𝑦 ′ − 𝑦 = 0
Habiamos encontrado que, usando el método de Frobenius (En torno a 𝑥 = 0) una solución es:
𝑦1 (𝑥) = 𝑒 𝑥

Para encontrar una segunda solución usando el método del Wronskiano identificamos 𝑃(𝑥) y 𝑄(𝑥).
𝑥 exp(− ∫𝑦 𝑃(𝑞)𝑑𝑞 )
𝑎
𝑦2 (𝑥) = 𝑦1 (𝑥) ∫ 2 (𝑦)
𝑑𝑦
𝑏 𝑦1

1−𝑥 1 1
𝑃(𝑥) = = − 1; 𝑄(𝑥) = −
𝑥 𝑥 𝑥
Construiremos ahora una segunda solución. Mostramos que
Calculamos
𝑦
𝑦
∫ 𝑃(𝑞)𝑑𝑞 = ln − (𝑦 − 𝑎) = ln 𝑦 − 𝑦 + 𝐶1
𝑎 𝑎

Con 𝐶1 una constante. De aquí

𝑦
𝑒𝑦
exp (− ∫ 𝑃(𝑞)𝑑𝑞 ) = 𝐶2
𝑎 𝑦

Escribir la constante 𝐶2 en función de 𝐶1 . Reemplazamos en (8.37)

𝑥
𝑒𝑦
𝑦2 (𝑥) = 𝑒 𝑥 ∫ 𝑑𝑦
𝑏 𝑦

Como estamos buscando una segunda solución, l.i. de 𝑦2 (𝑥), se puede elegir, sin pérdida de
generalidad, 𝐶2 = 1. Para resolver esta ecuación, escribimos el integrando en serie de Laurent. Se
tiene que

𝑒𝑦 1
= ∑ 𝑦 𝑛−1
𝑦 𝑛!
𝑛≥0

13
Reemplazamos en (24) e integramos término a término. Para integrar, separamos el primer término:
𝑒𝑦 1
= ∑ 𝑦 𝑛−1 =
𝑦 𝑛!
𝑛≥0

𝑒𝑦 1 1 1 1
= + ∑ 𝑦 𝑛−1 = + ∑ 𝑦𝑛
𝑦 𝑦 𝑛! 𝑦 (𝑛 + 1)!
𝑛≥1 𝑛≥0

Integramos

1
𝑦2 (𝑥) = 𝑒 𝑥 (ln 𝑥 − ln 𝑏 + ∑ 𝑥 𝑛+1 + 𝑐𝑡𝑒)
(𝑛 + 1)(𝑛 + 1)!
𝑛≥0

Notar que la constante aditiva al lado derecho corresponde a la suma de constantes que aparecen en
cada uno de los términos de integración.
Se encuentra

1
𝑦2 (𝑥) = 𝑒 𝑥 (ln 𝑥 + ∑ 𝑥 𝑛+1 ) + 𝐶′𝑦1 (𝑥) (8.38)
(𝑛 + 1)(𝑛 + 1)!
𝑛≥0

Sin pérdida de generalidad podemos elegir 𝑐𝑡𝑒 = 0. ¿Porqué?

Puede mostrarse que en general cuando la solución de la ecuación indicial es una raíz doble, se tiene
que

𝑦2 (𝑥) = 𝑦1 (𝑥) (ln 𝑥 + ∑ 𝐴𝑛 𝑥 𝑛 ) (8.39)

𝑛

La solución general es:

𝑦(𝑥) = 𝐴𝑦1 (𝑥) + 𝐵𝑦2 (𝑥)
Ejercicio.
Use el método de Frobenius para encontrar la (las) soluciones de la ecuación diferencial como una
serie en torno el punto 𝑥 = 0

𝑥 2 𝑦 ′′ − 2𝑥𝑦 ′ + 2(1 + 𝑥 2 )𝑦 = 0
Solución. El paso 1 es mostrar que 𝑥 = 0 es punto singular o singular regular. Hacerlo.
Postulamos una solución como en (8.28). Usamos (8.29) y reemplazamos

14
 𝑥 𝑠 (∑ 𝐴𝑛 (𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝑥 𝑛 − 2 ∑ 𝐴𝑛 (𝑛 + 𝑠)𝑥 𝑛 + 2 ∑ 𝐴𝑛 𝑥 𝑛 + 2 ∑ 𝐴𝑛 𝑥 𝑛+2 ) = 0
𝑛≥0 𝑛≥0 𝑛≥0 𝑛≥0

Para que esta igualdad se cumpla el término ( ) debe anularse. Ordenamos

∑ 𝐴𝑛 ((𝑛 + 𝑠)(𝑛 + 𝑠 − 3) + 2)𝑥 𝑛 + 2 ∑ 𝐴𝑛−2 𝑥 𝑛 = 0

𝑛≥0 𝑛≥2

Separamos los dos primeros términos de la serie de la izquierda.

𝐴𝑜 (𝑠(𝑠 − 3) + 2) + 𝐴1 ((𝑠 + 1)(𝑠 − 2) + 2)𝑥 + ∑(𝐴𝑛 ((𝑛 + 𝑠)(𝑛 + 𝑠 − 3) + 2) + 2𝐴𝑛−2 )𝑥 𝑛

𝑛≥2
=0
Los dos primeros términos corresponden a la EI. Los términos en la sumatoria dan origen a la RR.
Ecuación indicial:
𝐴𝑜 (𝑠(𝑠 − 3) + 2) = 0; 𝐴1 𝑠(𝑠 − 1) = 0 (𝐸𝐼)
2𝐴𝑛−2
𝐴𝑛 = − ; 𝑛≥2
(𝑛 + 𝑠)(𝑛 + 𝑠 − 3) + 2
Una solución de la EI (ambas) es: 𝑠 = 1.

Otras soluciones: 𝑠 = 0 y 𝐴𝑜 = 0; 𝑠 =? y 𝐴1 = 0.
Reemplazamos 𝑠 = 1 en la RR, y ordenando, se encuentra:
2𝐴𝑛−2
𝐴𝑛 = − ; 𝑛≥2 (𝑅𝑅)
𝑛(𝑛 − 1)
Ojo. Es importante verificar el que el denominador no se anula.
La RR nos permite encontrar dos soluciones l.i., una con términos pares de la serie y otra con términos
impares.
En este caso:
Se encuentras dos soluciones:

 𝐴𝑜 = 1; 𝐴1 = 0 → 𝑦1 (𝑥);
 𝐴𝑜 = 0; 𝐴1 = 1 → 𝑦2 (𝑥)
La solución general,
𝑦(𝑥) = 𝐴𝑦1 (𝑥) + 𝐵𝑦2 (𝑥)
Mostrar que si una condición de borde es 𝑦(0) = 0 entonces 𝐵 = 0.

Ejercicio. Use el método de Frobenius para resolver la ecuación diferencial

15
Expanda en torno a 𝑥 = 1.

𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ − 𝑦 = 0
Solución. Primero debemos verificar que el punto satisface las condiciones necesarias para que se
pueda usar el método de Frobenius. Identifica las funciones 𝑃(𝑥) y 𝑄(𝑥). Hacerlo.
Hacemos ahora el cambio de variables

𝑞 = 𝑥 − 1: 𝑥 2 = 𝑞 2 + 2𝑞 + 1
𝑑𝑦 𝑑𝑦 𝑑2 𝑦 𝑑2 𝑦
= : =
𝑑𝑥 𝑑𝑞 𝑑𝑥 2 𝑑𝑞 2
La ecuación diferencial se escribe

𝑞 2 𝑦 ′′ + 2𝑞𝑦 ′′ + 𝑦 ′′ + 𝑞𝑦 ′ + 𝑦 ′ − 𝑦 = 0: 𝑦 = 𝑦(𝑞)
Postulamos

𝑦(𝑞) = ∑ 𝐴𝑛 𝑞 𝑛+𝑠 (8.27)

𝑛≥0

Derivamos

𝑦 ′ (𝑞) = ∑(𝑛 + 𝑠)𝐴𝑛 𝑞𝑛+𝑠−1 ; 𝑦 ′′ (𝑞) = ∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑞𝑛+𝑠−2

𝑛≥0 𝑛≥0

Reemplazamos

∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑞𝑛+𝑠 + 2 ∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑞 𝑛+𝑠−1

𝑛≥0 𝑛≥0

+ ∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑞𝑛+𝑠−2 + ∑(𝑛 + 𝑠)𝐴𝑛 𝑞 𝑛+𝑠 + ∑(𝑛 + 𝑠)𝐴𝑛 𝑞𝑛+𝑠−1
𝑛≥0 𝑛≥0 𝑛≥0
𝑛+𝑠
− ∑ 𝐴𝑛 𝑞 =0
𝑛≥0

Factorizamos y ordenamos. El primer, cuarto y sexto término se juntan. Igualmente el segundo con
el quinto

𝑞 𝑠 (∑((𝑛 + 𝑠)2 − 1)𝐴𝑛 𝑞𝑛 + ∑(𝑛 + 𝑠)2 𝐴𝑛 𝑞𝑛−1 + ∑(𝑛 + 𝑠)(𝑛 + 𝑠 − 1)𝐴𝑛 𝑞𝑛−2 ) = 0
𝑛≥0 𝑛≥0 𝑛≥0

Hacemos ahora un cambio de índices y separamos los términos que acompañan a 𝑞 −2 y 𝑞 −1 .

16
𝑞 𝑠 (𝑠(𝑠 − 1)𝐴𝑜 𝑞 −2 + (𝑠(𝑠 + 1)𝐴1 + 𝑠 2 𝐴𝑜 )𝑞−1 + ∑(𝑛2 − 1)𝐴𝑛 𝑞𝑛 + ∑((𝑛 + 1 + 𝑠)2 )𝐴𝑛+1 𝑞𝑛
𝑛≥0 𝑛≥0

+ ∑(𝑛 + 2 + 𝑠)(𝑛 + 1 + 𝑠)𝐴𝑛+2 𝑞𝑛 ) = 0

𝑛≥0

Ecuación inidicial

𝑠(𝑠 − 1)𝐴𝑜 = 0: 𝑠(𝑠 + 1)𝐴1 + 𝑠 2 𝐴𝑜 = 0

Una solución es; 𝑠 = 0. 𝐴𝑜 y 𝐴1 arbitrrarios.
Relación de recurrencia (RR):
(𝑛 + 1)(𝑛 + 2)𝐴𝑛+2 + (𝑛 + 1)𝐴𝑛+1 + (𝑛2 − 1)𝐴𝑛 = 0
Notamos que:
(𝑛2 − 1) = (𝑛 − 1)(𝑛 + 1)
Como (𝑛 + 1) no se anula, se puede simplificar. La RR queda:
𝐴𝑛+1 + (𝑛 − 1)𝐴𝑛
𝐴𝑛+2 = . 𝑛≥0
𝑛+2
La solución se construye a partir de 𝐴𝑜 y 𝐴1 , que dependen de las CB.
𝐴1
𝑛 = 0: 𝐴2 =
2
𝐴2 +𝐴1 1
𝑛 = 1: 𝐴3 = : 𝐴3 = 𝐴1
3 2

Etcétera. Encontrar el término general. Demostrar su validez usando recurrencia. Notar que la
solución no depende de 𝐴𝑜 . Encuentre una segunda solución. Si no se puede encontrar una expresión
general, calcule los cinco primeros términos de la serie.

Ejercicio
Escriba las funciones 𝑃(𝑥) y 𝑄(𝑥) de la forma

𝑃(𝑥) = ∑ 𝑝𝑛 𝑥 𝑛 ; 𝑄(𝑥) = ∑ 𝑝𝑛 𝑥 𝑛
𝑛≥−1 𝑛≥−2

Luego escriba una solución de la ecuación diferencial de la forma

𝑦(𝑥) = ∑ 𝐴𝑛 𝑞𝑛+𝑠
𝑛≥0

Tome los primeros términos (cuatro términos) de la serie como primera solución.
Construya una segunda solución término por término, usando las expresiones vistas en esta sección

17
18