1.3.

ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 25

Suponiendo y ≥ 0 y C ≥ 0, se tiene
p
y0 = C 2 − y2,
ecuación de tipo separable, cuya solución es
arcsin(y/C) = x + D,
es decir
y = C sin(x + D),
que se convierte en
y(x) = C1 cos x + C2 sin x.

Ejemplo 1.48 (Caı́da libre) Recordamos la ecuación de la cáida libre (1.2): y” = −g, halleremos
su solució general mediante reducción de orden, utilizando las dos alternativas de reducción.
Comenzamos con la alterntiva A); para tal efecto, planteamos z(t) = y 0 (t), derivando y reemplazando
obtenemos
1
z 0 = −g (DR) ⇒ z = −gt + c, reemplazamos y 0 = −gt + c (DR) ⇒ y = − gt2 + ct + d,
2
Por lo tanto, la solución general es y = − 21 gt2 + ct + d .
Como en la ecuación diferencial no aparece t, la variable independiente de manera explı́cita, podemos
utilizar la alternativa B). Plantemos u(y) = y 0 (t), derivamos respecto a t, tenemos u0 y 0 = uu0 = y 00 ,
reemplazamos y obtenemos
1
uu0 = −g (S) ⇒ u2 = −gy + c, ⇒ u2 = −2gy + c,
2
Como u es la velocidad de caı́da libre, tenemos u ≤ 0, de donde
y0
u = − −2gy + c, reemplazamos y 0 = − −2gy + c (S) ⇒
p p
√ = 1,
− −2gy + c
integramos
Z
dy 1p
−2gy + c = t+d ⇒ −2gy + c = gt+d ⇒ −2gy+c = g 2 t2 +2gdt+d2 ,
p
√ = t+d ⇒
− −2gy + c g

1 (d2 − c) 1
de donde y = − gt2 − dt + = − gt2 + c1 t + c2 .
2 −2g 2

1.3.2. Ecuaciones diferenciales lineales de orden superior

Definición 1.49 Una ecuación diferencial lineal de orden n es una ecuación diferencial que
puede escribirse de la forma:

y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x)y = b(x) (L),

con a0 (x), . . . , an−1 (x) y b(x) funciones continuas que no dependen de la función incógnita y.
26 CAPÍTULO 1. ECUACIONES DIFERENCIALES ORDINARIAS

Ecuaciones diferenciales lineales homogéneas de orden superior

Se dirá que la ecuación (L) de orden n, es homogénea si b(x) es identicamente nula; es decir, si
la ecuación puede expresarse como
y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x)y = 0. (LH)
El estudio de las ecuaciones lineales homogéneas de orden n sigue el mismo camino que el seguido
para las ecuaciones lineales de primer orden.
Teorema 1.50 La solución general de una ecuación lineal homogénea de orden n, es un espacio
vectorial de dimensión n, para la adición de soluciónes y la multiplicación de las soluciónes por
constantes.

Demostración.- Al igual que en caso de primer orden, puede verificarse que la solución general de
(LH) es un subespacio vectorial. Probar que la dimensión es n, escapa de los objetivos del presente
curso.


Consecuencias Prácticas 1.51

1. Que la solución sea un espacio vectorial, significa que para toda familia de soluciones
{ϕ1 , ϕ2 , . . . , ϕm }, cualquier combinacioñ lineal es también solución; es decir:

c1 ϕ1 (x) + c2 ϕ2 (x) + · · · + ϕm (x) (1.26)

es solución de (LH)

2. Que la solución sea de dimensión n, significa que suficiente conocer una familia
{ϕ1 (x), . . . , ϕn (x)} de n soluciones linealmente independientes. Una tal familia se llama
Sistema fundamental de soluciones de (LH) (SF) y la solución general se escribe:

y = c1 ϕ1 (x) + c2 ϕ2 (x) + · · · + cn ϕn (x). (1.27)

Por lo tanto, nuestro propósito será determinar sistemas fundamentales de soluciones para expresar
la solución general de una ecuación diferencial (LH). La primera pregunta que nos hacemos es: ¿Si,
para una ecuación diferencial lineal homogénea de orden n, conocemos una familia de n solucio-
nes {ϕ1 (x), . . . , ϕn (x)}, cómo podemos saber si es o no es un sistema fundamental de soluciónes?
Razonemos:
{ϕ1 (x), . . . , ϕn (x)} es un (SF) ⇐⇒ {ϕ1 (x), . . . , ϕn (x)} es familia de soluciones y es linealmente
independiente, if f la única solución de la ecuación
c1 ϕ1 (x) + c2 ϕ2 (x) + · · · + cn ϕn (x) = 0,
con c1 , c2 ,...,cn como incógnitas es la trivial; es decir c1 = c2 = · · · = cn = 0, ⇐⇒ la única solución
del sistema lineal de n ecuaciónes y n incógnitas:

 c1 ϕ1 (x) + c2 ϕ2 (x) + · · · + cn ϕn (x) = 0,
 c1 ϕ01 (x) + c2 ϕ02 (x) + · · · + cn ϕ0n (x) = 0,


