106 Ampliación de Matemáticas

4.3.3. Ejercicios.

Determine la solución general de cada ecuación diferencial:

1. 4y 00 + y 0 = 0 7. y 000 − y = 0
2. y 00 + 9y = 0 8. 12y 00 − 5y 0 − 2y = 0
d4 y d3 y d2 y
3. y 00 − 36y = 0 9. dx4
+ dx 3 + dx2 = 0
˜
4. y 00 − y 0 − 6y = 0 10. y 00 + 4y 0 − y = 0
5. d2 y dy
+ 8 dx + 16y = 0 11. y 00 − 8y = 0
dx2
6. y − 4y − 5y 0 = 0
000 00 12. y 000 − 5y 00 + 3y 0 + 9y = 0
Resuelva cada ecuación diferencial, sujeta a las condiciones iniciales
indicadas.

13. y 00 + 16y = 0, y(0) = 2, y 0 (0) = −2

14. y 00 − y = 0, y(0) = 1, y 0 (0) = 1

15. y 00 + 6y 0 + 5y = 0, y(0) = 0, y 0 (0) = 3

16. y 00 + y 0 + 2y = 0, y(0) = 0, y 0 (0) = 0

17. y 00 + y = 0, y( π3 ) = 0, y 0 ( π3 ) = 2

18. y 000 + 12y 00 + 36y 0 = 0, y(0) = 0, y 0 (0) = 1, y 00 (0) = −7

4.4. Coeficientes indeterminados, método del anu-

lador.
Recuerda 4.4.1 Teniendo en cuenta la definición 4.1.4 de operador dife-
rencial, una ecuación lineal de coeficientes constantes

an y (n) + . . . + a0 y = g(x)

puede escribirse como

an Dn y + an−1 Dn−1 y + · · · + a1 Dy + a0 y = g(x)

,
(a Dn + an−1 Dn−1 + · · · + a1 D + a0 ) y = g(x)
| n {z }
P (D)
Antonio Baeza Salas 107

A P (D) se le denomina operador diferencial lineal de orden n. Este poli-

nomio podrá ser factorizado generalmente y los factores pueden conmutarse.
Ejemplo 4.4.1
y 00 − y = ex
puede escribirse
D 2 y − y = ex
(D2 − 1)y = ex
(D − 1)(D + 1)y = ex
•
Ejemplo 4.4.2 Una ecuación diferencial comoy 00 + 4y 0 + 4y = 0 se puede
escribir en la forma
(D2 + 4D + 4)y = 0
o sea
(D + 2)(D + 2)y = 0
o sea
(D + 2)2 y = 0
•

Nos interesa transformar una ecuación no homogénea en otra que si

lo sea. Para ello buscamos operadores que anulen el segundo miembro, la
función g(x). Estos reciben el nombre de operadores anuladores.

Definición 4.4.1 (Operador anulador.) Si L es un operador diferencial

con coeficientes constantes y f es una función suficientemente diferenciables
tal que
L(f (x)) = 0
se dice que L es un anulador de la función.

Veremos distintos tipos de operadores anuladores según la forma de la fun-

ción g; por ejemplo, la función constantes y = k es anulada por D porque
Dk = 0. La función y = x es anulada por el operador diferencial D 2 porque
la primera y la segunda derivada de x son 1 y 0, respectivamente. En forma
similar, D 3 x2 = 0, etcétera.
Recuerda 4.4.2 El operador diferencial D n anula cada una de las si-
guientes funciones:
1, x, x2 , x3 , . . . , xn−1
108 Ampliación de Matemáticas

Un polinomio se puede anular definiendo un operador que anule la po-

tencia máxima de x.
Recuerda 4.4.3 El operador diferencial (D − α)n anula cada una de las
siguientes funciones:
eαx , xeαx , . . . , xn−1 eαx

Recuerda 4.4.4 El operador diferencial [D 2 − 2αD + (α2 + β 2 )]n anula

