Ecuaciones Diferenciales Ordinarias

Lineares de Orden Superior

Raquel Serna
Objetivos

• Integración
• Ecuaciones cuadráticas

1/95
Método de coeficientes indeterminados.
Este método solo podrá aplicarse en el caso de tener una ecuación
ordinaria lineal se segundo orden de la forma

y 00 + ay 0 + by = r(x), (28)

donde a, b ∈ R y si la función r (x) toma las siguientes formas

• Exponencial: r(x) = eαx , con α ∈ R
• Polinomial: r(x) = p0 + p1 x + p2 x2 + · · · + pn xn , donde pi ∈ R y
n ∈ N∪ {0}
• Trigonométricos: r (x) = cos (ωx) o r (x) = sen (ωx) con ω ∈ R.
• Combinaciones: r (x) = eαx cos (ωx) o r (x) = eαx sen (ωx) o
r (x) = eαx (p0 + p1 x + p2 x2 + · · · + pn xn ) o
r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) cos (ωx) o
r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) sen (ωx) con ω, pi , α ∈ R.
41/95
Función exponencial r (x) = eαx , con α ∈ R
• Si λ = α no es una solución de la ecuación característica de la EDO
homogénea asociada entonces la solución particular de la EDO tendrá
la forma yp = Aeαx
• Si λ = α es una solución simple de la ecuación característica de la
EDO homogénea asociada entonces yp = Axeαx ,
• Si λ = α es una solución repetida de la ecuación característica de la
EDO homogénea asociada entonces yp = Ax2 eαx ,
para un A ∈ R indeterminado que posteriormente sera calculado
sustituyendo yp en la ecuación (28).
Ejemplo: Resolver y 00 + 3y 0 + 2y = 2e3x
Solución:
Observemos que la ecuación característica correspondiente a la EDO
homogénea asociada es
λ2 + 3λ + 2 = 0 =⇒ λ = −1, −2. 42/95
Luego y1 = e−x e y2 = e−2x son soluciones l.i. de la EDO homogénea
asociada y la solución particular tendrá la forma yp = Ae3x entonces
• yp0 = 3Ae3x
• yp00 = 9Ae3x
ya que se debe cumplir

yp00 + 3yp0 + 2yp = 2e3x

obtenemos(9Ae3x ) + 3 (3Ae3x ) + 2 (Ae3x ) = 2e3x , entonces

20Ae3x = 2e3x

e3x
=⇒ 20A = 2 =⇒ A = 101 , luego la solución particular es yp (x) = 10
.
Finalmente la solución general será

e3x
y (x) = c1 e−x + c2 e−2x +
10
43/95
donde c1 , c2 son constante indeterminadas.
Función polinomial r (x) = p0 + p1 x + p2 x2 + · · · + pn xn ; n ≥ 0
• Si λ = 0 no es una raíz de la ecuación característica de la EDO
homogénea asociada, entonces la solución particular de la EDO
tendrá la forma yp = q0 + q1 x + q2 x2 + · · · + qn xn .
• Si λ = 0 es una raíz simple de la ecuación característica de la EDO
homogénea asociada, entonces la solución particular de la EDO
tendrá la forma yp = x (q0 + q1 x + q2 x2 + · · · + qn xn ) .
• Si λ = 0 es una raiz repetida de la ecuación característica de la EDO
homogénea asociada, entonces la solución particular de la EDO
tendrá la forma yp = x2 (q0 + q1 x + q2 x2 + · · · + qn xn ) .
donde qi ∈ R son valores indeterminados que posteriormente serán
calculados sustituyendo yp en la ecuación (28).
Ejemplo: Resolver y 00 + 3y 0 = 4ex + 1 − 2x + e−3x
Solución:
Observemos que la ecuación característica correspondiente a la EDO
44/95
homogénea asociada es λ2 + 3λ = 0
entonces λ = 0, −3, luego dos soluciones linealmente
independientes de la EDO homogénea asociada serán
y1 = e0x = 1 e y2 = e−3x , seguidamente subdividimos el problema
en tres EDOs no homogéneas
• Subproblema 1: y 00 + 3y 0 = 4ex = r1 (x)
• Subproblema 2: y 00 + 3y 0 = 1 − 2x = r2 (x)
• Subproblema 3: y 00 + 3y 0 = e−3x = r3 (x)
Por tanto, la solución particular yp de la EDO dada tiene la forma

−3x
yp = ⏟Ae⏞x + x (B + Cx) + Dxe
⏟ ⏞
⏟ ⏞
yp1 yp2 yp3

donde yp1 , yp2 , yp3 representan las soluciones particulares de los

subproblemas 1,2 y 3 respectivamente, seguidamente derivamos yp
• yp0 = Aex + B + 2Cx + De−3x − 3Dxe−3x
• yp00 = Aex + 2C − 3De−3x − 3De−3x + 9Dxe−3x 45/95
dado que debe verificarse que

