Tema 5

Tema 5: Soluciones por medio de series y
funciones especiales
Ecuaciones diferenciales
Dr. Xavier Rivas

Grado en Fı́sica
Escuela Superior de Ingenierı́a y Tecnologı́a
Dr. Xavier Rivas PER3612 1 / 61

Índice
1 Introducción y objetivos
2 Puntos regulares
3 Puntos singulares y método de Frobenius
4 Ecuación de Legendre
5 Ecuación de Bessel
6 Otras ecuaciones

1
Introducción y objetivos

Introducción y objetivos
Identificar los puntos regulares, singulares y singulares-regulares

en una ecuación lineal de segundo orden.
Aprender a resolver ecuaciones lineales de segundo orden por
medio de series, cuando eso sea posible.
Reconocer las ecuaciones lineales más importantes en fı́sica y
conocer algunas de las funciones especiales que definen.

2
Puntos regulares

Funciones analı́ticas
Consideremos la ecuación
y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0 .
En principio, una ecuación de este tipo no admite soluciones en serie
para todo x ∈ R, sino que lo hace solamente en aquellos intervalos
en que las funciones a(x), b(x) son analı́ticas.
Recordamos que una función f (x) es analı́tica en un punto x0 si se
puede expresar como una serie de potencias (serie de Taylor) entorno
a dicho punto, es decir si
∞
f (n) (x0 )
(x − x0 )n ,
X
f (x) =
n=0 n!
en un cierto intervalo I 3 x0 . Para ello, la función tiene que ser
infinitamente derivable en x0 , y su serie de Taylor en x0 tiene que
tener un radio de convergencia R > 0.
Puntos regulares y singulares
Consideremos la EDO
y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0 .
Definición 2.1
Un punto x0 es un punto regular de la EDO si a(x) y b(x) son
analı́ticas en x0 . En caso contrario, x0 es un punto singular de la
EDO.
Si x0 es un punto regular, las funciones a(x), b(x) pueden escribirse
en serie de potencias,
∞ ∞
ak (x − x0 )k , bk (x − x0 )k ,
X X
a(x) = b(x) =
k=0 k=0
de modo que es razonable esperar que cualquier solución de la EDO

también sea una serie de potencias.
Teorema
Si x0 es un punto regular de y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0, su

solución general es
∞
ck (x − x0 )k = c0 y1 (x) + c1 y2 (x) ,
X
y(x) =
k=0
donde c0 , c1 son constantes arbitrarias y y1 (x), y2 (x) son soluciones

de la ecuación que admiten desarrollos en serie convergentes al
menos en el intervalo (x0 − R, x0 + R), donde R es el menor de los
dos radios de convergencia de a(x), b(x). Además, los coeficientes ck
para k ≥ 2 se pueden determinar de forma única en función de c0 y
c1 llevando la serie solución a la ecuación.

Puntos regulares y singulares
De esta manera obtenemos dos soluciones de la forma
y1 (x) = 1 + Σ , y2 (x) = (x − x0 ) + Σ ,
donde los sumatorios Σ contienen las potencias de x mayores o

iguales a 2. Estas soluciones satisfacen las condiciones iniciales
y1 (x0 ) = 1, y10 (x0 ) = 0, y2 (x0 ) = 0, y20 (x0 ) = 1, por lo que su
Wronskiano en x0 es no nulo y las dos soluciones son linealmente
independientes. Sabemos que las soluciones particulares existen y son
únicas, por ser a y b analı́ticas, y los valores y(x0 ) = y0 , y 0 (x0 ) = y00
determinan las constantes c0 = y0 y c1 = y00 .

3
Puntos singulares y método
de Frobenius

Puntos singulares regulares
Si x0 es un punto singular de la EDO
y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0 ,
es decir, si a(x), b(x) no son analı́ticas en x0 , entonces no podemos

aplicar el teorema anterior para encontrar una solución en dicho
punto. Sin embargo, en muchas ocasiones es precisamente cerca de
los puntos singulares donde nos interesa conocer las soluciones de la
EDO.
Afortunadamente existe un método para obtener las soluciones de la

