3.

Soluciones por medio de series

En el capítulo anterior vimos las escasas formas de resolver elementalmente la ecuación con coeficientes
variables [e] x"+ a(t)x' + b(t)x = 0 . Este capítulo trata una forma general de atacarla: suponer la solución
desarrollada en serie de potencias e introducir esta serie en la ecuación para determinar sus coeficientes.
En la sección 3.1 recordaremos la definición de función analítica (función descrita por una serie de
potencias convergente) y algunas manipulaciones matemáticas que se pueden hacer con ellas.
Cuando a y b son analíticas en un entorno de to (punto regular, sección 3.2) siempre se podrán
encontrar dos soluciones linealmente independientes de [e] en forma de serie de potencias por el siguiente
camino: llevando la serie a la ecuación conseguiremos expresar sus coeficientes ck en función de los dos
primeros c0 y c1 , que serán las dos constantes arbitrarias que aparecen en la solución de cualquier
ecuación de segundo orden (algunas veces podremos dar la expresión general del ck , pero otras nos
deberemos limitar a ir calculando coeficiente a coeficiente). Un teorema, que aceptaremos sin demostración,
nos asegurará que las series solución son convergentes al menos en el intervalo en que a y b lo eran.
Imponer datos iniciales en to será inmediato, pues tendremos que x(to)= c0 y x'(to)= c1 .
Si a y b no son analíticas, pero 'poco' (si to es singular regular), también se pueden utilizar series
para resolver [e] de una forma sólo algo más complicada (es el método de Frobenius de la sección 3.3).
Calcularemos primero una solución x1 que (si to = 0 ) será siempre de la forma t r ∑ (siendo r una de las
raíces del llamado polinomio indicial) y a continuación otra x2 , linealmente independiente de la anterior, que
unas veces será del mismo tipo y otras contendrá además un término incluyendo el ln t . De nuevo un
teorema no demostrado garantizará la convergencia de las series que vayan apareciendo.
Los cálculos de los coeficientes de las series son sencillos (aunque algo pesados). El problema básico
de utilizar series para resolver ecuaciones es la dificultad de obtener información sobre las soluciones que se
obtienen (y más cuando no podamos hallar siquiera su término general). Sin embargo ecuaciones del tipo [e]
surgen a menudo en problemas reales y las series son el único instrumento para resolverlas. Por eso existen
libros enteros (los de funciones especiales de la física) dedicados a estudiar las propiedades de las series
solución de algunas de estas ecuaciones (las de Legendre, Hermite, Bessel, Laguerre, Tchebycheff, …) .
Una mínima muestra de tales estudios son las propiedades de las series solución de las ecuaciones de
Legendre, Hermite y Bessel que se tratan en la sección 3.4. Para estudiar para grandes valores de t
esas ecuaciones, introduciremos de paso el llamado punto del infinito de una ecuación.

43
3.1 Funciones analíticas

Recordemos que una función f(t) es analítica en un punto t= t o si viene descrita por una serie de
potencias en un entorno de dicho punto:
∞
f(t) = ∑ ck(t–to)k = c0+c1(t–to) + c2(t–to)2+…
k=0
∞
A partir de ahora supondremos que to=0 (si no, haríamos t–to = s y estaríamos en ese caso): f(t) = ∑ c k t k .
k=0
Como se sabe, a cada serie de potencias está asociado un radio de convergencia R tal que:
si R=0 , la serie sólo converge en t=0 ; si R=∞ , converge para todo t ;
si 0<R<∞ , converge si |t|<R y diverge si |t|>R (en t=±R puede converger o divergir).
Además, si 0<to<R , la serie converge uniformemente en [–to,to] .
El R se puede calcular en muchas ocasiones aplicando el llamado criterio del cociente en valor absoluto:
∞ |ak+1|
Sea ∑ ak y sea p = klím . Entonces si p<1 , la ∑ converge, y si p>1 diverge.
k=0 → ∞ | a k|
Propiedad básica de las series de potencias es que, si R>0 , se pueden derivar e integrar término a término:
∞ ∞
f '(t) = ∑ kcktk-1 = c1+2c2 t +… , f"(t) = ∑ k(k–1)cktk–2 = 2c2+6c3 t +… , … ( ⇒ f(k)(0)= k! ck )
k=1 k=2
∞ ∞
∫ ∑ cktk k k+1c c
= C + ∑ k+1 t = C + c0 t + 21 t2 +… si |t|<R
k=0 k=0

También pueden sumarse, multiplicarse,… las series de potencias como si fuesen polinomios:
∞ ∞
f(t) = ∑ aktk si |t|<Rf y g(t) = ∑ bktk si |t|<Rg ⇒
∞ k=0 k=0
f(t)+g(t) = ∑ [ak+bk]tk , f(t)g(t)=a0b0+(a0b1+a1b0)t+ (a0b2+a1b1+a2b0)t2 + … , si |t|<min(Rf,Rg) .
k=0
(también el cociente f/g es analítico si tiene límite en t=0 ; ya veremos como hallar su desarrollo)
∞ (k)
f (0) k
Caso importante de estas series son las de Taylor: ∑ k!
t , para f con infinitas derivadas en 0 .
k=0
La mayoría de las funciones elementales coinciden con su serie de Taylor (y por tanto son analíticas) en todo
el intervalo de convergencia de esta serie. Por ejemplo, se tiene:
∞ ∞ ∞ ∞ ∞
tk [–1] k t2k+1 [–1] k t 2 k t2k+1 t2k
et = ∑ k ! , sen t = ∑ (2k+1)!
, cos t = ∑ (2k)!
, sht = ∑ , cht = ∑ ∀t∈R
k=0 k=0 k=0 k=0 (2k+1)! k=0 (2k)!

1 ∞
k
∞
[–1] k t k+1 ∞
[–1] k t2k+1 α α[α–1] 2
= ∑
1 – t k=0
t , ln(1+t) = ∑ k+1
, arctant = ∑ 2k+1
, [1+t ] = 1 + αt +
2!
t +… si |t |<1
k=0 k=0

Aunque una f sea C∞∈(R) y su serie de Taylor converja ∀t puede que ambas no coincidan, como le ocurre
a f(t)=e–1/t2, f(0)=0 que cumple f(k)(0)=0 ∀k (su serie de Taylor es Σ0.t k = 0 ). f , por tanto, no es analítica.
Para que una f sea analítica, es necesario que tenga infinitas derivadas en el punto, pero no es suficiente.
1
Ej 1. Hallemos de varias formas (algunas nada naturales) el desarrollo de f(t) = .
[1+t]2
∞ ∞ ∞
d 1 1 d
f(t) = –
dt 1+t
→
[1+t]2
=–
dt ∑ (–t)k = – k=1
∑ (–1)kktk–1 = ∑ (–1)k(k+1)tk = 1–2t+3t2–4t3+… si |t|<1 .
k=0 k=0