yp00 + 3yp0 = 4ex + 1 − 2x + e−3x

obtenemos
(Aex + 2C − 3De−3x − 3De−3x + 9Dxe−3x ) +
3 (Aex + B + 2Cx + De−3x − 3Dxe−3x ) = 4ex + 1 − 2x + e−3x entonces

2C + 3B + 6Cx + 4Aex − 3De−3x = 4ex + 1 − 2x + e−3x

consecuentemente 2C + 3B = 1; 6C = −2; 4A = 4 y −3D = 1 =⇒ A = 1;

D = − 31 ; C = − 31 y B = 95 , luego la solución particular será
5x x2 xe−3x
yp (x) = ex + 9
− 3
− 3
. Finalmente, la solución general será

5x x2 xe−3x
y (x) = c1 + c2 e−3x + ex + − −
9 3 3
donde c1 , c2 son constantes indeterminadas.
46/95
Función trigonométrica r (x) = cos (ωx) o r (x) = sen (ωx)
• Si λ = ±iω no es una raíz de la ecuación característica de la EDO
homogénea asociada, entonces la solución particular de la EDO
tendrá la forma yp = Acos (ωx) + Bsen (ωx).
• Si λ = ±iω es una raíz de la ecuación característica de la EDO
homogénea asociada, entonces la solución particular de la EDO
tendrá la forma yp = x (Acos (ωx) + Bsen (ωx)).
donde A, B ∈ R son valores indeterminado que posteriormente serán
calculados sustituyendo yp en la ecuación (28).
Ejemplo: y 00 + 2y 0 + y = 2senx − 3e−x + 1
Solución:
Observemos que la ecuación característica correspondiente a la EDO
homogénea asociada es λ2 + 2λ + 1 = 0 entonces λ = −1 (doble) . Luego

y1 = e−x e y2 = xe−x ,
47/95
después subdividimos el problema en tres EDOs no homogéneas
• Subproblema 1: y 00 + 2y 0 + y = 2senx = r1 (x)
• Subproblema 2: y 00 + 2y 0 + y = −3xsenx = r2 (x)
• Subproblema 3: y 00 + 2y 0 + y = 1= r3 (x)
Por tanto la solución particular de la EDO dada tiene la forma

yp = Asenx
⏟ +⏞ Bcosx + ⏟Cx2 e⏞−x + ⏟ D⏞
yp1 yp2 yp3

donde yp1 , yp2 , yp3 representan las soluciones particulares de los

subproblemas 1, 2 y 3 respectivamente, seguidamente derivamos yp
• yp0 = Acosx − Bsenx + 2Cxe−x − Cx2 e−x
• yp00 = −Asenx − Bcosx + 2Ce−x − 2Cxe−x − 2Cxe−x + Cx2 e−x

48/95
dado que debe verificarse que

yp00 + 2yp0 + yp = 2senx − 3e−x + 1

obtenemos
(−Asenx − Bcosx + 2Ce−x − 2Cxe−x − 2Cxe−x + Cx2 e−x ) +
2 (Acosx − Bsenx + 2Cxe−x − Cx2 e−x ) + Asenx + Bcosx + Cx2 e−x + D =
2senx − 3e−x + 1
entonces

(−2B) senx + 2Acosx + 2Ce−x + D = 2senx − 3e−x + 1

−3
luego −2B = 2, 2A = 0, 2C = −3 y D = 1 =⇒ B = −1;A = 0; C = 2
y
D=1
2 −x
por tanto la solución particular será yp (x) = −cos (x) − 3x 2e + 1
consecuentemente la solución general de la EDO dada será

3x2 e−x
y (x) = c1 e−x + c2 xe−x − cos (x) − +1
2
49/95
donde c1 , c2 son constantes indeterminadas.
En los casos anteriores se vio como encontrar la solución particular de la
EDO no homogénea en el caso de que r (x) tenga la forma de una función
exponencial, polinómica , de seno, coseno o la combinación lineal de ellos
(sumas, restas). Ahora veremos que este método también puede aplicarse
cuando r(x) es obtenido como multiplicación de estas funciones, es decir
• Si r (x) = eαx cos (ωx) o r (x) = eαx sen (ωx) y λ = α ± ωi no son
raíces de la ecuación característica de la EDO homogénea asociada,
entonces yp = Aeαx cos (ωx) + Beαx sen (ωx).
• Si r (x) = eαx cos (ωx) o r (x) = eαx sen (ωx) y λ = α ± ωi son raíces
de la ecuación característica de la EDO homogénea asociada,
entonces yp = x (Aeαx cos (ωx) + Beαx sen (ωx)).
• Si r (x) = eαx (p0 + p1 x + p2 x2 + · · · + pn xn ) y si λ = α no es raíz de
la ecuación característica de la EDO homogénea asociada, entonces
yp = eαx (q0 + q1 x + q2 x2 + · · · + qn xn ).
• Si r (x) = eαx (p0 + p1 x + p2 x2 + · · · + pn xn ) y si λ = α es raíz simple
de la ecuación característica de la EDO homogénea asociada,
entonces yp = xeαx (q0 + q1 x + q2 x2 + · · · + qn xn ). 50/95
 • Si r (x) = eαx (p0 + p1 x + p2 x2 + · · · + pn xn ) y si λ = α es raíz