EDO en un tipo particular de puntos singulares: los puntos singulares
regulares.

Definición 3.1
Un punto x0 es un punto singular regular si (x − x0 )a(x) y
(x − x0 )2 b(x) son analı́ticas en x0 .
Supongamos que x0 es un punto singular regular de la EDO
y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0 ,
y por simplicidad suponemos que x0 = 0. Para puntos singulares

regulares distintos de cero solo hay que considerar el cambio de
variable z = x − x0 . Multiplicamos la EDO por x2 :
x2 y 00 (x) + x2 a(x)y 0 (x) + x2 b(x)y(x) = 0 ,
y definimos
a∗ (x) = xa(x) , b∗ (x) = x2 b(x) .
Entonces, obtenemos
x2 y 00 (x) + xa∗ (x)y 0 (x) + b∗ (x)y(x) = 0 ,
donde a∗ (x), b∗ (x) son analı́ticas en x = 0:

∞ ∞
a∗ (x) = a∗k xk , b∗ (x) = b∗k xk .
X X
k=0 k=0
Estas expresiones son válidas en (−R, R), siendo R el menor de los

radios de convergencia de a∗ (x) y b∗ (x). Para establecer las
soluciones de la EDO, usaremos el método de Frobenius:

Método de Frobenius
La EDO x2 y 00 (x) + xa∗ (x)y 0 (x) + b∗ (x)y(x) = 0 tiene siempre una

solución y1 (x) de la forma
∞
y1 (x) = xr1 ck x k ,
X
c0 6= 0 ,
k=0
donde r1 es la mayor de las dos raı́ces reales r1 , r2 del polinomio

inicial definido como
q(r) = r2 + (a∗0 − 1)r + b∗0 .

La segunda solución linealmente independiente, y2 (x), es
Si r1 − r2 no es cero ni un entero positivo, entonces
∞
y2 (x) = xr2 bk x k ,
X
b0 6= 0 .
k=0
Si r1 = r2 , entonces
∞
y2 (x) = xr1 +1 bk xk + y1 (x) ln x .
X
k=0
Si n = r1 − r2 es un entero positivo, entonces

∞
r2
bk xk + ay1 (x) ln x ,
X
y2 (x) = x b0 6= 0 , a ∈ R .
k=0

Todas las soluciones están definidas al menos para x ∈ (0, R) y los

coeficientes ck , bk y la constante a pueden determinarse sustituyendo
cada una de las soluciones en la ecuación.
A partir de estas soluciones se obtienen otras que son válidas para

x ∈ (−R, 0) simplemente sustituyendo ln x por ln |x| y las
expresiones del tipo xr por |x|r .
A continuación veremos dos de las ecuaciones más relevantes que son

resolubles por medio de series. Estas ecuaciones serivarán como
ejemplos de los métodos descritos.

4
Ecuación de Legendre

La ecuación de Legendre es
(1 − x2 )y 00 − 2xy 0 + p(p + 1)y = 0 , donde p ∈ Z .
Esta ecuación aparece de forma natural al resolver la ecuación de

Laplace, que veremos en los temas dedicados a las EDPs, mediante
separación de variables en coordenadas esféricas. También aparece al
resolver la ecuación de Schrödinger para un potencial central en tres
dimensiones, por ejemplo, en el átomo de hidrógeno.
La ecuación de Legendre también puede escribirse como
2x 0 p(p + 1)
y 00 − y + y = 0.
1 − x2 1 − x2

De la ecuación anterior vemos que x = ±1 son puntos singulares

regulares. La ecuación anterior tiene soluciones analı́ticas al menos en
el intervalo |x| < 1. Vamos a buscar una solución en torno al punto
x = 0 de la forma ∞
ck x k .
X
y(x) =
k=0
Derivando la solución,
∞ ∞
y 0 (x) = ck kxk−1 = ck+1 (k + 1)xk ,
X X
k=1 k=0
∞ ∞
y 00 (x) = ck k(k − 1)xk−2 = ck+2 (k + 2)(k + 1)xk .
X X
k=2 k=0

