UNIVERSIDAD NACIONAL TECNOLÓGICA DE LIMA SUR

IS, IM, IE y IA

ASIGNATURA: INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES

Dr. Jesús Rule Flores Cruz

2025 - 0
UNIDAD DIDÁCTICA 2 SEMANA 7 SESIÓN 1

TEMA DE SESIÓN: ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

APRENDIZAJES ESPERADOS:

• Al finalizar la semana, el estudiante discute e interpreta los ciertos fenómenos de la naturaleza y los modela mediante
ecuaciones diferenciales de orden superior y encuentra su solución.

CAPACIDAD GENERAL:

• Al finalizar el curso, el estudiante aplica las series y ecuaciones diferenciales para resolver problemas relacionados a su
carrera profesional, apoyándose estratégicamente en herramientas matemáticas con responsabilidad y trabajo en equipo.

CAPACIDAD ESPECÍFICA:

• Resuelve situaciones problemáticas de contexto real referidas a analizar cambios discontinuos o regularidades, entre
fenómenos naturales, valores o expresiones; traduciéndolas a modelos matemáticos en ecuaciones diferenciales que mejor
se ajusten al comportamiento del fenómeno observado.
Ecuaciones lineales: teoría básica
Un problema de valor inicial de n-ésimo orden consiste en resolver la
EDO lineal:
n n −1
d y d dy
an ( x) n + an−1 ( x) n−1 +  + a1 ( x) + a0 ( x) y = g ( x)
dx dx dx
sujeta a las n condiciones iniciales:

y ( x0 ) = y0 , y( x0 ) = y1 ,  , y ( n−1) ( x0 ) = yn−1

Resolverlo consiste en encontrar una función y(x) en definida en un
intervalo I que contiene a x0, donde se cumplen la ecuación y las
condiciones iniciales.
3
La siguiente EDO lineal de orden n:
dny d n−1 y dy
an ( x) n + an−1 ( x) n−1 +  + a1 ( x) + a0 ( x) y = g ( x)
dx dx dx
Se dice que es no homogénea, si 𝑔 𝑥 ≠ 0.

dny d n−1 y dy
an ( x) n + an−1 ( x) n−1 +  + a1 ( x) + a0 ( x) y = g ( x)
dx dx dx
Si g(x) = 0 la ecuación es homogénea.

dny d n−1 y dy
an ( x) n + an−1 ( x) n−1 +  + a1 ( x) + a0 ( x) y = 0
dx dx dx
Veremos que para resolver una ecuación no homogénea tendremos que resolver
también la ecuación homogénea asociada.
4
Operadores diferenciales
• Sea Dy = dy/dx. Al símbolo D se le llama operador diferencial. Definimos
a un operador diferencial de n-ésimo orden u operador polinominal
como
L = an ( x) D n + an−1 ( x) D n−1 +  + a1 ( x) D + a0 ( x)

• El operador diferencial L es un operador lineal:

L{f ( x) + g ( x)} = L( f ( x)) + L( g ( x))

Podemos escribir las EDOs anteriores simplemente como

L(y) = 0 y L(y) = g(x)

5
Principio de superposición (ecuaciones homogéneas)
Sean y1, y2, …, yk soluciones de una ecuación diferencial homogénea de
n-ésimo orden en un intervalo I. Entonces la combinación lineal
y = c1y1(x) + c2y2(x) + …+ ckyk(x)
donde ci, i = 1, 2, …, k, son constantes arbitrarias, también es una
solución en el intervalo.
Ejemplo: Las funciones y1 = x2, y2 = x2 ln x son ambas soluciones en (0, ) de
x3 y − 2 xy + 4 y = 0
Luego y = x2 + x2 ln x también es una solución en (0, ).
Nota:
(A) y(x) = cy1(x) también es solución si y1(x) es una solución.
(B) Una ED lineal homogénea siempre posee la solución trivial y(x) = 0.
6
Dependencia e independencia lineal
Un conjunto de funciones f1(x), f2(x), …, fn(x) es linealmente dependiente en
un intervalo I, si existen ciertas constantes c1, c2, …, cn no todas nulas, tales
que:

c1f1(x) + c2f2(x) + … + cn fn(x) = 0

