UNIVERSIDAD AUTÓNOMA DE COAHUILA

FIME – UNIDAD NORTE

PROFESOR: FELIPE MORALES RODRÍGUEZ

Ecuaciones Diferenciales Lineales de Orden Superior.

(A) Ejemplos referentes al tema de las EDO ay’’ + by’ + cy = d(x)
1) y’’+ 4y = 2x + 3
Solución: La EDO homogénea es y’’+ 4y = 0. Se tiene a = 1, b = 0, c = 4, luego b2 – 4ac = –16<0.

b 4ac  b 2
r = 0,   = 2. Por lo tanto
2a 2a
yc = C1 erx cos(θx) + C2 erx sen(θx) = C1 cos(2x) + C2 sen(2x)
Sigue obtener yp con el método de variación de parámetros. En este ejemplo también se puede usar el otro
método (coeficientes indeterminados).
Se tiene y1 = cos(2x) y y2 = sen(2x). Entonces W=y1y2’–y1’y2=cos(2x)[2cos(2x)]–[–2sen(2x)]sen(2x)=2.
d ( x) 2 x  3
Por otro lado, g ( x)  = = 2x+3
a 1
y2 g ( x) y1 g ( x) sen(2 x)(2 x  3) cos(2 x)(2 x  3)
y p   y1  W
dx  y2  W
dx =  cos(2 x) 
2
dx  sen(2 x) 
2
dx

Haciendo las integrales y simplificando se obtiene

1 3
yp  x
2 4
La solución pedida es y = yc + yp, es decir,
1 3
y = C1 cos(2x) + C2 sen(2x) + x
2 4

2) y’’+ 3y’ +2y = cos(ex)

Solución: La EDO homogénea es y’’+ 3y’ + 2y = 0. Se tiene a = 1, b = 3, c = 2, luego b2 – 4ac = 1 > 0.

 b  b 2  4ac
Mediante el uso de r1, 2  se obtiene r1 = – 1 y r2 = – 2. Por lo tanto
2a
yc  C1e 1  C2e = C1e  x  C2e 2 x
r x r2 x

Sigue obtener yp con el método de variación de parámetros. En este ejemplo no se puede usar el otro método
(coeficientes indeterminados).
Se tiene y1 = e – x y y2 = e – 2x. Entonces W = y1y2’–y1’y2 = (e – x)(–2 e – 2x)–(– e – x )(e – 2x) = –e – 3x.

d ( x) cos(e x )
Por otro lado, g ( x)  = = cos(ex)
a 1
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
FIME – UNIDAD NORTE
PROFESOR: FELIPE MORALES RODRÍGUEZ

y2 g ( x) y1 g ( x) (e 2 x )[cos(e x )] x x
 2 x (e )[cos(e )]
y p   y1  W
dx  y2  W
dx =  e  x 
 e 3 x
dx  e   e 3 x dx
Haciendo las integrales y simplificando se obtiene
y p   e 2 x cos(e x )
La solución pedida es y = yc + yp, es decir,
y = C1e  x  C2e 2 x  e 2 x cos(e x )

(B) Ejemplo referente al tema de las EDO ax2y’’ + bxy’ + cy = d(x)

Ejemplo: x2y’’ – xy’ + y = ln(x)
Solución: La EDO homogénea es x2y’’ – xy’ + y = 0. Se tiene a = 1, b = –1, c = 1, luego (b – a)2 – 4ac = 0.
ba
Mediante el uso de r   se obtiene r = 1. Por lo tanto
2a
yc  C1 x r  C2 x r ln( x) = C1 x  C2 x ln( x)
Sigue obtener yp con el método de variación de parámetros. Para este tipo de EDO solamente se puede
emplear este método.
Se tiene y1 = x y y2 = xln(x). Entonces W = y1y2’–y1’y2 = (x)[1 + ln(x)]–[xln(x)](1) = x.
d ( x) ln( x)
Por otro lado, g ( x)  = 2
ax2 x
 ln( x)   ln( x) 
[ x ln( x)] 2  ( x)  2 
y2 g ( x) y1 g ( x)  x  dx  x ln( x)  x 
y p   y1  W
dx  y2  W
dx =  x 
x  x dx
Haciendo las integrales y simplificando se obtiene
y p  ln( x)  2
La solución pedida es y = yc + yp, es decir,
y = C1 x  C2 x ln( x) + ln( x )  2