–2[–2–1] 2 –2[–2–1][–2–2] 3
Otra forma: [1+ t ]–2 = 1 –2t + t + t +… = 1–2t+3t2–4t3+… si |t|<1 .
2! 3!

Multiplicando: f(t) = [1–t+t2–…][1–t+t2–…] = [1 + (–1–1) t + (1+1+1) t2 + … ] = 1–2t+3t2–… si |t|<1

También podemos 'dividir': buscar una serie Σcktk tal que [c0+c1t+c2t2+c3t3+…] [t2+2t+1] = 1 ;
así se va obteniendo: c0=1 ; 2c0+c1=0 → c1=–2c0=–2 ; c0+2c1+c2=0 → c2=–c0–2c1=3 ; …

[el radio R del desarrollo de un cociente P/Q , con P y Q polinomios, simplificados los factores
comunes, y Q(0)≠0 es la distancia al origen de la raíz (real o compleja) de Q más próxima].

44
3.2 Puntos regulares

Sea la ecuación [e] x" + a(t)x' + b ( t ) x = 0

Se dice que t= to es un punto regular de [e] si a y b son analíticas en t= to .

En caso contrario se dice que t= to es punto singular de [e].
Sea t=0 regular y llamemos R al menor de los radios de convergencia de a y b . Se podrá, pues, escribir :
∞ ∞
a(t) = ∑ a k tk , b(t) = k=0
∑ b k tk para |t |<R
k=0

Si sus coeficientes son analíticos, se puede esperar que cualquier solución de [e] también lo sea. De hecho
se tiene el siguiente teorema que no demostramos:
∞
Teor 1.
Si t=0 es regular, la solución general de [e] es x = ∑ cktk = c0 x1 + c1 x2 ,
k=0
donde c0 , c1 son arbitrarios y x1 , x2 son soluciones de [e] que admiten desarrollos
en serie convergentes, al menos, en el intervalo (–R,R) . Además, los coeficientes
ck (para k≥2 ) , se pueden determinar de forma única en función de c0 y c1 llevando
la serie de x a la ecuación [e] (con los coeficientes a(t) y b(t) desarrollados).
.

[Este último desarrollo, desde luego, no será necesario si a y b son polinomios].

Obtendremos así las soluciones en la forma x1 = 1+Σ , x2 = t+Σ , donde cada una de las series contienen
potencias de t mayores o iguales que 2 . Las soluciones halladas x1 y x2 satisfarán x1(0)=1, x1'(0)=0 ,
x2(0)=0, x2'(0)=1 , así que su wronskiano en t=0 es |W|(0)=1 y, por tanto, son linealmente independientes.
Si lo que queremos es la solución de [e] que satisface x(0)=xo , x'(0)=xo' (que por ser a y b analiticas
existe y es única) dada la forma de las series de x1 y x2 es inmediato que debe tomarse c0=xo , c1=xo' .
Para estudiar las soluciones de [e] en un entorno de otro to regular, el cambio de variable s= t–to lleva a
una ecuación en s para la que s=0 es regular. Aplicando lo anterior para hallar su solución obtenemos
∞ ∞
x= ∑ c k sk , es decir , x = k=0
∑ c k (t–to)k .
k=0

En la práctica operaremos como en los ejemplos siguientes:

2t 2
Ej 1. (1+t 2 ) x" +2tx ' – 2 x = 0 , es decir, x " + x' – x = 0 , a(t) y b(t) analíticas en para | t |<1 .
1+t2 1+t2
Existen, según el teorema, soluciones que son series de potencias convergentes al menos en (–1,1) .
Probamos una serie llevándola junto a sus derivadas a la ecuación inicial (mucho más corto que desarollar
a y b ) e intentamos escribir los ck en funcion de c0 y c1 . Para | t |<1 se tiene entonces:
∞ ∞ ∞ ∞ ∞ ∞
x = ∑ c k t k , x ' = ∑ kc k t k–1 , x" = ∑ k(k–1)c k t k–2 → ∑ [ k(k–1)c k t k–2+k(k–1)cktk] + ∑ 2kc k t k – ∑ 2c k tk = 0
k=0 k=1 k=2 k=2 k=1 k=0
Para que ocurra esto es necesario que el coeficiente de cada potencia de t se anule (hay que tener
cuidado con los primeros términos porque las series empiezan a aportar términos para distintos k ):
t0 : 2.1.c2 –2.c0 = 0 → c2 = c0
t1 : 3.2.c3 + [2–2]c1 = 0 → c3 = 0
………………………………………………
tk–2 : k(k–1)ck + [(k–2)(k–3)+2(k–2)–2]ck-2 = 0
De la última expresión deducimos la regla de recurrencia que nos proporciona el ck en función de los
anteriores ya conocidos (en este ejemplo, sólo en función de ck-2 ); factorizaremos, para facilitar los
cálculos, los polinomios que aparezcan calculando sus raíces:
k(k–3) k–3
ck = – c =– c , k=2,3,…
k(k–1) k-2 k–1 k-2
A partir de esta regla calculamos algunos ck más en la búsqueda de la expresión del término general:
1 1 3 1 k+1 1
c4 = – c2 = – c0 , c6= – c4 = c0 , … → c2k = (–1) c , k=1,2,…
3 3 5 5 2k–1 0
c5 =0 por estar en función de c3 que se anulaba. Análogamente c7 = c9 =…= 0

45
 Por tanto, la solución general es (los c0 y c1 quedan indeterminados como aseguraba el teorema):
∞
1 4 1 6 k+1 t2k
x = c0x1+c1x2 = c0 [1+ t2– t + t +…] + c1 t = c0 [1 + ∑ (–1) ] + c1 t
3 5 k=1 2k–1
El teorema predecía que las series solución convergerían al menos en (–1,1) . Esto es lo que sucede. La
serie de x1 (como se puede ver por el criterio del cociente) tiene R=1 . La 'serie' de x2 (truncada a partir
de su segundo término) converge para todo t .
[Esta ecuación se podía haber resuelto sin usar series. Era fácil darse cuenta de que
x2 = t era una solución, y podríamos haber calculado la x1 siguiendo la sección 2.2:
–∫ a 2t

∫ x2–2 e ∫ t–2 e ∫ ∫ t–2 (1 +t2)–1 dt = –1 – arctan t

–
x1 = x2 dt = t 1+t2 dt = t

cuyo desarrollo, salvo el signo, coincide con el obtenido anteriormente].