Si el conjunto no es linealmente dependiente, entonces es linealmente

independiente.
En otras palabras, si el conjunto es linealmente independiente, cuando:
c1f1(x) + c2f2(x) + … + cn fn(x) = 0
entonces necesariamente c1 = c2 = … = cn = 0.

7
¿Son estas funciones linealmente independientes?

c1f1(x) + c2f2(x) = 0
8
Ejemplo: Las funciones f1 = cos2 x, f2 = sin2 x, f3 = sec2 x, f4 = tan2 x son
linealmente dependientes en el intervalo (-/2, /2)
porque
c1 cos2 x +c2 sin2 x +c3 sec2 x +c4 tan2 x = 0
con
c1 = c2 = 1, c3 = -1, c4 = 1.

Ejemplo: Las funciones f1 = x½ + 5, f2 = x½ + 5x, f3 = x – 1, f4 = x2

son linealmente dependientes en el intervalo (0, ), porque

f2 = 1 f1 + 5 f3 + 0 f4

9
Wronskiano
Supongamos que cada una de las funciones f1(x), f2(x), …, fn(x) posee al menos
n – 1 derivadas. El determinante

f1 f2  fn
f1 ' f2 '  fn '
W ( f1 ,..., f n ) =
  
f 1( n−1) f 2( n−1)  f n( n−1)

se llama el Wronskiano de las funciones.

10
 TEOREMA Criterio para soluciones linealmente independientes

Sean y1(x), y2(x), …, yn(x) soluciones de una ED homogénea de n-ésimo orden

en un intervalo I. Este conjunto de soluciones es linealmente independiente
sí y sólo si W(y1, y2, …, yn)  0 para todo x en el intervalo.

DEFINICIÓN Conjunto fundamental de soluciones

Cualquier conjunto {y1(x), y2(x), …, yn(x)} de n soluciones linealmente

independientes de una ED homogénea de n-ésimo orden se llama conjunto
fundamental de soluciones.

11
 TEOREMA
Existencia de un conjunto fundamental

Existe un conjunto fundamental de soluciones para una ED lineal homogénea

de orden n en un intervalo I.

TEOREMA

Solución general (ecuaciones homogéneas)

Sea y1(x), y2(x), …, yn(x) un conjunto fundamental de soluciones de nuestra

ED lineal homogénea en un intervalo I. Entonces la solución general es
y = c1y1(x) + c2y2(x) + … + cnyn(x)
donde ci son constantes arbitrarias.
CH3_12x
12
• Las funciones y1 = e3x, y2 = e-3x son soluciones de y” – 9y = 0 en (-, )
Observa que e 3x
e −3 x
W ( e 3 x , e −3 x ) = −3 x
= −6  0
3e 3x
− 3e

• para todo x.
• Luego son independientes.
Así que y = c1y1 + c2y2 es la solución general.

Por ejemplo, la función y = 4 sinh(3x) - 5e3x es una solución.

Observemos que
−3 x −3 x −3 x  e3 x − e −3 x  −3 x
y = 2e − 7e
3x
= 2e − 2e
3x
− 5e = 4  − 5e
 2 
= 4 sinh 3x – 5e-3x
• Las funciones y1 = ex, y2 = e2x , y3 = e3x son
soluciones de y’’’ – 6y” + 11y’ – 6y = 0 en (-, ).
Como x 2x 3x
e e e
x
W (e , e 2x
,e )= e
3x x
2e 2x
3e = 2e  0
3x 6x

ex 4e2 x 9e3 x
para todo valor real de x.

y = c1ex + c2 e2x + c3e3x es la solución general en (-, ).

14
 TEOREMA
Solución General (Ecuaciones no homogéneas)

Sea yp cualquier solución particular de una EDO no homogénea en un intervalo I.

Y sea y1(x), y2(x), …, yk(x) un conjunto fundamental de soluciones de su EDO

homogénea asociada, entonces la solución general de la ecuación en el intervalo es
y= c1y1 + c2y2 +… + ckyk + yp
donde las ci , i= 1,2,….,n son constantes arbitrarias

y = c1y1 + c2y2 +… + ckyk + yp = yc + yp

