Operador Diferencial y Ecuaciones Diferenciales

1. Operador Diferencial
Un operador es un objeto matemático que convierte una funci´on en otra, por ejemplo, el opera-
dor derivada convierte una funci´on en una función diferente llamada la función derivada. Podemos
definir el operador derivada D que al actuar sobre una funci´on diferenciable produce la derivada
de esta, esto es:

D0 f(x) = f(x) , D1 f(x) = f ′ (x) , D2 f(x) = f ′′ (x) , . . . , Dn f(x) = f (n) (x) .

Es posible construir la siguiente combinación lineal con los operadores diferenciales:

P (D) = a0 + a1 D + a2 D2 + · · · + an Dn , an 6= 0 . (1)

donde a2 , a1 , a2 , . . . an son constantes. A este nuevo objeto lo podemos llamar el Operador Polino-
mial de orden n.

La utilidad de este objeto matemático quedar´a clara si hacemos la siguiente definici´on

P (D)y ≡ an Dn + an−1 Dn−1 + · · · + a2 D2 + a1 D + a0 y
= an Dn y + an−1 Dn−1 y + · · · + a2 D2 y + a1 Dy + a0 y
= an y(n) + an−1 y(n−1) + · · · + a2 y′′ + a1 y′ + a0 y (2)

Por otro lado, recordemos que una ecuaci´on diferencial lineal de orden n con coeficientes cons-
tantes es una ecuaci´on de la forma

an y(n) + an−1 y(n−1) + · · · + a2 y′′ + a1 y′ + a0 y = Q(x) , (3)

por lo tanto, (3) se puede escribir de una manera compacta como

P (D)y = Q(x) . (4)

El operador polinomial es lineal, esto significa que tiene las siguientes propiedades

Si f1 (x) y f2 (x) son dos funciones diferenciables de orden n, entones

P (D) [αf1 (x) + βf2 (x)] = αP (D)f1 (x) + βP (D)f2 (x)

donde α y β son constantes. Además:

Si y1 (x), y2 (x), . . . , yn (x) son n soluciones de la ecuación diferencial homogénea P (D)y = 0

entonces yh (x) = C1 y1 (x) + C2 y2 (x) + · · · + Cn yn (x) es también una solución.

Si yh (x) es una solución de P (D)y = 0 y yp (x) es una solución de P (D)y = Q(x) entonces
y(x) = yh (x) + yp (x) es una solución de P (D)y = Q(x).

Si yp1 (x), yp2 (x), . . . , ypn (x) son soluciones particulares de las respectivas n ecuaciones

P (D)y = Q1 (x), P (D)y = Q2 (x), . . . , P (D)y = Qn (x)

resulta entoneces que

P (D) [yp1 (x) + yp2 (x) + · · · + ypn (x)] = Q1 (x) + Q2 (x) + · · · + Qn (x)

implica que
yp (x) = yp1 (x) + yp2 (x) + · · · + ypn (x)
es una solución de
P (D)y = Q1 (x) + Q2 (x) + · · · + Qn (x)

Ejemplo Encuentre una ecuaci´on particular de

y′′ + y = x2 + xe2x + 3 ⇒ D2 + 1 y = x2 + xe2x + 3

Por alguno de los m´etodos anteriormente vistos podemos encontrar las soluciones particulares de
cada una de las siguientes ecuaciones




D2 + 1 y = x2 , D2 + 1 y = xe2x , D2 + 1 y = 3

las soluciones son, respectivamente

  
1 4
yp1 (x) = x2 − 2 , yp2 (x) = xe2x − e2x , yp3 (x) = 3
5 5

por lo tanto, una soluci´on particular de la ecuaci´on diferencial problema es

   
2 1 2x 4 2x 2 1 2x 4 2x
yp (x) = x − 2 + xe − e +3=x +1+ xe − e .
5 5 5 5

Al operador diferencial también se le pueden agregar las siguientes propiedades. Consideremos

los operadores

P1 (D) = an Dn + an−1 Dn−1 + · · · + a2 D2 + a1 D + a0 (5)

n n−1 2
P2 (D) = bn D + bn−1 D + · · · + b2 D + b 1 D + b 0 (6)

con n ≥ m.
 
(P1 +P2 )f(x) = P1 f(x)+P2 f(x) = an Dn + · · · + an Dn + · · · + (a2 + b2 )D2 + (a1 + b1 )D + a0 + b0 f(x)

[g(x)P (D)] f(x) = g(x) [P (D)f(x)]

g(x) [P1 (D) + P2 (D)] f(x) = g(x) [P1 (D)f(x)] + g(x) [P2 (D)f(x)]

[P1 (D)P2 (D)] f(x) = P1 (D) [P2 (D)f(x)]

Si el operador polinomial diferencial puede ser factorizado, es decir, si:

P (D) = an Dn + an−1 Dn−1 + · · · + a2 D2 + a1 D + a0 = an (D − m1 )(D − m2 ) · · · (D − mn ) ,

las cantidades: m1 , m2 , . . . , mn son las raices de la ecuación caracterı́stica de P (D)y = 0.

Si
P (D) = an Dn + an−1 Dn−1 + · · · + a2 D2 + a1 D + a0 ,
es un operador polinomial de orden n, entonces:

P (D + a) = an (D + a)n + an−1 (D + a)n−1 + · · · + a2 (D + a)2 + a1 (D + a) + a0

donde a es una constante.

Ejemplo
   
(D − 2D − 3) sen(x) + x2 = (D + 1)(D − 3) sen(x) + x2
  
= (D + 1) (D − 3) sen(x) + x2
 
= (D + 1) cos(x) + 2x − 3sen(x) − 3x2
= −sen(x) + 2 − 3 cos(x) − 6x + cos(x) + 2x − 3sen(x) − 3x2
= −4sen(x) − 2 cos(x) − 3x2 − 4x + 2

Teorema: Si
P (D) = an Dn + an−1 Dn−1 + · · · + a2 D2 + a1 D + a0
es un operador polinomial diferencial con coeficientes constantes y u(x) es una funci´on n veces
diferenciable, entonces:
P (D) [ueax ] = eax P (D + a)u ,

donde a es una constante.

Ejemplo Evaluar la siguiente expresi´on



D2 − D + 3 x3 e−2x

Al comparar la forma de esta expresión con la de teorema anterior, vemos que a = −2 y u(x) = x3 .
Por lo tanto, si hacemos

P (D − 2) = (D − 2)2 − (D − 2) + 3 = D2 − 5D + 9

se obtiene que





D2 − D + 3 x3 e−2x = e−2x (D − 2) x3 = e−2x D2 − 5D + 9 x3 = e−2x 9x3 − 15x2 + 6x .
2. Soluci´on de ecuaciones diferenciales lineales con el operador
polinomial
El m´etodo puede aplicarse a ecuaciones diferenciales de orden n, ya que si se tiene que

(D − m1 )(D − m2 ) · · · (D − mn )y = Q(x)

entonces se puede definir la función u como

u = (D − m2 ) · · · (D − mn )y

y lo que queda es resolver una ecuaci´on diferencial lineal de primer orden

(D − m1 )u = Q(x)

Una vez resuelta la ecuación para u se sustituye en:

(D − m2 )(D − m3 ) · · · (D − mn )y = u(x)

Repitiendo el proceso,

v = (D − m3 ) · · · (D − mn )y ⇒ (D − m2 )v = u(x)

Integrando:
(D − m3 ) · · · (D − mn )y = v(x)
Esta repetición nos llevará entonces hasta la soluci´on y(x).
Para resolver ecuaciones del tipo

an y (n) + an−1 y (n−1) + · · · + a2 y ′′ + a1 y ′ + a0 y = Q(x) , an 6= 0 ,

es decir, con coeficientes constantes, vamos a estudiar el siguiente ejemplo:
Ejemplo Resolver la ecuaci´on