Hemos tenido mucha suerte en el ejemplo 1 encontrando una regla de recurrencia tan sencilla y
pudiendo expresar el término general. Esto no sucede siempre como nos muestra el segundo ejemplo:
∞
Ej 2. x " + (t–2)x = 0 , x(0)=2, x '(0)=1 a(t)=0 y b(t) =t–2 analíticas en todo R → x = ∑ c k tk .
k=0

Llevando x , x' , x" a la ecuación e igualando a 0 los coeficientes de las diferentes potencias de t :
∞ ∞
∑ k(k–1)c k tk–2 + ∑[ ck tk+1 – 2c k tk ] = 0 →
k=2 k=0

t0 : 2.1.c2 –2.c0 = 0 → c2 = c0
1 1
t1 : 3.2.c3 + c0 – 2.c1 = 0 → c3 = – c0 + c1
6 3
…………………………………………………………………
1 2
tk–2 : k(k–1)ck + ck-3 – 2ck-2 = 0 → ck = – c + c , k=3,4,…
k(k–1) k-3 k(k–1) k-2
(regla de recurrencia de tres términos que suele traer muchos más problemas que las de dos).
Escribimos un par de términos más en función de c0 y c1 :
1 2 1 1 1 2 1 1
c4 = – c + c = c – c ; c5= – c + c =– c + c
12 1 12 2 6 0 12 1 20 2 20 3 15 0 30 1
No hay forma de encontrar la expresión del término general, aunque paso a paso podemos ir calculando
el número de términos que queramos.
1 1 1 5 1 1 4 1 5
La solución general es entonces: x = c0 [ 1+ t2 – t3 + t4 – t +… ] + c1 [ t + t3 – t + t +… ]
6 6 15 3 12 30
1 1
Y la particular pedida: x = 2 + t + 2t2 + t4 – t5 +… , serie convergente para todo t según el teorema.
4 10
Para calcular unos pocos términos (pero no para calcular muchos o para buscar la expresión del
término general) de una serie solución cerca de un punto regular se puede seguir el siguiente camino:
Haciendo t=0 en la ecuación: x" (0) + (0–2) x(0) = 0 → x" (0) = 4
Derivando la ecuación y volviendo a hacer t=0 : x"' +(t–2) x ' +x = 0 → x"'(0) = 2x '(0)–x(0) = 0
Derivando otra vez : x"" +(t–2) x" +2x ' = 0 → x""(0) = 2x"(0)–2x '(0) = 6 , ……
x"(0) 2 x"'(0) 3 4 2 0 3 4 4
Por tanto: x(t) = x(0) + x '(0) t + t + t +…= 2+t+ t + t + t +…
2 6 2 6 24

Si lo que buscásemos fuese la solución con x(2)=6, x '(2)=0 , la solución general de arriba, dada por
series en t=0 , claramente no sirve para imponer datos en t=2 . Debemos resolver en torno a este punto:
∞ ∞ ∞
d2x
s =t–2 → + s x = 0, s=0 regular → x = ∑ c k sk → ∑ k(k–1)c k sk–2 + ∑ c k sk+1 = 0 →
ds2 k=0 k=2 k=0
0 1 1 k–2 1
s : 2c2 = 0 ; s : 6c3 +c0 = 0 , c3 = – c0 ; … ; s : k(k–1)ck + ck-3 = 0 → ck = – c , k=3,4,…
6 k(k–1) k-3
1 1 1 1 1
→ c5 =c8 =… =0 ; c4 = – c ; c6 = – c = c ; c7 = – c = c →
4.3 1 6.5 3 6.5.3.2 0 7.6 4 7.6.4.3 1
1 3 1 6 1 4 1 6 1 6 1
s +…] → x =6 – s3 +
datos
x = c0 [ 1 – s + s +…] +c1 [ s – s + s +… = 6 – (t–2)3 + (t–2)6 +…
6 120 12 504 20 20
[ahora sí podríamos dar la expresión de su término general].

46
3.3 Puntos singulares regulares

Supondremos en esta sección que t= to es un punto singular de [e] x"+ a(t)x'+ b(t)x = 0 , es decir, que o
a(t) o b(t) o las dos no son analíticas en t= to , con lo que no es aplicable el método de la sección anterior.
Sin embargo, interesa precisamente en muchas ocasiones conocer el comportamiento de las soluciones de
[e] en las cercanías de sus puntos singulares. En general se podrá decir poco sobre este comportamiento,
salvo para un tipo particular de puntos sólo débilmente singulares: los singulares regulares.
to es punto singular regular de [e] si (t–to)a(t) y (t–to)2b(t) son analíticas en to .

1 1
Ej 1. t(t–1) 2 x" – t x ' +(t–1) x = 0 , es decir, x" –
2 x' + x=0.
(t–1) t(t–1)
t =0 y t =1 son puntos singulares de la ecuación (los demás son regulares).
Como además –t/(t–1)2 y t/(t–1) son analíticas en t=0 , este punto es singular regular.
Como –1/(t–1) no es analítica en 1 (aunque (t–1)/t sí lo sea), el punto t =1 es singular no regular.

Supondremos a partir de ahora que el punto singular regular de [e] que tratamos es t=0 . Sabemos que no
supone pérdida de generalidad ya que en el caso de que queramos estudiar las soluciones cerca de un to≠0
el cambio s= t–to nos traslada el problema al estudio de las soluciones cerca de 0 de la ecuación en s .
Multiplicando [e] por t2 y llamando a*(t) = ta(t) y b*(t) = t2b(t) obtenemos la ecuación:

[e*] t 2 x" + t a*(t) x' + b*(t) x = 0

Se trata, por tanto, de resolver [e*] cerca de t=0 suponiendo que a*(t) y b*(t) son analíticas en ese punto,
es decir, que admiten desarrollo en serie válido en |t |<R (mínimo de los radios de convergencia):
∞ ∞
a*(t) = ∑ a *k tk , b*(t) = k=0
∑ b *k tk
k=0

Para hallar las soluciones de [e*] se utiliza el llamado método de Frobenius, que detallaremos en el
teorema de esta sección. Aunque no lo demostramos, intentamos con las siguientes consideraciones hacer
creíbles sus hipótesis y conclusiones. La ecuación de la forma [e*] más sencilla es la conocida ecuación de
Euler (para ella a*(t) y b*(t) son 'series' que se reducen a su primer término). Viendo sus soluciones está
claro que no existen, en general, soluciones analíticas de [e*]. Pero ya que hay soluciones de Euler de la
forma t r se podría pensar que [e*] posee soluciones en forma de serie que comiencen por términos t r .
∞
Probemos por tanto en [e*] la solución x = t r ∑ cktk = c0 t r + c1 t r+1 + c2 t r+2 + … Debe ser entonces:
k=0
∞ ∞ ∞ ∞ ∞
∑ a *ktk ) ( k=0
∑ (k+r)(k+r–1)ck tk+r + ( k=0 ∑ (k+r)cktk+r ) + ( k=0
∑ b *ktk ) ( k=0
∑ cktk+r ) = 0
k=0

El coeficiente que acompaña a la potencia de menor orden ( t r ) debe anularse: [ r(r–1) + a*0 r + b*0 ] c0 = 0 .
Si la serie ha de empezar por términos de exponente r , debe ser c0 ≠0 . Por tanto, los únicos valores de r
para los que pueden existir soluciones no triviales de la forma t r ∑ son las raíces del polinomio:

q(r) = r(r–1) + a*0 r + b*0 , llamado polinomio indicial de [e*]

Lo anterior es coherente con lo conocido sobre ecuaciones de Euler. Para ellas, si el polinomio tenía dos
raíces distintas r1 y r2 dos soluciones linealmente independientes de la ecuación eran t r 1 y t r 2 . Si la raiz
era doble, sin embargo, sólo existía una solución de esa forma, y la segunda era la primera multiplicada por el
ln t ; por tanto, también es de esperar que en la solución general de [e*] aparezcan logaritmos.
Pero al resolver por series la ecuación [e*] pueden aparecer problemas que no se presentan en el caso
particular de las ecuaciones de Euler. Igualando a 0 el coeficiente que acompaña a t r+k tenemos:
[ (r+k)(r+k–1)+(r+k)a*0 + b*0 ] ck + [ (r+k–1)a*1+b*1 ] ck–1 + … = 0
donde los puntos representan términos que incluyen ck–2 , ck–3 , … De esta expresión podemos despejar
el ck en función de los anteriores ya calculados siempre que el corchete que le acompaña, que es q(r+k) (el
polinomio indicial evaluado en r+k ), no se anule. Si r1 es la mayor de las dos raices q(r1+k)≠0 para todo k .
Pero si r2 es la menor, y la diferencia r1–r2 es un entero positivo n , el q(r2+k)=0 si k=n , y, salvo que los
demás sumandos también se anulen (con lo que cn quedaría indeterminado), no hay forma de anular el
coeficiente de t r2+n y no pueden existir soluciones t r2 Σ .

47
 Enunciamos ya el teorema (aunque se podría considerar el caso de raíces complejas del polinomio
indicial, nos limitamos, por sencillez, a los casos reales):

Teor 1 . Supongamos que el polinomio indicial q(r)= r(r–1)+ a*0 r+ b*0 tiene raíces reales r1 y r2 con r1≥r2 .
r1 ∞
Entonces siempre existe una solución x1 de [e*] de la forma x1 =t ∑ c k tk , c0≠0 .
k=0

La segunda solución x2 linealmente independiente es, según los casos:

r2 ∞
a] Si r1–r2 no es cero ni entero positivo: x2 = t ∑ b k tk , b0≠0
k=0
r1+1 ∞
b] Si r1 = r2 , x2 = t ∑ b k tk + x1 ln t
k=0
r2 ∞
c] Si r1–r2 = n entero positivo, x2 = t ∑ b k tk + a x1 ln t , b0≠0 , a∈R
k=0