Nota:
yc =función complementaria
yp =una solución particular
15
• La función yp = -(11/12) – ½ x es una solución particular de

y − 6 y + 11y − 6 y = 3 x

La solución general es
11 1
y = yc + y p = c1e + c2e x 2x
+ c3e 3x
− − x
12 2

16
 TEOREMA
Principio de superposición
(ecuaciones no homogéneas)
Dadas k EDOs
an ( x) y ( n)
+ an−1 ( x) y ( n −1)
+  + a1 ( x) y + a0 ( x) y = gi ( x)
con i = 1, 2, …, k.
Si ypi denota una solución particular de la ED i-ésima correspondiente
a gi(x), tenemos que
y p = y p1 ( x ) + y p2 ( x ) +  + y pk ( x )

es una solución particular de

an ( x) y ( n ) + an−1 ( x) y ( n−1) +  + a1 ( x ) y + a0 ( x ) y
= g1 ( x ) + g 2 ( x ) +  + g k ( x )
17
• Observemos que
yp1 = -4x2 es una solución particular de
y"−3 y '+4 y = −16 x 2 + 24 x − 8
yp2 = e2x es una solución particular de
y"−3 y '+4 y = 2e 2x

yp3 = xex es una solución particular de

y"−3 y '+4 y = 2 xe x − e x
Entonces y = y p1 + y p2 + y p3 es una solución de

y − 3 y + 4 y = − 16 x 2
+ 24
  x − 8 + 2e 2x
+ 2 xe x
− e

x

g1 ( x ) g2 ( x ) g3 ( x )

18
Reducción de orden
Sabemos que la solución general de

a2 ( x) y + a1 ( x) y + a0 ( x) y = 0
es y = c1y1 + c2y2.
Supongamos que y1(x) denota una solución conocida (no trivial). Puesto
que la solución y2 es linealmente independiente, supongamos que
y2(x) = u(x) y1(x).
Nuestro objetivo será encontrar una tal u(x).
El método se conoce como reducción de orden.

19
Ejemplo:
Dada y1 = ex solución de y” – y = 0, hallar la segunda solución y2 por el método de
reducción de orden.
Solución Si y(x) = u(x)ex, entonces

y = uex + e xu , y = uex + 2e xu + e xu

que sustituyendo en la EDO: y"− y = e x (u"+2u ' ) = 0

Como ex  0, nuestra EDO se convierte en: u"+2u' = 0

Ahora "reduciremos" el orden de la ED gracias al cambio:
w = u’ w'+2w = 0
que integrando por separación de variables y deshaciendo el cambio, nos proporciona:

w = c1e −2 x = u u = −1/2 c1e −2 x + c2

20
Hemos hallado la segunda solución y2 por el método de reducción de orden:

c1 − x
y = u ( x)e = − e + c2e x
x
2

Recordemos que teníamos y1 = ex como primera solución de y” – y = 0.

Si tomamos c2 = 0, c1 = -2 para nuestra segunda solución, tenemos y2 = e-x.

Observa que W(ex, e-x)  0 para todo x, de modo que las soluciones son
independientes.

21
Caso general
• Escribimos la EDO en la forma estándar
y + P( x) y + Q( x) y = 0
Sea y1(x) una solución conocida de la EDO e y1(x)  0 para todo x en el
intervalo.
• Si definimos y(x) = u(x)y1(x), tenemos
y = uy1 + y1u  y = uy1 + y1u + y1u + y1u  y = uy1 + 2y1u + y1u
y + Py + Qy = uy(  
1 + 2y 1u + y 1u ) (
 + P uy1 + y 1u + Quy 1 )
= u[y    
1 + Py 1 + Qy 1 ] + y 1u + (2y 1 + Py 1 )u = 0
cero
 y1u + (2 y1 + Py1 )u = 0
empleando el cambio w = u’.
dw
y1w + (2 y1 + Py1 )w = 0 y1 + (2 y1 + Py1 )w = 0
dx
Dividiendo entre y1w y multiplicando por dx:
dw y1 dw y1
w
+ 2 dx = −Pdx 
y1
 w + 2 y1 dx = − Pdx + c