(C) Ejemplo referente al tema de las EDO any(n) + an–1y(n–1) + . . . + a2y’’ +a1y’+a0y=b(x)
Ejemplo: y’’’ – 4y’ = 3e– x
Solución: Primero se encuentra la solución yc = C1y1 + C2y2 + C3y3 de la EDO y’’’ – 4y’ = 0. Para esto se puede
usar lo indicado en el punto 3 del archivo de apoyo. Se tiene p(r) =r3–4r = r(r2–4) = r(r–2)(r+2). Las raíces y
multiplicidades son r1=0, m1=1, r2=2, m2=1, r3=– 2, m3=1 y por lo tanto yc = C1 + C2e2x + C3e–2x. Se usa lo que
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
FIME – UNIDAD NORTE
PROFESOR: FELIPE MORALES RODRÍGUEZ

acompaña a las constantes para escribir y1=1, y2= e2x, y3= e–2x. Enseguida se obtiene yp mediante el método de
variación de parámetros (solo este método se puede usar, fórmula en el punto 7 del archivo ya indicado). El
desarrollo para yp se muestra enseguida
y1 y2 y3 1 e2 x e 2 x 0 y2 y3 0 e2 x e 2 x
W  y1' y2' y3'  0 2e 2 x 2e 2 x =16, W1  0 y2' y3'  0 2e 2 x 2e 2 x = – 4
y1'' y2'' y3'' 0 4e2 x 4e 2 x 1 y2'' y3'' 1 4e 2 x 4e 2 x

y1 0 y3 1 0 e 2 x y1 y2 0 1 e2 x 0
2 x
W2  y '
1 0 y  0 0 2e
'
3 = 2e –2x
, W3  y '
1 y '
2 0  0 2e 2x
0 = 2e2x.
y ''
1 1 y''
3 0 1 4e 2 x y''
1 y ''
2 1 0 4e 2 x 1

En este caso g(x) = 3e—x (se obtiene dividiendo el lado derecho de la EDO entre lo que acompaña a y’’’).
Luego se usa la fórmula
W1 g ( x) W g ( x) W g ( x)
y p  y1  dx  y2  2 dx  y3  3 dx
W W W
(4)(3e x ) (2e2 x )(3e x ) (2e2 x )(3e x )
y p  (1)  dx  (e ) 
2x
dx  (e ) 
2 x
dx = e – x.
16 16 16
La solución general de la EDO es y = yc + yp = C1 + C2e2x + C3e – 2x + e – x.

(D) Ejemplo referente al tema de las EDO anxny(n) + an–1xn–1y(n–1) + . . . + a2x2y’’ + a1xy’ + a0y = b(x)
Ejemplo: x3y’’’ + xy’ – y = 54x4
Solución: Primero se encuentra la solución yc=C1y1+C2y2 +C3y3 de la EDO x3y’’’+xy’–y=0. Para esto se puede
usar lo indicado en el punto 4 del archivo de apoyo. Se inicia con y=xr y se obtiene y’=rxr – 1, y’’=r(r – 1)xr – 2,
y’’’=r(r – 1)(r – 2)xr – 3. Luego se sustituye en la EDO homogénea para tener x3[r(r–1)(r–2)xr–3]+x(rxr – 1)–xr=0,
así r(r–1)(r–2)+r–1=0, es decir, r3 – 3r2 + 3r – 1 = (r – 1)3 = 0 La raíz y multiplicidad son r1=1, m1=3 y por lo
tanto yc = C1x + C2xln(x) + C3xln2(x).
Se usa lo que acompaña a las constantes para escribir y1=x, y2= xln(x), y3=xln2(x). Enseguida se obtiene
yp mediante el método de variación de parámetros (solo este método se puede usar, fórmula en el punto 8 del
archivo ya indicado). El desarrollo para yp se muestra enseguida

y1 y2 y3 x x ln( x) x ln 2 ( x)
W  y1' y2' y3' = 1 1  ln( x) 2 ln( x)  ln 2 ( x) = 2
y1'' y2'' y3'' 0 1/ x 2 / x  2 ln( x) / x
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
FIME – UNIDAD NORTE
PROFESOR: FELIPE MORALES RODRÍGUEZ