Todas las soluciones están definidas al menos para 0<t<R y los coeficientes ck , bk y la
constante a se pueden determinar sustituyendo cada una de las soluciones en la ecuación.
.

Se comprueba sin dificultad que a partir de las soluciones anteriores obtenemos otras válidas para –R<t<0
sin más que sustituir ln t por ln|t| y las expresiones de la forma t r que preceden a las series por |t | r .
En el caso c] la constante a puede perfectamente ser 0 (como ocurre en las ecuaciones de Euler), con lo
que, a pesar de todo, hay dos soluciones independientes de la forma t r ∑ .

1 1
Ej 2. , o sea, t2 x"+ t x ' + t2 x = 0 .
2 t x " +x ' +t x = 0
2 2
1 1
t=0 es singular regular, puesto que a*(t) = y b*(t) = t2 son analíticas (en todo R).
2 2
1 1 1
Como a*0=a*(0) y b*0=b*(0) , el polinomio indicial es r(r–1)+ r +0 = r(r– ) y por tanto r1= y r2=0 ,
2 2 2
con r1–r2 no entero. Las dos soluciones linealmente independientes son, pues:
∞ ∞
x1 = ∑ ck tk+1/2 , c0≠0 y x2 = ∑ bk tk , b0≠0 , series convergentes para todo t .
k=0 k=0
Llevando la primera de ellas a la ecuación inicial:
∞ ∞ ∞
1 1 1
∑ 2 ( k+ )(k– ) ck tk–1/2 +
2 2 ∑ ( k+ 2 ) ck tk–1/2 + ∑ ck tk+3/2 = 0
k=0 k=0 k=0

(obsérverse que ahora todas las series comienzan por k=0 ya que, a diferencia de los
puntos regulares, no se anulan los primeros términos al derivar).
Igualando a 0 los coeficientes de las diferentes potencias de t :
t–1/2 : [ 2 (1/2)(–1/2) + 1/2 ] c0 = 0.c0 = 0 y c0 queda indeterminado como debía.
t1/2 : [ 2 (3/2)(1/2) + 3/2 ] c1 = 0 → c1 = 0 .
…………………………………………………
1 1 1 1
tk–1/2 : [2 (k+ )(k– ) + (k+ ) ] ck + ck-2 = 0 → ck = – c , k=2,3,…
2 2 2 k(2k+1) k-2
1 1 1
Por tanto: c3 = c5 = … = 0 , c2 = – c , c4 = – c = c , … y la primera solución es :
2.5 0 4.9 2 2.4.5.9 0
∞
1/2 (–1)m 2m
x1 = t ∑
[ 1 + m=1 2.4...2m.5.9...(4m+1)
t ] (eligiendo c0=1 )

∞ ∞ ∞
Para la segunda raiz del polinomio indicial: ∑ 2k(k–1) bk tk–1 + k=1
∑ k bk tk–1 + k=0
∑ bk tk+1 = 0 →
k=2
1
t0 : b1 = 0 , t1 : [4+2] b2 + b0 = 0 → b2 = – b0
6
…………………………………………………………
1
tk–1 : [2k(k–1)+k] bk + bk-2 = 0 → bk = – b , k=2,3,… →
k(2k–1) k-2
∞
1 1 (–1)m 2m
b3 = b5 = … = 0 , b4 = – b = b ,… →
4.7 2 2.4.3.7 0
x2 = 1 + ∑ 2.4...2m.3.7...(4m–1)
t
m=1

El criterio del cociente prueba que, como debían, ambas series convergen ∀t . La x2 es válida ∀t ,
pero x1 sólo si t>0 (en t=0 ni siquiera es derivable). Una x1 válida ∀t≠0 es x1 =|t|1/2 [1+∑] .

