Capı́tulo 6

EDOs de segundo orden lineales

A partir de este capı́tulo y hasta el final del capı́tulo 12 nos dedicaremos a desarrollar
métodos especı́ficos para encontrar la solución general de la ecuación lineal de segundo
orden:
y′′ + p(x)y′ + q(x)y = g(x). (6.1)

De hecho no partimos de cero, ya sabemos que la solución general tiene la forma

y(x) = c1 y1 (x) + c2 y2 (x) + y p (x), donde c1 , c2 son dos constantes arbitrarias que deben
determinarse con las condiciones auxiliares, tı́picamente condiciones iniciales de la forma
y(x0 ) = y0 , y′ (x0 ) = y′0 . Ya sabemos que si conocemos y1 , y2 podemos usar el método de
variación de los parámetros para encontrar la solución particular y p . De manera que nos
centraremos en la manera de encontrar dos soluciones linealmente independientes y1 , y2
de la ecuación diferencial homogénea y′′ + p(x)y′ + q(x)y = 0. Como comprobaremos,
una de las técnicas más útiles es la de desarrollo en serie de potencias.
Una pregunta que naturalmente puede surgir es la siguiente: si ya hemos visto en
capı́tulos anteriores que cualquier EDO lineal se puede escribir como un sistema de
ecuaciones diferenciales de primer orden, que necesidad hay de dedicar no uno sino
hasta siete capı́tulos a la resolución de ecuaciones de segundo orden? La respuesta es
BLA BLA BLA FALTA FALTA FALTA
Otro aspecto a tener en cuenta es que no siempre que veamos una EDO de segundo
orden debemos automáticamente intentar un desarrollo en serie para las soluciones. La
cabeza sirve para pensar, y hay que evaluar primero si existe algún otro método que
sepamos usar mediante el cual podamos resolver dicha EDO sin necesidad de recurrir al
desarrollo en serie de potencias que por lo general es largo y tedioso.
Ası́ pues, una vez hecha la advertencia, en este capı́tulo empezaremos repasando
primero las series de potencias y sus propiedades. Seguidamente veremos como podemos
hallar la solución mediante series de potencias en el caso de EDOs de segundo orden
homogéneas siempre y cuando el punto entorno al que pensamos desarrollar la serie sea
un punto ordinario1 . Finalmente, en este capı́tulo veremos como podemos manejar y
calcular series de potencias utilizando Sage.
1
Dejaremos para el capı́tulo 8 y siguientes el caso de desarrollos de las soluciones entorno a los
llamados puntos singulares.
116 EDOs de segundo orden lineales

6.1. Series de potencias

Una serie de potencias es una expresión del tipo:
∞
∑ an(x − x0)n (6.2)
n=0

dados x0 y la sucesión de números a0 , a1 , a2 , . . . . Recordemos algunos resultados de

estas series2 :

1) La serie converge en el punto x si converge la sucesión3 de sumas parciales

bm = ∑m n
n=0 an (x − x0 ) , es decir si existe lı́mm→∞ bm . Cuando una serie no converge
en un punto se dice que diverge en ese punto.
∞
2) Una serie converge absolutamente en el punto x si ∑ |an(x − x0)|n converge
n=0
en ese punto.

3) ¿Cómo se sabe si una serie converge o no en un punto x? El criterio del cociente

puede ayudar. Definamos l(x) como

an (x − x0 )n an
l(x) = lı́m n−1
= |x − x0 | lı́m ,
n→∞ an−1 (x − x0 ) n→∞ an−1

entonces,

si l(x) > 1, la serie no converge;

si l(x) < 1, la serie converge absolutamente;
si l(x) = 1, el criterio no decide.

4) Si la serie ∑∞ n
n=0 an (x−x0 ) converge en x = x1 entonces converge absolutamente
∀x tal que |x − x0 | < |x1 − x0 |; si diverge para x = x1 entonces diverge ∀x tal que
|x − x0 | > |x1 − x0 |.

