ECUACIONES DIFERENCIALES

LORENA MARITZA TERRIOS GUZMÁN

CAPÍTULO 4: ECUACIONES
DIFERENCIALES DE ORDEN SUPERIOR

Pág.

4.1 Teoría básica de las ecuaciones diferenciales homogéneas

de orden superior 46

4.2 Solución de ecuaciones diferenciales homogéneas con

coeficientes constantes 49

4.3 Teoría básica de las ecuaciones diferenciales no homogéneas 52

4.4 Métodos de solución para ecuaciones diferenciales no

homogéneas con coeficientes constantes

1. Método de los coeficientes indeterminados 53

2. Método del anulador 56

3. Método de variación de parámetros 61

4.5 Ecuaciones diferenciales homogéneas y no homogéneas con

coeficientes variables: Ecuación de Cauchy – Euler 63

Taller de ejercicios 67
 Lorena Maritza Terrios Guzmán 46
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

Capítulo 4 ECUACIONES DIFERENCIALES HOMOGÉNEAS DE ORDEN

SUPERIOR

Este capítulo inicia con la teoría básica de las ecuaciones diferenciales de

orden superior: definición de ecuación diferencial de orden superior,
clasificación y resultados básicos de la solución de una ecuación diferencial
de orden superior, además se presenta la definición de conjunto
fundamental de soluciones y de Wroskiano así como resultados que
relacionan estos dos conceptos. Posteriormente se describen métodos que
permiten hallar la solución general de ecuaciones diferenciales homogéneas
y no homogéneas con coeficientes constantes. El capítulo concluye con la
definición y método de solución de la ecuación de Cauchy-Euler que es, una
ecuación diferencial con coeficientes variables.

Para el desarrollo de la teoría y de los métodos que se abordan en este

capítulo se hace necesario el conocimiento y manejo de determinante de
una matriz, resultados que permiten obtener el determinante de una
matriz, propiedades respecto a dependencia e independencia lineal,
derivadas de orden superior, solución de ecuaciones polinomiales y solución
de sistemas de ecuaciones.

4.1 TEORÍA BÁSICA DE LAS ECUACIONES DIFERENCIALES

HOMOGÉNEAS DE ORDEN SUPERIOR

DEFINICIÓN 4.1: ECUACIONES DIFERENCIALES HOMOGÉNEAS. Una ecuación diferencial

ordinaria, lineal de orden n de la forma

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥 ) 𝑛
+ 𝑎𝑛−1 ( 𝑥 ) 𝑛−1
+ ⋯ + 𝑎1 (𝑥 ) + 𝑎0 (𝑥 )𝑦 = 0 (1)
𝑑𝑥 𝑑𝑥 𝑑𝑥

se llama ecuación diferencial homogénea.

NOTA:
 En este contexto la palabra homogénea no indica que los coeficientes sean funciones
homogéneas.

 Otra forma de denotar la ecuación diferencial de la definición 4.1 es

 Lorena Maritza Terrios Guzmán 47
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

𝐹(𝑥, 𝑦, 𝑦´, 𝑦 ′′, … , 𝑦 (𝑛−1) , 𝑦 (𝑛) ) = 0

TEOREMA 4.1: PRINCIPIO DE SUPERPOSICIÓN. Si 𝑦1 , 𝑦2 , … , 𝑦𝑘 son soluciones de una

ecuación diferencial lineal homogénea de orden n
𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥 ) 𝑛 + 𝑎𝑛−1 (𝑥 ) 𝑛−1 + ⋯ + 𝑎1 (𝑥 ) + 𝑎0 (𝑥 )𝑦 = 0 (1)
𝑑𝑥 𝑑𝑥 𝑑𝑥
entonces la combinación lineal 𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + ⋯ + 𝑐𝑘 𝑦𝑘 también es solución, donde cada
𝑐𝑖 , 𝑖 = 1,2, … , 𝑘, son constantes cualesquiera.

Demostración:
Como 𝑦1 , 𝑦2 , … , 𝑦𝑘 son soluciones de la E.D.L.H. (1), entonces

𝑎𝑛 (𝑥 )𝑦1 (𝑛) + 𝑎𝑛−1 (𝑥 )𝑦1(𝑛−1) + ⋯ + 𝑎1 (𝑥 )(𝑦1 )′ + 𝑎0 (𝑥 )𝑦1 = 0 (2)

𝑎𝑛 (𝑥 )𝑦2 (𝑛) + 𝑎𝑛−1 (𝑥 )𝑦2 (𝑛−1) + ⋯ + 𝑎1 (𝑥 )(𝑦2 )′ + 𝑎0 (𝑥 )𝑦2 = 0 (3)
⋮
𝑎𝑛 (𝑥 )𝑦𝑘 (𝑛) + 𝑎𝑛−1 (𝑥 )𝑦𝑘 (𝑛−1) + ⋯ + 𝑎1 (𝑥 )(𝑦𝑘 )′ + 𝑎0 (𝑥 )𝑦𝑘 = 0 (𝑘)

Ahora,

𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + ⋯ + 𝑐𝑘 𝑦𝑘 ⟺ 𝑎0 (𝑥 )𝑦 = 𝑎0 (𝑥 )[𝑐1 𝑦1 + 𝑐2 𝑦2 + ⋯ + 𝑐𝑘 𝑦𝑘 ]

𝑦 ′ = 𝑐1 (𝑦1 )′ + 𝑐2(𝑦2 )′ + ⋯ + 𝑐𝑘 (𝑦𝑘 )′ ⟺ 𝑎1 (𝑥 )𝑦 = 𝑎1 (𝑥 )[𝑐1 (𝑦1)′ + 𝑐2 (𝑦2)′ + ⋯ + 𝑐𝑘 (𝑦𝑘 )′ ]

⋮
𝑦 = 𝑐1(𝑦1 ) + 𝑐2(𝑦2 ) + ⋯ + 𝑐𝑘 (𝑦𝑘 ) ⟺ 𝑎𝑛 (𝑥 )𝑦 (𝑛) = 𝑎𝑛 (𝑥 )[𝑐1(𝑦1 )(𝑛) + ⋯ + 𝑐𝑘 (𝑦𝑘 )(𝑛)]
(𝑛) (𝑛) (𝑛) (𝑛)

Sea
𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑀 = 𝑎𝑛 (𝑥 ) 𝑛
+ 𝑎𝑛−1 ( 𝑥 ) 𝑛−1
+ ⋯ + 𝑎1 (𝑥 ) + 𝑎0 (𝑥 )𝑦
𝑑𝑥 𝑑𝑥 𝑑𝑥
Entonces

𝑀 = 𝑎𝑛 (𝑥)[𝑐1 (𝑦1 )(𝑛) + 𝑐2 (𝑦2 )(𝑛) + ⋯ + 𝑐𝑘 (𝑦𝑘 )(𝑛) ] + ⋯ + 𝑎1 (𝑥)[𝑐1 (𝑦1 )′ + 𝑐2 (𝑦2 )′ + ⋯ + 𝑐𝑘 (𝑦𝑘 )′ ]

Aplicando propiedades obtenemos:

𝑀 = 𝑐1 [𝑎𝑛 (𝑥)𝑦1 (𝑛) + ⋯ + 𝑎1 (𝑥)(𝑦1 )′ + 𝑎0 (𝑥)𝑦1 ] + ⋯ + 𝑐𝑘 [𝑎𝑛 (𝑥)𝑦𝑘 (𝑛) + ⋯ + 𝑎1 (𝑥)(𝑦𝑘 )′ + 𝑎0 (𝑥)𝑦𝑘 ]

Sustituyendo (1), (2),…, (k) en la última igualdad se tiene que

𝑀 = 𝑐1 (0) + ⋯ + 𝑐𝑘 (0) = 0

Por lo tanto, 𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + ⋯ + 𝑐𝑘 𝑦𝑘 es solución de la ecuación diferencial.

DEFINICIÓN 4.2: DEPENDENCIA E INDEPENDENCIA LINEAL DE FUNCIONES. Se dice que un

conjunto de funciones 𝑓1 (𝑥), 𝑓2 (𝑥), … , 𝑓𝑛 (𝑥) es linealmente dependiente en un intervalo I si
existen constantes 𝑐1 , 𝑐2 , … , 𝑐𝑛 no todas iguales a cero, tales que
𝑐1 𝑓1(𝑥 ) + 𝑐2𝑓2 (𝑥 ) + ⋯ + 𝑐𝑛 𝑓𝑛 (𝑥 ) = 0
 Lorena Maritza Terrios Guzmán 48
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

para toda 𝑥 en I. Si 𝑐1 = 𝑐2 = ⋯ = 𝑐𝑛 = 0 , entonces se dice que el conjunto de funciones es

linealmente independiente.
NOTA:
Para establecer si un conjunto es linealmente dependiente o independiente se pueden utilizar
además de la definición 4.2 propiedades que den a lugar, por ejemplo: si un conjunto de
funciones (establecido un orden) es linealmente dependiente, entonces existe al menos una
función que se puede expresar como combinación lineal de las que la preceden.

EJEMPLO 1
Demostremos que el conjunto de funciones {𝑠𝑒𝑛2𝑥, 𝑠𝑒𝑛𝑥𝑐𝑜𝑠𝑥 } es linealmente dependiente.

Solución:
Sean 𝑓1 (𝑥 ) = 𝑠𝑒𝑛2𝑥, 𝑓2 (𝑥 ) = 𝑠𝑒𝑛𝑥𝑐𝑜𝑠𝑥

Aplicando identidades trigonométricas tenemos que:

1 1
𝑓2 (𝑥 ) = 𝑠𝑒𝑛𝑥𝑐𝑜𝑠𝑥 = 𝑠𝑒𝑛2𝑥 = 𝑓1 (𝑥 )
2 2

Es decir, 𝑓2 se define como combinación lineal de 𝑓1 , por lo tanto 𝑓1 y 𝑓2 forman un

conjunto linealmente dependiente.

EJERCICIO 1
Determine si el conjunto {𝑥 2 + 1, 𝑥 } es linealmente dependiente o linealmente
independiente.

DEFINICIÓN 4.3: WRONSKIANO. Sean 𝑓1 (𝑥), 𝑓2 (𝑥), … , 𝑓𝑛 (𝑥) funciones que son al menos
(𝑛 − 1) veces derivables. El Wronskiano de 𝑓1 (𝑥 ), 𝑓2 (𝑥 ), … , 𝑓𝑛 (𝑥 ) se define como

𝑓1 (𝑥) 𝑓2 (𝑥) ⋯ 𝑓𝑛 (𝑥)

𝑊 (𝑥 ) = 𝑊 (𝑓1 , 𝑓2 , … , 𝑓𝑛 ) = || 𝑓1 ′(𝑥) 𝑓2 ′(𝑥) ⋯ 𝑓𝑛 ′(𝑥) |
|
⋮ ⋮ ⋮ ⋮
(𝑛−1) (𝑛−1)
𝑓1 (𝑥) 𝑓2 (𝑥) ⋯ 𝑓𝑛 (𝑛−1) (𝑥)