−  Pdx
ln|w| + ln|y |= −  Pdx + c  ln|wy |= −  Pdx + c  wy = c1e
2 2 2
1 1 1
Luego 1 −  Pdx e−  Pdx
 w = u´= c1 2 e  u = c1  2 dx + c2
y1 y1
Tomando c1 = 1, c2 = 0, obtenemos
e− P ( x ) dx
y2 = y1 ( x)  2 dx
y1 ( x)
23
La función y1= x2 es una solución de
x y"−3 xy '+4 y = 0
2

Hallar la solución general en (0, ).

Solución:
La forma estándar es 3 4
y − y + 2 = 0
x x
Dando los pasos anteriores, demuestra que:
3  dx / x
e
y2 = x  dx = x 2 ln x
2

La solución general es: x4

y = c1 x + c2 x ln x
2 2

24
Coeficientes indeterminados
Si queremos resolver
( n−1)
an y (n)
+ an−1 y +  + a1 y + a0 y = g ( x)
Tenemos que hallar y = yc + yp. Veamos cómo hacerlo,
en este caso general, mediante el método conocido
como de coeficientes indeterminados.

25
 Coeficientes indeterminados
Simplemente haremos una conjetura sobre la forma de la posible
solución particular a partir de la g(x) que deberá ser un polinomio, seno o
coseno, exponencial o combinación lineal de todas ellas...

Gracias a que las derivadas de las combinaciones lineales de estas

funciones vuelven a ser ellas mismas, parece razonable que busquemos
soluciones particulares de la misma forma...

Vamos a ilustrar la idea con algunos ejemplos

26
Resolver
y"+4 y '−2 y = 2 x − 3 x + 6
2
Solución:
Ya sabemos cómo obtener una solución yc de la
ecuación homogénea asociada. Ahora, queremos
hallar yp.
Como el lado derecho de la ED es un polinomio,
supondremos y p = Ax2 + Bx + C , entonces,
tras sustituir: y p ' = 2 Ax + B, y p "= 2 A

2A + 8Ax + 4B – 2Ax2 – 2Bx – 2C = 2x2 – 3x + 6

− 2 A = 2 , 8 A − 2B = −3 , 2 A + 4B − 2C = 6
A = −1, B = −5/2, C = −9 5
yp = −x − x − 9
2
2
27
Hallar una solución particular de
y"− y'+ y = 2sen (3x)
Solución:
Probemos yp = A cos(3x) + B sen(3x)
Tras sustituir,
(−8 A − 3B) cos(3x) + (3 A − 8B) sin( 3x) = 2 sin( 3x)
Luego
A = 6/73, B = −16/73
6 16
y p = cos(3x) − sen(3x)
73 73

28
Resolver y"−2 y '−3 y = 4 x − 5 + 6 xe 2 x
Solución:
yc = c1e− x + c2e3 x Solución homogénea
Pensando en el principio de superposición:

Probemos y p = Ax + B + Cxe + Ee
2x 2x

Tras sustituir,
−3 Ax − 2 A − 3B − 3Cxe 2 x + (2C − 3E )e 2 x = 4 x − 5 + 6 xe 2 x

Luego , A = −4/3, B = 23/9, C = −2, E = −4/3

4 23 4 2x
y p = − x + − 2 xe − e
2x
3 9 3
−x 4 23  4  2x
y = c1e + c2e − x + −  2 x +  e
3x
3 9  3
29
Determinar una yp de
y"−5 y '+4 y = 8e x

Solución:
Probemos: yp = Aex
Tras sustituir: 0 = 8ex (conjetura incorrecta)
El problema está en que la función complementaria es:
Y la suposición ya está presente en yc. yc = c1e + c2e
x 4x

Probemos como alternativa: yp = Axex.

Tras sustituir: -3Aex = 8ex
Entonces: A = -8/3,
yp = (−8/3)xe2x