y ′′′ + 2y ′′ − y ′ − 2y = e2x .

Al escribir esta ecuaci´on utilizando el operador diferencial resulta

D3 + 2D2 − D − 2 y = e2x ⇒ (D − 1) (D + 1) (D + 2) y = e2x

Si llamamos u = (D + 1) (D + 2) y entonces se puede ver que

(D − 1) u = e2x ⇒ u′ − u = e2x ⇒ u = e2x + C1 ex

Conocida la funci´on u entonces se tiene que

(D + 1) (D + 2) y = e2x + C1 ex

Repetimos el proceso, pero ahora hacemos v = (D + 2) y, por lo tanto

1 C1
(D + 1) v = e2x + C1 ex ⇒ v ′ + v = e2x + C1 ex ⇒ v = e2x + ex + c2 e−x
3 2
Conocida v, entonces
1 C1
(D + 2) y = e2x + ex + c2 e−x
3 2

Integrando una vez m´as:

1 C1 1
y ′ + 2y = e2x + ex + c2 e−x ⇒ y(x) = e2x + C1 ex + c2 e−x + C3 e−2x .
3 2 12
 Notemos que existe una manera bastante sencilla de resolver la ecuaci´on combinando varios
métodos. Si buscamos yh (x) a través de la ecuación caracterı́stica se obtiene que la solución es:
yh (x) = C1 ex + c2 e−x + C3 e−2x . Ahora bien, una soluci´on particular se obtiene al considerar
u(x) = e2x entonces v = 31 e2x y por lo tanto yp (x) = 121 e2x .
Ejercicio Resuelva la ecuación
y ′′ + y = ex .
2.1. El operador inverso
Los operadores suelen tener su inversa, esto significa que si

P (D)y = Q(x) ⇒ P −1 (D)P (D)y = P −1 (D)Q(x) ⇒ yp = P −1 (D)Q(x) ,

donde yp (x) es una solución particular de P (D)y = Q(x). Notemos que esto significa que
Z Z Z
−n
D Q(x) = ···
|{z} Q(x)dx
n veces

Ejemplo Evaluar
D−2 (2x + 3) .
Se tiene entonces lo siguiente:
Z Z
1 3
D−2 (2x + 3) = D−1 (2x + 3)dx = D−1 (x2 + 3x) = (x2 + 3x)dx = x3 + x2 .
3 2

De todos modos es fácil verificar que yp = 31 x3 + 23 x2 es una solución particular de D2 y = 2x + 3.

El operador inversa de P (D), es decir P −1 (D), puede tambi´en representarse como 1/P (D).
Debe quedar claro que su significado radica en el hecho de que al actuar sobre Q(x) produce una
soluci´on particular yp . Esta última notación resulta ser a veces más conveniente y dependiendo de
la forma de Q(x) puede resultar algunas veces de fácil solución.
En general, se tiene que sı́


P (D)y = an Dn + an−1 Dn−1 + · · · + a1 D + a0 y = Q(x) ,
 entonces
1 1
yp = [Q(x)] =   [Q(x)]
P (D) a1 a2 2 an−1 n−1 an n
a0 1 + a0 D + a0 D + · · · + a0 D + a0 D
1

= 1 + c1 D + c2 D2 + · · · + cn−1 Dn−1 + cn Dn [Q(x)] , a0 6= 0 ,
a0


donde el término 1 + c1 D + c2 D2 + · · · + cn−1 Dn−1 + cn Dn /a0 es la expansión en serie de 1/P (D).
Por otro lado, notemos que si a0 = 0, entonces



P (D)y = an Dn + an−1 Dn−1 + · · · + a1 D y = D an Dn−1 + an−1 Dn−2 + · · · + a1 y = Q(x) ,

Si ambos a0 = 0 y a1 = 0, entonces se puede seguir factorizando. Por lo tanto, en general:

P (D)y = Dr an Dn−r + · · · + ar+1 D + ar y = Q(x) .

esto significa que

1 1
yp = [Q(x)] = [Q(x)] ar 6= 0
P (D) Dr (an Dn−r + · · · + ar+1 D + ar )
 
1 1
= [Q(x)]
Dr (an Dn−r + · · · + ar+1 D + ar )

Muchos cálculos se simplificarán notablemente por la forma que pueda tener la funci´on Q(x),
veamos:
Si Q(x) = bxk , entonces
P (D)y = bxk
Note que si k = 0, entonces:
1 b
yp = [b] = , a0 6= 0 .
P (D) a0

porque b es constante.
 Ejemplos

1. Encontrar una soluci´on particular de

y ′′ − 2y ′ − 3y = 5 ⇒ (D2 − 2D − 3)y = 5 ,

aquı́, a0 = −3, a1 = −2, a2 = 1, b = 5 y k = 0, por lo tanto:

1 5
yp = [5] = − .
P (D) 3

2. Encontrar una soluci´on particular de

4y ′′ − 3y ′ + 9y = 5x2 ⇒ (4D2 − 3D + 9)y = 5x2 ,
ahora: a0 = 9, a1 = −3, a2 = 4, b = 5 y k = 2, esto significa que:
   
1  2 1  2 5 1 1 2  2 5 2 2 2
yp = 5x =
5x = 1+ D− D x = x + x− .
P (D) 9 1 − 93 D + 94 D2 9 3 3 9 3 3
Cuando
  se hace la expansión en serie de 1/P (D) no hace falta ir más alla de D3 , ya que
D3 x2 = 0.
3. Encontrar una solución particular de
y (5) − y (3) = 2x2 ⇒ (D5 − D3 )y = 2x2 ⇒ D3 (D2 − 1)y = 2x2
por lo tanto:  
1  2 1 1  2
yp = 2x = 3 2x
P (D) D D2 − 1
Trabajamos por partes,
1  2 1  2
2x = 2x
D2 − 1 (−1) (1 − D2 )
aquı́ a0 = −1, a1 = 0, a2 = 1, b = 2 y k = 2, por lo tanto
1  2 1  2 2

 2
2
2x = −2 2
x = −2 1 + D x = −2(x2 + 2)
D −1 (1 − D )
ahora tenemos que
  5 
1  2
 
−3 2
 x x3
yp = 3 −2(x + 2) = −2D x + 2 = −2 + .
D 60 3
Si Q(x) = beax , entonces
P (D)y = beax
Se puede demostrar (queda como ejercicio) que una soluci´on particular de esta ecuaci´on es
1 ax beax
yp = [be ] = , P (a) 6= 0 .
P (D) P (a)

Ejemplo Encuentre una solución particular de

y ′′′ − y ′′ + y ′ + y = 3e−2x
Entonces
(D3 − D2 + D + 1)y = 3e−2x ,
aquı́ b = 3 y a = −2, por lo tanto:
1  −2x  3e−2x 3e−2x 3 −2x
yp = 3e = = = − e ,
P (D) P (−2) (−2)3 − (−2)2 + (−2) + 1 13
Si Q(x) = b cos(ax) o Q(x) = b sen(ax). En este caso podemos utilizar la f´ormula de Euler y
volver al caso anterior

Ejemplo Encuentre una soluci´on particular de

y ′′ − 3y ′ + 2y = 3sen(2x)

entonces:
(D2 − 3D + 2)y = 3sen(2x) ,
sabemos que: e2ix = cos(2x) + isen(2x) y que la parte imaginaria de una solución particular
de
(D2 − 3D + 2)y = 3e2ix
será una solución de nuestro problema. Por lo tanto, con b = 3 y a = 2i se tiene:
 1  2ix  3e−2x 3e−2x 3e−2x 3(1 − 3i)e−2x
