3.

Ecuaciones diferenciales de orden

superior

( Chema Madoz, VEGAP, Madrid 2009)

n
n

1
d
y
d
d
y
any(x)0
(y0,x0)

a
(y1,y(n1x))0

a
(nx1)y
a
g(x)
n

1
1
0
Ecuaciones lineales: teora
bsica

Un problema de valor inicial de n-simo

orden consiste en resolver la EDO lineal:

sujeta a las n condiciones iniciales:

Resolverlo consiste en encontrar una funcin

y(x) definida en un intervalo I que contiene a
x0, donde se cumplen la ecuacin y las
condiciones iniciales.

Existencia de una solucin

nica
(Condicin suficiente)

Sea an(x), an-1(x), , a0(x), y g(x) continuas

en I,
con an(x) 0 para todo x de I. Si x = x0 es
cualquier punto de este intervalo, entonces
existe
una solucin
Ejemplo:
y del
3 y 5 yy(x)
7 y problema
0, y (1) 0 , yanterior
(1) 0, y(1)en
0 I
y es
nica.

posee la solucin trivial y(x) = 0. Como es una ED de

tercer orden lineal con coeficientes constantes, y(x) =
0 es la nica solucin en cualquier intervalo que
3
contenga a x = 1.

 Ejemplo: Comprueba que y = 3e2x + e2x 3x

es la nica solucin de
y"4 y 12 x, y (0) 4, y ' (0) 1
La ED es lineal, los coeficientes y g(x) son
todos funciones continuas, y a2(x) = 1 es
distinto de 0 en cualquier intervalo que
contenga x = 0. La solucin propuesta
cumple la EDO y es nica en I.
Comprueba que y = cx2 + x + 3 es solucin del PVI:

x 2 y 2 y 2 y 6,

y (0) 3 , y (0) 1

en toda la recta real. Este PVI tiene infinitas soluciones. Observa que el
coeficiente de la derivada a2(x) = x2 ms alta se hace cero en x = 0 y ese
punto necesariamente tiene que estar incluido en I porque lo imponen las
condiciones iniciales.
4

2
d
y
d
y
a2y(x)0
a,1(bx)
(1x)y
a
g(x)
0

Problemas de valores en la
frontera
Resolver:
sujeta a :

se llama problema de valor

en la frontera (PVF) y a las
restricciones se conocen
como condiciones de contorno
o condiciones en la frontera.
Nota: Las condiciones de contorno

pueden ser tambin sobre las derivadas.

5

Vimos que x = c1 cos 4t + c2 sin 4t era solucin de

x"16 x 0

x 16 x 0 , x(0) 0 , x 0
(a) Supongamos el PVF
2
Si x(0) = 0, entonces c1 = 0, y x(t) = c2 sen 4t.
Si x(/2) = 0, obtenemos 0 = 0 independientemente
de c2. De modo que tenemos infinitas soluciones.

x 16 x 0 , x(0) 0 , x 0
(b) Si
8
tenemos que c1 = 0, c2 = 0:
x(t) = 0, solucin nica.

(c) Si x 16 x 0 , x(0) 0 , x 1
2

tenemos que c = 0, y 1 = 0
1

(contradiccin). No hay solucin.

n
n

1
d
y
d
y
d
y
aann((xx))dny
(
a
)an1(x)dn1y
x
(

)
a
x

a
(
x
)
y

g
(
x
)
n

1
1
0
d
y
a1(x)
a0(x)y
0
La siguiente EDO lineal de orden n:

se dice que es no homognea.

si g(x) = 0 la ecuacin es homognea.

Veremos que para resolver una ecuacin no homognea

tendremos que resolver tambin la ecuacin homognea
asociada.

Operadores diferenciales
Sea Dy = dy/dx. Al smbolo D se le llama operador
diferencial. Definimos a un operador diferencial de
n-simo orden u operador polinominal como

L an ( x) D n an1 ( x) D n1 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)
8

Principio de superposicin
(ecuaciones homogneas)

Sean y1, y2, , yk soluciones de una ecuacin

