Tema 5

Ecuaciones diferenciales

Soluciones por medio de

series y funciones
especiales
 Índice

Esquema. . . . . . . . . . . . . . . . . . . . . . . 2

Ideas clave . . . . . . . . . . . . . . . . . . . . . . 3

5.1 Introducción y objetivos . . . . . . . . . . . . . 3

5.2 Puntos regulares . . . . . . . . . . . . . . . . 4

5.3 Puntos singulares y método de Frobenius . . . . . . 6

5.4 Ecuación de Legendre . . . . . . . . . . . . . . 8

5.5 Ecuación de Bessel . . . . . . . . . . . . . . . 13

5.6 Otras ecuaciones . . . . . . . . . . . . . . . . 19

5.7 Cuaderno de ejercicios . . . . . . . . . . . . . . 24

5.8 Referencias bibliográficas . . . . . . . . . . . . . 25

 Esquema

Soluciones en forma de serie para ecuaciones de segundo orden

con coeficientes variables en torno a:
- Puntos regulares
- Puntos singulares regulares (Método de Frobenius)

Ejemplos: Funciones especiales:

- Ecuación de Legendre - Polinomios de Legendre
- Ecuación de Bessel - Funciones de Bessel
- Polinomios de Hermite
- Polinomios de Laguerre

Otras ecuaciones notables

Ecuaciones diferenciales
2
Tema 5. Esquema
 Ideas clave

5.1 Introducción y objetivos

Hasta ahora hemos podido resolver sistemas de ecuaciones lineales y ecuaciones li‐
neales de orden 𝑛 siempre que sus coeficientes fueran constantes. El único caso en el
que pudimos resolver elementalmente una ecuación de segundo orden con coeficien‐
tes variables fue en el de las ecuaciones de Euler, y solamente porque un afortunado
cambio de variable podía transformarla en una ecuación con coeficientes constantes.

En este tema vamos a ver una forma de atacar ecuaciones lineales homogéneas de
segundo orden con coeficientes variables. Es decir, ecuaciones del tipo:

𝑦″ (𝑥) + 𝑎(𝑥)𝑦′ (𝑥) + 𝑏(𝑥)𝑦(𝑥) = 0 . (1)

Si 𝑎(𝑥) y 𝑏(𝑥) son analíticas en un cierto intervalo, podemos esperar que la solución
también lo sea. De hecho, un teorema nos dirá que se pueden encontrar dos solucio‐
nes linealmente independientes en forma de serie. El método para resolver la Ecua‐
ción (1) consiste pues en suponer una solución en forma de serie:

𝑦(𝑥) = Σ∞ 𝑘
𝑘=0 𝑐𝑘 (𝑥 − 𝑥0 ) ,

y llevar la solución a la Ecuación (1) para obtener los coeficientes 𝑐𝑘 de la serie. De ese
modo encontraremos una regla de recurrencia que nos dará todos los coeficientes
𝑐𝑘 en función de los dos primeros 𝑐0 y 𝑐1 , que son son las dos constantes arbitrarias
que aparecen en la solución de cualquier ecuación de segundo orden. Estas series
convergen al menos en el mismo intervalo que 𝑎(𝑥) y 𝑏(𝑥).

Para ciertos puntos en que 𝑎(𝑥) y 𝑏(𝑥) no son analíticas se puede aplicar el llamado
método de Frobenius. Esta técnica también nos permite obtener unas soluciones en
forma de series que daremos en un teorema sin demostración.

Ecuaciones diferenciales
3
Tema 5. Ideas clave
El cálculo de los coeficientes 𝑐𝑘 es muy tedioso. Sin embargo, ecuaciones del tipo Ecua‐
ción (1) aparecen en muchos problemas de física. Algunas de estas ecuaciones tienen
nombre propio, como la ecuaciones de Legendre, Bessel, Hermite, Laguerre, etc. que
definen una serie de funciones, llamadas funciones especiales. Las funciones espe‐
ciales tienen gran importancia en las matemáticas puras y aplicadas y existen libros
enteros dedicados a ellas (Berry, 2001). La definición y propiedades de las funciones
más comunes se pueden encontrar en cualquier libro de fórmulas y tablas, por ejem‐
plo en (M.R. Spigel, 2005). Este artículo (González & Manotas, 2016) contiene una in‐
teresante discusión sobre las funciones especiales que veremos en este tema, y otras
adicionales, incluyendo su origen, propiedades, contexto histórico, y aplicaciones.

Los objetivos de este tema son:

▸ Identificar los puntos regulares, singulares y singulares‐regulares en una ecua‐

ción lineal de segundo orden.

▸ Aprender a resolver ecuaciones lineales de segundo orden por medio de series,

cuando eso es posible.

▸ Reconocer las ecuaciones lineales más importantes en física y conocer algunas

propiedades de las funciones especiales que definen.

5.2 Puntos regulares

Consideremos la ecuación:

𝑦″ (𝑥) + 𝑎(𝑥)𝑦′ (𝑥) + 𝑏(𝑥)𝑦(𝑥) = 0 . (2)

En principio, una ecuación de este tipo no admite soluciones en forma de serie en

todo 𝑥 ∈ ℝ, sino que lo hace solamente aquellos intervalos en que las funciones 𝑎(𝑥)
y 𝑏(𝑥) son analíticas. Recordemos que una función 𝑓 (𝑥) es analítica en un punto 𝑥0
si se puede expresar como una serie de potencias (serie de Taylor) en torno a dicho

Ecuaciones diferenciales
4
Tema 5. Ideas clave
punto:
1 𝑑 𝑛𝑓
𝑓 (𝑥) = Σ∞
𝑛=0 (𝑥 − 𝑥0 )𝑛 .
𝑛
𝑛! 𝑑𝑥 𝑥|
0

Para ello la función debe ser derivable infinitas veces en 𝑥0 , y la serie debe ser con‐
vergente en un entorno de 𝑥0 de radio 𝑅, llamado radio de convergencia.

Definición 1: Puntos regulares y singulares

Un punto 𝑥0 es un punto regular de la Ecuación (2) si 𝑎(𝑥) y 𝑏(𝑥) son analíticas

en dicho punto. En el caso contrario se dice que 𝑥0 es un punto singular de la
Ecuación (2).

Si 𝑥0 es un punto regular de Ecuación (2) entonces 𝑎(𝑥) y 𝑏(𝑥) se pueden escribir como
series de potencias:

𝑎(𝑥) = Σ∞
𝑘=0 𝑎𝑘 (𝑥 − 𝑥0 )
𝑘

𝑏(𝑥) = Σ∞ 𝑘
𝑘=0 𝑏𝑘 (𝑥 − 𝑥0 ) ,