cada una de las siguientes funciones:
eαx cos βx, xeαx cos βx, . . . , xn−1 eαx cos βx
eαx sen βx, xeαx sen βx, . . . , xn−1 eαx sen βx

Recuerda 4.4.5 Supongamos que L(y) = g(x) es una ecuación diferencial

lineal con coeficientes constantes, y que la entrada g(x) consiste en sumas y
productos finitos de las funciones mencionadas anteriormente, esto es, que
g(x) es una combinación lineal de las funciones
k, xm , xm eαx , xm eαx cos βx, xm eαx sen βx
en donde m es entero no negativo y α y β son números reales. Ya sabemos
que esa función g se puede anular con un operador diferencial, L1 , de orden
mı́nimo, formado por un producto de los operadores D n , (D − α)n y [D2 −
2αD + (α2 + β 2 )]n . Aplicamos L1 a ambos lados de la ecuación L(y) = g(x)
y obtenemos L1 L(y) = L1 (g(x)) = 0 al resolver la ecuación homogénea
y de orden superior L1 L(y) = L1 (g(x)) = 0, descubrimos la forma de una
solución particular, yp , de la ecuación original no homogénea L(y) = g(x). A
continuación sustituimos esa forma supuesta en L(y) = g(x) para determinar
una solución particular explı́cita. Este procedimiento de determinación de
yp se llama método de los coeficientes indeterminados.
Recordemos que la solución general de una ecuación diferencia lineal no
homogénea L(y) = g(x) es y = yc +yp donde yc es la función complementaria,
esto es, la solución general de la ecuación homogénea asociada L(y) = 0. La
solución general de cada ecuación L(y) = g(x) esta definida en el intervalo
(−∞, +∞).

Ejemplo 4.4.3 Resuelve

y 00 + 3y 0 + 2y = 4x2 (4.23)
Solución:
Antonio Baeza Salas 109

Paso 1 Primero resolvemos la ecuación homogénea y 00 + 3y 0 + 2y = 0. A

partir de la ecuación auxiliar m2 + 3m + 2 = (m + 1)(m + 2) = 0, se
tiene que la función complementaria es

yc = c1 e−x + c2 e−2x

Paso 2 Como el operador diferencial D 3 anula a 4x2 , vemos que D 3 (D2 +

3D + 2)y = 4D 3 x2 es lo mismo que

D3 (D2 + 3D + 2)y = 0

La ecuación auxiliar de la ecuación anterior, de quinto orden es

m3 (m2 + 3m + 2) = 0

m3 (m + 1)(m + 2) = 0
tiene por raı́ces m1 = m2 = m3 = 0, m4 = −1, m5 = −2. Ası́ la
solución general debe ser

y = c1 + c2 x + c3 x2 + c4 e−x + c5 e−2x

Los términos cuarto y quinto constituyen la función complementaria de

la ecuación original (4.23). Podemos decir que una solución particular,
yp , también deberı́a satisfacer la ecuación

D3 (D2 + 3D + 2)y = 0

. Esto significa que los términos restantes en la ecuación

y = c1 + c2 x + c3 x2 + c4 e−x + c5 e−2x

han de tener la forma básica de yp :

y = A + Bx + Cx2 (4.24)

Para que la ecuación (4.24) sea una solución particular de la (4.23), se

necesita determinar los coeficientes especı́ficos A, B y C. Derivamos
la función (4.24) para obtener

yp0 = B + 2Cx, yp00 = 2C

110 Ampliación de Matemáticas

y sustituidos en (4.23) para llegar a

yp00 + 3yp0 + 2yp = 2C + 3B + 6Cx + 2A + 2Bx + 2Cx2 = 4x2

Como se supone que esta última ecuación tiene que er una identidad,
los coeficientes de las potencias de igual grado en x deben ser iguales.
Esto es:
2C = 4, 2B + 6C = 0, 2A + 3B + 2C = 0
de donde se tiene que

yp = 7 − 6x + 2x2