TEOREMA 4.2: CRITERIO PARA SOLUCIONES LINEALMENTE INDEPENDIENTES. Sean

𝑦1 , 𝑦2 , … , 𝑦𝑛 soluciones de una ecuación diferencial lineal homogénea de orden n, en un
intervalo I. El conjunto de soluciones es linealmente independiente en I si y solo si
𝑊 (𝑦1 , 𝑦2 , … , 𝑦𝑛 ) ≠ 0, para toda 𝑥 en I.

EJEMPLO 2
Utilicemos el teorema 4.2 para analizar el conjunto de funciones {𝑒 𝑥 , 𝑒 2𝑥 , 𝑒 3𝑥 } .

Solución:
𝑒𝑥 𝑒 2𝑥 𝑒 3𝑥
𝑊 (𝑒 𝑥 , 𝑒 2𝑥 , 𝑒 3𝑥 ) = |𝑒 𝑥 2𝑒 2𝑥 3𝑒 3𝑥 |
𝑒𝑥 4𝑒 2𝑥 9𝑒 3𝑥
Resolvamos el determinante por la primera fila.
 Lorena Maritza Terrios Guzmán 49
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

2𝑥
𝑊 (𝑒 𝑥 , 𝑒 2𝑥 , 𝑒 3𝑥 ) = 𝑒 𝑥 |2𝑒 2𝑥 3𝑒 3𝑥 | − 𝑒 2𝑥 |𝑒 𝑥 3𝑒 3𝑥 | + 𝑒 3𝑥 |𝑒 𝑥 2𝑒 2𝑥 |
4𝑒 9𝑒 3𝑥 𝑒𝑥 9𝑒 3𝑥 𝑒𝑥 4𝑒 2𝑥
𝑊 (𝑒 𝑥 , 𝑒 2𝑥 , 𝑒 3𝑥 ) = 𝑒 𝑥 (18𝑒 5𝑥 − 12𝑒 5𝑥 ) − 𝑒 2𝑥 (9𝑒 4𝑥 − 3𝑒 4𝑥 ) + 𝑒 3𝑥 (4𝑒 3𝑥 − 2𝑒 3𝑥 )

𝑊 (𝑒 𝑥 , 𝑒 2𝑥 , 𝑒 3𝑥 ) = 6𝑒 6𝑥 − 6𝑒 6𝑥 + 2𝑒 6𝑥 2𝑒 6𝑥 ≠ 0, ∀𝑥 ∈ ℝ

Por el teorema 4.2 se concluye que {𝑒 𝑥 , 𝑒 2𝑥 , 𝑒 3𝑥 } es linealmente independiente.

EJERCICIO 2
Analice la independencia lineal de cada conjunto de funciones utilizando el teorema 4.2.
2
a) {𝑥𝑒 𝑥 , 𝑒 𝑥 } b) {𝑥 2 , 2𝑥, 2𝑥 2 + 4𝑥 }

DEFINICIÓN 4.4: CONJUNTO FUNDAMENTAL DE SOLUCIONES. Cualquier conjunto

𝑦1 , 𝑦2 , … , 𝑦𝑛 de soluciones linealmente independientes de una ecuación diferencial lineal
homogénea de orden n en un intervalo I se denomina conjunto fundamental de soluciones en
el intervalo.

TEOREMA 4.3: EXISTENCIA DE UN CONJUNTO FUNDAMENTAL. Existe un conjunto

fundamental de soluciones para la ecuación diferencial lineal homogénea de orden n en un
intervalo I.

TEOREMA 4.4: SOLUCIÓN GENERAL. Si 𝑦1 , 𝑦2 , … , 𝑦𝑛 es un conjunto fundamental de

soluciones de la ecuación diferencial lineal homogénea de orden n en un intervalo I, entonces
la solución general en el intervalo es
𝑦 = 𝑐1𝑦1 (𝑥 ) + 𝑐2 𝑦2 (𝑥 ) + ⋯ + 𝑐𝑛 𝑦𝑛 (𝑥 )
donde 𝑐1 , 𝑐2, … , 𝑐𝑛 son constantes arbitrarias.

TAREA 1: Determine si el conjunto de funciones dado es linealmente dependiente o

linealmente independiente.

a) {𝑥, 𝑥𝑙𝑛𝑥 } b) {𝑒 𝑥 , 𝑒 𝑥+2 , 𝑒 −𝑥 } c) {1, 𝑠𝑒ℎ𝑛𝑥, 𝑐𝑜𝑠ℎ𝑥 }

4.2 SOLUCIÓN DE ECUACIONES DIFERENCIALES HOMOGÉNEAS

DE ORDEN SUPERIOR CON COEFICIENTES CONSTANTES

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

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 𝑛 + 𝑎𝑛−1 𝑛−1 + ⋯ + 𝑎1 + 𝑎0 𝑦 = 0 (3)
𝑑𝑥 𝑑𝑥 𝑑𝑥
donde 𝑎0 , 𝑎1 , … , 𝑎𝑛−1 , 𝑎𝑛 son números reales con 𝑎𝑛 ≠ 0
 Lorena Maritza Terrios Guzmán 50
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

NOTA:
 Como los coeficientes 𝑎0 , 𝑎1 , … , 𝑎𝑛−1 , 𝑎𝑛 , son funciones constantes, entonces son continuas
para todo número real y 𝑎𝑛 ≠ 0, entonces la ecuación diferencial (3) tiene soluciones
definidas para número real.
𝑑𝑦 𝑑2 𝑦 𝑑𝑛 𝑦
 Si 𝑦 = 𝑒 𝑚𝑥 , entonces = 𝑚𝑒 𝑚𝑥 , 2
= 𝑚2 𝑒 𝑚𝑥 , … , = 𝑚𝑛 𝑒 𝑚𝑥 .
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑛

Luego en (3):
𝑎𝑛 𝑚𝑛 𝑒 𝑚𝑥 + 𝑎𝑛−1 𝑚𝑛−1 𝑒 𝑚𝑥 + ⋯ + 𝑎1 𝑚𝑒 𝑚𝑥 + 𝑎0 𝑒 𝑚𝑥 = 𝑒 𝑚𝑥 (𝑎𝑛 𝑚𝑛 + 𝑎𝑛−1 𝑚𝑛−1 + ⋯ + 𝑎1 𝑚 + 𝑎0 )

= 𝑒 𝑚𝑥 𝑃(𝑚) = 0 (4)
𝑛 𝑛−1
donde 𝑃(𝑚) = 𝑎𝑛 𝑚 + 𝑎𝑛−1 𝑚 + ⋯ + 𝑎1 𝑚 + 𝑎0

De (4) se infiere que

a) 𝑃 (𝑚) = 0
b) 𝑦 = 𝑒 𝑚𝑥 es solución de la ecuación (3)

 𝑎𝑛 𝑚𝑛 + 𝑎𝑛−1 𝑚𝑛−1 + ⋯ + 𝑎1 𝑚 + 𝑎0 = 0 se llama ecuación auxiliar de la ecuación

diferencial lineal homogénea de orden n (1).

 El método de solución de una ecuación diferencial lineal homogénea con coeficientes

constantes sigue los siguientes pasos:

PASO 1: Plantear la ecuación auxiliar asociada a la ecuación diferencial lineal homogénea y

hallar las raíces de la misma.

PASO 2: Determinar las soluciones de la ecuación diferencial homogénea basándose en los

siguientes casos:

Caso 1: Si 𝑚1 , 𝑚2 , … . , 𝑚𝑛 son raíces reales y distintas, entonces 𝑒 𝑚1 𝑥 , 𝑒 𝑚2 𝑥 , … , 𝑒 𝑚𝑛 𝑥 son

soluciones linealmente independientes de (3) y la solución general es de la forma
𝑦 = 𝑐1𝑒 𝑚1 𝑥 + 𝑐2 𝑒 𝑚2 𝑥 + ⋯ + 𝑐𝑛 𝑒 𝑚𝑛 𝑥

Caso 2: Si 𝑚1 = 𝑎 + 𝑏𝑖 y 𝑚2 = 𝑎 − 𝑏𝑖 son raíces complejas, entonces 𝑒 (𝑎+𝑏𝑖)𝑥 y 𝑒 (𝑎−𝑏𝑖)𝑥 son

soluciones linealmente independientes de (3), pero son funciones complejas de valor real. Sin
embargo, veamos que es posible, a partir de ellas, obtener soluciones reales linealmente
independientes.

Aplicando la fórmula de Euler, tenemos que:

𝑒 𝑚1 𝑥 = 𝑒 (𝑎+𝑏𝑖)𝑥 = 𝑒 𝑎𝑥 (𝑐𝑜𝑠𝑏𝑥 + 𝑖𝑠𝑒𝑛𝑏𝑥 ) = 𝑒 𝑎𝑥 𝑐𝑜𝑠𝑏𝑥 + 𝑖𝑒 𝑎𝑥 𝑠𝑒𝑛𝑏𝑥 (𝑎)

𝑒 𝑚2 𝑥 = 𝑒 (𝑎−𝑏𝑖)𝑥 = 𝑒 𝑎𝑥 (𝑐𝑜𝑠𝑏𝑥 − 𝑖𝑠𝑒𝑛𝑏𝑥 ) = 𝑒 𝑎𝑥 𝑐𝑜𝑠𝑏𝑥 − 𝑖𝑒 𝑎𝑥 𝑠𝑒𝑛𝑏𝑥 (𝑏)

 Lorena Maritza Terrios Guzmán 51
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

Sumando (a) y (b) obtenemos:

1 𝑚1𝑥 1
2
𝑒 + 2𝑒 𝑚2 𝑥 = 𝑒 𝑎𝑥 𝑐𝑜𝑠𝑏𝑥
Restando (a) y (b) obtenemos:
1 𝑚1 𝑥 1 𝑚2 𝑥
2𝑖
𝑒 − 2𝑖
𝑒 = 𝑒 𝑎𝑥 𝑠𝑒𝑛𝑏𝑥

Luego, 𝑒 𝑎𝑥 𝑐𝑜𝑠𝑏𝑥, 𝑒 𝑎𝑥 𝑠𝑒𝑛𝑏𝑥 son solución de la ecuación diferencial dada, por ser combinación
lineal de dos soluciones de dicha ecuación y además son linealmente independiente.

Por lo tanto, la solución general contiene una función de la forma

𝑒 𝑎𝑥 (𝑐1𝑐𝑜𝑠𝑏𝑥 + 𝑐2𝑠𝑒𝑛𝑏𝑥)

Caso 3: Si 𝑚 es una raíz de multiplicidad k, entonces las k soluciones linealmente

independientes correspondientes son 𝑒 𝑚𝑥 , 𝑥𝑒 𝑚𝑥 , … , 𝑥 𝑘−1 𝑒 𝑚𝑥 . La solución general contiene
una expresión de la forma

𝑐1𝑒 𝑚𝑥 + 𝑐2𝑥𝑒 𝑚𝑥 + ⋯ + 𝑐𝑘 𝑥 𝑘−1 𝑒 𝑚𝑥 = 𝑒 𝑚𝑥 (𝑐1 + 𝑐2 𝑥 + ⋯ + 𝑐𝑘 𝑥 𝑘−1 )

