240 Ecuaciones diferenciales

4. y 00 C 4y D .5x 2 C 4x 3/e x .
5. 4y 00 4y 0 C y D . 3x C 4/ sen x 2.2x 1/ cos x .
6. y 00 C y 0 D 3 .

7. y 00 C 4y 0 C 4y D 2.x C 3/ .
8. y 00 y0 12y D e 4x .
9. y 00 C 25y D 6 sen x .
10. y 00 2y 0 C 5y D e x sen x .
11. y 00 y D x2ex .
12. y 00 C 25y D 20 sen 5x C x cos x .
13. y 00 y 0 D .3x 2 C 6x C 1/e x .
14. y 00 2y 0 8y D 2 sen x C .12x 10/e 4x .
15. y 00 4y 0 C 4y D .18x 4/e 2x .
16. y 00 C 12y 0 C 100y D 48 sen 10t .
17. y 00 5y 0 C 6y D 2 sen 2x cos 2x C .3 2x/e 2x .

4.7 Variacin de parmetros

4.7.1 Variacin de parmetros para ED de orden 2
El mtodo que presentamos en esta seccin, llamado de variacin de parmetros, es un procedimiento til
y ms general que el de coeficientes indeterminados para la obtencin de una solucin particular yp .x/ de
la ED lineal no homognea y se basa en el conocimiento de la solucin general de la ED lineal homognea
asociada a la ED original.
Presentamos primeramente este mtodo para las ED lineales de segundo orden. Diremos que el mtodo de
variacin de parmetros se usa para obtener una solucin 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 solucin general de la ED lineal homognea asociada

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

Si suponemos que la solucin general de la ED lineal homognea (4.22) est dada por la combinacin 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 adems 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 solucin general de y 00 C p.x/y 0 C q.x/y D 0.
El mtodo de variacin de parmetros propone que la solucin 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.

4.7 Variacin de parmetros 241

Para determinarlas, buscaremos qu condiciones deben cumplir para que la funcin propuesta yp sea solu-
cin de la ED (4.21). Lo primero que necesitamos hacer es calcular las derivadas yp0 , yp00 para utilizarlas en
(4.21) y ver cules 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 expresin de yp00 no aparezcan u100 & u200 , ya que la inclusin de
estas segundas derivadas en yp00 hara mucho ms compleja la obtencin de las funciones u1 & u2 .
Luego entonces, usando la condicin 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 solucin de la ED lineal y 00 C p.x/y 0 C q.x/y D g.x/ slo si se cumple que

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

Esta condicin, 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 homognea 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 incgnitas:
(
u10 1 C u20 2 D 0I
(4.24)
u10 10 C u20 20 D g.x/ :

El determinante s del sistema es  

 2 
s D  10 D W .1 ; 2/:
1 20 
Y debido a que el wronskiano W .1 ; 2 /.x/ 0, entonces s 0, por lo que el sistema de ecuaciones tiene
una nica solucin. Dicha solucin nica la calculamos usando la regla de Cramer:
 
 0 2 
g.x/  0 

2 g.x/2 g.x/2 .x/
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//
242 Ecuaciones diferenciales

De donde obtenemos u1 & u2 mediante integracin

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 solucin particular de la lineal es
Z Z
g.x/2 .x/ g.x/1 .x/
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 solucin 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 mtodo de variacin de parmetros, calcular una solucin particular y escribir la solu-
cin 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 homognea 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 solucin 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 2 0
8 C 6/ C 3x u2 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
0 2 0
2xu1 C 3x u2 D x 3 :

El determinante del sistema es  2

x 3 

x
W D  D 3x 4 2x 4 D x 4 :
2x 3x 2
4.7 Variacin de parmetros 243

La solucin del sistema es

x 3 
 
 0
 3
x 3x 2  1
u10 D D D x 4
:
 2 W  x4
x 0 
x 3

2x x 1
u20 D D Dx 5
:
W x4
Integrando:
Z 3
4 x 1 3
u1 D x dx D C C1 D x C C1 :
3 3
Z 4
5 x 1 4
u2 D x dx D C C2 D x C C2 :
4 4
1 1
Tomando C1 D 0 y C2 D 0, obtenemos, u1 D x 3
y u2 D x 4
. Por ello, una solucin particular es
3 4
1 1 1 1 1
yp D u1 x 2 C u2 x 3 D x 3 2
x x 4 3
x D x 1
x 1
D x 1
:
3 4 3 4 12
Entonces la solucin 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 mtodo de variacin de parmetros, calcular una solucin particular y escribir la solu-
cin 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 homognea
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 solucin particular propuesta por este mtodo:
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 condicin

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
244 Ecuaciones diferenciales

(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
0 0
u1 C .1 C ln x/u2 D 4x 1 ln x:

El determinante del sistema es

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

La solucin del sistema es

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

dx 4
Z Z
1
u1 D 4 .ln x/2 dx D 4 .ln x/2
x D .ln x/3 C C1 :
x 3
dx
Z Z
1
u2 D 4 x 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 solucin 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 solucin 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 mtodo de variacin de parmetros, calcular una solucin particular y escribir la solu-
cin general de la ED lineal
y 00 C y D sec2 x:

H Primero se obtiene un conjunto fundamental de soluciones para la ED homognea asociada

y 00 C y D 0:

Proponiendo y D e r x para resolver la ED homognea:

p
r2 C 1 D 0 ) r D 1 D 0 1i:
4.7 Variacin de parmetros 245

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 homognea asociada.
Se propone como solucin 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 solucin 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 solucin 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:
Entonces la solucin general es
y D .sen x/ ln.sec x C tan x/ 1 C c1 sen x C c2 cos x:

246 Ecuaciones diferenciales

Ejemplo 4.7.4 Utilizando el mtodo de variacin de parmetros, calcular una solucin particular y escribir la solu-
cin 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 homognea asociada
y 00 3y 0 C 2y D 0I
proponemos y D e r x 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 solucin 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 C u1 e x C 2u20 e 2x C 4u2 e 2x / 3.u1 e x C 2u2 e 2x / C 2.u1 e x C u2 e 2x / D )
1 C ex
e 3x
) u10 e x C u1 e x .1 3 C 2/ C 2u20 e 2x C u2 e 2x .4 6 C 2/ D )
1 C ex
e 3x
) u10 e x C 2u20 e 2x D : (4.33)
1 C ex
Entonces u10 & u20 satisfacen el sistema formado por (4.32) y (4.33):