..

 .
(n−1) (n−1) (n−1)

(x) + · · · + cn ϕn

c1 ϕ1 (x) + c2 ϕ2 (x) = 0,
1.3. ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 27

es la trivial, ⇐⇒ la matriz asociada al sistema lineal de ecuaciones:

 ϕ1 (x) ϕ2 (x) ··· ϕn (x) 
 ϕ01 (x) ϕ02 (x) ··· ϕ0n (x) 
W (x) =  .. .. 
 . . 
(n−1) (n−1) (n−1)
ϕ1 (x) ϕ2 (x) · · · ϕn (x)

es inversible ⇐⇒ det W (x) 6= 0.

La matriz W (x) se llama matriz wronskiana asociada a la familia {ϕ1 (x), . . . , ϕn (x)} y
det W (x) se llama wronskiano asociado a la familia {ϕ1 (x), . . . , ϕn (x)}. En consecuencia, podemos
formular el teorema:

Teorema 1.52 Sea

y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x)y = 0. (LH)

y {ϕ1 (x), . . . , ϕn (x)} una familia de soluciones de (LH), entonces !son equivalentes los siguientes
enunciados:

(i) {ϕ1 (x), . . . , ϕn (x)} es un (SF) de (LH),

(ii) det W (x) 6= 0,

(iii) y = c1 ϕ1 (x) + c2 ϕ2 (x) + · · · + cn ϕn (x) es solución general de (LH).

Ejemplo 1.53 Consideremos la ecuación diferencial lineal homogénea de segundo orden

y 00 + y = 0,

cos x y sin x son soluciones de la ecuación en cuestión;

cos00 (x) + cos(x) = − cos(x) + cos(x) = 0,

sin00 (x) + sin(x) = − sin(x) + sin(x) = 0.

Ahora veamos si {cos x, sin x} forman un (SF), el wronskiano está dado por
 
cos x sin x
det W (x) = det = 1,
sinx cos x

de donde cos x y sin x es un sistema fundamental, por lo que la solución general está dada por
y(x) = c1 cos x + c2 sin x .

La siguiente pregunta que nos formulamos es la siguiente:¿Cómo podemos encontrar familias de

soluciones, mejor si son sistemas fundamentales de soluciones?
Lastimosamente, no existe un método general para encontrar sistemas fundamentales de solu-
ciones para una ecuación diferencial (LH); sino más bien, existen métodos para tipos especı́ficos de
ecuaciones lineales homogéneas de orden superior. En este texto veremos dos situaciones donde es
posible determinar sistemas fundamentales.
28 CAPÍTULO 1. ECUACIONES DIFERENCIALES ORDINARIAS

A) Ecuaciones diferenciales lineales de orden 2, donde se conoce una solución no nula

Sea
y 00 + p(x)y 0 + q(x)y = 0, LH
de orden 2 y ϕ(x) una solución no nula. Suponemos que la otra solución del sistema funda-
mental a determinar es de la forma
y(x) = c(x)ϕ(x),
donde c(x) es una función no constante, caso contrario tendrı́amos dos soluciones linealmente
dependientes. Derivando, se obtiene
y 0 (x) = c0 (x)ϕ(x) + c(x)ϕ0 (x),
y 00 (x) = c00 (x)ϕ(x) + 2c0 (x)ϕ0 (x) + c(x)ϕ00 (x);
reemplazando en la ecuación diferencial, se tiene
ϕ(x)c00 (x) + (2ϕ0 (x) + ϕ(x)p(x))c0 (x) + c(x)(ϕ00 (x) + p(x)ϕ0 (x) + q(x)ϕ(x)) = 0,
| {z }
=0

de donde c0 (x) es solución de la ecuación de primer orden lineal homogénea, que ya sabemos
resolver,
2ϕ0 (x) + ϕ(x)p(x)
z 0 (x) + z(x) = 0,
ϕ(x)
por consiguiente
c0 (x) = eA(x) ,
2ϕ0 (x) + ϕ(x)p(x)
donde A(x) es una primitiva de − , por lo tanto
ϕ(x)
Z x
c(x) = eA(s) ds.
x0