0 y2 y3 0 x ln( x) x ln 2 ( x)
W1  0 y2' y3' = 0 1  ln( x) 2 ln( x)  ln 2 ( x) = xln2(x)
1 y2'' y3'' 1 1/ x 2 / x  2 ln( x) / x

y1 0 y3 x 0 x ln 2 ( x) y1 y2 0 x x ln( x) 0
W2  y1' 0 y3' = 1 0 2 ln( x)  ln 2 ( x) = –2 xln(x), W3  y1' y2' 0 = 1 1  ln( x) 0 = x
y ''
1 1 y ''
3 0 1 2 / x  2 ln( x) / x y''
1 y ''
2 1 0 1/ x 1

54𝑥 4
En este caso 𝑔(𝑥) = = 54𝑥 (se obtiene dividiendo el lado derecho de la EDO entre lo que acompaña a
𝑥3

y’’’). Luego se usa la fórmula

W1 g ( x) W g ( x) W g ( x)
y p  y1  dx  y2  2 dx  y3  3 dx
W W W
[ xln 2 (x)](54x) [2xln(x)](54x) (x)(54x)
Sustituyendo se tiene y p  x  dx + x ln( x)  dx + x ln 2 ( x)  dx
2 2 2
Haciendo las integrales y lo demás se obtiene yp = 2x4. Por lo tanto, la solución general de la EDO es
y = yc + yp = C1x + C2xln(x) + C3xln2(x) + 2x4

(E) Aplicaciones
1) Un cuerpo que pesa 64 libras al sujetarse al extremo de un resorte lo estira 0.32 pie. Si el cuerpo ocupa una
posición que está 8 pulgadas sobre la posición de equilibrio y desde ahí se le comunica una velocidad
dirigida hacia abajo de 5 pie/segundo. (a) Encontrar la ecuación del movimiento; (b) ¿Cuál es la posición,
velocidad y aceleración del cuerpo en t = 3 s?
Solución: Este es un ejemplo correspondiente a los denominados “sistema masa resorte”. Por lo tanto, la
EDO que se usa es my’’ + cy’ + ky =f(t). En este problema, el medio no ofrece resistencia al movimiento
y no se tiene la actuación de fuerza exterior, entonces c = 0 y f(t) = 0. La constante k del resorte se obtiene
con la fórmula k = W/d = 64/0.32 = 200 lb/pie. La masa m se halla con la igualdad m = W/g = 64/32 = 2
slugs.
Puesto que el cuerpo está a 8 pulgadas sobre la posición de equilibrio, se tiene la condición y(0) = –8 pulg.
= –8/12 pie = –2/3 pie. Además, al proporcionársele una velocidad dirigida hacia debajo de 5 pies/seg., se
tiene la condición y’(0) = 5 pies/seg. Así, se necesita resolver el problema
2y’’ + 200y = 0 ó y’’ + 100y = 0, y(0) = –2/3, y’(0) = 5
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
FIME – UNIDAD NORTE
PROFESOR: FELIPE MORALES RODRÍGUEZ

La EDO tiene la presentación ay’’ + by’ + cy = 0 con a = 1, b = 0, c = 100. Se obtiene b2 – 4ac = – 400<0.
𝑏 √4𝑎𝑐−𝑏 2
Se usan las fórmulas 𝑟 = − 2𝑎, 𝜃 = para obtener r = 0 y θ = 10. Por lo tanto, la solución general
|2𝑎|

es
y = C1 ert cos(θt) + C2 ert sen(θt) = C1cos(10t) + C2sen(10t)
De esta solución se puede encontrar
y’ = –10C1sen(10t) + 10C2cos(10t)
Al sustituir la condición y(0) = –2/3, se escribe la igualdad –2/3 = C1cos(0) + C2sen(0). Luego C1= –2/3.
Al sustituir la condición y’(0) = 5, se escribe la igualdad 5 = –10C1sen(0) + 10C2cos(0). Luego C2= 0.5
Entonces: (a) La ecuación de movimiento es y = – (2/3) cos(10t) + 0.5 sen(10t); (b) Si t = 3, posición = y
= – (2/3) cos(30) + 0.5 sen(30) = – 0.59685 pie (el signo negativo indica que el cuerpo está por encima
de la posición de equilibrio), velocidad = y’ = (20/3) sen(10t) + 5 cos(10t) = (20/3) sen(30) + 5 cos(30) =
–5.81562 pies/seg (signo negativo indica cuerpo va hacia arriba), aceleración = y’’ = (200/3) cos(10t) –
50 sen(10t) = 59.685 pies/seg2 .