30
• Si ninguna función en la supuesta yp es parte de yc
En la siguiente tabla se muestran soluciones particulares de prueba.
g (x) Forma de yp
1. 1(una constante) A
2. 5x + 7 Ax + B
3. 3x 2 − 2 Ax2 + Bx + C
4. x3 − x + 1 Ax3 + Bx 2 + Cx + E
5. sen 4 x A cos 4 x + B sen 4 x
6. cos 4 x A cos 4 x + B sen 4 x
7. e5 x Ae5 x
8. (9 x − 2)e5 x ( Ax + B )e5 x
9. x 2 e5 x ( Ax 2 + Bx + C )e5 x
10. e3 x sen 4 x Ae 3 x cos 4 x + Be3 x sen 4 x
11. 5 x 2 sen 4 x ( Ax 2 + Bx + C ) cos 4 x + ( Ex2 + Fx + G) sen 4 x
12. xe3 x cos 4 x ( Ax + B)e3 x cos 4 x + (Cx + E )e3 x sen 4 x

31
Hallar la forma de yp de
(a) y"−8 y '+25 y = 5 x3e − x − 7e − x
Solución:
Tenemos que g ( x) = (5 x3 − 7)e − x y probamos con
−x
y p = ( Ax + Bx + Cx + E )e
3 2

No hay duplicación entre los términos yp e yc

(b) y” + 4y = x cos x
Solución:
Probamos con x p = ( Ax + B) cos x + (Cx + E ) sin x
Tampoco hay duplicidad entre los términos yp y yc .

32
Hallar la forma de yp de
y − 9 y + 14 y = 3x 2 − 5sen 2 x + 7 xe 6 x
Solución:
Para 3x2:
y p1 = Ax 2 + Bx + C
Para -5 sen 2x:
y p2 = E cos 2 x + Fsen2 x
Para 7xe6x:
y p3 = (Gx + H )e6 x
Ningún término de y p = y p1 + y p2 + y p3 duplica un
término de yc

33
Así que la regla formal en este caso es que la solución
particular es una combinación lineal de las funciones
linealmente independientes que se generan mediante
diferenciaciones repetidas de g(x).

¿Y cuál es la regla si la solución particular así propuesta

es también una solución de la ecuación homogénea
asociada?

Si alguna yp contiene términos que duplican los términos de

yc, entonces esa yp se debe multiplicar por xn, donde n es el
entero positivo más pequeño que elimina esa duplicación.

34
Resolver y"+ y = 4 x + 10 senx, y( ) = 0, y' ( ) = 2
Solución: yc = c1 cos x + c2 senx
Primero probamos: yp = Ax + B + C cos x + E sen x
Pero hay una duplicación.
Entonces probamos con
yp = Ax + B + Cx cos x + Ex sen x

Tras sustituir y simplificar,

A = 4, B = 0, C = -5, E = 0

Luego y = c1 cos x + c2 sen x + 4x – 5x cos x

Como y() = 0, y’() = 2, tenemos

y = 9 cos x + 7 sen x + 4x – 5x cos x

35
Resolver y"−6 y '+9 y = 6 x 2 + 2 − 12e3 x
Solución:
yc = c1e3x + c2xe3x Este término está
duplicado, aparece ya
yp = 
Ax
+
2
+
Bx C + Ee

3x en yc.

y p1 y p2

Debemos probar con: yp = 

Ax
+
Bx+
2
C+
Exe 2 3x

Tras sustituir y simplificar, y p1 y p2

A = 2/3, B = 8/9, C = 2/3, E = -6

Luego
2 2 8 2
y = c1e + c2 xe + x + x + − 6 x e
3x 3x 2 3x
3 9 3

36
Resolver y + y"= e cos x
x

Solución:
m3 + m2 = 0; m = 0, 0, -1
yc = c1+ c2x + c3e-x
Probamos como solución particular:
yp = Aex cos x + Bex sen x
Tras sustituir y simplificar,
A = -1/10, B = 1/5
Luego,
−x 1 x 1 x
y = yc + y p = c1 + c2 x + c3e − e cos x + e senx
10 5

37
Hallar la forma de yp de y ( 4 ) + y = 1 − x 2e − x
Solución:
yc = c1+ c2x + c3x 2 + c4e-x
2 −x −x −x
Prueba: y p = A + 
Bxe +Cxe
 + Ee
y p1 y p2
Como aparece repetido en la solución homogénea, necesitaremos
multiplicar A por x3 y (Bx2e-x + Cxe-x + Ee-x) por x. Prueba ahora:

3 −x 2 −x −x
yp = 
Ax 3 + 
Bxe +
Cxe +
Exe


y p1 y p2

38
Método de variación de parámetros
a2 ( x) y + a1 ( x) y + a0 ( x) y = g ( x)
y + P( x) y + Q( x) y = f ( x)
donde P(x), Q(x) y f(x) son continuas en I.
Conocidas y1(x) e y2(x) soluciones l. i. de
la ec. homogénea asociada, probaremos
como solución particular:
y p = u1 ( x) y1 ( x) + u2 ( x) y2 ( x)

39
 y p = u1 ( x) y1 ( x) + u2 ( x) y2 ( x)
Sustituimos yp’, yp” en la EDO:

yp + P( x) yp + Q( x) y p 0 0

= u1[ y1 + Py1 + Qy1 ] + u2[ y2 + Py2 + Qy2 ]

+ y1u1 + u1 y1 + y2u2 + u2 y2 + P[ y1u1 + y2u2 ] + 2 y1u1 + 2 y2 u2

d d
= [ y1u1 ] + [ y2u2 ] + P[ y1u1 + y2u2 ] + y1u1 + y2 u2
dx dx
d
= [ y1u1 + y2u2 ] + P[ y1u1 + y2u2 ] + y1u1 + y2 u2 = f ( x)
dx
40
 d
[ y1u1 + y2u2 ] + P[ y1u1 + y2u2 ] + y1u1 + y2 u2 = f ( x)
dx

Necesitaremos dos ecuaciones para encontrar valores de u1 y u1. Exijamos

que: y1u1’ + y2u2’ = 0, para obtener una ecuación adicional y de paso que
la EDO se reduzca a: y1’u1’ + y2’u2’ = f(x).
De modo que nos queda el sistema de ecuaciones:

y1u1’ + y2u2’ = 0

y1’u1’ + y2’u2’ = f(x)

41
Expresado en términos de determinantes

W1 y2 f ( x ) y W2 y1 f ( x )
u1 = =− u 2 = =
W W W W
donde

y1 y2 0 y2 y1 0
W= , W1 = , W2 =
y1 y2 f ( x) y2 y1 f ( x)
De donde encontraremos, por integración, las soluciones.

42
 y"−4 y '+4 y = ( x + 1)e 2x

Solución:
m2 – 4m + 4 = 0, m = 2 (cero doble)
y1 = e2x, y2 = xe2x,
Resolver e 2x
xe2 x
W (e 2 x , xe2 x ) = = e 4x
0
2e 2 x 2 xe + e
2x 2x

Como f(x) = (x + 1)e2x, entonces:

0 xe2 x e2 x 0
W1 = = −( x + 1) xe , W2 = 2 x
4x
= ( x + 1) e 4x

( x + 1)e2 x 2 xe 2x
2e ( x + 1)e2 x

W1 y2 f ( x ) W2 y1 f ( x)
u1 = =− u2 = =−
W W W W
( x + 1) xe 4x
( x + 1)e 4x
u1 = − 4x
= − x − x , u2 = −
2
4x
= x +1
e e 43
 u1 = − x − x ,
2
u2 = x + 1
Luego
u1 = (-1/3)x3 – ½ x2, u2 = ½ x2 + x

Recordemos que: y p = u1 ( x) y1 ( x) + u2 ( x) y2 ( x)
y1 = e2x, y2 = xe2x

 1 3 1 2  2x  1 2  1 3 2x 1 2 2x
x p =  − x − x e +  x + x  xe = x e + x e
2x
 3 2  2  6 2
1 3 2x 1 2 2x
y = yc + y p = c1e + c2 xe + x e + x e
2x 2x
6 2 44
 Resolver 4 y"+36 y = csc 3x
Solución:
y” + 9y = (1/4) csc 3x
m2 + 9 = 0, m = 3i, -3i
y1 = cos 3x, y2 = sin 3x, f(x) = (1/4) csc(3x)
Como

cos 3x sin 3x
W (cos 3x , sin 3x) = =3
− 3sin 3x 3 cos 3x
0 sin 3x 1 cos 3x 0 1 cos 3x
W1 = = − , W2 = =
1/4 csc 3x 3 cos 3x 4 − 3 sin 3x 1/4 csc 3x 4 sin 3x

