4.

Ecuaciones diferenciales lineales de orden superior

4.1 Definición y teorema de existencia

Un problema de valor inicial para una ecuación diferencial lineal de orden 𝑛

𝑑𝑛𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
Resolver: 𝑎𝑛 (𝑥) 𝑑𝑥 𝑛 + 𝑎𝑛−1 (𝑥) 𝑑𝑥 𝑛−1 + . . . +𝑎1 (𝑥) 𝑑𝑥 + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥)

(1)
(𝑛−1)
Sujeta a: 𝑦(𝑥0 ) = 𝑦0 , 𝑦 ′ (𝑥0 ) = 𝑦0′ . . . 𝑦 (𝑛−1) (𝑥0 ) = 𝑦0
(𝑛−1)
en donde 𝑦0 , 𝑦0′ , . . . 𝑦0 , son constantes arbitrarias. Se busca una solución en algún intervalo
𝐼 que contenga al punto 𝑥 = 𝑥0 .

Teorema 4.1: Sean 𝑎𝑛 (𝑥), 𝑎𝑛−1 (𝑥),. . . , 𝑎1 (𝑥), 𝑎0 (𝑥) 𝑦 𝑔(𝑥) continuas en un intervalo 𝐼y sea
𝑎𝑛 (𝑥) ≠ 0 para todo x en este intervalo. Si 𝑥 = 𝑥0 es cualquier punto de este intervalo, entonces
existe una solución 𝑦(𝑥) del problema del valor inicial (1) en el intervalo y esta solución es única.

4.2 Teoría preliminar

4.2.1 Problemas de valor en la frontera

Otro tipo de problemas consiste en resolver una ecuación diferencial de orden dos o de orden
mayor que dos en la cual la variable dependiente y (o sus derivadas) se especifica en dos puntos
diferentes. Un problema como
𝑑2 𝑦 𝑑𝑦
Resolver: 𝑎2 (𝑥) + 𝑎1 (𝑥) + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥)
𝑑𝑥 2 𝑑𝑥

Sujeta a: 𝑦(𝑎) = 𝑦0 , 𝑦(𝑏) = 𝑦1

se llama un problema de valores de frontera de dos puntos o, simplemente, un problema de

valores en la frontera.

Ejemplo 4.1: El problema resolver la ecuación diferencial

𝑥 2 𝑦 ′′ − 2𝑥𝑦 ′ + 2𝑦 = 6
Sujeta a:

𝑦(1) = 0, 𝑦(2) = 3

es un problema de valor en la frontera.

4.2.2 Dependencia e independencia lineal

Definición 4.1: Se dice que un conjunto de funciones 𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥) es linealmente
dependiente en un intervalo 𝐼 si existen constantes 𝑐1 , 𝑐2 , . . . , 𝑐𝑛 no todas nulas, tales que

𝑐1 𝑓1 (𝑥) + 𝑐2 𝑓2 (𝑥)+ . . . +𝑐𝑛 𝑓𝑛 (𝑥) = 0

Elaborado por Javier E. Herrera C. pág. 1

para toda 𝑥 en el intervalo.

Definición 4.2: Se dice que un conjunto de funciones 𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥) es linealmente
independiente en un intervalo 𝐼 si no es linealmente dependiente en el intervalo.

En otras palabras, un conjunto de funciones es linealmente independiente en un intervalo si las

únicas constantes para las cuales

𝑐1 𝑓1 (𝑥) + 𝑐2 𝑓2 (𝑥)+ . . . +𝑐𝑛 𝑓𝑛 (𝑥) = 0

Para toda x en el intervalo, son 𝑐1 = 𝑐2 = . . . = 𝑐𝑛 = 0.

Ejemplo 4.2: Las funciones 𝑓1 (𝑥) = 𝑥, 𝑓2 (𝑥) = 𝑥 2 𝑦 𝑓3 (𝑥) = 4𝑥 − 3𝑥 2 son linealmente

dependientes ya que

𝑐1 𝑥 + 𝑐2 𝑥 2 + 𝑐3 (4𝑥 − 3𝑥 2 ) = 0
cuando 𝑐1 = −4, 𝑐2 = 3, 𝑐3 = 1

Ejemplo 4.3: Las funciones 𝑓1 (𝑥) = 5, 𝑓2 (𝑥) = 𝑐𝑜𝑠 2 𝑥 𝑦 𝑓3 (𝑥) = 𝑠𝑒𝑛2 𝑥 son linealmente
dependientes ya que

𝑐1 5 + 𝑐2 𝑐𝑜𝑠 2 𝑥 + 𝑐3 𝑠𝑒𝑛2 𝑥 = 0
1
cuando 𝑐1 = − , 𝑐2 = 1, 𝑐3 = 1
5

4.2.3 El Wronskiano

Para determinar de forma mecánica si un conjunto de funciones es linealmente dependiente (ld) o

linealmente independiente (li), podemos utilizar el siguiente determinante

Definición 4.3: El Wronskiano

Suponga que cada una de las funciones 𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥) tienen al menos 𝑛 − 1
derivadas. El determinante

𝑊൫𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥)൯ =

las primas denotan derivadas se llama el Wronskiano de las funciones.

Elaborado por Javier E. Herrera C. pág. 2

Teorema 4.2: Si el determinante

𝑊൫𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥)൯ =

no es cero por lo menos en un punto del intervalo 𝐼, entonces 𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥) son
linealmente independientes en el intervalo.

Ejemplo 4.4: Las funciones 𝑓1 (𝑥) = 𝑒 3𝑥 𝑦 𝑓2 (𝑥) = 𝑒 −3𝑥 son linealmente independiente, ya que

Corolario 4.1: Si 𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥) tienen al menos 𝑛 − 1 derivadas y son linealmente
dependientes en 𝐼, entonces

𝑊൫𝑓1 (𝑥), 𝑓2 (𝑥), . . . , 𝑓𝑛 (𝑥)൯ = 0

para toda 𝑥 del intervalo.

4.3 Solución de ecuaciones diferenciales lineales

Una ecuación diferencial lineal de la forma

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

Se llama homogénea, en tanto que a

𝑑𝑛𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥) 𝑑𝑥 𝑛 + 𝑎𝑛−1 (𝑥) 𝑑𝑥 𝑛−1 + . . . +𝑎1 (𝑥) 𝑑𝑥 + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥) (3)

en donde 𝑔(𝑥) no es idénticamente nula, se llama no homogénea.

Ejemplo 4.5: La ecuación

𝑥 2 𝑦 ′′ − 2𝑥𝑦 ′ + 2𝑦 = 6
Es una ecuación diferencial lineal no homogénea de segundo orden y la ecuación

𝑥 3 𝑦 ′′′ + 6𝑥𝑦 ′′ − 4𝑦 = 0
es una ecuación diferencial lineal homogénea de tercer orden.