Ejemplo 1.54 Consideremos la ecuación de segundo orden

x2 y 00 − xy 0 − 3y = 0,
se verifica que y(x) = x3 es una solución particular no nula de la ecuación diferencial de se-
gundo orden en cuestion. Determinemos otra solución linealmente independiente. Planteamos
y(x) = c(x)x3 ,
derivando y remplazando en lae cuación diferencial, se obtiene
x5 c00 + (6x4 − x4 )c0 (x) = 0,
es decir c0 (x) es solución de la ecuación diferencial
5
z 0 = − z,
x
por lo que c0 (x) = 1/x5 , de donde
1
c(x) = − x−4 ,
4
1
La otra solución linealmente independiente encontrada es y = − x−4 · x, utilizando la parte 1
4
1
de las consecuencias prácticas 1.51, y = es solución La solución de la ecuación diferencial,
x
está dada por
c2
y(x) = c1 x3 + .
x
1.3. ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 29

B) Ecuaciones diferenciales lineales homogéneas a coeficientes constantes Una ecuación

diferencial lineal homogénea a coeficientes constantes de orden n es una ecuación diferencial
de la forma
y (n) + an−1 y (n−1) + · · · + a1 y 0 + a0 y = 0. (LHC)
Ahora intentemos encontrar sistemas fundamentales para este tipo de ecuaciones. Comenza-
mos recordando el caso de primer orden

y 0 + ay = 0 (LHC) ⇒ solución general y = ce−ax ⇒ SF = {e−ax }. (1.28)

Observamos que erx con r ∈ R

es una solución, (con r = −a). Busquemos soluciónes de la
rx
forma y = e para (LHC) de orden n. Consideramos la ecuació diferencial

y (n) + an−1 y (n−1) + · · · + a1 y 0 + a0 y = 0, (LHC)

derivamos y: y 0 = rerx , y 00 = r2 erx ,...,y (n) = rn erx , reemplazamos en (LHC), lo que da

rn erx + an−1 rn−1 erx + · · · + a1 rerx + a0 erx = (rn + an−1 rn−1 + · · · + a1 r + a0 )erx = 0,

llamamos p(λ) = λn + an−1 rλn−1 + · · · + a1 λ + a0 polinomio caracterı́stico de (LHC).

Hemos mostrado la:

Proposición 1.55 erx , con r ∈ R, es solución de la ecuación diferencial

y (n) + an−1 y (n−1) + · · · + a1 y 0 + a0 y = 0, (LHC)

si y solamente si r es raı́z del polinomio caracterı́stico de (LHC)

p(λ) = λn + an−1 rλn−1 + · · · + a1 λ + a0

Ejemplo 1.56 Consideremos la ecuación diferencial (LHC)

y 00 − 3y 0 + 2y = 0,

El polinomio caratecterı́stico de esta ecuación es

p(λ) = λ2 − 3λ + 2 = (λ − 2)(λ − 1). (1.29)

Las raı́ces del polinomio caracterı́stico son λ1 = 2, λ2 = 1. Por la proposición 1.55 e2x , ex
son soluciónes de la ecuación diferencial (1.29). Utilizando el wronskiano, deducimos que son
linealmente independientes y como el orden de la ecuación es 2, tenemos que la solución general
es y = c1 e2x + c2 ex .

Remarca 1.57 Raı́ces diferentes dan soluciones linealmente independientes; en efecto, la ma-
triz wronskiana corresponde casi a una matriz de Vandermonde.

Corolario 1.58 Sea

y (n) + an−1 y (n−1) + · · · + a1 y 0 + a0 y = 0, (LHC)
cuyo polinomio caracterı́stico tiene todas sus raı́ces reales y simples, que notamos r1 , r2 ,...,rn .
Entonces {er1 x , er2 x , . . . , ern x } es un sistema fundamental y

y = c1 er1 x + c2 er2 x + · · · cn ern x

es solución general.
30 CAPÍTULO 1. ECUACIONES DIFERENCIALES ORDINARIAS

La pregunta natural que uno se hace, después de haber visto la relación entre raı́ces del
polinomio caracterı́stico de una (LHC) y sus soluciónes, ¿qué sucede si las raı́ces del polino-
mio caracterı́stico son complejas no reales? y ¿qué sucede si existe alguna raı́z del polinomio
caracterı́stico que se repite?

Ejemplo 1.59

1. y 00 − 2y 0 + 2y = 0 tiene como polinomio caracterı́stico p(λ) = λ2 − 2λ + 2 cuyas raı́ces

