Análisis III B - Turno mañana 1

Teorema de Fuchs
Hemos visto que si P (z) y Q(z) son funciones analı́ticas, entonces las soluciones de la
ecuación y 00 + P (z) y 0 + Q(z) y = 0 son analı́ticas. Parecerı́a razonable esperar que si P (z) y
Q(z) tienen un polo en z0 , las soluciones de la ecuación diferencial también tengan un polo en
ese punto.
Es fácil ver que esto no es cierto en general ya que, por ejemplo, una solución de la ecuación
y + z12 y 0 − z13 y = 0 es la función −z e1/z . El desarrollo de Laurent de dicha función en un
00

entorno a z = 0 tiene infinitos términos con potencias negativas, o sea tiene una singularidad
esencial en z = 0.

Sea
A(z) y 00 + B(z) y 0 + C(z) y = 0
una ecuación diferencial homogénea de segundo orden cuyos coeficientes son funciones analı́ticas
en una cierta región. Dividiendo por A(z) ambos miembros de la ecuación anterior se tiene una
ecuación de la forma

y 00 + p(z) y 0 + q(z) y = 0. (1)

Decimos que z0 ∈ C es un punto singular de la ecuación (1) si lo es de p y/o de q.

En muchos casos interesa obtener soluciones de la ecuación (1) definidas en 0 < |z −zo | < R,
donde z0 es un punto singular, y de la forma
∞
X
zr an z n , (2)
n=0

con a0 6= 0, y r ∈ C.

Un punto z0 tal que sea una singularidad de la ecuación y tal que exista una solución de la
forma 2 se dice que es un punto singular regular o un punto de Fuchs.

Teorema de Fuchs: z0 es un punto singular regular o un punto de Fuchs de la ecuación

(1) si y sólo si

1. z0 es un punto singular de p y/o q y

2. (z − z0 ) p(z) y (z − z0 )2 q(z) son analı́ticas en z0 .

Consideremos, por simplicidad, z0 = 0. De lo anterior, p(z) tiene a la sumo un polo simple

y q(z) a lo sumo un polo de orden dos en z0 , por la tanto podemos escribirlas como:

∞ ∞
A X B C X
p(z) = + αn z n y q(z) = + + βn z n . (3)
z n=0
z2 z n=0
2 Análisis III B - Turno mañana

∞
P
Buscamos soluciones de la ecuación (1) de la forma z r an z n , con a0 6= 0, y r ∈ C.
n=0

Podemos considerar, por la linealidad de la ecuación, que a0 = 1 y ası́ la solución toma la

forma

y(z) = z r + a1 z r+1 + · · · + an z r+n + . . . (4)

Si reemplazamos (3) y (4) en la ecuación diferencial (1) tenemos:

∞
! 
P∞ A X
n
∞
P
(n + r) (n + r − 1) an z n+r−2 + + αn z (n + r) an z n+r−1 +
n=0 z n=0 n=0
∞
! 
B C X ∞
P
+ 2
+ + βn z n an z n+r = 0
z z n=0 n=0

o equivalentemente
[r (r − 1) z r−2 + a1 (r − 1) r z r−1 + · · · + an (r + n) (r + n − 1) z r+n−2 + . . . ]+
 
A
+ + α0 + α1 z + . . . [r z r−1 + a1 (r + 1) z r + · · · + an (r + n) z r+n−1 + . . . ] +
z
 
B C
+ 2 + + β0 + β1 z + . . . [z r + a1 z r+1 + · · · + an z r+n + . . . ] = 0.
z z

Multiplicando las series tenemos

r (r − 1) z r−2 + a1 (r − 1) r z r−1 + · · · + an (r + n) (r + n − 1) z r+n−2 + · · · +

+A r z r−2 +A a1 (r +1) z r−1 +· · ·+A an (r +n) z r+n−2 +· · ·+α0 r z r−1 +α0 a1 (r +1) z r +· · · +

+an (r + n) z r+n−1 + · · · + α1 r z r + α1 a1 (r + 1) z r+1 + · · · + α1 an (r + n) z r+n + · · · + B z r−2 +

+B a1 z r−1 +B a2 z r +· · ·+B an z r+n−2 +· · ·+C z r−1 +C a1 z r +C a2 z r+1 +· · ·+C an z r+n−1 +

+ · · · + β0 z r + β0 a1 z r+1 + β0 a2 z r+2 + · · · + β0 an z r+n + · · · = 0.

Reodenando adecuadamente llegamos a

r (r − 1) z r−2 + a1 (r − 1) r z r−1 + · · · + an (r + n) (r + n − 1) z r+n−2 + · · · + A r z r−2 +