y es de esperar que cualquier solución de Ecuación (2) también tenga esta forma. De
hecho podemos enunciar el siguiente teorema:

Teorema 1

Si 𝑥0 es un punto regular de 𝑦″ (𝑥) + 𝑎(𝑥)𝑦′ (𝑥) + 𝑏(𝑥)𝑦(𝑥) = 0, su solución general

es:

𝑦(𝑥) = Σ∞
𝑘=0 𝑐𝑘 (𝑥 − 𝑥0 )
𝑘
(3)

= 𝑐0 𝑦1 (𝑥) + 𝑐1 𝑦2 (𝑥) ,

donde 𝑐0 y 𝑐1 son constantes arbitrarias y 𝑦1 (𝑥), 𝑦2 (𝑥) son soluciones de la ecua‐

ción que admiten desarrollos en serie convergentes al menos en el intervalo 𝑥 ∈
(−𝑅, 𝑅), donde 𝑅 es el menor de los dos radios de convergencia de 𝑎(𝑥) y 𝑏(𝑥).
Además, los coeficientes 𝑐𝑘 para 𝑘 ≥ 2 se pueden determinar de forma única en

Ecuaciones diferenciales
5
Tema 5. Ideas clave
 función de 𝑐0 y 𝑐1 llevando la serie Ecuación (3) a la Ecuación (2) con los coeficien‐
tes 𝑎(𝑥) y 𝑏(𝑥) desarrollados.

De esta manera obtenemos dos soluciones de la forma:

𝑦1 (𝑥) = 1 + Σ ,

𝑦2 (𝑥) = (𝑥 − 𝑥0 ) + Σ ,

donde los sumatorios Σ contienen las potencias de 𝑥 mayores o iguales que 2. Estas
soluciones satisfacen las condiciones iniciales 𝑦1 (𝑥0 ) = 1, 𝑦′1 (𝑥0 ) = 0, 𝑦2 (𝑥0 ) = 0,
𝑦′2 (𝑥0 ) = 1, por tanto su wronskiano en 𝑥0 es no‐nulo y las dos soluciones son lineal‐
mente independientes. Sabemos que la soluciones particulares existen y son únicas,
por ser a y b analíticas, y los valores iniciales 𝑦(𝑥0 ) = 𝑦0 , 𝑦′ (𝑥0 ) = 𝑦′0 determinan las
constantes 𝑐0 = 𝑦0 y 𝑐1 = 𝑦′0 .

5.3 Puntos singulares y método de Frobenius

Si 𝑥0 es un punto singular de la Ecuación (2), es decir, si 𝑎(𝑥) y 𝑏(𝑥) no son analíticas
en 𝑥0 , entonces no se puede aplicar el Teorema 1 para encontrar una solución en
dicho punto. Sin embargo en muchas ocasiones es precisamente en las proximidades
de los puntos singulares donde nos interesa conocer las soluciones de la ecuación.
Afortunadamente existe un método para obtener las soluciones de la Ecuación (2) en
una clase de puntos singulares llamados puntos singulares regulares.

Definición 2: Puntos singulares regulares

Un punto 𝑥0 es un punto singular regular de la Ecuación (2) si (𝑥 − 𝑥0 )𝑎(𝑥) y

(𝑥 − 𝑥0 )2 𝑏(𝑥) son analíticas en dicho punto.

Suponemos entonces que 𝑥0 es un punto singular de la Ecuación (2), y por simplicidad

suponemos que 𝑥0 = 0. Para puntos singulares regulares distintos de cero no hace

Ecuaciones diferenciales
6
Tema 5. Ideas clave
falta más que hacer el cambio de variable 𝑧 ≡ 𝑥 − 𝑥0 para llevar el punto singular
regular al origen. Si multiplicamos la Ecuación (2) por 𝑥2 :

𝑥2 𝑦″ (𝑥) + 𝑥2 𝑎(𝑥)𝑦′ (𝑥) + 𝑥2 𝑏(𝑥)𝑦(𝑥) = 0 ,

y definimos:

𝑎∗ (𝑥) ≡ 𝑥𝑎(𝑥) ,

𝑏∗ (𝑥) ≡ 𝑥2 𝑏(𝑥) ,

obtenemos la ecuación:

𝑥2 𝑦″ (𝑥) + 𝑥𝑎∗ (𝑥)𝑦′ (𝑥) + 𝑥2 𝑏∗ (𝑥)𝑦(𝑥) = 0 , (4)

donde 𝑎∗ (𝑥) y 𝑏∗ (𝑥) son analíticas en 𝑥 = 0:

𝑎∗ (𝑥) = Σ∞ ∗ 𝑘
𝑘=0 𝑎𝑘 𝑥

𝑏∗ (𝑥) = Σ∞ ∗ 𝑘
𝑘=0 𝑏𝑘 𝑥 .

Estas expresiones son válidas para 𝑥 ∈ (−𝑅, 𝑅), siendo 𝑅 el menor de los rangos de
convergencia de 𝑎∗ (𝑥) y 𝑏∗ (𝑥). Para establecer las soluciones de la Ecuación (4) en
forma de serie vamos a emplear el llamado método de Frobenius que expresaremos
en forma de teorema sin demostración:

Teorema 2: Método de Frobenius

La ecuación Ecuación (4) tiene siempre una solución 𝑦1 (𝑥) de la forma:

𝑦1 (𝑥) = 𝑥𝑟1 Σ∞ 𝑘
𝑘=0 𝑐𝑘 𝑥 , 𝑐0 ≠ 0 ,

donde 𝑟1 es la mayor de las dos raíces reales del polinomio indicial de Ecuación (4),
definido como:
𝑞(𝑟) = 𝑟2 + (𝑎∗0 − 1)𝑟 + 𝑏∗0 .

La segunda solución linealmente independiente 𝑦2 (𝑥) es, según los casos:

Ecuaciones diferenciales
7
Tema 5. Ideas clave
 Si 𝑟1 − 𝑟2 no es cero 𝑦2 (𝑥) = 𝑥𝑟2 Σ∞ 𝑘
𝑘=0 𝑏𝑘 𝑥 , 𝑏0 ≠ 0
ni entero positivo
Si 𝑟1 = 𝑟2 𝑦2 (𝑥) = 𝑥𝑟1 +1 Σ∞ 𝑘
𝑘=0 𝑏𝑘 𝑥 + 𝑦1 (𝑥) ln 𝑥
Si 𝑟1 − 𝑟2 = 𝑛 𝑦2 (𝑥) = 𝑥𝑟2 Σ∞ 𝑘
𝑘=0 𝑏𝑘 𝑥 + 𝑎𝑦1 (𝑥) ln 𝑥 , 𝑏0 ≠ 0, 𝑎 ∈ ℝ
entero positivo