repetida de la ecuación característica de la EDO homogénea asociada,
entonces yp = x2 eαx (q0 + q1 x + q2 x2 + · · · + qn xn ).
• Si r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) cos (ωx) o
r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) sen (ωx) y λ = ± ωi no son
raíces de la ecuación característica de la EDO homogénea asociada,
entonces

yp = q0 + q1 x + q2 x2 + · · · + qn xn cos (ωx) + r0 + r1 x + r2 x2 + · · · + rn xn sen (ωx) .

(︀ )︀ (︀ )︀

• Si r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) cos (ωx) o

r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) sen (ωx) y λ = ± ωi son
raíces de la ecuación característica de la EDO homogénea asociada,
entonces

yp = x q0 + q1 x + q2 x2 + · · · + qn xn cos (ωx) + x r0 + r1 x + r2 x2 + · · · + rn xn sen (ωx)

(︀ )︀ (︀ )︀

51/95
Ejemplo: Resolver y 00 + 2y 0 + 2y = e−x sen(x) con
Solución:

52/95
53/95
Ejemplo: Resolver y 00 + y = (x + 1) senx
Solución:
En este caso observamos que la ecuación característica correspondiente es

λ2 + 1 = 0 =⇒ λ = ±i

por otro lado, observamos que el termino independiente de la EDO se

obtiene como multiplicación de un polinomio de primer grado y la función
seno, entonces usando el metodo de coeficientes indeterminados tenemos
que la solución particular yp tendra la siguiente estructura

yp = x (Ax + B) sen (x) + x (Dx + E) cos (x)

luego
yp0 = (2Ax + B) senx + (Ax2 + Bx) cosx + (2Dx + E) cosx − (Dx2 + Ex) senx
= ((2A − E) x − Dx2 + B) senx + (Ax2 + (B + 2D) x + E) cosx

54/95
seguidamente
yp00 = (2A − E − 2Dx) senx + ((2A − E) x − Dx2 + B) cosx +
(2Ax + B + 2D) cosx − (Ax2 + (B + 2D) x + E) senx
= (−Ax2 − (4D + B) x − 2E) senx + (−Dx2 + (4A − E) x + 2B + 2D) cosx
luego, ya que

yp00 + yp = (x + 1) senx

obtenemos
(−4Dx − 2E) senx + (4Ax + 2B + 2D) cosx = (x + 1) senx
entonces
−4D = 1;−2E = 1 ; 4A = 0; 2B + 2D = 0
luego D = − 41 ; E = − 21 ; A = 0 ; B = 41
x2 cosx
luego la solución particular será yp = xsenx
4
− 4
− xcosx
2
luego la solución general de la EDO sera

xsenx x2 cosx xcosx

y(x) = c1 cosx + c2 senx + − −
4 4 2
donde c1 y c2 son constante indeterminadas. 55/95
Ejemplo: Resolver y 00 + 3y 0 + 2y = xe−x + 2e−x senx
Solución
En este caso la ecuación característica a la EDO homogenea asociada es:
λ2 + 3λ + 2 = 0
resolviendo obtenemos λ = −1, −2. Luedo dos soluciones l.i. de la EDO
homogénea es y1 = e−x e y2 = e−2x . Para resolver la EDO no homogénea
dada, tenemos que dividirla en dos EDOs no homogéneas
y 00 + 3y 0 + 2y = xe−x = r1 (x) y y 00 + 3y 0 + 2y = 2e−x sen (x) = r2 (x),
Usando el método de coeficiente indeterminados, determinaremos la
estructura de las respectivas soluciones particulares,es decir
yp1 (x) = x (Ax + B) e−x
yp2 (x) = De−x cos (x) + Ee−x sen (x)
consecuentemente la solución particular yp de la EDO dada será
yp = x (Ax + B) e−x + De−x cos (x) + Ee−x sen (x)

56/95
derivando obtenemos
yp0 = e−x −Ax2 + (2A − B) x + B
[︀ ]︀
1
yp00 = e−x Ax2 + (B − 4A) x + 2A − 2B
[︀ ]︀
1
yp0 = e−x [(E − D) cos (x) − (D + E) sen (x)]
2
yp00 = e−x [2Dsen (x) − 2Ecos (x)]
2
dado que yp00 + 3yp0 + 2yp1 = xe−x obtenemos
1 1
e−x Ax2 + (B − 4A) x + 2A − 2B +
[︀ ]︀

3e−x −Ax2 + (2A − B) x + B + 2x (Ax + B) e−x = xe−x