Introduciendo estas series en la ecuación, llegamos a
∞
ck+2 (k + 2)(k + 1)xk − ck+2 (k + 2)(k + 1)xk+2
X
k=0

k+1 k
− 2ck+1 (k + 1)x + p(p + 1)ck x = 0.
A partir de esta ecuación, queremos encontrar una recurrencia para

calcular los coeficientes ck . Para ello, está claro que los coeficientes
de cada potencia xk deben ser nulos. Si nos fijamos en el término en
xk , vemos que
(ck+2 (k + 2)(k + 1) − ck k(k − 1) − 2ck k + p(p + 1)ck ) xk = 0 ,
de donde obtenemos la recurrencia
k(k + 1) − p(p + 1)
ck+2 = ck .
(k + 2)(k + 1)
La recurrencia anterior puede reescribirse como
(p − k + 2)(p + k − 1)
ck = − ck−2 ,
k(k − 1)
de donde podemos obtener todos los coeficientes ck a partir de dos
constantes arbitrarias c0 , c1 :
p(p + 1)
c2 = − c0 ,
2
(p − 1)(p + 2)
c3 = − c1 ,
6
(p − 2)(p + 3) p(p − 2)(p + 1)(p + 3)
c4 = − c2 = c0 ,
12 24
(p − 3)(p + 4) (p − 1)(p − 3)(p + 2)(p + 4)
c5 = − c3 = c1 ,
20 120
..
.
de donde podemos obtener la fórmula general para los coeficientes

pares e impares de la serie:
p(p − 2) · · · (p − 2n + 2)(p + 1)(p + 3) · · · (p + 2n − 1)
c2n = (−1)n c0 ,
(2n)!
(p − 1)(p − 3) · · · (p − 2n + 1)(p + 2)(p + 4) · · · (p + 2n)
c2n+1 = (−1)n c1 .
(2n + 1)!
Ası́ pues, la solución general de la ecuación de Legendre es

∞ ∞ ∞
k 2n
c2n+1 x2n+1 ,
X X X
y(x) = ck x = c0 + c1 x + c2n x +
k=0 n=1 n=1

Es decir,
y(x) = c0 y1 (x) + c1 y2 (x) ,
donde
y1 (x) =
∞
X p(p − 2) · · · (p − 2n + 2)(p + 1)(p + 3) · · · (p + 2n − 1) 2n
1+ (−1)n x ,
n=1
(2n)!
y2 (x) =
∞
X (p − 1)(p − 3) · · · (p − 2n + 1)(p + 2)(p + 4) · · · (p + 2n) 2n+1
x+ (−1)n x ,
n=1
(2n + 1)!
que convergen en el disco |x| < 1.

Las series y1 (x), y2 (x) se truncan para valores de p pares positivos e
impares positivos respectivamente. Si p = 2m, y1 (x) se convierte en
un polinomio de grado 2m:
p = 0: y1 (x) = 1,
p = 2: y1 (x) = 1 − 3x2 ,
p = 4: y1 (x) = 1 − 10x2 + 35 3
x4 ,
..
.
mientras que si p = 2m + 1, y2 (x) se convierte en un polinomio de
grado 2n + 1:
p = 1: y2 (x) = x,
p = 3: y2 (x) = x − 53 x3 ,
p = 5: y2 (x) = x − 143
x3 + 215
x5 ,
..
.
Polinomios de Legendre
Estos polinomios Pn (x), normalizados de modo que Pn (1) = 1, son
los llamados polinomios de Legendre:
P0 (x) = 1 ,
P1 (x) = x ,
3 1
P2 (x) = x2 − ,
2 2
5 3 3
P3 (x) = x − ,
2 2
35 4 15 2 3
P4 (x) = x − x + ,
8 4 8
63 5 35 3 15
P5 (x) = x − x + x,
8 4 8
..
.

Los polinomios de Legendre tienen muchas propiedades interesantes y
numerosas aplicaciones, como la expansión en serie de un potencial
1
de la forma y expansiones multipolares.
r
El polinomio Pn tiene n ceros reales, todos en el intervalo (−1, 1).
Como Pn (−x) = (−x)n Pn (x),
los polinomios Pn con n par tienen simetrı́a par,
los polinomios Pn con n impar tienen simetrı́a impar.