+[A a1 (r + 1) + α0 r] z r−1 + [A a2 (r + 2) + α0 a1 (r + 1) + α1 r] z r + · · · + B z r−2 +

(B a1 + C) z r−1 + (B a2 + C a1 + β0 ) z r + · · · = 0.
Análisis III B - Turno mañana 3

O sea, podemos escribir la ecuación de la forma

∞
X
r−2
z bn z n = 0,
n=0

lo que implica que todos los coeficientes bn son iguales a cero. Entonces

b0 = r (r − 1) + A r + B = 0. Luego si llamamos polinomio indicial al polinomio I(r) =

r (r − 1) + A r + B, obtenemos a partir de sus raı́ces los valores de r que permiten obtener las
soluciones de la forma (3).

b1 = (r+1) r a1 +A (r+1) a1 +α0 r+B a1 +C = [r (r + 1) + A (r + 1) + B ] a1 +(α0 r+C) = 0

| {z }
I(r+1)

Luego, I(r + 1) a1 + α0 r + C = 0.

En general, si se anula el coeficiente de z n+r−2 se obtiene:

bn = an (r + n) (r + n − 1) + an (r + n) A + Pn + an B + Qn =,

= [(r + n) (r + n − 1) + (r + n) A + B ] an + Pn + Qn = 0
| {z }
I(r+n)

donde Pn es función de a1 , ..., an−1 y de α0 , α1 ..., αn−1 y Qn es función de a1 , ..., an−1 y de

C, β0 , β1 , ..., βn−2 .

Luego la fórmula de recurrencia puede escribirse en la forma I(r + n) an + Tn = 0 (4). A

partir de esta ecuación se pueden determinar los coeficientes por recurrencia, siempre que la
ecuación se pueda resolver para n ≥ 1 y ello será posible si I(r + n) 6= 0, n ≥ 1.

Tenemos tres casos distintos:

1. Primer caso: I(r) tiene dos raı́ces r1 y r2 distintas tales que su diferencia no es un número
entero.
Entonces, ∀ n ≥ 1, r1 + n 6= r2 y r2 + n 6= r1 , lo que implica que para n ≥ 1 se tiene
que I(r1 + n) 6= 0 , I(r2 + n) 6= 0; luego para ambos valores de ecuación siempre puede
resolverse y se obtiene dos soluciones:

y1 (z) = z r1 [1 + a1 z + ... + an z n + ...] = z r1 f1 (z)

y2 (z) = z r2 [1 + b1 z + ... + bn z n + ...] = z r2 f2 (z)

4 Análisis III B - Turno mañana

siendo f1 y f2 analı́ticas en 0 y f1 (0) = f2 (0) = 1.

Estas dos soluciones son linealmente independientes y por lo tanto la solución general de
la ecuación es
y(z) = C1 z r1 f1 (z) + C2 z r2 f2 (z).

Se dice entonces que, en ese punto, la ecución es del primer tipo de Fuchs.

2. Segundo caso: I(r) tiene una raı́z doble r.

En este caso hay una solución de la forma

y1 (z) = z r f (z)

Una segunda solución es de la forma

y2 (z) = y1 (z) g(z),

donde g(z) se puede determinar reemplazando y2 en la ecuación diferencial.

La ecución se dice que es del segundo tipo de Fuchs.
Ecuación del segundo tipo de Fuchs

3. Tercer caso: I(r) tiene dos raı́ces distintas cuya diferencia es un entero, es decir, las raı́ces
son r1 y r2 = r1 + p, para algún p ∈ N. En este caso I(r2 + n) 6= 0 para todo n ≥ 1.
Luego la ecuación (4) tiene solución para r = r2 . Por lo tanto, la ecuación admite una
solución de la forma

y2 (z) = z r2 (1 + a1 z + ... + an z n + . . . ).

Ahora, I(r1 +n) = 0 para n = p y no nulo para cualquier otro número natural. La ecuación
(4), toma la forma 0 ap + Tp = 0. Si Tp 6= 0, la determinación de ap es imposible y la
∞
P
ecuación diferencial tiene una sola solución de la forma y(z) = an z n+r y la ecuación
n=0
es del segundo tipo de Fuchs.
En cambio, si se tiene que Tp = 0, la ecuación admite como solución cualquier valor de
ap . Fijando uno de ellos arbitrariamente, se pueden determinar todos los valores de an y
P∞
se obtienen dos soluciones de la forma y(z) = an z n+r , o sea, la ecuación es del primer
n=0
tipo de Fuchs.

Veamos algunos ejemplos de soluciones alrededor de un punto singular:

Análisis III B - Turno mañana 5

Ejemplo 1

Resolver la ecuación x2 y 00 + (x2 + x) y 0 − y = 0.

Esta ecuación tiene al origen como punto singular. Escribamos la ecuación de la forma
   
00 1 0 1
y + 1+ y + − 2 y = 0.
x x
1 1
Entonces p(x) = + 1 y q(x) = − 2
x x
De acuerdo con (2) tenemos A = 1, α0 = 1 αi = 0 ∀i ≥ 1, B = −1, C = 0 y βj = 0 ∀j ≥ 0.

Luego el polinomio indicial es I(r) = r (r − 1) + r − 1 = r 2 − 1, cuyas raı́ces son r1 = −1,

r2 = 1.
r2 = r1 + p, ⇒ p = 2.
La raı́ces difieren en un entero, y como

Tp = [α0 (r1 + 1) + C] a1 + α1 r1 + β0 = [1 (−1 + 1) + 0] a1 + 0 (−1) + 0 = 0

entonces, la ecuación tiene una solución de la forma

∞
! ∞
X X
n
−1
y1 (x) = x 1+ an x =x −1
+ an xn−1
n=1 n=1

y otra de la forma

∞
! ∞
X X
n
y2 (x) = x 1+ bn x =x+ bn xn+1
n=1 n=1

Hallemos y10 y y100

∞ ∞
1 X 2 X
y10 (x) =− 2 + (n − 1) an xn−2 y100 (x) = 3+ (n − 1) (n − 2) an xn−3 .
x n=1
x n=1

Ahora reemplacemos en la ecuación diferencial

∞
! ∞
! ∞
2 X 1 X 1 X
x2 + (n − 1) (n − 2) an xn−3 +(x2 +x) − 2+ (n − 1) an xn−2
− − an xn−1 = 0
x3 n=1 x n=1
x n=1

Multiplicando y reaordenando adecuadamente obtenemos

6 Análisis III B - Turno mañana

∞ ∞ ∞ ∞
2 X n−1
X
n 1 X n−1 1 X
+ (n−1) (n−2) an x −1+ (n−1) an x − + (n−1) an x − − an xn−1 =
x n=1 n=1
x n=1
x n=1

∞ ∞ ∞ ∞
2 X n−1
X
n−1 1
X
n−1 1
X
= + (n−1) (n−2) an x −1+ (n−2) an−1 x − + (n−1) an x − − an xn−1 =
x n=1 n=2
x n=1 x n=1

∞ ∞ ∞ ∞
2 X n−1
X
n−1 1
X
n−1 1
X
= + (n−1) (n−2) an x −1+ (n−2) an−1 x − + (n−1) an x − −a1 − an xn−1 =
x n=2 n=2
x n=2 x n=2

∞
P
= −1 − a1 + [(n − 1) (n − 2) an + (n − 2) an−1 + (n − 1) an − an ] xn−1 =
n=2

∞
P
= −1 − a1 + [n (n − 2) an + (n − 2) an−1] xn−1 = 0.
n=2

Luego,

−1 − a1 = 0 y n (n − 2) an + (n − 2) an−1 = 0, entonces,
an−1
a1 = −1 y an = −
, n ≥ 2.
n
Si tomamos a2 = 0 tenemos que an = 0, n ≥ 2 y la solución es
1
y1 (x) = − 1.
x

Hacemos lo mismo para la otra solución:

∞
P ∞
P
y20 (x) = 1 + (n + 1) bn xn y200 (x) = (n + 1) n bn xn−1
n=1 n=1

Reemplazando en la ecuación diferencial tenemos:

∞ ∞
! ∞
X X X
2 n−1 2 n
x (n + 1) n bn x + (x + x) 1+ (n + 1) bn x −x− bn xn+1 =
n=1 n=1 n=1
∞
X ∞
X ∞
X ∞
X
n+1 2 n+2 n+1
= (n + 1) n bn x +x + (n + 1) bn x +x+ (n + 1) bn x −x− bn xn+1 =
n=1 n=1 n=1 n=1
∞
X ∞
X
= (n + 2) n bn xn+1 + x2 + (n + 1) bn xn+2 =
n=1 n=1
∞
X ∞
X
2 n+1 2
= 3 b1 x + (n + 2) n bn x +x + (n + 1) bn xn+2 =
n=2 n=1
Análisis III B - Turno mañana 7

∞
X ∞
X
2 n+2
= (3 b1 + 1) x + (n + 3) (n + 1) bn+1 x + (n + 1) bn xn+2 =
n=1 n=1
∞
X
2
= (3 b1 + 1) x + [(n + 3) (n + 1) bn+1 + (n + 1) bn ] xn+2 = 0.
n=1