[︀ ]︀

luego
e−x [2Ax + 2A + B] = xe−x
=⇒ 2A = 1; 2A + B = 0 =⇒ A = 21 y B = −1. Por tanto,
(︁ 2 )︁
yp1 (x) = x2 − x e−x . Por otro lado dado que
yp00 + 3yp0 + 2yp2 = 2e−x sen (x)
2 2
obtenemos
57/95
e−x [2Dsen (x) − 2Ecos (x)] +
3e−x [ (E − D) cos (x) − (D + E) sen (x)] +
2 De−x cos (x) + Ee−x sen (x) = 2e−x sen (x)
[︀ ]︀

=⇒ e−x [(E − D) cos (x) − (D + E) sen (x)] = 2e−x sen (x)

=⇒ E − D = 0 y −D − E = 2 =⇒ D = E = −1. Por tanto,
yp2 (x) = e−x [−cos (x) − sen (x)], luego la solución particular de la
EDO no(︁ 2homogénea
)︁ dada es
x
yp = 2 − x e + e−x [−cos (x) − sen (x)] . Finalmente la
−x

solución general será (︁ 2 )︁

yp (x) = c1 e−x + c2 e−2x + e−x x2 − x − cos (x) − sen (x)
donde c1 y c2 son constantes indeterminadas.

58/95
Ejercicios Propuestos

1. Resolver los siguientes PVI

a) y 00 + 3y 0 + 2y = ex + e2x , y (0) = 0, y 0 (0) = 0
b) y 00 + 4y 0 + 4y = e−2x , y (0) = 1, y 0 (0) = 0
c) y 00 + 4y = 2sen (2x) , y (0) = 0, y 0 (0) = −1
2. Use el método de variación de parámetros para resolver los
siguientes PVI
a) y 00 + y = sec (x) , y (0) = 0, y 0 (0) = 0
b) y 00 + 3y 0 + 2y = sen (ex ) , y (0) = 1, y 0 (0) = 1
ex
c) y 00 − y = , y (0) = y 0 (0) = 0
1 + e2x

59/95
 ⎧
⎨y 00 + y = r (x)
3. Considere el problema
⎩y (0) = y 0 (0) = 0
∫︀ x
a) Verifique que la función ϕ (x) = 0 sen (x − t) r (t) dt es una
solución del problema dado.
b) Obtenga esta formula usando el método de variación de
parámetros.
4. Muestre que la EDO x2 y 00 + xy 0 + x2 − 1/ 4 y = 0 para x > 0
(︀ )︀

puede ser reducida a la EDO v00 + v = 0 mediante la sustitución

y = x−1/ 2 v (x).
Determine la solución general de la EDO original.

60/95
Ecuaciones diferenciales ordinarias lineales de
orden superior
Una EDO lineal de orden superior es una ecuación de la forma

d(n) y d(n−1) y dy
a0 (x) + a1 (x) + · · · + an−1 (x) + an (x)y = F(x) (29)
dxn dxn−1 dx
donde a0 , a1 , · · · , an−1 , an , F son funciones continuas en el intervalo I ⊂ R,
con a0 (x) 6= 0 para todo x ∈ I.

Problemas de valores iniciales (PVIs)

Un PVI es un problema formado por una ecuación diferencial de orden n
junto con los valores fijos que debe tomar la solución y sus derivadas en un
mismo punto x = x0 .

y (x0 ) = y0 , y 0 (x0 ) = y1 , y 00 (x0 ) = y2 , · · · , y (n−1) (x0 ) = yn−1 ,

61/95
Problema de valores en la frontera
Es el formado por la ecuación diferencial junto con los valores fijados
que debe tomar la solución de la EDO en distintos x( x = x0 ,
x = x1 , . . . , x = xn−1 ), es decir

y (x0 ) = y0 , y (x1 ) = y1 , y (x2 ) = y2 , · · · , y (xn−1 ) = yn−1 .

62/95
Teorema (Existencia de una solución única)
Sean an (x) , an−1 (x) , · · · , a1 (x) , a0 (x) y g (x) funciones continuas
en el intervalo abierto I y sea a0 (x) 6= 0 para toda x en este
intervalo. Si x = x0 es cualquier punto en este intervalo, entonces
para cualesquiera valores de y0 , y1 , y2 , · · · , yn−1 ∈ R existe una
solución única y(x) de la EDO (31) en el intervalo I tal que

y(x0 ) = y0 , y 0 (x0 ) = y1 , y 00 (x0 ) = y2 , · · · , y (n−1) (x0 ) = yn−1 ,

63/95
Ecuaciones diferenciales ordinarias lineales de
orden superior con coeficientes constantes

Una EDO lineal de orden superior con coeficientes constantes es una

ecuación de la forma

y (n) + a1 y (n−1) + · · · + an y = r(x) (30)

donde a1 , · · · , an son constantes y r es una función que depende de la

variable independiente x.

