Ecuaciones Diferenciales De Orden n.

Gómez Ávalos, G. (2022) Ecuaciones Diferenciales de Orden n.

Universidad Andrés Bello, Chile.
Transformaciones lineales Página 2

Unidad III: Ecuaciones diferenciales de orden n.

1 Ecuaciones diferenciales lineales: Caso no homogéneas. 3
1.1 Aniquiladores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Método de coecientes indeterminados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3 Método de variación de parámetros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 Reducción de orden. 12
2.1 Fórmulas de Abel y Liouville. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

2
Transformaciones lineales Página 3

1 Ecuaciones diferenciales lineales: Caso no homogéneas.

1.1 Aniquiladores
Deseamos resolver ecuaciones diferenciales del tipo:
f (D)y = h(x).

Donde h(x) es alguna función diferenciable asociada a una variable x en un intervalo I ⊆ R.

La idea es encontrar un operador diferencial g(D), si es que existe, tal que aniquile a h(x), osea
g(D)h(x) = 0.

Luego,
g(D)f (D)y = g(D)h(x) = 0.

La cual es una ecuación homogénea que podemos resolver con las técnicas aprendidas, cuya solución llamamos solución
general ampliada, y a partir de ella determinamos la forma de la solución particular.
Ejemplos 1.

1. Resolver la ecuación diferencial D2 (D − 1)y = 3ex + sin(x).

Solución: Comenzamos determinando la solución de la ecuación diferencial homogénea asociada. Para ello,
determinamos las raíces de la ecuación característica: m2 (m − 1) = 0 ⇒ m = 0 ∨ m = 1.
Luego, la solución homogénea está dada por:
yh = c1 + c2 x + c3 ex .

Ahora, deseamos determinar una solución particular del problema. Así, tenemos que un aniquilador de la función
h(x) = 3ex + sin(x) es g(D) = (D − 1)(D2 + 1), entonces aplicamos dicho operador a la ecuación diferencial:

D2 (D − 1)y = 3ex + sin(x) ⇒ D2 (D − 1)2 (D2 + 1)y = 0.

Esta última ecuación diferencial, denominada general ampliada, es de tipo homogénea, por lo que la solución de
la general ampliada es:
yga = c1 + c2 x + c3 ex + c4 xex + c5 cos(x) + c6 sin(x) .
| {z } | {z }
yh yp

Por tanto, la forma de la solución particular es: yp = c4 xex + c5 cos(x) + c6 sin(x), con c4 , c5 y c6 constantes a
determinar. Para ello, derivamos tres veces yp :
yp = c4 xex + c5 cos(x) + c6 sin(x) ⇒ yp0 = c4 (ex + xex ) − c5 sin(x) + c6 cos(x)
⇒ yp00 = c4 (2ex + xex ) − c5 cos(x) − c6 sin(x)
⇒ yp000 = c4 (3ex + xex ) + c5 sin(x) − c6 cos(x).

y reemplazamos en la ecuación diferencial:

D2 (D − 1)y = 3ex + sin(x) ⇔ (D3 − D2 )y = 3ex + sin(x) ⇔ y 000 − y 00 = 3ex + sin(x).

Obteniendo:
yp000 − yp00 = c4 (3ex + xex ) + c5 sin(x) − c6 cos(x) − (c4 (2ex + xex ) − c5 cos(x) − c6 sin(x))
= c4 ex + (c5 − c6 ) cos(x) + (c5 + c6 ) sin(x)
= 3ex + sin(x).

Comparando términos semejantes formamos el sistema:

c4 =3 c4 =3
1
c5 − c6 =0 ⇔ c5 = c6 ⇒ 2c6 = 1 ⇔ c6 = = c5 .
c5 + c6 =1 c5 + c6 =1 2

3
Transformaciones lineales Página 4

1
Así, yp = 3xex + (cos(x) + sin(x)).
2
Concluyendo que la solución general del problema es:
1
y = yh + yp = c1 + c2 x + c3 ex + 3xex + (cos(x) + sin(x)).
2

2. Resolver la ecuación diferencial (D2 + D − 2)y = 2x − 40 cos(2x).