son λ1 = 1 + i y λ2 = 1 − i raı́ces complejas no reales.
2. y 00 − 2y 0 + y = 0 tiene como polinomio caracterı́stico p(λ) = λ2 − 2λ1 + 1 = (λ − 1)2 .
λ = 1 es una raı́z que se repite dos veces, por la proposición 1.55, y = ex es una solución
no nula de la ecuación diferencial, falta una para tener un (SF).

Comencemos con la situación de raı́ces complejas.

Un poco de variable compleja
R
Recordemos que C el plano complejo es 2 provisto de una adición y una multiplicación que
hacen de C un cuerpo conmutativo. Un elemento z ∈ C puede escribirse de la manera siguiente

z = x + iy con x, y ∈ R.
Una función f : I ⊂ R → C admite la representación siguiente
f (x) = u(x) + iv(x),

donde u, v : I → . R
Ahora bien, si u y v son derivables, f es derivable y

f 0 (x) = u0 (x) + v 0 (x).

Proposición 1.60 Sean f, g : I ⊂ R → C, a ∈ C, entonces, se tiene:

i) (f + g)0 (x) = f 0 (x) + g 0 (x).
ii) (af )0 (x) = af 0 (x).
iii) (f (x)g(x))0 = f 0 (x)g(x) + f (x)g 0 (x)

Demostración.- Ejercicio.

Remarca 1.61 Las reglas de cálculo para derivadas son las mismas que para las funciones
reales.

Definición 1.62 Sea a = α + iβ ∈ C, x ∈ R, se define

eax = eαx (cos βx + i sin βx). (1.30)

Proposición 1.63 Se tiene:

i) Dados a, b ∈ C, se verifica
e(a+b)x = eax ebx . (1.31)
1.3. ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 31

ii) Para a ∈ C, se verifica

(eax )0 = aeax . (1.32)

Demostración.- Ejercicio.

Reconsideremos nuevamente la ecuación diferencial lineal homogénea de orden n a coeficientes

constantes
y (n) + an−1 y (n−1) + · · · + a1 y 0 + a0 y = 0. (LHC)
y planteamos
L(y) = y (n) + an−1 y (n−1) + · · · + a1 y 0 + a0 y. (1.33)
Remarcamos que L es una aplicación lineal y resolver la ecuación (LHC), es encontrar y(x)
tal que L(y) = 0.
Ahora bien, en lugar de considerar solamente soluciones reales, podemos prolongar nuestras
R
soluciones a soluciones y : → C.
Utilizando las reglas de cálculo para derivadas, se tiene para r ∈ C

L(erx ) = (erx )(n) + an−1 (erx )(n−1) + · · · + a1 (erx )0 + a0 erx

= rn erx + an−1 rn−1 erx + · · · + a1 rerx + a0 erx
= erx (rn + an−1 rn−1 + · · · + a1 r + a0 )
= erx p(r).

donde p(λ) es el polinomio caracterı́stico de (LHC) y 1.55 es válida para raı́ces complejas.
Además, si r = α + iβ, con β 6= 0, es una raı́z de p(λ), como los coeficientes de p(λ) son reales,
el conjugado de r, r̄ = α − iβ es también raiz de p(λ) y por consiguiente aportan al (SF):
{erx , er̄x }. Ahora bien, por la consecuencia 1 de 1.51, tenemos que:

1 αx+iβx 1 αx−iβx
eαx cos(βx) = e + e , (1.34)
2 2
1 αx+iβx 1
eαx sin(βx) = e − eαx−iβx (1.35)
2i 2i

son también soluciónes de (LHC) y como generan erx y er̄x , son linealmente independientes.

En sı́ntesis, las raı́ces complejas r = α + iβ y r̄ = α − iβ, con β 6= 0 aportan al (SF) de

(LHC) con eαx cos(βx) y eαx sin(βx).

Ejemplo 1.64 Reconsideremos la ecuación diferencial y 00 −2y 0 +2y = 0, tiene como polinomio
caracterı́stico p(λ) = λ2 − 2λ + 2 cuyas raı́ces son λ1 = 1 + i y λ2 = 1 − i raı́ces complejas
no reales. El aporte al sistema fundamental de este par de raı́ces es ex cos x y ex sin x. Como
el orden de la ecuación es 2, estás soluciones forman un sistema fundamental y la solución
general de la ecuación diferencial (LCH) es

y = c1 ex cos +c2 ex sin x.

32 CAPÍTULO 1. ECUACIONES DIFERENCIALES ORDINARIAS

Polinomio caracterı́stico con raı́ces múltiples