Todas las soluciones están definidas al menos para 𝑥 ∈ (0, 𝑅) y los coeficientes
𝑐𝑘 ,𝑏𝑘 ,y la constante 𝑎 se pueden determinar sustituyendo cada una de las solu‐
ciones en la ecuación.

A partir de estas soluciones se obtienen otras que son válidas para 𝑥 ∈ (−𝑅, 0) sin
más que sustituir ln 𝑥 por ln |𝑥| y las expresiones del tipo 𝑥𝑟 por |𝑥|𝑟 .

A continuación, vamos a ver dos de las ecuaciones más relevantes que son resolu‐
bles por medio de series. Estas ecuaciones servirán como ejemplo ilustrativo de los
métodos descritos anteriormente.

5.4 Ecuación de Legendre

La ecuación de Legendre es:

(1 − 𝑥2 )𝑦″ − 2𝑥𝑦′ + 𝑝(𝑝 + 1)𝑦 = 0 , (5)

donde 𝑝 es un número entero. Esta ecuación aparece de forma natural al resolver la

ecuación de Laplace, que veremos en los temas sobre ecuaciones en derivadas par‐
ciales, mediante separación de variables en coordenadas esféricas. También aparece
al resolver la ecuación de Schrödinger para un potencial central en tres dimensiones,
por ejemplo, en el átomo de hidrógeno. La Ecuación (5) se puede escribir como:

2𝑥 𝑝(𝑝 + 1)
𝑦″ − 𝑦 ′
+ 𝑦 = 0,
(1 − 𝑥2 ) (1 − 𝑥2 )

Ecuaciones diferenciales
8
Tema 5. Ideas clave
en donde se ve claramente que 𝑥 = ±1 son puntos singulares regulares. La Ecua‐
ción (5) por tanto tiene soluciones analíticas al menos en el intervalo|𝑥| < 1. Vamos
a buscar una solución en torno al punto regular 𝑥 = 0 de la forma:

𝑦(𝑥) = Σ∞ 𝑘
𝑘=0 𝑐𝑘 𝑥 . (6)

Sus derivadas respecto a 𝑥 son:

𝑦′ (𝑥) = Σ∞
𝑘=1 𝑐𝑘 𝑘𝑥
𝑘−1
= Σ∞ 𝑘
𝑘=0 𝑐𝑘+1 (𝑘 + 1)𝑥 ,

𝑦″ (𝑥) = Σ∞
𝑘=2 𝑐𝑘 𝑘(𝑘 − 1)𝑥
𝑘−2
= Σ∞ 𝑘
𝑘=0 𝑐𝑘+2 (𝑘 + 2)(𝑘 + 1)𝑥 .

Introduciendo estas series en Ecuación (5) obtenemos:

Σ∞ 𝑘
𝑘=0 [𝑐𝑘+2 (𝑘 + 2)(𝑘 + 1)𝑥 − 𝑐𝑘+2 (𝑘 + 2)(𝑘 + 1)𝑥
𝑘+2
− 2𝑐𝑘+1 (𝑘 + 1)𝑥𝑘+1 + 𝑝(𝑝 + 1)𝑐𝑘 𝑥𝑘 ] = 0 .

A partir de esta ecuación tenemos que encontrar la regla de recurrencia para los co‐
eficientes 𝑐𝑘 . Para ello nos damos cuenta de que los términos en cada potencia de 𝑥
se deben anular. Si nos fijamos en el término en 𝑥𝑘 vemos que es:

𝑘
[𝑐𝑘+2 (𝑘 + 2)(𝑘 + 1) − 𝑐𝑘 𝑘(𝑘 − 1) − 2𝑐𝑘 𝑘 + 𝑝(𝑝 + 1)𝑐𝑘 ] 𝑥 = 0 .

Para que se cumpla la igualdad para todo 𝑥 el coeficiente debe ser cero, por tanto:

𝑐𝑘+2 (𝑘 + 2)(𝑘 + 1) − 𝑐𝑘 (𝑘(𝑘 + 1) − 𝑝(𝑝 + 1)) = 0 .

De donde obtenemos la regla de recurrencia:

𝑘(𝑘 + 1) − 𝑝(𝑝 + 1)
𝑐𝑘+2 = 𝑐𝑘 ,
(𝑘 + 2)(𝑘 + 1)

o expresada de otra manera que nos resultará más práctica:

(𝑝 − 𝑘 + 2)(𝑝 + 𝑘 − 1)
𝑐𝑘 = − 𝑐𝑘−2 ,
𝑘(𝑘 − 1)

Ecuaciones diferenciales
9
Tema 5. Ideas clave
que determina todos los coeficientes 𝑐𝑘 de la serie Ecuación (6) a partir de dos cons‐
tantes arbitrarias 𝑐0 y 𝑐1 :

𝑝(𝑝 + 1)
𝑐2 = − 𝑐 ,
2×1 0
(𝑝 − 1)(𝑝 + 2)
𝑐3 = − 𝑐1 ,
3×2
(𝑝 − 2)(𝑝 + 3) 𝑝(𝑝 − 2)(𝑝 + 1)(𝑝 + 3)
𝑐4 = − 𝑐2 = 𝑐0 ,
4×3 4×3×2×1
(𝑝 − 3)(𝑝 + 4) (𝑝 − 1)(𝑝 − 3)(𝑝 + 2)(𝑝 + 4)
𝑐5 = − 𝑐3 = 𝑐1 ,
5×4 5×4×3×2
⋯⋯

De estas expresiones podemos obtener la fórmula general para los coeficientes pares
e impares de la serie:

𝑝(𝑝 − 2) ⋯ (𝑝 − 2𝑛 + 2)(𝑝 + 1)(𝑝 + 3) ⋯ (𝑝 + 2𝑛 − 1)

𝑐2𝑛 = (−1)𝑛 𝑐0 ,
(2𝑛)!
(𝑝 − 1)(𝑝 − 3) ⋯ (𝑝 − 2𝑛 + 1)(𝑝 + 2)(𝑝 + 4) ⋯ (𝑝 + 2𝑛)
𝑐2𝑛+1 = (−1)𝑛 𝑐1 .
(2𝑛 + 1)!

Así que la solución general de la Ecuación (5) es:

𝑦(𝑥) = Σ∞
𝑘=0 𝑐𝑘 𝑥
𝑘

= 𝑐0 + 𝑐1 𝑥 + Σ∞ 2𝑛 ∞
𝑛=1 𝑐2𝑛 𝑥 + Σ𝑛=1 𝑐2𝑛+1 𝑥
2𝑛+1
,

esto es:
𝑦(𝑥) = 𝑐0 𝑦1 (𝑥) + 𝑐1 𝑦2 (𝑥) ,

donde:

𝑝(𝑝 − 2) ⋯ (𝑝 − 2𝑛 + 2)(𝑝 + 1)(𝑝 + 3) ⋯ (𝑝 + 2𝑛 − 1) 2𝑛

𝑦1 (𝑥) = 1 + Σ∞
𝑛=1 (−1)
𝑛
𝑥 ,
(2𝑛)!
𝑛 (𝑝 − 1)(𝑝 − 3) ⋯ (𝑝 − 2𝑛 + 1)(𝑝 + 2)(𝑝 + 4) ⋯ (𝑝 + 2𝑛) 2𝑛+1
𝑦2 (𝑥) = 𝑥 + Σ∞
𝑛=1 (−1) 𝑥 ,
(2𝑛 + 1)!

que convergen en |𝑥| < 1. Las series 𝑦1 (𝑥) y 𝑦2 (𝑥) se truncan para valores de 𝑝 pares
positivos e impares positivos respectivamente. Si 𝑝 = 2𝑚, 𝑦1 (𝑥) se convierte en un

Ecuaciones diferenciales
10
Tema 5. Ideas clave
polinomio de grado 2𝑚:

𝑝=0→ 𝑦1 (𝑥) = 1 ,

𝑝=2→ 𝑦1 (𝑥) = 1 − 3𝑥2 ,

35 4
𝑝=4→ 𝑦1 (𝑥) = 1 − 10𝑥2 + 𝑥 ,
3
⋯ ⋯

Mientras que si 𝑝 = 2𝑚 + 1, 𝑦2 (𝑥) se convierte en un polinomio de grado 2𝑚 + 1:

𝑝=1→ 𝑦2 (𝑥) = 𝑥 ,
5
𝑝=3→ 𝑦2 (𝑥) = 𝑥 − 𝑥3 ,
3
14 21
𝑝=5→ 𝑦2 (𝑥) = 𝑥 − 𝑥3 + 𝑥5 ,
3 5
⋯ ⋯

Estos polinomios 𝑃𝑛 (𝑥), debidamente normalizados de forma que 𝑃𝑛 (1) = 1, definen

los llamados polinomios de Legendre de grado 𝑛:

𝑃0 (𝑥) = 1 ,

𝑃1 (𝑥) = 𝑥 ,
3 1
𝑃2 (𝑥) = 𝑥2 − ,
2 2
5 3 3
𝑃3 (𝑥) = 𝑥 − 𝑥 ,
2 2
35 4 15 2 3
𝑃4 (𝑥) = 𝑥 − 𝑥 + ,
8 4 8
63 5 35 3 15
𝑃5 (𝑥) = 𝑥 − 𝑥 + 𝑥,
8 4 8
⋯

Los polinomios de Legendre tienen muchas propiedades interesantes y numerosas

aplicaciones, como la expansión en serie de un potencial de la forma 1/𝑟 y expansiones
multipolares.

El polinomio 𝑃𝑛 (𝑥) tiene 𝑛 ceros reales, todos ellos en el intervalo (−1, 1). Como 𝑃𝑛 (−𝑥) =
(−𝑥)𝑛 𝑃𝑛 (𝑥), los 𝑃𝑛 (𝑥) con 𝑛 par tienen simetría par, y aquellos con 𝑛 impar tienen sime‐

Ecuaciones diferenciales
11
Tema 5. Ideas clave
tría impar (ver la Figura 1). Los polinomios de Legendre se pueden definir de muchas
maneras, por ejemplo mediante la fórmula de Rodrigues:

1 𝑑𝑛 2
𝑃𝑛 (𝑥) = 𝑛 𝑛
(𝑥 − 1)𝑛 ,
2 𝑛! 𝑑𝑥

o mediante una fórmula de recurrencia:

(𝑛 + 1)𝑃𝑛+1 (𝑥) = (2𝑛 + 1)𝑥𝑃𝑛 (𝑥) − 𝑛𝑃𝑛−1 (𝑥) .

Los polinomios de Legendre forman un sistema completo de polinomios ortogonales.

El que son ortogonales significa que la integral:

1
𝑃𝑚 (𝑥)𝑃𝑛 (𝑥)𝑑𝑥
∫−1

2
es 0 cuando 𝑛 ≠ 𝑚. Si 𝑛 = 𝑚 entonces la integral es 2𝑛+1
. Utilizando la delta de Kro‐
necker podemos expresar esta propiedad como:

1
2
𝑃𝑚 (𝑥)𝑃𝑛 (𝑥)𝑑𝑥 = 𝛿 .
∫−1 2𝑛 + 1 𝑚𝑛

El que formen un sistema completo significa que cualquier función 𝑓 (𝑥) continua por
partes en el intervalo [−1, 1] se puede expandir en serie de polinomios de Legendre
de la forma:
𝑓 (𝑥) = Σ∞
𝑘=0 𝑎𝑘 𝑃𝑘 (𝑥) .

En el próximo tema profundizaremos sobre el significado de estas expansiones.

y
y

1 1

P0(x)
P3(x) P1(x)
0.5
P4(x) 0.5 P5(x)

0 0

P2(x)
-0.5 -0.5

-1 (a) (b)
-1

-1 -0.5 0 0.5 1 x x

Figura 1: Polinomios de Legendre. (a) 𝑃0 (𝑥), 𝑃2 (𝑥), 𝑃4 (𝑥); (b) 𝑃1 (𝑥), 𝑃3 (𝑥), 𝑃5 (𝑥).

Ecuaciones diferenciales
12
Tema 5. Ideas clave
5.5 Ecuación de Bessel
La ecuación de Bessel es:

𝑥2 𝑦″ + 𝑥𝑦′ + (𝑥2 − 𝑝2 )𝑦 = 0 . 𝑝 ≥ 0. (7)

Esta ecuación también aparece al resolver la ecuación de Laplace en coordenadas cilín‐

dricas y la ecuación de Helmholtz en coordenadas esféricas (La ecuación de Helmholtz
es una forma independiente del tiempo de la ecuación de ondas que veremos más
adelante). Reescribiendo la ecuación de Bessel en la forma:

1 𝑥2 − 𝑝2
𝑦″ + 𝑦′ + 𝑦 = 0,
𝑥 𝑥2

vemos claramente que 𝑥 = 0 es un punto singular regular de la Ecuación (7). Además,

comparando la Ecuación (7) con Ecuación (4) vemos que 𝑎∗ (𝑥) = 1 y 𝑏∗ (𝑥) = 𝑥2 − 𝑝2 ,
y el primer término de su desarrollo es 𝑎∗0 = 1 y 𝑏∗0 = −𝑝2 . Podemos aplicar el método
de Frobenius y proponer una solución de la forma:

𝑦1 (𝑥) = 𝑥𝑟1 Σ∞ 𝑘
𝑘=0 𝑐𝑘 𝑥 ,