Paso 3 La solución general de la ecuación es y = yc + yp

y = c1 e−x + c2 e−2x + 7 − 6x + 2x2

•
Ejemplo 4.4.4 Resuelve

y 00 − 3y 0 = 8e3x + 4 sen x (4.25)

Solución:

Paso 1 Primero resolvemos la ecuación homogénea. A partir de la ecuación

auxiliar m2 − 3m = m(m − 3) = 0, se tiene que la función complemen-
taria es
yc = c1 + c2 e3x
.

Paso 2 En vista de que (D − 3)e3x = 0 y (D 2 + 1) sen x = 0, aplicamos el

operador diferencial (D − 3)(D 2 + 1) a ambos lados de (4.25):

(D − 3)(D 2 + 1)(D 2 − 3D)y = 0

La ecuación auxiliar de esta ecuación es

(m − 3)(m2 + 1)(m2 − 3m) = 0

m(m − 3)2 (m2 + 1) = 0

De modo que

y = c1 + c2 e3x + c3 xe3x + c4 cos x + c5 sen x

Antonio Baeza Salas 111

Después de excluir la combinación lineal de términos c1 + c2 e3x que

corresponde a yc , llegamos a la forma de yp :

yp = Axe3x + B cos x + C sen x

Sustituimos yp en (4.25), simplificamos y obtenemos

yp00 − 3yp0 = 3Ae3x + (−B − 3C) cos x + (3B − C) sen x = 8e3x + 4 sen x

Igualamos coeficientes:

 3A = 8
−B − 3C = 0

3B − C = 4

Vemos que A = 8/3, B = 6/5, C = −2/5 y en consecuencia

8 6 2
yp = xe3x + cos x − sen x
3 5 5

Paso 3 La solución general es:

8 6 2
y = c1 + c2 e3x + xe3x + cos x − sen x
3 5 5
•
Ejemplo 4.4.5 Resuelve

y 00 + y = x cos x − cos x (4.26)

Solución:
Paso 1 Primero resolvemos la ecuación homogénea. A partir de la ecua-
ción auxiliar m2 + 1 = (m − i)(m + i) = 0, se tiene que la función
complementaria es
yc = c1 cos x + c2 sen x
.

Paso 2 (D 2 + 1)2 es un anulador del lado derecho de la ecuación (4.26):

(D2 + 1)(D 2 + 1)y = 0 (D 2 + 1)3 y = 0

Como i y −i son raı́ces complejas de multiplicidad 3 de la ecuación

auxiliar de la última ecuación diferencial, se tiene

y = c1 cos x + c2 sen x + c3 x cos x + c4 x sen x + c5 x2 cos x + c6 x2 sen x

112 Ampliación de Matemáticas

Después de excluir la combinación lineal de términos c1 cos x + c2 sen x

que corresponde a yc , llegamos a la forma de yp :

yp = Ax cos x + Bx sen x + Cx2 cos x + Ex2 sen x

Sustituimos yp en (4.26), simplificamos y obtenemos

yp00 +yp = 4Ex cos x−4Cx sen x+(2B+2C) cos x+(−2A+2E) sen x = x cos x−cos x

Igualamos coeficientes:


 4E = 1

−4C = 0
 2B + 2C = −1


−2A + 2E = 0

Vemos que A = 1/4, B = −1/2, C = 0, E = 1/4 y en consecuencia

Paso 3 La solución general es:

1 1 1
y = c1 cos x + c2 sen x + x cos x − x sen x + x2 sen x
4 2 4

•
Recuerda 4.4.6 [Coeficientes indeterminados, método del anula-
dor.] La ecuación diferencial L(y) = g(x) tiene coeficientes constantes y la
función g(x) consiste en sumas y productos finitos de constantes, polinomios,
funciones exponenciales eax , senos y cosenos.

1. Se determina la solución complementaria, yc , de la ecuación homogénea

L(y) = 0.

2. Ambos lados de la ecuación no homogénea L(y) = g(x) se someten a