5) Definimos el radio de convergencia ρ = lı́mn→∞ aan−1 n

. Es tal que si |x − x0 | < ρ
la serie (6.2) converge absolutamente en x; si |x − x0 | > ρ la serie diverge en x. Si
ρ > 0, la serie (6.2) converge a una función
∞
f (x) = ∑ an(x − x0)n (6.3)
n=0

válida en el intervalo de convergencia x ∈ (x0 − ρ , x0 + ρ ).

2
Todos los enunciados siguientes son válidos también para números complejos.
3
Una sucesión bm converge al valor b si ∀ε > 0, ∃M tal que |bm − b| < ε, ∀m > M
6.1 Series de potencias 117

6a) Si la serie ∑∞ n
n=0 an (x − x0 ) converge a f (x) con radio de convergencia ρ1 y la
∞
serie ∑n=0 bn (x − x0 )n converge a g(x) con radio de convergencia ρ2 , entonces la
serie ∑∞ n
n=0 (an + bn )(x − x0 ) converge a f (x) + g(x) con un radio de convergencia
ρ ≥ mı́n(ρ1 , ρ2 ).
6b) En las condiciones anteriores, f (x)g(x) viene definido por la serie ∑∞
n=0 cn (x −
x0 ) n

⋆ Problema - 6.1. con cn = ∑nm=0 am bn−m .

El radio de convergencia de esta última serie es ρ ≥ mı́n(ρ1 , ρ2 ).

7) Si ρ > 0, la función definida por (6.3) es continua y tiene derivadas de cualquier

orden también continuas válidas dentro del intervalo de convergencia |x − x0 | < ρ ,
o x ∈ (x0 − ρ , x0 + ρ ) para números reales. Las derivadas se calculan como:

⋆ Problema - 6.2.
∞ ∞
f ′ (x) = ∑ nan(x − x0)n−1 = ∑ (n + 1)an+1(x − x0)n, (6.4)
n=0 n=0
∞ ∞
′′ n−2
f (x) = ∑ n(n − 1)an(x − x0) = ∑ (n + 2)(n + 1)an+2(x − x0)n, (6.5)
n=0 n=0
∞
(n + k)!
f (k) (x) = ∑ an+k (x − x0 )n . (6.6)
n=0 n!

Los coeficientes de la serie son:

⋆ Problema - 6.3.
f (n) (x0 )
an =
n!

8) Si ∑∞ n ∞ n
n=0 an (x − x0 ) = ∑n=0 bn (x − x0 ) , entonces an = bn , ∀n ≥ 0. En particular
∞ n
si ∑n=0 an (x − x0 ) = 0, entonces an = 0, ∀n ≥ 0.

9) Una función f (x) que tiene un desarrollo en serie (serie de Taylor) :

∞
f (x) = ∑ an(x − x0)n
n=0

válido para |x − x0 | < ρ , o x ∈ (x0 − ρ , x0 + ρ ) para números reales con ρ > 0 (esta
condición es esencial), se llama una función analı́tica en el punto x = x0 . Por
extensión, llamaremos “función analı́tica”(sin más) a la función que es analı́tica
(es decir, desarrollable en serie de Taylor con radio de convergencia no nulo) en
todos los puntos x0 de la recta real R.
118 EDOs de segundo orden lineales

Una√función que no es analı́tica no es una función “rara”. Por ejemplo, la función

f (x) = x no es analı́tica en x = 0pporque no tiene derivadas de ningún orden en el
punto x = 0. Tampoco lo es f (x) = |x| aunque es continua en x = 0. Otro caso clásico
es el de la función ( 2
e−1/x , si x 6= 0,
f (x) =
0, si x = 0.
que es continua y tiene todas las derivadas continuas en x = 0, pero al ser f (n) (x = 0) =
(n)
f (0) n
0, ∀n ≥ 0, la serie resultante es ∑∞
n=0 n! x = 0 que solo converge a la función f (x)
en x = 0 y tiene, por tanto, radio de convergencia nulo. Visualmente, ninguna de las
funciones anteriores tiene ningún aspecto “raro”.
10) Un punto donde una función f (x) no es analı́tica, se llama un punto singular
o singularidad de f (x).