45
 W1 1 W2 1 cos3x
u1 = =− u2 = =
W 12 W 12 sen3x

Entonces u1 = −1/12 x, u2 = −1/36 ln | sen3x |

1 1
y p = − x cos3x + (sen3x) ln | sen3x |
12 36
1 1
y = yc + y p = c1 cos3x + c2 sen3x − x cos3x + (sen3x) ln | sen3x |
12 36

46
Resolver 1
y"− y =
Solución: x
m2 – 1 = 0, m = 1, -1
y1 = ex, y2 = e-x, f(x) = 1/x, y W(ex, e-x) = -2
Luego
e− x (1/ x) 1 x e−t
u1 = − , u1 =  dt
−2 2 x0 t
x t
e (1/ x) 1 ex

u2 = − , u2 = −  dt
−2 2 x0 t
1 x x e−t 1 − x x et
yp = e  dt − e  dt
2 x0 t 2 x0 t

−t t
1 x e 1 x e
y = yc + y p = c1e x + e x  dt − e− x  dt
2 x0 t 2 x0 t
47
Ecuación de Cauchy-Euler
Forma de ecuación de Cauchy-Euler
n n −1
d y n −1 d y dy
an x n
n
+ an−1x n −1
+  + a1x + a0 y = g ( x)
dx dx dx
• Método de solución
Probamos y(x) = xm, donde debemos determinar m, para
resolver la ecuación homogénea asociada: Observa que:
k
d y m−k
ak xk
k
= ak x k
m ( m − 1)( m − 2)  ( m − k + 1) x
dx
= ak m(m − 1)(m − 2)(m − k + 1) x m

(an m(m −1)(m − 2)(m − n + 1) + ... + a1m + a0 )xm = 0

48
 2
Ecuación auxiliar ax 2 d y
2
dy
+ bx + cy = g ( x)
dx dx
Para n = 2, y = xm, tenemos Observa que tenemos
(am(m – 1) + bm + c)xm = 0, o que ax2 es igual a
am2 + (b – a)m + c = 0 cero en x = 0. Para
Caso 1: Raíces reales y distintas asegurar existencia y
unicidad, tomaremos
y = c1 x + c2 x
m1 m2
I = (0, ).

2
d
Resolver x 2 y − 2 x dy − 4 y = 0
2
Solución: dx dx
Tenemos a = 1, b = -2 , c = -4
m2 – 3m – 4 = 0, m = -1, 4,
y = c1x-1 + c2x4
49
Caso 2: Raíces reales repetidas
• Dedujimos y2 = xm 1
ln x
Luego
y = c1 x + c2 x ln x
m1 m1

Resolver d 2
y dy
4x 2
2
+ 8x + y = 0
Solución: dx dx
(am(m – 1) + bm + c)xm = 0, o am2 + (b – a)m + c = 0
Tenemos a = 4, b = 8, c = 1
4m2 + 4m + 1 = 0, m = -½ , -½
−1/ 2 −1/ 2
y = c1 x + c2 x ln x
50
Caso 3: Raíces complejas conjugadas
• Orden superior: multiplicidad k
k −1
x , x ln x , x (ln x) ,  , x (ln x)
m1 m1 m1 2 m1

• Caso 3: raíces complejas conjugadas

m1 =  + i, m2 =  – i,
y = C1x( + i) + C2x( - i)
Como
xi = (eln x)i = ei ln x = cos( ln x) + i sen( ln x)
x-i = cos ( ln x) – i sen ( ln x)
Luego
y = c1x cos( ln x) + c2x sen( ln x)
= x [c1 cos( ln x) + c2 sen( ln x)]

51
 1
Resolver 4 x y + 17 y = 0, y (1) = −1, y ' (1) = −
2

Solución: (am(m – 1) + bm + c)xm = 0, o am2 + (b – a)m + c2= 0

Tenemos a = 4, b = 0 , c = 17
4m2 − 4m + 17 = 0, m = ½ + 2i
y = x1/ 2 [c1 cos( 2 ln x) + c2 sin( 2 ln x)]
Aplicando y(1) = -1, y’(1) = 0, tenemos que c1 = -1, c2 = 0,