yp = 3e = = = =
P (D) P (2i) (2i)2 − 3(2i) + 2 −2(1 + 3i) −20
3 3
= − (1 − 3i)(cos(2x) + isen(2x)) = − [cos(2x) + 3sen(2x) + i(sen(2x) − 3 cos(2x))]
20 20
la parte imaginaria de esta ultima
´ ecuación ser´a nuestra soluci´on particular
3
yp = − [sen(2x) − 3 cos(2x)] .
20

Si Q(x) = Q1 (x) + Q2 (x) + · · · + Qn (x), entonces

1 1 1 1
yp = [Q(x)] = [Q1 (x)] + [Q2 (x)] + · · · + [Q2 (x)]
P (D) P (D) P (D) P (D)

Ejercicio Encuentre la solución general de

y ′′ − y = x3 + 3x − 4 .

3. Un método para resolver ecuaciones diferenciales lineales con

el operador polinomial inverso y utilizando fracciones parciales
La mejor manera de ilustrar el método es con un ejemplo. Consideremos la siguiente ecuaci´on
diferencial:

y ′′′ − 5y ′′ + 8y ′ − 4y = 2e4x ⇒ D3 − 5D2 + 8D − 4 y = 2e4x .
Busquemos su solución general.
Podemos notar que
D3 − 5D2 + 8D − 4 = (D − 2)2 (D − 1) ,
por lo tanto
1  4x  1  4x  1  4x 
yp = P −1 (D) [Q(x)] = 2e = 3 2e = 2e
P (D) D − 5D2 + 8D − 4 (D − 2)2 (D − 1)
 on en fracciones parciales
podemos hacer una expansi´
1 1 1 1
2
= − +
(D − 2) (D − 1) D−1 D−2 (D − 2)2

por lo tanto:
1  4x  1  4x  1  4x 
yp = 2e − 2e + 2e .
D−1 D−2 (D − 2)2
Ahora estamos en capacidad de resolver cada t´
ermino por separado:

2e4x 2e4x 2e4x e4x

yp = − + = .
3 2 4 6
Par otro lado, la familia de soluciones de la ecuación homogénea es

yh = (C1 + C2 x)e2x

on general resulta ser entonces:

y la soluci´

e4x
y(x) = (C1 + C2 x)e2x + .
6

Ejercicios Encuentre las soluciones particulares de las siguientes ecuaciones

1. y ′′ + y ′ = x2 + e2x

2. y ′′ − 2y ′ − 3y = 3eix

3. y ′′ − 2y ′ − 3y = 3 sen(x)

4. y ′′ − 2y ′ − 3y = 3 cos(x)

5. y ′′′ − y ′′ + y ′ − y = 4e−3x + 2x4

6. y (5) + 2y ′′′ + y ′ = 2x + sen(x) + cos(x)

 Variación de parámetros
Variación de parámetros para ED de orden 2
El método que presentamos en esta sección, llamado de variación de parámetros, es un procedimiento útil
y más general que el de coeficientes indeterminados para la obtención de una solución particular yp .x/ de
la ED lineal no homogénea y se basa en el conocimiento de la solución general de la ED lineal homogénea
asociada a la ED original.
Presentamos primeramente este método para las ED lineales de segundo orden. Diremos que el método de
variación de parámetros se usa para obtener una solución particular yp .x/ de la ED lineal

y 00 C p.x/y 0 C q.x/y D g.x/; (4.21)

a partir del conocimiento de la solución general de la ED lineal homogénea asociada

y 00 C p.x/y 0 C q.x/y D 0: (4.22)

Si suponemos que la solución general de la ED lineal homogénea (4.22) está dada por la combinación lineal

.x/ D c1 1 .x/ C c2 2 .x/;

esto significa que 1 .x/ & 2 .x/ son soluciones de la ED (4.22) y además su wronskiano W Œ1 .x/; 2 .x/ ¤ 0
en todo el intervalo .˛; ˇ/ donde las funciones p.x/ & q.x/ son continuas. Es decir, 1 .x/ & 2 .x/ forman
un conjunto fundamental de soluciones para la ED (4.22).
Supongamos pues que .x/ D c1 1 .x/ C c22 .x/ sea la solución general de y 00 C p.x/y 0 C q.x/y D 0.
El método de variación de parámetros propone que la solución particular yp .x/ tenga la misma forma que
.x/, pero permitiendo cambiar c1 y c2 . Esto es, propone que yp .x/ sea

yp .x/ D u1 1 .x/ C u2 2 .x/;

donde u1 & u2 son funciones de x desconocidas.

Para determinarlas, buscaremos qué condiciones deben cumplir para que la función propuesta yp sea solu-
ción de la ED (4.21). Lo primero que necesitamos hacer es calcular las derivadas yp0 , yp00 para utilizarlas en
(4.21) y ver cuáles propiedades deben cumplir u1 , u2 :

yp D u1 1 C u2 2 ) yp0 D u10 1 C u1 10 C u20 2 C u2 20 : (4.23)

Antes de obtener yp00, suponemos que

u10 1 C u20 2 D 0:
Esto se hace con la finalidad de que en la expresión de yp00 no aparezcan u100 & u200, ya que la inclusión de
estas segundas derivadas en yp00 haría mucho más compleja la obtención de las funciones u1 & u2 .
Luego entonces, usando la condición anterior en (4.23):

yp0 D u1 10 C u2 20 ;

por lo cual, al volver a derivar, obtenemos:

yp00 D u10 10 C u1 100 C u20 20 C u2 200 :

Ahora bien, yp es solución de la ED lineal y 00 C p.x/y 0 C q.x/y D g.x/ sólo si se cumple que

yp00 C p.x/yp0 C q.x/yp D g.x/:

÷
Esta condición, al usar las expresiones obtenidas de yp , yp0 & yp00 nos da:

Œu10 10 C u1 100 C u20 20 C u2 200  C p.x/Œu1 10 C u2 20  C q.x/Œu1 1 C u2 2  D g.x/:
š š
yp00 yp0 yp

Al factorizar u1 , u2 :

u10 10 C u1 Œ100 C p.x/10 C q.x/1  C u20 20 C u2 Œ200 C p.x/20 C q.x/2  D g.x/:

Pero por ser 1 & 2 soluciones de la ED homogénea se cumple:

100 C p.x/10 C q.x/1 D 0 & 200 C p.x/20 C q.x/2 D 0:

Entonces debe cumplirse que

u10 10 C u20 20 D g.x/:
Hemos obtenido un sistema de dos ecuaciones con dos incógnitas:
(
u10 1 C u20 2 D 0I
(4.24)
u10 10 C u20 20 D g.x/ :

El determinante s del sistema es  

1 2 
s D  0 0  D W.1 ; 2/:
1 2
Y debido a que el wronskiano W.1 ; 2 /.x/ ¤ 0, entonces s ¤ 0, por lo que el sistema de ecuaciones tiene
una única solución. Dicha solución única la calculamos usando la regla de Cramer:
 
 0 2 
 
g.x/  0  g.x/2 g.x/2 .x/
2
u10 D D ) u10 .x/ D :
s W.1 ; 2 / W.1 .x/; 2 .x//
 
 1 0 
 0 
 g.x/ g.x/1 g.x/1 .x/
u20 D 1 D ) u20 .x/ D :
s W.1 ; 2 / W.1 .x/; 2 .x//

De donde obtenemos u1 & u2 mediante integración

g.x/2 .x/ g.x/1 .x/

Z Z
u1 D dx & u2 D dx:
W Œ1 .x/; 2 .x/ W Œ1 .x/; 2 .x/

Sustituyendo u1 .x/ & u2 .x/ en yp .x/ se tiene que la solución particular de la lineal es

g.x/2 .x/ g.x/1 .x/

Z Z
yp .x/ D 1 .x/ dx C 2 .x/ dx: (4.25)
W Œ1.x/; 2 .x/ W Œ1 .x/; 2 .x/

Finalmente, podemos escribir la solución general de la ED lineal como