diferencial homognea de n-simo orden en
un intervalo I. Entonces la combinacin lineal
y = c1y1(x) + c2y2(x) + + ckyk(x)
donde ci, i = 1, 2, , k, son constantes
2
2
Ejemplo:
Las
funciones
y
=
x
,
y
=
x
ln x
1
2
arbitrarias, tambin 3es una solucin
enson
el ambas
x y 2 xy 4 y 0
soluciones
en
(0,
)
de
intervalo. 2
2
Luego y = x + x ln x tambin es una solucin en (0,
).

Nota:
(A) y(x) = cy1(x) tambin es solucin si y1(x) es una solucin.

9
(B) Una ED lineal homognea siempre posee la solucin trivial y(x)
=

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
10
entonces necesariamente c = c = = c = 0.

Son estas funciones linealmente independientes?

c1f1(x) + c2f2(x) = 0
11

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
12

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

W ( f1 ,..., f n )

f1

f2

fn

f1 '

f2 '

fn '

f 1( n1)

f 2( n1)

f n( n1)

se llama el Wronskiano de las funciones.

13

TEOREMA

Criterio para soluciones

linealmente independientes

Sean y1(x), y2(x), , yn(x) soluciones de una

ED homognea de n-simo orden en un intervalo

Este conjunto de soluciones es linealmente
independiente si y slo si W(y1, y2, , yn) 0 para
todo x en el intervalo.
DEFINICIN

Conjunto fundamental de soluciones

Cualquier conjunto y1(x), y2(x), , yn(x) de n

soluciones linealmente independientes de una ED

homognea de n-simo orden se llama conjunto
14
fundamental de soluciones.

TEOREMA

Existencia de un conjunto fundamental

Existe un conjunto fundamental de solucion

para una ED lineal homognea de orden n en un
intervalo I.
TEOREMA

Solucin general (ecuaciones homogneas)

Sea y1(x), y2(x), , yn(x) un conjunto fundament

de soluciones de nuestra ED lineal homognea e

un intervalo I. Entonces la solucin general es
y = c1y1(x) + c2y2(x) + + cnyn(x)
CH3_15x
donde ci son constantes
arbitrarias.
15

 Las funciones y1 = e3x, y2 = e-3x son

soluciones de
y 9y = 0 en (-, )
3x
3 x
e
e
Observa que
3x
3 x
W (e , e )
6 0
3x
3 x
3e

3e

para todo x. Luego son independientes.

As que y = c1y1 + c2y2 es la solucin
general.
Por ejemplo, la funcin y = 4 sinh(3x) 5e3x
es una solucin.
que
3 x Observemos
3 x

y 2e 2e
3x

3 x

5e

3 x

e e
4
2

5e 3 x

16

= 4 sinh 3x

ex

e2 x

e3 x

x, y = e2x , y = e3x son

Las funciones
y
=
e
1
x
2x
3x
x 2 2x
3 x3
6x
W
(
e
,
e
,
e
)

e
2
e
3
e

2
e
0
soluciones de y 6y + 11y 6y = 0 en(-,
).
e x 4e 2 x 9e3 x
Como

para todo valor real de x.

y = c1ex + c2 e2x + c3e3x es la solucin general
en (-, ).
17

TEOREMA

Solucin General

(Ecuaciones no homogneas)
Sea yp cualquier solucin particular de una EDO no
homognea en un intervalo I. Y sea y1(x), y2(x), , yk(x)
un conjunto fundamental de soluciones de su EDO
homognea asociada, entonces la solucin general de
la ecuacin 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

= funcin complementaria + una solucin particular
18

 La funcin yp = -(11/12) x es una

solucin particular de

y 6 y 11 y 6 y 3 x
La solucin general es

11 1
y yc y p c1e c2e c3e x
12 2
x

2x

3x

19

TEOREMA

Dadas k EDOs

Principio de superposicin
(ecuaciones no homogneas)

an ( x) y ( n ) an 1 ( x) y ( n 1) a1 ( x) y a0 ( x) y g i ( x)

con i = 1, 2, , k.
Si ypi denota una solucin 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 solucin particular de
an ( x) y ( n ) an1 ( x) y ( n1) a1 ( x) y a0 ( x) y

g1 ( x) g 2 ( x) g k ( x)
20

y"3 y '4 y 16 x 2 24 x 8
Observemos que
yp1 = -4x2 es una solucin particular de

y"3 y '4 y 2e 2 x