EJEMPLO 3
Utilicemos el método descrito para hallar el conjunto fundamental de soluciones de la
siguiente ecuación diferencial.
𝑦 (4) − 13𝑦 ′′ + 36𝑦 = 0 (1)
Solución:

PASO 1: Plantear la ecuación auxiliar asociada a la ecuación diferencial lineal homogénea y

hallar las raíces de la misma.
𝑦 = 𝑒 𝑚𝑥 (2) 𝑦 ′ = 𝑚𝑒 𝑚𝑥 (3) 𝑦 ′′ = 𝑚2 𝑒 𝑚𝑥 (4)
3 𝑚𝑥 (4) 4 𝑚𝑥
𝑦′′′ = 𝑚 𝑒 (5) 𝑦 =𝑚 𝑒 (6)

Sustituyamos (6), (4) y (2) en (1):

𝑚4 𝑒 𝑚𝑥 − 13𝑚2 𝑒 𝑚𝑥 + 36𝑒 𝑚𝑥 = 𝑒 𝑚𝑥 (𝑚4 − 13𝑚2 + 36) = 0 (7)

De (7) se concluye que
⏟𝑚4 − 13𝑚2 + 36 = 0 (8)
𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 𝐴𝑢𝑥𝑖𝑙𝑖𝑎𝑟
Hallemos las raíces (7).

𝑚4 − 13𝑚2 + 36 = (𝑚2 − 9)(𝑚2 − 4) = (𝑚 − 3)(𝑚 + 3)(𝑚 − 2)(𝑚 + 2) = 0

𝑚1 = 3, 𝑚2 = −3, 𝑚3 = 2, 𝑚3 = −2

PASO 2: Como las raíces de la ecuación auxiliar son todas reales y distintas estamos en el caso
1 que describe el método, por lo tanto la solución general de la ecuación diferencial dada es:

𝑦 = 𝑐1 𝑒 3𝑥 + 𝑐2 𝑒 −3𝑥 + 𝑐3 𝑒 2𝑥 + 𝑐4 𝑒 −2𝑥
 Lorena Maritza Terrios Guzmán 52
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

EJERCICIO 4
1. Obtenga la solución general de cada ecuación diferencial.

a) 𝑦 (6) − 𝑦 (4) = 0 b) 𝑦 ′′′ + 4𝑦 ′′ + 5𝑦 ′ = 0

2. Resuelva el siguiente problema con condición inicial.

4𝑦 ′′ + 4𝑦 ′ + 17𝑦 = 0, 𝑦(0) = −1, 𝑦 ′ (0) = 2

3. Encuentre la solución del problema de valor frontera dado.

𝑦 ′′′ − 2𝑦 ′′ + 𝑦 ′ = 0, 𝑦(0) = 0, 𝑦(−1) = −1, 𝑦 ′′(0) = 3

4.3 ECUACIONES DIFERENCIALES NO HOMOGÉNEAS

DEFINICIÓN 4.5: ECUACIONES DIFERENCIALES NO HOMOGÉNEAS. Una ecuación

diferencial, ordinaria, lineal de orden n de la forma
𝑑𝑛 𝑦 𝑑𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥 ) 𝑑𝑥𝑛 + 𝑎𝑛−1 (𝑥 ) 𝑑𝑥𝑛−1 + ⋯ + 𝑎1 (𝑥 ) 𝑑𝑥 + 𝑎0 (𝑥 )𝑦 = 𝑔(𝑥 ) (2)

donde 𝑔(𝑥) ≠ 0, se llama ecuación diferencial no homogénea.

NOTA:
 La ecuación (1) de la definición 4.1 se llama ecuación diferencial homogénea asociada a la
ecuación diferencial no homogénea de la definición 4.5.

 Otra forma de denotar la ecuación diferencial de la definición 4.5 es

𝐹(𝑥, 𝑦, 𝑦´, 𝑦 ′′, … , 𝑦 (𝑛−1) , 𝑦 (𝑛) ) = 𝑔(𝑥 )

 Cualquier función 𝑦𝑝 libre de parámetros arbitrarios, que satisface una ecuación diferencial
lineal no homogénea (2) se dice que es una solución particular de la ecuación.

TEOREMA 4.5: EXISTENCIA DE UNA SOLUCIÓN ÚNICA. Sean 𝑎0 (𝑥 ), 𝑎1 (𝑥 ), … , 𝑎𝑛−1 (𝑥 ), 𝑎𝑛 (𝑥 )

y 𝑔(𝑥) funciones continuas en un intervalo I, y sea 𝑎𝑛 (𝑥 ) ≠ 0 para todo 𝑥 de I. Si 𝑥0 ∈ 𝐼,
entonces existe una solución única 𝑦 = 𝑓 (𝑥 ) del problema con condiciones iniciales definido
como
𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥 ) 𝑛 + 𝑎𝑛−1 (𝑥 ) 𝑛−1 + ⋯ + 𝑎1 (𝑥 ) + 𝑎0 (𝑥 )𝑦 = 𝑔(𝑥 )
𝑑𝑥 𝑑𝑥 𝑑𝑥
𝑦(𝑥0 ) = 𝑦0 , 𝑦 ′ (𝑥0 ) = 𝑦1 , 𝑦 ′′(𝑥0 ) = 𝑦2 , … , 𝑦 (𝑛−1) (𝑥0 ) = 𝑦𝑛−1
 Lorena Maritza Terrios Guzmán 53
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

NOTA:
Las hipótesis del teorema 4.5 son importantes dado que si una de ellas no se cumple podría
ocurrir el problema con condiciones iniciales no tenga solución única o que no tenga solución.

EJERCICIO 3
Analizar la solución dada a cada problema con condiciones iniciales, según el teorema 4.5.
a) 𝑦 = 3𝑒 2𝑥 + 𝑒 −2𝑥 − 3𝑥: solución de (1)

𝑦 ′′ − 4𝑦 = 12𝑥
} (1)
𝑦(0) = 4, 𝑦 ′(0) = 1

b) 𝑦 = 𝑐𝑥 2 + 𝑥 + 3: solución de (2)

𝑥 2 𝑦 ′′ − 2𝑥𝑦 ′ + 2𝑦 = 6
} (2)
𝑦(0) = 3, 𝑦 ′(0) = 1

TEOREMA 4.6: SOLUCIÓN GENERAL DE UNA ECUACIÓN DIFERENCIAL NO HOMOGÉNEA.

Sea 𝑦𝑝 una solución particular de una ecuación diferencial lineal no homogénea de orden n en
un intervalo I, y sea {𝑦1 , 𝑦2 , … , 𝑦𝑛 } un conjunto fundamental de soluciones de la ecuación
diferencial lineal homogénea asociada en I. Entonces la solución general de la ecuación en el
intervalo es
𝑦 = 𝑐1 𝑦1 (𝑥 ) + 𝑐2𝑦2 (𝑥 ) + ⋯ + 𝑐𝑛 𝑦𝑛 (𝑥 ) + 𝑦𝑝
donde 𝑐1 , 𝑐2, … , 𝑐𝑛 son constantes arbitrarias.

NOTA:
Del teorema 4.6 a la combinación lineal 𝑦𝑐 = 𝑐1𝑦1 (𝑥 ) + 𝑐2𝑦2 (𝑥 ) + ⋯ + 𝑐𝑛 𝑦𝑛 (𝑥 ) se le llama
solución complementaria para la ecuación lineal no homogénea, en consecuencia la solución
general de una ecuación diferencial lineal no homogénea es 𝑦 = 𝑦𝑐 + 𝑦𝑝 .

4.4 MÉTODOS DE SOLUCIÓN PARA UNA ECUACIÓN DIFERENCIAL

LINEAL NO HOMOGÉNEA CON COEFICIENTES CONSTANTES

1. MÉTODO DE LOS COEFICIENTES INDETERMINADOS

Este método permite encontrar la solución de una ecuación diferencial lineal no homogénea
con coeficientes constantes. Del teorema 4.6 se tiene que la solución general de la ecuación no
homogénea (2) es de la forma 𝑦 = 𝑦𝑐 + 𝑦𝑝 , donde 𝑦𝑐 es la solución de la ecuación diferencial
lineal homogénea asociada.
La idea fundamental de este método se basa en lanzar una conjetura acerca de la forma de 𝑦𝑝,
motivada por las clases de funciones que definen a la función de entrada 𝑔(𝑥). El método
general se limita a E.D. lineales no homogéneas tales que los coeficientes 𝑎𝑖 , 𝑖 = 1,2, … , 𝑛 son
 Lorena Maritza Terrios Guzmán 54
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

constantes y 𝑔(𝑥) es una función constante o polinomial o exponencial de la forma 𝑒 𝑎𝑥 o

trigonométrica de la forma 𝑠𝑒𝑛𝑏𝑥 o 𝑐𝑜𝑠𝑏𝑥 o sumas y productos finitos de estas funciones.
Estas funciones tienen la notable propiedad de que las derivadas de sus sumas y productos
son de nuevo sumas y productos de funciones con estas mismas características.

(𝑛) (𝑛−1)
Dado que la combinación lineal de las derivadas 𝑎𝑛 𝑦𝑝 + 𝑎𝑛−1 𝑦𝑝 + 𝑎1 𝑦𝑝′ … + 𝑎0 𝑦𝑝 (siendo
𝑦𝑝 solución de la E.D.L. no homogénea), deber ser idéntica a la función de entrada 𝑔(𝑥),
parece razonable suponer que 𝑦𝑝 tiene la misma forma que 𝑔(𝑥).

Los pasos a seguir para utilizar el método son los siguientes:

PASO 1: Hallar la solución de la ecuación diferencial homogénea asociada.

PASO 2: Establecer la conjetura respecto a la solución particular considerando los siguientes

aspectos:
- La solución particular debe tener la forma de la función de entrada.
- Analizar las derivadas sucesivas de 𝑦𝑝 que se deben sustituir en la ecuación
diferencial. Es decir si no aparecen en la conjetura establecida, entonces se adicionan.
- No debe contener ningún término de la solución complementaria de la ecuación
diferencial. Si contiene un término de la solución complementaria de la ecuación
diferencial, entonces se multiplica a dicho término de la solución particular elegida por
𝑥 𝑘 , donde k es el menor entero que elimina esta duplicación.

PASO 3: Establecer y resolver el sistema de ecuaciones que permita hallar las constantes
determinadas en el paso 2.

PASO 4: Describir la solución general de la ecuación diferencial dada como 𝑦 = 𝑦𝑐 + 𝑦𝑝 .

EJEMPLO 4
Resolver la siguiente ecuación diferencial aplicando el método de los coeficientes
indeterminados.
𝑦 ′′ − 3𝑦 ′ + 2𝑦 = ⏟
2𝑥 − 3𝑒 𝑥 + 12𝑥𝑒 −𝑥 (1)
𝑔(𝑥)
Solución:

PASO 1: Hallar la solución de la ecuación diferencial homogénea asociada a (1):

𝑦 ′′ − 3𝑦 ′ + 2𝑦 = 0 (2)
La ecuación auxiliar de (2) es:
𝑚2 − 3𝑚 + 2 = (𝑚 − 1)(𝑚 − 2) = 0 (3)

Las soluciones de (3) son: 𝑚1 = 1, 𝑚2 = 2

Luego la solución complementaria se define como

𝑦𝑐 = 𝑐1 𝑒⏟𝑥 + 𝑐2 𝑒⏟
2𝑥

𝑦1 𝑦2
 Lorena Maritza Terrios Guzmán 55
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

PASO 2: Planteamiento de la conjetura.

𝑦𝑝 : debe tener la misma forma de la función de entrada, 𝑔(𝑥)

𝑦𝑝 = 𝐴𝑥 + 𝐵𝑒 𝑥 + 𝐶𝑥𝑒 −𝑥
𝑦𝑝 : Adicionar (si es necesario) las derivadas sucesivas que se obtienen de 𝑦𝑝 .

𝑦𝑝 = 𝐴𝑥 + 𝐵𝑒 𝑥 + 𝐶𝑥𝑒 −𝑥 + 𝐷
⏟ 𝐹𝑒 −𝑥
+ ⏟
↑
𝑑𝑢𝑝𝑙𝑖𝑐𝑎 𝑟𝑒𝑠𝑢𝑙𝑡𝑎 𝑎𝑙 𝑟𝑒𝑠𝑢𝑙𝑡𝑎
𝑎 𝑦1 𝑑𝑒𝑟𝑖𝑣𝑎𝑟 𝑥 𝑎𝑙
𝑑𝑒𝑟𝑖𝑣𝑎𝑟
𝑥𝑒 −𝑥
𝑦𝑝 : Eliminar la duplicación con la solución complementaria y con los términos que describen
a la solución particular
𝑦𝑝 = 𝐴𝑥 + 𝐵𝑥𝑒 𝑥 + 𝐶𝑥𝑒 −𝑥 + 𝐷 + 𝐹𝑒 −𝑥 (4)
La ecuación (4) describe de manera correcta la forma de la solución particular.
De (4) tenemos que:
𝑦′𝑝 = 𝐴 + 𝐵(𝑥 + 1)𝑒 𝑥 + 𝐶 (1 − 𝑥 )𝑒 −𝑥 − 𝐹𝑒 −𝑥 (5)
𝑦′′𝑝 = 𝐵(𝑥 + 2)𝑒 𝑥 + 𝐶 (𝑥 − 2)𝑒 −𝑥 + 𝐹𝑒 −𝑥 (6)

PASO 3: Hallemos A, B, C y F.
- Sustituyendo (4), (5) y (6) en (1) tenemos:

𝐵(𝑥 + 2)𝑒 𝑥 + 𝐶(𝑥 − 2)𝑒 −𝑥 + 𝐹𝑒 −𝑥 − 3 [𝐴

⏟ ⏟ + 𝐵(𝑥 + 1)𝑒 𝑥 + 𝐶(1 − 𝑥)𝑒 −𝑥 − 𝐹𝑒 −𝑥 ] + 2 [𝐴𝑥
⏟ + 𝐵𝑥𝑒 𝑥 + 𝐶𝑥𝑒 −𝑥 + 𝐷 + 𝐹𝑒 −𝑥 ] = 𝑔(𝑥)
𝑦′′𝑝 𝑦′𝑝 𝑦𝑝

- Para obtener el sistema de ecuaciones en las variables A, B, C y F se igualan los coeficientes

de los términos semejantes de cada lado de la última igualdad.

𝑥: 2𝐴 = 2 ⟺ 𝐴 = 1
𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒𝑠: − 3𝐴 + 2𝐷 = 0 ⟺ 𝐷 = 3/2
𝑒𝑥: 2𝐵 − 3𝐵 = −3 ⟺ 𝐵=3
𝑥
𝑥𝑒 : 𝐵 − 3𝐵 + 2𝐵 = 0 ⟺ 0 = 0
𝑥𝑒 −𝑥 : 𝐶 + 3𝐶 + 2𝐶 = 12 ⟺ 𝐶=2
−𝑥
𝑒 : − 2𝐶 + 𝐹 − 3𝐶 + 3𝐹 + 2𝐹 = 0 ⟺ −5𝐶 + 6𝐹 = 0 ⟺ 𝐹 = 5/3
Luego,
3 5
𝑦𝑝 = 𝑥 + 3𝑥𝑒 𝑥 + 2𝑥𝑒 −𝑥 + + 𝑒 −𝑥
2 3
PASO 4: La solución general de la ecuación diferencial dada es
3 5 −𝑥
𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1𝑒 𝑥 + 𝑐2 𝑒 2𝑥 + 𝑥 + 3𝑥𝑒 𝑥 + 2𝑥𝑒 −𝑥 + + 𝑒
2 3
EJERCICIO 5
Obtenga la solución de cada ecuación diferencial aplicando el método de los coeficientes
indeterminados.

a) 𝑦 ′′ − 8𝑦´ + 25𝑦 = 5𝑥 2 𝑒 −𝑥 b) 𝑦 ′′′ + 𝑦 ′′ = 𝑒 𝑥 𝑐𝑜𝑠𝑥

 Lorena Maritza Terrios Guzmán 56
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

TAREA 2: Plantee y halle la solución particular de cada ecuación diferencial según el método
de los coeficientes indeterminado.

a) 4𝑦 ′′′ − 20𝑦 ′′ − 𝑦 ′ + 5𝑦 = 5𝑥 3 𝑒 5𝑥 + 𝑠𝑒𝑛𝑥 b) 𝑦 (4) − 4𝑦 ′′′ + 6𝑦 ′′ − 4𝑦 ′ + 𝑦 = 𝑥 + 𝑒 𝑥

2. MÉTODO DEL ANULADOR

El método del anulador es un procedimiento que permite determinar una solución particular
𝑦𝑝 , de una ecuación diferencial lineal no homogénea 𝐿(𝑦) = 𝑔(𝑥 ), con coeficientes constantes
donde la función de entrada consiste en sumas y productos finitos de las funciones constantes,
polinomiales, exponenciales y trigonométricas de la forma

𝑘, 𝑥 𝑛 , 𝑥 𝑛 𝑒 ∝𝑥 , 𝑥 𝑛 𝑒 ∝𝑥 𝑠𝑒𝑛𝛽𝑥, 𝑥 𝑛 𝑒 ∝𝑥 𝑐𝑜𝑠𝛽𝑥

Para abordar el método del anulador es necesario considerar los siguientes aspectos.

DEFINICIÓN 4.6: OPERADOR DIFERENCIAL. Sea 𝑦 = 𝑓(𝑥) 𝑛 veces derivable. Un operador

diferencial lineal de orden 𝑛 es una función 𝐿 definida como

𝐿(𝑦) = 𝑎𝑛 𝐷 𝑛 𝑦 + 𝑎𝑛−1 𝐷 𝑛−1 𝑦 + ⋯ + 𝑎1 𝐷𝑦 + 𝑎0 𝑦 = (𝑎𝑛 𝐷 𝑛 + 𝑎𝑛−1 𝐷 𝑛−1 + ⋯ + 𝑎1 𝐷 + 𝑎0 )𝑦 (1)

𝑑𝑛 𝑦 𝑑𝑛−1 𝑦 𝑑𝑦
Donde 𝐷 𝑛 𝑦 = 𝑑𝑥𝑛 , 𝐷 𝑛−1 𝑦 = 𝑑𝑥𝑛−1 , … , 𝐷𝑦 = 𝑑𝑥

NOTA:
𝑑𝑛 𝑦 𝑑𝑛−1 𝑦 𝑑𝑦
 Si 𝑎𝑛 𝑑𝑥𝑛 + 𝑎𝑛−1 𝑑𝑥𝑛−1 + ⋯ + 𝑎1 𝑑𝑥 + 𝑎0 𝑦 = 𝑔(𝑥 ) (2), entonces por la definición 4.6
tenemos que (2) se puede expresar como

𝑎𝑛 𝐷 𝑛 𝑦 + 𝑎𝑛−1 𝐷 𝑛−1 𝑦 + ⋯ + 𝑎1 𝐷𝑦 + 𝑎0 𝑦 = 𝐿(𝑦) = 𝑔(𝑥 ) (3)

 Si cada 𝑎𝑖 en (3) es constante, entonces un operador diferencial lineal se puede factorizar,

siempre que el polinomio 𝑎𝑛 𝑚𝑛 + 𝑎𝑛−1 𝑚𝑛−1 + ⋯ + 𝑎1 𝑚 + 𝑎0 sea factorizable, es decir si 𝑟 es
una raíz de la ecuación auxiliar, entonces
𝐿 = (𝐷 − 𝑟)𝑃(𝐷)
donde 𝑃(𝐷) es un operador diferencial lineal de orden 𝑛 − 1.

Por ejemplo: sea 𝑚2 + 5𝑚 + 6 = (𝑚 + 3)(𝑚 + 2)

Luego si 𝑦 = 𝑓(𝑥) es dos veces derivable, entonces (𝐷 2 + 5𝐷 + 6)𝑦 = (𝐷 + 3)(𝐷 + 2)𝑦

PROPIEDADES DE LOS OPERADORES DIFERENCIALES: Sean 𝑦, 𝑦1 y 𝑦2 funciones

suficientemente derivables en una variable y, 𝐿, 𝐿1 y 𝐿2 operadores diferenciales lineales,
entonces
a) 𝐿2 (𝑦) + 𝐿1 (𝑦) = 𝐿1 (𝑦) + 𝐿2 (𝑦)
 Lorena Maritza Terrios Guzmán 57
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

b) 𝐿2 𝐿1 (𝑦) = 𝐿1 𝐿2 (𝑦)
c) 𝐿(𝑘𝑦) = 𝑘𝐿(𝑦), 𝑘 ∈ ℝ
d) 𝐿(𝑦1 + 𝑦2 ) = 𝐿(𝑦1 ) + 𝐿(𝑦2 )

DEFINICIÓN 4.7: OPERADOR DIFERENCIAL ANULADOR. Si L es un operador diferencial

lineal con coeficientes constantes y 𝑓 es una función suficientemente derivable tal que
𝐿(𝑓(𝑥 )) = 0
entonces se dice que L es un anulador de la función f.

Por ejemplo D anula a 𝑓 (𝑥 ) = 𝑘; 𝐷 2 anula a 𝑓 (𝑥 ) = 2𝑥 + 3; 𝐷 3 anula a 𝑓 (𝑥 ) = 𝑥 2 .

NOTA:
 Sea 𝑃𝑛−1 (𝑥 ) = 𝑎𝑛−1 𝑥 𝑛−1 + 𝑎𝑛−2 𝑥 𝑛−2 + ⋯ + 𝑎1 𝑥 + 𝑎0 .
Se distinguen los siguientes tipos de operadores diferenciales anuladores:

TIPO 1: 𝐷 𝑛 [𝑃𝑛−1 (𝑥 )] = 0

TIPO 2: (𝐷−∝)𝑛 [𝑃𝑛−1 (𝑥 )𝑒 ∝𝑥 ] = 0

TIPO 3: [𝐷 2 − 2 ∝ 𝐷 + (∝2 + 𝛽2 )]𝑛 [𝑃𝑛−1 (𝑥 ) 𝑒 ∝𝑥 𝑐𝑜𝑠𝛽𝑥 ] = 0

[𝐷 2 − 2 ∝ 𝐷 + (∝2 + 𝛽2 )]𝑛 [𝑃𝑛−1 (𝑥 ) 𝑒 ∝𝑥 𝑠𝑒𝑛𝛽𝑥 ] = 0

PROPIEDADES DE LOS OPERADORES DIFERENCIALES ANULADORES:

a) Si 𝐿 es un operador diferencial anulador tal que 𝐿(𝑦1 ) = 0 y 𝐿(𝑦2 ) = 0, entonces
𝐿(𝑐1𝑦1 + 𝑐2 𝑦2) = 0.
b) Si 𝐿1 y 𝐿2 son operadores diferenciales anuladores tales que 𝐿1 (𝑦1 ) = 0 y 𝐿2 (𝑦2 ) = 0,
entonces 𝐿1 𝐿2 (𝑦1 + 𝑦2 ) = 0

NOTA:
 Si 𝑦1 𝑦 𝑦2 tienen el mismo anulador L, entonces L es un anulador de menor orden de
𝑦1 + 𝑦2 .

 Si L es un anulador de 𝑦1 entonces 𝐿𝑛 , 𝑛 ∈ ℕ, 𝑛 > 1 es también un anulador de 𝑦1 .

 (𝐿1 )𝑛 (𝐿2 )𝑚 (𝑦1 ) = 0 y (𝐿1 )3𝑛 (𝐿2 )𝑚 (𝐿3)𝑘 (𝑦2 ) = 0 ⇒ (𝐿1 )3𝑛 (𝐿2 )𝑚 (𝐿3 )𝑘 (𝑦1 + 𝑦2 ) = 0

EJEMPLO 5
Determinemos el operador diferencial anulador de la función 𝑓(𝑥 ) = 4𝑒 2𝑥 − 10𝑥𝑒 2𝑥

Solución:
Dado que
𝑓(𝑥 ) = 4𝑒 2𝑥 − 10𝑥𝑒 2𝑥 = 𝑒 2𝑥 (4 − 10𝑥 ) (1)
 Lorena Maritza Terrios Guzmán 58
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

Vemos que ∝= 2 y la función polinomial es de grado 1, entonces el operador anulador para f

es tipo 2. Es decir es de la forma (𝐷 − 2)2

(𝐷 − 2)2 [𝑒 2𝑥 (4 − 10𝑥)] = 0 (2)

Verifiquemos que se satisface la igualdad anterior.

(𝐷 − 2)2 [𝑒 2𝑥 (4 − 10𝑥)] = (𝐷 2 − 4𝐷 + 4)[𝑒 2𝑥 (4 − 10𝑥)] = 0

(𝐷 − 2)2 [𝑒 2𝑥 (4 − 10𝑥)] = 𝐷 2 [𝑒 2𝑥 (4 − 10𝑥)] − 4𝐷 [𝑒 2𝑥 (4 − 10𝑥 )] + 4[𝑒 2𝑥 (4 − 10𝑥 )] = 0 (3)

𝐷[𝑒 2𝑥 (4 − 10𝑥 )] = 𝑒 2𝑥 (−10) + (4 − 10𝑥 )2𝑒 2𝑥 = −2𝑒 2𝑥 (1 + 10𝑥 ) (4)

𝐷 2 [𝑒 2𝑥 (4 − 10𝑥)] = −2[(1 + 10𝑥 )2𝑒 2𝑥 + 𝑒 2𝑥 (10)] = −8𝑒 2𝑥 (3 + 5𝑥 ) (5)

Sustituyamos (1), (4) y (5), en (2):

(𝐷 − 2)2 [𝑒 2𝑥 (4 − 10𝑥)] = −8𝑒 2𝑥 (3 + 5𝑥 ) + 8𝑒 2𝑥 (1 + 10𝑥 ) + 8[𝑒 2𝑥 (2 − 5𝑥 )]

= 8𝑒 2𝑥 (−3 − 5𝑥 + 1 + 10𝑥 + 2 − 5𝑥 ) = 𝑒 2𝑥 (0) = 0

EJERCICIO 6
Determinemos el operador diferencial que anule a la función 𝑓(𝑥 ) = 5𝑐𝑜𝑠2𝑥 + 9𝑒 𝑥 𝑠𝑒𝑛3𝑥

TAREA 3: Demuestre que:

a) Si L es un anulador de 𝑦1, también lo es 𝐿𝑛 , 𝑛 ∈ ℕ, 𝑛 > 1.
b) Si 𝐿1 𝐿2 (𝑦1 + 𝑦2 ) = 0, entonce 𝐿1 𝐿2 (𝑦1 ) = 0

Los pasos a seguir para utilizar el método del anulador son los siguientes:

PASO 1: Hallar la solución de la ecuación diferencial lineal homogénea asociada 𝐿(𝑦) = 0

PASO 2: Determinar el operador diferencial de menor orden que anula a 𝑔(𝑥 ), 𝐿1 (𝑔(𝑥)) = 0 y
aplicarlo a la ecuación diferencial no homogénea, esto es

𝐿1 𝐿(𝑦) = 𝐿1 (𝑔(𝑥)) = 0

PASO 3: Resolver la E.D.L.H. 𝐿1 𝐿(𝑦) = 0

En este paso se descubre la forma de una solución particular 𝑦𝑝 para la ecuación diferencial
lineal no homogénea.

PASO 4: Eliminar de la solución del paso 3 los términos que se duplican en la solución
complementaria. Formar una combinación lineal 𝑦𝑝 de los términos restantes.

PASO 5: Se sustituye esta forma supuesta en 𝐿(𝑦) = 𝑔(𝑥 ), para encontrar 𝑦𝑝 explícitamente.
 Lorena Maritza Terrios Guzmán 59
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

EJEMPLO 6
Apliquemos el método del anulador para resolver la siguiente ecuación diferencial no
homogénea.
𝑦 ′′ − 2𝑦 ′ + 𝑦 = 10𝑒 −2𝑥 𝑐𝑜𝑠𝑥 ⟺ (𝐷 2 − 2𝐷 + 1)𝑦 = 10𝑒 −2𝑥 𝑐𝑜𝑠𝑥 (1)

Solución:

PASO 1: Hallar la solución de

𝑦 ′′ − 2𝑦 ′ + 𝑦 = (𝐷 2 − 2𝐷 + 1)𝑦 = 0

Establezcamos la ecuación auxiliar y hallemos sus raíces:

𝑚2 − 2𝑚 + 1 = (𝑚 − 1)2 = 0

La raíz es 𝑚 = 1 de multiplicidad algebraica 2.

La solución complementaria es entonces:

𝑦𝑐 = 𝑐1 𝑒 𝑥 + 𝑐2 𝑥𝑒 𝑥

PASO 2: Identificar el operador anulador de la función de entrada.

𝑛 ∝ 𝛽

𝑐𝑜𝑠 ⏞
⏞ ⏞𝑥
Como 𝑔(𝑥 ) = 10𝑒 −2𝑥 𝑐𝑜𝑠𝑥 = 10𝑥 𝑒 0 −2
1∙𝑥

Por lo tanto el operador anulador es tipo 3, esto es

[𝐷 2 − 2(−2)𝐷 + ((−2)2 + (1)2 )]1 [10𝑒 −2𝑥 𝑐𝑜𝑠𝑥 ] = 0

(𝐷 2 + 4𝐷 + 5)[10𝑒 −2𝑥 𝑐𝑜𝑠𝑥 ] = 0

PASO 3: Resolver la ecuación diferencial homogénea que se obtiene de aplicar en (1), el

anulador encontrado en el paso 2.
(𝐷 2 + 4𝐷 + 5)(𝐷 2 − 2𝐷 + 1)𝑦 = (𝐷 2 + 4𝐷 + 5)[10𝑒 −2𝑥 𝑐𝑜𝑠𝑥 ]

(𝐷 2 + 4𝐷 + 5)(𝐷 2 − 2𝐷 + 1)𝑦 = 0 (2)

Hallemos y resolvamos la ecuación auxiliar de (2):

(𝑚2 + 4𝑚 + 5)(𝑚2 − 2𝑚 + 1) = (𝑚2 + 4𝑚 + 5)(𝑚 − 1)2 = 0

𝑚1 = −2 + 𝑖, 𝑚2 = −2 − 𝑖 , 𝑚3 = 1 ∶ 𝑚𝑢𝑙𝑡𝑖𝑝𝑙𝑖𝑐𝑖𝑑𝑎𝑑 𝑎𝑙𝑔𝑒𝑏𝑟𝑎𝑖𝑐𝑎 2

La solución de (2) es entonces:

𝑐1 𝑒 𝑥 + 𝑐2 𝑥𝑒 𝑥 + ⏟
𝑦=⏟ 𝑒 −2𝑥 (𝑐3𝑐𝑜𝑠𝑥 + 𝑐4 𝑠𝑒𝑛𝑥 )
𝑦𝑐 𝑦𝑝
PASO 4: Al eliminar los términos que duplican la solución complementaria obtenemos:
 Lorena Maritza Terrios Guzmán 60
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

𝑦𝑝 = 𝑒 −2𝑥 (𝐴𝑐𝑜𝑠𝑥 + 𝐵𝑠𝑒𝑛𝑥 ) (3)

PASO 5: Hallemos las constantes

𝑦′𝑝 = 𝑒 −2𝑥 (−𝐴𝑠𝑒𝑛𝑥 + 𝐵𝑐𝑜𝑠𝑥 ) − 2𝑒 −2𝑥 (𝐴𝑐𝑜𝑠𝑥 + 𝐵𝑠𝑒𝑛𝑥 )

𝑦′𝑝 = 𝑒 −2𝑥 [(−𝐴 − 2𝐵)𝑠𝑒𝑛𝑥 + (𝐵 − 2𝐴)𝑐𝑜𝑠𝑥 ] (4)

𝑦′′𝑝 = 𝑒 −2𝑥 [(3𝐴 − 4𝐵)𝑐𝑜𝑠𝑥 + (3𝐵 + 4𝐴)𝑠𝑒𝑛𝑥 ] (5)

Sustituyendo (3), (4) y (5) en (1):

𝑦′′𝑝 − 2𝑦 ′𝑝 + 𝑦𝑝 = ⏟
10𝑒 −2𝑥 𝑐𝑜𝑠𝑥
𝑔(𝑥)
−2𝑥 −2𝑥 −2𝑥
(8𝐴 − 6𝐵)𝑒 𝑐𝑜𝑠𝑥 + (6𝐴 + 8𝐵)𝑒 𝑠𝑒𝑛𝑥 = 10𝑒 𝑐𝑜𝑠𝑥 (6)