y.x/ D yp .x/ C .x/I

y.x/ D yp .x/ C Œc11 .x/ C c2 2 .x/:

con la yp .x/ obtenida previamente.

Ejemplo 4.7.1 Utilizando el método de variación de parámetros, calcular una solución particular y escribir la
solu-ción general de la ED lineal
1
x 2y 00 4xy 0 C 6y D ;
x
considerando que y1 D x 2 & y2 D x 3 forman un conjunto fundamental de soluciones para la ED homogénea asociada

x 2 y 00 4xy 0 C 6y D 0:

H Sea yp .x/ D u1 .x/y1 .x/ C u2 .x/y2 .x/ una solución particular.

Entonces:
yp D u1 x 2 C u2 x 3 ) yp0 D u10 x 2 C 2u1 x C u20 x 3 C 3u2 x 2 :
Suponiendo que

u10 x 2 C u20 x 3 D 0; (4.26)

se tiene que

yp0 D 2u1 x C 3u2 x 2 & yp00 D 2u10 x C 2u1 C 3u20 x 2 C 6u2 x:

Sustituyendo en la ED normalizada
4 0 6 1
y 00 y C 2y D 3;
x x x
se obtiene:
4 6 1
.2u10 x C 2u1 C 3u20 x 2 C 6u2 x/ .2u1 x C 3u2 x 2 / C 2 .u1 x 2 C u2 x 3 / D 3 )
x x x
1
) 2xu10 C u1 .2 8 C 6/ C 3x 2 u20 C u2 .6x 12x C 6x/ D 3 )
x
1
) 2xu10 C 3x 2u20 D 3 : (4.27)
x
Entonces u10 & u20 satisfacen el sistema formado por las ecuaciones (4.26) y (4.27)
(
x 2 u10 C x 3 u20 D 0I
2xu10 C 3x 2u20 D x 3 :
 El determinante del sistema es  2 3

x x 4
W D  2  D 3x
 2x 4 D x 4 :
2x 3x

La solución del sistema es

 0 x3 
 
 3 2

x 3x 1
u10 D D 4 D x 4:
 2 W  x
x 0 
2x x 3  x 1

0
u2 D D 4 D x 5:
W x
Integrando:
3
1 x
Z
4
u1 D x C C1 D x 3 C C1 :
dx D
3 3
x 4 1
Z
5
u2 D x dx D C C2 D x 4 C C2 :
4 4
1 1 4
Tomando C1 D 0 y C2 D 0, obtenemos, u1 D x 3
y u2 D x . Por ello, una solución particular es
3 4
1 1 4 3 1 1 1 1
yp D u1 x 2 C u2 x 3 D x 3 x 2 x x D x 1
x 1
D x :
3 4 3 4 12
Entonces la solución general de la ED lineal es 

1 1
y D yp .x/ C c1 y1 .x/ C c2 y2 .x/ ) y D x C c1 x 2 C c2 x 3 :
12

Ejemplo 4.7.2 Utilizando el método de variación de parámetros, calcular una solución particular y escribir la solu-
ción general de la ED lineal
 x 2 y 00 xy 0 C y D 4x ln x;
considerando que y1 D x & y2 D x ln x forman un conjunto fundamental de soluciones para la ED homogénea
asociada
x 2y 00 xy 0 C y D 0:

H Sea yp .x/ D u1 .x/y1 .x/ C u2 .x/y2 .x/ una solución particular propuesta por este método:
Entonces:
yp D u1 x C u2 x ln x & yp0 D u10 x C u1 C u20 x ln x C u2 .1 C ln x/:
Imponiendo la condición

u10 x C u20 x ln x D 0; (4.28)

resulta  
1
yp0 D u1 C u2 .1 C ln x/ & yp00 D u10 C u20 .1 C ln x/ C u2 :
x
Sustituyendo en la ED normalizada

1 0 1 4
y 00 y C 2 y D ln x;
x x x
se obtiene:
 
1 1 1 1 4
u10 C u20 .1 C ln x/ C u2 Œu1 C u2 .1 C ln x/ C u1 C u2 ln x D ln x )
x x x x x

(aquí factorizamos u1 & u2 )

   
1 1 1 1 1 1 4
) u10 C u20 .1 C ln x/ C u2 ln x C ln x C u1 C D ln x )
x x x x x x x
 4
) u10 C .1 C ln x/u20 D ln x: (4.29)
x
Entonces u10 & u20 satisfacen el sistema conformado por las ecuaciones (4.28) y (4.29).
(
xu10 C .x ln x/u20 D 0I
u10 C .1 C ln x/u20 D 4x 1 ln x:

El determinante del sistema es

x x ln x 
 
W D   D x C x ln x x ln x D x:
1 1 C ln x 

La solución del sistema es

ln
 
 0 x x 
4x ln x 1 C ln x
 1 
4.ln x/2
0
u1 D D D 4x 1 .ln x/2 :
 W  x
x 0 
 1 4x 1 ln x 4 ln x

0
u2 D D D 4x 1 ln x:
W x
Integrando:

dx 4
Z Z
u1 D 4 x 1 .ln x/2 dx D 4 .ln x/2 D .ln x/3 C C1 :
x 3
dx
Z Z
u2 D 4 x 1 ln x dx D 4 .ln x/ D 2.ln x/2 C C2 :
x
4
Tomando simplemente u1 D .ln x/3 & u2 D 2.ln x/2 , se tiene que una solución particular es
3
4 2
yp D u1 x C u2 x ln x D .ln x/3 x C 2.ln x/2 x ln x ) yp .x/ D x.ln x/3 :
3 3
Y ahora, la solución general de la ED está dada por

2
y D yp .x/ C c1 y1 .x/ C c2 y2 .x/ ) y D x.ln x/3 C c1x C c2 x ln x:
3
Ejemplo 4.7.3 Utilizando el método de variación de parámetros, calcular una solución particular y escribir la solu-

ción general de la ED lineal
y 00 C y D sec2 x:

H Primero se obtiene un conjunto fundamental de soluciones para la ED homogénea asociada

y 00 C y D 0:

Proponiendo y D e rx para resolver la ED homogénea:

p
r2 C 1 D 0 ) r D ˙ 1 D 0 ˙ 1i:

Entonces: (
y1 D e 0x sen 1x D sen x:
y2 D e 0x cos 1x D cos x:
Estas dos funciones forman un conjunto fundamental de soluciones de la ED homogénea asociada.
Se propone como solución particular:
yp D u1 sen x C u2 cos x;
se obtiene:
yp0 D u10 sen x C u1 cos x C u20 cos x u2 sen x:
Suponiendo:
u10 sen x C u20 cos x D 0; (4.30)
se tiene:
yp0 D u1 cos x u2 sen x & yp00 D u10 cos x u1 sen x u20 sen x u2 cos x:
Sustituyendo en
yp00 C yp D sec2 x;
se obtiene:
 .u10 cos x u1 sen x u20 sen x u2 cos x/ C .u1 sen x C u2 cos x/ D sec2 x )
) u10 cos x u20 sen x D sec2 x: (4.31)
De este modo u10 & u20 satisfacen el sistema formado por (4.30) y (4.31):
(
u10 sen x C u20 cos x D 0:
u10 cos x u20 sen x D sec2 x:

El determinante del sistema es

sen x cos x 
 
W D  D sen 2 x cos 2 x D 1 ) W D 1:
cos x sen x 

La solución del sistema es

cos x 
 
 0
sec x sen x 
 2
sec2 x cos x
u10 D D D sec x:
W 1
sen x
 
0 
cos x sec2 x  sen x sec2 x
 
u20 D D D sen x sec2 x:
W 1
Integrando:
Z
u1 D sec x dx D ln.sec x C tan x/ C C I
Z Z
2
u2 D sen x sec x dx D sec x tan x dx D sec x C C:

Tomando u1 D ln.sec x C tan x/ & u2 D sec x, se tiene que una solución particular es
yp D u1 sen x C u2 cos x ) yp D .sen x/ ln.sec x C tan x/ .sec x/ cos x )
) yp D .sen x/ ln.sec x C tan x/ 1:
1 sen x C c2 cos x:
Entonces la solución general es

y D .sen x/ ln.sec x C tan x/ 1 C c
Ejemplo 4.7.4 Utilizando el método de variación de parámetros, calcular una solución particular y escribir la solu-
ción general de la ED lineal
e 3x
y 00 3y 0 C 2y D :
1 C ex
H Primero se obtiene un conjunto fundamental de soluciones para la ED homogénea asociada
y 00 3y 0 C 2y D 0I
proponemos y D e rx y se obtiene:
r2 3r C 2 D 0 ) r1 D 1 & r2 D 2;
entonces:
y1 D e x & y2 D e 2x ;
son funciones que forman un conjunto fundamental de soluciones.
Se propone como solución particular:
yp D u1 e x C u2 e 2x ;
y se obtiene:
yp0 D u10 e x C u1 e x C u20 e 2x C 2u2 e 2x :
Suponiendo
u10 e x C u20 e 2x D 0; (4.32)
se tiene:
yp0 D u1 e x C 2u2 e 2x & yp00 D u10 e x C u1 e x C 2u20 e 2x C 4u2 e 2x :
Sustituyendo en
e 3x
yp00 3yp0 C 2yp D ;
1 C ex
se obtiene:
e 3x
.u10 e x x
C u1 e C 2u20 e 2x 2x x 2x
C 4u2 e / 3.u1 e C 2u2 e / C 2.u1 e C u2 e / D ) x 2x
1 C ex
0 x x 0 2x 2x e 3x
) u1 e C u1 e .1 3 C 2/ C 2u2 e C u2 e .4 6 C 2/ D )
1 C ex
0 x 0 2x e 3x
) u1 e C 2u2 e D : (4.33)
1 C ex
Entonces u10 & u20 satisfacen el sistema formado por (4.32) y (4.33):

 x u10 C e 2x u20 D 0:
e
e 3x
e x u10 C 2e 2x u20 D :
1 C ex
El determinante del sistema es
  x
e 2x 

e
W D  x D 2e 3x e 3x D e 3x ) W D e 3x :
e 2e 2x 
La solución del sistema es
 
 0
 e 2x 
 e 3x 
2x 
 2e 
 1 C ex e 2x e 3x e 2x
u10 D D 3x D :
W e .1 C e x / 1 C ex
 
e x 0 


e x e 3x 
 1 C ex  e x e 3x ex
u20 D D 3x D :
W e .1 C e x / 1 C ex

e 2x ex
Z
Integrando:
Z1 C e x
dx & u2 D dx:
1 C ex
u1 D
Utilizando el cambio de variable t D 1 C e x :

u1 D ln.1 C e x / .1 C e x / C C1 :
u2 D ln.1 C e x / C C2 :

Tomando u1 D ln.1 C e x / .1 C e x / & u2 D ln.1 C e x /, se obtiene la solución particular

yp D u1 e x C u2 e 2x D Œln.1 C e x / .1 C e x /e x C Œln.1 C e x /e 2x D

D e x ln.1 C e x / C e 2x ln.1 C e x / e x .1 C e x / D Œe x ln.1 C e x /Œ1 C e x  e x .1 C e x / )
) yp .x/ D e x .1 C e x /Œln.1 C e x / 1:

Por lo tanto, la solución general de la ED es

y D yp .x/ C c1y1 .x/ C c2 y2 .x/ ) y D e x .1 C e x /Œln.1 C e x / 1 C c1e x C c2 e 2x :

 

De acuerdo a lo discutido en esta sección, una solución particular yp de la ED lineal normalizada

y 00 C p.x/y 0 C q.x/y D g.x/

tiene la forma:
yp .x/ D u1 1 .x/ C u2 2 .x/;
donde f 1 .x/; 2 .x/ g es un conjunto fundamental de soluciones de la ED homogénea asociada:

y 00 C p.x/y 0 C q.x/y D 0I

y las funciones u1 & u2 se obtienen integrando, respectivamente, las funciones u10 & u20 que son soluciones
del sistema ˚
1 u10 C 2 u20 D 0I
10 u10 C 20 u20 D g.x/I

el cual resolvemos por el método de Cramer. Aún más:

W1 W2
u10 D & u20 D ;
W W
donde W D W Œ1 .x/; 2 .x/ 6D 0 [para toda x en el intervalo donde p.x/ & q.x/ sean continuas] y W1
& W2 están dadas mediante un determinante en el que se sustituye la columna 1 & 2, respectivamente, de
W por la columna:
 
0
:
g.x/
Es conveniente señalar las fortalezas y debilidades de este método en comparación con el método de coefi-
cientes indeterminados, tratado en la sección anterior.

1. El método de coeficientes indeterminados es en muchos casos más sencillo y fácil de aplicar que el
de variación de parámetros, pero tiene esta limitación: sólo es aplicable cuando g.x/ tiene la forma
de un polinomio, exponencial, combinación lineal de senos y cosenos o una suma de las funciones
mencionadas. Además la ED lineal debe tener coeficientes constantes.
 2. El método de variación de parámetros en cambio es aplicable en pricipio para cualquier función g.x/
y, en este sentido, es mucho más general. Además se puede usar para ED lineales aun cuando no
sean de coeficientes constantes. Lo único que se necesita es contar con un conjunto fundamental
de soluciones f 1 ; 2 g. Esta generalidad del método tiene un costo, puesto que para obtener una
solución particular yp es necesario calcular las integrales:

2 .x/g.x/ 1 .x/g.x/
Z Z
u1 D dx & u2 D dx:
W.x/ W.x/

(La resolución de estas integrales es dificultosa.)

Ejercicios 4.7.1 Variación de parámetros para ED de orden 2. Soluciones en la página 467

Utilizando variación de parámetros, calcular una solución particular y escribir la solución general de la ecuación
diferencial dada. Considerar que las funciones y1 D y1 .x/ & y2 D y2 .x/ forman un conjunto fundamental de
soluciones para la ecuación homogénea asociada.

1. x 2 y 00 6xy 0 C 10y D 8x 3 I y1 D x 2 & y2 D x 5 .

p 1
2. x 2 y 00 xy 0 3y D 30 xI y1 D x 3 & y2 D .
x
65
3. x 2 y 00 C xy 0 C 8y D p
3
I y1 D x 4 & y2 D x 2 .
x
6
4. x 2 y 00 C 8xy 0 C 12y D 2
I y1 D x 3 & y2 D x 4 .
x
5. x 2 y 00 6xy 0 C 10y D 4x ln x 5xI y1 D x 5 & y2 D x 2 .

Utilizando variación de parámetros, determinar una solución particular y escribir la solución general de la ecuación
diferencial dada.
6. y 00 y D ex . 14. y 00 C 9y D 9 sec 3x tan 3x .

7. y 00 yDe x
. 15. y 00 yDe 2x
sen e x
.

8. y 00 C y D sen x . 16. y 00 C 4y D sen 2 2x .

9. y 00 C y D cos x . 17. y 00 C 4y D cos 2 2x .

ex
10. y 00 2y 0 C y D 6xe x . 18. y 00 2y 0 C y D .
x
11. y 00 C 2y 0 C y D 12xe x
. e x
19. y 00 C 2y 0 C y D .
12. y 00 C y D tan x . x
1
13. y 00 C 4y D 4 sec 2x . 20. y 00 C 3y 0 C 2y D .
1 C e 2x
 Variación de parámetros para ED de orden n