11) Los polinomios son funciones analı́ticas para todo x.

6.2. Buscando la solución a las EDO de segundo or-

den homogéneas
Después de esta introducción, vayamos ahora al problema que nos ocupa. Se trata de
encontrar dos soluciones linealmente independientes, y1 (x), y2 (x), de la ecuación lineal
de segundo orden homogénea:

R(x)y′′ + P(x)y′ + Q(x)y = 0

o escrita en la forma
y′′ + p(x)y′ + q(x)y = 0 (6.7)
con p(x) = P(x)/R(x), q(x) = Q(x)/R(x). Recordemos que el teorema de existencia y
unicidad de la solución nos dice que si p(x) y q(x) son funciones continuas en un intervalo
(α , β ) entonces existe una solución única y(x) válida en todo ese intervalo y satisfaciendo
condiciones iniciales y(x0 ) = y0 , y′ (x0 ) = y′0 con x0 ∈ (α , β ). Esta solución se encuentra
escribiendo la solución general y(x) = c1 y1 (x) + c2 y2 (x) y eligiendo las constantes c1 , c2
de manera que se satisfagan las condiciones iniciales anteriores.
El método de solución para encontrar y1 (x), y2 (x) consiste en buscar soluciones en
forma de serie de potencias y(x) = ∑∞ n
n=0 an (x − x0 ) . Para que este desarrollo en serie
pueda conducir a encontrar y1 (x), y2 (x) es necesario imponer la condición de que p(x) y
q(x) sean analı́ticas en x0 , es decir que se pueda escribir

p(x) = ∑∞ n
n=0 pn (x − x0 ) , radio de convergencia ρ1 > 0,

q(x) = ∑∞ n
n=0 qn (x − x0 ) , radio de convergencia ρ2 > 0.

Si se da esta condición se dice que x0 es un punto ordinario de la ecuación diferencial.

Un punto que no es ordinario se llama un punto singular. Lo usual es que el punto
6.2 Buscando la solución a las EDO de segundo orden homogéneas 119

x0 coincida con el punto donde se fija la condición inicial, pero esto no es necesario y
podemos buscar las soluciones y1 (x), y2 (x) sin hacer, de momento, ninguna referencia a
la condición inicial.

⋆ Problema - 6.4. La ecuación de Bessel x2 y′′ + xy′ + (x2 − ν 2 ) = 0 tiene un punto de

singularidad en x = 0, todos los demás son puntos ordinarios. La ecuación de Legendre
(1 − x2 )y′′ − 2xy′ + α (α + 1)y = 0 tiene como puntos singulares x = −1, +1.

Teorema.- Si x0 es un punto ordinario de la ecuación diferencial (6.7), entonces es

posible encontrar dos soluciones y1 (x), y2 (x) de la forma y(x) = ∑∞ n
n=0 an (x − x0 ) con un
radio de convergencia ρ ≥ mı́n(ρ1 , ρ2 ).
El teorema tiene dos partes: (1) demostrar que es posible encontrar dos conjuntos de
coeficientes an de manera que en cada caso y(x) = ∑∞ n
n=0 an (x − x0 ) sea efectivamente
una solución de la ecuación diferencial y (2) demostrar que el radio de convergencia
de esta serie satisface ρ ≥ mı́n(ρ1 , ρ2 ). Demostraremos aquı́ solo la primera parte de
este teorema. La demostración es constructiva, es decir, vamos a encontrar expresiones
explı́citas para los coeficientes an . Lo mejor es empezar con un ejemplo. Sea la ecuación
y′′ +y = 0. Tomamos x0 = 0 (es lo más sencillo) y nos preguntamos si es posible encontrar
una solución de la forma y(x) = ∑∞ n
n=0 an x .