Principio de superposición

Elaborado por Javier E. Herrera C. pág. 3

Teorema 4.3: Sean 𝑦1 , 𝑦2 , . . . , 𝑦𝑘 soluciones de la ecuación diferencial lineal homogénea de
orden 𝑛 (2) en un intervalo 𝐼. Entonces la combinación lineal

𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + . . . +𝑐𝑘 𝑦𝑘
en donde los 𝑐𝑖 , 𝑖 = 1,2, … , 𝑘 son constantes arbitrarias, también es una solución en el intervalo.

Corolario 4.2:

a. Si 𝑦1 (𝑥) es una solución de una ecuación diferencial homogénea, entonces un múltiplo

constante de ella, 𝑦 = 𝑐1 𝑦1 (𝑥), también es una solución.

b. Una ecuación diferencial lineal homogénea siempre tiene la solución trivial 𝑦 = 0

Ejemplo 4.6: Las funciones 𝑦1 = 𝑒 𝑥 𝑦 𝑦2 = 𝑒 −𝑥 son soluciones de la ecuación diferencial

𝑦 ′′ − 𝑦 = 0
en el intervalo (−∞, ∞). Por el principio de superposición, la combinación lineal

𝑦 = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 −𝑥
también es una solución de la ecuación en el intervalo.

Teorema 4.4: Sean 𝑦1 , 𝑦2 , . . . , 𝑦𝑛 , 𝑛 soluciones de la ecuación (2) en un intervalo 𝐼. Entonces el

conjunto de soluciones es linealmente independiente si y sólo si

𝑊(𝑦1 , 𝑦2 , . . . , 𝑦𝑛 ) ≠ 0
para toda 𝑥 en el intervalo.

Definición 4.4: Se llama conjunto fundamental de soluciones en un intervalo 𝐼 a cualquier

conjunto 𝑦1 , 𝑦2 , . . . , 𝑦𝑛 de 𝑛 soluciones linealmente independientes de la ecuación diferencial
lineal homogénea (2) en el intervalo 𝐼.

Teorema 4.5: Sea 𝑦1 , 𝑦2 , . . . , 𝑦𝑛 un conjunto fundamental de soluciones de la ecuación (2). Si

𝑌(𝑥) es cualquier solución de (2) en 𝐼, entonces es posible encontrar constantes 𝐶1 , 𝐶2 , . . . , 𝐶𝑛
tales que

𝑌 = 𝐶1 𝑦1 (𝑥) + 𝐶2 𝑦2 (𝑥) + . . . +𝐶𝑛 𝑦𝑛 (𝑥)

Teorema 4.6: Existe un conjunto fundamental de soluciones para la ecuación diferencial (2) en un
intervalo 𝐼.

Definición 4.5: Sea 𝑦1 , 𝑦2 , . . . , 𝑦𝑛 un conjunto fundamental de soluciones de la ecuación (2) en

un intervalo 𝐼. Se define como solución general de la ecuación en el intervalo a

𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + . . . +𝑐𝑛 𝑦𝑛
en donde los 𝑐𝑖 𝑐𝑜𝑛 𝑖 = 1,2, . . . , 𝑛 son constantes arbitrarias.

Ejemplo 4.7: Las funciones 𝑦1 = 𝑒 𝑥 𝑦 𝑦2 = 𝑒 −𝑥 son soluciones de la ecuación diferencial

𝑦 ′′ − 𝑦 = 0

Elaborado por Javier E. Herrera C. pág. 4

en el intervalo (−∞, ∞). Observemos que estas soluciones son linealmente independientes en el
eje 𝑥, en efecto

para toda 𝑥, por lo tanto 𝑦1 , 𝑦 𝑦2 forman un conjunto fundamental de soluciones con lo cual

𝑦 = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 −𝑥
es la solución general de la ecuación en el intervalo.

Ejemplo 4.8: Compruebe que las funciones 𝑦1 = 𝑥, 𝑦2 = 𝑥 −2 𝑦 𝑦3 = 𝑥 −2 𝑙𝑛𝑥 forman un conjunto

fundamental de soluciones de la ecuación diferencial

𝑥 3 𝑦 ′′′ + 6𝑥 2 𝑦 ′′ + 4𝑥𝑦 ′ − 4𝑦 = 0
en el intervalo (0, ∞). Forme la solución general.

Solución: Estas funciones satisfacen la ecuación diferencial (compruébelo) y además son

linealmente independientes ya que

𝑊(𝑥, 𝑥 −2 , 𝑥 −2 𝑙𝑛𝑥) = 9𝑥 −6 ≠ 0
para (0, ∞). La solución general en este intervalo es

𝑦 = 𝑐1 𝑥 + 𝑐2 𝑥 −2 + 𝑐3 𝑥 −2 𝑙𝑛𝑥
Cualquier función 𝑦𝑝 libre de parámetros arbitrarios, que satisfaga la ecuación (3) se dice que es
una solución particular o integral particular de la ecuación.

Teorema 4.7: Sea𝑦𝑝 cualquier solución particular de la ecuación diferencial lineal no homogénea
de n-ésimo orden (3) en un intervalo 𝐼, y sea 𝑦1 , 𝑦2 , . . . , 𝑦𝑛 un conjunto fundamental de
soluciones de la ecuación diferencial homogénea asociada (2) en 𝐼. Entonces la solución general de
la ecuación en el intervalo es

𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + . . . +𝑐𝑛 𝑦𝑛 + 𝑦𝑝 ,

donde las 𝑐𝑖 , 𝑖 = 1,2, … , 𝑛 son constantes arbitrarias.

En este teorema vemos que la solución general de una ecuación lineal no homogénea está
compuesta por la suma de dos funciones

𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + . . . +𝑐𝑛 𝑦𝑛 + 𝑦𝑝 = 𝑦𝑐 + 𝑦𝑝

La combinación lineal 𝑦𝑐 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + . . . +𝑐𝑛 𝑦𝑛 , que es la solución general de la ecuación

(2), se llama función complementaria para la ecuación (3), entonces quiere decir que para resolver
una ecuación diferencial lineal no homogénea, debemos resolver primero la ecuación homogénea
asociada y luego se encuentra una solución particular de la ecuación no homogénea, por lo que la
solución general de la ecuación no homogénea es entonces

𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑓𝑢𝑛𝑐𝑖ó𝑛 𝑐𝑜𝑚𝑝𝑙𝑒𝑚𝑒𝑛𝑡𝑎𝑟𝑖𝑎 + 𝑐𝑢𝑎𝑙𝑞𝑢𝑖𝑒𝑟 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑝𝑎𝑟𝑡𝑖𝑐𝑢𝑙𝑎𝑟

Elaborado por Javier E. Herrera C. pág. 5

4.2.4 Reducción de orden