yp2 = e2x es una solucin particular de

y"3 y '4 y 2 xe x e x

yp3 = xex es una solucin particular de

y y p1 y p2 y p3
Entonces

es una solucin de

2
2x
x
x
y 3 y 4 y
16
x

24
x

2
e

2
xe

g1 ( x )

g2 ( x )

g3 ( x )
21

Reduccin de orden
Sabemos que la solucin general de

a2 ( x) y a1 ( x) y a0 ( x) y 0

es y = c1y1 + c2y1.
Supongamos que y1(x) denota una solucin
conocida (no trivial). Puesto que la solucin y2
es linealmente independiente, supongamos
que y2(x) = u(x) y1(x). Nuestro objetivo ser
encontrar una tal u(x). El mtodo se conoce
como reduccin de orden.
22

Dada y1 = ex solucin de y y = 0, hallar la segunda solucin y2 por el

mtodo de reduccin de orden.
Solucin
Si y(x) = u(x)ex, entonces

que sustituyendo en la EDO:x

x
x
x
x
Como ex 0, nuestra EDO se convierte en:
Ahora "reduciremos" el orden de la ED gracias al cambio:
x
y
"

e
(u"2u ' ) 0
w = u
que integrando por separacin
u"2u '
de variables y deshaciendo el
cambio, nos proporciona:

y ue e u , y ue 2e u e u

w'2 w 0

w c1e 2 x u
u 1/2 c1e 2 x c2

23

Hemos hallado la segunda solucin y2

por el mtodo de reduccin de orden:
c1 x
y u ( x)e e c2e x
2
x

Recordemos que tenamos y1 = ex

como primera solucin de y y = 0. Si
tomamos c2 = 0, c1 = -2 para nuestra
segunda solucin, tenemos y2 = e-x.
Observa que W(ex, e-x) 0 para todo x,
de modo que las soluciones son
independientes.
24

yP
(x)yQ
(x)y0

Caso general
Escribimos la EDO en la forma
estndar

Sea y1(x) una solucin conocida de la

EDO e y1(x) 0 para todo x en el
intervalo.

u
y

y
u
,
y

u
y

2
y
u

y
u
1
1
1
1
1
Si definimos y(x) = u(x)y1(x), tenemos

y Py Qy
u[ y1 Py1 Qy1 ] y1u (2 y1 Py1 )u 0

cero

25

y1u ( 2 y1 Py1 )u 0
dw
y1
(2 y1 Py1 ) w 0
dx
Dividiendo
dw
y1
entre y w
2 dx Pdx
y1
y multiplicando w
dw
y1
por dx:
2 dx Pdx c

empleando el cambio w =
w u. y1

y1w (2 y1 Py1 ) w 0
1

ln | wy12 | Pdx c

Luego

u c1

wy12 c1e Pdx

Pdx

y12

dx c2

Tomando c1 = 1, c2 = 0, obtenemos
y2 y1 ( x)

P ( x ) dx

y12 ( x )

dx
26

La funcin y1= x2 es una solucin de

2

x y"3 xy '4 y 0
Hallar la solucin general en (0, ).
Solucin:
La forma estndar es

3
4
y y 2 0
x
x
Dando los pasos anteriores, demuestra
que:
3 dx / x

La solucin
general 2es:
2 e

y2 x

dx x ln x
2

y c1x c2 x ln x
27

 La ecuacin diferencial ay + by = 0 se resuelve

ya sea mediante separacin de variables o
mediante la ayuda de un factor integrante.
Observa que si despejamos y de la ecuacin
diferencial ay + by = 0 se obtiene y = ky,
donde k es una constante.

Esto nos revela la "naturaleza" de la solucin: la

nica funcin elemental no trivial cuya derivada
es una mltiplo de si misma es la funcin
exponencial, y(x) = emx. Lo que resta ser
determinar el valor de m...

28

ay
bcy
0am
2
bc0

Ecuaciones lineales
homogneas con coeficientes
constantes

an y

(n)

an1 y

( n 1)

a2 y a1 y a0 y 0

donde ai son constantes, an 0.

Ecuacin o polinomio auxiliar :

Para n = 2,
Si probamos y(x) = emx,
e mx (am 2 bm c) 0

obtenemos la ecuacin auxiliar.

29

2
am

bm
c
0

Las dos races del polinomio auxiliar son:

2

m1 (b b 4ac ) / 2a
2

m2 (b b 4ac ) / 2a

(1) b2 4ac > 0: reales y distintas, m1 m2 .

(2) b2 4ac = 0: reales e iguales, m1 = m2 = -b/
(2a).
(3) b2 4ac < 0: complejas conjugadas,

m1 i , m2 i

30

 Caso 1: Races reales y distintas

La solucin general es

y c1e

m1 x

c2 e

m2 x

Por qu?

Caso 2: Races reales repetidas

y1 e

m1x

Para obtener la segunda solucin utilizamos el

mtodo de reduccin de orden, recordando que

m1 = m2 = -b/(2a).

2 m1x
e
La solucinm1general
es m1x
x
y2 e 2 m1x dx e dx xem1x
e

y c1e

m1x

c2 xe

m1x
31

ie

cos
in

Caso 3: Races complejas conjugadas

Escribimos m1 i,, una
m2 solucin
i
( i ) x
( i ) x
general es
y C1e
C2 e
Usando la frmula de Euler:

i x

i x

cos x i sin x e
e

i x

2 cos x e

i x

i x

cos x i sin x

i x

2i sin x

32

y
ceos(
x)
cesin(
x)

x
112

( i ) x
( i ) x
y

C
e

C
e
Como
es solucin
general, tomando
1
2

C1 = C2 = 1 y C1 = 1, C2 = -1 , tenemos dos
soluciones:

y1 ex (eix e ix ) 2ex cos x

y2 ex (eix e ix ) 2iex sin x

As, ex cos x y ex sen x son un conjunto

fundamental de soluciones y la solucin general
es

33

 Resolver las EDs siguientes:

(a) 2 y"5 y '3 y 0

2m 2 5m 3 (2m 1)(m 3) , m1 1/2 , m2 3

y c1e x / 2 c2 e 3 x

(b) y"10 y '25 y 0

m 2 10m 25 (m 5) 2 , m1 m2 5
y c1e5 x c2 xe5 x

(c) y"4 y '7 y 0

m 2 4m 7 0 , m1 2 3i , m2 2 3i
2 , 3 , y e 2 x (c1 cos 3 x c2 sen 3 x)

34

Resolver 4 y"4 y '17 y 0, y (0) 1, y ' (0) 2

Solucin:
4m 2 4m 17 0, m1 1/2 2i
y e x / 2 (c1 cos 2 x c2 sin 2 x)
y (0) 1, c1 1,

e y ' (0) 2,

c2 3/4

35

c
o
s
k
x

c
s
i
n
k
x
1e
2
k
x

k
x
e
yc1osh(kx)c2sinh(kx)

Resolver las ecuaciones:

y k 2 y 0,

y k 2 y 0, k 0

Para la primera ecuacin :

Para la segunda ecuacin :

Como
Luego

y1 1/2(e kx e kx ) cosh(kx)
y2 1/2(e kx e kx ) sinh( kx)

36

(a
)y
n
(
n

1
)

a
y

a
y

a
y
y

0
n

1
2
1
0
n
n

1
2
m
m
a0
n
1m
1m
Ecuaciones de orden
superior

Dada la EDO:

La ecuacin asociada

se llama su ecuacin auxiliar .

37

Resolver

y 3 y 4 y 0

Solucin:

m3 3m 2 4 (m 1)(m 2 4m 4) (m 1)(m 2) 2
m2 m3 2

y c1e x c2e 2 x c3 xe2 x

Resolver
Solucin:

d4y
d2y
2 2 y 0
4
dx
dx

m 4 2m 2 1 (m 2 1) 2 0
m1 m3 i, m2 m4 i
y C1eix C2e ix C3 xeix C4 xeix
c1 cos x c2 sin x c3 x cos x c4 x sin x38

Races complejas repetidas

Si m1 = + i es una raz compleja
de multiplicidad k, entonces m2 =
i es tambin una raz compleja
de multiplicidad k. Las 2k soluciones
linealmente independientes son :
ex cos x , xex cos x , x 2ex cos x , , x k 1ex cos x

ex sen x , xex senx , x 2 ex sen x , , x k 1ex sen x

39

Coeficientes
indeterminados

Si queremos resolver

an y

(n)

an1 y

( n 1)

a1 y a0 y g ( x)

Tenemos que hallar y = yc + yp. Veamos

cmo hacerlo, en este caso general,
mediante el mtodo conocido como de
coeficientes indeterminados.

40

Coeficientes
indeterminados
Simplemente haremos una conjetura sobre la
forma de la posible solucin particular a partir de la
g(x) que deber ser un polinomio, seno o coseno,
exponencial o combinacin 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
41

y"4 y '2 y 2 x 3 x 6
Resolver
Solucin:
Ya sabemos cmo obtener una solucin yc de la
ecuacin homognea asociada. Ahora, queremos
hallar yp.
Como el lado derecho de la ED es un polinomio,
supondremos
y p Ax 2 Bx C ,
entonces,
y p ' 2 Ax B, y p " 2 A
tras sustituir:
2A + 8Ax + 4B 2Ax2 2Bx 2C = 2x2 3x + 6

2 A 2 , 8 A 2 B 3 , 2 A 4 B 2C 6
A 1, B 5/2, C 9

5
yp x x 9
42
2
2

Hallar una solucin particular de

y" y ' y 2 sen(3x)

Solucin:
Probemos yp = A cos(3x) + B sen(3x)
Tras sustituir,
(8 A 3B ) cos(3 x) (3 A 8 B ) sin(3 x) 2 sin(3 x)

Luego

A 6/73, B 16/73
6
16
y p cos(3x) sen(3 x)
73
73
43

Resolver
Solucin:

y"2 y '3 y 4 x 5 6 xe2 x

yc c1e x c2e3 x

Solucin homognea

Pensando en el principio de superposicin:

Probemos y p Ax B Cxe 2 x Ee 2 x
Tras sustituir,
3 Ax 2 A 3B 3Cxe2 x (2C 3E )e 2 x

Luego

4 x 5 6 xe2 x

A 4/3, B 23/9, C 2, E 4/3

4
23
4 2x
2x
y p x 2 xe e
3
9
3
4
23
4 2x
x
3x
y c1e c2e x 2 x e
3
9
3

44

Determinar una yp de

y"5 y '4 y 8e x

Solucin:
Probemos: yp = Aex
Tras sustituir: 0 = 8ex (conjetura incorrecta)
El problema est en que la funcin complementaria es:
Y la suposicin ya est presente en yc.

yc c1e c2 e
x

Probemos como alternativa: yp = Axex.

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

45

4x

 Si ninguna funcin en la supuesta y p es parte de yc

En la siguiente tabla se muestran soluciones particulares de prueba.

46

Hallar la forma de yp de
3 x
x
(a) y"8 y '25 y 5 x e 7e
Solucin:
Tenemos que g ( x) (5 x 3 7)e x y probamos con
3

y p ( Ax Bx Cx E )e

No hay duplicacin entre los trminos yp e yc

(b) y + 4y = x cos x
Solucin:
Probamos con x ( Ax B ) cos x (Cx E ) sin x
p
Tampoco hay duplicidad entre los trminos yp y
yc .
47

y 9 y 14 y 3 x 2 5sen2 x 7 xe 6 x
Hallar la forma de2 yp de

y p1 Ax Bx C

Solucin:
Para 3x2:

y p2 E cos 2 x Fsen2 x
y p3 (Gx H )e6 x

Para -5 sen 2x:

Para 7xe :
6x

y p y p1 y p2 y p3

Ningn trmino de

duplica un trmino de yc
48

As que la regla formal en este caso es que la solucin

particular es una combinacin lineal de las
funciones linealmente independientes que se
generan mediante diferenciaciones repetidas de
g(x).
Y cul es la regla si la solucin particular as propuesta
es tambin una solucin de la ecuacin homognea
asociada?
Si alguna yp contiene trminos que
duplican los trminos de yc, entonces
esa yp se debe multiplicar por xn,
donde n es el entero positivo ms
pequeo que elimina esa duplicacin.

49

y" y 4 x 10 senx, y ( ) 0, y ' ( ) 2

yc c1 cos x c2 senx
Resolver
Solucin:
Primero probamos: yp = Ax + B + C cos x + E sen x
Pero hay una duplicacin.
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
50

2
3
x
y
x
B
E
e
p

A
yp1xC
y
p
2
2
2
3
x
pyp1yp2

y
A
xB
xe
C
E

2
3x
y
"

6
y
'

9
y

6
x

12
e
Resolver

Solucin:
yc = c1e3x + c2xe3x

Debemos probar con:

Este trmino est

duplicado, aparece ya
en yc.

Tras sustituir y simplificar,

A = 2/3, B = 8/9, C = 2/3, E = -6
Luego
2 2 8
2
3x
3x
2 3x
y c1e c2 xe x x 6 x e
3
9
3

51

Resolver y y" e cos x

x

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

52

x
y

B
x
e

C
x
e

E
e
pyp1yp2
3
ypA
xyp13B
ex
xyp22exE
C
xex
y

( 4)

2 x

y 1 x e

Hallar la forma de yp de

Solucin:
yc = c1+ c2x + c3x2 + c4e-x
Prueba:

Como aparece repetido en la solucin homognea,

necesitaremos multiplicar A por x3 y (Bx2e-x + Cxe-x
+ Ee-x) por x. Prueba ahora:

53

Mtodo del anulador

Sigue los apuntes de Jose Olarrea.

54

Mtodo de variacin de
parmetros

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

Conocidas y1(x) e y2(x) soluciones l. i. de
la ec. homognea asociada, probaremos
como solucin particular:

y p u1 ( x) y1 ( x) u2 ( x) y2 ( x)

55

y p u1 ( x) y1 ( x) u2 ( x) y2 ( x)

yp P ( x) yp Q ( x) y p

Sustituimos yp, yp en la EDO:

u1[ y1 Py1 Qy1 ] u2 [ y2 Py2 Qy2 ]

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

d
d
[ y1u1 ] [ y2u2 ] P[ y1u1 y2u2 ] y1u1 y2 u2
dx
dx
d
[ y1u1 y2u2 ] P[ y1u1 y2u2 ] y1u1 y2 u2 f ( x)
dx
56

d
[ y1u1 y2u2 ] P[ y1u1 y2u2 ] y1u1 y2 u2 f ( x)
dx
Necesitaremos dos ecuaciones para encontrar
valores de u1 y u1. Exijamos que: y1u1 + y2u2 = 0,
para obtener una ecuacin adicional y de paso que
la EDO se reduzca a: y1u1 + y2u2 = f(x).
De modo que nos queda el sistema de ecuaciones:
y1u1 + y2u2 = 0
y1u1 + y2u2 = f(x)

57

W
y
f
(
x
)
W
y
f
(
x
)
2
1
1
2
u2W
u1W
Expresado en trminos de determinantes
y

y1
W
y1

donde

y2
0
, W1
y2
f ( x)

y2
y1
, W2
y2
y1

0
f ( x)

De donde encontraremos, por integracin, las

soluciones.
58

W
y
f
(
x
)
W
y
f
(
x
)
2
1
1
2
u2W
u1W
y"4 y '4 y ( x 1)e

2x

Solucin:
m2 4m + 4 = 0, m = 2 (cero doble)
y1 = e2x, y2 = xe2x,

Resolver
W (e 2 x

2x
e
, xe2 x )
2e 2 x

xe2 x
4x

e
0
2x
2x
2 xe e

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

0
W1
( x 1)e 2 x

xe2 x
e2 x
4x
( x 1) xe , W2 2 x
2x
2 xe
2e

( x 1) xe
u1
4x
e

4x

0
4x

(
x

1
)
e
( x 1)e 2 x

( x 1)e
x x , u2
4x
e
2

4x

x 1

59

u1 x x ,

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

2x
x p x x e x x xe x e x e
2
6
2
3
2

1 3 2x 1 2 2x
y yc y p c1e c2 xe x e x e
6
2
2x

2x

60

4 y"36 y csc 3 x
Resolver
Solucin:
y + 9y = (1/4) csc 3x
m2 + 9 = 0, m = 3i, -3i
y1 = cos 3x, y2 = sin 3x, f(x) = (1/4) csc(3x)

cos 3 x
sin 3 x
Como
W (cos 3 x , sin 3 x)
3
3 sin 3 x 3 cos 3 x
0

sin 3 x

cos 3 x
0
1
1 cos 3 x
W1
, W2

4
4 sin 3 x
1/4 csc 3 x 3 cos 3 x
3 sin 3 x 1/4 csc 3 x

61

W1
1
u1

W
12

W2 1 cos 3x
u2

W 12 sen3 x

Entonces u1 1/12 x,

u2 1/ 36 ln | sen3 x |

1
1
y p x cos 3 x ( sen3 x) ln | sen3 x |
12
36
1
1
y yc y p c1 cos 3 x c2 sen3 x x cos 3 x ( sen3 x) ln | sen3 x |
12
36

62

1
y" y
x
Resolver
Solucin:
x
t
x e
e
(
1
/
x
)
1
u1 -1
dt
m2u1
1 = 0, m =, 1,
x0 t

2
2
x
-x
y1 = e , yx2 = e , f(x) = 1/x,t y W(ex, e-x) = -2
Luego u e (1 / x) , u 1 x e dt
2

2 x 0 t

1 x x e t
1 x x et
yp e
dt e
dt
x0 t
2 x0 t
2
t
t
x
x
1
e
1
e
y yc y p c1e x e x
dt e x
dt
x0 t
2 x0 t
2
63

Ecuaciones de orden superior

Para las EDs de n-simo orden de la forma

y ( n ) Pn1 ( x) y ( n1) P1 ( x) y P0 ( x) y f ( x)

tomamos yp = u1y1 + u2y2 + + unyn, donde yi , i = 1, 2, , n, son la

familia de soluciones independientes que forman yc. As:

Suposiciones
para simplificar
la EDO:

y1u1 y2u2 ynun 0

y1u1 y2 u2 yn un 0

Que nos lleva a las ecuaciones solucin uk = Wk/W con k = 1, 2, , n. Donde

( n 1) de la y's
( ny
1W
) es el determinante
( n 1) que se obtiene de
W es el wronskiano
k
1
1
2
2
n
n
sustituir en W la k-sima columna por (0, 0,..., f(x)).

u y

u y

u f ( x )

64

Ecuacin de CauchyEuler
Forma de ecuacin de Cauchy-Euler
n
n 1
d
y
d
y
dy
n
n 1
an x
an1 x
a1 x a0 y g ( x)
n
n 1
dx
dx
dx

Mtodo de solucin
Probamos y(x) = xm, donde debemos determinar m,
para
resolver la ecuacin homognea asociada: Observa que:

k
d
y
k
mk
k
ak x m(m 1)(m 2)(m k 1) x
ak x
k
dx
ak m(m 1)(m 2)(m k 1) x m

an m(m 1)(m 2) (m n 1) ... a1m a0 x m 0

65

m
m
1
y
c1x
c2x2

Ecuacin auxiliar

2
d
y
dy
2
ax
bx cy g ( x)
2
dx
dx

Observa que
Para n = 2, y = xm, tenemos
tenemos que ax
(am(m 1) + bm + c)xm = 0,
o a cero en
es igual
x = 0. Para
2
am + (b a)m + c = 0
asegurar
existencia y
Caso 1: Races reales y distintas
unicidad,
2

tomaremos
I = (0, ).

2
d
y
dy
2
x
2x 4 y 0
2
dx
dx

Resolver
Solucin:
Tenemos a = 1, b = -2 , c = -4
m2 3m 4 = 0, m = -1, 4,
y = c1x-1 + c2x4

66

m
m
1
yycx
c1/21xc2
cx1/22xln1xln

Caso 2: Races reales repetidas

m1

Dedujimos y2 x ln x
Luego

2
d
y
dy
2
Resolver 4 x 2 8 x y 0
dx
dx

Solucin:
Tenemos a = 4, b = 8, c = 1
4m2 + 4m + 1 = 0, m = - , -

67

Caso 3: Races complejas

conjugadas

Orden superior: multiplicidad k

m1

m1

m1

m1

x , x ln x , x (ln x) , , x (ln x)

k 1

Caso 3: races 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)]