Igualando los coeficientes de los términos semejantes de cada lado de (6) se obtiene:

𝑒 −2𝑥 𝑐𝑜𝑠𝑥: 8𝐴 − 6𝐵 = 10 (7)

𝑒 −2𝑥 𝑠𝑒𝑛𝑥: 6𝐴 + 8𝐵 = 0 (8)
4𝐵
De (8), 𝐴 = − 3

Sustituyendo A en (7) obtenemos

4𝐵 50 3 4
8 (− ) − 6𝐵 = 10 ⟺ − 𝐵 = 10 ⟺ 𝐵 = − ⟹ 𝐴 =
3 3 5 5

Por lo tanto la solución general es:

𝑒 −2𝑥
𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1𝑒 𝑥 + 𝑐2𝑥𝑒 𝑥 + (4𝑐𝑜𝑠𝑥 − 3𝑠𝑒𝑛𝑥 )
5
EJERCICIO 7
Aplique el método del anulador para resolver cada ecuación diferencial no homogénea.

a) 𝑦 ′′ − 3𝑦 ′ = 8𝑒 3𝑥 + 4𝑠𝑒𝑛𝑥 b) 𝑦 ′′′ − 𝑦 ′′ = 2 + 𝑥 + 𝑒 2𝑥

TAREA 4: Describa y halle la solución particular de cada ecuación diferencial, según el método
del anulador.
a) 𝑦 ′′′ − 3𝑦 ′′ + 2𝑦 = 5𝑥 3 𝑒 𝑥 b) 𝑦 (4) − 4𝑦 ′′′ + 6𝑦 ′′ − 4𝑦 ′ + 𝑦 = 𝑠𝑒𝑛𝑥 − 2𝑥𝑒 −𝑥
 Lorena Maritza Terrios Guzmán 61
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

3. MÉTODO DE VARIACIÓN DE PARÁMETROS

El método de variación de parámetros es un procedimiento para encontrar la solución

particular 𝑦𝑝 de una ecuación diferencial lineal no homogénea de orden n de la forma:

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 𝑛
+ 𝑎𝑛−1 𝑛−1
+ ⋯ + 𝑎1 + 𝑎0 𝑦 = 𝑔(𝑥 ) (1)
𝑑𝑥 𝑑𝑥 𝑑𝑥
En (1) 𝑔(𝑥) es una función continua en un intervalo I.

El método sigue los siguientes pasos:

PASO 1: Hallar la solución 𝑦𝑐 , de la ecuación homogénea asociada. Sea esta

𝑦𝑐 = 𝑐1𝑦1 + 𝑐2 𝑦2 + ⋯ + 𝑐𝑛 𝑦𝑛 (2)

PASO 2: Si 𝑎𝑛 ≠ 1, dividir entre 𝑎𝑛 en (1) para obtener la forma estándar de la ecuación (1):

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
+ 𝑃𝑛−1 + ⋯ + 𝑃1 + 𝑃0 𝑦 = 𝑓 (𝑥 ) (3)
𝑑𝑥 𝑛 𝑑𝑥 𝑛−1 𝑑𝑥
PASO 3: Una solución particular se obtiene así:

𝑦𝑝 = 𝑦1 𝑢1 (𝑥 ) + 𝑦2 𝑢2 (𝑥 ) + ⋯ + 𝑦𝑛 𝑢𝑛 (𝑥 ) (4)

donde
𝑊𝑘 (𝑥 )
𝑢𝑘 (𝑥 ) = ∫ 𝑑𝑥 , 𝑘 = 1,2, … , 𝑛 (5)
𝑊 (𝑥 )
W es el Wronskiano de 𝑦1 , 𝑦2 , … , 𝑦𝑛, y 𝑊𝑘 es el determinante que se obtiene al reemplazar la
k-ésima columna del Wronskiano por la columna cuyas entradas son (0,0, … ,0, 𝑓(𝑥))

Al resolver cada integral en (5) se omite la constante de integración.

PASO 4: La solución general de (1) es 𝑦 = 𝑦𝑐 + 𝑦𝑝

NOTA:
 El método de variación de parámetros tiene una ventaja particular sobre el método de los
coeficientes indeterminados en cuanto a que es aplicable a ecuaciones no homogéneas donde
la función de entrada no se limita a una función que es combinación lineal de funciones
constantes, polinomiales, exponenciales de la forma 𝑒 𝑎𝑥 o seno y coseno, de a forma
𝑠𝑒𝑛𝑏𝑥 o 𝑐𝑜𝑠𝑏𝑥.

 El método de variación de parámetros siempre produce una solución particular 𝑦𝑝 , siempre

y cuando se pueda resolver la ecuación homogénea asociada.
 Lorena Maritza Terrios Guzmán 62
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

EJEMPLO 7
Apliquemos el método de variación de parámetros para obtener la solución de la siguiente
ecuación diferencial.

2𝑦 ′′′ − 6𝑦 ′′ + 4𝑦′ = ⏟
2(𝑥 − 1)𝑒 𝑥 (1)
𝑔(𝑥)
Solución:

PASO 1: Hallemos la solución complementaria.

2𝑚3 − 6𝑚2 + 4𝑚 = 2𝑚(𝑚2 − 3𝑚 + 2) = 2𝑚(𝑚
⏟ − 1)(𝑚 − 2) = 0
𝑒𝑐𝑢𝑎𝑐𝑖ó𝑛 𝑎𝑢𝑥𝑙𝑖𝑎𝑟

𝑚1 = 1, 𝑚2 = 2, 𝑚3 = 0 ⇔ 𝑦𝑐 = 𝑐1 𝑒⏟𝑥 + 𝑐2 𝑒⏟
2𝑥
+ 𝑐3 ∙ ⏟
1
𝑦1 𝑦2 𝑦3

PASO 2: Dividamos entre 2 para llevar (1) a la forma estándar.

𝑦 ′′′ − 3𝑦 ′′ + 2𝑦′ = ⏟
(𝑥 − 1)𝑒 𝑥
𝑓(𝑥)
PASO 3: describamos y obtengamos la solución particular.

𝑦𝑝 = 𝑦1 𝑢1 (𝑥 ) + 𝑦2 𝑢2 (𝑥 ) + 𝑦3 𝑢3 (𝑥 ) + 𝑦4 𝑢4 (𝑥 ) = 𝑒 𝑥 𝑢1 (𝑥 ) + 𝑒 2𝑥 𝑢2 (𝑥 ) + 1 ∙ 𝑢3 (𝑥 )

𝑒𝑥 𝑒 2𝑥 1
𝑥 2𝑥
𝑊 (𝑥 ) = 𝑊 (𝑒 , 𝑒 , 1) = |𝑒 𝑥 2𝑒 2𝑥 0| = 2𝑒 3𝑥
𝑒𝑥 4𝑒 2𝑥 0
0 𝑒 2𝑥 1
𝑊1 (𝑥) = | 0 2𝑒 2𝑥 0| = 2(1 − 𝑥 )𝑒 3𝑥
(𝑥 − 1)𝑒 𝑥 4𝑒 2𝑥 0

𝑒𝑥 0 1
𝑊2 (𝑥) = |𝑒 𝑥 0 0| = (𝑥 − 1)𝑒 2𝑥
𝑒𝑥 (𝑥 − 1)𝑒 𝑥 0

𝑒𝑥 𝑒 2𝑥 0
𝑊3 (𝑥) = |𝑒 𝑥 2𝑒 2𝑥 0 | = (𝑥 − 1)𝑒 4𝑥
𝑒𝑥 4𝑒 2𝑥 (𝑥 − 1)𝑒 𝑥

𝑊1 (𝑥 ) 2(1 − 𝑥 )𝑒 3𝑥 𝑥2 1
𝑢1 (𝑥 ) = ∫ 𝑑𝑥 = ∫ 𝑑𝑥 = ∫(1 − 𝑥 )𝑑𝑥 = 𝑥 − = (2𝑥 − 𝑥 2 )
𝑊 (𝑥 ) 2𝑒 3𝑥 2 2

𝑊2 (𝑥 ) (𝑥 − 1)𝑒 2𝑥 1 −𝑥
1 −𝑥
𝑢2 (𝑥 ) = ∫ 𝑑𝑥 = ∫ 𝑑𝑥 = ∫(𝑥 − 1)𝑒 𝑑𝑥 = − 𝑥𝑒
𝑊 (𝑥 ) 2𝑒 3𝑥 2 2

𝑊3 (𝑥 ) (𝑥 − 1)𝑒 4𝑥 1 1
𝑢3 (𝑥 ) = ∫ 𝑑𝑥 = ∫ 𝑑𝑥 = ∫(𝑥 − 1)𝑒 𝑥 𝑑𝑥 = (𝑥 − 2)𝑒 𝑥
𝑊 (𝑥 ) 2𝑒 3𝑥 2 2
 Lorena Maritza Terrios Guzmán 63
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

Por lo tanto,
𝑥
2𝑥 − 𝑥 2 2𝑥 (−
1 −𝑥 1 𝑥
−𝑒 𝑥 2
𝑦𝑝 = 𝑒 ( )+𝑒 𝑥𝑒 ) + 1 ∙ [ (𝑥 − 2)𝑒 ] = (𝑥 − 2𝑥 + 2)
2 2 2 2

PASO 4: La solución general de (1) es

𝑒𝑥 2
𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 2𝑥 + 𝑐3 − (𝑥 − 2𝑥 + 2)
2
EJERCICIO 8
Aplicar el método de variación de parámetros para obtener la solución de cada ecuación
diferencial.
𝑒 3𝑥
a) 𝑦 ′′ − 3𝑦 ′ + 2𝑦 = b) 3𝑦 ′′′ − 12𝑦′ = 9𝑡𝑎𝑛𝑥
1+𝑒 𝑥

TAREA 5: Obtenga la solución general de cada ecuación diferencial aplicando el método de

variación de parámetros.

a) 𝑦 ′′ − 𝑦 = 𝑒 −2𝑥 𝑠𝑒𝑛(𝑒 −𝑥 ) b) 𝑦 ′′′ − 𝑦 ′′ = 2 + 𝑥 + 𝑒 2𝑥

4.5 ECUACIONES DIFERENCIALES CON COEFICIENTES VARIABLES:

ECUACIÓN DE CAUCHY – EULER

DEFINICIÓN4.7: ECUACIÓN DE CAUCHY – EULER. Una ecuación diferencial lineal de la forma

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 𝑥 𝑛 𝑛 + 𝑎𝑛−1 𝑥 𝑛−1 𝑛−1 + ⋯ + 𝑎1 𝑥 + 𝑎0 𝑦 = 𝑔(𝑥 ) (1)
𝑑𝑥 𝑑𝑥 𝑑𝑥
donde 𝑎𝑛 , 𝑎𝑛−1 , … , 𝑎0 son constantes, se denomina Ecuación de Cauchy – Euler.

NOTA:
 La característica observable de este tipo de ecuación es que los coeficientes son monomios