64/95
Ecuaciones homogéneas Si en la ecuación anterior r(x) = 0 para
toda x, entonces la ecuación se llamará homogénea, es decir

y (n) + a1 y (n−1) + · · · + an y = 0 (31)

donde a1 , · · · , an son constantes. El polinomio característico p (λ)

correspondiente a esta ecuación es

P (λ) = λn + a1 λn−1 + · · · + an . (32)

luego la ecuación característica asociada a la EDO homogénea sera

p(λ) = 0

es decir
λn + a1 λn−1 + · · · + an = 0
65/95
Proposición

• Si λ es una raíz simple de la ecuación característica p (λ) = 0

entonces y = eλx es una solución de la EDO (31).
• Si λ es una raíz repetida de multiplicidad r, entonces todas las
funciones y = xk−1 eλx , para k = 1, 2, · · · , r son soluciones de la
ecuación (31).
• Si λ = σ + iω es una raíz simple compleja de la ecuación
característica p (λ) = 0 entonces y = eσx cos (ωx) e
y = eσx sen (ωx) son dos soluciones de la ecuación(31)
• si λ = σ + iω es una raíz compleja con multiplicidad r, entonces
todas la funciones y = xk−1 eσx cos (ωx) e y = xk−1 eσx sen (ωx),
para k = 1, 2, · · · , r son soluciones de la EDO (31).