Cuando una raı́z del polinomio caracterı́stico, sea real o compleja, se repite más de una vez,
con los resultados vistos, faltan soluciones para tener completo un (SF). Se dirá que r es una
raı́z de multiplicidad m o que se repite m veces de p(λ), si podemos escribir p(λ):

p(λ) = (λ − r)m q(λ), con q(λ) polinomio y q(r) 6= 0. (1.36)

Ahora bien, planteando los operadores D(y) = y 0 e I(y) = y, el operador (1.33) asociado a la
ecuación diferencial (LHC) puede escribirse como

L(y) = (Dn + an−1 Dn−1 + · · · + a1 D + a0 I)(y), (1.37)

y si r es una raı́z de multiplicidad m del polinomio caracterı́stico p(λ) de (LHC), L(y) se

escribe también
L(y) = Q̃ ◦ (D − rI)m , (1.38)
donde Q(y) es el operador ecuación diferencial (LHC) asociado al polinomio caracterı́stico
q(λ). Observamos inmediatamente que si ϕ(x) es solución de (D − rI)m (y) = 0, entonces
también es solución de L(y) = O.
Por consiguiente, para ver el aporte de la raı́z r de multiplicidad m, es suficiente ver las
soluciones linealmente independientes, con las que aporta (D − rI)m (y) = 0. Tenemos

Proposición 1.65 Son soluciones linealmente independientes de (D − rj I)m (y) = 0 las fun-
ciones
erx , xerx , . . . , xm−1 erx . (1.39)

Demostración.- Por induccción sobre m. Para m = 1, el resultado ya sido mostrado. Supo-

nemos cierto para m − 1, a mostrar que es cierto para m. Solo debemos verificar que xm−1 erx
es solución, en efecto

(D − rI)m (xm−1 erx ) = (D − rI)m−1 (D − rI)(xm−1 erx ),

m−1 rx
= (D − rI)m−1 ((m − 1)xm−2 erx + rx| e {z− rxm−1 erx}),
=0
h.i.
= (m − 1)(D − rI)m−1 (xm−2 erx ) = 0

Remarca 1.66 Si r es una raiz no real de multiplicidad m, entonces r̄ es una raiz de mul-
tiplicidad m y contribuyen ambas raices con 2m soluciones linealmente al (SF), utilizando la
misma idea de (1.34) y (1.35); es decir

1 k αx+iβx
xk eαx cos(βx) = (x e + xk eαx−iβx ), (1.40)
2
1 k αx+iβx
xk eαx sin(βx) = (x e − xk eαx−iβx ) (1.41)
2i
para k = 0, 1, . . . , m − 1 y r = α + βi, con β 6= 0.

Finalmente enunciamos:
1.3. ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 33

Teorema 1.67 Consideremos la ecuación diferencial lineal homogénea a coeficientes

constantes
y (n) + an−1 y (n−1) + · · · + a1 y 0 + a0 y = 0. (LHC)
y sea p(λ) = λn +an−1 λn−1 +· · ·+a1 λ+a0 el respectivo polinomio caracterı́stico. Entonces:

a) Si r ∈ R es una raı́z que se repite m en p(λ), la contribución de r al (SF) es:

erx , xerx , . . . , xm−1 erx , (1.42)

un total de m soluciones.
b) Si α + iβ, α − iβ, con β 6= 0 es un par de raı́ces no reales que se repiten m veces, la
contribución del par de raı́ces al (SF) es:

eαx cos βx, xeαx cos βx, ..., xm−1 eαx cos βx
(1.43)
eαx sin βx, xeαx sin βx, ..., xm−1 eαx sin βx

Ejemplo 1.68 Consideremos la ecuacion diferencial (LHC)

y (6) − 2y (5) + 3y (4) − 4y (3) + 3y 00 − 2y 0 + y = 0,

el polinomio caracterı́stico de la ecuación está dado por

l(λ) = λ6 − 2λ5 + 3λ4 − 4λ3 + 3λ2 − 2λ + 1 = (λ − 1)2 (λ2 + 1)2 .

Deducimos que 1, i y −i son raices del polinomio caracterı́stico, las tres de multiplicidad 2,
por lo que el sistema fundamental buscado está dado por las soluciones:

ex , xex , cos x, sin x, x cos x, x sin x.

La solución general, será por consiguiente

y(x) = c1 ex + c2 xex + c3 cos x + c4 sin x + c5 x cos x + c6 x sin x.

C) Utilización de Series Enteras.- Recomendable para determinar la solución que satisfage un

problema a valor inicial. Recordemos los conceptos necesarios para desarrollar este procedi-
miento.