En la sección 2, estudiamos algunos casos de reducción de orden, en esta se estudia la reducción

de orden de una ecuación diferencial lineal de segundo orden., el método consiste en que dada
una solución 𝑦1 (𝑥) de la ecuación diferencial lineal de segundo orden
𝑑2 𝑦 𝑑𝑦
𝑑𝑥 2
+ 𝑝(𝑥) 𝑑𝑥 + 𝑞(𝑥)𝑦 = 0 (4)

puede determinarse una segunda solución 𝑦2 (𝑥) que sea linealmente independiente con 𝑦1 (𝑥),
como lo describe el siguiente ejemplo.

Ejemplo 4.9: Dado que 𝑦1 (𝑥) = 𝑥𝑒 −𝑥 es una solución de la ecuación diferencial

𝑦 ′′ + 2𝑦 ′ + 𝑦 = 0
encuentre una segunda solución.

Solución: Sea 𝑦 = 𝑢(𝑥)𝑦1 (𝑥) = 𝑢(𝑥)𝑒 −𝑥 , por lo tanto,

𝑦 ′ = (1 − 𝑥)𝑒 −𝑥 + 𝑥𝑒 −𝑥 𝑢′ , 𝑦 ′′ = 𝑥𝑒 −𝑥 𝑢′′ + 2(1 − 𝑥)𝑒 −𝑥 𝑢′ − (2 − 𝑥)𝑒 −𝑥 𝑢

sustituyendo en la ecuación dada tenemos:
2
𝑢′′ + 𝑢′ = 0
𝑥
Si 𝑢′ = 𝑝, obtenemos la ecuación lineal de primer orden
2
𝑝′ + 𝑝 = 0
𝑥
Cuya solución es
𝑐
𝑝=
𝑥2
y como

𝑢′ = 𝑝
tenemos que,
𝑐1
𝑢=
𝑥
Tomando 𝑐1 = 1, una segunda es
1
𝑦2 = 𝑥𝑒 −𝑥 = 𝑒 −𝑥
𝑥
Este procedimiento se reduce utilizando el siguiente teorema:

Teorema 4.7: Si 𝑦1 (𝑥) es una solución de la ecuación diferencial

Elaborado por Javier E. Herrera C. pág. 6

 𝑑2 𝑦 𝑑𝑦
2
+ 𝑝(𝑥) + 𝑞(𝑥)𝑦 = 0
𝑑𝑥 𝑑𝑥
entonces una segunda solución 𝑦2 (𝑥) de (4) linealmente independiente con 𝑦1 (𝑥) es
𝑒 − ∫ 𝑝(𝑥)𝑑𝑥
𝑦2 (𝑥) = 𝑦1 (𝑥) ∫ 𝑑𝑥 (5)
𝑦12 (𝑥)

Si aplicamos la ecuación (5), al ejemplo 4.9, tenemos lo siguiente

𝑒 − ∫ 2𝑑𝑥 𝑒 −2𝑥
𝑦2 (𝑥) = 𝑥𝑒 −𝑥 ∫ 𝑑𝑥 = 𝑥𝑒 −𝑥
∫ 𝑑𝑥
(𝑥𝑒 −𝑥 )2 (𝑥𝑒 −𝑥 )2
= −𝑒 −𝑥
4.3 Ecuaciones lineales homogéneas con coeficientes constantes

La ecuación lineal de primer orden

𝑦 ′ + 𝑎𝑦 = 0,
𝑑𝑦
= −𝑎𝑦
𝑑𝑥
𝑑𝑦
= −𝑎𝑑𝑥
𝑦
𝑙𝑛𝑦 = −𝑎𝑥 + 𝑐
𝑦 = 𝑐1 𝑒 −𝑎𝑥
en donde 𝑎 es una constante tiene como solución como solución la función exponencial

𝑦 = 𝑐1 𝑒 −𝑎𝑥 𝑒𝑛 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜 (−∞, ∞).

Por lo que, al tratar de determinar la solución de la ecuación de orden superior

𝑑𝑛𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 𝑑𝑥 𝑛 + 𝑎𝑛−1 𝑑𝑥 𝑛−1 + . . . +𝑎1 𝑑𝑥 + 𝑎0 𝑦 = 0

en donde 𝑎𝑖 , 𝑖 = 1,2, . . . , 𝑛 son constantes, en el intervalo (−∞, ∞), sería natural pensar en
que si existen soluciones exponenciales en el intervalo.

La respuesta es que sí son funciones exponenciales sus soluciones o se construyen a partir de

funciones exponenciales.

4.3.1 Ecuación auxiliar. Raíces reales distintas, raíces reales iguales y raíces complejas
conjugadas

Comenzaremos con el caso particular de la ecuación de segundo orden

𝑎𝑦 ′′ + 𝑏𝑦 ′ + 𝑐𝑦 = 0 (6)

Supongamos que una solución tiene la forma

𝑦 = 𝑒 𝑚𝑥

Elaborado por Javier E. Herrera C. pág. 7

Luego,

𝑦 ′ = 𝑚𝑒 𝑚𝑥 , 𝑦 ′′ = 𝑚2 𝑒 𝑚𝑥
Ahora sustituyendo en (6) se tiene: 𝑎𝑦 ′′ + 𝑏𝑦 ′ + 𝑐𝑦 = 0

(𝑎𝑚2 + 𝑏𝑚 + 𝑐)𝑒 𝑚𝑥 = 0; (7)

y puesto que 𝑒 𝑚𝑥 ≠ 0, se concluye que (7) se cumple si y sólo si 𝑚 satisface la ecuación

(𝑎𝑚2 + 𝑏𝑚 + 𝑐) = 0 (8)

Esta última ecuación se le llama ecuación auxiliar o ecuación característica de la ecuación

diferencial (6), observemos que la ecuación (8) es una ecuación cuadrática, por lo que al resolverla
debemos considerar tres casos, ya que las raíces deben ser

−𝑏 + √𝑏 2 − 4𝑎𝑐 −𝑏 − √𝑏 2 − 4𝑎𝑐
𝑚1 = , 𝑚2 =
2𝑎 2𝑎
Caso 1:

Si 𝒃𝟐 − 𝟒𝒂𝒄 > 𝟎: Bajo esta suposición la ecuación auxiliar tiene dos raíces reales distintas
𝑚1 𝑦 𝑚2 , así encontramos dos soluciones de la ecuación (6) 𝑦1 = 𝑒 𝑚1 𝑥 y 𝑦2 = 𝑒 𝑚2 𝑥 , las cuales
son linealmente independientes en (−∞, ∞), por lo tanto, forman un conjunto fundamental de
soluciones y como consecuencia la solución general de (6) es:

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