66/95
Sistema fundamental de soluciones
Un conjunto de n funciones {y1 , y2 , · · · , yn } es llamado sistema
fundamental de soluciones de la EDO (31) en ]a, b[ si son
linealmente independientes en el intervalo ]a, b[ .
NOTA IMPORTANTE:Dado un conjunto fundamental de la EDO(31)
{y1 , y2 , · · · , yn } la solución general de la EDO (31) tendrá la forma

y (x) = c1 y1 + c2 y2 + · · · + cn yn

donde c1 , c2 , · · · , cn son constanetes indeterminadas.

67/95
Ejemplo: Resolver la EDO y 000 − 6y 00 + 11y 0 − 6y = 0
Solución:

68/95
Ejemplo: Resolver la EDO y (4) − y 0 = 0
Solución:
En este caso la ecuación característica correspondiente seria
λ4 − λ = 0 =⇒ λ (λ3 − 1) = 0 =⇒ λ p ( λ − 1) ( λ + λ + 1) = 0
2
−1± 3i
resolviendo obtenemos λ = 0, 1, 2
luego obtenemos las soluciones:
y1 (x) = 1
y2 (x) = ex (︁ p )︁
x
y3 (x) = e− 2 cos 23x
x
(︁ p )︁
y4 (x) = e− 2 sen 23x
esta soluciones son linealmente independientes(l.i.), consecuentemente la
solución general de la EDO será
(︃ p )︃ (︃ p )︃
− 2x
3x − 2x
3x
y (x) = c1 + c2 e + c3 e
x
cos + c4 e sen
2 2
donde c1 , c2 , c3 y c4 son constantes indeterminadas.
69/95
Ecuaciones no Homogéneas
Una EDO lineal de orden superior con coeficientes constantes no
homogénea es una ecuación de la forma

y (n) + a1 y (n−1) + · · · + an y = r(x) (33)

donde a1 , · · · , an son constantes.

Teorema
Cualquier solución de la EDO (33) tiene la forma

y = c1 y1 + c2 y2 + · · · + Cn yn + yp , (34)

donde c1 , c2 , · · · , cn ∈ R e y1 , y2 , · · · yn son soluciones linealmente

independientes de la ecuación homogénea asociada a la EDO (31) e yp es
una solución particular de la EDO (33). La solución (34) será llamada
solución general de la ecuación (33). 70/95
A continuación presentamos dos métodos para encontrar una
solución particular yp de la EDO no homogénea.

71/95
Método de Variación de parámetros
Asumamos que y1 , y2 , · · · , yn son n soluciones linealmente independientes
de EDO (31). Definamos

yp = c1 (x) y1 + c2 (x) y2 + c3 (x) y3 + · · · + cn (x) yn ,

donde ci (x) son funciones por determinar, además para simplificar cuentas
supondremos que

y1 c10 + y2 c20 + y3 c30 + · · · + yn cn0 = 0

y10 c10 + y20 c20 + y30 c30 + · · · + yn0 cn0 = 0
y1 c1 + y200 c20 + y300 c30 + · · · + yn00 cn0 = 0
00 0

..
.
y1(n−2) c10 + y2(n−2) c20 + y3(n−2) c30 + · · · + yn(n−2) cn0 = 0

Luego sustituyendo yp , yp0 , yp00 , · · · , ypn en la ecuación no homogénea (29) tal

que 72/95
 yp(n) + a1 yp(n−1) + · · · + an yp = r(x) (35)

se verificará que

y1(n−1) c10 + y2(n−1) c20 + y3(n−1) c30 + · · · + yn(n−1) cn0 = r (x) .

entonces se tiene un sistema de n ecuaciones y n incógnitas ci0 (x) por

resolver. Luego de hacer diversos cálculos obtenemos que

(−1)n−i W (y1 , y2 , · · · , yi−1 , yi+1 , · · · , yn ) r(x)

ci0 (x) =
W (y1 , y2 , · · · , yn )

para cada i ∈ {1, 2, · · · , n}, donde W (y1 , y2 , · · · , yn ) es el wronskiano de

las funciones y1 , y2 , · · · , yn . Por tanto,
∫︁
(−1)n−i W (y1 , y2 , · · · , yi−1 , yi+1 , · · · , yn ) r(x)
ci (x) = dx
W (y1 , y2 , · · · , yn )

para cada i ∈ {1, 2, · · · , n} 73/95

NOTA IMPORTANTE: La determinante de una matriz de 3 × 3 se calcula
como sigue
⃒ ⃒
⃒ a
⃒ 11 a12 a13 ⃒
⃒ ⃒ ⃒ ⃒ ⃒ ⃒ ⃒
⃒ a
⃒ 22 a23 ⃒
⃒ a
⃒ 21 a23 ⃒
⃒ a
⃒ 21 a22
⃒ ⃒ ⃒
⃒ a21 a22 a23 ⃒ = a11 ⃒ ⃒ − a12 ⃒ ⃒ + a13 ⃒
⃒ ⃒ ⃒
⃒
⃒ ⃒ ⃒ a32 a33 ⃒ ⃒ a31 a33 ⃒ ⃒ a31 a32 ⃒
⃒ a31 a32 a33 ⃒

74/95
 4x
2e
Ejemplo: Resolver la EDO y 000 − 6y 00 + 11y 0 − 6y = 1+e x

Solución:
Notemos que la ecuación característica correspondiente a la EDO
homogénea asociada es

λ3 − 6λ2 + 11λ − 6 = 0,

luego λ = 1, 2, 3 , entonces obtenemos que y1 = ex , y2 = e2x e y3 = e3x son

tres soluciones l.i. de la EDO homogénea asociada, luego buscamos la
solución particular yp haciendo uso del método de variación de parámetros,
es decir
yp = c1 (x) y1 + c2 (x) y2 + c3 (x) y3 ,

por tanto tendremos que resolver el sistema

y1 c10 + y2 c20 + y3 c30 = 0

y10 c10 + y20 c20 + y30 c30 = 0
y100 c10 + y200 c20 + y300 c30 = r (x)
75/95
luego
)r(x) (−1)3−2 W(y1 ,y3 )r(x) (−1)3−3 W(y1 ,y2 )r(x)
c10 (x) = (−1)W(yW(y
3−1
2 , y3
, c0 (x) = y c30 (x) =
1 , y2 , y3 ) 2 W(y1 ,y2 ,y3 ) W(y1 ,y2 ,y3 )
0btenemos
⃒ ⃒
⃒ ex e2x e3x ⃒
⃒ ⃒
W (y1 , y2 , y3 ) = ⃒ ex 2e2x 3e3x ⃒ = ex 18e5x − 12e5x −e2x 9e4x − 3e4x +e3x 4e3x − 2e3x =
⃒ ⃒ (︀ )︀ (︀ )︀ (︀ )︀
⃒ x
4e2x 9e3x
⃒
⃒ e ⃒

= 6e6x − 6e6x + 2e6x = 2e6x

y ⃒ ⃒
⃒ e2x e3x ⃒⃒
W(y2 , y3 ) = ⃒ ⃒ = 3e5x − 2e5x = e5x
⃒
⃒ 2e2x 3e3x ⃒
⃒ ⃒
⃒ ex e3x ⃒⃒
W(y1 , y3 ) = ⃒ ⃒ = 3e4x − e4x = 2e4x
⃒
⃒ ex 3e3x ⃒
⃒ ⃒
⃒ ex e2x ⃒⃒
W(y1 , y2 ) = ⃒ ⃒ = 2e3x − e3x = e3x
⃒
⃒ ex 2e2x ⃒
entonces

76/95
 W (y2 , y3 ) r (x) 1 e3x e2x
∫︁ ∫︁ ∫︁
c1 (x) = dx = e−x r (x) dx = dx = − ex + ln (1 + ex ) + ,
W (y1 , y2 , y3 ) 2 1+e x 2

W (y1 , y3 ) r (x) 2e2x

∫︁ ∫︁ ∫︁
c2 (x) = − dx = − e−2x r (x) dx = − dx = 2ln (1 + ex ) − 2ex
W (y1 , y2 , y3 ) 1 + ex

y
W (y1 , y2 ) r (x) 1 ex
∫︁ ∫︁ ∫︁
c3 (x) = dx = e−3x r (x) dx = dx = ln (1 + ex ) .
W (y1 , y2 , y3 ) 2 1 + ex
luego
yp = c1 (x) ex + c2 (x) e2x + c3 (x) e3x

e2x
(︁ )︁
yp = 2
− ex + ln (1 + ex ) ex + (2ln (1 + ex ) − 2ex ) e2x + ln (1 + ex ) e3x

3e2x
yp = ex (1 + ex )2 ln (1 + ex ) − − ex
2

Por tanto,
yp = ex (1 + ex )2 ln (1 + ex ) .
77/95
Método de coeficientes indeterminados.
Este método solo podrá aplicarse en el caso de tener una ecuación
ordinaria lineal de orden n de la forma

y (n) + a1 y (n−1) + · · · + an y = r(x) (36)

donde ai ∈ R y si la función r (x) toma las siguientes formas

• Exponencial: r(x) = eαx , con α ∈ R
• Polinomial: r(x) = p0 + p1 x + p2 x2 + · · · + pn xn , donde pi ∈ R y
n ∈ N∪ {0}
• Trigonométricos: r (x) = cos (ωx) o r (x) = sen (ωx) con ω ∈ R.
• Combinaciones: r (x) = eαx cos (ωx) o r (x) = eαx sen (ωx) o
r (x) = eαx (p0 + p1 x + p2 x2 + · · · + pn xn ) o
r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) cos (ωx) o
r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) sen (ωx) con ω, pi , α ∈ R.
78/95
Función exponencial r (x) = eαx , con α ∈ R