Definición 1.69 Sea f : I ⊂ R una función indefinidamente derivable, se dice que f es

análitica en el punto x0 , si
∞
X
f (x) = cn (x − x0 )n , |x − x0 | < ρ,
n=0

con 0 < ρ ≤ +∞.

Damos algunas de las propiedades, sin demostración de las funciones análiticas.

34 CAPÍTULO 1. ECUACIONES DIFERENCIALES ORDINARIAS

Proposición 1.70 Supongamos que f, g son análiticas en el punto x0 , entonces

∞
X
(αf + βg)(x) = (αcn + βdn )(x − x0 )n ,
n=0
X∞
f 0 (x) = ncn (x − x0 )n−1 ,
n=1
(n)
f (x0 ) = n!cn ,
∞ n
!
X X
f (x)g(x) = ck dn−k (x − x0 )n ,
n=0 k=0

donde
∞
X ∞
X
n
f (x) = cn (x − x0 ) , g(x) = dn (x − x0 )n .
n=0 n=0

Consideremos la ecuación lineal homogénea de orden n

y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x) = 0. (LH)

x0 es un punto ordinario, si cada una de las funciones ak (x) es análitica en x = x0 , caso

contrario se dirá que x0 es un punto singular. Enunciamos sin demostración el:

Teorema 1.71 Si x0 es un punto ordinario de la ecuación (LH), entonces cualquier solución

es análitica en x = x0 .

La determinación de un sistema fundamental para la ecuación (LH), se la hace eligiendo un

punto x0 que sea ordinario y se considera los n problemas a valor inicial

y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x) = 0,

y (k) (x0 ) = δkj , k = 1, . . . , n − 1; (j = 1, . . . , n),

donde δkj = 0 si k 6= j y δkj = 1 si j = k. Luego se resuelve cada uno de los problemas

planteando
∞
X
y(x) = cn (x − x0 )n ,
n=0

desarrollando en series de potencias cada una de las aj (x) y utilizando las reglas de cálculo
dadas en la proposición 1.70.

Ejemplo 1.72 Determinemos la solución general de

y 00 + xy = 0.

Observamos que la ecuación es (LH) y los coeficientes de las derivadas de y son análiticas,
(a1 (x) = 0, a0 (x) = x), para cualquier punto de la recta real. Elegimos por comodidad x0 = 0
ya que x ya está desarrollada en serie de potencias. Planteando
∞
X
y(x) = cn xn ,
n=0
1.3. ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 35

se tiene Y cn = y (n) (0)/n! y en particular y 00 (0) = 0y(0) da c2 = 20.

Por otro lado para n ≥ 3, se tiene
n−2
X n−2

y n (x) = (−xy(x))(n−2) = − x(k) y (n−k−2) (x) = −xy (n−2) (x) − (n − 2)y (n−3) (x),
k
k=0

de donde
(n − 2)(n − 3)!
cn = − cn−3 . (1.44)
n!
Se observa inmediatamente que c3k+2 = 0.
Planteando c0 = y0 y c1 = y1 , donde y0 y y1 son los valores iniciales y utilizando la fórmula
(1.44), se obtiene los valores de los restantes coeficientes.

Puntos Singulares Regulares

Definición 1.73 Se dice que el punto x0 singular de la ecuación (LH) de orden n es un punto
singular regular si las funciones

an−k (x − x0 )k , k = 1, . . . , n

son funciones anaáliticas.

Lo ideal es trabajar en puntos ordinarios; sin embargo la descripción de algunos fenómenos

pasa por la solución de problemas diferenciales a valor inicial, en el que x0 es un punto singular
regular.
Para resolver utilizando series enteras, se plantea
∞
X ∞
X
γ n
y(x) = (x − x0 ) cn (x − x0 ) = cn (x − x0 )n+γ
n=0 n=0

R
donde γ ∈ , (se deberá determinar su valor) y se desarrolla en serie de potencias las funciones
an−k (x)(x − x0 )k para obtener una ecuación para γ y relaciones recursivas para los ck .

Ejemplo 1.74 Resolvamos la ecuación diferencial

x2 y 00 + xy 0 − m2 y = 0.

Para que el problema tenga sentido, cualquier problema a valor inicial en x = 0, debe satisfacer
la condición que y0 = 0, dejando libre y1 . Rescribiendo la ecuación en la forma requerida, se
tiene
1 m2
y 00 + y 0 − 2 y = 0.
x x
Se observa inmediatamente que 1 y m2 constantes son análiticas.
Planteando
X∞
y(x) = cn xn+γ ,
n=0