la acción de un operador diferencial, L1 , que anule la función g(x).

3. Se determina la solución general de la ecuación diferencial homogénea

de orden superior L1 L(y) = 0.

4. De la solución obtenida en el paso anterior, se eliminan todos los térmi-

nos duplicados en la solución complementaria que se determino en el
paso 1). Se forma una combinación lineal, yp , con los términos restan-
tes. esta será la forma de una solución particular de L(y) = g(x).
Antonio Baeza Salas 113

5. Se sustituye yp que se determinó en el paso 4) en L(y) = g(x). Se

igualan los coeficientes de las diversas funciones a cada lado de la
igualdad y se despejan los coeficientes desconocidos en yp del sistema
de ecuaciones resultante.

6. Con la solución particular que se determinó en el paso anterior, se

forma la solución general y = yc + yp de la ecuación diferencial dada.

4.4.1. Ejercicios.
1. En los problemas que siguen escriba la ecuación diferencial dada en
la forma L(y) = g(x), donde L es un operador diferencial lineal con
coeficientes constantes. Si es posible, factorice L.

a) 9y 00 − 4y = sen x
b) y 00 − 4y 0 − 12y = x − 6
c) y 000 + 10y 00 + 25y 0 = ex

2. Determine un operador lineal que anule la función dada:

a) 1 + 6x − 2x3
b) 1 + 7e2x
c) cos 2x

3. Resuelva la respectiva ecuación diferencial por el método de los coefi-

cientes indeterminados:

a) y 00 − 9y = 54
b) y 00 + 4y 0 + 4y = 2x + 6
c) y 00 − y 0 − 12y = e4x
d ) y 00 + 25y = 6 sen x

4.5. Variación de parámetros.

Se trata de un método por el que obtenemos una solución particular de
la no homogénea.
Se resuelve la ecuación homogénea y se impone que la solución particular
que busco tiene la misma forma cambiando las constantes por funciones que
si pueden depender de x.
114 Ampliación de Matemáticas

Para aplicar el método de variación de parámetros a una ecuación dife-

rencial de segundo orden,

a2 (x)y 00 + a1 (x)y 0 + a0 (x)y = g(x) (4.27)

llevamos la ecuación diferencial a su forma reducida

y 00 + P (x)y 0 + Q(x)y = f (x) (4.28)

Suponemos que P, Q, f son continuas en algún intervalo I. Para la ecuación

lineal de segundo orden (4.27) se busca una solución de la forma

yp = u1 (x)y1 (x) + u2 (x)y2 (x) (4.29)

en que y1 y y2 formen un conjunto fundamental de soluciones, en I, de

la forma homogénea asociada de (4.27). Aplicamos dos veces la regla del
producto para diferenciar yp y obtenemos:

yp0 = u1 (x)y10 (x) + u01 (x)y1 (x) + u2 (x)y20 (x) + u02 (x)y2 (x)
yp00 = u1 (x)y100 (x) + u01 (x)y10 (x) + u001 (x)y1 (x) + u01 (x)y10 (x)+
u2 (x)y200 (x) + u02 (x)y20 (x) + u002 (x)y2 (x) + u02 (x)y20 (x)

Sustituimos (4.29)en la ecuación (4.27)

yp00 + P (x)yp0 + Q(x)yp = f (x)