• Si λ = α no es una solución de ecuación característica asociada,

entonces la solución particular de la EDO tendrá la forma yp = Aeαx ,
donde A es una constante a determinar.
• Si λ = α es una solución de multiplicidad k de ecuación característica
asociada , entonces la solución particular de la EDO tendrá la forma
yp = Axk eαx ,
donde A es una constante a determinar sustituyendo yp en la ecuación (36).

79/95
Función polinomial r (x) = p0 + p1 x + p2 x2 + · · · + pn xn ; n ≥ 0

• Si λ = 0 no es una raíz de la ecuación característica de la EDO

homogénea asociada, entonces la solución particular de la EDO
tendrá la forma yp = q0 + q1 x + q2 x2 + · · · + qn xn .
• Si λ = 0 es una raíz de multiplicidad k de la ecuación característica de
la EDO homogénea asociada, entonces la solución particular de la
EDO tendrá la forma yp = xk (q0 + q1 x + q2 x2 + · · · + qn xn )
donde qi ∈ R son valores indeterminados que posteriormente serán
calculados sustituyendo yp en la ecuación (36).

80/95
Función trigonométrica r (x) = cos (ωx) o r (x) = sen (ωx)

• Si λ = ±iω no es una raíz de la ecuación característica de la EDO

homogénea asociada, entonces la solución particular de la EDO
tendrá la forma yp = Acos (ωx) + Bsen (ωx).
• Si λ = ±iω es una raíz de multiplicidad k de la ecuación característica
de la EDO homogénea asociada, entonces la solución particular de la
EDO tendrá la forma yp = xk (Acos (ωx) + Bsen (ωx)).
donde A, B ∈ R son valores indeterminado que posteriormente serán
calculados sustituyendo yp en la ecuación (36).

81/95
En los casos anteriores se vio como encontrar la solución particular de la
EDO no homogénea en el caso de que r (x) tenga la forma de una función
exponencial, polinómica , de seno, coseno o la combinación lineal de ellos
(a través de sumas, restas). Ahora veremos que este método también
puede aplicarse cuando r(x) es obtenido como multiplicación de estas
funciones, es decir
• Si r (x) = eαx cos (ωx) o r (x) = eαx sen (ωx) y λ = α ± ωi no son
raíces de la ecuación característica de la EDO homogénea asociada,
entonces yp = Aeαx cos (ωx) + Beαx sen (ωx).
• Si r (x) = eαx cos (ωx) o r (x) = eαx sen (ωx) y λ = α ± ωi son raíces
de multiplicidad k de la ecuación característica de la EDO homogénea
asociada, entonces yp = xk (Aeαx cos (ωx) + Beαx sen (ωx)).
• Si r (x) = eαx (p0 + p1 x + p2 x2 + · · · + pn xn ) y si λ = α no es raíz de
la ecuación característica de la EDO homogénea asociada, entonces
yp = eαx (q0 + q1 x + q2 x2 + · · · + qn xn ).
• Si r (x) = eαx (p0 + p1 x + p2 x2 + · · · + pn xn ) y si λ = α es una raíz de
multiplicidad k de la ecuación característica de la EDO homogénea
asociada, entonces yp = xk eαx (q0 + q1 x + q2 x2 + · · · + qn xn ). 82/95
 • Si r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) cos (ωx) o
r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) sen (ωx) y λ = ± ωi no son
raíces de la ecuación característica de la EDO homogénea asociada,
entonces