2) Un objeto que pesa 1/4 de libra, se deja caer desde lo alto de un edificio. Determinar el desplazamiento
en términos del tiempo si la resistencia del aire es fr = v/160. Calcular la distancia recorrida cuando han
transcurrido 4 segundos.
Solución: El movimiento es vertical, luego se usa la ecuación my’’ = ∑Fy . Las fuerzas que actúan sobre
el cuerpo son el peso y la resistencia del aire, entonces ∑Fy = W – fr = 1/4 – v/160. Por otro lado m =
W/g = (1/4)/32=1/128 slug. Por lo tanto, la EDO a considerar es (1/128)y’’ = 1/4 – v/160 = 1/4 – y’/160
ó sea y’’ = 32 – (4/5)y’, es decir, y’’ + (4/5)y’ = 32.
Las condiciones son: y(0) = 0 (la distancia se empieza a medir cuando se deja caer el objeto) y y’(0) = 0
(el objeto se deja caer y entonces la velocidad inicial es cero). El problema es
y’’ + (4/5)y’ = 32, y(0) = 0, y’(0) = 0
Nota: La variable independiente es t (tiempo). En el desarrollo siguiente viene siendo la x, es decir, x
indica tiempo.
La EDO tiene la presentación ay’’ + by’ + cy =d(x). La EDO homogénea es y’’+ (4/5)y’ = 0. Se tiene a=1,

 b  b 2  4ac
b = 4/5, c = 0, luego b2 – 4ac = 16/25 > 0. Mediante el uso de r1, 2  se obtiene r1 = 0 y
2a
r2 = – 4/5. Por lo tanto
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
FIME – UNIDAD NORTE
PROFESOR: FELIPE MORALES RODRÍGUEZ

4
 x
yc  C1e 1  C2e = C1  C2e
r x r2 x 5

Sigue obtener yp con el método de variación de parámetros (también se puede usar el método de
coeficientes indeterminados).
4 4 4 4
4 4
Se tiene y1 = 1 y y2 = 𝑒 − 5 𝑥 . Entonces W = y1y2’–y1’y2 = (1)( − 5 𝑒 − 5 𝑥 .) – (0 )( 𝑒 − 5 𝑥 .) = − 5 𝑒 − 5 𝑥 .

d ( x) 32
Por otro lado, g ( x)  = = 32
a 1
4
 x
y2 g ( x) y1 g ( x)
 
4
y p   y1 dx  y2 dx =  (1) (e )(32) dx  e  5 x (1)(32) dx
5

W W  4  4x  4  4x
5
 e 5
 e
5 5

Haciendo las integrales y simplificando se obtiene

yp = 40x – 50
La solución general es y = yc + yp, es decir,
4 4
 x  x
y = C1  C2e 5
+ 40x – 50 = K1  C2e 5
+ 40x con K1 = C1 – 50
De esta solución se puede encontrar
4
4  x
y’ =  C2 e 5 + 40
5
4
 ( 0)
Al sustituir la condición y(0) = 0, se escribe la igualdad 0 = K1  C2e 5
+ 40(0). Luego, se forma la
ecuación K1 + C2 = 0.
4
4  (0)
Al sustituir la condición y’(0) = 0, se escribe la igualdad 0 =  C2 e 5 + 40 . Luego C2 = 50. Al
5
sustituir este valor de C2 en la ecuación K1 + C2 = 0 se obtiene K1 = –50. Por lo tanto
4
 x
y =  50  50e 5
+ 40x
El desplazamiento en términos del tiempo viene dado por
4
 t
y =  50  50e 5
+ 40t
La distancia recorrida, cuando han transcurrido 4 segundos, es
4
 ( 4)
y =  50  50e 5
+ 40(4) = 112.03811 pies
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
FIME – UNIDAD NORTE
PROFESOR: FELIPE MORALES RODRÍGUEZ