en variable 𝑥, cuyo grado coincide con el orden de la derivada a la que multiplican.

 Analicemos una ecuación de Cauchy – Euler de orden 2, esto es

𝑎2 𝑥 2 𝑦 ′′ + 𝑎1 𝑥𝑦 ′ + 𝑎0 𝑦 = 𝑔(𝑥) (1)
Cuya ecuación de Cauchy – Euler de orden 3, homogénea asociada es

𝑎2 𝑥 2 𝑦 ′′ + 𝑎1 𝑥𝑦 ′ + 𝑎0 𝑦 = 0 (2)
- Si
𝑦 = 𝑥 𝑚 ⟹ 𝑦 ′ = 𝑚𝑥𝑚−1 ⟹ 𝑦 ′′ = 𝑚(𝑚 − 1)𝑥 𝑚−2

- Sustituyendo 𝑦, 𝑦 ′ y 𝑦′′ en (2) se tiene que

𝑎2 𝑥 2 𝑚(𝑚 − 1)𝑥 𝑚−2 + 𝑎1 𝑥𝑚𝑥 𝑚−1 + 𝑎0 𝑥 𝑚 = 0
 Lorena Maritza Terrios Guzmán 64
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

𝑎2 𝑚(𝑚 − 1)𝑥 𝑚 + 𝑎1 𝑚𝑥 𝑚 + 𝑎0 𝑥 𝑚 = 0
[𝑎2 𝑚(𝑚 − 1) + 𝑎1 𝑚 + 𝑎0 ]𝑥 𝑚 = 0 (3)

Así 𝑦 = 𝑥 𝑚 es una solución de la ecuación diferencial homogénea asociada a (2) siempre que
m sea una solución de la ecuación auxiliar

𝑎2 𝑚(𝑚 − 1) + 𝑎1 𝑚 + 𝑎0 = 0

𝑎2 𝑚2 + (−𝑎2 + 𝑎1 )𝑚 + 𝑎0 = 0

El proceso anterior permite describir un método de solución para una ecuación de Cauchy-
Euler de orden n.

MÉTODO DE SOLUCIÓN DE UNA ECUACIÓN DE ECUACIÓN DE CAUCHY – EULER

Dada la ecuación de Cauchy-Euler

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 𝑥 𝑛 + 𝑎𝑛−1 𝑥 𝑛−1
+ ⋯ + 𝑎1 𝑥 + 𝑎0 𝑦 = 𝑔(𝑥 ) (1)
𝑑𝑥 𝑛 𝑑𝑥 𝑛−1 𝑑𝑥
y la ecuación homogénea asociada a (1):

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 𝑥 𝑛 𝑛
+ 𝑎𝑛−1 𝑥 𝑛−1
𝑛−1
+ ⋯ + 𝑎1 𝑥 + 𝑎0 𝑦 = 0 (2)
𝑑𝑥 𝑑𝑥 𝑑𝑥
El método de solución de (1) sigue los siguientes pasos:

PASO 1: Sustituir 𝑦 = 𝑥 𝑚 y sus derivadas de orden superior necesarias, en la ecuación de

Cauchy – Euler homogénea asociada.

PASO 2: Obtener la ecuación auxiliar y determinar su solución

𝑃 (𝑚 ) = 0 (3)

PASO 3: La solución de la ecuación de Cauchy – Euler homogénea se obtiene considerando los

siguientes casos.

- Caso 1: Si 𝑚1 , 𝑚2 , … , 𝑚𝑛 son raíces reales y distintas, entonces las soluciones de (2) son de
la forma
𝑥 𝑚1 , 𝑥 𝑚2 , … , 𝑥 𝑚𝑛

- Caso 2: Si 𝑚 es una raíz de multiplicidad algebraica k, entonces las soluciones de (2)

correspondiente a esta raíz son de la forma

𝑥 𝑚 , 𝑥 𝑚 𝑙𝑛𝑥, … , 𝑥 𝑚 (𝑙𝑛𝑥 )𝑘−1

- Caso 3: Si 𝑚1 = 𝑎 + 𝑏𝑖 y 𝑚2 = 𝑎 − 𝑏𝑖 son raíces complejas, entonces las soluciones de (2)

correspondiente a esta raíz son de la forma

𝑥 𝑎 𝑐𝑜𝑠(𝑏𝑙𝑛𝑥 ), 𝑥 𝑎 𝑠𝑒𝑛(𝑏𝑙𝑛𝑥 )
 Lorena Maritza Terrios Guzmán 65
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

PASO 4: Hallar la solución particular utilizando el método de variación de parámetros.

PASO 5: Expresar la solución general de (1) como 𝑦 = 𝑦𝑐 + 𝑦𝑝 .

NOTA:
Si la ecuación de Cauchy-Euler a resolver es homogénea, entonces solo se aplican los tres
primeros pasos del procedimiento anterior.

EJEMPLO 8
Obtengamos la solución de
2𝑥 2 𝑦 ′′ − 5𝑥𝑦 ′ + 5𝑦 = 2𝑥 2 𝑙𝑛𝑥 (1)

Solución:
PASO 1: 2𝑥⏟ 2 𝑦 ′′ − 5𝑥𝑦 ′ + 5𝑦 = 0 (2)
𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝐶𝑎𝑢𝑐ℎ𝑦−𝐸𝑢𝑙𝑒𝑟 𝑎𝑠𝑜𝑐𝑖𝑎𝑑𝑎

𝑦 = 𝑥 𝑚 (3) 𝑦 ′ = 𝑚𝑥 𝑚−1 (4) 𝑦 ′′ = 𝑚(𝑚 − 1)𝑥 𝑚−2 (5)

Sustituyamos (3), (4) y (5) en (2):

2𝑥 2 𝑚(𝑚 − 1)𝑥 𝑚−2 − 5𝑥𝑚𝑥 𝑚−1 + 5𝑥 𝑚 = 2𝑚(𝑚 − 1)𝑥 𝑚 − 5𝑚𝑥 𝑚 + 5𝑥 𝑚 = 0

= [2𝑚(𝑚 − 1) − 5𝑚 + 5]𝑥𝑚
= (2𝑚2 − 7𝑚 + 5)𝑥 𝑚 = 0
PASO 2: La ecuación auxiliar es:
2𝑚2 − 7𝑚 + 5 = (2𝑚 − 5)(𝑚 − 1) = 0

Las raíces de la ecuación auxiliar son:

5
𝑚1 = 1, 𝑚2 =
2

PASO 3: Como las soluciones de la ecuación son reales y distintas aplicamos el caso 1 para
describir la solución de la ecuación diferencial, esto es
5⁄ 2
𝑦𝑐 = 𝑐1 𝑥 𝑚1 + 𝑐2𝑥 𝑚2 = 𝑐1 ⏟
𝑥 1 + 𝑐2 𝑥⏟
𝑦1 𝑦2

PASO 4: Hallar 𝑦𝑝 , aplicando el método de variación de parámetros:

- Llevemos (1) a la forma estándar dividiendo entre 𝑥 2 :

5 ′ 5
𝑦 ′′ − 𝑦 + 2 𝑦 = 𝑙𝑛𝑥 ⏟
2𝑥 2𝑥 𝑓(𝑥)
- La solución particular, según el método de variación de parámetros, es de la forma:

𝑦𝑝 = 𝑦1 𝑢1 (𝑥 ) + 𝑦2 𝑢2 (𝑥 ) = 𝑥𝑢1 (𝑥 ) + 𝑥 5⁄2 𝑢2 (𝑥 ) (6)

- Hallemos 𝑢1 (𝑥) y 𝑢2 (𝑥 ):
 Lorena Maritza Terrios Guzmán 66
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

𝑥 𝑥 5⁄ 2 3
5⁄ 2 ) 5 3⁄ 2 | = 𝑥 5⁄2
𝑊 (𝑥 ) = 𝑊(𝑥, 𝑥 =|
1 𝑥 2
2

0 𝑥 5⁄2
𝑥 0
𝑊1 (𝑥 ) = | 5 3⁄2| = −𝑥 5⁄2 𝑙𝑛𝑥, 𝑊2 (𝑥 ) = | | = 𝑥𝑙𝑛𝑥
𝑙𝑛𝑥 𝑥 1 𝑙𝑛𝑥
2

𝑊1 (𝑥 ) −𝑥 5⁄2 𝑙𝑛𝑥 2 2 2𝑥
𝑢1 (𝑥 ) = ∫ 𝑑𝑥 = ∫ 𝑑𝑥 = − ∫ 𝑙𝑛𝑥𝑑𝑥 = − (𝑥𝑙𝑛𝑥 − 𝑥 ) = (1 − 𝑙𝑛𝑥) (7)
𝑊 (𝑥 ) 3 5⁄ 2 3 3 3
2𝑥
𝑊2 (𝑥 ) 𝑥𝑙𝑛𝑥 2 2
𝑢2 (𝑥 ) = ∫ 𝑑𝑥 = ∫ 𝑑𝑥 = ∫ 𝑥 −3⁄2𝑙𝑛𝑥𝑑𝑥 = − (2 + 𝑙𝑛𝑥 ) (8)
𝑊 (𝑥 ) 3 5⁄ 2 3 √ 𝑥
2𝑥
Sustituyendo (7) y (8) en (6) tenemos que
2𝑥 −2 2𝑥 2
𝑦𝑝 = 𝑥 [ (1 − 𝑙𝑛𝑥)] + 𝑥 5⁄2 [ (2 + 𝑙𝑛𝑥 )] = (7 + 2𝑙𝑛𝑥 )
3 √𝑥 3
PASO 5: La solución general de (1) es

5⁄ 2
2𝑥 2
𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1 𝑥 + 𝑐2𝑥 + (7 + 2𝑙𝑛𝑥 )
3
EJERCICIO 9
Obtener la solución general de cada ecuación de Cauchy –Euler.
a) 𝑥 2 𝑦 ′′ − 3𝑥𝑦 ′ + 3𝑦 = 2𝑥 4 𝑒 𝑥 b) 𝑥 3 𝑦 ′′′ + 5𝑥 2 𝑦 ′′ + 7𝑥𝑦 ′ + 8𝑦 = 0

TAREA 6: Obtenga la solución general de cada ecuación diferencial.

a) 𝑥 2 𝑦 ′′ − 2𝑥𝑦 ′ = 𝑥 3 𝑙𝑛𝑥 b) 𝑥 2 𝑦 ′′ + 𝑥𝑦´ − 9𝑦 = 𝑥 1⁄2 + 𝑥 −1⁄2

 Lorena Maritza Terrios Guzmán 67
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

TALLER 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

1. Determine si cada conjunto de funciones es linealmente independiente

1) 𝑓1 (𝑥 ) = 𝑥, 𝑓2 (𝑥 ) = 𝑥 2 , 𝑓3 (𝑥 ) = 4𝑥 − 3𝑥 2 2) 𝑓1 (𝑥 ) = 0, 𝑓2 (𝑥 ) = 𝑥, 𝑓3 (𝑥 ) = 𝑒 𝑥