yp = q0 + q1 x + q2 x2 + · · · + qn xn cos (ωx) + r0 + r1 x + r2 x2 + · · · + rn xn sen (ωx) .

(︀ )︀ (︀ )︀

• Si r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) cos (ωx) o

r (x) = (p0 + p1 x + p2 x2 + · · · + pn xn ) sen (ωx) y λ = ± ωi son
raíces de multiplidad k de la ecuación característica de la EDO
homogénea asociada, entonces

yp = xk q0 + q1 x + q2 x2 + · · · + qn xn cos (ωx) + xk r0 + r1 x + r2 x2 + · · · + rn xn sen (ωx) .

(︀ )︀ (︀ )︀

83/95
Ejemplo: Resuelva y 000 + y 00 = ex cos (x)
Solución:
Notemos que la ecuación característica correspondiente a la EDO
homogénea asociada es
λ3 + λ2 = 0,

luego λ = 0(doble) y λ = −1, entonces

y1 = 1 , y2 = x e y3 = e−x son tres soluciones l.i. de la EDO homogénea
asociada, luego yh = c1 + c2 x + c3 e−x . Por otro lado la solución particular
tomará la forma
yp = Aex cos(x) + Bex sen(x)

luego
yp0 = ex [(A + B) cos (x) + (B − A) sen (x)]
yp00 = ex [2Bcos (x) − 2Asen (x)]
yp000 = ex [ (2B − 2A) cos (x) − (2A + 2B) sen (x)]

84/95
Dado que yp debe satisfacer la EDO dada tenemos que

yp000 +yp00 = (−2A + 4B) ex cos(x)+(−4A − 2B) ex sen(x) = ex cos(x)

entonces −2A + 4B = 1, −4A − 2B = 0, luego A = − 101 y B = 1

5
consecuentemente
1 1
yp = − ex cos(x) + ex sen(x).
10 5
Por tanto la solución general de la ecuación es

1 1
y = yh + yp = c1 + c2 x + c3 e−x − ex cos(x) + ex sen(x)
10 5

85/95
Ejemplo: Resuelva y (4) + y 000 = 1 − x2 e−x
Solución: Notemos que la ecuación característica correspondiente a la EDO
homogénea asociada es
λ4 + λ3 = 0.

Para encontrar la solución particualar de la EDO subdividimos el problema

en dos como sigue:
y (4) + y 000 = 1 = r1 (x) y y (4) + y 000 = −x2 e−x = r2 (x)
luego las soluciones particulares yp1 , yp2 de las ultimas dos EDO tendran la
siguiente estructura
yp1 = Ax3 e yp2 = x (B + Cx + Dx2 ) e−x
luego la solución particualr de la EDO dada tendrá la forma

yp = yp1 + yp2

entonces
yp (x) = Ax3 + Bx3 + Cx2 + Dx e−x
(︀ )︀

yp0 (x) = 3Ax2 + (3Bx2 + 2Cx + D) e−x − (Bx3 + Cx2 + Dx) e−x =
3Ax2 + (D + (2C − D)x + (3B − C) x2 − Bx3 ) e−x 86/95
yp00 = 6Ax + 2C − 2D + (6B − 4C + D) x + (C − 6B) x2 + Bx3 e−x
[︀ ]︀

yp000 = 6A+ 6B − 6C + 3D + (6C − 18B − D) x + (9B − C) x2 − Bx3 e−x

[︀ ]︀

yp(4) = 12C − 24B − 4D + (36B − 8C + D) x + (C − 12B) x2 + Bx3 e−x

[︀ ]︀

entonces como

yp(4) + yp000 = 6A + 6C − 18B − D + (18B − 2C) x − 3Bx2 e−x

[︀ ]︀

87/95
tenemos que

6A + 6C − 18B − D + (18B − 2C) x − 3Bx2 e−x = 1 − x2 e−x

[︀ ]︀

entonces 6A = 1, 6C − 18B − D = 0, 18B − 2C = 0 y −3B = −1

luego A = 1/ 6 , B = 1/ 3, C = 3 y D = 12. entonces
(︃ )︃
x3 x3
yp = + + 3x + 12x e−x
2
6 3

luego la solución general tiene la forma

(︃ )︃
x3 x3
y (x) = c1 + c2 x + c3 x2 + c4 e−x + + + 3x2 + 12x e−x .
6 3

88/95
Ejercicios Propuestos

Resolver los siguientes PVI

1. y 000 + y 00 − 2y = 0 , y (0) = y 0 (0) = 0, y 00 (0) = 1
2. y (4) − 3y 00 + 2y = 0, y (0) = y 0 (0) = y 00 (0) = 0, y 000 (0) = 1
Resolver las siguientes EDO
1. y (5) + 32y = x6 + x4 + x2 + 1
2. y (4) + 4y 00 + 4y = 3cos (2x)
3. y 000 − y 00 + y 0 − y = ex senx + ex
4. y 000 − y 00 + y 0 − y = ex − xcosx

89/95