Solución: Comenzamos determinando la solución de la ecuación diferencial homogénea asociada. Para ello,
determinamos las raíces de la ecuación característica: m2 − m − 2 = 0 ⇔ (m + 2)(m − 1) ⇒ m = −2 ∨ m = 1.
Luego, la solución homogénea está dada por:
yh = c1 e−2x + c2 ex .

Ahora, deseamos determinar una solución particular del problema. Así, tenemos que un aniquilador de la función
h(x) = 2x − 40 cos(2x) es g(D) = D2 (D2 + 4), entonces aplicamos dicho operador a la ecuación diferencial:

(D2 + D − 2)y = 2x − 40 cos(2x) ⇒ (D − 1)(D + 2)D2 (D2 + 4)y = 0.

Esta última ecuación diferencial, denominada general ampliada, es de tipo homogénea, por lo que la solución de
la general ampliada es:
yga = c1 e−2x + c2 ex + c3 + c4 x + c5 cos(2x) + c6 sin(2x) .
| {z } | {z }
yh yp

Por tanto, la forma de la solución particular es: yp = c3 +c4 x+c5 cos(2x)+c6 sin(2x), con c3 , c4 , c5 y c6 constantes
a determinar. Para ello, derivamos dos veces yp :
yp = c3 + c4 x + c5 cos(2x) + c6 sin(2x) ⇒ yp0 = c4 − 2c5 sin(2x) + 2c6 cos(2x)
⇒ yp00 = −4c5 cos(2x) − 4c6 sin(2x).

y reemplazamos en la ecuación diferencial:

(D2 + D − 2)y = 2x − 40 cos(2x) ⇔ y 00 + y 0 − 2y = 2x − 40 cos(2x).

Obteniendo:
y 00 + y 0 − 2y = −4c5 cos(2x) − 4c6 sin(2x) + c4 − 2c5 sin(2x) + 2c6 cos(2x)
− 2(c3 + c4 x + c5 cos(2x) + c6 sin(2x))
= (−2c3 + c4 ) − 2c4 x + (−4c5 + 2c6 − 2c5 ) cos(2x) + (−4c6 − 2c5 − 2c6 ) sin(2x)
= 2x − 40 cos(2x).

Comparando términos semejantes formamos el sistema:

−2c3 + c4 =0 c3 = − 21 c4
−2c4 =2 c4 = −1
⇔ ⇒ 20c6 = −40 ⇔ c6 = −2 ⇒ c5 = 6.
−6c5 + 2c6 = −40 −6c5 + 2c6 = −40
−6c6 − 2c5 =0 c5 = −3c6

Así, yp = 21 − x + 6 cos(2x) − 2 sin(2x).

Concluyendo que la solución general del problema es:
1
y = yh + yp = yh = c1 e−2x + c2 ex + − x + 6 cos(2x) − 2 sin(2x).
2

4
Transformaciones lineales Página 5

1.2 Método de coecientes indeterminados

Existe una alternativa para el método anterior, aniquiladores, el cual propone la forma de la solución particular de
una ecuación diferencial de coecientes constantes:
(an Dn + an−1 Dn−1 + · · · + a1 D + a0 )y = h(x),

di
donde para todo i ∈ {1, 2, . . . , n}, con Di = i y a0 , ai ∈ R constantes, a partir de un resumen que se presenta en la
dx
siguiente tabla:
Función Forma de Multiplicidad k
h(x) = yp = de la raíz en
polinomio caract.

p(x) xk (a0 + a1 x + · · · + am xm ) 0
grado(p) = m

p(x)eax xk (a0 + a1 x + · · · + am xm )eax a

grado(p) = m

c1 cos(bx) + c2 sin(bx) xk (a0 cos(bx) + b0 sin(bx)) bi

p(x) cos(bx) + q(x) sin(bx) xk (a0 + a1 x + · · · + ar xr ) cos(bx)

grado(p) = m +xk (b0 + b1 x + · · · + br xr ) sin(bx) bi
grado(p) = n r = max{n, m}

c1 eax cos(bx) + c2 eax sin(bx) xk (a0 cos(bx) + b0 sin(bx))eax a + bi