68

Resolver
Solucin:

/y
1
2

xcos(2lnx)

1
4 x y 17 y 0, y (1) 1, y ' (1)
2
2

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,

69

3
2
d
y
d
y
dy
3
2
x
5x
7x 8y 0
3
2
dx
dx
dx

Resolver
Solucin: dy
2
m 1 d y
m2
m

mx
,

m
(
m

1
)
x
,
Sea y = x ,
2
dx
dx

d3y
m 3

m
(
m

1
)(
m

2
)
x
3
dx
Luego tenemos 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)
70

x 2 y"3 xy '3 y 2 x 4e x
Resolver
Solucin:
Tenemos (m 1)(m 3) = 0, m = 1, 3
yc = c1x + c2x3
Usando variacin de parmetros,
yp = u1y1 + u2y2, donde y1 = x, y2 = x3
Escribimos la ED como

3
3
2 x

y y 2 y 2x e
x
x

Luego P(x) = -3/x, Q(x) = 3/x2, f(x) = 2x2ex

71

x x3
3
W
2x ,
2
1 3x
0
W1
As
2 x 2e x

x
0
x3
5 x
3 x

2
x
e
,
W

2
x
e
2
2 x
2
1 2x e
3x

5 x
2 x 5e x
2
x
e
2 x
x