48
Más complicadas son las cuentas para los otros dos casos del teorema:

1 1
Ej 3. t 2 x " +2t2 x ' + (t2+ )x = 0 t=0 singular regular; a*(t)=2t , b*(t) = t2+ analíticas en R .
4 4
∞ ∞
1 1
r(r–1)+
4
=0 → r = doble → x1 =
2 ∑ ck tk+1/2 → ∑ [ k2cktk+1/2+(2k+1)cktk+3/2+cktk+5/2 ] = 0 →
k=0 k=0
2k–1 1 1 1 k1
c1=–c0 , ck=– c – c , k=2,3,… → c2= c0 , c3=– c0 , … , ck=(–1) c → x1 = t1/2 e–t
k 2 k-1 k 2 k-2 2 6 k! 0
∞
Como la raíz es doble, la otra solución necesariamente contiene un logaritmo: x2 = t 3/2 ∑ b ktk + x1 ln t →
k=0
∞ ∞ 1 2
3 k+1/2 1 3 1
k–1/2
x2' = ∑ (k+ 2 ) bkt + x + x1' ln t , x2" =
t 1 ∑ (k+ 2 )(k+ 2 ) bkt – t2 x1 + t x1' + x1" ln t →
k=0 k=0
∞
1
∑ [(k2+2k+1) bktk+3/2 +(2k+3) bktk+5/2 +bktk+7/2] + ln t [ t2x1"+2t2x1'+ (t2+ )x ] = 0
4 1
k=0

El último corchete es 0 por ser x1 solución (lo que acompaña al logaritmo siempre se anula).

→ b0 = b1 = b2 =…= 0 → x2 = t1/2 e–t ln t

1
(Para comprobar el resultado podríamos hacer y = et x → t2 y "+ y = 0 (Euler) → y = c1 t1/2 + c2 t1/2 ln t ;
4
también podríamos, una vez hallada la x1 , haber calculado otra solución utilizando la fórmula de 2.2:
e–2t
x2 = t1/2 e–t ∫ –2t dt = t1/2 e–t ln t , exactamente la misma solución hallada con las series].
te

Ej 4. t 2 x " +2t2 x ' – 2 x = 0 t=0 es singular regular, a*(t)=2t y b*(t)=–2 analíticas en R .

El polinomio indicial r(r–1) +0.r – 2 tiene por raíces r1=2 y r2=–1 . Así pues:
∞ ∞ ∞ ∞
x1 = ∑ ck tk+2 , c0≠0 → ∑ (k+2)(k+1) ck tk+2 + ∑ 2(k+2) ck tk+3 – k=0
∑ 2ck tk+2 = 0
k=0 k=0 k=0

2(k+1)
→ c0 indeterminado , ck = – c , k=1,2,… →
k(k+3) k–1
3 4 2(k+1) 2k 2(k–1) … (–2)k(k+1)
c1=– c0 , c2= c0 , c3=– c0 , … , ck = (–1)k c0 = 6c0
5 15 k(k+3) (k-1)(k+2) (k-2)(k+1) (k+3)!
∞ (–2) k (k+1) ∞
1 (–2)k(k+1)(k+2) k+1
Por tanto, eligiendo c0 = , x1 = ∑ (k+3)! tk+2 → x1' = ∑ t
6 k=0 k=0 (k+3)!
∞
La segunda solución (caso c] del teorema) es x2 = ∑ bktk–1 + a x1 ln t , b0≠0 , a constante (quizás nula) →
k=0
∞
∑ [ (k–1)(k–2)bk tk–1+2(k–1)bktk–2bktk–1 ] + a [ (–1+2t) x1 + 2t x1' ] + a ln t [ t2x1"+2t2x1'–2x1 ] = 0
k=0

Como siempre, el tercer corchete se anula, por ser x1 solución. Sustituyendo las series de x1 y x1' escritas
arriba en el segundo corchete y agrupando potencias de t :
a 2a 2
–2b0 –2b1 –2b2 t + [2b3+2b2–2b3 – + ] t + … = 0 → b1=–b0 , b2=0 , a=0 ; b0 , b3 indeterminados.
6 3
Como a =0 , en la expresión de x2 no aparece el ln t . Sabíamos que debía ser b0≠0 . El hecho de que también
b3 quede indeterminado se debe a que proporciona potencias t2 , comienzo de la serie de x1 . Elegimos b0=1
y b3=0 (para no recalcular x1 ). Como en la regla de recurrencia cada bk depende de bk–1 es b4=b5=…=0 .
1
Concluimos que: x2 = (1– t) [es fácil comprobar que satisface la ecuación].
t

∫ (1–t)
1–t t2
e –2t
[De esta x2 obtendríamos otra con: x1 = 2 dt , pero ahora la primitiva no es calculable].
t

[En este ejemplo, si, en vez de partir de la raíz mayor de la ecuación indicial, hubiésemos sustituido la x2 ,
habríamos obtenido las dos series de un tirón; pero esto ocurriría por que casualmente resulta ser a=0 ; si
fuera a≠0 sólo obtendríamos la solución trivial y=0 y deberíamos empezar de nuevo desde el principio].

49
 Para distinguir en el caso c] del teorema de Frobenius si aparecen logaritmos o no (es decir, si es o no a≠0 ) no es
necesario hallar la expresión del término general, bastan con los primeros términos de la x1 . Por ejemplo:

Ej 5. Estudiemos cuántas soluciones analíticas linealmente independientes tiene: 2 t x" + ( t –2) x ' + 3 x = 0 .
2 t 3t
t x" + t ( –1) x ' + x = 0 → t =0 es singular regular con r (r –1)–r +0 = 0 → r1=2 , r2=0 .
2 2
∞
Seguro que es analítica x1 = ∑ cktk+2 →
k=0
∞ ∞
∑ [ 2(k+2)(k+1)cktk+1 – 2(k+2)cktk+1 +(k+2)cktk+2 +3cktk+2 ] = k=0
∑ [ 2(k+2)k cktk+1+ (k+5) cktk+2 ] =0.
k=0