Descripción del método general
Una vez discutido el método de variación de parámetros para ecuaciones diferenciales de orden 2, en esta
sección extenderemos dicho método a ecuaciones diferenciales de orden n para n > 2. Así, consideraremos
el caso de la ED lineal no homogénea normalizada:

y .n/ C an 1 .x/y .n 1/
C


C a0 .x/y D g.x/I (4.34)

cuya ED homogénea asociada es

y .n/ C an 1 .x/y .n 1/
C


C a0 .x/y D 0: (4.35)

Suponemos conocido el conjunto fundamental de soluciones:

f1 .x/; 2 .x/;


; n.x/g;

con el cual podemos formar la solución general de la ecuación (4.35),

.x/ D c1 1 .x/ C c2 2 .x/ C


C cnn.x/: (4.36)

En lo que sigue, supondremos también que las funciones an 1 .x/;


; a0 .x/ & g.x/ son continuas en el
intervalo .˛; ˇ/, donde el conjunto de funciones f1.x/; 2 .x/;


; n.x/g satisface:
 
 1 .x/ 2 .x/


n.x/ 
 0 0 0 
 1 .x/  2 .x/


 n .x/ 
W Œ1 .x/; 2 .x/;


; n.x/ D  :: :: :: ::  ¤ 0, para todo x 2 .˛; ˇ/:
 
 : : : : 
 .n 1/

1 .x/ 2.n 1/ .x/


n.n 1/ .x/

Como vimos en la sección anterior, el método se apoya en la idea de que las constantes c1 ; c2;


; cn de la
ecuación (4.36) las podemos reemplazar por un conjunto de funciones indeterminadas fu1 .x/; u2 .x/;


; un.x/g
por conocer, que permiten generar una solución particular

yp .x/ D u1 .x/1 .x/ C u2 .x/2 .x/ C


C un .x/n .x/ (4.37)

de la ecuación (4.34). Como queremos determinar n funciones, es de esperarse que debemos imponer n
condiciones a las funciones uj .x/; j D 1; 2; : : :; n.
Es claro que una de las condiciones debe ser sin lugar a dudas el cumplimiento de la ecuación (4.34).
La idea central de este método radica en no resolver para las funciones uj .x/ ecuaciones diferenciales que
sean de orden mayor que 1, es decir, buscaremos las restantes n 1 condiciones de manera que nunca
.k/
tengamos que considerar alguna relación en la que intervenga alguna derivada uj .x/ para k > 1.
Con esto en mente, a partir de (4.37) obtenemos al derivar:

yp0 D Œu1 10 C u2 20 C


C un n0  C Œu10 1 C u20 2 C


C un0 n :

Para asegurar que no aparezca alguna uj00.x/, al calcular la segunda derivada, requerimos que en el intervalo
.˛; ˇ/:

u10 1 C u20 2 C


C un0 n D 0: (4.38)

Entonces, tenemos ahora:

yp0 D u1 10 C u2 20 C


C unn0 :
Por lo tanto,
yp00 D Œu1 100 C u2 200 C


C un n00 C Œu10 10 C u20 20 C


C un0 n0 :
Por la misma razón que expusimos anteriormente, ahora planteamos la condición:

u10 10 C u20 20 C


C un0 n0 D 0: (4.39)

De donde se desprende que

yp00 D u1 100 C u2 200 C


C unn00: (4.40)

Si proseguimos de la misma manera, determinamos que se deben cumplir las relaciones:

.h/ .h/
u10 1 C u20 2 C


C un0 n.h/ D 0; h D 0; 1; 2; : : :; n 2; (4.41)

yp.k/ D u1 1.k/ C u2 2.k/ C


C un n.k/ para k D 0; 1; 2; : : :; n 1: (4.42)

Finalmente, si para k D 0; 1; 2; : : :; n 1 usamos las relaciones (4.42) en (4.34), hallamos:
 ™
u1 1.n/ C


C un n.n/ C u10 1.n
.n/
yp
1/
C


C un0 n.n 1/
Can 1 Œu1 1.n
C


Cu 
 
1/
C


C
.n 1/
yp
.n 1/
n n

C a0 Œu1 1 C


C unn  D g.x/: (4.43)


yp

Ahora reacomodamos la ecuación (4.43) factorizando u1 , u2


un :

u1 Œ1.n/ C an 1 1.n 1/
C


C a0 1  C


C un Œn.n/ C an 1 n.n 1/
C


C a0 n C
.n 1/
C Œu10 1 C


C un0 n.n 1/  D g.x/: (4.44)
Observamos que como cada una de las funciones del conjunto f1 .x/; 2 .x/;


; n.x/g satisface a la ecuación
(4.35), todos los términos en (4.44) son iguales a cero con excepción del último; esto permite llegar a la con-
clusión de que
.n 1/
u10 1 C


C un0 n.n 1/
D g.x/: (4.45)
Al reunir todas las condiciones indicadas en (4.41) junto con la ecuación (4.45), concluimos que las funciones
incógnitas u

0 0
1 ; : : : ; un satisfacen a las condiciones:



 u10 1 C u20 2 C


C un0 n D 0I
u 0  0 C u 0  0 C


C un0 n0 D 0I


 1 1 2 2


::
: (4.46)

 .n 2/ .n 2/ .n 2/



 u10 1 C u20 2 C


C un0 n D 0I
u 0  .n 1/ C u 0  .n 1/ C


C u 0  .n
 1/
D g.x/:
1 1 2 2 n n

Ahora bien, como indicamos anteriormente:

W Π1.x/; 2 .x/;


; n.x/ ¤ 0 , para todo x 2 .˛; ˇ/ :

Por lo tanto, considerando (4.46) como un sistema de n ecuaciones con n incógnitas u10 ; : : : ; un0 , podemos
utilizar la regla de Cramer para obtener la solución única para u10 ; : : : ; un0 . Obtenemos:
 Wk Wk
uk0 D D I k D 1; 2; : : : ; n: (4.47)
W Œ1.x/; 2 .x/;


; n .x/ W
Donde Wk difiere de W Œ1 .x/; 2 .x/;


; n .x/ D W en que tiene, como columna k-ésima, a la columna:
 
0
 0 
 ::  :
 
 : 
g.x/
Wk
Las funciones W D W.x/ y Wk D Wk .x/ son continuas; en consecuencia, las expresiones uk0 D son
W
integrables; por lo tanto, integrando se determinan las funciones incógnitas u1 .x/; : : : ; un .x/.
Finalmente observe que, si en lugar de la ecuación diferencial (4.34), tuviéramos
an .x/y .n/ C an 1 .x/y .n 1/
C


C a0 .x/y D b.x/;
b.x/
–es decir, si an .x/ no fuese idénticamente igual a 1–, entonces habría que poner g.x/ D en lugar de
an .x/
g.x/ en la última ecuación de (4.46). En otras palabras, habría que normalizar la ED.

Ejemplo 4.7.5 Resolver la ED y .3/ C 4y 0 D cot 2x.

H Primero observamos que esta ecuación diferencial no puede ser resuelta por medio del método de
coeficientes indeterminados debido a la presencia de la función cotangente. La ED está normalizada, ya
que el coeficiente de la mayor derivada de la ecuación es igual a 1.
Para la solución, determinamos en primer lugar la ecuación característica correspondiente a la ED asociada:

r 3 C 4r D r.r 2 C 4/ D 0:

Las raíces de esta ecuación algebraica son r1 D 0I r2;3 D ˙2i. Por lo tanto, el conjunto fundamental de
soluciones está integrado por las 1funciones
2 3

 D 1;  D cos 2x &  D sen 2x:

El wronskiano W D W .1 ; 2; 3 / es

 cos 2x sen 2x 
 
1
2 sen 2x 2 cos 2x  D 8 sen 2 2x C 8 cos 2 2x D

W D W.1 ; 2; 3 / D 0
0 4 cos 2x 4 sen 2x 
D 8Œsen 2 2x C cos 2 2x D 8:

De acuerdo a lo discutido en el procedimiento general, requerimos hallar W1 ; W2 y W3 . Tenemos:

cos sen
 
 0 2x 2x 
2 sen 2x 2 cos 2x  D cot 2xŒ2 cos 2 2x C 2 sen 2 2x D 2 cot 2xI
 
W1 D  0
cot 2x 4 cos 2x 4 sen 2x

sen
 
1 0 2x 
2 cos 2x  D cot 2xŒ2 cos 2x D 2 cot 2x cos 2xI
 
W2 D 0 0
0 cot 2x 4 sen 2x

y también:

cos 2x
 
1 0 
2 sen 2x 0  D 2 sen 2x cot 2x D 2 cos 2x:

W3 D 0
0 4 cos 2x cot 2x

De esta manera, obtenemos para u1 ; u2 ; u3 las siguientes expresiones:

Para u1 ,
2 cot 2x
 
W1 1 1 1
Z Z
u1 D dx D dx D ln.sen 2x/ D ln.sen 2x/:
W 8 4 2 8
Para u2 ,

W2 2 cot 2x cos 2x 1 cos 2 2x 1 sen 2 2x

1
Z Z Z Z
u2 D dx D dx D dx D dx D
W 8 4 sen 2x 4 sen 2x
Z   
1 1 1 1
Z
D csc 2x dx sen 2x dx D ln.csc 2x cot 2x/ C cos 2x D
4 4 2 2
1 1
D ln.csc 2x cot 2x/ cos 2x:
8 8
Finalmente para u3 ,

2 cos 2x
 
W3 1 1 1
Z Z
u3 D dx D dx D sen 2x D sen 2x:
W 8 4 2 8
Estamos ahora listos para obtener la solución general la cual, como sabemos, resultará sumando la solución
complementaria y la solución particular que calculamos con el método de variación de parámetros. De esta
manera la solución general de la ED es

y D Œc11 C c2 2 C c3 3  C Œu1 1 C u2 2 C u3 3  D
D c1 C c2 cos 2x C c3 sen 2xC
1 1 1 2 1 2
C ln.sen 2x/ ln.csc 2x cot 2x/ cos 2x cos 2x sen 2x D
8 8 ā8 8
1
D
8
1 1
D c1 C c2 cos 2x C c3 sen 2x C ln.sen 2x/ ln.csc 2x cot 2x/ cos 2x:
8 8

La ecuación diferencial de Cauchy-Euler

El ejemplo anterior muestra cómo generar una solución particular y la consecuente solución general de
la ED por el método de variación de parámetros si tan sólo se conoce un conjunto fundamental de solu-
ciones de la ecuación diferencial homogénea asociada. Hay que decir que tuvimos a nuestro favor que la
ecuación homogénea es de coeficientes constantes, característica que facilitó la determinación del conjunto
fundamental.
Si los coeficientes no son constantes, la tarea puede resultar más compleja con excepción de algunos casos
particulares, como el que discutimos a continuación.

 La ED de Cauchy-Euler de orden n tiene la siguiente forma general:

n n 1
nd y n 1d y dy
an x C an 1 x C


C a1 x C a0 y D g.x/; donde a0 ; a1 ; : : : ; an son constantes.
dx n dx n 1 dx

La ED de Cauchy-Euler de orden 3 tiene la siguiente forma:

d 3y 2
2d y dy
a3 x 3 C a2 x C a1 x C a0 y D g.x/: (4.48)
dx 3 dx 2 dx
Este tipo de ED se puede reducir a una ED de coeficientes constantes si se realiza el siguiente cambio de
variable:
x D et :
Derivamos usando la regla de la Cadena:
dy dy dy
dy dy
y0 D D dt D dt D dtt D e t :
dx dx d t e dt
.e /
dt dt
Es decir,
dy dy
y0 D De t : (4.49)
dx dt 

Al calcular la segunda derivada:

 
dy d 0 dy 2
e t td y t dy
2
d y dy 0
dt dt e e 
d 2y dy

y 00 D 2 D D dt D D dt 2 dt D e 2t
:
dx dx dx d t et dt 2 dt
.e /
dt dt

Es decir,

d 2y d 2y
 
2t dy
y 00 D 2 D e : (4.50)
dx dt 2 dt

Para la tercera derivada encontramos que

  2 
dy 00 d d y dy
3 e 2t
d y dy 00
dt dt 2 dt
y 000 D 3 D D dt D D
dx dx dx d t
.e /
dt dt
 3 2   2 
d y d y d y dy
e 2t 2e 2t  3
d 2y

3 2 dt 2 dt
dt dt 3t d y dy
D D e 3 2 C2 :
et dt 3 dt dt
Es decir,

d 3y d 3y d 2y
 
3t dy
000
y D 3 De 3 2 C2 : (4.51)
dx dt 3 dt dt

Al usar (4.49), (4.50) y (4.51) en (4.48), y tomando en consideración que x 1

D e t , obtenemos una ED con
coeficientes constantes y variable independiente t.

Ejemplo 4.7.6 Resolver la ED x 3 y .3/ C 3x 2y 00 3xy 0 D x ln x.

H Si usamos (4.49), (4.50) y (4.51) en la ED, obtenemos:

d 3y d 2y d 2y
       
3t 3t dy 2t 2t dy t t dy
e e 3 2 C2 C 3e e 3e e D e t ln e t :
dt 3 dt dt dt 2 dt dt

Que se simplifica en
d 3y d 2y dy d 2y dy dy
3 2 C2 C3 2 3 3 D te t ;
dt 3 dt dt dt dt dt
o bien en
d 3y dy
4 D te t :
dt 3 dt

Ésta es una ED con coeficientes constantes; para resolverla aplicamos el procedimiento conocido. En este
caso, la ecuación característica asociada a la ED homogénea es

r 3 4r D r.r 2 4/ D r.r 2/.r C 2/ D 0;

cuyas raíces son r1 D 0, r2 D 2 & r3 D 2. Deducimos entonces que un conjunto fundamental de soluciones
está integrado por las funciones:

1 D 1; 2 D e 2t & 3 D e 2t
:

El wronskiano correspondiente para este caso es

 1 e 2t e 2t 
 

W D W .1 ; 2; 3 / D 0 2e 2t 2e 2t  D 8 C 8 D 16:


0 4e 2t 4e 2t 


Y las funciones W1 , W2 y W3 son

e 2t e 2t 
 
0
2e 2t 2e 2t  D te t Π2 4te t :

W1 D  0 2 D
te t 4e 2t 4e 2t 

e 2t 
 
1 0
2e 2t  D 2te t :

W2 D 0 0
0 te t 4e 2t 

e 2t
 
1 0 
2e 2t 0  D 2te 3t :

W3 D 0
0 4e 2t te t 

Por último, integramos por partes para determinar las funciones incógnitas u1 ; u2 & u3 :

W1 4te t 1 1 t
Z Z Z
u1 D dt D dt D te t dt D e .t 1/:
W 16 4 4
W2 2te t 1 1 t
Z Z Z
u2 D dt D dt D te t dt D e .t C 1/:
W 16 8 8
W3 2te 3t 1 1
Z Z Z
u3 D dt D dt D te 3t dt D e 3t .3t 1/:
W 16 8 72
En conclusión, la solución de la ED está dada por

y D Œc1 1 C c22 C c33  C Œu1 1 C u2 2 C u3 3  D

1 t 1 t 1
D Œc1 C c2e 2t C c3 e 2t  e .t 1/ e .t C 1/e 2t C e 3t .3t 1/e 2t
D
4  8 72 
1 1 1 1 1 1
D Œc1 C c2e 2t C c3 e 2t  C e t tC t C t D
4 4 8 8 24 72
1
D Œc1 C c2e 2t C c3 e 2t  C e t .1 3t/:
9
Sin embargo, en la ecuación inicial, no es t la variable independiente, sino x.
De x D e t , hallamos que t D ln x. Si lo utilizamos en el resultado anterior, encontramos:
1
y D c1 C c2 x 2 C c3 x 2
C x.1 3 ln x/:
9