u1

x
e
,
u

e
2
3
2x
2 x3
x
2 x
x
x
u

e
u1 x e 2 xe 2e , 2

Hallamos

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
3

2 x

y yc y p c1 x c2 x 2 x e 2 xe

x
72

Una ecuacin 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 y xy y ln x
2

x et
t ln x
dy dy dt 1 dy

dx dt dx x dt
d 2 y d 1 dy
1 dy 1 d dy
1 dy 1 d dy



2

2
dx
dx x dt
x dt x dx dt
x dt x dt dx
1 dy 1 d 1 dy
2

x dt x dt x dt

1 d 2 y dy
2

2
x dt
dt

73

x 2 y xy y ln x
d2y
dy
2 y t
2
dt
dt

y c1e t c2te t 2 t

y c1 x c2 x ln x 2 ln x
74

Unos ejemplos de
ecuaciones no lineales
2

y" 2 x( y ' )
Resolver
Solucin:
Esta ecuacin no lineal carece de trmino en y.
Sea u(x) = y, entonces du/dx = y,
du
2xu 2
dx

du
1
2
2

2
x
dx
1
2
u
(Se escribe en esta forma solo por conveniencia para luego
integrar)

dy
1
2 2
dx
x c1
Entonces,
dx
1
1 x
y 2 2 tan
c2
c1
c1
x c1
Como u-1 = 1/y,