Caso 2: 𝒃𝟐 − 𝟒𝒂𝒄 = 𝟎: En este, se tiene que 𝑚1 = 𝑚2 , y por lo tanto se obtiene sólo una solución
exponencial 𝑦1 = 𝑒 𝑚1 𝑥 de (6), pero una segunda solución de la ecuación (6) viene dada por

𝑒 2𝑚1 𝑥
𝑦2 = 𝑒 𝑚1 𝑥 ∫ 𝑑𝑥 = 𝑒 𝑚1 𝑥 ∫ 𝑑𝑥 = 𝑥𝑒 𝑚1 𝑥
𝑒 2𝑚1 𝑥
entonces la solución general es

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

Caso 3: 𝒃𝟐 − 𝟒𝒂𝒄 < 𝟎: Con esta suposición se obtiene que 𝑚1 𝑦 𝑚2 son raíces complejas
(números complejos conjugados), entonces se puede escribir 𝑚1 = 𝛼 + 𝛽𝑖 y 𝑚2 = 𝛼 − 𝛽𝑖, donde
𝛼 𝑦 𝛽 > 0 son reales, de manera formal no hay diferencia en este caso y el caso 1, por lo que

𝑦 = 𝑐1 𝑒 (𝛼+𝛽𝑖 )𝑥 + 𝑐2 𝑒 (𝛼−𝛽𝑖) 𝑥

En donde 𝑦1 = 𝑒 (𝛼+𝛽𝑖 )𝑥 y 𝑦2 = 𝑒 (𝛼−𝛽𝑖) 𝑥 son dos soluciones linealmente independientes.

Por conveniencia se prefiere trabajar con funciones reales en lugar de exponenciales, entonces
gracias a la fórmula de Euler: 𝑒 𝑖𝜃 = 𝑐𝑜𝑠𝜃 + 𝑖𝑠𝑒𝑛𝜃 donde 𝜃 es cualquier número real, se tiene que

𝑒 (𝛼+𝛽𝑖)𝑥 = 𝑒 𝛼𝑥 (𝑐𝑜𝑠𝛽𝑥 + 𝑖𝑠𝑒𝑛𝛽𝑥) (9)

𝑒 (𝛼−𝛽𝑖)𝑥 = 𝑒 𝛼𝑥 (𝑐𝑜𝑠𝛽𝑥 − 𝑖𝑠𝑒𝑛𝛽𝑥) (10)

Elaborado por Javier E. Herrera C. pág. 8

Y como estamos interesados sólo en soluciones que sean funciones con valores reales, podemos
sumar (9) y (10)y dividir por 2, o restar y dividir por 2i, obteniendo así

𝑒 𝛼𝑥 𝑐𝑜𝑠𝛽𝑥 𝑦 𝑒 𝛼𝑥 𝑠𝑒𝑛𝛽𝑥 𝑒 𝑖𝜋 = −1, 𝑒 𝑖𝜋 + 1 = 0

Estas son soluciones linealmente independientes, de modo que la solución general en este caso es

𝑦 = 𝑒 𝛼𝑥 (𝑐1 𝑐𝑜𝑠𝛽𝑥 + 𝑐2 𝑠𝑒𝑛𝛽𝑥)

Ejemplo 4.10: Resolver 4𝑦 ′′ + 𝑦 ′ = 0

Solución: La ecuación auxiliar de esta ecuación diferencial es

4𝑚2 + 𝑚 = 0
Con lo cual,

4𝑚2 + 𝑚 = 𝑚(4𝑚 + 1) = 0, 𝑚 = 0 𝑦 4𝑚 + 1 = 0
1
1
Así, 𝑚1 = 0 y 𝑚2 = − 4, por lo tanto, 𝑦1 = 𝑒 0𝑥 = 1, 𝑦2 = 𝑒 −4𝑥 ,la solución general es
1
𝑦 = 𝑐1 + 𝑐2 𝑒 −4𝑥
𝑑2 𝑦 𝑑𝑦
Ejemplo 4.11: Resolver 𝑑𝑡 2 − 4 𝑑𝑡 − 5𝑦 = 0, 𝑦(1) = 0, 𝑦 ′ (1) = 2

Solución: Como la ecuación auxiliar es

𝑚2 − 4𝑚 − 5 = 0
Entonces, 𝑚 = −1 y 𝑚 = 5, por lo que

𝑦 = 𝑐1 𝑒 −𝑡 + 𝑐2 𝑒 5𝑡
Ahora de 𝑦(1) = 0, 𝑦 ′ (1) = 2, se tiene

𝑐1 𝑒 −1 + 𝑐2 𝑒 5 = 0

−𝑐1 𝑒 −1 + 5𝑐2 𝑒 5 = 2
𝑒 𝑒 −5
Con lo cual, 𝑐1 = − y 𝑐2 = , luego solución particular es
3 3

1 1−𝑡 1 5𝑡−5
𝑦=− 𝑒 + 𝑒
3 3
Ejemplo 4.12: Resuelva la ecuación diferencial 𝑦 ′′ + 8𝑦 ′ + 16𝑦 = 0

Solución: Como la ecuación auxiliar de esta ecuación diferencial es

𝑚2 + 8𝑚 + 16 = 0

Entonces, 𝑚1 = 𝑚2 = −4, por lo tanto, 𝑦 = 𝑐1 𝑒 −4𝑥 + 𝑐2 𝑥𝑒 −4𝑥

Ejemplo 4.13: Resolver 𝑦 ′′ − 4𝑦 ′ + 5𝑦 = 0

Elaborado por Javier E. Herrera C. pág. 9

Solución: La ecuación auxiliar de esta ecuación es

𝑚2 − 4𝑚 + 5 = 0
Entonces, 𝑚1 = −2 + 𝑖 y 𝑚2 = −2 − 𝑖; por lo tanto,
𝑦 = 𝑒 𝛼𝑥 (𝑐1 𝑐𝑜𝑠𝛽𝑥 + 𝑐2 𝑠𝑒𝑛𝛽𝑥)

𝑦 = 𝑒 −2𝑥 (𝑐1 𝑐𝑜𝑠𝑥 + 𝑐2 𝑠𝑒𝑛𝑥)

En general, para resolver una ecuación diferencial de 𝑛 − é𝑠𝑖𝑚𝑜 orden
𝑑𝑛𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 𝑑𝑥 𝑛 + 𝑎𝑛−1 𝑑𝑥 𝑛−1 + . . . +𝑎1 𝑑𝑥 + 𝑎0 𝑦 = 0 (2)

Se debe resolver una ecuación polinomial de 𝑛 − é𝑠𝑖𝑚𝑜 grado

𝑎𝑛 𝑚𝑛 + 𝑎𝑛−1 𝑚𝑛−1 + . . . +𝑎2 𝑚2 + 𝑎1 𝑚 + 𝑎0 = 0 (11)