donde 𝑟1 es la mayor de las raíces del polinomio indicial, que en este caso es:

𝑞(𝑟) = 𝑟(𝑟 − 1) + 𝑟 − 𝑝2

= 𝑟 2 − 𝑝2 ,

con raíces 𝑟 = ±𝑝. Por tanto, debemos probar la solución:

𝑦1 (𝑥) = Σ∞
𝑘=0 𝑐𝑘 𝑥
𝑝+𝑘
, (8)

que converge en todo ℝ. Sus derivadas son:

𝑦′1 (𝑥) = Σ∞
𝑘=0 𝑐𝑘 (𝑝 + 𝑘)𝑥
𝑝+𝑘−1
,

𝑦″1 (𝑥) = Σ∞
𝑘=0 𝑐𝑘 (𝑝 + 𝑘)(𝑝 + 𝑘 − 1)𝑥
𝑝+𝑘−2
.

Ecuaciones diferenciales
13
Tema 5. Ideas clave
Llevando estas expresiones a la Ecuación (7) obtenemos:

Σ∞
𝑘=0 [𝑐𝑘 (𝑝 + 𝑘)(𝑝 + 𝑘 − 1)𝑥
𝑝+𝑘
+ 𝑐𝑘 (𝑝 + 𝑘)𝑥𝑝+𝑘 + 𝑐𝑘 𝑥𝑝+𝑘+2 − 𝑝2 𝑐𝑘 𝑥𝑝+𝑘 ] = 0 ⇒

Σ∞ 𝑝2 + 𝑝𝑘 − 𝑝
+ 𝑝𝑘 + 𝑘2 − 
𝑘=0 [𝑐𝑘 ( 𝑘 + 𝑝
+  𝑝2 )𝑥𝑝+𝑘 + 𝑐𝑘 𝑥𝑝+𝑘+2 ] = 0 ⇒
𝑘−

Σ∞
𝑘=0 [𝑐𝑘 (2𝑝 + 𝑘)𝑘𝑥
𝑝+𝑘
+ 𝑐𝑘 𝑥𝑝+𝑘+2 ] = 0 . (9)

Como los términos en cada potencia de 𝑥 se deben anular, tenemos que:

𝑐𝑘 (2𝑝 + 𝑘)𝑘 + 𝑐𝑘−2 = 0 ,

lo cual nos da la regla de recurrencia:

1
𝑐𝑘 = − 𝑐 .
(2𝑝 + 𝑘)𝑘 𝑘−2

Además, si escribimos los primeros términos de la Ecuación (9):

𝑐0 𝑥𝑝+2 + (2𝑝 + 1)𝑐1 𝑥𝑝+1 + 𝑐1 𝑥𝑝+3 + 2(2𝑝 + 2)𝑐2 𝑥𝑝+2 + 𝑐2 𝑥𝑝+4 + ⋯ = 0 ,

vemos que la única manera de que el término en 𝑥𝑝+1 se anule es haciendo 𝑐1 = 0.

Por tanto, todos los coeficientes 𝑐𝑘 con 𝑘 impar son cero, mientras que los coeficientes
con 𝑘 par son:

1
𝑐2 = − 𝑐 ,
2 × 2(𝑝 + 1) 0
1 1
𝑐4 = − 𝑐 =− 4 𝑐 ,
4 × 2(𝑝 + 2) 2 2 × 2(𝑝 + 2)(𝑝 + 1) 0
⋯

es decir:

𝑐2𝑛+1 = 0 ,
(−1)𝑛
𝑐2𝑛 = 𝑐 .
22𝑛 𝑛!(𝑝 + 𝑛)(𝑝 + 𝑛 − 1) ⋯ (𝑝 + 1) 0

Ecuaciones diferenciales
14
Tema 5. Ideas clave
Así que la primera solución de la ecuación de Bessel es:

(−1)𝑛
𝑦1 (𝑥) = 𝑐0 𝑥𝑝 1 + Σ∞
𝑛=1 2𝑛 𝑥2𝑛 . (10)
( 2 𝑛!(𝑝 + 𝑛)(𝑝 + 𝑛 − 1) ⋯ (𝑝 + 1) )

Esta solución suele expresarse en términos de la función gamma, que se define como:

Definición 3: Función gamma

La función gamma viene definida:

∞
Γ(𝑠) = 𝑥𝑠−1 𝑒−𝑥 𝑑𝑥 , 𝑥 > 0.
∫0

La función gamma tiene la siguiente importante propiedad:

Γ(𝑠 + 1) = 𝑠Γ(𝑠) ,

por tanto:
Γ(𝑠 + 𝑛 + 1) = (𝑠 + 𝑛)(𝑠 + 𝑛 − 1) ⋯ (𝑠 + 1)𝑠Γ(𝑠) ,

y si 𝑠 ∈ ℕ, entonces Γ(𝑠 + 1) = 𝑠!

Utilizando la función gamma, la solución Ecuación (10) se puede escribir como:

𝑛
𝑝 ∞ (−1) Γ(𝑝 + 1) 𝑥 2𝑛
𝑦1 (𝑥) = 𝑐0 𝑥 Σ𝑛=0 (2) .
𝑛!Γ(𝑝 + 𝑛 + 1)

1
Esta solución, con 𝑐0 = 2𝑝 Γ(𝑝+1)
, define la llamada función de Bessel de primera especie
y orden 𝑝:

Definición 4: Función de Bessel de primera 1a y orden 𝑝

𝑥 𝑝 (−1)𝑛 𝑥 2𝑛
𝐽𝑝 (𝑥) = ( ) Σ∞
𝑛=0 𝑛!Γ(𝑝 + 𝑛 + 1) ( 2 ) .
2

Como ejemplo, veamos la forma de las funciones de Bessel de primera especie con

Ecuaciones diferenciales
15
Tema 5. Ideas clave
𝑝 = 0 y 𝑝 = 1:

(−1)𝑛 𝑥 2𝑛
𝐽0 (𝑥) = Σ∞ ,
𝑛=0
(𝑛!)2 ( 2 )
(−1)𝑛 𝑥 2𝑛+1
𝐽1 (𝑥) = Σ∞
𝑛=0 𝑛!(𝑛 + 1)! ( 2 ) .

Al igual que 𝐽0 y 𝐽1 todas la 𝐽𝑝 son funciones oscilatorias (Figura 2), de hecho, para 𝑥
suficientemente grande se tiene que:

2 𝜋
𝐽𝑝 ∼ √ cos(𝑥 − (2𝑝 + 1) ) . (11)
𝜋𝑥 4