⋆ Problema - 6.5. Sustitución en la ecuación diferencial lleva a:

∞
∑ [(n + 2)(n + 1)an+2 + an] xn = 0
n=0

De donde deducimos que (n + 2)(n + 1)an+2 + an = 0, n = 0, 1, 2, . . . . De aquı́ se deduce

la relación de recurrencia:
−an
an+2 = , n = 0, 1, 2, . . .
(n + 2)(n + 1)

⋆ Problema - 6.6. De donde se obtiene

(−1)k
a2k = (2k)! a0 , k = 0, 1, 2, 3, . . .
(−1)k
a2k+1 = (2k+1)! a1 , k = 0, 1, 2, 3, . . .

O sea, que a0 y a1 no se determinan por el método, pero dados a0 y a1 se conocen

todos los coeficientes an , n ≥ 2. Veremos que esto está relacionado con las constantes
de integración de la solución general.

⋆ Problema - 6.7. Reordenando los términos de la serie, podemos escribir:

∞ ∞ ∞
(−1)k 2k (−1)k 2k+1
y(x) = ∑ anxn = a0 ∑ (2k)!
x + a1 ∑
(2k + 1)!
x
n=0 k=0 k=0
120 EDOs de segundo orden lineales

Lo bueno es que ya hemos acabado. Esta fórmula da la solución. Para entenderla un

poco más, definamos las siguientes funciones que llamaremos seno y coseno

(−1) 2k k
cos(x) ≡ ∑∞
k=0 (2k)! x ,

(−1) k
sin(x) ≡ ∑∞
k=0 (2k+1)! x
2k+1 .

De manera que la solución se escribe y(x) = a0 cos(x) + a1 sin(x), siendo a0 y a1 con-

stantes arbitrarias.

⋆ Problema - 6.8. Esta solución se obtiene fácilmente si vemos que y′′ + y = 0 es

una ecuación diferencial de coeficientes constantes probando y(x) = erx .

¿Pero son estas funciones las que hemos llamado seno y coseno de toda la vida?

⋆ Problema - 6.9. Las dos series que definen a sin(x) y cos(x) tienen radio de
convergencia infinito.

Esta definición de las funciones sin x y cos x es consistente con todas las propiedades
conocidas de dichas funciones definidas de otra manera.

⋆ Problema - 6.10. Derivando las series, se puede ver que:

d cos(x)
= − sin x,
dx
d sin(x)
= cos x.
dx
También se puede ver que cos(0) = 1, sin(0) = 0.

Quizá la propiedad más tı́pica sea cos(x)2 + sin(x)2 = 1.

⋆ Problema - 6.11. Lo podemos demostrar si verificamos que f (x) = cos(x)2 +

f (x)
sin(x)2 satisface d dx = 0, por lo que f (x) es constante igual a f (0) = 1.

⋆ Problema - 6.12. También podemos demostrar que sin(−x) = − sin(x), cos(−x) =

cos(x).

⋆ Problema - 6.13. Es muy útil la relación de Euler eiz = cos(z) + i sin(z) válida para
cualquier número complejo z.
6.3 Series de potencias con Sage 121

6.3. Series de potencias con Sage

Si deseamos calcular la serie de Taylor de una cierta función en Sage basta proceder
como en cualquiera de los siguientes ejemplos mediante la función taylor( f , xo , nterminos ),

sage: var(’x’); f(x)=sin(x);g(x)=cos(x)