p(x)eax cos(bx) + q(x)eax sin(bx) xk (a0 + a1 x + · · · + ar xr )eax cos(bx)

grado(p) = m +xk (b0 + b1 x + · · · + br xr )eax sin(bx) a + bi
grado(p) = n r = max{n, m}

5
Transformaciones lineales Página 6

Ejemplos 2.

1. Resolver la ecuación diferencial y 00 − 3y 0 + 2y = (x2 + x)e3x .

Solución: Comenzamos resolviendo la ecuación homogénea asociada: y 00 − 3y 0 + 2y = 0.

Raíces del polinomio característico: f (m) = m2 − 3m + 2 = 0.
√
2 3± 9−8
m − 3m + 2 = 0 ⇒ m = ⇔ m1 = 2 ∧ m2 = 1.
2
Luego, la solución de la homogénea asociada: yh = c1 e2x + c2 ex .
Aplicamos el método de coecientes constantes para determinar la solución particular de la ecuación diferencial
completa:
h(x) = (x2 + x)e3x ,

donde observando la tabla:

Función Forma de Multiplicidad k
h(x) = yp = de la raíz en
polinomio caract.

p(x)eax xk (a0 + a1 x + · · · + am xm )eax a

grado(p) = m

identicamos p(x) = x2 + x, grado(p) = m = 2 y a = 3, entonces k = 0, pues 3 no es raíz del polinomio

característico. Así, la forma de la solución particular es yp = x0 (a0 + a1 x + a2 x2 )e3x = (a0 + a1 x + a2 x2 )e3x ,
con a0 , a1 y a2 constantes a determinar.
Derivamos hasta obtener la derivada de tercer orden para luego reemplazar en la ecuación diferencial:
yp = (a0 + a1 x + a2 x2 )e3x ⇒ yp0 = ((a1 + 3a0 ) + (2a2 + 3a1 )x + 3x2 )e3x
⇒ yp00 = ((9a0 + 6a1 + 2a2 ) + (12a2 + 9a1 )x + 9a2 x2 )e3x .

Entonces,
yp00 − 3yp0 + 2yp = ((2a0 + 3a1 + 2a2 ) + (2a1 + 6a2 )x + 2a2 x2 )e3x
= (x2 + x)e3x .

Comparando formamos el sistema:

2a0 + 3a1 + 2a2 =0
1
2a1 + 6a2 = 1 ⇔ a0 = 1, a1 = −1 y a2 = .
2a2 =1 2

x2
 
Así, una solución particular de la ecuación diferencial es: yp = 1−x− e3x .
2
Por lo tanto, la solución general está dada por:
x2
 
2x x
y = yh + yp = c1 e + c2 e + 1 − x − e3x .
2

6
Transformaciones lineales Página 7

2. Resolver la ecuación diferencial y 00 + y = 3 sin(x).

Solución: Comenzamos resolviendo la ecuación homogénea asociada: y 00 + y = 0.

Raíces del polinomio característico: f (m) = m2 + 1 = 0:
m2 + 1 = 0 ⇒ m = ±i.

Luego, la solución de la homogénea asociada: yh = c1 cos(x) + c2 sin(x).

Aplicamos el método de coecientes constantes para determinar la solución particular de la ecuación diferencial
completa:
h(x) = 3 sin(x),

donde observando la tabla:

Función Forma de Multiplicidad k
h(x) = yp = de la raíz en
polinomio caract.

c1 cos(bx) + c2 sin(bx) xk (a0 cos(bx) + b0 sin(bx)) bi

identicamos y b = 1, entonces k = 1, pues bi = i es raíz del polinomio característico con multiplicidad 1. Así,
la forma de la solución particular es yp = x(a0 cos(x) + b0 sin(x)), con a0 y b0 constantes a determinar.
Derivamos hasta obtener la derivada de tercer orden para luego reemplazar en la ecuación diferencial:
yp = a0 x cos(x) + b0 x sin(x) ⇒ yp0 = a0 cos(x) + b0 sin(x) − a0 x sin(x) + b0 x cos(x)
⇒ yp00 = 2b0 cos(x) − 2a0 sin(x) − a0 x cos(x) − b0 x sin(x).