Los polinomios de Legendre se pueden definir de muchas maneras,
por ejemplo mediante la fórmula de Rodrigues:
1 dn 2
Pn (x) = n (x − 1) ,
2 n! dn x
o mediante la recurrencia
(n + 1)Pn+1 (x) = (2n + 1)xPn (x) − nPn−1 (x) .
Los polinomios de Legendre forman un sistema completo de
polinomios ortogonales, es decir que

Z 1 0
 n 6= m ,
Pm (x)Pn (x)dx = 2
−1 
 n = m,
2n + 1
que puede escribirse como
Z 1
2
Pm (x)Pn (x)dx = δnm .
−1 2n + 1
El hecho de que formen un sistema completo significa que cualquier

función f (x) continua a trozos en el intervalo [−1, 1] se puede
escribir como una serie de polinomios de Legendre de la forma
∞
X
f (x) = ak Pk (x) .
k=0
Más adelante veremos más sobre estas expansiones.

5
Ecuación de Bessel

La ecuación de Bessel es
x2 y 00 + xy 0 + (x2 − p2 )y = 0 , para p ≥ 0 .
Esta ecuación también aparece al resolver la ecuación de Laplace en

coordenadas cilı́ndricas y la ecuación de Hemholtz en coordenadas
esféricas. La ecuación de Hemholtz es una forma independiente del
tiempo de la ecuación de ondas que veremos más adelante.
Si reescribimos la ecuación de Bessel como

1 x2 − p 2
y 00 + y 0 + y = 0,
x x2
vemos claramente que x = 0 es un punto singular regular. Además,
vemos que a∗ (x) = 1 y b∗ (x) = x2 − p2 , y el primer término del
desarrollo es a∗0 = 1 y b∗0 = −p2 .
Podemos aplicar el método de Frobenius y proponer una solución de
la forma ∞
y1 (x) = xr1 ck x k ,
X
k=0
donde r1 es la mayor de las raı́ces del polinomio inicial, que en este

caso es
q(r) = r(r − 1) + r − p2 = r2 − p2 = (r + p)(r − p) ,
que tiene raı́ces r = ±p. Por tanto, debemos probar la solución

∞
ck xp+k ,
X
y1 (x) =
k=0
que converge en todo R.

Las derivadas de y1 (x) son
∞
y10 (x) = ck (p + k)xp+k−1 ,
X
k=0
∞
y100 (x) = ck (p + k)(p + k − 1)xp+k−2 .
X
k=0
Sustituyendo en la ecuación,
∞
X
ck (p + k)(p + k − 1)xp+k + ck (p + k)xp+k + ck xp+k+2 − p2 ck xp+k = 0

k=0
∞
X
ck (p2 + pk − p + pk + k 2 − k + p + k − p2 )xp+k + ck xp+k+2 = 0

k=0
∞
X
ck (2p + k)kxp+k + ck xp+k+2 = 0

(A)
k=0

Como los términos se tienen que anular, tenemos que
ck (2p + k)k + ck−2 ,
de modo que obtenemos la recurrencia

1
ck = − ck−2 .
(2p + k)k
Si escribimos los primeros términos de la ecuación (A),
c0 xp+2 +(2p+1)c1 xp+1 +c1 xp+3 +2(2p+2)c2 xp+2 +c2 xp+4 +· · · = 0 ,
vemos que la única manera de que el término en xp+1 se anule es

haciendo c1 = 0.

Por tanto, todos los coeficientes ck con k impar son cero, mientras
que los coeficientes con k par son
1
c2 = − c0 ,
2 · 2(p + 1)
1 1
c4 = − 2 c2 = 4 c0 ,
2 · 2(p + 2) 2 · 2(p + 2)(p + 1)
..
.
es decir,
c2n+1 = 0 ,
(−1)n
c2n = c0 .
22n n!(p + n)(p + n − 1) · · · (p + 1)

Ası́ que la primera solución de la ecuación de Bessel es