En el caso de que todas las raíces sean reales y distintas la solución de la ecuación (2) es

𝑦 = 𝑐1 𝑒 𝑚1 𝑥 + 𝑐2 𝑒 𝑚2 𝑥 + . . . +𝑐𝑛 𝑒 𝑚𝑛 𝑥 .

Resumir los otros dos casos no es tan sencillo que el primer caso pues surgen diferentes
combinaciones, pero cuando una raíza tiene una multiplicidad 𝑘, por ejemplo; 𝑚1 entonces se
puede demostrar que las soluciones linealmente independientes son

𝑒 𝑚1 𝑥 , 𝑥𝑒 𝑚1 𝑥 , 𝑥 2 𝑒 𝑚1 𝑥 , . . . , 𝑥 𝑘−1 𝑒 𝑚1 𝑥
Y la solución general es la combinación lineal

𝑐1 𝑒 𝑚1 𝑥 + 𝑐2 𝑥𝑒 𝑚1𝑥 + 𝑐3 𝑥 2 𝑒 𝑚1 𝑥 + . . . +𝑐𝑘 𝑥 𝑘−1 𝑒 𝑚1 𝑥

Ejemplo 4.14: Resuelva 𝑦 ′′′ − 4𝑦 ′′ − 5𝑦 ′ = 0

Solución: La ecuación auxiliar es 𝑚3 − 4𝑚2 − 5𝑚 = 0 y sus raíces son

𝑚1 = 0, 𝑚2 = 5, 𝑦 𝑚3 = −1
Por lo que,

𝑦 = 𝑐1 + 𝑐2 𝑒 5𝑥 + 𝑐3 𝑒 −𝑥
Ejemplo 4.15: Resuelva 𝑦 ′′′ + 𝑦 ′′ − 2𝑦 = 0

Solución: La ecuación auxiliar es 𝑚3 + 𝑚2 − 2 = 0 y sus raíces son 𝑚1 = 1 , 𝑚2 = −1 + 𝑖

y 𝑚3 = −1 − 𝑖 por lo que,

𝑦 = 𝑐1 𝑒 𝑡 + 𝑒 −𝑡 (𝑐2 𝑐𝑜𝑠𝑡 + 𝑐3 𝑠𝑒𝑛𝑡)

𝑑5 𝑦 𝑑4 𝑦 𝑑3 𝑦 𝑑2 𝑦 𝑑𝑦
Ejemplo 4.16: Resuelva 𝑑𝑥 5 + 5 𝑑𝑥 4 − 2 𝑑𝑥 3 − 10 𝑑𝑥 2 + 𝑑𝑥 + 5𝑦 = 0

Solución: La ecuación auxiliar relacionada a esta ecuación diferencial es

𝑚5 + 5𝑚4 − 2𝑚3 − 10𝑚2 + 𝑚 + 5 = 0

Elaborado por Javier E. Herrera C. pág.

10
y sus raíces son 𝑚1 = −1, 𝑚2 = −1, 𝑚3 = 1, 𝑚4 = 1 𝑦 𝑚5 = −5, entonces la solución general
de la ecuación diferencial es:

𝑦 = 𝑐1 𝑒 −𝑥 + 𝑐2 𝑥𝑒 −𝑥 + 𝑐3 𝑒 𝑥 + 𝑐4 𝑥𝑒 𝑥 + 𝑐5 𝑒 −5𝑥

4.4 Ecuaciones no homogéneas con coeficientes constantes

4.4.1 Coeficientes indeterminados
El método de los coeficientes indeterminados se aplica a la ecuación
𝑑𝑛𝑦 𝑑 𝑛−1 𝑦 𝑑𝑦
𝑎𝑛 (𝑥) 𝑑𝑥 𝑛 + 𝑎𝑛−1 (𝑥) 𝑑𝑥 𝑛−1 + . . . +𝑎1 (𝑥) 𝑑𝑥 + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥) (3)

donde 𝒈(𝒙) contiene un polinomio, términos de la forma 𝒔𝒆𝒏𝒃𝒙, 𝒄𝒐𝒔𝒃𝒙, 𝒆𝒃𝒙 donde 𝑏 es una
constante o combinaciones de sumas y productos de estos.

Debido a que la combinación lineal de derivadas

(𝑛) (𝑛−1)
𝑎𝑛 𝑦𝑝 + 𝑎𝑛−1 𝑦𝑝 + . . . +𝑎1 𝑦𝑝′ + 𝑎0 𝑦𝑝

debe ser idéntica a 𝑔(𝑥), parece razonable suponer que 𝑦𝑝 tiene la misma forma que 𝑔(𝑥).

Los siguientes ejemplos ilustran el método.

1
Ejemplo 4.17: Obtenga la solución general de la ecuación diferencial 𝑦 ′′ + 𝑦 ′ + 𝑦 = 𝑥 2 − 2𝑥
4
1
Solución: Primero se resuelve la ecuación homogénea asociada 4 𝑦 ′′ + 𝑦 ′ + 𝑦 = 0, la ecuación
1
auxiliar es 4 𝑚2 + 𝑚 + 1 = 0 cuyas raíces son 𝑚1 = 𝑚2 = −2, entonces 𝑦𝑐 = 𝑐1 𝑒 −2𝑥 + 𝑐2 𝑥𝑒 −2𝑥 ,
ahora debido a que 𝑔(𝑥) es un polinomio cuadrático supongamos que una solución particular es
de la forma de un polinomio cuadrático, es decir, 𝑦𝑝 = 𝐴𝑥 2 + 𝐵𝑥 + 𝐶, luego 𝑦𝑝′ = 2𝐴𝑥 + 𝐵,
1
𝑦𝑝′′ = 2𝐴 sustituyendo en la ecuación diferencial 4 𝑦 ′′ + 𝑦 ′ + 𝑦 = 𝑥 2 − 2𝑥, obtenemos

1
(2𝐴) + 2𝐴𝑥 + 𝐵 + 𝐴𝑥 2 + 𝐵𝑥 + 𝐶 = 𝑥 2 − 2𝑥
4
Es decir,
1
𝐴𝑥 2 + (2𝐴 + 𝐵)𝑥 + 𝐴 + 𝐵 + 𝐶 = 𝑥 2 − 2𝑥
2
Como se supone que esta última ecuación es una identidad, los coeficientes de los términos
semejantes deben ser iguales, así
1
𝐴 = 1, 2𝐴 + 𝐵 = −2, 𝐴+𝐵+𝐶 =0
2
Por lo tanto,

2(1) + 𝐵 = −2, 𝑒𝑠 𝑑𝑒𝑐𝑖𝑟 𝐵 = −4

1 7
(1) − 4 + 𝐶 = 0, 𝑐𝑜𝑛 𝑙𝑜 𝑐𝑢𝑎𝑙 𝐶 =
2 2