Entonces,
yp00 + yp = 2b0 cos(x) − 2a0 sin(x)
= 3 sin(x).

Comparando formamos el sistema:

2b0 =0 3
⇔ a0 = − ∧ b0 = 0.
−2a0 =3 2

3
Así, una solución particular de la ecuación diferencial es: yp = − x cos(x).
2
Por lo tanto, la solución general está dada por:
3
y = yh + yp = c1 cos(x) + c2 sin(x) − x cos(x).
2

7
Transformaciones lineales Página 8

1.3 Método de variación de parámetros

Consideremos una ecuación diferencial de orden n completa:
an (x)y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x)y = b(x), an (x) 6= 0, ∀x ∈ I = (a, b). (1)
Siendo b(x) y ai (x), con i ∈ {1, 2, . . . , n}, funciones continuas en un intervalo I = (a, b).
El método propone como solución particular:
yp = c1 (x)y1 + · · · + cn (x)yn , (2)
donde yh = c1 y1 + · · · + cn yn es la solución general de la ecuación diferencial homogénea asociada:
an (x)y (n) + an−1 (x)y (n−1) + · · · + a1 (x)y 0 + a0 (x)y = 0.

Para encontrar las funciones ci (x), i ∈ {1, · · · , n} de (1) debemos obtener un sistema de n ecuaciones. Así, debemos
imponer n condiciones a través de un sistema de n ecuaciones de la forma:
yp0 = c01 (x)y1 + · · · + c0n (x)yn +c1 (x)y10 + · · · + cn (x)yn0 .
| {z }
=0

Luego,
yp00 = c01 (x)y10 + · · · + c0n (x)yn0 +c1 (x)y100 + · · · + cn (x)yn00 .
| {z }
=0

Así, sucesivamente realizamos este proceso (n − 1) veces:

(n−2) (n−1)
yp(n−1) = c01 (x)y1 + · · · + c0n (x)yn(n−2) +c1 (x)y1 + · · · + cn (x)yn(n−1) .
| {z }
=0

Por lo tanto, para que yp sea una solución particular de (1) debemos tener que:
(n−1)
c01 (x)y1 + · · · + c0n (x)yn(n−1) = b(x).

Finalmente, hemos formado el sistema:

c01 (x)y1 + · · · + c0n (x)yn =0
c01 (x)y10 + · · · + c0n (x)yn0 =0
..
.
(n−2) (n−2)
c01 (x)y1 + · · · + c0n (x)yn =0
(n−1) (n−1)
c01 (x)y1 + · · · + c0n (x)yn = b(x)

8
Transformaciones lineales Página 9

Ejemplo 3. Resolver la ED y 00 + y = tan(x).

Solución: Comenzamos resolviendo la ecuación diferencial homogénea asociada: y 00 + y = 0, determinando las raíces
del polinomio caracaterístico:
m2 + 1 = 0 ⇔ m = ±i.

Luego, la solución homogénea está dada por: yh = c1 cos(x) + c2 sin(x).

Ahora, por el método de variación de parámetros se propone una solución de la forma yp = c1 (x) cos(x) + c2 (x) sin(x),
con c1 (x) y c2 (x) dos funciones diferenciables a determinar. Para ello, debemos resolver el sistema:
c01 cos(x) + c02 sin(x) = 0
.
−c01 sin(x) + c02 cos(x) = tan(x)

A partir de la primera ecuación del sistema se tiene que:

sin(x)
c01 = −c02
cos(x)

Así, reemplazando en la segunda ecuación:

sin2 (x)
−c01 sin(x) + c02 cos(x) = tan(x) ⇒ c02 + c02 cos(x) = tan(x)
cos(x)
 2
sin (x) + cos2 (x)

⇒ c02 = tan(x)
cos(x)
 
1 sin(x)
⇒ c02 =
cos(x) cos(x)
⇒ c02 = sin(x)
Z
⇒ c2 (x) = sin(x)dx

⇒ c2 (x) = − cos(x).