Entonces,

−bn
b1 = −1/3 y (n + 3) (n + 1) bn+1 + (n + 1) bn = 0, o sea, bn+1 = , n ≥ 1.
n+3
−b1 1 2
b2 = = =
4 43 4!
−b2 2 2
b3 = =− =−
5 5 4! 5!
−b3 2 2
b4 = = =
6 6 5! 6!
−b4 2 2
b5 = = =−
7 7 6! 7!
..
.

(−1)n 2
bn = , n≥1
(n + 2)!

Luego, la otra solución nos queda

∞
X (−1)n 2 n+1 2 2 2
y2 (x) = x + x = x − x2 + x3 − x4 + .... =
n=1
(n + 2)! 3! 4! 5!
   
1 2 2 2 3 2 4 1 2 2 2 3 2 4
= x − x + x + ... = −2 + 2 x + 2 − 2 x + x − x + x + ... =
x 2! 3! 4! x 2! 3! 4!
∞
!
1 X (−1)n xn
= 2x− 2+2 .
x n=0
n!

Entonces,
2

y2 (x) = x − 1 + e−x .
x
Por lo tanto

x − 1 + e−x 1−x e−x 1−x

y(x) = C1 + C2 = A1 + A2 .
x x x x
8 Análisis III B - Turno mañana

Ejemplo 2

Resolver la ecuación x2 y 00 − 3 x y 0 + 4 y = 0.

El punto singular es el origen. Escribamos la ecuación de la forma

3 4
y 00 − + 2 y = 0.
x x
3 4
Entonces p(x) = − y q(x) = − 2
x x
De acuerdo con (2) tenemos A = −3, αi = 0 ∀i ≥ 0, B = 4, C = 0 y βj = 0 ∀j ≥ 0.

Luego el polinomio indicial es I(r) = r (r − 1) − 3 r + 4 = (r − 2)2 . r = 2 es raı́z doble.

Entonces, es del segundo tipo de Fuchs y tiene una solución de la forma

∞
X ∞
X
2 n 2
y1 (x) = x (1 + an x ) = x + an xn+2
n=1 n=1

Hallemos y10 y y100

∞
P ∞
P
y10 (x) = 2 x + (n + 1) an xn+2 y100 (x) = 2 + (n + 2) (n + 1) an xn .
n=1 n=1

Ahora reemplacemos en la ecuación diferencial

∞
X ∞
X ∞
X
2 n n+1 2
x (2 + (n + 2) (n + 1) an x ) − 3 x (2 x + (n + 2) an x ) + 4 (x + an xn+2 ) = 0
n=1 n=1 n=1

∞
X ∞
X ∞
X
2 n+2 2 n+2 2
2x + (n + 2) (n + 1) an x − 6x − 3 (n + 2) an x +4x +4 an xn+2 =
n=1 n=1 n=1

∞
X
2
= (2 − 6 + 4) x + [(n + 2) (n + 1) − 3 (n + 2) + 4] an xn+2 = 0
n=1

Luego,

[(n + 2) (n + 1) − 3 (n + 2) + 4] an = n2 an = 0,

entonces, an = 0, n ≥ 1.
Análisis III B - Turno mañana 9

Luego la solución es
y1 (x) = x2 .

Busquemos la otra solución

y2 (x) = u(x) y1(x) = u(x) x2

y20 (x) = u0 (x) x2 + u(x) 2 x y y200 (x) = u00 (x) x2 + 2 u0 (x) 2 x + u(x) 2

Reemplacemos en la ecuación diferencial

x2 (u00 (x) x2 + 2 u0(x) 2 x + u(x) 2) − 3 x (u0(x) x2 + u(x) 2 x) + 4 u(x) x2 = 0

u00 (x) x4 + 4 u0(x) x3 + u(x) 2 x2 − 3 u0(x) x3 − 6 u(x) x2 + 4 u(x) x2 = 0 ⇒

u00 (x) x2 + 4 u0(x) x + u(x) 2 − 3 u0(x) x − 6 u(x) + 4 u(x) = 0 ⇒

u00 (x) x2 + u0 (x) x = 0 ⇒ u00 (x) x + u0(x) = 0

k
Sea z = u0 , entonces z 0 x + z = 0, que tiene como solución z = .
x
k
Luego, u0 = ⇒ u(x) = k ln|x|.
x
Entonces la solución es y2 (x) = x2 ln|x|.

Por lo tanto,

y(x) = C1 x2 + C2 x2 ln|x|.