Cada 𝐽𝑝 tienen un número infinito de ceros en (0, ∞) que se pueden encontrar tabu‐
lados en cualquier libro de funciones especiales. Las funciones de Bessel son especial‐
mente importantes en la propagación de ondas y en los problemas con potenciales
estáticos. Al resolver problemas en sistemas de coordenadas cilíndricas se obtienen
funciones de Bessel de orden entero (𝑝 = 𝑛); en coordenadas esféricas se obtienen
ordenes semi‐enteros (𝑝 = 𝑛 + 1/2).

Las funciones de Bessel también se pueden definir de forma recursiva mediante la

fórmula:
2𝑝
𝐽𝑝−1 (𝑥) + 𝐽𝑝+1 (𝑥) = 𝐽 (𝑥) .
𝑥 𝑝
Para valores no‐enteros de 𝑝 las funciones 𝐽𝑝 (𝑥) y 𝐽−𝑝 (𝑥) son linealmente indepen‐
dientes. Por otro lado, si 𝑝 es entero (𝑝 = 𝑛) se cumple que:

𝐽−𝑛 (𝑥) = (−1)𝑛 𝐽𝑛 (𝑥) .

Además, las funciones de Bessel con 𝑝 semi‐entero se pueden expresar de manera

exacta por medio de senos y cosenos:

2 2
𝐽 1 (𝑥) = √ sin 𝑥 𝐽− 1 (𝑥) = √ cos 𝑥
2 𝜋𝑥 2 𝜋𝑥
2 1 2 1
𝐽 3 (𝑥) = √ ( sin 𝑥 − cos 𝑥) 𝐽− 3 (𝑥) = −√ (sin 𝑥 + cos 𝑥)
2 𝜋𝑥 𝑥 2 𝜋𝑥 𝑥
⋯ ⋯

Ecuaciones diferenciales
16
Tema 5. Ideas clave
Ahora vamos a tratar de encontrar una segunda solución de Ecuación (7)linealmente
independiente de la primera. El Teorema 2 nos da tres posibles soluciones dependien‐
do de si la diferencia entre las raíces del polinomio indicial 𝑟1 − 𝑟2 es entera, no entera,
o cero. En este caso, 𝑟1 − 𝑟2 = 2𝑝. Veamos las diferentes posibilidades:

▸ Si 2𝑝 no es cero ni entero. La segunda solución a la ecuación de Bessel es de la

forma:
𝑦2 (𝑥) = 𝑥−𝑝 Σ∞ 𝑘
𝑘=0 𝑏𝑘 𝑥 , 𝑥 > 0,

que es idéntica a la Ecuación (8) excepto por el signo de 𝑝. Por tanto, para en‐
contrar los coeficientes 𝑏𝑘 no hay más que repetir los cálculos anteriores susti‐
tuyendo 𝑝 por −𝑝. Así llegamos a la solución 𝑦2 (𝑥) = 𝐽−𝑝 (𝑥), con:

𝑥 −𝑝 (−1)𝑛 𝑥 2𝑛
𝐽−𝑝 (𝑥) = ( ) Σ∞𝑛=0 𝑛!Γ(−𝑝 + 𝑛 + 1) ( 2 ) .
2

▸ Si 𝑝 es cero. Entonces el Teorema 2 nos propone la siguiente solución:

𝑦2 (𝑥) = Σ∞
𝑘=0 𝑏𝑘 𝑥
𝑘+1
+ 𝐽0 (𝑥) ln 𝑥 . 𝑥>0

Introduciendo esta solución en la Ecuación (7) se determinan los coeficientes

𝑏𝑘 , dando lugar a:

(−1)𝑛+1 1 1 𝑥 2𝑛
𝑦2 (𝑥) = Σ∞ ( 1 + + ⋯ + ) ( + 𝐽0 (𝑥) ln 𝑥 , 𝑥 > 0,
𝑛=1
(𝑛!)2 2 𝑛 2)

la cual define la función de Bessel de segunda especie y orden 0:

Definición 5: Función de Bessel de segunda especie y orden 0

Esta función se expresa matemáticamente como:

(−1)𝑛+1 1 1 𝑥 2𝑛
𝑌0 (𝑥) = Σ∞ 1 + + ⋯ + + 𝐽0 (𝑥) ln 𝑥 , 𝑥 > 0.
𝑛=1
(𝑛!)2 ( 2 𝑛) (2)

▸ Si 2𝑝 es entero. Entonces la segunda solución a la ecuación de Bessel es:

𝑦2 (𝑥) = 𝑥−𝑝 Σ∞ 𝑘
𝑘=0 𝑏𝑘 𝑥 + 𝑎𝐽𝑝 (𝑥) ln 𝑥 , 𝑥 > 0,

Ecuaciones diferenciales
17
Tema 5. Ideas clave
 1 3 5
Si 2𝑝 es impar (𝑝 = , , , … ),
2 2 2
se puede demostrar que 𝑎 = 0 y la solución
tiene la misma forma que en el caso a), es decir, 𝑦2 (𝑥) = 𝐽−𝑝 (𝑥). Si 2𝑝 es par
(𝑝 = 1, 2, 3, … ) la solución es una función de Bessel de segunda especie y orden
𝑝, definida como:

Definición 6: Función de Bessel de segunda especie y orden 𝑝 (𝑝 ∈ ℕ)

1 𝑥 −𝑝 1 1 1 𝑥 2𝑝 𝑝−1 (𝑝 − 𝑛 − 1)! 𝑥 2𝑛
𝑌𝑝 (𝑥) = − 1 + + ⋯ + + Σ
2 ( 2 ) ( 𝑝! ( 2 𝑝) (2) 𝑛=0 𝑛! (2) )
1 𝑥 𝑝 (−1)𝑛+1 1 1 1 1 𝑥 2𝑛
− ( ) Σ∞ 𝑛=0 𝑛!(𝑛 − 𝑝)! (( 1 + + ⋯ + + 1 + + ⋯ +
2 2 2 𝑛) ( 2 𝑛 + 𝑝 )) ( 2 )
+ 𝐽𝑝 (𝑥) ln 𝑥 , 𝑥>0
y
y

1 1

J0(x) Y0(x)
0.5 Y1(x)
J1(x)
0.5

-0.5

x x
-0.5 -1
0 5 10 15 20 25 0 5 10 15 20 25

Figura 2: Funciones de Bessel de primera especie y funciones de Bessel de segunda

especie. Elaboración propia.

Las funciones de Bessel de segunda especie están relacionadas con las de pri‐
mera especie mediante la siguiente fórmula:

𝐽𝑝 (𝑥) cos(𝑝𝜋) − 𝐽− 𝑝(𝑥)

𝑌𝑝 (𝑥) = .
sin(𝑝𝜋)

En el caso de orden entero (𝑝 = 𝑛) las funciones se definen tomando el límite:

𝑌𝑛 (𝑥) = lím 𝑌𝑝 (𝑥) ,

𝑝→𝑛

Ecuaciones diferenciales
18
Tema 5. Ideas clave
 de lo cual se puede deducir que cuando 𝑝 es entero (𝑝 = 𝑛) se cumple que:

𝑌−𝑛 (𝑥) = (−1)𝑛 𝑌𝑛 (𝑥) .

5.6 Otras ecuaciones

A continuación, vamos a ver brevemente otras ecuaciones con nombre propio que
tienen relevancia en física y que definen otras series de funciones especiales. En es‐
te vídeo demostramos algunas de las propiedades de las funciones especiales más
importantes.

Ŏ Accede al vídeo: Funciones especiales.

Ecuación general de Legendre

La ecuación general de Legendre es:

2 ″ ′ 𝑚2
(1 − 𝑥 ) 𝑦 − 2𝑥𝑦 + [ℓ(ℓ + 1) − 𝑦 = 0, (12)
1 − 𝑥2 ]

cuyas soluciones son las funciones asociados de Legendre 𝑃ℓ𝑚 (𝑥). Los indices ℓ y 𝑚
(que son enteros) se llaman grado y orden respectivamente de la función asociada
de Legendre. La Ecuación (12) tiene soluciones no‐nulas y no‐singulares en [−1, 1]
solo si ℓ y 𝑚 son enteros no‐negativos con 0 ≤ 𝑚 ≤ ℓ o alguno de ellos es negativo
con 0 ≤ |𝑚| ≤ |ℓ|. Además, cuando 𝑚 es par, la función 𝑃ℓ𝑚 (𝑥) es un polinomio.
Cuando 𝑚 = 0 y ℓ es un entero, las funciones 𝑃ℓ𝑚 (𝑥) son idénticas a los polinomios de
Legendre.

Las funciones asociadas de Legendre forman parte de la definición de otras funciones

especiales, los esféricos armónicos, que tienen un papel fundamental en el estudio del
momento angular en mecánica cuántica.

Ecuaciones diferenciales
19
Tema 5. Ideas clave
 ▸ Definición para ℓ y 𝑚 enteros no‐negativos. Las funciones 𝑃ℓ𝑚 (𝑥), están relacio‐
nadas con los polinomios de Legendre mediante la fórmula:

𝑑𝑚
𝑃ℓ𝑚 (𝑥) = (−1)𝑚 (1 − 𝑥2 )𝑚/2 𝑃 (𝑥) ,
𝑑𝑥𝑚 ℓ

Las funciones asociadas de Legendre no son mutuamente ortogonales en gene‐

ral. Por ejemplo, 𝑃11 no es ortogonal a 𝑃22 . Sin embargo, algunos sub‐conjuntos
sí son ortogonales. Suponiendo que 0 ≤ 𝑚 ≤ ℓ , para 𝑚 fijo las funciones satis‐
facen la condición de ortogonalidad:

1
2(ℓ + 𝑚)!
𝑃𝑘𝑚 𝑃ℓ𝑚 𝑑𝑥 = 𝛿 .
∫−1 (2ℓ + 1)(ℓ − 𝑚)! 𝑘ℓ

▸ Definición para 𝑚 negativo y/o ℓ negativo. La Ecuación (12) es claramente inva‐

riante bajo un cambio de signo de 𝑚. Las funciones con 𝑚 negativo son propor‐
cionales a la de 𝑚 positivo:

(ℓ − 𝑚)! 𝑚
𝑃ℓ−𝑚 = (−1)𝑚 𝑃 .
(ℓ + 𝑚)! ℓ

A partir de su definición, se puede comprobar que las funciones asociadas de

Legendre son pares o impares de acuerdo a la relación:

𝑃ℓ𝑚 (−𝑥) = (−1)ℓ+𝑚 𝑃ℓ𝑚 (𝑥) ,

y cumplen la fórmula de recurrencia:

𝑚
(ℓ − 𝑚 + 1)𝑃ℓ+1 (𝑥) = (2ℓ + 1)𝑥𝑃ℓ𝑚 (𝑥) − (ℓ + 𝑚)𝑃ℓ−1
𝑚
(𝑥) .

Ecuación de Hermite
La ecuación de Hermite de orden 𝜆 es:

𝑦″ − 2𝑥𝑦′ + 2𝜆𝑦 = 0 ,

Ecuaciones diferenciales
20
Tema 5. Ideas clave
y aparece en el problema del oscilador armónico cuántico entre otros. La ecuación de
Hermite admite como soluciones linealmente independientes en series de potencias
las siguientes:

𝜆(𝜆 − 2) ⋯ (𝜆 − 2𝑛 + 2) 2𝑛
𝑦1 (𝑥) = 1 + Σ∞ 𝑛 𝑛
𝑛=1 (−1) 2 𝑥 ,
(2𝑛)!

(𝜆 − 1)(𝜆 − 3) ⋯ (𝜆 − 2𝑛 + 1) 2𝑛+1
𝑦2 (𝑥) = 𝑥 + Σ∞ 𝑛 𝑛
𝑛=1 (−1) 2 𝑥 .
(2𝑛 + 1)!
Para 𝜆 = entero par, la función 𝑦1 (𝑥) es polinómica, y cuando 𝜆 = entero impar, es la
función 𝑦2 (𝑥) la que es un polinomio. Los polinomios se suelen escoger de forma que
el coeficiente dominante para el de grado 𝑛 sea 2𝑛 . Estos son los llamados polinomios
de Hermite 𝐻𝑛 (𝑥). Pueden obtenerse mediante la fórmula general:

2 𝑑 𝑛 −𝑥2
𝐻𝑛 (𝑥) = (−1)𝑛 𝑒𝑥 𝑒 .
𝑑𝑥𝑛

Los primeros son:

𝐻0 (𝑥) = 1 ,

𝐻1 (𝑥) = 2𝑥 ,

𝐻2 (𝑥) = 4𝑥2 − 2 ,

𝐻3 (𝑥) = 8𝑥3 − 12𝑥 ,

𝐻4 (𝑥) = 16𝑥4 − 48𝑥2 + 12 ,

Los polinomios de Hermite tienen simetría par o impar dependiendo del orden 𝑛:

𝐻𝑛 (−𝑥) = (−1)𝑛 𝐻𝑛 (𝑥) ,

verifican la ley de recurrencia:

𝐻𝑛+1 (𝑥) = 2𝑥𝐻𝑛 (𝑥) − 2𝑛𝐻𝑛−1 (𝑥).

Ecuaciones diferenciales
21
Tema 5. Ideas clave
 2
Los polinomios de Hermite son ortogonales respecto a la función peso 𝑤(𝑥) = 𝑒−𝑥
en (−∞, ∞), es decir:

∞
2
𝐻𝑚 (𝑥)𝐻𝑛 (𝑥) 𝑒−𝑥 𝑑𝑥 = √𝜋 2𝑛 𝑛! 𝛿𝑛𝑚 ,
∫−∞

y además forman un sistema completo.

Ecuación de Laguerre
La ecuación de Laguerre es:

𝑥𝑦″ + (1 − 𝑥)𝑦′ + 𝑛𝑦 = 0 , (13)

tiene soluciones no‐singulares solo si 𝑛 es un entero no‐negativo. Esta ecuación apa‐

rece en mecánica cuántica en la parte radial de la ecuación de Schrödinger para un
átomo hidrogenoide, entre otros problemas.

Las soluciones de la Ecuación (13) son los polinomios de Laguerre 𝐿𝑛 (𝑥), cuya fórmula
general es:
𝑛
𝑛 (−1)𝑘 𝑘
𝐿𝑛 (𝑥) = 𝑥 .
∑ (𝑘) 𝑘!
𝑘=0

Los polinomios de Hermite también se pueden definir mediante la fórmula de Rodri‐

gues:
𝑒𝑥 𝑑 𝑛 −𝑥 𝑛
𝐿𝑛 (𝑥) = 𝑒 𝑥 ),
𝑛! 𝑑𝑥𝑛 (
o mediante la fórmula de recurrencia:

(2𝑘 + 1 − 𝑥)𝐿𝑘 (𝑥) − 𝑘𝐿𝑘−1 (𝑥)

𝐿𝑘+1 (𝑥) = ,
𝑘+1

a partir de los dos primeros polinomios:

𝐿0 (𝑥) = 1 ,

𝐿1 (𝑥) = 1 − 𝑥 .

Además, los polinomios de Laguerre son ortogonales en [0, ∞) respecto a la función

Ecuaciones diferenciales
22
Tema 5. Ideas clave
peso 𝑤(𝑥) = 𝑒−𝑥 , es decir:

∞
𝐿𝑚 (𝑥)𝐿𝑛 (𝑥)𝑒−𝑥 𝑑𝑥 = 𝛿𝑚𝑛 .
∫0

Ecuación generalizada de Laguerre

La ecuación de Laguerre es un caso particular de la ecuación generalizada de Laguerre:

𝑥𝑦″ + (𝛼 + 1 − 𝑥)𝑦′ + 𝑛𝑦 = 0 ,

(𝛼)
cuyas soluciones son los polinomios generalizados (o asociados) de Laguerre 𝐿𝑛 (𝑥).
Los polinomios de Laguerre son un caso especial de los generalizados con 𝛼 = 0, es
(0)
decir, 𝐿𝑛 (𝑥) = 𝐿𝑛 (𝑥).

Los polinomios generalizados de Laguerre se pueden definir mediante la fórmula ge‐

neral:
(𝛼) 𝑛 𝑛 + 𝛼 (−1)𝑘 𝑘
𝐿𝑛 (𝑥) = Σ𝑘=0 𝑥 ,
(𝑛 − 𝑘) 𝑘!

mediante la fórmula de Rodrigues:

(𝛼) 𝑥−𝛼 𝑒𝑥 𝑑 𝑛 −𝑥 𝑛+𝛼

𝐿𝑛 (𝑥) = 𝑒 𝑥 ),
𝑛! 𝑑𝑥𝑛 (

o de manera recursiva, utilizando la fórmula de recurrencia:

(𝛼) (𝛼)
(𝛼) (2𝑘 + 1 + 𝛼 − 𝑥)𝐿𝑘 (𝑥) − (𝑘 + 𝛼)𝐿𝑘−1 (𝑥)
𝐿𝑘+1 (𝑥) = ,
𝑘+1

a partir de los dos primeros polinomios:

(𝛼)
𝐿0 (𝑥) = 1 ,
(𝛼)
𝐿1 (𝑥) = 1 + 𝛼 − 𝑥 .

Los polinomios generalizados de Laguerre son ortogonales en [0, ∞) respecto a la fun‐

ción peso 𝑤(𝑥) = 𝑥𝛼 𝑒−𝑥 :

∞
(𝛼) (𝛼) Γ(𝑛 + 𝛼 + 1)
𝐿𝑚 (𝑥)𝐿𝑛 (𝑥)𝑥𝛼 𝑒−𝑥 𝑑𝑥 = 𝛿𝑚𝑛 .
∫0 𝑛!

Ecuaciones diferenciales
23
Tema 5. Ideas clave
5.7 Cuaderno de ejercicios

Ejercicio 1. Estudiar si las siguientes ecuaciones tienen puntos singulares o no, y

clasificarlos:

1. 𝑥2 𝑦″ − 2𝑦 = 0.

2. 𝑥𝑦″ + 2𝑦′ + ln 𝑥 𝑦 = 0.

3. 𝑥 sin 𝑥𝑦″ + 3𝑦′ + 𝑥𝑦 = 0.

Ejercicio 2. Resolver por medio de series en torno a 𝑥 = 0:

1. 𝑦″ + 𝑥2 𝑦 = 0.

2. 𝑥𝑦″ + (1 − 𝑥)𝑦′ + 𝑦 = 0.

3. 𝑦″ − 𝑒𝑥 𝑦 = 0.

Ejercicio 3. Deducir a partir de la fórmula de Rodrigues para los polinomios de

Legendre las siguientes fórmulas de recurrencia:

▸ 𝑛𝑃𝑛 = 𝑥𝑃𝑛′ − 𝑃𝑛−1

′
.

▸ (𝑛 + 1)𝑃𝑛+1 − (2𝑛 + 1)𝑥𝑃𝑛 + 𝑛𝑃𝑛−1 = 0.

Ejercicio 4. Comprobar que la ecuación de Bessel tiene un punto singular no‐

regular en el infinito.

Ecuaciones diferenciales
24
Tema 5. Ideas clave
 Ejercicio 5. Comprobar:

2𝑝
1. 𝐽𝑝−1 + 𝐽𝑝+1 = 𝑥 𝑝
𝐽 .

2. 𝐽𝑝′ = 12 (𝐽𝑝−1 − 𝐽𝑝+1 ).

Ejercicio 6. Hallar 𝐽5/2 y 𝐽−5/2 .

5.8 Referencias bibliográficas

Berry, M. (2001). Why are special functions special? Physics Today, 54(4), 11–12.

González, J. A. A. & Manotas, A. T. (2016). Una breve discusión sobre algunas funciones
especiales. Sigma, 12(2), 1–30.

M.R. Spigel, L. A. (2005). Fórmulas y tablas de matemática aplicada.

Ecuaciones diferenciales
25
Tema 5. Ideas clave