Además,
sin(x) sin2 (x) 1 − cos2 (x) 1
c01 = −c02 =− =− = cos(x) − = cos(x) − sec(x).
cos(x) cos(x) cos(x) cos(x)

Integrando tenemos:
Z
c1 (x) = cos(x) − sec(x)dx ⇒ c1 (x) = sin(x) − ln(| sec(x) + tan(x)|).

Por lo tanto, la solución general del problema está dada por:

y = c1 cos(x) + c2 sin(x) + (sin(x) − ln(| sec(x) + tan(x)|)) cos(x) − cos(x) sin(x)
= c1 cos(x) + c2 sin(x) − ln(| sec(x) + tan(x)|) cos(x).

9
Transformaciones lineales Página 10

Observación 4. Ahora, presentaremos una técnica útil que nos ayudará a resolver el sistema:
c01 (x)y1 + · · · + c0n (x)yn =0
c01 (x)y10 + · · · + c0n (x)yn0 =0
..
. . (3)
(n−2) (n−2)
c01 (x)y1 + ··· + c0n (x)yn =0
(n−1) (n−1)
c01 (x)y1 + ··· + c0n (x)yn = b(x)

Al ser un sistema lineal de n ecuaciones y n incógnitas podemos resolver utilizando matrices. Así, si denimos la
matriz de coecientes:
···
 
y1 y2 yn
 y10 y20 ··· 0
yn 
A =  .. .. .. ..  ∈ Mn×n ,
 
 . . . . 
(n−1) (n−1) (n−1)
y1 y2 ··· yn

el vector columna:
c01 (x)
 
 c02 (x) 
X =  .  ∈ Mn×1 .
 
 .. 
c0n (x)

y el vector:

0
 0 
X =  .  ∈ Mn×1 .
 
 .. 
b(x)

Podemos reescribir el sistema (3) como la ecuación matricial:

A · X = B.

El cual tiene única solución si |A| =

6 0, donde |A| representa el determinante de la matriz de coecientes. Recordar el
método de Cramer para resolver sistemas de ecuaciones lineales.
Además, notamos que |A| = 6 0 es equivalente a:

y1 y2 ··· yn
y10 y20 ··· yn0
W [y1 , y2 , . . . , yn ] = .. .. .. .. 6= 0.
. . . .
(n−1) (n−1) (n−1)
y1 y2 ··· yn

Es decir, el sistema tendrá única solución si el Wronskiano de las funciones y1 , y2 , . . . , yn es no nulo, el cual por
simplicidad denotaremos por W , y las soluciones serán dadas por el método de Cramer. Por ejemplo,
0 y2 ··· yn
1 0 y20 ··· yn0
c01 (x) = .. .. .. .. .
W . . . .
(n−1) (n−1)
b(x) y2 ··· yn

10
Transformaciones lineales Página 11

Ahora, veamos un ejemplo:

Ejemplo 5. Resolver la ecuación diferencial y 00 + y = csc(x).
Solución: Comenzamos resolviendo la ecuación diferencial homogénea asociada: y 00 + y = 0, determinando las raíces
del polinomio caracaterístico:
m2 + 1 = 0 ⇔ m = ±i.

Luego, la solución homogénea está dada por: yh = c1 cos(x) + c2 sin(x).

Ahora, por el método de variación de parámetros se propone una solución de la forma yp = c1 (x) cos(x) + c2 (x) sin(x),
con c1 (x) y c2 (x) dos funciones diferenciables a determinar. Para ello, debemos resolver el sistema:
c01 cos(x) + c02 sin(x) = 0
.
−c01 sin(x) + c02 cos(x) = csc(x)

Donde:
cos(x) sin(x)
W [cos(x), sin(x)] = W = = cos2 (x) + sin2 (x) = 1.
− sin(x) cos(x)

Luego,
1 0 sin(x)
c01 (x) = = − sin(x) csc(x) = −1.
W csc(x) cos(x)

Integrando:
Z
c1 (x) = − dx = −x + k1 .

Además,
1 cos(x) 0 cos(x)
c02 (x) = = = cot(x).
W − sin(x) csc(x) sin(x)