¿Cómo se responde lo siguiente? Calcular el tiempo que tarda en llegar al suelo si el edificio tiene 120
pies de altura.

3) Find the charge on the capacitor in an RLC series circuit when L=0.5 h, R=10 Ω, C=0.01 f, E(t)=150 V,
q(0)=1 C, and i(0) = 0 A. What is the charge on the capacitor after a long time?
1
Solución: La EDO a usar es Lq ' ' Rq ' q  E (t ) , es decir, 0.5q’’ + 10q’ + 100q = 150, equivalente a
C
q’’ + 20q’ + 200q = 300.
Las condiciones son: q(0) = 1 y q’(0) = 0. El problema es
q’’ + 20q’ + 200q = 300, q(0) = 1, i(0) = q’(0) = 0
Nota: La variable independiente es t (tiempo) y la dependiente es q (carga). En el desarrollo siguiente, el
tiempo se representa por x y la carga por la variable y.
La EDO tiene la presentación ay’’ + by’ + cy =d(x). La EDO homogénea es y’’ + 20y’ + 200y = 0. Se
𝑏 √4𝑎𝑐−𝑏 2
tiene a=1, b = 20, c = 200, luego b2 – 4ac = – 400 < 0. Se usan las fórmulas 𝑟 = − 2𝑎, 𝜃 = para
|2𝑎|

obtener r = – 10 y θ = 10. Por lo tanto, la solución yc es

yc = C1 erx cos(θx) + C2 e rx sen(θx) = C1 e – 10x cos(10x) + C2 e – 10x sen(10x)
Sigue obtener yp. Vamos a usar el método de coeficientes indeterminados
Se tiene y1=e –10xcos(10x) y y2=e –10xsen(10x). El lado derecho de la EDO es 300, entonces α = 0, β = 0,
n = 0. Sustituyendo estos valores en la fórmula indicada por el método
yp=eαx(Anxn+An–1xn–1+...+A2x2+A1x+A0)cos(βx) + eαx(Bnxn+Bn–1xn–1+...+B2x2+B1x+B0)sen(βx)
se obtiene yp = A0. Sigue analizar si la yp se modifica. Como β = 0 y las funciones eαx = 1, xeαx = x no son soluciones
de la homogenea (no son ni y1 ni y2), la yp no cambia (inciso (c) del teorema de la página 33 del archivo de los
apuntes). Por lo tanto yp = A0 , yp’ = 0, yp’’ = 0. Sustituyendo en la EDO y’’ + 20y’ + 200y = 300 se obtiene 0 +
0 + 200A0 = 300, es decir, A0 = 300/200 = 3/2.
Así, yp = 3/2, y = yc + yp = C1 e – 10x cos(10x) + C2 e – 10x sen(10x) + 3/2. De aquí se deducen las igualdades
q = C1 e – 10t cos(10t) + C2 e – 10t sen(10t) + 3/2
q = e – 10t [C1 cos(10t) + C2 sen(10t)] + 3/2
q’ = e – 10t [–10C1 sen(10t) + 10C2 cos(10t)] + [C1 cos(10t) + C2 sen(10t)](– 10 e – 10t)
Al sustituir la condición q(0) = 1, se escribe la igualdad 1 = (1)(C1 + 0) +3/2. Luego, C1 = – 1/2.
Al sustituir la condición q’(0) = 0, se escribe la igualdad 0 = (1)(0 + 10C2)+ (C1 + 0)(–10). Entonces, 0 =
10C2 – 10C1. De esto se obtiene C2 = C1 = – 1/2. Por lo tanto, la carga en el capacitor está dada por
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
FIME – UNIDAD NORTE
PROFESOR: FELIPE MORALES RODRÍGUEZ

q = e – 10t [– 0.5 cos(10t) – 0.5 sen(10t)] + 3/2

En la forma amplitud – fase se tiene
q = 0.7071067812 e – 10t cos(10t + 2.35619449) + 3/2
Puesto que: (a) las funciones seno y coseno están acotadas por los números – 1 y 1, (b) la función
exponencial tiende a cero cuando t tiende al infinito; la carga en el capacitor después de un tiempo largo
(tiempo tendiendo al infinito) es
lim 𝑞 = 3/2
𝑡→∞