Elaborado por Javier E. Herrera C. pág.

11
Así,

𝑦𝑝 = 𝐴𝑥 2 + 𝐵𝑥 + 𝐶
7
𝑦𝑝 = 𝑥 2 − 4𝑥 +
2
Por último, la solución general de la ecuación dada es
7
𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1 𝑒 −2𝑥 + 𝑐2 𝑥𝑒 −2𝑥 + 𝑥 2 − 4𝑥 +
2
En la siguiente tabla se presentan ejemplos específicos de 𝑔(𝑥) junto con una correspondiente
forma de la solución particular.

Ejemplo 4.18: Resuelva la siguiente ecuación diferencial 𝑦 ′′ + 3𝑦 ′ + 2𝑦 = 6

Solución: Primero resolvemos la ecuación homogénea asociada 𝑦 ′′ + 3𝑦 ′ + 2𝑦 = 0

Ecuación auxiliar es

𝑚2 + 3𝑚 + 2 = 0
Por lo tanto.

𝑚1 = −2 𝑚2 = −1
Con lo cual,

Elaborado por Javier E. Herrera C. pág.

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

La solución particular sugerida es

𝑦𝑝 = 𝐴, 𝑙𝑢𝑒𝑔𝑜 𝑦𝑝′ = 0, 𝑦𝑝′′ = 0

Sustituimos en 𝑦 ′′ + 3𝑦 ′ + 2𝑦 = 6, entonces 0 + 3(0) + 2𝐴 = 6

𝐴=3
Luego 𝑦𝑝 = 3

𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1 𝑒 −2𝑥 + 𝑐2 𝑒 −𝑥 + 3

Ejemplo 4.19: Resolver la ecuación 𝑦 ′′ − 10𝑦 ′ + 25𝑦 = 30𝑥 + 3

Primero resolvemos la ecuación homogénea

𝑦 ′′ − 10𝑦 ′ + 25𝑦 = 0
Ecuación auxiliar es

𝑚2 − 10𝑚 + 25 = 0
Luego,

𝑚1 = 𝑚2 = 5

Con lo cual

𝑦𝑐 = 𝑐1 𝑒 5𝑥 + 𝑐2 𝑥𝑒 5𝑥
Ahora la solución particular sugerida es

𝑦𝑝 = 𝐴𝑥 + 𝐵, 𝑦𝑝′ = 𝐴. 𝑦𝑝′′ = 0

Por lo tanto,

𝑦𝑝′ = 𝐴. 𝑦𝑝′′ = 0

Sustituyendo en 𝑦 ′′ − 10𝑦 ′ + 25𝑦 = 30𝑥 + 3, se tiene

0 − 10𝐴 + 25(𝐴𝑥 + 𝐵) = 30𝑥 + 3

−10𝐴 + 25𝐴𝑥 + 25𝐵 = 30𝑥 + 3
25𝐴𝑥 − 10𝐴 + 25𝐵 = 30𝑥 + 3
25𝐴 = 30

Elaborado por Javier E. Herrera C. pág.

13
 6
𝐴=
5
−10𝐴 + 25𝐵 = 3
3
𝐵=
5
6 3
𝑦𝑝 = 𝑥 +
5 5
Por lo que la solución general es
6 3
𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1 𝑒 5𝑥 + 𝑐2 𝑥𝑒 5𝑥 + 𝑥 +
5 5
Si alguna 𝑦𝑝𝑖 , contiene soluciones que duplican los términos de 𝑦𝑐 , se debe multiplicar por 𝑥 𝑛 ,
donde 𝑛 desde el entero positivo que elimina esa duplicación; Ilustremos esta situación con los
siguientes ejemplos

Ejemplo 4.20: Determine la forma de la solución particular de la ecuación

a. 𝑦 ′′ − 𝑦 ′ = −3.
𝑥
1
b. 𝑦 ′′ − 𝑦 ′ + 𝑦 = 3 + 𝑒 2
4

Solución:

a. La ecuación auxiliar es

𝑚2 − 𝑚 = 0,

por lo que

𝑚 = 0 𝑦 𝑚 = 1,

Entonces

𝑦1 = 1 y 𝑦2 = 𝑒 𝑥

Así

𝑦𝑐 = 𝑐1 𝑒 𝑥 + 𝑐2 .

Ahora, la 𝑦𝑝 sugerida sería 𝑦𝑝 = 𝐴,

Pero existe una duplicidad en 𝑦𝑐 = 𝑐1 𝑒 𝑥 + 𝑐2 , pero multiplicando 𝑦𝑝 = 𝐴 por 𝑥, se elimina la

duplicidad, es decir,

𝑦𝑝 = 𝐴𝑥
𝑥
1
b. La ecuación auxiliar de esta ecuación 𝑦 ′′ − 𝑦 ′ + 𝑦 = 3 + 𝑒 2 es
4

1
𝑚2 − 𝑚 + =0
4
Elaborado por Javier E. Herrera C. pág.
14
Por lo tanto,
1
𝑚1 = 𝑚2 = 2.

Con lo cual,
𝑥 𝑥
𝑦𝑐 = 𝑐1 𝑒 2 + 𝑐2 𝑥𝑒 2
1 𝑥
𝑦 ′′ − 𝑦 ′ + 𝑦 = 3 + 𝑒 2
4
Observado la tabla 4.1, se sugiere que
𝑥
𝑦𝑝 = 𝐴 + 𝐵𝑒 2
𝑥
pero existe una duplicidad 𝑦𝑐 , por lo que multiplicamos 𝑦𝑝2 = 𝐵𝑒 2 por 𝑥 con lo cual
𝑥
𝑦𝑝2 = 𝐵𝑥𝑒 2

Pero nuevamente existe duplicación, procedemos nuevamente a multiplicar por 𝑥, así

𝑥
𝑦𝑝2 = 𝑥 2 𝐵𝑒 2 .

Por lo tanto, finalmente

𝑥
𝑦𝑝 = 𝐴 + 𝐵𝑥 2 𝑒 2

4.4.2 Variación de parámetros

Este procedimiento sin mencionarlo se introdujo en la sección 2.7 (ecuaciones lineales de primer
orden), por lo que procedemos a generalizar el mismo con ecuaciones de orden superior, primero
lo desarrollaremos en una ecuación de segundo orden.

Dada

𝑎2 𝑦 ′′ + 𝑎1 𝑦 ′ + 𝑎0 𝑦 = 𝑔(𝑥)
Primero se encuentra la solución complementaria 𝑦𝑐 , luego se escribe la ecuación diferencial en su
forma estándar

𝑦 ′′ + 𝑃𝑦 ′ + 𝑄𝑦 = 𝑓(𝑥)
Para esta ecuación se busca una solución particular de la forma