∞
(−1)n
!
p
x2n
X
y1 (x) = c0 x 1+ 2n n!(p + n)(p + n − 1) · · · (p + 1)
.
n=1 2
Esta solución se suele expresar en términos de la función Gamma:

La función Gamma
Definición 5.1
La función Gamma es la función
Z ∞
Γ(s) = xs−1 e−x dx , x > 0.
0
La función Gamma satisface la propiedad
Γ(s + 1) = sΓs ,
de modo que
Γ(s + n + 1) = (s + n)(s + n − 1) · · · (s + 1)sΓ(s) ,
y si n ∈ N, entonces Γ(s + 1) = s!.

Funciones de Bessel
Utilizando la función Gamma, podemos escribir la solución y1 como

∞
(−1)n Γ(p + 1) x
2n
p
X
y1 (x) = c0 x .
n=0 n!Γ(p + n + 1) 2
1
Esta solución, con c0 = , define la función de Bessel de
2p Γ(p + 1)
primera especie y orden p:
Definición 5.2
La función de Bessel de primera especie y orden p es
∞
(−1)n
p X 2n
x x
Jp (x) = .
2 n=0 n!Γ(p + n + 1) 2

Funciones de Bessel
Como ejemplo, veamos la forma de las funciones de Bessel de

primera especie con p = 0 y p = 1:
∞
(−1)n x 2n
X
J0 (x) = 2
,
n=0 (n!) 2
∞
(−1)n
2n+1
X x
J1 (x) = .
n=0 n!(n + 1)! 2
Del mismo modo que J0 y J1 , todas las Jp son funciones oscilatorias.

De hecho, para x suficientemente grande,
s
2 π

Jp (x) ∼ cos x − (2p + 1) .
πx 4

Funciones de Bessel
Cada función Jp tiene un número infinito de ceros en (0, ∞) que se

pueden encontrar tabulados en cualquier libro de funciones especiales.
Las funciones de Bessel son especialmente importantes en la
propagación de ondas y en los problemas con potenciales estáticos.
Al resolver problemas en sistemas de coordenadas cilı́ndricas se

obtienen funciones de Bessel de orden entero (p = n), en coordenadas
esféricas se obtienen órdenes semi-enteros (p = n + 1/2).

Funciones de Bessel
Las funciones de Bessel también se pueden definir de forma recursiva

mediante la fórmula
2p
Jp−1 (x) + Jp+1 (x) = Jp (x) .
x
Para valores no-enteros de p las funciones Jp (x) y J−p (x) son
linealmente independientes. Por otro lado, si p es entero (p = n) se
cumple que
J−n (x) = (−1)n Jn (x) .

Funciones de Bessel
Además, las funciones de Bessel con p semi-entero se pueden expresar
de manera exacta mediante senos y cosenos:
r r
2 2
J 1 (x) = sin x , J− 1 (x) = cos x ,
r
2 πx 2
rπx
2 1 2 1

J 3 (x) = sin x − cos x , J− 3 (x) = − sin x + cos x ,
2 πx x 2 πx x
.. ..
. .
Ahora vamos a buscar una segunda solución de la ecuación de Bessel

linealmente independiente de la primera. Por un teorema visto
anteriormente, tenemos tres soluciones dependiendo de si la
diferencia entre las raı́ces del polinomio inicial r1 − r2 es entera, no
entera, o cero. En este caso, r1 − r2 = 2p. Estudiemos las diferentes
posibilidades:
Funciones de Bessel
2p no es cero ni entero. La segunda solución de la ecuación de
Bessel es de la forma
∞
−p
bk x k ,
X
y2 (x) = x x > 0,
k=0
que coincide con la primera expresión de la solución y1 que

obtuvimos, salvo por el signo de p. Por tanto, para encontrar los
coeficientes bk no hay más que repetir los cálculos anteriores
sustituyendo p por −p. Ası́, llegamos a la solución
y2 (x) = J−p (x), con
∞
(−1)n
−p X 2n
x x
J−p (x) = .
2 n=0 n!Γ(−p + n + 1) 2

Funciones de Bessel
p es cero. Entonces, proponemos una solución del tipo