Calculemos, por ejemplo, los tres primeros términos de esta primera serie:
5 3 5
t1 : 0 c0 =0 , c0 indeterminado ; t2 : 6c1+5c0 =0 , c1 =– c0 ; t3 : 16c2+6c1 =0 , c2 =– c1 = c
6 8 16 0
5 5 4 … 5 5
→ x1 = t2 – t3 + t + [ → x1' = 2t – 2 t2 + 4 t3 + … ]
6 16
∞
La segunda solución x2 = ∑ bktk + a x1 ln t será analítica si a=0 y no lo será si a≠0 . Hay que trabajar:
k=0
∞ ∞
a 2a a
x2' = ∑ k bktk–1 + a x1' ln t + x
t 1
; x2" = ∑ k(k–1) bktk–2 + a x1" ln t + x '– x →
t 1 t2 1
k=1 k=2
∞ ∞ ∞
4a
∑ k(k–1) bktk–1 + k=1
∑ [ k b k tk – 2k bktk–1] + k=0
∑ 3bktk + 4ax1' – x + ax1 = 0 →
t 1
k=2

3 3
t0 : –2b1+3b0 = 0 → b1 = b0 ; t1 : 4b2+b1–4b2+3b1+8a – 4a = 0 → a = –b1 = – b0 ≠ 0
2 2
3 3
Por tanto, la segunda solución contiene el ln t y no es analítica en t=0 : x2 = 1 + t+…– x ln t .
2 2 1

Ej 6. t x" + 2 e t x ' = 0 Sabemos resolver la ecuación siguiendo la sección 2.2, pero utilicemos Frobenius:

t=0 es singular regular [ a*(t) =2 e t y b*(t)≡0 analíticas en todo R ]. λ(λ–1)+2λ=0 → r1=0 , r2=–1 .
∞ ∞
La x1 = ∑ ck tk se ve (¡a ojo!) que es x1 ≡ 1 . La otra es: x2 = ∑ bk tk–1 + a ln t →
k=0 k=0
∞ ∞ a
a
x2' = ∑ (k–1)bktk–2 + t
, x2" = ∑ (k–1)(k–2)bktk–3 – t2 →
k=0 k=0
1
2b0 t–2 +2b3 t +… – at–1 + [2+2t+t2+ t3+ …][ at–1–b0 t–2 +b2 +2 b3 t +… ] = 0 →
3
t–2 : 2b0–2b0 = 0 → b0 indeterminado como debía.

t–1 : –a+2a–2b0 = 0 → a = 2b0 (aparecen, pues, logaritmos).

1 3
t0 : 2a–b0+2b2 = 0 → b2 = b0 – a = – b0
2 2
1 2
t1 : 2b3+a– b0+2b2+4b3 = 0 → b3 = b0
3 9
1 3 2 2
……………………………………………………… → x2 = 2 ln t + – t + t +…
t 2 9
2et
Resolvamos la ecuación ahora sin series: x ' =v → v' =– v →
t
–∫ 2 t–1et dt ∫ (–2/t – 2 – t – t2 /3 +…) dt 2 3
v=Ce =Ce = C t–2 e– 2t – t /2 – t /9 +… =
1 1 1 1 2 1 3
= C t–2 [ 1 + (–2t– t2– t3+…) + (–2t– t2–…) + (–2t–…) + … ] →
2 9 2 2 6

∫ [t–2 – 2t–1 + 23 – 9 t + … ] dt = K – C [ 2 ln t + t
4 1 3 2
x = K+ C – t + t2 +… ]
2 9

50
3.4 Ecuaciones de Legendre, Hermite y Bessel.

La ecuación de Legendre es la ecuación [L] (1–t 2 ) x" – 2tx' + p(p+1)x = 0

Resolvámosla primero en torno a t=0 que es punto regular. Como a(t)= –2t/(1–t2) y b(t)= p(p+1)/(1–t2) son
analíticas en |t |<1 la ecuación tiene soluciones analíticas al menos en ese intervalo. Probamos pues:
∞ ∞ ∞ ∞
(p–k+2)(p+k–1)
x = ∑ c k tk → ∑ [ k ( k–1)cktk–2–k(k–1)cktk] – ∑ 2kc k tk + ∑ p(p+1)cktk = 0 → ck = –
k(k–1)
ck-2 , k=2,3,…
k=0 k=2 k=1 k=0

p(p+1) (p–1)(p+2) p(p–2)(p+1)(p+3) (p–1)(p–3)(p+2)(p+4)

→ c2 = – 2.1 c0 , c3 = – 3.2 c1 , c4 = 4! c0 , c5 = 5! c1 ,… →
∞
n p(p–2)…(p–2n+2)(p+1)(p+3)…(p+2n–1) 2n
x1 = 1 + ∑ (–1) (2n)! t
n=1
∞ (p–1)(p–3)…(p–2n+1)(p+2)(p+4)…(p+2n) 2n+1
x2 = t + ∑ (–1)n (2n+1)! t
n=1

Si p es un entero par positivo, p=2m , x1 se convierte en un polinomio de grado 2m :

35
p=0 → x1 = 1 , p=2 → x1 = 1–3t2 , p=4 → x1 = 1–10t2 + 3 t4 ,…

Si p es impar, p=2m+1 , es x2 quien se reduce a un polinomio de grado 2m+1 :

5 14 21
p=1 → x2 = t , p=3 → x2 = t – 3 t3 , p=5 → x2 = t – 3 t3 + 5 t5 , …

Se llama polinomio de Legendre de grado n al polinomio Pn solución de [L] con p=n∈N , Pn(1)=1 :
3 1 5 3 35 15 3 63 35 15
P0 =1 , P1 = t , P2 = 2 t2 – 2 , P3 = 2 t3 – 2 t , P4 = 8 t4 – 4 t2+ 8 , P5 = 8 t5 – 4 t3+ 8 t , …

n P2
Como Pn(–t) = (–1) Pn(t) , los P2m tienen simetría par y los P2m+1 impar. P0
Observemos que los P2m+1 y las derivadas de los P2m se anulan en 0 .
Se pueden probar además las siguientes propiedades de los Pn : P1
P3
1 dn 2
Pn(t) = n (t –1) n , fórmula de Rodrigues.
2 n! dtn
-1 1
Pn tiene n ceros reales, todos en (–1,1) . Los Pn son ortogonales:
1 1 2
∫–1 PnPm dt = 0 ∫–1 Pn
2
, si m≠n ; dt = 2n+1

Para las aplicaciones de la ecuación [L] a las EDPs se necesitará saber qué
soluciones están acotadas en ±1 . Para estudiarlo resolvemos la ecuación en torno a t=1 , haciendo s= t–1 :
[L1] s(s+2)x"+ 2(s+1)x' – p(p+1)x = 0
Para ella s=0 es punto singular regular, con a*(s)=2(s+1)/(s+2) , b*(s)=–p(p+1)s/(s+2) analíticas para |s|<2 .
El polinomio indicial tiene raíz r=0 doble ∀p . Dos soluciones linealmente independientes de [L1] son:
∞ ∞
x1 = ∑ cksk y x2 = |s| ∑ bksk + x1 ln|s| , c0=1
k=0 k=0
y las series convergen al menos para |s|<2 . Sin necesidad de calcular ningún coeficiente podemos ya
afirmar que la x1 siempre está acotada en s= 0 ( t= 1 ), mientras que la x2 no lo está (tiende a –∞ si s→0 ).
Vamos a calcular x1 y comprobar que si p=n obtenemos de nuevo los Pn [ pues x1(1)=1 ] . Debe ser:
∞ ∞ ∞
p(p+1)–k(k–1)
∑ [ k ( k–1)cksk+2k(k–1)cksk–1] + k=1
∑ [2kcksk+2kcksk–1] – k=0
∑ p(p+1)cksk = 0 → ck = 2k
2 ck-1 , k=1,2,…
k=2
p(p+1) p(p+1)[p(p+1)–2.1] 2 p(p+1)[p(p+1)–2.1]…[p(p+1)–n(n–1)] n
→ x1(s) = 1 + s + s +…+ s +…
2 16 2n(n!)2
Si p=n la regla de recurrencia nos dice que el cn+1 y todos lo siguientes se anulan. En particular:
6[6-2] 3 1
p=0 → x1 = 1 ; p=1 → x1 = 1+s = t ; p=2 → x1 = 1+3s+ 16 s2 = 2 t2 – 2 ; …

Se prueba que si p≠n la x1 no está acotada cuando s → –2 ( t → –1 ) con lo que, salvo constantes, las
únicas soluciones de [L] acotadas a la vez en t=1 y t=–1 son los polinomios de Legendre.

51
Las series solución de [L] hablan sólo de lo que ocurre para t∈(–1,1) . Las de [L1] de t∈(–1,3) . Para analizar
lo que ocurre 'en el infinito' hacemos en [L] el cambio s=1/t y estudiamos s=0 (punto del infinito de [L]):
dx 2 2
= –s 2 dx ; d x = s4 d x +2s3 dx → [L ] s2 (s2 –1)x"+ 2s3 x' + p(p+1)x = 0
dt ds dt 2 ds2 ds ∞

s=0 es singular regular, con a*(s)=2s2/(s2–1) y b*(s)=p(p+1)/(s2–1) analíticas en |s|<1 . Las series solución
de [L∞] convergen al menos en ese intervalo y por tanto de ellas podemos extraer información sobre las
soluciones de [L] para |t |>1 . Por ejemplo, como el polinomio indicial de [L∞] es r2 – r– p(p+1) (de raices 1+p
y –p ) y para todo p es r1=1+p>0 deducimos que siempre hay soluciones de [L] que tienden a 0 si t→∞.
∞ ∞
(Pues x1∞(s) = s1+p ∑ cksk → 0 si s→0+ ; o sea, x1∞(t) = t–(1+p) ∑ ckt–k → 0 si t→∞ ).
k=0 k=0

Resolvamos por series [L∞] en el caso más sencillo: p=0 (único valor para el que s=0 es regular):
∞
k–2 1 1 1 1
x = ∑ cksk → ck = ck-2 , k=2,3,… → x = c0 +c1 [ s + s3 + s5 +… ] = c0 +c1 [ t–1 + t–3 + t–5+ … ] ,
k=0 k 3 5 3 5
serie (no de potencias) que nos describe las soluciones de [L] para |t |>1 , valores para los que converge.
–1 1 1+t 1 1+s
[De otra forma: (1–t2)x" +2tx ' =0 → x' =c1(1–t2) → x = c0 +c1
2
ln|
1–t
| = c0 +c1 2 ln| 1–s | , válidas si t,s≠1 ].

Otra ecuación que aparece al abordar problemas físicos es la de Hermite: [H] x " – 2tx' + 2px = 0 .
∞
Como t=0 es regular, existe solución analítica x = ∑ cktk , convergente además en todo R. Entonces
k=0
∞ ∞ ∞
(k–2–p)
∑ k ( k–1)cktk–2 – k=1 ∑ 2pcktk = 0
∑ 2kc k tk + k=0 → ck = 2 c
k(k–1) k–2
, k=2,3,… →
k=2

∞ ∞
n (–p)(2–p)…(2n–2–p) 2n n (1–p)(3–p)…(2n–1–p) 2n+1
x = c1 x1 + c2 x2 = c1 [ 1 + ∑ 2 (2n)! t ] + c2 [ t + ∑ 2 (2n+1)! t ]
n=1 n=1

Como ocurría con Legendre, [H] posee solución polinómica cuando p= n∈N . Si p=2m , la x1 pasa a ser un
polinomio de grado 2m , y si p=2m+1 es la x2 la que se convierte en un polinomio de ese mismo grado:
2
p=0 → x1 = 1 ; p=1 → x2 = t ; p=2 → x1 = 1 –2t2 ; p=3 → x2 = t – 3 t3 ; …

Se llaman polinomios de Hermite Hn(t) a las soluciones polinómicas de [H] tales que los términos que
contienen la potencia más alta de t son de la forma 2n tn , es decir, a:

H0 = 1 ; H1 = 2 t ; H2 = 4 t2 – 2 ; H3 = 8 t3 – 12t ; …

Citemos ahora, también sin demostrarlas, una serie de propiedades de los Hn que serán útiles, por ejemplo,
cuando aparezcan en física cuántica. Una forma de generarlos todos nos lo proporciona el hecho de que:
2 ∞
1
e2ts–s = ∑ Hn(t) sn (a la exponencial de la izquierda se le llama función generatriz de los Hn )
k=0 n!
[Nos limitamos a comprobarlo para los 4 que ya hemos calculado:
4t3 3 1 4t3
( 1 +2ts +2t2s2 +
3
s +… ) (1– s2 + 2 s4 – … ) = 1 + 2t s + (2t2–1) s2 + ( 3 –2t) s3 + … ].
2 dn 2
De la función generatriz se deduce una fórmula de Rodrigues similar a la de [L]: Hn(t) = (–1)net n e–t .
dt
∂ 2ts–s2 t 2 ∂ –(t–s)2 ∂ ∂ n t2 d –z2
[Pues Hn(t) = ∂sn e |s=0 = e ∂sn e |s=0 = ( t–s=z , ∂s = – ∂z ) = (–1) e dzn e |z=t ].
2
[H] no aparece inicialmente en cuántica, sino u"+ (2p+1–t2 )u = 0 . Haciendo u= xe–t /2 en ella se llega a [H].
No es fácil probar (no nos lleva a nada el análisis del punto del infinito de [H], por ser singular no regular) que
2
las únicas soluciones u de la inicial que tienden a 0 si |t |→∞ son las de la forma un(t)= e –t /2 Hn(t) ,
llamadas funciónes de Hermite de orden n. Sólo estas un interesarán físicamente.
Como los Pn de Legendre, se puede ver que también las un son ortogonales, ahora en (–∞,∞) :
∞ ∞ ∞ ∞
∫–∞ un um dt = ∫–∞ Hn Hm e–t ∫–∞ un2 dt = ∫–∞ Hn2 e–t dt = 2n n! √
π
2 2
dt = 0 , si m≠n ;

∞ ∞ ∞ ∞ ∞ 2 2 ∞ ∞
[Lo comprobamos si n=0,1 : ∫–∞ u0 u1 = ∫–∞ 2te–t = 0 , ∫–∞ u02 = ∫–∞ e–t =
√ π , ∫–∞ u1 = –2te–t |–∞ + ∫–∞ 2e–t =2
√π ]
2 2 2

52
 La última ecuación de interés físico que vemos es la de Bessel: [B] t 2 x" + t x' + (t2 –p 2 ) x = 0 , p≥0 .

Para ella, t=0 es punto singular regular con polinomio indicial r2–p2 de raíces r1=p , r2=–p .
∞
Hay pues una solución (acotada en t=0 ∀p ) de la forma x1 = tp ∑ c k tk , t>0 , serie convergente en todo R .
k=0
Llevándola a [B]:
∞ c
∑ [k(2p+k)cktp+k+cktp+k+2] = 0 → c1=0 ; ck = – k(2p+k) ,
k–2
k=2,3,… →
k=0
∞
c0 c0 (–1)m t2m
c1 = c3 = …=0 ; c2 = – 2 ; c4 = 4 ; … → x1 = c0 tp [ 1 + ∑ ]
2 (p+1) 2 2(p+1)(p+2) m=1 22mm! (p+1)…(p+m)

6 Podemos dar una expresión más corta de la x1 utilizando las propiedades de la

función gamma (función que generaliza el factorial para números no enteros)
Γ definida por la integral impropia convergente:
∞
Γ(s) = ∫0 e–x xs–1dx si s>0 , y extendida a los s<0 mediante:
Γ(s+n)
2 Γ(s) = (s+n–1)...(s+1)s si –n<s<–n+1 , n∈N .
1
Se cumplen para la Γ las siguientes igualdades:
– 2 -1 1 2 3 4
∞ ∞
 π ; Γ(s+1)= –e–x xs|∞0 +sΓ(s) = sΓ(s)
Γ(1)= ∫0 e–x dx =1 ; Γ(1/2)=2∫ 0 e–u du = √
2

→ Γ(s+n) = (s+n–1)…(s+1)s Γ(s) → Γ(n+1) = n! , n∈N

Eligiendo c0 = 1 / [2p Γ(p+1)] obtenemos entonces la siguiente solución de [B]:

∞
(–1)m
J p (t) ≡ [ 2t ]p ∑ m! Γ(p+m+1)
[ 2t ] 2m función de Bessel de primera especie y orden p
m=0

En particular se tiene:
∞ ∞
(–1) m t 2m (–1)m
J0(t) = ∑ 2 [ 2
] , J1(t) = ∑ [ t ]2m+1
m=0 (m!) m=0 m! (m+1)! 2

cuyás gráficas resultan ser las de la derecha. Se prueba

que, al igual que J0 y J1 , todas las Jp son oscilatorias y
que para t grande se parecen a:

√ πt cos[ t –(2p+1)4 ]
π
Jp∼
2

[el punto del infinito es no regular y de nuevo no es útil].

Cada Jp tiene un número infinito de ceros en (0,∞) [que deben conocerse para resolver algunas EDPs]:
los de J0 son: 2.4048, 5.5201, 8.6532, … ; los de J1 : 3.8317, 7.0156, 10.1735, …

Para calcular una segunda solución linealmente independiente de [B] (necesariamente no acotada en t=0 ),
≠ 0,1,… existe una x2 de la forma:
el método de Frobenius nos dice que si r1–r2 = 2 p≠
∞ ∞
t –p (–1)m t 2m
x2 = t–p ∑ b k tk , t>0 (llevándola a [B] se tiene J–p(t) ≡ [ 2 ] ∑ [ ] )
k=0 m=0 m! Γ(–p+m+1) 2

Si p∉Z , pero 2p∈Z ( p= 1 / 2 ,3 / 2 ,5 / 2 , … ) , en principio la x2 podría contener un lnt pero se puede ver
que no es así (caso c] de Frobenius, pero con a=0 ). De hecho, haciendo p=1/2 en J±p se obtienen:

√
∞

√ √
2 (–1)m t2m+1
∑
2 2
J1/2(t) = = s e n t , J–1/2(t) = … = cos t ,
22m+1
t m! (m+1/2)…1/2 Γ(1/2)
m=0 πt πt
1
que son soluciones independientes de [B] con p= 2 (la expresión asintotica de arriba es exacta para J1/2 ).
2p
Además probaremos una relación que expresa Jp+1 en función de los anteriores: J p+1 = t J p – J p–1 ,

con lo que todas las J(2n+1)/2 , n∈Z son funciones elementales (las demás Jp no lo son).

53
Para p= n = 0,1,… el atajo anterior no resulta, pues es fácil comprobar que si cambiamos n por –n la J–n
resultante no es linealmente independiente de la Jn [se tiene que J–n = (–1)nJn ]. Así que nosotros nos
veríamos obligados a calcular las x2 de Frobenius (y comprobaríamos que aparecen logaritmos en su larga
expresión). Por ejemplo, para p=0 (que seguro contiene logaritmos) se acaba obteniendo:
∞
(–1)m+1 1 1 t 2m
x2(t) = ∑ 2 [1 + + … + ][ ]
2 m 2
+ J0(t) ln t ≡ K0(t) , t >0
m=0 (m!)
función de Bessel de segunda especie y orden 0
Pero en muchos problemas físicos en los que surge la ecuación [B] es necesario que las soluciones estén
acotadas, y para ellos no servirá de nada el conocimiento de estas complicadas segundas soluciones.

Lo que sí puede ser útil en el futuro será conocer las siguientes propiedades de las derivadas de las Jp :
d d
dt
[ tp Jp(t) ] = tp Jp–1(t) , dt [ t– p Jp(t) ] = – t– p Jp+1(t) (En particular es [ t J 1 ] ' =J0 , [ J 0 ] ' =– J1 ) .

∞ ∞
d (–1)m t2m+2p (–1)m t2m+p–1
(Son inmediatas: dt ∑ 2m+p = tp ∑ 2m+p–1 y similar la otra).
m=0 2 m! Γ(p+m+1) m=0 2 m! Γ(p+m)
p p
De ellas se deduce, derivando y despejando la Jp' en las dos, que: Jp' = Jp–1 – t Jp = –Jp+1 + t Jp ,

y de las dos últimas igualdades deducimos la relación de recurrencia de las Jp ya citada.

54