𝑦𝑝 = 𝑢1 (𝑥)𝑦1 (𝑥) + 𝑢2 (𝑥)𝑦2

Y para obtener 𝑢1 𝑦 𝑢2 se utilizan las fórmulas

𝑊1 𝑦2 𝑓(𝑥) 𝑊2 𝑦1 𝑓(𝑥)
𝑢1′ = 𝑊
=− 𝑊
y 𝑢2′ = 𝑊
= 𝑊

donde

Elaborado por Javier E. Herrera C. pág.

15
 𝑦1 𝑦2 0 𝑦2 𝑦1 0
𝑊 = |𝑦 ′ ′ |, 𝑊1 = | ′ |, 𝑊2 = | |
1 𝑦2 𝑓(𝑥) 𝑦2 𝑦1′ 𝑓(𝑥)
Por último, la solución general sigue siendo 𝑦 = 𝑦𝑐 + 𝑦𝑝 .

Ejemplo 4.21: Utilizando el método de variación de parámetros resuelva

a. 𝑦 ′′ + 𝑦 = 𝑠𝑒𝑐𝑥

b. 3𝑦 ′′ − 6𝑦 ′ + 6𝑦 = 𝑒 𝑥 𝑠𝑒𝑐𝑥

Solución:

𝑦 ′′ + 𝑦 = 𝑠𝑒𝑐𝑥
a. Como la ecuación auxiliar es 𝑚2 + 1 = 0 entonces 𝑚1 = 𝑖 y 𝑚2 = −𝑖; luego

𝑦𝑐 = 𝑐1 𝑐𝑜𝑠𝑥 + 𝑐2 𝑠𝑒𝑛𝑥, 𝑦1 = 𝑐𝑜𝑠𝑥, 𝑦2 = 𝑠𝑒𝑛𝑥

Luego,
𝑐𝑜𝑠𝑥 𝑠𝑒𝑛𝑥
𝑊=| | = 𝑐𝑜𝑠 2 𝑥 + 𝑠𝑒𝑛2 𝑥 = 1
−𝑠𝑒𝑛𝑥 𝑐𝑜𝑠𝑥
𝑦 ′′ + 𝑦 = 𝑠𝑒𝑐𝑥
Identificando 𝑓(𝑥) = 𝑠𝑒𝑐𝑥, obtenemos
0 𝑦2
| |
𝑊1 𝑓(𝑥) 𝑦2′ 0 𝑠𝑒𝑛𝑥
𝑢1′ = = =| | = −𝑠𝑒𝑐𝑥 𝑠𝑒𝑛𝑥 = −𝑡𝑎𝑛𝑥,
𝑊 1 𝑠𝑒𝑐𝑥 𝑐𝑜𝑠𝑥
𝑢1′ = −𝑡𝑎𝑛𝑥

𝑢1 = − ∫ 𝑡𝑎𝑛𝑥 𝑑𝑥 = 𝑙𝑛|𝑐𝑜𝑠𝑥|

𝑊2 𝑦 0 𝑐𝑜𝑠𝑥 0
𝑢2′ = = | 1′ |=| | = 𝑐𝑜𝑠𝑥 𝑠𝑒𝑐𝑥 = 1
𝑊 𝑦1 𝑓(𝑥) −𝑠𝑒𝑛𝑥 𝑠𝑒𝑐𝑥
𝑢2′ = 1
Integrando tenemos,

𝑢2 = ∫ 𝑑𝑥 = 𝑥

𝑢1 = 𝑙𝑛|𝑐𝑜𝑠𝑥| y 𝑢2 = 𝑥

Así,

𝑦𝑝 = 𝑢1 𝑦1 + 𝑢2 𝑦2 = 𝑐𝑜𝑠𝑥 𝑙𝑛|𝑐𝑜𝑠𝑥| + 𝑥 𝑠𝑒𝑛𝑥

𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1 𝑐𝑜𝑠𝑥 + 𝑐2 𝑠𝑒𝑛𝑥𝑐𝑜𝑠𝑥 + 𝑐𝑜𝑠𝑥 𝑙𝑛|𝑐𝑜𝑠𝑥| + 𝑥 𝑠𝑒𝑛𝑥

Elaborado por Javier E. Herrera C. pág.

16
 3𝑦 ′′ − 6𝑦 ′ + 6𝑦 = 𝑒 𝑥 𝑠𝑒𝑐𝑥

b. La ecuación auxiliar es 3𝑚2 − 6𝑚 + 6 = 0, 𝑚1 = 1 + 𝑖, 𝑚2 = 1 − 𝑖 por lo que

𝑦𝑐 = 𝑒 𝑥 (𝑐1 𝑐𝑜𝑠𝑥 + 𝑐2 𝑠𝑒𝑛𝑥)
y
𝑒 𝑥 𝑐𝑜𝑠𝑥 𝑒 𝑥 𝑠𝑒𝑛𝑥
𝑊=| | = 𝑒 2𝑥
𝑒 𝑐𝑜𝑠𝑥 − 𝑒 𝑥 𝑠𝑒𝑛𝑥
𝑥
𝑒 𝑐𝑜𝑠𝑥 + 𝑒 𝑥 𝑠𝑒𝑛𝑥
𝑥

Por otro lado, escribiendo la ecuación diferencial en su forma estándar se tiene

1
𝑦 ′′ − 2𝑦 ′ + 2𝑦 = 𝑒 𝑥 𝑠𝑒𝑐𝑥
3
Por lo tanto,
0 𝑒 𝑥 𝑠𝑒𝑛𝑥
|1 𝑥 | 1
𝑒 𝑠𝑒𝑐𝑥 𝑒 𝑥 𝑐𝑜𝑠𝑥 + 𝑒 𝑥 𝑠𝑒𝑛𝑥 ( 𝑒 𝑥 𝑠𝑒𝑐𝑥) (𝑒 𝑥 𝑠𝑒𝑛𝑥) 1
𝑢1′ = − 3 =− 3 = − 𝑡𝑎𝑛𝑥
𝑒 2𝑥 𝑒 2𝑥 3
Luego
1
𝑢1 = 𝑙𝑛|𝑐𝑜𝑠𝑥|,
3
𝑒 𝑥 𝑐𝑜𝑠𝑥 0
| 𝑥 1 𝑥 | 1
𝑒 𝑐𝑜𝑠𝑥 − 𝑒 𝑥 𝑠𝑒𝑛𝑥 𝑒 𝑠𝑒𝑐𝑥 (𝑒 𝑥 𝑐𝑜𝑠𝑥) ( 𝑒 𝑥 𝑠𝑒𝑐𝑥) 1
𝑢2′ = 3 = 3 =
𝑒 2𝑥 𝑒 2𝑥 3
Luego,
𝑥
𝑢2 =
3
Por lo que,
1 𝑥
𝑦𝑝 = 𝑒 𝑥 𝑐𝑜𝑠𝑥 𝑙𝑛|𝑐𝑜𝑠𝑥| + 𝑒 𝑥 𝑠𝑒𝑛𝑥
3 3
Así,
1 𝑥
𝑦 = 𝑒 𝑥 (𝑐1 𝑐𝑜𝑠𝑥 + 𝑐2 𝑠𝑒𝑛𝑥) + 𝑒 𝑥 𝑐𝑜𝑠𝑥𝑙𝑛|𝑐𝑜𝑠𝑥| + 𝑒 𝑥 𝑠𝑒𝑛𝑥
3 3
El método que acabamos de describir se puede generalizar a ecuaciones de orden superior luego
que sean escrito en forma estándar.