se obtiene
y 0 (x) = P∞ n+γ−1 ,
P
n=0 (n + γ)cn x
00 ∞
y (x) = n=0 (n + γ)(n + γ − 1)cn xn+γ−2 .
36 CAPÍTULO 1. ECUACIONES DIFERENCIALES ORDINARIAS

Remplazando en la ecuación diferencial, se tiene

P∞
− 1) + (n + γ) − m2
cn xn+γ−2 = 0,


n=0 (n + γ)(n + γP
∞ 2 2 n+γ−2 = 0
n=0 (n + γ) − m cn x

de donde obtenemos las ecuaciones

(n + γ)2 − m2 cn = 0.

De estas ecuaciones deducimos que cn = 0 para todo n excepto para solamente uno de los cn .
Si suponemos c0 6= 0, se tiene
γ 2 − m2 = 0,
por lo que γ = m o γ = −m. De aquı́ obtenemos dos solucioens diferentes:

y(x) = xm , y(x) = x−m .

Ecuaciones diferenciales lineales (no homogéneas)

Teorema 1.75 La solución general de la ecuación diferencial lineal (L) de orden n,

y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x)y = b(x) (L)

está dada por la siguiente regla memotécnica:

Solución general de (LH)

Solución general de (L) = asociada: + Una solución (parti- (1.45)
(n) (n−1)
cular) de (L)
y +an−1 (x)y +···+a0 (x)y=0

Demostración.- Sea ψ una solución particular de (L) y ϕ cualquier solución de (LH), entonces

ψ+ϕ

es una solución de (L); en efecto

(ϕ + ψ)0 (x) = ϕ0 (x) + ψ 0 (x)

= a(x)ϕ(x) + (a(x)ψ(x) + b(x))
= a(x)(ϕ + ψ)(x) + b(x).

Sea ψ una solución particular de (L), entonces cualquier solución de (L) es de la forma

ψ + ϕ,

donde ϕ es solución de (LH). En efecto, sea φ una solución de (L), suficiente ver que ψ − φ es
solución de (LH):

(ψ − φ)0 (x) = ψ 0 (x) − φ0 (x)

= (a(x)ψ(x) + b(x)) − (a(x)φ(x) + b(x))
= a(x)(ψ(x) − φ(x)).

1.3. ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 37

Determinación de una solución particular de (L)

La determinación de la solución general de una ecuación lineal de primer orden pasa: por la
determinación de la solución general de la ecuación lineal homogenea asociada (algo que ya sabemos
hacer) y por la determinación de una solución particular de (L). Al igual que en el caso de primer
orden, para el caso de orden superior veremos las dos mismas maneras de encontrar una solución
particular:
A) Por tanteo o adivinando una solución particular
Este procedimiento es producto de la experiencia y la genialidad de la persona que lo utiliza;
sin embargo, en algunos casos se puede encontrar la solución, sobre todo, cuando b(x) la parte
no homogénea de (L) es un polinomio, una función sin o cos y cuando es exponencial y la
forma de proceder es análoga al caso de primer orden.

B) Variación de constantes
Cuando no es posible adivinar las soluciones particulares, se tiene el método de variación de
constantes.
Sea {ϕ1 , ϕ2 , . . . , ϕn } un sistema fundamental de soluciones de (LH). Al igual que en el caso
de primer orden se supone que la solución particular es de la forma
n
X
ψ(x) = ck (x)ϕk (x).
k=1

Derivando ψ, se obtiene
n
X n
X
ψ 0 (x) = c0k (x)ϕk (x) + ck (x)ϕ0k (x),
k=1 k=1

suponiendo
n
X
ck (x)ϕ0k (x) = 0,
k=1
se tiene
n
X
0
ψ (x) = c0k (x)ϕk (x).
k=1

Derivemos ψ0 y obtenemos
n
X n
X
00
ψ (x) = c0k (x)ϕ0k (x) + ck (x)ϕ00k (x),
k=1 k=1

suponiendo
n
X
c0k (x)ϕ0k (x) = 0,
k=1
se tiene
n
X
ψ 00 (x) = ck (x)ϕ0k (x).
k=1

Repetimos el proceso de derivación hasta determinar ψ (n−1) , suponiendo por lo tanto

n
(j)
X
c0k (x)ϕk (x) = 0, j = 0, 1, . . . , n − 2.
k=1
38 CAPÍTULO 1. ECUACIONES DIFERENCIALES ORDINARIAS

y obteniendo por consiguiente

n
(j)
X
(j)
ψ (x) = ck (x)ϕk (x), j = 0, 1, . . . , n − 1.
k=1