∞
bk xk+1 + J0 (x) ln x ,
X
y2 (x) = x > 0,
k=0
que define la función de Bessel de segunda especie y orden

0:
Definición 5.3
La función de Bessel de segunda especie y orden 0 es
∞
(−1)n+1
2n
1 1 x
X
Y0 (x) = 2
1 + + ··· + +J0 (x) ln x , x > 0.
n=1 (n!) 2 n 2

Funciones de Bessel
2p es entero. Entonces, la segunda solución a la ecuación de

Bessel es
∞
y2 (x) = x−p bk xk + aJp (x) ln x ,
X
x > 0.
k=0
Si 2p es impar, es decir si p = 12 , 32 , 25 , . . . , se puede demostrar

que a = 0 y la solución es y2 (x) = J−p (x).
Si 2p es par, es decir si p = 1, 2, 3, . . . , la solución es una

función de Bessel de segunda especie y orden p:

Funciones de Bessel
Definición 5.4
La función de Bessel de segunda especie y orden p, p ∈ N, es
p−1
!
−p 2p
1 x 1 1 1 x (p − n − 1)! x 2n
X
Yp (x) = − 1 + + ··· + +
2 2 p! 2 p 2 n! 2
n=0
∞
p X 2n
1 x (−1)n+1 1 1 1 1 x

− 1+ + ··· + + 1+ + ··· +
2 2 n!(n − p)! 2 n 2 n+p 2
n=0
+ Jp (x) ln x , x > 0.

Funciones de Bessel
Figura: Funciones de Bessel de primera y segunda especies
Las funciones de Bessel de segunda especie están relacionadas con las

de primera especie:
Jp (x) cos pπ − J−p (x)
Yp (x) = .
sin pπ
Si p = n es entero, Y−n (x) = (−1)n Yn (x).
6
Otras ecuaciones

Otras ecuaciones
A continuación veremos brevemente otras ecuaciones que tienen

relevancia en fı́sica y que definen otras series de funciones especiales:
ecuación general de Legendre,
ecuación de Hermite,
ecuación de Laguerre,
ecuación general de Laguerre.

Ecuación general de Legendre
La ecuación general de Legendre es

!
2 00 0 m2
(1 − x )y − 2xy + `(` + 1) − y = 0,
1 − x2
cuyas soluciones son las funciones asociadas de Legendre P`m (x).

Los ı́ndices `, m ∈ Z se llaman grado y orden de la función asociada
de Legendre.
La ecuación general de Legendre tiene soluciones no-nulas y

no-singulares en [−1, 1] solo si ` y m son enteros no negativos con
0 ≤ m ≤ ` o alguno de ellos es negativo con 0 ≤ |m| ≤ |`|.

Cuando m es par, la función P`m (x) es un polinomio. Cuando m = 0

y ` ∈ Z, las funciones P`m (x) coinciden con los polinomios de
Legendre.
Las funciones asociadas de Legendre forman parte de la definición de

otras funciones especiales, los esféricos armónicos, que tienen un
papel fundamental en el estudio del momento angular en mecánica
cuántica.

Caso `, m enteros no negativos. Las funciones P`m (x) están

relacionadas con los polinomios de Legendre:
m dm
P`m (x) = (−1)m (1 − x2 ) 2 P` (x) .
dxm
Las funciones asociadas de Legendre no son mutuamente
ortogonales en general. Por ejemplo, P11 no es ortogonal a P22 .
Sin embargo, algunos subconjuntos sı́ son ortogonales. Si
0 ≤ m ≤ `, para m fijo, las funciones satisfacen
Z 1
2(` + m)!
Pkm (x)P`m (x)dx = δk` .
−1 (2` + 1)(` − m)!