𝑦 (𝑛) + 𝑃𝑛−1 (𝑥)𝑦 (𝑛−1) + . . . +𝑃1 (𝑥)𝑦 ′ + 𝑃0 (𝑥)𝑦 = 𝑓(𝑥) (*)

Si 𝑦𝑐 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + . . . +𝑐𝑛 𝑦𝑛 es la solución complementaria de (*) entonces una solución

particular es

𝑦𝑝 = 𝑢1 (𝑥)𝑦1 (𝑥) + 𝑢2 𝑦2 (𝑥) + . . . +𝑢𝑛 (𝑥)𝑦𝑛 (𝑥)

Elaborado por Javier E. Herrera C. pág.

17
En este caso,
𝑊𝑘
𝑢𝑘′ = , 𝑘 = 1,2,3, . . . , 𝑛
𝑊
donde 𝑊 es el Wrosnkiano de 𝑦1 , 𝑦2 , . . . , 𝑦𝑛 y 𝑊𝑘 es el determinante que se obtiene al
reemplazar la 𝑘 − é𝑠𝑖𝑚𝑎 columna del Wrosnkiano por la columna que consta de (0,0,0, . . . , 𝑓(𝑥)).

Ejemplo 4.22: Resuelva la ecuación 𝑦 ′′′ + 𝑦 ′ = 𝑡𝑎𝑛𝑥

Solución: La ecuación auxiliar es 𝑚3 + 𝑚 = 0,𝑚1 = 0, 𝑚2 = −𝑖, 𝑚3 = 𝑖

por lo tanto; 𝑦𝑐 = 𝑐1 + 𝑐2 𝑐𝑜𝑠𝑥 + 𝑐3 𝑠𝑒𝑛𝑥 , 𝑦1 = 1, 𝑦2 = 𝑐𝑜𝑠𝑥, 𝑦3 = 𝑠𝑒𝑛𝑥, y

1 𝑐𝑜𝑠𝑥 𝑠𝑒𝑛𝑥
𝑊 = |0 −𝑠𝑒𝑛𝑥 𝑐𝑜𝑠𝑥 | = 1
0 −𝑐𝑜𝑠𝑥 −𝑠𝑒𝑛𝑥
Como 𝑓(𝑥) = 𝑡𝑎𝑛𝑥, obtenemos

𝑊1 0 𝑐𝑜𝑠𝑥 𝑠𝑒𝑛𝑥
𝑢1′ = =| 0 −𝑠𝑒𝑛𝑥 𝑐𝑜𝑠𝑥 | = 𝑡𝑎𝑛𝑥
𝑊
𝑡𝑎𝑛𝑥 −𝑐𝑜𝑠𝑥 −𝑠𝑒𝑛𝑥

𝑊2 1 0 𝑠𝑒𝑛𝑥
𝑢2′ = = |0 0 𝑐𝑜𝑠𝑥 | = −𝑠𝑒𝑛𝑥
𝑊
0 𝑡𝑎𝑛𝑥 −𝑠𝑒𝑛𝑥

𝑊3 1 𝑐𝑜𝑠𝑥 0 𝑐𝑜𝑠 2 𝑥 − 1
𝑢3′ = = |0 −𝑠𝑒𝑛𝑥 0 | = −𝑠𝑒𝑛𝑥 𝑡𝑎𝑛𝑥 = = 𝑐𝑜𝑠𝑥 − 𝑠𝑒𝑐𝑥
𝑊 𝑐𝑜𝑠𝑥
0 −𝑐𝑜𝑠𝑥 𝑡𝑎𝑛𝑥
Como,

𝑢1′ = 𝑡𝑎𝑛𝑥, 𝑢2′ = −𝑠𝑒𝑛𝑥 y 𝑢3′ = 𝑐𝑜𝑠𝑥 − 𝑠𝑒𝑐𝑥, entonces

𝑢1 = −𝑙𝑛|𝑐𝑜𝑠𝑥|
𝑢2 = 𝑐𝑜𝑠𝑥
𝑢3 = 𝑠𝑒𝑛𝑥 − 𝑙𝑛|𝑠𝑒𝑐𝑥 + 𝑡𝑎𝑛𝑥 |
𝑦𝑝 = 𝑢1 𝑦1 + 𝑢2 𝑦2 + 𝑢3 𝑦3 = −𝑙𝑛|𝑐𝑜𝑠𝑥| + 𝑐𝑜𝑠 2 𝑥 + 𝑠𝑒𝑛2 𝑥 − 𝑠𝑒𝑛𝑥 𝑙𝑛|𝑠𝑒𝑐𝑥 + 𝑡𝑎𝑛𝑥|

𝑦 = 𝑦𝑐 + 𝑦𝑝 = 𝑐1 + 𝑐2 𝑐𝑜𝑠𝑥 + 𝑐3 𝑠𝑒𝑛𝑥 − 𝑙𝑛|𝑐𝑜𝑠𝑥| + 𝑐𝑜𝑠 2 𝑥 + 𝑠𝑒𝑛2 𝑥 − 𝑠𝑒𝑛𝑥 𝑙𝑛|𝑠𝑒𝑐𝑥 + 𝑡𝑎𝑛𝑥|

𝑦 = 𝑐1 + 𝑐2 𝑐𝑜𝑠𝑥 + 𝑐3 𝑠𝑒𝑛𝑥 − 𝑙𝑛|𝑐𝑜𝑠𝑥| + 1 − 𝑠𝑒𝑛𝑥 𝑙𝑛|𝑠𝑒𝑐𝑥 + 𝑡𝑎𝑛𝑥|

𝑦 = 𝑐4 + 𝑐2 𝑐𝑜𝑠𝑥 + 𝑐3 𝑠𝑒𝑛𝑥 − 𝑙𝑛|𝑐𝑜𝑠𝑥| − 𝑠𝑒𝑛𝑥 𝑙𝑛|𝑠𝑒𝑐𝑥 + 𝑡𝑎𝑛𝑥|
donde 𝑐4 = 𝑐1 + 1

Elaborado por Javier E. Herrera C. pág.

18