3) 𝑓1 (𝑥 ) = 5, 𝑓2 (𝑥 ) = 𝑐𝑜𝑠 2 𝑥, 𝑓3 (𝑥 ) = 𝑠𝑒𝑛2 𝑥 4)𝑓1 (𝑥 ) = 𝑥, 𝑓2 (𝑥 ) = 𝑥 + 1, 𝑓3 (𝑥 ) = 𝑥 2

5) 𝑓1 (𝑥 ) = 𝑐𝑜𝑠2𝑥, 𝑓2 (𝑥 ) = 1, 𝑓3 (𝑥 ) = 𝑐𝑜𝑠 2 𝑥 6) 𝑓1 (𝑥 ) = 𝑒 𝑥 , 𝑓2 (𝑥 ) = 𝑒 −𝑥 , 𝑓3 (𝑥 ) = 𝑠𝑒𝑛ℎ𝑥

2. Compruebe que las funciones dadas forman un conjunto fundamental de soluciones de la

ecuación diferencial dada.

1) 𝑦1 = 𝑒 −3𝑥 , 𝑦2 = 𝑒 4𝑥 ; 𝑦 ′′ − 𝑦 ′ − 12𝑦 = 0
2) 𝑦1 = 𝑥 3 , 𝑦2 = 𝑥 4 ; 𝑥 2 𝑦 ′′ − 6𝑥𝑦 ′ + 12𝑦 = 0
3) 𝑦1 = 𝑥, 𝑦2 = 𝑥 −2 , 𝑦3 = 𝑥 −2 𝑙𝑛𝑥; 𝑥 3 𝑦 ′′′ + 6𝑥 2 𝑦 ′′ + 4𝑥𝑦 ′ − 4𝑦 = 0

3. Muestre que la función dada es la solución general de la ecuación diferencial indicada.

1) 𝑦 = 𝐴𝑐𝑜𝑠𝑥 + 𝐵𝑠𝑒𝑛𝑥 + 𝑥𝑠𝑒𝑛𝑥 + 𝑐𝑜𝑠𝑥𝑙 𝑛(𝑐𝑜𝑠𝑥 ) , 𝑦 ′′ + 𝑦 = 𝑠𝑒𝑐𝑥

1 1
2) 𝑦 = 𝐴𝑥 −1⁄2 + 𝐵𝑥 −1 + 𝑥 2 − 𝑥, 2𝑥 2 𝑦 ′′ + 5𝑥𝑦 ′ + 𝑦 = 𝑥 2 − 𝑥
15 6

4. Encontrar la solución general de cada ecuación diferencial lineal homogénea.

1) 4𝑦 ′′ + 𝑦 ′ = 0 2) 𝑦 ′′ − 36𝑦 = 0
3) 𝑦 ′′ + 8𝑦 ′ + 16𝑦 = 0 4) 12𝑦 ′′ − 5𝑦′ − 2𝑦 = 0
5) 𝑦 ′′ + 9𝑦 = 0 6) 𝑦 ′′′ − 4𝑦 ′′ − 5𝑦 ′ = 0
7) 𝑦 ′′′ − 𝑦 = 0 8) 𝑦 ′′′ − 5𝑦 ′′ + 3𝑦 ′ + 9𝑦 = 0
9) 𝑦 ′′′ + 𝑦 ′′ − 2𝑦 = 0 10) 𝑦 (4) + 𝑦 ′′′ + 𝑦 ′′ = 0
11) 𝑦 (5) + 5𝑦 (4) − 2𝑦 ′′′ − 10𝑦 ′′ + 𝑦 ′ + 5𝑦 = 0 12) 2𝑦 (5) − 7𝑦 (4) + 12𝑦 ′′′ + 8𝑦 ′′ = 0

5. Obtener la solución de cada problema de valor frontera.

1) 𝑦 ′′ − 10𝑦 ′ + 25𝑦 = 0, 𝑦(0) = 1, 𝑦(1) = 0
𝜋
2) 𝑦 ′′ + 𝑦 = 0, 𝑦(0) = 1, 𝑦 ′ (3 ) = −1

3) 𝑦 ′′ − 2𝑦 ′ + 𝑦 = 0, 𝑦(0) = 1, 𝑦′(𝑙𝑛2) = 1
 Lorena Maritza Terrios Guzmán 68
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

6. Aplicar el método de los coeficientes indeterminados para resolver cada ecuación

diferencial.
1
1) 𝑦 ′′ + 3𝑦 ′ + 2𝑦 = 6 2) 4 𝑦 ′′ + 𝑦 ′ + 𝑦 = 𝑥 2 − 2𝑥

3) 𝑦 ′′ − 8𝑦 ′ + 20𝑦 = 100𝑥 2 − 26𝑥𝑒 𝑥 4) 𝑦 ′′ + 3𝑦 = −48𝑥 2 𝑒 3𝑥

5) 4𝑦 ′′ − 4𝑦 ′ − 3𝑦 = 𝑐𝑜𝑠2𝑥 6) 𝑦 ′′ + 2𝑦 ′ = 2𝑥 + 5 − 𝑒 2𝑥

7) 𝑦 ′′ + 2𝑦 ′ − 24𝑦 = 16 − (𝑥 + 2)𝑒 4𝑥 8) 𝑦 ′′′ − 𝑦 ′′ − 4𝑦 ′ + 4𝑦 = 5 − 𝑒 𝑥 + 𝑒 2𝑥

9) 𝑦 (4) + 2𝑦 ′′ + 𝑦 = (𝑥 − 1)2 10) 𝑦 (4) − 𝑦 ′′ = 4𝑥 + 2𝑥𝑒 −𝑥

11) 𝑦 ′′′ − 2𝑦 ′′ − 4𝑦 ′ + 8𝑦 = 6𝑥𝑒 2𝑥 12) 𝑦 ′′ + 2𝑦 ′ + 𝑦 = 𝑠𝑒𝑛𝑥 + 3𝑐𝑜𝑠2𝑥

13) 𝑦 (4) + 𝑦 ′′′ = 1 − 𝑥 2 𝑒 −𝑥

7. Resolver cada ecuación diferencial usando el método del anulador.

1) 𝑦 ′′ − 9𝑦 = 54 2) 2𝑦 ′′ − 7𝑦 ′ + 5𝑦 = −29

3) 𝑦 ′′′ + 2𝑦 ′′ + 𝑦 ′ = 10 4) 𝑦 ′′′ + 𝑦 ′′ = 8𝑥 2

5) 𝑦 ′′ − 2𝑦 ′ + 2𝑦 = 𝑥 3 + 4𝑥 6) 𝑦 ′′ + 6𝑦 ′ + 8𝑦 = 2𝑥 + 3𝑒 −2𝑥

7) 𝑦 ′′ − 𝑦 = 5 + 𝑥 2 𝑒 𝑥 8) 𝑦 ′′ − 2𝑦 ′ + 5𝑦 = 𝑒 𝑥 𝑠𝑒𝑛𝑥

1
9) 𝑦 ′′ + 3𝑦 ′ − 10𝑦 = 𝑥 (𝑒 𝑥 + 1) 10) 𝑦 ′′ + 𝑦 ′ + 4 𝑦 = 𝑒 𝑥 (𝑠𝑒𝑛3𝑥 − 𝑐𝑜𝑠3𝑥)

11) 𝑦 ′′ + 4𝑦 = −8 + 4𝑐𝑜𝑠𝑥 + 3𝑠𝑒𝑛𝑥 12) 𝑦 ′′ − 2𝑦 ′ − 3𝑦 = 4𝑒 𝑥 − 9

13) 𝑦 ′′′ − 4𝑦 ′′ + 4𝑦 ′ = 5𝑥 2 − 6𝑥 + 4𝑥 2 𝑒 2𝑥 + 3𝑒 5𝑥

8. Aplicar el método de variación de parámetros para resolver cada ecuación diferencial.

1) 𝑦 ′′ + 𝑦 = 𝑠𝑒𝑐 𝑥 2) 𝑦 ′′ + 𝑦 = 𝑡𝑎𝑛 𝑥

3) 𝑦 ′′ + 𝑦 = 𝑠𝑒𝑛 𝑥 4) 𝑦 ′′ + 𝑦 = 𝑠𝑒𝑐 𝑥𝑡𝑎𝑛 𝑥

5) 𝑦 ′′ + 𝑦 = 𝑐𝑜𝑠 2 𝑥 6) 𝑦 ′′ + 𝑦 = 𝑠𝑒𝑐 2 𝑥

𝑒 2𝑥
7) 𝑦 ′′ − 𝑦 = 𝑐𝑜𝑠ℎ 𝑥 8) 𝑦 ′′ − 4𝑦 =
𝑥
 Lorena Maritza Terrios Guzmán 69
CAPÍTULO 4: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

9𝑥 1
9) 𝑦 ′′ − 9𝑦 = 10) 𝑦 ′′ + 3𝑦 ′ + 2𝑦 =
𝑒 3𝑥 1+𝑒 𝑥

𝑒𝑥
11) 𝑦 ′′ − 2𝑦 ′ + 𝑦 = 12) 𝑦 ′′ + 3𝑦 ′ + 2𝑦 = 𝑠𝑒𝑛𝑒 𝑥
1+𝑥 2

13) 𝑦 ′′ − 2𝑦 ′ + 𝑦 = 𝑒 𝑥 𝑎𝑟𝑐𝑡𝑎𝑛 𝑥 14) 2𝑦 ′′ + 2𝑦′ + 𝑦 = 4√𝑥

15) 𝑦 ′′ + 2𝑦′ + 𝑦 = 𝑒 −𝑥 𝑙𝑛 𝑥 16) 3𝑦 ′′ − 6𝑦′ + 6𝑦 = 𝑒 𝑥 𝑠𝑒𝑐 𝑥

17) 4𝑦′′ − 4𝑦′ + 𝑦 = 𝑒𝑥⁄2 √1 − 𝑥2 18) 𝑦 ′′ − 𝑦 ′ − 2𝑦 = 𝑠𝑒𝑛 2𝑥

9. Obtener la solución general de cada ecuación de Cauchy-Euler.

1) 𝑥 2 𝑦 ′′ + 3𝑥𝑦 ′ = 0 2) 𝑥 2 𝑦 ′′ − 5𝑥𝑦 ′ + 8𝑦 = 0

3) 𝑥𝑦 ′′ + 𝑦 ′ = 𝑥 4) 𝑥 2 𝑦 ′′ − 5𝑥𝑦 ′ + 8𝑦 = 8𝑥 6

5) 𝑥 2 𝑦 ′′ − 4𝑥𝑦 ′ + 6𝑦 = 𝑙𝑛𝑥 2 6) 𝑥 2 𝑦 ′′ + 10𝑥𝑦 ′ + 8𝑦 = 𝑥 2

7) 𝑥 3 𝑦 ′′′ + 5𝑥 2 𝑦 ′′ + 7𝑥𝑦 ′ + 8𝑦 = 0 8) 𝑥 2 𝑦 ′′ + 𝑥𝑦´ − 9𝑦 = 𝑥 1⁄2 + 𝑥 −1⁄2