Integrando:
Z
c1 (x) = cot(x)dx

= ln(| sin(x)|) + k2 .

Por lo tanto, la solución general del problema está dada por:

y = c1 cos(x) + c2 sin(x) + ln(| sin(x)|) cos(x) − x sin(x).

11
Transformaciones lineales Página 12

2 Reducción de orden.
Hasta el momento solo hemos estudiado ecuaciones diferenciales lineales de orden n con coecientes constantes. Para
el caso de orden 2, donde los coecientes son funciones no constantes, puede resultar muy complejo determinar una
solución general, pero si conocemos una solución podemos determinar otra linealmente independiente a través de la
técnica de reducción de orden. Así, determinar la solución general del sistema.

Dada una ecuación diferencial de segundo orden de la forma:

y 00 + p(x)y 0 + q(x)y = 0, (4)
supongamos que conocemos una solución particular y1 (x). Entonces, es posible encontrar una segunda solución
linealmente independiente de (4), la cual denotaremos por y2 (x), mediante el siguiente método:
Supongamos que y2 = v(x)y1 , donde v(x) debe ser determinada.
y2 = v(x)y1 ⇒ y20 = v 0 y1 + vy10
⇒ y200 = v 00 y1 + 2v 0 y10 + vy100 .

Luego, sustituyendo en la ecuación (4) tenemos:

y200 + py20 + qy2 = 0
(v 00 y1 + 2v 0 y10 + vy100 ) + p(v 0 y1 + vy10 ) + qvy1 = 0
y1 v 00 + (2y10 + py1 )v 0 + (y100 + py10 + qy1 )v = 0.

Así, como y1 es solución de (4) se tiene:

y100 + p(x)y10 + q(x)y1 = 0,

entonces
y1 v 00 + (2y10 + py1 )v 0 = 0. (5)
Ahora, para resolver la ecuación diferencial anterior (5) hacemos el cambio de variable w = v 0 , entonces:
y1 v 00 + (2y10 + py1 )v 0 = 0 ⇒ y1 w0 + (2y10 + py1 )w = 0
dw 2y 0 + py1
⇒ =− 1 dx
w y1
2y10 + py1
Z
⇒ ln(w) = − dx
y1
0 +py
2y1 1
R
⇒ w = e− y1 dx

0 +py
2y1 1
R
⇒ v 0 = e− y1 dx
Z 0 +py
2y1 1
R
⇒ v = e− y1 dx
dx.

Notando que:
R
0 +py 0 p(x)dx
−
R 2y1 1
dx −2
R y1 R
dx− pdx ln(y1−2 )
R
p(x)dx e
e y1
=e y1
=e e = .
y12

Por lo tanto, tenemos que:

R
0 +py
e− p(x)dx
Z 2y1
Z
1
R
− dx
y2 = vy1 = y1 e y1
dx = y1 dx.
y12

12
Transformaciones lineales Página 13

Ejemplo 6. Si y1 = x es solución de la ecuación diferencial x2 y 00 + 3xy 0 − 3y = 0. Determine su solución general

para x ∈ (0, +∞).
Solución. Supongamos que y2 = v(x)x es otra solución de x2 y 00 + 3xy 0 − 3y = 0. Luego, derivando tenemos:
y2 = v(x)x ⇒ y20 = v 0 x + v
⇒ y200 = v 00 x + 2v 0 .

Reemplazando en la ecuación diferencial:

x2 y200 + 3xy20 − 3y2 = 0
v 00 x3 + 2v 0 x2 + 3x(v 0 x + v) − 3vx = 0
v 00 x3 + 5x2 v 0 = 0

Haciendo w = v 0 , entonces tenemos que w0 = v 0 y