Ejemplo 4.7.7 Resolver el siguiente PVI:

ln x
x 2 y .3/ xy 00 C y 0 D ; con y.1/ D 0; y 0 .1/ D 1 & y 00 .1/ D 1:
x
H Primero multiplicamos la ecuación por x a fin de llevarla a una ecuación del tipo Cauchy-Euler:

x 3 y .3/ x 2 y 00 C xy 0 D ln x:

Si hacemos ahora el cambio de variable x D e t e incorporamos los resultados (4.49), (4.50) y (4.51):
 3
d 2y
    2   
3t 3t d y dy 2t 2t d y dy t t dy
e e 3 2 C2 e e Ce e D t:
dt 3 dt dt dt 2 dt dt
Al simplificar, hallamos:
d 3y d 2y dy d 2 y dy dy
3 2 C2 C C D t;
dt 3 dt dt dt 2 dt dt

o bien:
d 3y d 2y dy
4 2 C 4 D t:
dt 3 dt dt
Ésta es una ED con coeficientes constantes. La ecuación característica correspondiente a la ED homogénea
asociada es
r 3 4r 2 C 4r D r.r 2 4r C 4/ D r.r 2/2 D 0:
Por lo tanto las raíces de la ecuación son r1 D 0 y r2 D r3 D 2. En consecuencia, las funciones que integran
el conjunto fundamental de soluciones de la ED homogénea asociada son

1 D 1; 2 D e 2t & 3 D te 2t :

Con ellas podemos calcular el wronskiano W D W.1 ; 2 ; 3/:

1 e 2t 2t
 
te 
W D W.1 ; 2; 3 / D 0 2e  .2t C 1/e D .8t C 8/e 4t  .8t C 4/e 4t D 4e 4t :
2t  2t
 

0 4e 2t  .4t C 4/e 2t 

De manera similar:
e 2t te 2t
 
0 
2e 2t .2t C 1/e 2t  D tŒ.2t C 1/e 4t 2te 4t  D te 4t ;
 
W1 D 0
t 4e 2t .4t C 4/e 2t 

te 2t
 
1 0 
0 .2t C 1/e 2t  D t.2t C 1/e 2t D . 2t 2 t/e 2t ;
 
W2 D 0
0 t .4t C 4/e 2t 

e 2t 0
 
1
2e 2t 0 D 2te 2t :

W3 D 0
0 4e 2t t 

Así, por el método de variación de parámetros, las funciones incógnitas u1 ; u2 ; u3:

te 4t 1 t2
 
W1 1
Z Z
u1 D dt D 4t
dt D D t 2:
W 4e 4 2 8
W2 . 2t 2 t/e 2t 1
Z Z Z
u2 D dt D 4t
dt D .2t 2 C t/e 2t dt:
W 4e 4
Si en la última integral aplicamos integración por partes, obtenemos:

1 2t
u2 D e .3 C 6t C 4t 2 /:
16
Otra integración por partes nos da como resultado:

W3 2t e 2t 1 1
Z Z Z
2t 2t
u3 D dt D 4t
dt D te dt D e .2t C 1/:
W 4e 2 8

De x D e t , obtenemos y D ln x. Por lo tanto, la solución general de la ecuación diferencial es

1 2 1 1
y D Œc1 C c2 e 2t C c3 t e 2t  C t C e 2t .3 C 6t C 4t 2 /e 2t e 2t
.2t C 1/t e 2t D
8 16 8
1 3 3 1 1 2 1
D Œc1 C c2 e 2t C c3 t e 2t  C t 2 C C t C t2 t t D
8 16 8 4 4 8
1 1
D Œc1 C c2 e 2t C c3 t e 2t  C t 2 C t:
8 4

De x D e t , obtenemos t D ln x, por lo tanto, la solución general de la ED en la variable original x es

1 1
y D c1 C c2x 2 C c3 x 2 ln x C ln 2 x C ln x:
8 4
Ahora bien, de las condiciones iniciales:

y.1/ D 0 ) c1 C c2 D 0:

Para aplicar la segunda condición requerimos la primera derivada; ésta es

1 ln x 1
y 0 D 2c2 x C c3 Œx C 2x ln x C C :
4 x 4x
Por lo tanto, y 0 .1/ D 1 produce:

1 3
1 D 2c2 C c3 C o bien 2c2 C c3 D :
4 4
Finalmente, para usar la tercera condición, requerimos la segunda derivada; al calcularla encontramos:
 1 1 ln x
 
1
y D 2c2 C c3 Œ3 C 2 ln x C
00
:
4 x2 4x 2

De donde, por la condición y 00.1/ D 1, hallamos:

1 1
1 D 2c2 C 3c3 C o bien 1 D 2c2 C 3c3:
4 4
Resolviendo el sistema, obtenemos: 
c Cc D 0I
1 2


3
2c2 C c3 D I

 4
2c2 C 3c3 D 1:


13 13 7
Encontramos que c1 D ; c2 D & c3 D . Al utilizar estos valores en la solución general, hallamos
16 16 8
la siguiente solución particular:

13 13 2 7 2 1 1
yD C x x ln x C ln 2 x C ln x D
16 16 8 8 4
1
D Œ 13 C 13x 2 14x 2 ln x C 2 ln 2 x C 4 ln x:
16


Ejercicios 4.7.2 Variación de parámetros para ED de orden n. Soluciones en la página 468

Utilice el método de variación de parámetros o el correspondiente a la ED de Cauchy-Euler para proporcionar la
solución de cada una de las siguientes ecuaciones diferenciales. Si se indica, utilice la información proporcionada.

1. y .3/ y 00 D 12x 2 C 6x .

2. y .4/ C y 00 D x 2 C x .

3. y .3/ C 3y 00 C 3y 0 C y D e x
.
2x 3 C x 2 4x 6
4. y .3/
2y 00 0
y C 2y D .
x4
5. y 000 D 2y 00 C 1 .
 6. y .4/ C 16y 00 D 64 cos.4x/ .

7. y .3/ 4y 00 C 4y 0 D 12e 2x C 24x 2 .

8. y .4/ 2y 00 C y D 100 cos 3x .

9. y .3/ 6y 00 C 11y 0 6y D e x .
24.x C y/
10. y .3/ D .
x3
11. x 3 y .3/ x 2 y 00 C 2xy 0 2y D x 3 .

12. x 3 y .3/ C 5x 2 y 00 C 2xy 0 2y D x 4 .

13. x 3 y .3/ 4x 2y 00 C 8xy 0 8y D 4 ln x .

14. x 3 y .3/ C x 2 y 00 6xy 0 C 6y D 30x .

15. xy .3/ C 2xy 00 xy 0 2xy D 1, si el conjunto fundamental de soluciones está integrado por
1 D e x ; 2 D e x
; 3 D e 2x
.

16. x 2 y .3/ 2y 0 D 5 ln x, si el conjunto fundamental de soluciones está integrado por

1 D 1; 2 D ln x & 3 D x 3 .

17. y .3/ y 0 D 2x; con y.0/ D 0; y 0 .0/ D 1; y 00 .0/ D 2 .

18. y .4/ y D 8e x ; con y.0/ D 1; y 0 .0/ D 0; y 00 .0/ D 1; y 000 .0/ D 0 .

19. y .3/ C 3y 00 C 3y 0 C y D 12e x

; con y.0/ D 1; y 0 .0/ D 0; y 00 .0/ D 3 .

20. y .4/ y D cos x; con y.0/ D 1; y 0 .0/ D 1; y 00 .0/ D y 000 .0/ D 0 .