75

yy" ( y ' )

Resolver
Solucin:
Esta ecuacin no lineal carece de trmino en x.
Sea u(x) = y, entonces
du = (du/dy)(dy/dx)
du dy
2
y y= du/dx
= u du/dy
u u

dy

c2 e

ln|u| = ln|y| + c1, u = c2y (donde

Como u = dy/dx = c2y, dy/y = c2 dx
ln|y| = c2x + c3,
c2 x

c1

y c4e

76

y x y y 2 ,

y (0) 1 ,

y(0) 1

Supongamos que existe solucin para:

y(0)
y(0) 2 y(0) 3 y ( 4) (0) 4
y ( x) y (0)
x
x
x
x
Si adems suponemos
que
y(x)
admite
1!
2!
3!
4!
desarrollo en serie de Taylor centrado en 0:

d
2

y ( x) ( x y y ) 1 y 2 yy
Como y(0) = -1,
dx y(0) = 1, de la ED original:
y(0) = 0 + y(0) y(0)2 = 2.
Derivando sucesivamente la ED original:

77

d
y ( x) (1 y 2 yy) y 2 yy 2( y) 2
dx
( 4)

... podemos utilizar el mismo mtodo para

obtener y(3)(0) = 4, y(4)(0) = 8, etc.
Y encontrar una aproximacin en Taylor de
la solucin:

2 3 1 4
y ( x) 1 x x x x
3
3
2

78

Una ltima observacin: La ED de este

ejemplo: y x y y 2 , y (0) 1 , y(0) 1
es equivalente (mediante cambio de
variable) al sistema de ecuaciones
diferenciales:
dy

u
dx

du
2
x y y

dx
y (0) 1, u (0) 1

79