sage: [Link](x,0,6) #6 primeros términos de sin(x) entorno a x=0
sage: taylor(g,2,4) #4 primeros términos de cos(x) entonro a x=2
Para calcular la función asociada a una serie de potencias podemos sumar sus térmi-
nos utilizando la función sum,
sage: var(’x k c n’)
sage: assume(abs(k*x)<1)
sage: sum(c*(k*x)ˆn, n,0,oo)
que dará como resultado −c/(k ∗ x − 1). En el anterior ejemplo conviene notar el uso
de la función assume() para fijar condiciones que se presupone deben obedecer nuestras
variables, en nuestro caso |kx| < 1. Por otra parte conviene también señalar que el infinito
en Sage se puede representar bien mediante dos letras o consecutivas oo, o mediante la
palabra in f ty o in f .
Como hemos visto en las secciones anteriores, con bastante frecuencia estaremos
interesados en calcular el radio de convergencia de una cierta función, supongamos por
ejemplo que queremos calcular el radio para la serie,
∞
xn
∑ .
n=0 (3n)!

Para ello introduciremos en Sage el siguiente código,

sage: var(’n’)
sage: a(n)=1/factorial(3*n)
sage: rho=limit(abs(a(n+1)/a(n)),n=+oo)
sage: rho
que nos dará como solución rho = 0. En este ejemplo conviene resaltar el uso de la
función f actorial(m) que nos calcula m!, ası́ como la función limit que podemos utilizar
siempre que queramos evaluar un lı́mite.
Otra forma muy potente para calcular series de potencias asociadas funciones con-
siste en usar los anillos matemáticos para series de potencias implementados en Sage.
Para nuestros propósitos, el siguiente ejemplo será suficiente: supongamos que queremos
calcular la serie de potencias asociada a f (t) = 1/(1 + t 2 ), en este caso escribiremos,
sage: R.<t> = PowerSeriesRing(QQ)
sage: f=1/(1+tˆ2) #No poner f(t), solo f
sage: [Link](5)
sage: f
122 EDOs de segundo orden lineales

El uso del anillo PoweSeriesRing fuerza la representación de f en serie de potencias. El

primer resultado mostrará los cinco primeros terminos de la serie t 4 −t 2 + 1 mientras que
el segundo resultado nos mostrará todos los terminos de la serie hasta la precisión que
tengamos por defecto establecida en nuestro Sage. Existen diversos anillos matemáticos
en Sage, para una mejor comprensión basta mencionar que ZZ corresponde a los números
racionales4 .
Como hemos visto, el proceso de buscar la solución a las EDO de segundo orden
mediante series de potencias nos conduce hacia una ecuación de recurrencia, como por
ejemplo (n+2)(n+1)an+2 +an = 0 para n = 0, 1, 2, ... Sage también nos permite obtener
fácilmente la solución a la ecuación de recurrencia como sigue

sage: maxima(’load(solve_rec)’)
sage: sol=maxima(’solve_rec((n+2)*(n+1)*a[n+2]+a[n]=0,a[n])’)
sage: show(sol)

la solución será,
k2 (−1)n/2 k1 (−1)3n/2
an = −
n! n!
donde k1 y k2 són constantes relacionadas con el valor que queremos que tomen los
primeros terminos. En el anterior ejemplo hemos usado maxima que es uno de los pa-
quetes especializados integrados en Sage. Si deseáramos fijar de antemano el valor de los
primeros términos, [Link]. a0 = 1 i a1 = 1, nada más sencillo, añadiremos dichas condiciones
al final de la instrucción solver ec, es decir: solver ec((n + 2) ∗ (n + 1) ∗ a[n + 2] + a[n] =
0, a[n], a[0] = 1, a[1] = 1).

4
Otros anillos implementados en Sage són por ejemplo ZZ para los enteros, RR para los reales, CC
para los números complejos
6.4 Ejercicios 123

6.4. Ejercicios
1. Para cuales de las ecuaciones diferenciales siguientes es posible encontrar una
solución en forma de serie de potencias ∑∞ n
n=0 an x ?

(a) y′′ + xex y′ + 2y = 0