cero cero
z }| { z }| {
u1 (y1 + P (x)y1 + Q(x)y1 ) +u2 (y2 + P (x)y20 + Q(x)y2 )
00 0 00

+y1 u001 + u01 y10 + y2 u002 + u02 y20 + P [y1 u01 + y2 u02 ] + y10 u01 + y20 u02
= f (x)
d 0 d 0 0 0 0 0 0 0
dx [y1 u1 ] + dx [y2 u2 ] + P [y1 u1 + y2 u2 ] + y1 u1 + y2 u2
= f (x)
d 0 0 0 0 0 0 0 0
dx [y1 u1 + y2 u2 ] + P [y1 u1 + y2 u2 ] + y1 u1 + y2 u2
= f (x)
(4.30)
Dado que buscamos determinar dos funciones desconocidas, u1 y u2 , es
de esperar que necesitemos dos ecuaciones. Las podemos obtener si estable-
cemos la hipótesis adicional de que las funciones u1 y u2 satisfacen y1 u01 +
y2 u02 = 0. Esta hipótesis es pertinente porque si pedimos que y1 u01 +y2 u02 = 0
la ecuación (4.30) se reduce a y10 u01 + y20 u02 = f (x). Con ello ya tenemos las
dos ecuaciones que deseábamos, aunque sea para determinar las derivadas
u01 y u02 . Aplicamos la regla de Cramer y la solución del sistema

y1 u01 + y2 u02 = 0
y10 u01 + y20 u02 = f (x)
Antonio Baeza Salas 115

se puede expresar en términos de los determinates

W1 W2
u01 = y u02 = (4.31)
W W
en donde
     
 y1 y2   0 y 2
  y 0 
W =  0  W1 =   W2 =  10  (4.32)
y1 y20  0
f (x) y2  y1 f (x) 

Las funciones u1 y u2 se determinan integrando los resultados en (4.31). Se

que el determinante W es el wronskiano de y1 y y2 . Sabemos que por la
independencia de y1 y y2 en I, W 6= 0 para todo x en el intervalo I.
Recuerda 4.5.1 Para resolver (4.27),

Paso 1 Se halla la función complementaria yc = c1 y1 + c2 y2 ,

Paso 2 Se calcula el wronskiano W (y1 , y2 ).

Paso 3 Se lleva la ecuación a su forma reducida (4.28) para hallar f (x).

W1 W2
Paso 4 Se determinan u1 y u2 integrando u01 = W y u02 = W .

Paso 5 Una solución particular es yp = u1 y1 + u2 y2 .

Paso 6 La solución general de la ecuación es y = yc + yp

Ejemplo 4.5.1 Resuelve la ecuación y 00 − 4y 0 + 4y = (x + 1)e2x

Solución: Partimos de la ecuación auxiliar m2 − 4m + 4 = (m − 2)2 = 0,
y tenemos que yc = c1 e2x + c2 xe2x . Identificamos y1 = e2x y y2 = xe2x y
calculamos el wronskiano
 2x 
 e xe2x 
 = e4x
W =  2x 2x

2x
2e e + 2xe 

Como la ecuación diferencial dada está en la forma reducida, vemos que

f (x) = (x + 1)e2x Ası́ tenemos que:
 2x
  2x 
 0 xe  4x
 e 0 
 = (x+1)e4x
W1 =  2x 2x 2x
 = −(x+1)xe W2 =  2x
(x + 1)e 2xe + e  2e (x + 1)e2x 

de (4.31) se tiene que

u01 = −x2 − x u02 = x + 1

116 Ampliación de Matemáticas

, por tanto
u1 = −x3 /3 − x2 /2 u2 = x2 /2 + x
3 x2 2x
y se tiene yp = ( x6 + 2 )e y

x3 x2 2x
y = yc + yp = c1 e2x + c2 xe2x + ( + )e
6 2
•

4.5.1. Ejercicios.
Resolver cada una de las ecuaciones diferenciales por variación de paráme-
tros:

1. y 00 + y = sen x
1
2. y 00 + 3y 0 + 2y = 1+ex

e2x
3. y 00 − 4y = x

4. y 000 + y 00 = 8x2
1
5. y 00 − y = x

4.6. Ecuación de Cauchy-Euler.

Toda ecuación diferencial lineal de la forma
dn y n−1 d
n−1 y dy
an xn + a n−1 x + a1 x + a0 y = g(x)
dxn dxn−1 dx
donde los coeficientes an , · · · , a0 son constantes, tiene los nombres de ecua-
ción de Cauchy-Euler, ecuación de Euler Cauchy, ecuación de Euler o ecua-
ción equidimensional.
Examinaremos el caso de la ecuación homogénea de segundo orden:

d2 y dy
ax2 2
+ bx + cy = 0
dx dx
La solución de ecuaciones orden superior será análoga. Una vez deter-
minada la función complementaria yc (x) también podemos resolver la ecua-
ción no homogénea ax2 y 00 + bxy 0 + cy = g(x) con el método de variación de
parámetros.
Antonio Baeza Salas 117

4.6.1. Método de solución.

Intentaremos una solución de la forma y = xm , donde m está por deter-
minar. La primera y segunda derivada son,

dy d2 y
= mxm−1 = m(m − 1)xm−2
dx dx2
En consecuencia:
2
d y dy
ax2 dx 2 + bx dx + cy = ax2 m()m − 1)xm−2 + bxmxm−1 + cxm
= am(m − 1)xm + bmxm + cxm
= xm (am(m − 1) + bm + c)