e x u10 C e 2x u20 D 0:
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 solucin del sistema es
 
 0
 e 2x 
 e 3x 
 2e 2x 
 1 C ex e 2x e 3x e 2x
u10 D D D :
W e 3x .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 x
D :
W e .1 C e / 1 C ex
4.7 Variacin de parmetros 247

Integrando:
e 2x ex
Z Z
u1 D dx & u2 D dx:
1 C ex 1 C ex
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 solucin 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 solucin 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 seccin, una solucin 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 homognea 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 mtodo de Cramer. An ms:

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 estn dadas mediante un determinante en el que se sustituye la columna 1 & 2, respectivamente, de
W por la columna:
 
0
:
g.x/
Es conveniente sealar las fortalezas y debilidades de este mtodo en comparacin con el mtodo de coefi-
cientes indeterminados, tratado en la seccin anterior.

1. El mtodo de coeficientes indeterminados es en muchos casos ms sencillo y fcil de aplicar que el

de variacin de parmetros, pero tiene esta limitacin: slo es aplicable cuando g.x/ tiene la forma
de un polinomio, exponencial, combinacin lineal de senos y cosenos o una suma de las funciones
mencionadas. Adems la ED lineal debe tener coeficientes constantes.
248 Ecuaciones diferenciales

2. El mtodo de variacin de parmetros en cambio es aplicable en pricipio para cualquier funcin g.x/
y, en este sentido, es mucho ms general. Adems 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 mtodo tiene un costo, puesto que para obtener una
solucin 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 resolucin de estas integrales es dificultosa.)

Ejercicios 4.7.1 Variacin de parmetros para ED de orden 2. Soluciones en la pgina 467

Utilizando variacin de parmetros, calcular una solucin particular y escribir la solucin general de la ecuacin
diferencial dada. Considerar que las funciones y1 D y1 .x/ & y2 D y2 .x/ forman un conjunto fundamental de
soluciones para la ecuacin homognea 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 I y1 D x 3
& y2 D x 4
.
x2
5. x 2 y 00 6xy 0 C 10y D 4x ln x 5xI y1 D x 5 & y2 D x 2 .

Utilizando variacin de parmetros, determinar una solucin particular y escribir la solucin general de la ecuacin
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

4.7.2 Variacin de parmetros para ED de orden n

Descripcin del mtodo general
Una vez discutido el mtodo de variacin de parmetros para ecuaciones diferenciales de orden 2, en esta
seccin extenderemos dicho mtodo a ecuaciones diferenciales de orden n para n > 2. As, consideraremos
el caso de la ED lineal no homognea normalizada:

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