(b) sin(x)y ′′ + y′ + y = 0
√
(c) y′′ + xy′ + cos(x)y = 0
y
(d) (1 − x)y′′ + xy′ + 1+x 2 =0

Determinar en cada caso qual serı́a el intervalo mı́nimo de validez de la solución

en serie.

2. Determinar el radio de convergencia y el intervalo de validez las siguientes series

de potencias.
∞
(−3)n
(a) ∑ (x − 5)n
n=1 n 7n+1
∞
(−1)n n
(b) ∑ (x + 3)n
n=0 4n
∞
2n
(c) ∑ (4x − 8)n
n=1 n
∞
x2n
(d) ∑ (−3)n
n=1
∞
(−1)n x2n
(e) ∑ (n!)2 22n
n=0

3. Demostar, haciendo uso del desarrollo en serie que ex+y = ex ey . Utilizar este resul-
tado y la relación d’Euler eix = cos(x) + i sin(x) para demostrar que: sin(x + y) =
sin(x) cos(y) + cos(x) sin(y).

4. Definiendo π como la primera raı́z positiva de sin(x) = 0 probar que sin(x + 2π ) =

sin(x). Calcular el número π haciendo uso de las aproximaciones polinómicas a
sin(x).

5. Probar que los coeficientes an tales que verifican la relación ∑∞

n=1 n an x
n−1 =

∑∞ n
n=0 an x vienen dados por an = a0 /n! donde n = 1, 2, 3, ....

6. Probar las siguientes afirmaciones:

(a) Si el radio de convergencia de una serie ∑ an xn es R, entonces el radio de
convergencia de la serie ∑ an x2n es R1/2 .
(b) Si una serie de potencias ∑ an xn converge para x = c, entonces tenemos
garantizado que convergirá absolútamente para todo valor x tal que |x| < |c|.

7. En la ecuación
y′′ + p(x)y′ + q(x)y = g(x)
124 EDOs de segundo orden lineales

hacemos el cambio de variable

 Zx 
1
z(x) = y(x) exp dt p(t)
2 x0
Demostar que la ecuación resultante no tiene termino en derivada primera de z(x).
A esta estructura se la llama la forma normal de la ecuación diferencial.
8. Resolver las siguientes ecuaciones hallando primero su forma normal:
(a) x2 y′′ + xy′ + 2− 1 y=0


 x 4
(b) y′′ − 2x y′ + a2 + x22 y = 0
(c) x2 y′′ + 4x2 y′ + (4x2 − 6)y = 0
¿Cuales de estas soluciones se pueden desarrollar en serie de Taylor entorno de
x = 0?
9. Encontrar una representación en serie de potencias para ln(1 + x) para |x| < 1
usando la relación existente entre ln(1 + x) y su derivada. Calcular su radio de
convergencia
10. El test de la integral de Cauchy o Maclaurin dice que si tenemos una serie ∑n f (n)
siendo fR(n) una función monótona decreciente entonces la serie convergirá si la
integral 1∞ f (s) ds es finita, en otro caso la serie será divergente.
(a) Probar dicho test de la integral.
(b) Si definimos la función zeta de Riemann como ξ (p) = ∑∞ −p
n=1 n , decir para que
valores de p dicha función tendrá un valor finito y para que valores divergirá hacia
el infinito.
11. Investigar la convergencia de las siguientes series:
∞
(a) ∑ (xk − xk+1).
k=0
∞
(−1)k x
(b) ∑ √k ( 2 − 1)k .
k=1

12. Encontrar una representación en serie de potencias para las siguientes funciones
entorno del punto xo = 0 y determinar el radio de convergencia de dichas series:
x2
(a) f (x) =
1 − 2x
1
(b) f (x) =
1 − x2
x2
(c) f (x) =
1+x
6x
(d) f (x) = 2
Z5xx − 4x − 1
ln(1 + t)
(e) f (x) = dt
0 t
(f) f (x) = sinh(x)