Finalmente introducimos los valores de ψ (j) en la respectiva ecuación diferencial, obteniendo

n n n−1 n
!
(n−1) (n) (j)
X X X X
0
ck (x)ϕk (x) + ck (x)ϕk (x) + aj (x) ck (x)ϕk (x) = b(x).
k=1 k=1 j=0 k=1

Reagrupando los terminos, se obtiene

 
n n n−1
(n−1) (n) (j)
X X X
c0k (x)ϕk (x) + ck (x) ϕk (x) + aj (x)ϕk (x) = b(x)
k=1 k=1 j=0
| {z }
=0 solución de (LH)

Resumiendo hemos obtenido n ecuaciones lineales para c01 , c02 , . . . , c0n ,

c01 (x)ϕ1 (x) + c02 (x)ϕ2 (x) + ··· + c0n (x)ϕn (x) = 0
c01 (x)ϕ01 (x) + c02 (x)ϕ02 (x) + ··· + c0n (x)ϕ0n (x) = 0
.. ..
. .
(n−1) (n−1) (n−1)
c01 (x)ϕ1 (x) + c02 (x)ϕ2 (x) + ··· + c0n (x)ϕn (x) = b(x)

que escrito de manera matricial

ϕ1 (x) ϕ2 (x) ··· ϕn (x)   c0 (x) 

0
  
1
 ϕ01 (x) ϕ02 (x) ··· ϕ0n (x)   c02 (x) 
 0 
.. ..   ..   ... ;
. =
    
 . .
(n−1) (n−1) (n−1)
ϕ1 (x) ϕ2 (x) · · · ϕn (x) c0 (x) b(x)
| {z } | n{z } | {z }
W (x) C 0 (x) B(x)

es decir
W (x)C 0 (x) = B(x).
Ahora bien, como {ϕ1 , ϕ2 , . . . , ϕn } es un sistema fundamental de (LH), W (x) es inversible,
por lo que
C 0 (x) = (W (x))−1 B(x),
obtenemos las funciones ck (x) integrando las respectivas c0k (x)

Ejemplo 1.76 Hallemos la solución general de

1
y 00 + y = .
sin x
La ecuación (LH) asociada, admite como sistema fundamental {cos x, sin x}. La aplicación de
variación de constantes conduce a considerar el sistema de ecuaciones
  0   
cos x sin x c1 (x) 0
= .
− sin x cos x c02 (x) 1/over sin x
1.3. ECUACIONES DIFENCIALES DE ORDEN SUPERIOR 39

Resolviendo el sistema, por ejemplo por determinantes obtenemos:

cos x
c01 (x) = −1, c02 (x) = ,
sin x
de donde
c1 (x) = −x, c2 (x) = ln(sin x).
La solución general está dada por consiguiente

y(x) = c1 cos x + c2 sin x − x cos x + ln(sin x) sin x.

Proposición 1.77 El problema a valor inicial

y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x)y = b(x)




 y(x0 ) = y0 ,



y 0 (x0 ) = y1 ,
 ..



 .
y (n−1) (x0 ) = yn−1


tiene solución única.

Demostración.- En clases.

Ejemplo 1.78 Hallar el valor de y(ln 2), sabiendo que y es la solución del problema a valor inicial:
 00
 y + y 0 − 2y = 4e2t ,
y(0) = 6,
 0
y (0) = −5.

Hallamos la solución general de la ecuación (L) del problema a valor inicial; para tal efecto, comen-
zamos resolviendo la ecuación (LH) asociada:

y 00 + y 0 − 2y = 0, (LHC) ⇒ p(λ) = λ2 + λ − 2 = (λ + 2)(λ − 1),

Las raı́ces del polinomio caracterı́stico son λ = −2 y λ = 1, de donde SF = {e−2t , et } y la solución

general de (LH) asociada es y = c1 e−2t + c2 et .
Pasamos a encontrar una solución particular de (L) por tanteo, planteando y = αe2t , derivamos y
reemplazamos, lo que da:
4αe2t + 2αe2t − 2αe2t = 4e2t ⇒ α = 1,
solución particular encontrada y = e2t y solución general de (L):

y = c1 e−2t + c2 et + e2t .

Ahora, encontramos los valores de c1 y c2 reemplazando las condiciones iniciales en la solución

general:  
y(0) = c1 + c2 + 1 = 6, c1 + c2 = 5,
⇒ ⇒ c1 = 4, c2 = 1.
−2c1 + c2 + 2 = −5. −2c1 + c2 = −7

La solución del problema a valor inicial es y = 4e−2t + et + e2t e y(ln 2) = 1 + 2 + 4 = 7.