Caso m negativo y/o ` negativo. La ecuación general de
Legendre es invariante bajo un cambio de signo de m. Las
funciones con m negativo son proporcionales a las de m positivo:
(` − m)! m
P`−m = (−1)m P
(` + m)! `
A partir de la definición, se puede ver que las funciones

asociadas de Legendre son pares o impares de acuerdo a
P`m (−x) = (−1)`+m P`m (x) ,
y cumplen la recurrencia
m
(` − m + 1)P`+1 (x) = (2` + 1)xP`m (x) − (` + m)P`−1
m
(x) .

Ecuación de Hermite
La ecuación de Hermite de orden λ es
y 00 − 2xy 0 + 2λy = 0 ,
y aparece en el problema del oscilador armónico cuántico entre otros.

La ecuación de Hermite admite como soluciones linealmente
independientes en series de potencias las siguientes:
∞
λ(λ − 2) · · · (λ − 2n + 2) 2n
(−1)n 2n
X
y1 (x) = 1 + x ,
n=1 (2n)!
∞
(λ − 1)(λ − 3) · · · (λ − 2n + 1) 2n+1
(−1)n 2n
X
y2 (x) = x + x .
n=1 (2n + 1)!

Para λ un entero par, la función y1 (x) es polinómica. Cuando λ es un

entero impar, la función y2 (x) es polinómica. Los polinomios se
suelen escoger de forma que el coeficiente dominante para el de
grado n sea 2n . Estos son lo llamados polinomios de Hermite
Hn (x). Pueden obtenerse mediante la fórmula general
2 dn −x2
Hn (x) = (−1)n ex e .
dxn
Los primeros polinomios de Hermite son
H0 (x) = 1 , H1 (x) = 2x , H2 (x) = 4x2 − 2 ,

H3 (x) = 8x3 − 12x , H4 (x) = 16x4 − 48x2 + 12 .

Los polinomios de Hermite tienen simetrı́a par o impar en función de

n:
Hn (−x) = (−1)n Hn (x) ,
y verifican la recurrencia
Hn+1 (x) = 2xHn (x) − 2nHn−1 (x) .
Los polinomios de Hermite son ortogonales respecto a la función peso

2
w(x) = e−x en (−∞, +∞), es decir
Z +∞
2 √ n
Hm (x)Hn (x)e−x dx = π2 n!δnm ,
−∞
y forman un sistema completo.

Ecuación de Laguerre
La ecuación de Laguerre es
xy 00 + (1 − x)y 0 + ny = 0 ,
y tiene soluciones no-singulares solo si n es un entero no-negativo.

Esta ecuación aparece en mecánica cuántica en la parte radial de la
ecuación de Schrödinger para un átomo de hidrógeno, entre otros
problemas.
Las soluciones de la ecuación de Laguerre son los polinomios de

Laguerre Ln (x), cuya fórmula general es
n
n (−1)k k
!
X
Ln (x) = x .
k=0 k k!

Ecuación de Laguerre
Los polinomios de Laguerre también se pueden definir mediante la
fórmula de Rodrigues:
ex dn −x n
Ln (x) = (e x ) ,
n! dxn
(2k + 1 − x)Lk (x) − kLk−1 (x)
Lk+1 (x) = ,
k+1
a partir de los polinomios
L0 (x) = 1 , L1 (x) = 1 − x .
Además, los polinomios de Laguerre son ortogonales [0, +∞)
respecto a la función peso w(x) = e−x , es decir
Z ∞
Lm (x)Ln (x)e−x dx = δmn .
0

Ecuación generalizada de Laguerre
La ecuación de Laguerre que hemos visto es un caso particular de la

ecuación generalizada de Laguerre
xy 00 + (α + 1 − x)y 0 + ny = 0 ,
cuyas soluciones son los polinomios generalizados (o asociados)

de Laguerre L(α) n (x). Los polinomios de Laguerre son un caso
especial de los generalizados con α = 0, es decir
L(0)
n (x) = Ln (x) .

Los polinomios generalizados de Laguerre se pueden definir mediante

la fórmula general
n
n + α (−1)k k
!
L(α)
X
n (x) = x ,
k=0 n−k k!
mediante la fórmula de Rodrigues
x−α ex dn −x n+α
L(α)
n (x) = (e x ) ,
n! dxn

(α) (α)
(α) (2k + 1 + α − x)Lk (x) − (k + α)Lk−1 (x)
Lk+1 (x) = ,
k+1
a partir de los dos primeros polinomios
(α) (α)
L0 (x) = 1 , L1 (x) = 1 + α − x .
Los polinomios generalizados de Laguerre son ortogonales en [0, +∞)

respecto a la función peso w(x) = xα e−x :
Z ∞
α −x Γ(n + α + 1)
L(α) (α)
m (x)Ln (x)x e dx = δmn .
0 n!


Tema 5

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

Tema 5

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tema 5

Cargado por

Copyright:

Formatos disponibles

Tema 5: Soluciones por medio de series y

Dr. Xavier Rivas

Dr. Xavier Rivas PER3612 1 / 61

3 Puntos singulares y método de Frobenius

Dr. Xavier Rivas PER3612 2 / 61

Dr. Xavier Rivas PER3612 3 / 61

Identificar los puntos regulares, singulares y singulares-regulares

Dr. Xavier Rivas PER3612 4 / 61

Dr. Xavier Rivas PER3612 5 / 61

de modo que es razonable esperar que cualquier solución de la EDO

Si x0 es un punto regular de y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0, su

donde c0 , c1 son constantes arbitrarias y y1 (x), y2 (x) son soluciones

Dr. Xavier Rivas PER3612 8 / 61

De esta manera obtenemos dos soluciones de la forma

donde los sumatorios Σ contienen las potencias de x mayores o

Dr. Xavier Rivas PER3612 9 / 61

Dr. Xavier Rivas PER3612 10 / 61

Si x0 es un punto singular de la EDO

y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0 ,

es decir, si a(x), b(x) no son analı́ticas en x0 , entonces no podemos

Afortunadamente existe un método para obtener las soluciones de la

Dr. Xavier Rivas PER3612 11 / 61

Supongamos que x0 es un punto singular regular de la EDO

y 00 (x) + a(x)y 0 (x) + b(x)y(x) = 0 ,

y por simplicidad suponemos que x0 = 0. Para puntos singulares

x2 y 00 (x) + x2 a(x)y 0 (x) + x2 b(x)y(x) = 0 ,

x2 y 00 (x) + xa∗ (x)y 0 (x) + b∗ (x)y(x) = 0 ,

donde a∗ (x), b∗ (x) son analı́ticas en x = 0:

Estas expresiones son válidas en (−R, R), siendo R el menor de los

Dr. Xavier Rivas PER3612 13 / 61

La EDO x2 y 00 (x) + xa∗ (x)y 0 (x) + b∗ (x)y(x) = 0 tiene siempre una

donde r1 es la mayor de las dos raı́ces reales r1 , r2 del polinomio

q(r) = r2 + (a∗0 − 1)r + b∗0 .

Dr. Xavier Rivas PER3612 14 / 61

Si n = r1 − r2 es un entero positivo, entonces

Dr. Xavier Rivas PER3612 15 / 61

Todas las soluciones están definidas al menos para x ∈ (0, R) y los

A partir de estas soluciones se obtienen otras que son válidas para

A continuación veremos dos de las ecuaciones más relevantes que son

Dr. Xavier Rivas PER3612 16 / 61

Dr. Xavier Rivas PER3612 17 / 61

(1 − x2 )y 00 − 2xy 0 + p(p + 1)y = 0 , donde p ∈ Z .

Esta ecuación aparece de forma natural al resolver la ecuación de

La ecuación de Legendre también puede escribirse como

Dr. Xavier Rivas PER3612 18 / 61

De la ecuación anterior vemos que x = ±1 son puntos singulares

Dr. Xavier Rivas PER3612 19 / 61

A partir de esta ecuación, queremos encontrar una recurrencia para

de donde podemos obtener la fórmula general para los coeficientes

Ası́ pues, la solución general de la ecuación de Legendre es

Dr. Xavier Rivas PER3612 22 / 61

que convergen en el disco |x| < 1.

Dr. Xavier Rivas PER3612 25 / 61

Dr. Xavier Rivas PER3612 26 / 61

El hecho de que formen un sistema completo significa que cualquier

Más adelante veremos más sobre estas expansiones.

Dr. Xavier Rivas PER3612 28 / 61

Dr. Xavier Rivas PER3612 29 / 61