w0 x3 + 5x2 w = 0
w0 x3 = −5x2 w /dividiendo por x3 , x > 0.
5
w0 = − w.
x
Esta última ecuación es de variables separables, entonces:
5 dw 5
w0 = − w ⇒ = − dx
x w x
⇒ ln(w) = ln(x−5 )
⇒ w = x−5
⇒ v 0 = w = x−5
1
⇒ v = − 4.
4x
1
Entonces, y2 = v(x) · x = − . Así, la solución general de la ecuación diferencial es:
4x3
−1
y = c1 x + c2 · , con c1 y c2 constantes reales arbitrarias.
4x3
1
Además, observemos que c2 · − sigue siendo una constante arbitraria, entonces podemos decir que la solución general
4
de la ecuación diferencial está dada por:
1
y = c1 x + c2 · , con c1 y c2 constantes reales arbitrarias.
x3

13
Transformaciones lineales Página 14

2.1 Fórmulas de Abel y Liouville.

Los siguientes resultados nos permiten resumir la técnica de reducción de orden para ecuaciones diferenciales de la
forma (4).

Teorema 1 (fórmula de Abel).

El Wronskiano de dos soluciones y1 e y2 , denotado por W [y1 , y2 ], de la ecuación diferencial homogénea

y 00 + p(x)y 0 + q(x)y = 0 es una solución deR la ecuación lineal de primer orden W 0 + p(x)W = 0. Por tanto,
existe una constante c tal que W (x) = ce p(x)dx .

Demostración. Por denición:

y1 y2
W = W [y1 , y2 ] = = y1 y20 − y10 y2 ⇒ W 0 = y10 y20 + y1 y200 − y100 y2 − y10 y20
y10 y20
⇒ W 0 = y1 y200 − y100 y2 .

Entonces,
W 0 + p(x)W = y1 y200 − y100 y2 + p(x)[y1 y20 − y10 y2 ]
= y1 [y200 + p(x)] − y20 [y100 + p(x)y10 ]
= y1 [−q(x)y2 ] − y2 [q(x)y1 ]
= 0.

El siguiente resultado resume la técnica de resolución de orden.

Teorema 2 (fórmula de Liouville).

Si conocemos una función y1 de la ecuación diferencial:

y 00 + p(x)y 0 + q(x)y = 0, (6)
R
e− p(x)dx
Z
Entonces, la función y2 = y1 dx es una solución de (6). Además, la solución general de (6) es:
y12

y = c1 y1 + c2 y2 , c1 , c2 ∈ R constantes arbitrarias.

Ejemplos 7.

1. Si y1 = x es solución de la ecuación diferencial x2 y 00 − 3xy 0 + 4y = 0. Determinar su solución general para

x ∈ (0, +∞).

Solución. Como x > 0, entonces:

3 0 4
x2 y 00 − 3xy 0 + 4y = 0 ⇒ y 00 − y + 2 y = 0.
x x
Aquí, vamos a utilizar la fórmula de Liouville para determinar una solución y2 . Para ello, identicamos p(x) =
3
− , entonces integrando:
x
Z Z
3
p(x) dx = − dx = ln(x−3 ),
x

14
Transformaciones lineales Página 15

donde hemos considerado la constante de integración igual a cero.

Luego,
R
e− p(x)dx
Z
y2 = y1 dx
y12
−3
e− ln(x )
Z
=x dx
x2
3
eln(x )
Z
=x dx
x2
Z 3
x
=x dx
x2
x2 x3
Z
= x x dx = x · = .
2 2
Por lo tanto, la solución general de la ecuación diferencial está dada por:
y = c1 x + c2 x3 , c1 , c2 ∈ R constantes arbitrarias.

15
Transformaciones lineales Página 16

Referencias bibliográcas.
1. Blanchard, P., Devaney, R.& Hall, G. (1998). Ecuaciones Diferenciales. International Thompson Editores,
Argentina.
2. Bronson, R. (1988). Ecuaciones diferenciales modernas: Teoría y problemas. MacGraw-Hill, México.
3. Campbell, S. & Haberman, R. (1998). Introducción a las Ecuaciones Diferenciales con problemas de valor en la
frontera. MacGraw-Hill, México.
4. Nagle, R. & et. al (2005). Ecuaciones diferenciales y problemas de valor en la frontera (4a edición). Person,
México.
5. Zill, D. & Cullen, M. (2008). Matemáticas avanzadas para ingeniería, Vol. 1: Ecuaciones diferenciales (3a
edición). MacGraw-Hill, México.

16