y = − x1/2 cos( 2 ln x)

52
 3 2
Resolver d y d y dy
x 3
3
+ 5x 2
2
+ 7x + 8y = 0
dx dx dx
Solución:
Sea y = xm,
2
dy d y
m−1
= mx , 2 = m(m − 1) x m−2 ,
dx dx
d3y m−3
= m ( m − 1)( m − 2) x
dx3

Luego temenos: xm (m + 2)(m2 + 4) = 0

m = -2, m = 2i, m = -2i
y = c1x -2 + c2 cos(2 ln x) + c3 sin(2 ln x)
53
Resolver x 2 y"−3xy '+3 y = 2 x 4e x
Solución:
Tenemos (m – 1)(m – 3) = 0, m = 1, 3
yc = c1x + c2x3
Usando variación de parámetros,

yp = u1y1 + u2y2, donde y1 = x, y2 = x3

3 3
 
Escribimos la ED como y − y + 2 y = 2 x e
2 x
x x
Luego P(x) = -3/x, Q(x) = 3/x2, f(x) = 2x2ex

54
Así x x3
W= 2
= 2 x 3
,
1 3x
0 x3 x 0
W1 = = −2 x 5 x
e , W2 = 2 x = 2 x 3 x
e
2 x 2e x 3x 2
1 2x e

Hallamos
2 x 5e x 2 x 5 x
e
u1 = − 3
= − x 2 x
e , u 
2 = = e x
2x 2 x3
u1 = − x 2e x + 2 xe x − 2e x , u2 = e x

y p = u1 y1 + u2 y2 = (− x 2e x + 2 xe x − 2e x ) x + e x x3
= 2 x 2e x − 2 xe x
y = yc + y p = c1x + c2 x + 2 x e − 2 xe
3 2 x x
55
Una ecuación de Cauchy-Euler siempre se puede escribir como un lineal
de coeficientes constantes haciendo el cambio de variable: x = et.
Por ejemplo: Resuelve así: x 2
y − xy + y = ln x
dt 1
x = e  t = ln x 
t
=
dx x
dy dy dt 1 dy
= =
dx dt dx x dt
d 2 y d  dy  d  1 dy  1 dy 1  d  dy  
=  =  =− 2 +   
dx 2
dx  dx  dx  x dt  x dt x  dx  dt  
1 dy 1  d  dy   1 dy 1  d  1 dy  
=− 2 +    = − 2 +   
x dt x  dt  dx   x dt x  dt  x dt  
1  d 2 y dy 
= 2 2 − 
x  dt dt 
56
x 2 y − xy + y = ln x
dy 1 dy d 2 y 1  d 2 y dy 
=  2 = 2 2 − 
dx x dt dx x  dt dt 
d 2 y dy dy
2
− − + y=t
dt dt dt

d2y dy
2
−2 + y =t
dt dt
y = c1et + c2te t + 2 + t

y = c1 x + c2 x ln x + 2 + ln x
57
BIBLIOGRAFIA
https://www.google.com.pe/search?q=(%C2%A9+Chema+Madoz,+VEGAP,+Madrid+2009)&tbm=isch&source=iu&ictx=1
&fir=SNleSXNRch0W3M%252C5G9ylbhrbgBBkM%252C_&vet=1&usg=AI4_-
kQtXl22a749gwvMV_E0zwlOGYpnqQ&sa=X&ved=2ahUKEwjq-
K_f1r3qAhVuILkGHULtARYQ9QEwAnoECAkQGQ&biw=1020&bih=613#imgrc=-Cti5MZ_LU5w_M
https://www.google.com.pe/search?q=(%C2%A9+Chema+Madoz,+VEGAP,+Madrid+2009)&tbm=isch&source=iu&ictx=1
&fir=SNleSXNRch0W3M%252C5G9ylbhrbgBBkM%252C_&vet=1&usg=AI4_-
kQtXl22a749gwvMV_E0zwlOGYpnqQ&sa=X&ved=2ahUKEwjq-
K_f1r3qAhVuILkGHULtARYQ9QEwAnoECAkQGQ&biw=1020&bih=613#imgrc=Qkbzjq79FkYIcM