Ası́, y = xm es una solución de la ecuación diferencial siempre que m sea

una solución de la ecuación auxiliar

am(m − 1) + bm + c = 0 ó am2 + (b − a)m + c = 0 (4.33)

Hay tres casos distintos por considerar que dependen de si las raı́ces de
esta ecuación cuadrática son reales y distintas, reales e iguales o complejas.

CASO 1: raı́ces reales distintas.

Sean m1 y m2 las raı́ces reales de (4.33), tales que m1 6= m2 . Enton-
ces y1 = xm1 y y2 = xm2 forman un conjunto fundamental de soluciones.
Ası́ pues, la solución general es

y = c 1 x m1 + c 2 x m2 (4.34)

CASO 2: raı́ces reales iguales.

Si las raı́ces de (4.33) son repetidas sólo llegaremos a una solución, que es
y = xm1 . Cuando las raı́ces de la ecuación cuadrática am2 + (b − a)m + c = 0
son iguales, el discriminante de los coeficientes tiene que ser cero. De acuerdo
con la fórmula cuadrática, la raı́z debe ser m1 = −(b − a)/2a.
Podemos formar ahora una segunda solución, y2 , aplicando una reduc-
ción de orden. Primero escribimos la ecuación de Cauchy-Euler en la forma:

d2 y b dy c
2
+ + 2y = 0
dx ax dx ax
b
R b
tenemos por tanto P (x) = ax e ( ax )dx = ab ln x.
Por tanto
118 Ampliación de Matemáticas

R − ab ln x
y2 = xm1 e x2m dx
m
R − b 1 −2m
= x 1 x ax 1 dx
R b b−a
= xm1 R x− a x a dx
= xm1 dx x =x
m1 ln x

La solución general es:

y = c1 xm1 + c2 xm1 ln x (4.35)

CASO 3: raı́ces complejas.

Si las raı́ces de (4.33) son el par conjugado m1 = α + iβ, m2 = α − iβ,
la solución general es

y = xα (c1 cos(β ln x) + c2 sen(β ln x)) (4.36)

4.6.2. Ejercicios
Resuelva las siguientes ecuaciones diferenciales:
1. x2 y 00 − 2y = 0

2. x2 y 00 + xy 0 + 4y = 0

3. x2 y 00 − 3xy 0 − 2y = 0

4. x2 y 00 + 6xy 0 + y = 0
Resuelve cada ecuación diferencial sujeta a las condiciones iniciales
indicadas.

5. x2 y 00 + 3xy 0 = 0, y(1) = 0, y 0 (1) = 4

6. x2 y 00 + xy 0 + y = 0, y(1) = 1, y 0 (1) = 2

4.7. Modelado con ecuaciones diferenciales de or-

den superior.
4.7.1. Sistemas de resorte y masa: movimiento libre no amor-
tiguado.
Ley de Hooke. Supongamos que, como en la figura 4.1(b), una masa
m1 esta unida a un resorte flexible colgado de un soporte rı́gido. Cuando