C a0 .x/y D g.x/I (4.34)
4.7 Variacin de parmetros 249

cuya ED homognea 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 solucin general de la ecuacin (4.35),

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


C cn n .x/: (4.36)

En lo que sigue, supondremos tambin 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/ 
 10 .x/ 20 .x/


n0 .x/ 
 
W 1 .x/; 2 .x/;


; n .x/ D  :: :: :: ::  0, para todo x 2 .; / :

 .n :1/ : : :
 
.n 1/ .n 1/


1 .x/  2 .x/


 n.x/

Como vimos en la seccin anterior, el mtodo se apoya en la idea de que las constantes c1 ; c2;


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


; un .x/g
por conocer, que permiten generar una solucin particular

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


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

de la ecuacin (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 ecuacin (4.34).
La idea central de este mtodo 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 relacin 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 un n0 :
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 razn que expusimos anteriormente, ahora planteamos la condicin:

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 un n00 : (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)
250 Ecuaciones diferenciales

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 1/ C


C un0 n.n 1/

Can 1 u1 1.n
C


Cu  1/
C


C .n 1/
n n

.n/ .n 1/
yp yp

C a0 u1 1 C


C un n D g.x/: (4.43)

yp

Ahora reacomodamos la ecuacin (4.43) factorizando u1 , u2


un :

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


C a0 1 C


C un n.n/ C an .n 1/
1 n 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 ecuacin
(4.35), todos los trminos en (4.44) son iguales a cero con excepcin del ltimo; esto permite llegar a la con-
clusin 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 ecuacin (4.45), concluimos que las funciones
incgnitas u10 ; : : : ; un0 satisfacen a las condiciones:



u10 1 C u20 2 C


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


C D

u
1 1
 u 
2 2 u 
n n 0I
::

: (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 incgnitas u10 ; : : : ; un0 , podemos
utilizar la regla de Cramer para obtener la solucin 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 incgnitas u1 .x/; : : : ; un .x/.
Finalmente observe que, si en lugar de la ecuacin diferencial (4.34), tuviramos
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 idnticamente igual a 1, entonces habra que poner g.x/ D en lugar de
an .x/
g.x/ en la ltima ecuacin de (4.46). En otras palabras, habra que normalizar la ED.
4.7 Variacin de parmetros 251

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

H Primero observamos que esta ecuacin diferencial no puede ser resuelta por medio del mtodo de
coeficientes indeterminados debido a la presencia de la funcin cotangente. La ED est normalizada, ya
que el coeficiente de la mayor derivada de la ecuacin es igual a 1.
Para la solucin, determinamos en primer lugar la ecuacin caracterstica correspondiente a la ED asociada:

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

Las races de esta ecuacin algebraica son r1 D 0I r2;3 D 2i . Por lo tanto, el conjunto fundamental de
soluciones est integrado por las funciones

1 D 1; 2 D cos 2x & 3 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 8sen 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 2x sen 2x 
 
 0
2 sen 2x 2 cos 2x  D cot 2x2 cos 2 2x C 2 sen 2 2x D 2 cot 2xI

W1 D  0
cot 2x 4 cos 2x 4 sen 2x 

sen 2x 
 
1 0
2 cos 2x  D cot 2x2 cos 2x D 2 cot 2x cos 2xI

W2 D 0 0
0 cot 2x 4 sen 2x 

y tambin:

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
Z Z  
W1 1 1 1
u1 D dx D dx D ln.sen 2x/ D ln.sen 2x/:
W 8 4 2 8
Para u2 ,

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

Z Z Z Z
W2 1 1
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
252 Ecuaciones diferenciales

Estamos ahora listos para obtener la solucin general la cual, como sabemos, resultar sumando la solucin
complementaria y la solucin particular que calculamos con el mtodo de variacin de parmetros. De esta
manera la solucin 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 1
C ln.sen 2x/ ln.csc 2x cot 2x/ cos 2x cos 2 2x sen 2 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 ecuacin diferencial de Cauchy-Euler

El ejemplo anterior muestra cmo generar una solucin particular y la consecuente solucin general de
la ED por el mtodo de variacin de parmetros si tan slo se conoce un conjunto fundamental de solu-
ciones de la ecuacin diferencial homognea asociada. Hay que decir que tuvimos a nuestro favor que la
ecuacin homognea es de coeficientes constantes, caracterstica que facilit la determinacin del conjunto
fundamental.
Si los coeficientes no son constantes, la tarea puede resultar ms compleja con excepcin de algunos casos
particulares, como el que discutimos a continuacin.

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

d ny n 1d
n 1
y dy
an x n C an 1x 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 3
C a2 x 2
C a1 x C a0 y D g.x/: (4.48)
dx dx 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 0 d
e t dy t d 2y t dy
d 2y dy 0 e e 
d 2y dy

y 00 D D D dt D dt dt
D dt 2 dt D e 2t
:
dx 2 dx dx d t et dt 2 dt
.e /
dt dt
4.7 Variacin de parmetros 253

Es decir,

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

Para la tercera derivada encontramos que

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

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

Es decir,

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

Al usar (4.49), (4.50) y (4.51) en (4.48), y tomando en consideracin 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 C 2 C 3 3 3 D te t ;
dt 3 dt 2 dt dt 2 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 ecuacin caracterstica asociada a la ED homognea es

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

cuyas races 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 
254 Ecuaciones diferenciales

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 

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

Por ltimo, integramos por partes para determinar las funciones incgnitas 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 3t
Z Z Z
u3 D dt D dt D te 3t dt D e .3t 1/:
W 16 8 72
En conclusin, la solucin 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 ecuacin 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 ecuacin por x a fin de llevarla a una ecuacin 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 2y dy dy
3 C2 C C D t;
dt 3 dt 2 dt dt 2 dt dt
4.7 Variacin de parmetros 255

o bien:
d 3y d 2y dy
4 C4 D t:
dt 3 dt 2 dt
sta es una ED con coeficientes constantes. La ecuacin caracterstica correspondiente a la ED homognea
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 races de la ecuacin son r1 D 0 y r2 D r3 D 2. En consecuencia, las funciones que integran
el conjunto fundamental de soluciones de la ED homognea 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 te 2t
 

2t
.2t C 1/e  D .8t C 8/e 4t
2t 
.8t C 4/e 4t D 4e 4t :
 
W D W .1 ; 2; 3 / D 0 2e

0 4e 2t .4t C 4/e 2t 

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

te 2t
 
1 0 
0 .2t C 1/e  D t.2t C 1/e 2t D . 2t 2 t/e 2t ;
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 mtodo de variacin de parmetros, las funciones incgnitas u1 ; u2 ; u3:

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

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

2te 2t
Z Z Z
W3 1 2t 1 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 solucin general de la ecuacin diferencial es

1 1 2t 1
y D c1 C c2 e 2t C c3 te 2t C t 2 C e .3 C 6t C 4t 2 /e 2t e 2t
.2t C 1/te 2t D
8 16 8
1 3 3 1 1 2 1
D c1 C c2 e 2t C c3 te 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 te 2t C t 2 C t:
8 4
256 Ecuaciones diferenciales

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

1 2 1
y D c1 C c2x 2 C c3 x 2 ln x C ln 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 condicin 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 condicin, requerimos la segunda derivada; al calcularla encontramos:

1 1 ln x
 
1
y 00 D 2c2 C c3 3 C 2 ln x C 2
:
4 x 4x 2

De donde, por la condicin 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 C c2 D 0I
1


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 solucin general, hallamos
16 16 8
la siguiente solucin particular:

13 13 7 2 1 1
yD C x2 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 Variacin de parmetros para ED de orden n. Soluciones en la pgina 468

Utilice el mtodo de variacin de parmetros o el correspondiente a la ED de Cauchy-Euler para proporcionar la
solucin de cada una de las siguientes ecuaciones diferenciales. Si se indica, utilice la informacin 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 y 0 C 2y D .
x4
5. y 000 D 2y 00 C 1 .
4.8 Variacin de parmetros y reduccin de orden 257

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
x
1 D e ; 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 .

4.8 Variacin de parmetros y reduccin de orden

Con lo anteriormente tratado hasta aqu, podemos hacer las afirmaciones siguientes:

1. Con el mtodo de variacin de parmetros podemos resolver la ED lineal no homognea:

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

siempre y cuando conozcamos la solucin general .x/ D c1 1 .x/ C c2 2 .x/ de la ED lineal ho-
mognea asociada y 00 C p.x/y 0 C q.x/y D 0.

2. Con el mtodo de variacin de parmetros podemos resolver la ED lineal con coeficientes constantes

ay 00 C by 0 C cy D g.x/:

3. Aplicando primero el mtodo de reduccin de orden y luego el mtodo de variacin de parmetros,

podemos resolver la ED lineal
y 00 C p.x/y 0 C q.x/y D g.x/;
con el conocimiento de slo una solucin y D 1 .x/ de la ED lineal homognea asociada

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

Sobre las dos primeras afirmaciones ya hemos ejemplificado; y ahora lo haremos sobre la ltima afirmacin.