Ecuacionesdiferenciales LHV 0

Universidad Tecnica Federico Santa Mara
Departamento de Matematica
Ecuaciones Diferenciales Ordinarias
1. Introduccion
Definicion 1.1 Una ecuacion diferencial ordinaria (EDO) es una relacion del siguiente tipo

F x, y, y (x), y (x), , y (n) (x) = 0, (1)
donde la incognita y = y(x) es una funcion de la variable independiente x.
La funcion F establece una relacion entre las derivadas de la incognita y. El adjetivo ordinaria
hace referencia a que solo se deriva con respecto a una variable.
Definicion 1.2 El orden de una ecuacion diferencial es el numero que indica las veces que
aparece derivada la variable dependiente y con respecto a la variable independiente x.
Por ejemplo, y + sen(x + y) = x2 2xy, es de orden 2.
Definicion 1.3 El grado de una EDO es la potencia a la cual esta elevada la derivada de
mayor orden.
Por ejemplo, xy + cos(x2 + y 3 ) = x2 2x3 y, es de orden 2 y de grado 1.
Definicion 1.4 La ecuacion 1 es lineal en la variable y si ella se puede escribir de la forma:
an (x)y (n) + an1 (x)y (n1) + + a1 (x)y + a0 (x)y = q(x), (2)
Donde las funciones ai (x) R, i = 1, . . . , n son llamados los coeficientes de la EDO. Si

estos coeficientes dependen explcitamente de la variable independiente, se dice que la EDO es
de coeficientes variables. Si no, se dice que la EDO es de coeficientes constantes.
Por ejemplo, la ecuacion
cos(x)y + x2 y (x) + 2xy = tan x,
es lineal, de segundo orden, de primer grado y con coeficientes variables. Mientras que la
ecuacion
2y + 3y + y + 5y = ex ,
es lineal, de tercer orden, primer grado y con coeficientes constantes.
Una EDO no lineal, como su nombre lo indica, es una ecuacion que no es lineal(ya que tiene
una gran cantidad de formas distintas). Por ejemplo, la ecuacion
y(2 + (y )3 ) = 1,
es no lineal, de orden uno y grado tres.
2. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
2. Ecuaciones Diferenciales de Primer Orden
En esta seccion estudiaremos algunas tecnicas que nos permitiran determinar soluciones de
una gran cantidad de EDO.
2.1. Integracion Directa
El problema tipo es
y = f (x). (3)
Si se tiene que f es integrable, podemos usar el Teorema Fundamental del Calculo (TFC) para
obtener soluciones dadas por Z
y(x) = f (x) dx + C, (4)
Esta es una solucion general de la ecuacion 3 incluyendo la constante arbitraria C R. El valor
de C se puede determinar mediante una condicion inicial.
Definicion 2.1 Una condicion inicial es un punto cualquiera (x0 , y0 ) por donde pasa la
solucion de la EDO, la que permite calcular el valor de C cuando x = x0 e y = y0 .
Teorema 2.1 Todo
Ejemplo 2.1 las soluciones de la ecuacion y = cos x son
Z
y = cos x dx = sen x + C, C R.
Supongamos ahora que queremos encontrar la solucion de 3 definida sobre un intervalo I y

que ademas pase por un punto (x0 , y0 ) dado. Esto es, queremos resolver el siguiente problema
con una condicion inicial : (
y = f (x), x I
y(x0 ) = y0 .
Si integramos lo anterior entre x0 y x I, se obtiene:
Z x Z x

y (z)dz = f (z)dz,
x0 x0
usando el TFC sigue que Z x
y(x) = y(x0 ) + f (z)dz. (5)
x0
Comparando 4 con 5 se deduce que la constante arbitraria en 4 queda determinada en forma
unica como C = y(x0 ).
2
2.2. VARIABLES SEPARABLES
2.2. Variables Separables
Definicion 2.2 Estas EDOs se escriben como
y = f (x)g(y). (6)
Note que si en algun punto y0 ocurre que g(y0) = 0, entonces la funcion constante y(x) = y0
es solucion de 6. En aquellos puntos donde g(y) 6= 0 tenemos que
y
= f (x).
g(y)
Integrando y gracias al cambio de variable y = y(x):
Z
1 dy dy
dx = (y = y(x) dy = dx)
g(y) dx dx
Z Z
dy
= f (x) dx + C, (7)
g(y)
donde C R es una constante. Como veremos, no siempre es sencillo encontrar una expresion
explcita para la solucion y(x) desde 7.
Ejemplo 2.2
y = xy.
En este caso, f (x) = x y g(y) = y. Note que y(x) = 0 es una solucion. Si y 6= 0, entonces
Z Z
dy
= x + C.
y
De esta manera
x2
ln |y| = + C,
2
y luego
x2
|y| = ke 2 ,
de esto se desprende que k = eC > 0 y luego las soluciones son positivas. Ademas, si k = 0 se
recupera la solucion nula. De esta manera
x2
y = ke 2 .
Ejemplo 2.3
2y+3
y = .
x y
3
2 y+3
En este caso f (x) = x
y g(y) = y
. En este caso, la funcion constante y(x) = 3 es una
solucion de la EDO. Ademas, existe la restriccion y 6= 0. Si y 6= 3, entonces
Z Z
3 2
1 dy = dx
y+3 x
y + ln(y + 3)3 = ln x2 + C
(y + 3)3 ey = kx2 ,
donde k = eC > 0. Como vemos, en este caso la solucion es implcita.
Ejemplo 2.4 Se ha establecido que la velocidad de la desintegracion del atomo del radio es
directamente proporcional a su masa en cada instante. Determinar la ley de variacion de la
masa del radio en funcion del tiempo, si para t = 0 la masa es de m0 .
Sea m = m(t) la masa en el instante t. Luego, por la desintegracion del elemento se tiene que
la masa en el instante t + t es m m. Entonces
m(t + t) m(t) m
= .
t t
En el lmite cuando t 0 se tiene dm
dt
. Ademas, de acuerdo a la hipotesis:
dm
= km, k > 0.
dt
Integrando, se obtiene m(t) = cekt . Si imponemos la condicion inicial se obtiene que c = m0
y entonces
m(t) = m0 ekt .
2
Ejemplo 2.5 Resolver la EDO xy 3 dx + ex dy = 0 separando las variables queda:
Z
x2 3 1 2 1
xe dx + y dy = 0 2xex dx y 2 = C
2 2
1 2 1 2
ex y 2 = C ex + y 2 = C
2 2
Ejemplo 2.6 Resolver la EDO sen(x)sen(y)dx + cos(x)cos(y)dy = 0 separando las variables
queda:
senx cosy
dx + dy = 0 lncosx + lnseny = lnC
cosx seny
seny
= C seny = Ccosx
cosx
4
2.2.1. Ecuaciones Diferenciales Homogeneas de Primer Orden
Como veremos, este tipo de ecuaciones no es de variables separables, pero usando alguna
tecnica algebraica (una sustitucion adecuada) podremos reescribirlas de dicha manera. Primero,
definimos
Definicion 2.1. Sea f : R2 R una funcion de dos variables y k N. Se dice que la funcion
f (x, y) es homogenea de grado k si
f (x, y) = k f (x, y), R \ {0}.
Notemos que si contamos con dos funciones f, g homogeneas de grado k, g 6= 0. Entonces,

para x 6= 0
y y y
f (x, y) f (x, x ) xk f (1, ) f (1, )
= x = x = x =h y .
g(x, y) y y y x
g(x, x ) xk g(1, ) g(1, )
x x x
Es decir, que el cuociente de dos funciones homogeneas es una funcion que depende solo de la
razon entre x e y.
Una EDO de primer orden es homogenea si puede ser escrita como
y
y = h . (8)
x
Observacion 2.1 Es importante no confundir con las EDOs lineales homogeneas, que estu-
diaremos mas adelante.
Para resolver una ecuacion homogenea, introducimos el cambio de variable y = vx. Notamos
que
y = xv + v.
Si introducimos esto en 8 se obtiene
xv + v = h(v),
y luego
1
v = (h(v) v) , (9)
x
que es una ecuacion diferencial de variables separables.
5
Ejemplo 2.7 Consideremos la EDO

x+y
y = .
xy
Note que las funciones del lado derecho son homogeneas de grado 1. Si introducimos el cambio
de variable y = vx, entonces
x + xv
xv + v =
x xv
1 + v2
xv = , x 6= 0
1v
1 1 + v2
v = ,
x 1v
que como vemos es de variables separables. Resolviendo
1v dx
dv =
1 + v2 x
1
arctan v ln(1 + v 2 ) = ln |x| + C.
2
Volviendo a las variables originales, la solucion en forma implcita queda expresada como
y 1 y 2
arctan ln 1 + = ln |x| + C.
x 2 x
Ejemplo 2.8
x2 + y 2
y = .
2xy
En este caso, las funciones del lado derecho son homogeneas de grado 2. Nuevamente, intro-
duciendo el cambio de variable y = vx, se tiene
x2 + x2 v 2
xv + v =
2x2 v
1 v2
xv = , x 6= 0
2v
1 1 v2
v = ,
x 2v
que es de variables separables. Resolviendo
2v dx
dv =
1 v2 x
2
ln(1 v ) = ln |x| + C
(1 v 2 )1 = k|x|.
6
Volviendo a las variables originales
y 2 = x2 k1 |x|.
Ejemplo 2.9
x2 xy + y 2
y =
xy
En este caso, las funciones del lado derecho son homogeneas de grado 2. Nuevamente, intro-
duciendo el cambio de variable y = vx, se tiene
x2 x2 v + x2 v 2
xv + v =
x2 v
1v
xv = , x 6= 0
v
11v
v = ,
x v
v dx
dv =
1v x
v ln(v 1) = ln |x| + C
Volviendo a las variables originales

y yx
ln = ln|x| lnC
x x
y
ln(y x) = lnC
x
y
(y x)e x = C
2.2.2. Ejercicios Propuestos
1. (x2 y 2 )y = xy 7. (y a)dt + t2 dy = 0
2. (2x + 3y)dx + (x 2y)dy = 0 8. xdt (t2 a2 )dx = 0

dx 1+x2
3. ydx xdy = 0 9. dy
= 1+y 2

4. (1 + u)vdu + (1 v)udv = 0 10. (1 + s2 )dt tds = 0
5. (1 + y)dt (1 t)dy = 0 11. d + tg d = 0
6. (t2 xt2 ) dx
dt
+ x2 + tx2 = 0 12. sec2 tg d + sec2 tg d = 0
7
p
13. (1 + x2 )dy 1 y 2 dx = 0 15. xy 2 dy = (x3 + y 3 )dx
14. 3ex tgy dx + (1 ex )sec2 y dy = 0 16. tgx sen2 y dx + cos2 x ctgy dy = 0
2.2.3. Ecuaciones Reducibles a homogeneas.
Si se tiene que la ecuacion posee la forma

ax + by + c
y = , (10)
x + y +
donde el segundo miembro es la razon entre dos rectas. Se tienen dos posibilidades.
a) a/b 6= /. Es decir, rectas son no paralelas. Se busca el punto de interseccion (x0 , y0 )

y se reescribe 10 como
a(x x0 ) + b(y y0 )
y = . (11)
(x x0 ) + (y y0 )
Ahora se introduce el cambio de variable
(
u = x x0
v = y y0 ,
dy dv dv
se tiene que du = dx y ademas dx
= dx
= du
. Reemplazando en 11 se obtiene
dv au + bv
= ,
du u + v
que es homogenea de grado 1.
Ejemplo 2.10 Consideremos la siguiente ecuacion
dy 2x + y + 2
= .
dx xy+1
Evidentemente, las rectas del lado derecho no son paralelas. El punto de interseccion es
(x0 , y0 ) = (1, 0). Reescribimos la ecuacion como
dy 2(x (1)) + (y 0) 2(x + 1) + y
= = .
dx (x (1)) (y 0) (x + 1) y
Ahora introducimos el cambio de variable u = x + 1, v = y y repitiendo el procedimiento
anterior llegamos a
dv 2u + v
= ,
du uv
que es homogenea de
grado 1.Se deja como ejercicio la resolucion.
y
Respuesta: 2 arc tg = ln[C(x2 + 2x + 1 + y 2 )]
x+1
8
2.3. ECUACIONES DIFERENCIALES EXACTAS
b) a/b = /. Las rectas son paralelas, pero esto implica que

a
= = = a = b.
b
Luego
dy (ax + by) + c
= .
dx (ax + by) +
Hacemos el cambio de variable
dz dy dy 1 dz a
z = ax + by =a+b = .
dx dx dx b dx b
Reemplazando
dz b(z + c)
= + a,
dx z +
que es de variables separables.
Ejemplo 2.11 Consideramos

2x + y 1 2x + y 1
y = y = .
4x + 2y + 5 2(2x + y) + 5
Hacemos z = 2x + y. Luego z = 2 + y . Reemplazando
z1 5z + 9
z = + 2 z = ,
2z + 5 2z + 5
Z Z
2z + 5
dz = dx + C.
5z + 9
Mediante el cambio de variable u = 5z + 9 es posible integrar para obtener
2 7
(5z + 9) + ln |5z + 9| = x + C.
25 25
Finalmente
2 7
(5(2x + y) + 9) + ln |5(2x + y) + 9| = x + C.
25 25
2.3. Ecuaciones Diferenciales Exactas
Dada una funcion de dos variables F (x, y), su diferencial total, dF , esta dado por
F F
dF = dx + dy,
x y
9
F F
siempre y cuando ,
x y
existan. Puesto que las derivadas parciales de F tambien son funciones
de x e y es posible escribir lo anterior como
dF = M(x, y)dx + N(x, y)dy.
Consideremos ahora la curva de nivel F (x, y) = c R definida sobre una region del plano R.
Si suponemos que F C 1 (R), entonces
df = M(x, y)dx + N(x, y)dy = 0, (x, y) R. (12)

F
Si ademas, exigimos que y
6= 0 R, entonces el Teorema de la Funcion Implcita implica
que para cada punto de R, existe una funcion y = y(x) de forma unica y diferenciable donde
ademas
M(x, y)
y = ,
N(x, y)
de donde y es una solucion de 12.
Se dice a la expresion M(x, y)dx + N(x, y)dy que es una diferencial exacta si existe una funcion
F tal que dF = M(x, y)dx + N(x, y)dy. En tal caso, la ecuacion 12 se denomina ecuacion
diferencial exacta y la expresion F (x, y) = c, c arbitraria, define la solucion general de 12 en
toda la region R, llamada la integral general de esta ecuacion.
Ejemplo 2.12 La ecuacion

xdx + ydy = 0,
es exacta, ya que
x2 + y 2 x2 + y 2
= x, = y.
x 2 y 2
De donde la integral de esta ecuacion es
F (x, y) = x2 + y 2 = 2c > 0.
Sin embargo, no contamos con una tecnica que nos permita reconocer una ecuacion dife-
rencial exacta para abordar su resolucion. Esto lo podemos hacer simplemente exigiendo a F
que sea una funcion dos veces diferenciable con segundas derivadas parciales continuas en R,
es decir, F C 2 . Si esto ocurre, entonces
2F 2F
= ,
yx xy
10
F F
donde hemos definido a M(x, y) = x
y N(x, y) = y
. Luego, esta ultima condicion puede
escribirse como
M N
= . (13)
y x
Este criterio sera efectivo para determinar la exactitud de una ecuacion si el recproco de 13 es
verdadera. Para tal fin, es necesario anadir una hipotesis adicional. Esto es, se debe exigir que
R sea una region simplemente conexa, concepto que sera profundizado en el curso siguiente.
Informalmente hablando, diremos que R es simplemente conexa si no posee agujeros en su
interior. Enunciamos el siguiente resultado:
Teorema 2.1. Sean M(x, y), N(x, y) continuamente diferenciables en una region simplemente
M N
conexa R. Entonces M(x, y) dx + N(x, y) dy es un diferencial exacta si y solo si y
= x
.
Ejemplo 2.13 La ecuacion (x2 + y)dx + (x 2y)dy = 0 es exacta pues

M 2 N
= (x + y) = 1 = (x 2y) = 1.
y y x y
Buscamos una curva solucion de la forma F (x, y) = c. Para resolver, recordamos que por
F
definicion: x
= M y luego, integrando con respecto a x considerando a la variable y constante
F
= (x2 + y)
x Z
F (x, y) = (x2 + y) dx + A(y)
x3
F (x, y) = + xy + A(y),
3
donde A(y) es una funcion que solo depende de y. Derivando la ultima expresion y comparando
con N:
F dA
= x+
y dy
dA
N(x, y) = x +
dy
dA
x 2y = x + A(y) = y 2 + C.
dy
Entonces la curva solucion es
x3
F (x, y) = + xy y 2 = c.
3
Ejemplo 2.14 Resolver la ecuacion:
(3x2 y + 2xy)dx + (x3 + x2 + 2y)dy = 0
11
Sera exacta? si pues
M(x, y) = 3x2 y + 2xy, N(x, y) = x3 + x2 + 2y

M N
= 3x2 + 2x = 3x2 + 2x
y x
Por lo tanto, podemos ver que corresponde a una ecuacion diferencial exacta
M N
= = 3x2 + 2x
y x
F
= 3x2 y + 2xy F (x, y) = x3 y + x2 y + A(y)
x
F dA dA
= 3x3 + 2x2 + = x3 + x2 + 2y = 2y A(y) = y 2
y dy dy
Por lo tanto,
F (x, y) = x3 y + x2 y + y 2
y la solucion es
x3 y + x2 y + y 2 = C
Ejemplo 2.15 Resolver

2x y 2 3x2
dx + dy = 0
y3 y4
La EDO es exacta pues:
d 2x 6x d y 2 3x2
= 4 =
dy y3 y dx y4
dF 2x x2
= 3 F (x, y) = 3 + A(y)
dx y y
2 2
dF d x 3x dA y 2 3x2 dA 1
= + A(y) = + = =
dy dy y 3 y4 dy y4 dy y2
1 x2 1
A(y) = F (x, y) = 3
y y y
Luego la solucion general de la EDO exacta es
x2 1
=C
y3 y
12
Ejemplo 2.16 Resolver la EDO siguiente

y2 1 1 x2
dx + dy = 0
(x y)2 x y (x y)2
| {z } | {z }
M (x,y) N (x,y)
La EDO es exacta pues:

dM d y2 1 2xy d 1 x2 dN
= 2
= 3
= 2
=
dy dy (x y) x (x y) dx y (x y) dx
Por lo tanto
Z
dF y2 1 y2 1 y2
= 2
2
dx = ln|x| + A(y)
dx (x y) x (x y) x xy

dF d y2 2xy y 2 dA 1 x2
= ln|x| + A(y) = + =
dy dy xy (x y)2 dy y (x y)2
dA 1 y2 y2 y
= 1 + A(y) = ln|y| y F (x, y) = ln|x| + ln|y| y = + ln| | y
dy y xy xy x
Luego la solucion general de la EDO exacta es
y y2
ln| | y = C
x xy

1. (x2 + y)dx + (x 2y)dy = 0 1 3y 2 2y
6. + dx = dy
x2 x4 x3
2. (y 3x2 )dx (4y x)dy = 0
x2 dy y 2 dx
3. (y 3 x)y = y 7. =0
(x y)2
4. 2(3xy 2 + 2x3 )dx + 3(2x2 y + y 2)dy = 0 8. (3x2 y 2 )dy 2xydx = 0
xdx + (2x + y)dy
5. =0
(x + y)2 9. 2xydx + (1 + x2 )dy = 0
2.3.2. Factor de Integracion
Supongamos que la ecuacion 12 no es exacta. Para poder escribirla en dicha forma, se

multiplica la ecuacion por una funcion = (x, y) 6= 0 (llamando a usualmente factor
13
integrante) de modo que la ecuacion
(M) dx + (N) dy = 0, (14)
resulte ser exacta. Al realizar un calculo sencillo se llega a la ecuacion

N M
M N = ,
y x x y
que es una ecuacion diferencial parcial de primer orden, posiblemente mas compleja que el
problema inicial. Sin embargo, si por ejemplo, el factor integrante lo suponemos que depende
solo de y, entonces

d N M
M =
dy x y

1 d 1 N M
=
dy M x y

d 1 N M
(ln ) =
dy M x y
Z
1 N M
ln = dy
M x y
Puesto que en este caso se ha supuesto que = (y), para que eso suceda se debe cumplir
que el lado derecho de la ultima ecuacion no depende de x y la ecuacion resultante es una EDO
exacta al ser multiplicada por el factor integrante encontrado.
De forma similar, si = (x), debe verificarse que la expresion

1 M N
,
N y x
solo dependa de x, por lo tanto, se tiene que
Z
1 M N
ln = dx
N y x
Ejemplo 2.17 Dada la ecuacion (y + xy 2 ) dx x dy = 0, se tiene que
M N
= 1 + 2xy, = 1
y x
y la ecuacion no es exacta, pero si consideramos la razon

1 N M 1 1 2xy 2(1 + xy) 2
= 2
= = .
M y x y + xy y(1 + +xy) y
14
Luego el factor integrante depende de y. Este factor confirma la EDO

d 2
(ln ) = ,
dy y
1
de donde ln = 2 ln y y entonces = es un factor integrante. Multiplicando la ecuacion
y2
original por este factor se obtiene

1 x
+ x dx 2 dy = 0,
y y
que es exacta (siempre es conveniente verificarlo)

1 1 x
+x = 2 = 2 .
y y y x y
x x2
Donde la solucion es + =C
y 2
Ejemplo 2.18 Resolver la ecuacion

2 y3
2xy + x y + dx + (x2 + y 2 )dy = 0
3
Aqu vemos que
1 M N 2x + x2 + y 2 2x
= =1
N y x x2 + y 2

d 1 M N
Por lo tanto, = dx = dx ln() = x = ex multiplicando por
N y x
= ex se obtiene:
x 2 y3
e 2xy + x y + dx + ex (x2 + y 2)dy = 0
3
que resulta ser exacta y cuya solucion es :

x 2 y2
ye x + =C
3
15

x x x
1. x + e y dx + e y 1 dy = 0 de modo que satisfaga la condicion inicial y(0) = 2.
y

ln(lny) 2xy 3 lnx 2 2
2. + dx + + x y dy = 0
x 3 ylny
3. (xcosy yseny)dy + (xseny + ycosy)dx = 0

1 y x 1
4. dx + dy = 0
x x2 + y 2 x2 + y 2 y
5. (x + y 2 )dx (2xy)dy = 0 11. 2y 2dx + (2x + 3xy)dy = 0
6. y(1 + xy)dx = xdy

12. (x2 + 2y)dx + xdy = 0
y
7. dx + (y 3 lnx)dy = 0
x
13. (x2 + y 2 )dx + 3xydy = 0
9. ydx + (y x)dy = 0
10. (x2 + y 2 + x)dx + ydy = 0 14. xydx + (x2 + 2y 2 + 2)dy = 0
Hay otros casos de factor integrantes, por ejemplo, si = (x2 + y 2), se efectua un cambio
de variables z = x2 + y 2 en

N M
M N = ,
y x x y
quedando
z z N M
M N = ,
z y z x x y

d d N M
M 2y N 2x = ,
dz dz x y
de donde la expresion !
N M
x
y
2yM 2xN
solo dependa de z, o sea, es funcion de z, donde
Z N M
!
x
y
ln = dz
2yM 2xN
16
2.4. ECUACIONES DIFERENCIALES LINEALES
(x3 + xy 2 y)dx + (x2 y + y 3 + x)dy = 0
si sabe que el factor integrante es = (x2 + y 2 ), podemos ver que

N M
!
y
x 2xy + 1 (2xy 1) 1 1
= = =
2yM 2xN 2y(x3 + xy 2 y) 2x(x2 y + y 3 + x) x2 + y 2 z
R 1 1 1
Por lo tanto, = e z
dz
= elnz = = 2
z x + y2
1
Al multiplicar la EDO por = , y ordenandola convenientemente queda
x2 + y 2
ydx xdy
xdx + ydx =
x2 + y 2
arreglandola queda
ydxxdy dx xdy
y2 y
y2
xdx + ydx = x2 +y 2
= x2
y2
1+ y2
Lo que da
x2 y 2 x
+ = arctg +C
2 2 y
2.4. Ecuaciones Diferenciales Lineales
Una ecuacion de la forma

a1 (x) y + a0 (x) y = h(x), (15)
donde a1 (x), a0 (x) y h(x) son funciones continuas en un intervalo I, a1 (x) 6= 0 x I se
denomina Ecuacion Diferencial lineal de Primer Orden, y se dice normal si en 15
a1 (x) 6= 0 x I, dividiendo por a1 (x) 6= 0 queda:
y + p(x) y = q(x), (16)
donde p(x) = a0 (x)/a1 (x) y q(x) = h(x)/a1 (x).

A la solucion de 16 la llamaremos solucion general de la ecuacion . Para caracterizarla conta-
mos con el siguiente teorema:
Teorema 2.2. Principio de Linealidad. La solucion general de 16 esta dada por
y(x) = yh (x) + yp (x), (17)
17
donde yh = yh (x) es la solucion del problema homogeneo
y + p(x) y = 0 (18)
llamada solucion homogenea e yp = yp (x) es una solucion particular de
yp + p(x) yp = q(x).
Demostracion. Primero verificamos que y = y(x) dada por 17 es una solucion general de
16. En efecto, al reemplazar
y + p(x) y = q(x)
yh + yp + p(x)yh + p(x)yp = q(x)
(yh + p(x)yh ) + (yp + p(x)yp ) = q(x)
0 + (yp + p(x)yp ) = q(x)
yp + p(x)yp = q(x).
Ya que yh es solucion del problema homogeneo. Entonces, yp es la solucion particular y se tiene

que:
yp + p(x)yp = q(x),
es valida trivialmente y por lo tanto, y(x) dada por 17 es solucion.
Para verificar que todas las soluciones estan dadas por 17 nos daremos otra solucion y 6= y.
Puesto que la ecuacion es lineal, entonces la suma y + y tambien es solucion y entonces
(y + y) + p(x) (y + y) = q(x)
y + yh + yp + p(x)y + p(x)yh + p(x)yp = q(x)
(y + p(x)y) + (yh + p(x)yh ) + (yp + p(x)yp ) = q(x)
(y + p(x)y) + 0 + (yp + p(x)yp ) = q(x)
(y + p(x)y) + (yp + p(x)yp ) = q(x).
El segundo termino del lado izquierdo se anula pues es solucion del problema homogeneo, en
tanto que el tercer termino vale q(x) pues yp es solucion particular. Luego
y + p(x)y = 0,
lo que es una contradiccion al supuesto de que y sea solucion de 16. Por lo tanto, todas las
soluciones de 16 estan dadas por la identidad 17
18
El teorema anterior es de gran utilidad, pues ademas de ser aplicable a ecuaciones lineales
de orden superior, permite desarrollar un algoritmo para resolver 16. Siempre resolveremos
primero el problema homogeneo, notando que la ecuacion 18:
y + p(x) y = 0,
es de variables separables, con solucion dada por

R
yh (x) = ce p(x) dx
, c R. (19)
2.4.1. Variacion de Parametros
En segundo lugar, buscaremos la solucion particular, que asumiremos estara dada por
yp (x) = k(x)yh (x) (un multiplo funcional de la solucion homogenea). Si esto es as, solo resta
encontrar a la funcion k(x). Si reemplazamos a yp en la ecuacion se tiene que
(k(x)yh (x)) + p(x)(k(x)yh (x)) = q(x)

k yh (x) + kyh (x) + p(x)k(x)yh (x) = q(x)
k yh (x) + k(x)(yh (x) + p(x)yh (x)) = q(x)
k yh (x) + k(x)(0) = q(x)
k yh (x) = q(x)
q(x)
k = ,
yh (x)
recordando que yh (x) es solucion homogenea. Una solucion de la ultima ecuacion es
Z
q(x)
k(x) = dx.
yh (x)
Y luego la solucion particular es
Z R
R
p(x) dx
yp (x) = k(x)yh (x) = q(x)e dx e p(x) dx . (20)
Entonces, la solucion general de 16 es
y(x) = yh (x) + yp (x)

R
Z R R
p(x) dx p(x) dx
= ke + q(x)e dx e p(x) dx (21)
R
Z R

p(x) dx p(x) dx
= e c + q(x)e dx , c R. (22)
19
Ejemplo 2.20 La ecuacion (1 + x2 )y xy = 1, es lineal. En este caso, a1 (x) = (1 + x2 ),

a0 (x) = x y h(x) = 1. Note que 1 + x2 > 0 x R. Luego, podemos normalizar y resolver
en todo R. Al hacerlo, obtenemos
x 1
y y = .
1 + x2 1 + x2
La ecuacion homogenea asociada es
x
y y = 0.
1 + x2
Separando variables y resolviendo
x
y = 2
y
Z Z1 + x
dy x
= dx
y 1 + x2
ln y = ln(1 + x2 )1/2 + k

yh (x) = c 1 + x2 .
En tanto, la solucion particular es yp (x) = k(x)yh (x), donde k(x) esta dado por
Z Z Z
q(x) 1/(1 + x2 ) dx
k(x) = dx = dx = .
yh (x) 1+x 2 (1 + x2 )3/2
Mediante el cambio de variable x = tg , se llega a
x
k(x) = sen = . x
1 + x2
Luego 1
x
yp (x) = 1 + x2 = x.
1 + x2
La solucion general es

y(x) = c 1 + x2 + x.
(x4 + 2y)dx = xdy
Ordenandola
2
xy 2y = x4 y y = x3
x
20
R
Z R

(2/x) dx 3 (2/x) dx
y(x) = e C+ x e dx
Z Z 3
2 lnx 3 2 lnx 2 x
y(x) = e C+ xe dx = x C + dx
x2

x2 2
y(x) = x C + 2y = x2 C + x4
2
(3xy + 3y 4)dx + (1 + x)2 dy = 0
Ordenandola
3 4
(1 + x)2 y + 3(1 + x)y = 4 y + y=
1+x (1 + x)2
Z
R
3/(1+x) dx 4 R
3/(1+x) dx
y(x) = e C+ e dx
(1 + x)2
Z Z
3 ln(1+x) 4 3 ln(1+x) 1 (1 + x)3
y(x) = e C+ e dx = C +4 dx
(1 + x)2 (1 + x)3 (1 + x)2
C 4 (1 + x)2
y(x) = + = C(1 + x)3 + 2(1 + x)1
(1 + x)3 2 (1 + x)3
(1 + 2x tg y) dy = dx
Ordenandola
dx
2x tg y = 1
dy
R
Z R

2 tg y dy 2 tg y dy
x(y) = e C+ e dy
Z Z
2 ln|cosy| 2 ln|cos(y)| 2 2
x(y) = e C+ e dy = sec y C + cos y dy

2 1 sen(2y)
x(y) = sec y C + y+ 4xcos2 y = C + 2y + sen(2y)
2 2
21

di
L + Ri = E,
dt
donde L, R, E son constantes con la c.i. i(0) = 0
Ordenandola
R Z R
di R E t E t
+ i= i(t) = e L C + e L dt
dt L L L

R R R
t E L t t E
i(t) = e L C + e L = Ce L +
LR R

R
E E E t
i(0) = 0 0 = C + C = i(t) = 1 e L
R R R
Ejemplo 2.25 Encontrar la velocidad inicial mnima que debe tener un cuerpo que se dispara
para que escape de la atraccion de la tierra. Desprecie la resistencia del aire.
k
a(r) = , k<0
r2
como en
gR2
r = R, a(R) = g a(r) = 2
r
por otro lado
dv dv dr dv gR2 dv
a= = = v a(r) = 2 = v
dt dr dt dr r dr
gR2 v2
+C =
r 2
pero en
v02 2gR2
r = R, v = v0 C = gR v 2 = + v02 2gR
2 r
Ejemplo 2.26 Un embudo con angulo de salida de 60o y una seccion transversal con area
de 0, 5cm2 , contiene agua. En t = 0, se abre la salida. Determinar el tiempo que tardara en
vaciarse, suponiendo que la altura inicial del agua es de h(0) = 10cm.

Dato: v = 0, 6 2gh donde h =altura instantanea, v =velocidad del lquido.
dV =volumen de agua que sale.

dV = 0, 5 v dt con v = 0, 6 2gh, por otro lado
22

dV = r 2 dh con r = h tg(30)= h/ 3
h2 p 0, 9 2g
dV = dh = 0, 3 2ghdt dt = h3/2 dh
3

0, 9 2g 2 2
t = h5/2 + C h(t = 0) = 10 0 = 105/2 + C
5 5
2 2
t= 105/2 h5/2 h = 0 t = 105/2 = 99, 7seg 1 40
5 0, 9 2g 5 0, 9 2g
Ejemplo 2.27 Considere un estanque de 10.000 metros cubicos de agua. En el instante t = 0
el agua esta limpia, pero cuenta con una entrada de lquido A y una salida C. Por la entrada A
fluyen 500 metros cubicos al da; en tanto, por C salen 1.300 metros cubicos al da. Considere
ademas, que en t = 0, por A comienza a entrar agua salada con una concentracion de 5 kilos por
cada 1.000 metros cubicos, que se mezcla dentro del estanque de forma homogenea. Ademas,
por la contaminacion ambiente, ingresan deshechos que provocan la disminucion del volumen
del estanque a razon de 50 metros cubicos por da.
dS
Sea S(t) la cantidad de sal dentro del estanque en el instante t. Entonces, dt
es la diferencia
instantanea entre la sal entrante por A y la saliente por C. Para formular una ecuacion dife-
rencial que determine a S(t) debemos considerar dos factores. En primer lugar, la sal entrante
al estanque, es el producto entre la concentracion de sal entrante con la velocidad de entrada,
es decir
5 5
sal entrante = 500 = kilos por da.
1,000 2
Por otro lado, la sal que sale del estanque es el producto de su velocidad de salida con la
concentracion de sal en el instante t, esto es
S(t) 26
sal saliente = 1300 = S(t) kilos por da,
10,000 50t 200 t
donde V (t) = 10,00050t es el volumen del estanque en el instante t. El problema tendra sentido
entonces para 0 t < 200. Luego, la ecuacion diferencial que gobierna la cantidad de sal dentro
del estanque es
dS
= sal entrante sal saliente
dt
5 26
= S(t)
2 200 t
dS 26 5
+ S(t) = .
dt 200 t 2
23
La ecuacion homogenea asociada es

dS 26
+ S(t) = 0,
dt 200 t
con solucion Sh (t) = c(200 t)26 . La solucion particular esta dada por Sp (t) = k(t)Sh (t), donde
k(t) satisface: Z
5 dt 1 1
k(t) = 26
= .
2 (200 t) 10 (200 t)25
Luego,
200 t
Sp (t) = .
10
La solucion general es
200 t
S(t) = c(200 t)26 +
.
10
Para determinar c usamos el dato inicial: en t = 0 no hay sal en el estanque, luego S(0) = 0.
As
20
0 = 20026 c + 20 c = .
20026
Con lo cual, la cantidad de sal en el estanque S Ht L
20
esta dada por:

15
26 10
200 t 200 t
S(t) = 20 + ,
200 10 5
50 100 150 200 t

Figura 2.13. Concentracion de sal en el estanque. Cuando
cuya grafica se presenta en la figura 2.27. t 200, la solucion tiende asintoticamente a la recta 20 t/10.
1 + cosx y 2y
1. dy 1 + cosx dx = 0 4. y + = x3
sen2 x senx x
di
2. L + Ri = Esen(wt), con i(0) = 0 5. y 2 dx (2xy + 3)dy = 0
dt
y
3. y =x 6. xy + y ex = 0, con y(a) = b
x
24
y p
7. y = 1 + x, con y(0) = 0 9. (1 + y 2 )dx = 1 + y 2 seny xy dy
1 x2
8. y y tgx = secx, con y(0) = 0
25
2.5. ECUACIONES QUE PUEDEN REDUCIRSE A LA FORMA LINEAL
2.5. Ecuaciones que pueden reducirse a la forma lineal
Hay ciertas E.D. no lineales de primer orden que pueden reducirse a E.D.L. mediante una
adecuada sustitucion.
2.5.1. Ecuacion de Bernoulli
Una E.D. no L. de Bernoulli es de la forma
y + p(x) y = q(x) y n , (23)
donde n 6= 0, conduce a una E.D.L y n 6= 1 a una E.D. de variables separables.

Si la E.D. 23 se divide por y n queda
dy
y n + p(x) y 1n = q(x)
dx
Si se efectua la sustitucion
dv dy 1 dv dy
v = y 1n = (1 n)y n = y n
dx dx (1 n) dx dx
quedando
1 dv
+ p(x) v = q(x)
(1 n) dx
que es una E.D.L. en la variable v
y y + x y2 = x
Ordenandola
Z
x 2 1(1) dv R
2x dx
R
2x dx
y + xy = v = y = y + 2xv = 2x v(x) = e C + 2x e dx
y dx
Z
x2 x2 2 2 2
v(x) = e C + 2x e dx = ex C + ex = 1 + Cex
por lo tanto
2
y 2 = 1 + Cex
26

5
x y y = x3 y 3
2
Ordenandola
y 5 y 2 5 dv 2
y = x2 y 3 y 3 y = x2 v = y 2 = y 1(3) + v = 5x2
x 2 x 2 dx x
Z Z
R
2/x dx 2
R
2/x dx 2 lnx 2 2 lnx 1
v(x) = e C + 5x e dx = e C + 5x e dx = 2 C + x5
x
por lo tanto
y 2 = Cx2 + x3

x y 4y = x2 y
Ordenandola

4y 2 1 y dv 2 x
y =x y y 4 = x v = y 1/2 = y 1(1/2) v =
x y x dx x 2
Z Z
R
2/x dx x R 2/x dx 2 lnx x 2 lnx 2 1
v(x) = e C+ e dx = e C+ e dx = x C + ln|x|
2 2 2
por lo tanto 2
2 1 4 1
y = x C + ln|x| y = x C + ln|x|
2 2
1. x y + y = x2 y 2 3. 2y dx + (2x y x3 )dy = 0

2. 2x y y y 2 + x = 0 4. 3x dy = y 1 + x senx 3y 3 senx dx
2.5.3. Ecuacion de Ricatti
Una E.D. NO L. de Ricatti es de la forma
y + a2 (x) y 2 + a1 (x)y + a0 (x) = 0 (24)
donde ai (x) son funciones continuas en el intervalo I i = 0, 1, 2 y a2 (x) 6= 0 x I. Para

1
convertirla en una E.D.L. de primer orden se utiliza la sustitucion y(x) = y1 (x) + donde
v
27
y1 (x) es una solucion particular de la ecuacion de Ricatti dada en 24, (esta solucion particular
y1 (x) se encuentra por inspeccion).
Demostracion
Dada
1 dy dy1 1 dv
y(x) = y1 (x) + = 2
v dx dx v dx
reemplazando en 24 se obtiene
2
dy1 1 dv 1 1
2 + a2 (x) y1 (x) + + a1 (x) y1 (x) + + a0 (x) = 0
dx v dx v v
dy1 1 dv 2 1 a1 (x)
2 + a2 (x)y12 (x) + a2 (x) y1 (x) + a2 (x) 2 + a1 (x)y1 (x) + + a0 (x) = 0 (25)
dx v dx v v v
pero y1 (x) es solucion particular de la ecuacion de Ricatti 24, o sea, que
y1 (x) + a2 (x)y12 (x) + a1 (x)y1 (x) + a0 (x) = 0
por lo tanto, en 25
dy1 1 dv 2 1 a1 (x)
2 + a2 (x)y12 (x) + a2 (x) y1 (x) + a2 (x) 2 + a1 (x)y1 (x) + + a0 (x) = 0
dx v dx v v v
queda
1 dv 2 1 a1 (x)
+ a 2 (x) y 1 (x) + a2 (x) + =0
v 2 dx v v2 v
que multiplicandola por v 2 , se obtiene
dv
2a2 (x) v y1 (x) a2 (x) a1 (x) v = 0
dx
ordenandola queda
dv
[2a2 (x) y1 (x) + a1 (x)] v = a2 (x) (26)
dx
obteniendose una E.D.L. en la variable v
dv
+ [p(x)]v = q(x)
dx
28

y = x + 1 xy 2 + 2x2 y x3
si se conoce que y1 (x) = 1 + x, ordenandola
y + xy 2 2x2 y x 1 + x3 = 0
1
aplicando la sustitucion y(x) = 1 + x + , se obtiene aplicando 26
v
dv dv
2x(1 + x) 2x2 v = x 2xv = x
dx dx
Z Z
x2 x2
R R
2x dx 2x dx
v(x) = e C + xe dx = e C + xe dx

x2 1 2 2 1
v(x) = e C ex = Cex
2 2
luego obtenemos que la solucion es
2
y(x) = 1 + x +
Cex2 1
1
y sen2 (x) y 2 + y + cos2 x = 0
senx cosx
si se conoce que y1 (x) = cotx
1
aplicando la sustitucion y(x) = cotx + , se obtiene aplicando 26
v

dv 2 1 2 dv 1
2cotx sen x + v = sen (x) 2senx cosx + v = sen2 (x)
dx senx cosx dx senx cosx

dv 2
sen(2x) + v = sen2 (x)
dx sen(2x)
R
Z R

[sen(2x)+2cosec(2x)] dx 2 [sen(2x)+2cosec(2x)] dx
v(x) = e C sen (x)e dx
como Z
cos(2x)
[sen(2x) + 2cosec(2x)] dx = ln|cosec(2x) + ctg(2x)|
2
y Z Z
cos(2x) cos(2x)
2 +ln|cosec(2x)+ctg(2x)|
sen (x)e 2 dx = sen2 (x)e 2 (cosec(2x) + ctg(2x))dx
Z Z
cos(2x) cos(2x)
2 2
sen (x)e (cosec(2x) + ctg(2x))dx = sen2 (x)e 2 (ctgx)dx =
29
Z Z Z
2 cos(2x) cos(2x) sen(2x) cos(2x) 1 1
sen (x)e 2 (ctgx)dx = senx cosxe 2 dx = e 2 dx = e 2 cos(2x)
2 2
1
e 2 cos(2x) 1 1 1 1 1
v(x) = C e 2 cos(2x)
= tgxe 2
cos(2x)
C e 2 cos(2x)
cosec(2x) + ctg(2x) 2 2
1
cos(2x) 1
v(x) = C tgx e 2 tgx
2

2
y(x) = ctgx +
1
cos(2x)
C tgx e 2 tgx
Nota tambien puede obtenerse como solucion
2
y(x) = ctgx +
C tgx esen2 x tgx
o tambien
2
y(x) = ctgx +
C tgx ecos2x tgx
Observacion En los siguientes problemas para la busqueda de la solucion particular en las

ecuaciones de Ricatti, se pueden suponer que son polinomios del tipo a0 + a1 x + a2 x2 + a3 x3 .......
y se deben buscar los ai intentando con a0 , a0 +a1 x, a0 +a1 x+a2 x2 +a3 x3 , etc., reemplazando
en la E.D. de Ricatti dada.

2x2 y y 2 x2 = 0
intentamos con y = a vemos que 2x2 0 a2 x2 = 0 vemos que a depende de x (contradiccion

pues se supuso que era constante).
intentamos con y = a+bx vemos que 2x2 b(a+bx)2 x2 = 0 x2 (2bb2 1)2abxa2 = 0
de donde debe cumplirse simultaneamente que
2b b2 1 = 0
b = 1
2ab = 0 y1 = 0 + 1 x = x
a = 0
a2 = 0
30
1
efectuando la sustitucion y(x) = x + , se obtiene aplicando 26
v
1 1 v 1
v (2 2
x)v = 2 v + = 2
2x 2x x 2x
Z Z

R 1
1/x dx 1 R 1/x dx lnx 1 1 lnx
v(x) = e C e dx = e C e dx
2 x2 2 x2
Z
1 1 1 1 1 C lnx
v(x) = C x dx = C lnx =
x 2 x2 x 2 x 2x
2x
y(x) = x +
C ln|x|
1. y + y 2 1 = 0 4. y xy 2 + (2x 1)y = x 1
2. y y 2 + 1 = 0 5. xy y 2 + 2xy + 2x + 4 = 0
3. y + y 2 + 3y + 2 = 0 6. y + y 2 2x2 y + x4 2x = 0
7. y + (1 + ex )y 2 2x2 (1 + ex )y + (1 + ex )x4 2x = 0
1
8. y + y 2 1 = 0, tal que y(0) =
3
2.5.5. Sustituciones varias
Existen muchas E.D.N.L. que mediante una sustitucion adecuada la transforma en una E.D.
conocida.
p
Ejemplo 2.34 Resolver y = 3x + 5y + 9
Veamos el caso mas general y = f (ax + by + c) con a, b, c R constantes.

dv dy dy 1 dv
Sea v = ax + by + c =a+b = a
dx dx dx b dx

1 dv dv dv
a = f (v) = b f (v) + a = dx
b dx dx b f (v) + a
31
3. EXISTENCIA Y UNICIDAD DE LAS SOLUCIONES DE UNA EDL DE PRIMER ORDEN
Z
dv 2 2u du 2 u du
obtenemos = dx haciendo v = u = dx =x+C
5 v+3 5u + 3 5 u + 0, 6
0, 4(u 0, 6 ln|u + 0, 6|) = x + C 0, 4(3x + 5y + 9 0, 6 ln|3x + 5y + 9 + 0, 6|) = x + C

y + p(x)y = q(x)ylny
dv 1 dy
Sea v = lny = y = yv reemplazando
dx y dx
yv + p(x)y = q(x)yv v + p(x) = q(x)v v q(x)v = p(x)
que es una E.D.L. en la variable v.
Utilice lo obtenido anteriormente para resolver el siguiente problema de Cauchy con x > 0.

xy 4x2 y + 2y ln(y) = 0
y(1) = 1
Solucion: Usando el cambio de variables v = ln(y), la EDO
xy 4x2 y + 2y ln(y) = 0
se transforma en
2
v + v = 4x
x
la cual es lineal. La solucion es x2 v = x4 + c con c R, equivalentemente
x2 ln(y) = x4 + c
. Como y(1) = 1, entonces c = 1 y as la solucion del problema de Cauchy es
x2 ln(y) = x4 1
3. Existencia y unicidad de las soluciones de una EDL de primer

orden
Problemas de valor inicial

Se puede asegurar que cada EDL de primer orden en un intervalo I, tiene soluciones. De hecho
tiene una infinidad de soluciones, una para cada valor de C en la expresion
R
Z R

p(x)dx p(x)dx
y(x) = e C + q(x) e dx (27)
32
4. ECUACIONES DIFERENCIALES LINEALES DE ORDEN SUPERIOR
y la solucion general de una tales ecuaciones es, por lo tanto, una familia de curvas paralelas
parametricas (cuyo parametro es C) en el plano xy determinada por el intervalo I. Es facil
ver que hay una curva solucion que pasa por un punto cualquiera (x0 , y0), y que 27 puede
resolverse para determinar C cuando x = x0 e y = y0 . El problema de encontrar una funcion
y = f (x) que sea una solucion de una EDL normal de primer orden y + p(x) y = q(x), y que
tambien satisfaga la condicion inicial y(x0 ) = y0 es llamado un problema de valor inicial para
la ecuacion dada.
Teorema 3.1 Si p(x), q(x) C(I), y sea x0 I, entonces la EDL normal de primer orden
y + p(x) y = q(x),
tiene una y solamente una solucion que satisface la condicion inicial y(x0 ) = y0 .
Por lo tanto,todo problema de valor inicial referente a una EDL normal de primer orden, tiene
precisamente una solucion.
4. Ecuaciones Diferenciales Lineales de Orden Superior
Teorema 4.1 Sean

y1 (x), y2 (x), y3 (x), ........yn (x) C n1 (I).
Suponiendo que x0 I de modo que el conjunto

n on
(n1)
yi (x0 ), yi (x0 ), yi (x0 ), yi (x0 ), ......, yi (x0 )
i=1
es linealmente independiente(L.I.) en Rn las funciones {yj (x)}j=1,2,3,....,n son linealmente

independientes.
Ejemplo 4.1 {ex , xex , x2 ex } es L.I. en R3 ya que:

y1 (x) = ex (ex , ex , ex )x=0 = (1, 1, 1)

y2 (x) = xex (xex , ex (1 + x), ex (2 + x))x=0 = (0, 1, 2)

y3 (x) = x2 ex (x2 ex , xex (2 + x), ex (2 + 4x + x2 )x=0 = (0, 0, 2)
vemos claramente que el conjunto {(1, 1, 1), (0, 1, 2), (0, 0, 2)} es L.I.
33
Ejemplo 4.2 {1, x, x2 } es L.I. en R3 ya que:

y1 (x) = 1 (1, 0, 0)x=0 = (1, 0, 0)

y2 (x) = x (x, 1, 0)x=0 = (0, 1, 0)

y3 (x) = x2 (x2 , 2x, 2)x=0 = (0, 0, 2)
Ejemplo 4.3 {ex , x, cosx} es L.I. en R3 ya que:

y1 (x) = ex (ex , ex , ex )x=0 = (1, 1, 1)

y2 (x) = x (x, 1, 0)x=0 = (0, 1, 0)

y3 (x) = cosx (cosx, senx, cosx)x=0 = (1, 0, 1)
Definicion 4.1 Sean
y1 (x), y2 (x), y3 (x), ........yn (x) C n1 (I). x I
el determinante
y1 (x) y2 (x) y3 (x) ... ... yn (x)
y1 (x) y2 (x) y3 (x) ... ... yn (x)
y1 (x) y2(x) y3 (x) ... ... yn (x)
W [y1 (x), y2 (x), y3 (x), ........yn (x)] = (28)
... ... ... ... ... ...
... ... ... ... ... ...
(n1) (n1) (n1) (n1)
y1 (x) y2 (x) y3 (x) ... ... yn (x)
se llama el Wronskiano de las funciones
y1 , y2 , y3 , ........yn
Ejemplo 4.4
ex xex x2 ex
W [ex , xex , x2 ex ] = xex (1 + x)ex (2 + x)ex = 2e3x
x2 ex (x2 + 2x)ex (x2 + 4x + 2)ex
34
Ejemplo 4.5
1 x x2
W [1, x, x2 ] = 0 1 2x = 2
0 0 2
Ejemplo 4.6
ex x cosx
x
W [e , x, cosx] = e x
1 senx = e3x (xcosx xsenx 2cosx)
ex 0 cosx
Ejemplo 4.7
1 x cosx
W [1, x, cosx] = 0 1 senx = cosx
0 0 cosx
Ejemplo 4.8
1 senx cosx
W [1, senx, cosx] = 0 cosx senx = 1
0 senx cosx
1. W [x, senx, cosx] = 5. W [x, senx, cos(2x)] =
2. W [senx, sen(2x), sen(3x)] = 6. W [x + 1, 2 + senx, 4 + cosx] =
3. W [senx, sen2 x, sen3 x] = 7. W [2x + 3, x + senx, x + cosx] =
4. W [x, sen(2x), cosx] = 8. W [1, x, x2 , x3 ] =
Observacion Si
W [y1 , y2, y3 , ........yn ] 6= 0 para algun x y1 , y2 , y3, ........yn
son funcionalmente independientes.
35
4.1. SOLUCION GENERAL DE UNA EDL HOMGENEA DE ORDEN N
4.1. Solucion general de una EDL homgenea de orden n
Aunque no hay metodos explcitos para resolver ED arbitrarias, existen metodos sistematicos
para ciertas ED., aqu veremos ED para la cual hay soluciones unicas.
Sean
y1 , y2 , y3 , ........yn soluciones L.I. de la E.D.L homogenea de orden n
d(n) y d(n1) y d2 y dy
an (x) n
+ an1 (x) n1
+ ...... + a2 (x) 2
+ a1 (x) + a0 (x)y = 0 (29)
dx dx dx dx
entonces la solucion general (o solucion homogenea) de la ecuacion es:
yh (x) = C1 y1 (x) + C2 y2 (x) + C3 y3 (x) + ........ + Cn yn (x) (30)
donde los Ci son constantes que dependen de las n condiciones iniciales (solucion propia o
intrnseca del sistema) y que se determinan solo en el caso de que h(x) = 0 en (31) que
corresponde a (29).
4.2. Solucion general de una EDL no homgenea de orden n
Sea yp una solucion particular(solucion forzante del sistema) de la E.D.L no homogenea de

orden n
an (x) n
+ an1 (x) n1
+ ...... + a2 (x) 2
+ a1 (x) + a0 (x)y = h(x) (31)
dx dx dx dx
y sea yh (x) dada por (30) que es la solucion de la ecuacion homogenea dada en (29), entonces
la solucion general de la ecuacion es:
yG (x) = yh (x) + yp (x), (32)
o sea,
yG (x) = C1 y1 (x) + C2 y2 (x) + C3 y3 (x) + ........ + Cn yn (x) +yp (x) (33)
| {z }
yh (x)
donde los Ci son constantes que dependen de las n condiciones iniciales que se determinan
despues de haber encontrado la solucion particular yp .
Nota Debemos generalizar el teorema 3.1 que era para EDL normales de primer orden a EDL
normales de orden n.
36
4.3. OPERADORES DIFERENCIALES
Teorema 4.2 Teorema de la existencia y unicidad para EDL normales de orden n
Sean a0 (x), a1 (x), a2 (x), ......, an1 (x), an (x), h(x) C(I), y sea x0 I. Entonces , si
y0 , y1 , y2 , ......., yn1 son n numeros arbitrarios, existe una y solamente una solucion y(x) que
satisface la EDL normal de orden n
an (x) n
+ an1 (x) n1
+ ...... + a2 (x) 2
+ a1 (x) + a0 (x)y = h(x)
dx dx dx dx
con la propiedad de que
y(x0 ) = y0 , y (x0 ) = y1 , y (x0 ) = y2 , ......, y (n1) (x0 ) = yn1.
Como en el caso de ecuaciones de primer orden, el problema de encontrar una solucion de

la EDL anterior que satisfaga la n condiciones dadas es llamado un problema de valor inicial
con las condiciones
y(x0 ) = y0 , y (x0 ) = y1 , y (x0 ) = y2 , ......, y (n1) (x0 ) = yn1.
4.3. Operadores Diferenciales
Sea D el operador que denota diferenciacion con respecto a una variable (por ejemplo x),
D el operador que indica una doble diferenciacion, D3 el operador que indica una triple dife-
2
renciacion, etc.
dy dk y
Ejemplo 4.9 El operador D aplicado a y es Dy , generalizando Dk y k indicando la
dx dx
diferenciacion kesima con respecto a la variable x.
Pueden formarse polinomios en D, los que tienen expresiones de la forma
L(D) = an Dn + an1 Dn1 + ...... + a2 D2 + a1 D + a0 Id ,
donde los a0 , a1 , a2 , ...., an son numeros reales, y Id corresponde a la funcion identidad.(Tales

expresiones son conocidas como operadores diferenciales lineales de orden n de coeficientes
constantes. Este operador aplicado a y queda

L(D)(y) = an Dn + an1 Dn1 + ...... + a2 D2 + a1 D + a0 I (y) =
L(D)(y) = an y (n) + an1 y (n1) + ...... + a2 y + a1 y + a0 y
37
4.3. OPERADORES DIFERENCIALES
Ejemplo 4.10
d2 y dy
(D2 + D 2I)y = (D2 + D 2)y = + 2y
dx2 dx
se puede ver que el polinomio (D2 + D 2) = (D + 2)(D 1) = (D 1)(D + 2) y en este caso
(coeficientes constantes) es facil verificar que:

dy
(D + 2)(D 1)y = (D + 2)[(D 1)y] = (D + 2) y =
dx

d dy dy d2 y dy dy
(D + 2)(D 1)y = y +2 y = 2 + 2 2y
dx dx dx dx dx dx
d2 y dy
(D + 2)(D 1)y = + 2y
dx2 dx
por otro lado

dy
(D 1)(D + 2)y = (D 1)[(D + 2)y] = (D 1) + 2y =
dx

d dy dy d2 y dy dy
(D 1)(D + 2)y = + 2y + 2y = 2 + 2 2y
dx dx dx dx dx dx
d2 y dy
(D + 2)(D 1)y = 2 + 2y
dx dx
Lo que se concluye que los operadores diferenciales de coeficientes constantes son conmuta-
tivos, lo que no sucede con los operadores de coeficientes variables.
Se puede verificar que:

Dk (y1 (x) + y2 (x)) = Dk y1 (x) + Dk y2 (x) , R
L(D)(y1 (x) + y2 (x)) = L(D)(y1 (x)) + L(D)(y2 (x)), R
38
5. EDL NORMALES HOMOGENEAS
5. EDL normales homogeneas
5.1. EDL normal homogeneas de segundo orden
Sean a0 (x), a1 (x), a2 (x) C(I), y sea a2 (x) 6= 0 con x I. Una EDL normal de segundo
orden homogenea es de la forma
d2 y dy
a2 (x) 2
+ a1 (x) + a0 (x)y = 0
dx dx
dado que a2 (x) 6= 0 la EDL anterior se puede reescribir como:
y + p(x)y + q(x)y = 0 (34)
donde
a1 (x) a0 (x)
p(x) = , q(x) =
a2 (x) a2 (x)
5.2. Teorema de Abel
El teorema de Abel propone encontrar una segunda solucion y2 (x) de una EDL homogenea
de segundo orden si se conoce previamente una primera solucion y1 (x), para ello propone que
la segunda solucion es un multiplo funcional de la primera y2 (x) = K(x) y1 (x)
Teorema 5.1 Sean p(x), q(x) C(I) Si y1 (x) es una solucion no trivial de la EDL
homogenea
y + p(x)y + q(x)y = 0
entonces y2 (x) = K(x) y1 (x) donde

Z R
e p(x)dx
K(x) = dx
y12 (x)
Teorema 5.2 Sean y1 (x) e y2 (x) dos soluciones no triviales de la EDL (34). Entonces cualquier
combinacion lineal de ellas: y(x) = C1 y1 (x) + C2 y2 (x) es tambien una solucion de (34).
Volviendo al teorema de Abel si y2 (x) es una solucion L.I. de y1 (x) de la EDL (34), entonces
la solucion general de (34) es
yG (x) = yh (x) = C1 y1 (x) + C2 y2 (x) = y1 (x) [C1 + C2 K(x)]
39
5.2. TEOREMA DE ABEL
Ejemplo 5.1 Dada la EDL de segundo orden
(1 2x x2 )y + 2(1 + x)y 2y = 0
a) Encontrar la segunda solucion de la EDL dada, si se sabe que y1 (x) = x + 1
b) Determine la solucion homogenea de la EDL.
Solucion: a) Utilizando Abel

Z R 2(1 + x)
Z
dx 2 Z
e 1 2x x2 eln(x + 2x 1) (x2 + 2x 1)
K(x) = dx = dx = dx =
(1 + x)2 (1 + x)2 (1 + x)2
2
=x+ .
1+x
Tenemos que y2 (x) = 2 + x(x + 1) .
b) La solucion general es de la forma yh = C1 (x + 1) + C2 (2 + x(x + 1))
Ejemplo 5.2 Dada la EDL de segundo orden
(1 + 4x x2 )y + 2(x 2)y 2y = 0
a) Encontrar la segunda solucion de la EDL dada, si se sabe que y1 (x) = 2 x

Z R 2(x 2)
Z
dx 2 Z 2
e 1 + 4x x2 eln(x 4x 1) x 4x 1
K(x) = dx = dx = dx
(2 x)2 (2 x)2 (x 2)2
Z
5 5
K(x) = 1 dx = x +
(x 2)2 x2
Tenemos que y2 (x) = 2x x2 5
b) La solucion general es de la forma yh = C1 (2 x) + C2 (2x x2 5).
2
Ejemplo 5.3 Dada la EDL de segundo orden y + y + y = 0
x
40
5.2. TEOREMA DE ABEL
sen(x)
a) Encontrar la segunda solucion de la EDL dada, si se sabe que y1 (x) =
x

Z R 2
Z
dx Z
e x e2ln(x) 2
K(x) = 2 dx = 2 dx = cosec (x) dx = cotan(x)
sen(x) sen(x)
x x
cos(x)
Tenemos que y2 (x) =
x
sen(x) cos(x)
b) La solucion general es de la forma yh = C1 + C2
x x
En los problemas siguientes se da la EDL y una solucion y1 (x)
1. y 4y + 4y = 0, donde y1 (x) = e2x
2. y 2ay + a2 y = 0, donde y1 (x) = eax
3. 3xy y = 0, donde y1 (x) = 4
4. (1 x2 )y 2xy = 0, donde y1 (x) = 7
5. (1 x2 )y 2xy + 2y = 0, donde y1 (x) = x
6. 2xy ex y = 0, donde y1 (x) = 7
7. x2 y 3xy + 4y = 0, donde y1 (x) = x2
8. y + tgx y 6ctgx y = 0, donde y1 (x) = sen3 x

2x 2
9. y + p(x)y + y = 0, donde y1 (x) = ex
x2 2x
10. xy y + 4x3 y = 0, donde y1 (x) = sen(x2 )
11. y 2y + 2y = 0, donde y1 (x) = ex cosx
41
6. OPERADOR DIFERENCIAL LINEAL DE ORDEN N DE COEFICIENTES CONSTANTES
12. x2 y 3xy + 4y = 0, donde y1 (x) = x2 lnx
13. y + (tgx 2ctgx)y = 0, donde y1 (x) = sen3 x

2 x2
14. y + y + q(x)y = 0, donde y1 (x) = ex
x2 2x

2x 2
15. y + p(x)y + 2
y = 0, donde y1 (x) = x2
x 2x

2 x2
16. y + y + q(x)y = 0, donde y1 (x) = x2
x2 2x
17. x2 y + xy + 4y = 0, donde y1 (x) = sen(2lnx)
6. Operador diferencial lineal de orden n de coeficientes constantes
Sea
L(D) = an Dn + an1 Dn1 + ...... + a2 D2 + a1 D + a0 I,
donde los a0 , a1 , a2 , ...., an son numeros reales. Este operador aplicado a y queda

L(D)(y) = an Dn + an1 Dn1 + ...... + a2 D2 + a1 D + a0 I (y) =
| {z }
L(D)
L(D)(y) = an y (n) + an1 y (n1) + ...... + a2 y + a1 y + a0 y

d0 y
donde I(y) = D0 (y) = =y correponde a la funcion identidad, que usualmente se omite
dx0
por simplicidad.
Ejercicio Apliquemos L(D) a la funcion emx
L(D)(emx ) = an mn emx + an1 mn1 emx + ...... + a2 memx + a1 y + a0 emx =
L(D)(emx ) = emx (an mn + an1 mn1 + ...... + a2 m + a1 y + a0 ) = emx L(m)

| {z }
L(m)
Si m es una raz de la ecuacion L(m) = 0 (llamada ecuacion caracterstica)

L(m) = 0 L(D)(emx ) = 0, pero una ecuacion L(m) = 0 de grado n en la variable m es
un polinomio que tiene n soluciones reales. Por lo tanto, para cada mi tal que L(mi ) = 0
L(D)(emi x ) = 0; luego, si m1 , m2 , m3 , ......., mi , ......mn son soluciones de L(m) = 0
L(D) (C1 em1 x + C2 em2 x + C3 em3 x + ...... + Cn emn x ) = 0
42
L(D)(y) = 0 yh (x) = C1 em1 x + C2 em2 x + C3 em3 x + ...... + Cn emn x
Los valores de las constantes Ci R se determinan mediante evaluacion con las condiciones
iniciales, cuando se trate de un problema de valor inicial.
Ejemplo 6.1 Resolver el siguiente problema de valor inicial

2y + 5y 12y = 0

y(0) = 2

y (0) = 3
Solucion:
3
(2D2 + 5D 12)(y) = 0 L(m) = 2m2 + 5m 12 = 0 m1 = 4 m2 =
2
3
yh (x) = C1 em1 x + C2 em2 x = C1 e4x + C2 e 2 x
con las c.i.
( (
y(0) = 2 C1 + C2 = 2 C1 = 0 3
yh (x) = 2e 2 x
y (0) = 3 4C1 + 23 C2 = 3 C2 = 2
Nota
Consideremos la ED de segundo orden de coeficientes constantes general:
ay + by + cy = 0
de lo cual obtenemos la ecuacion caracterstica:
am2 + bm + c = 0
Las soluciones pueden ser:
a) Dos races reales y distintas m1 y m2 , yh (x) = C1 em1 x + C2 em2 x como en el ejemplo

6.1.
b
b) Dos races reales e iguales m1 = m2 = m = , por lo tanto, se tiene solo una
2a
bx

solucion y1 (x) = emx = e 2a y la segunda solucion se obtiene por el teorema de Abel
bx

bx Z a bx
e
y2 (x) = e 2a dx = xe 2a
2bx

e 2a
43
luego:
bx bx

yh (x) = C1 e 2a + C2 xe 2a
Observacion
Al factor x que se obtiene por el teorema de Abel se le conoce como factor de desacoplo,
el cual se aplica cada vez que sucedan races repetidas, como se vera en el ejemplo 6.2
c) Dos races complejas conjugadas donde m1 y m2 son de la forma:
m1 = + i
m2 = i
usando las relaciones de Euler
ei = cos + i sin
ei = cos i sin
se obtiene
yh (x) = ex (C1 cos(x) + C2 sin(x))

(D 3)4 (y) = 0 yh (x) = e3x C1 + C2 x + C3 x2 + C4 x3
y en general

(D a)n (y) = 0 yh (x) = eax C1 + C2 x + C3 x2 + .... + Cn xn1
y para a = 0
Dn (y) = 0 yh (x) = C1 + C2 x + C3 x2 + .... + Cn xn1
((D 1)2 + 4)(y) = 0 m = 1 2i yh (x) = ex (C1 sen(2x) + C2 cos(2x))
y en general
((D a)2 + b2 )(y) = 0 yh (x) = eax (C1 sen(bx) + C2 cos(bx))
y para a = 0
(D2 + b2 )(y) = 0 yh (x) = C1 sen(bx) + C2 cos(bx)
44
7. OPERADOR ANULADOR
(D4 +2D3 +D2 )(y) = 0 m1 = m2 = 1, m3 = m4 = 0 yh (x) = ex (C1 + C2 x)+(C3 + C4 x)
(D2 + 32 )(D 2)(D + 7)(y) = 0 yh (x) = C1 sen(3x) + C2 cos(3x) + C3 e2x + C4 e7x
((D 1)2 + 4)2 (y) = 0 yh (x) = ex (C1 sen(2x) + C2 cos(2x)) + ex (C1 xsen(2x) + C2 xcos(2x))
y en general
((D a)2 + b2 )n (y) = 0

yh (x) = eax sen(bx) C1 + C2 x + ... + Cn xn1 + eax cos(bx) Cn+1 + Cn+2 x + ... + C2n xn1
y para a = 0
(D2 + b2 )n (y) = 0

yh (x) = sen(bx) C1 + C2 x + ... + Cn xn1 + cos(bx) Cn+1 + Cn+2 x + ... + C2n xn1
7. Operador Anulador
Definicion 7.1 Si
an Dn + an1 Dn1 + + a1 D + a0
es tal que:

an Dn + an1 Dn1 + + a1 D + a0 (f (x)) = 0
se dice que
an Dn + an1 Dn1 + + a1 D + a0
anula a f.
Ejemplo 7.1 Dn anula a 1, x, x2 , x3 , . . . , xn1
Ejemplo 7.2 (D a)n anula a eax , xeax , x2 eax , . . . , xn1 eax
Ejemplo 7.3 Obtener un operador que anule a:
e3x + xex
45
Solucion:
(D + 3)(e3x ) = 0
(D 1)2 (xex ) = 0
luego el operador es (D + 3)(D 1)2 pues
(D + 3)(D 1)2 (e3x + xex ) = 0
Ejemplo 7.4 Si y son numeros reales, de modo que la ecuacion (m2 2m+(2 + 2 ))n = 0,
tiene races complejas: + i, y i, ambas de multiplicidad n, luego:
n
D2 2D + (2 + 2 )
anula a las funciones:
ex cos(x), xex cos(x), x2 ex cos(x), . . . , xn1 ex cos(x),
ex sin(x), xex sin(x), x2 ex sin(x), . . . , xn1 ex sin(x)
Ejemplo 7.5 Obtener un operador diferencial que anule a:
7 x + 6sen(3x)
los anuladores son:

D2 (7 x) = 0
(D2 + 9)(sen(3x)) = 0
D2 (D2 + 9)(7 x + 6sen(3x)) = 0
por lo tanto el operador es D2 (D2 + 9)
46
En los siguientes ejercicios encontrar la solucion general de cada EDL y cuando se den c.i.,
de la solucion particular.
1. y 4y = 0 7. x + 5x + 6x = 0, x(0) = 1, x (0) = 2
2. x 3x + 2x = 0 8. z z 6z = 0, z(0) = 1, z (0) = 1
3. y + 6y + 9y = 0 9. x x x + x = 0
4. y 5y = 0 10. y (4) 9y = 0
5. y + 2y + 2 y = 0 11. z (5) 2z + y = 0
6. x + x 6x = 0, x(0) = 0, x (0) = 5 12. y (4) 5y + 4y = 0
13. z 3z + 3z z = 0, z(0) = 1, z (0) = 2, z (0) = 3
14. z (4) 4z = 0, z(0) = 1, z (0) = 3, z (0) = 0, z (0) = 16
15. z (4) 10z + 9z = 0, z(0) = 3, z (0) = 9, z (0) = 11, z (0) = 81
En los siguientes ejercicios encontrar el operador anulador L(D)(f (x)) = 0
1. f (x) = senx cos(2x) 7. f (x) = x2 senx + xe4x
2. f (x) = xsenx + 4cos(3x) 8. f (x) = sen2 x + cos2 (2x)
3. f (x) = 4x 6cos2 (2x) 9. f (x) = e2x + e2x
4. f (x) = xsen2 x + 5cos(2x) 10. f (x) = senhx + cosh(2x)
5. f (x) = ex senx + ex cos(2x) 11. f (x) = 1 + x + x2 + senx + cos(2x)
6. f (x) = xex senx + x2 12. f (x) = e2x senx + xcos(2x)
47
8. EDL DE COEFICIENTES CONSTANTES NO HOMOGENEOS
8. EDL de coeficientes constantes no homogeneos
Una EDL de CC no homogenea, puede escribirse como:
L(D)(y) = h(x) (35)
Sea yp una solucion particular(solucion forzante del sistema) de la E.D.L no homogenea de

orden n
an n
+ an1 n1
+ ...... + a2 2
+ a1 + a0 y = h(x) (36)
dx dx dx dx
y sea yh (x) la solucion de la EDL de CC homogenea dada en la seccion anterior, entonces la
solucion general de la ecuacion es:
yG (x) = yh (x) + yp (x), (37)
o sea,
yG (x) = C1 y1 (x) + C2 y2 (x) + C3 y3 (x) + ........ + Cn yn (x) +yp (x) (38)
| {z }
yh (x)
donde los Ci son constantes que dependen de las n condiciones iniciales que se determinan
despues de haber encontrado la solucion particular yp .
Ahora debemos encontrar yp existen varios metodos para determinar yp , los dos mas conocidos
son:
1. El Metodo de los Coeficientes Indeterminados (MCI) y
2. El Metodo de Variacion de Parametros (MVP).
8.1. Metodo de los Coeficientes Indeterminados
El MCI tiene la ventaja de ser un metodo sencillo, pero para poder aplicarlo es necesario
conocer la forma que debera tener la solucion particular.
Es posible encontrar una solucion si

una constante k

un polinomio en x

h(x) = una funcion ex

sen(x), cos(x)

suma y/o producto de las funciones anteriores
48
8.1. METODO DE LOS COEFICIENTES INDETERMINADOS
Sea L(y) = h(x) y sea L2 el operador anulador de h(x) de modo que
L2 (h(x)) = 0 L2 L(y) = L2 (h(x)) = 0 L (y) = 0 donde L = L2 L
Sea y (x) = yh + yp tal que L (y (x)) = 0 yp = y yh
Se sustituye yp en la ED L(yp ) = h(x) y de este modo se pueden determinar las constantes que
aparecen en yp .

y + 4y = 7
Solucion:
(D2 + 4)(y) = 0 yh = C1 sen(2x) + C2 cos(2x)
L2 (7) = D(7) = 0 L2 = D L = D(D2 + 4)
y = C1 sen(2x) + C2 cos(2x) + C3 yp = C3
7
yp + 4yp = 0 + 4C3 7 C3 =
4
7
yG (x) = C1 sen(2x) + C2 cos(2x) +
4
Ahora se procede a calcular las constantes Ci mediante las c.i.
Ejemplo 8.2
Resolver
y + 8y = 5x + ex
Solucion:
D2 (x) = 0,
(D + 1)ex = 0
L = D2 (D + 1)(D2 + 8) = 0
La solucion es de la forma:

y (x) = C1 cos(2 2x) + C2 sen(2 2x) + C3 + C4 x + C5 ex
| {z } | {z }
yh yp
Como yp = A + Bx + Cex , sustituyendo en la ecuacion:
yp + 8yp = 5x + ex
49
8.1. METODO DE LOS COEFICIENTES INDETERMINADOS
luego:
Cex + 8A + 8Bx + 8Cex = 5x + ex
as:
A = 0
5
B = 8
1
C = 9
por lo tanto, la solucion es:

5 1
yG (x) = C1 cos(2 2x) + C2 sen(2 2x) + x + ex
| {z } |8 {z9 }
yh
yp

(D2 + 2D + 1)(y) = 8ex + 25cos(2x)
Solucion:
yh (x) = C1 ex + C2 xex
L2 = (D 1)(D2 + 4) yp = Aex + Bsen(2x) + Ccos(2x)
(D2 + 2D + 1)(yp ) = (D2 + 4 + 2D 3)(yp ) = (D2 + 4)(Aex ) + (2D 3)(yp ) =
(D2 1 + 5)(Aex ) + (2(D 1) 1)(yp ) =
= 5Aex Aex + 4Bcos(2x) 4Csen(2x) 3Bsen(2x) 3Ccos(2x) =
= 4Aex + (4B 3C)cos(2x) (3B + 4C)sen(2x) 2ex + 25cos(2x)
as:
4A = 8 A = 2
4B 3C = 25 B = 4
3B + 4C = 0 C = 3
por lo tanto, la solucion particular es:
yp = 2ex + 4sen(2x) 3cos(2x)
y la solucion general es
yG (x) = C1 ex + C2 xex + 2ex + 4sen(2x) 3cos(2x)

| {z } | {z }
yh yp
50
8.2. METODO DE VARIACION DE PARAMETROS
Empleando el MCI, encontrar yp en:
1. y + 6y + 10y = x4 4x2 + 2 8. y + 4y + 4y = xex 4senx
2. 6y + 2y y = 7xex (1 + x) 9. y (4) 2y + y = ex
3. y + y 2y = 3ex senx 10. y (4) 2y + y = x3
4. y 4y + 8y = e2x (1 + sen(2x) 11. y (4) 2y + y = ex + x3
5. y (4) 8y + 16y = xsenhx 12. y (4) + y = cos(4x)
6. y + 4y = 16xsenh(2x) 13. y (4) 2y + y = sen(2x)
7. y 4y + 5y = 2e2x cosx 14. y (4) 2y + y = ex + x3 + sen(2x)
8.2. Metodo de Variacion de Parametros
Primero consideremos la EDL normal de segundo orden definida en un intervalo I
y + a1 (x)y + a0 (x)y = h(x) (39)
y la solucion homogenea asociada es
yh (x) = C1 y1 (x) + C2 y2 (x)
Se busca una solucion particular yp (x) de (39) tal que yp (x0 ) = 0, yp (x0 ) = 0 donde x0 I
La construccion de yp (x) se empieza con la suposicion de que cualquier solucion particular
de (39) debe estar relacionada con yh (x) y para ello se intenta alterar esta ultima tal que se
convierta en una solucion particular de (39). Una forma de hacerlo es permitir que las constantes
C1 y C2 , sean funciones de x, (se le llama parametros), o sea, K1 (x) y K2 (x), de modo que
yp (x) = K1 (x)y1 (x) + K2 (x)y2 (x) (40)
Si (40) se deriva una, dos veces, simplificando la notacion omitiendo la mencion de la variable
x
yp (x) = K1 y1 + K1 y1 + K2 y2 + K2 y2
yp(x) = K1 y1 + 2K1 y1 + K1 y1 + K2 y2 + 2K2 y2 + K2 y2
51
y se sustituye en (39), donde se obtiene:
K1 (y1 + a1 (x)y1 + a0 (x)y1 ) + K2 (y2 + a1 (x)y2 + a0 (x)y2 )+

| {z } | {z }
=0 =0
+K1 y1 + 2K1 y1 + K2 y2 + 2K2 y2 + a1 K1 y1 + a1 K2 y2 = h(x)
pues y1 e y2 eran soluciones homogeneas,
K1 y1 + 2K1 y1 + K2 y2 + 2K2 y2 + a1 K1 y1 + a1 K2 y2 = h(x)
reagrupando
(K1 y1 + K2 y2 ) + a1 (K1 y1 + K2 y2 ) + (K1 y1 + K2 y2 ) = h(x)
Esta identidad debe mantenerse si (40) ha de ser una solucion particular de (39) y escogiendo
K1 (x) y K2 (x), tal que
(K1 y1 + K2 y2 ) +a1 (K1 y1 + K2 y2 ) + (K1 y1 + K2 y2 ) = h(x)

| {z } | {z } | {z }
=0 =0 =h(x)
K1 y1 + K2 y2 = 0

K1 y1 + K2 y2 = h(x)
Resolviendo por Cramer:
h(x)y2 h(x)y1
K1 = y K2 =
W W
donde W es el Wronskiano de las funciones y1 (x), y2 (x)
y1 y2
W = W [y1(x), y2 (x)] =
y1 y2
Z x
y2 (x)y1 (t) y1 (x)y2 (t)
yp (x) = h(t) dt
x0 W [y1 (t), y2 (t)]
y + y = secx
en este caso no se puede aplicar MCI, pues no existe un operador L tal que L(secx) = 0
Solucion:
yh (x) = C1 senx + C2 cosx
yp (x) = K1 (x)senx + K2 (x)cosx
52
K1 (x)senx + K2 (x)cosx = 0
K1 (x)cosx K2 (x)senx = secx
senx cosx
W = W [y1 (x), y2 (x)] = = 1
cosx senx
K1 = secx cosx = 1 K2 = secx senx = tgx
Z
K1 = x K2 = tgxdx = ln|cosx|
yp = K1 senx + K2 cosx = xsenx + cosx ln|cosx|

yG = C1 senx + C2 cosx + xsenx + cosx ln|cosx|

y y 2y = ex senx
Este ejemplo se resolvera por MVP, aun cuando se puede aplicar MCI, pues existe un operador
L tal que L(ex senx) = 0 L = (D + 1)2 + 1 (se deja como ejercicio comprobarlo)
Solucion:
yh (x) = C1 e2x + C2 ex
yp (x) = K1 (x)e2x + K2 ex
K1 (x)e2x + K2 ex = 0
2K1 (x)e2x K2 ex = ex senx
e2x ex 1 1
W = W [y1 (x), y2 (x)] = 2x x
= e2x ex = 3ex
2e e 2 1
1
K1 = e3x senx K2 = senx
3
1 3x 1
K1 = e (3senx cosx) K2 = cosx
10 3
1 1
yp = K1 e2x + K2 ex = ex (3senx + cosx) + ex cosx
30 3
1 1
yG = C1 e2x + C2 ex ex (3senx + cosx) + ex cosx
30 3
y + y = tgx
en este caso no se puede aplicar MCI, pues no existe un operador L tal que L(tgx) = 0
Solucion:
yh (x) = C1 senx + C2 cosx
53
yp (x) = K1 (x)senx + K2 (x)cosx

K1 (x)senx + K2 (x)cosx = 0
K1 (x)cosx K2 (x)senx = tgx
senx cosx
W = W [y1 (x), y2 (x)] = = 1
cosx senx
K1 = senx K2 = senx tgx
Z Z
sen2 x
K1 = cosx K2 = dx = (secx cosx) dx = senx ln|secx + tgx|
cosx
yp = K1 senx+K2 cosx = cosxsenx+cosxsenxcosxln|secx+tgx| = cosxln|secx+tgx|
yG = C1 senx + C2 cosx + xsenx cosx ln|secx + tgx|
xy + y = x + 1
en este caso no se puede aplicar MCI, pues no es una EDL de CC es de coeficientes variables,
por inspeccion vemos que una solucion homogenea es y1 = 1
Solucion:
Por el teorema de Abel se puede encontrar que y2 = ln|x| x R {0}, luego tenemos que
yh (x) = C1 + C2 ln|x|
yp (x) = K1 (x) + K2 (x)lnx
K1 (x) + K2 (x)lnx = 0
h(x)
K2 (x) x1 = x
= x+1
x
Se divide por x la funcion h(x) pues la ED debe ser normal y el termino x multiplica a y , o
sea, estamos resolviendo y + x1 y = x+1
x
1 lnx 1
W = W [y1(x), y2 (x)] = 1
=
0 x
x
K1 = (1 + x)lnx K2 = 1 + x
2
x2 x x2
K1 = +x + x ln|x| K2 = x +
4 2 2
2
2 2
x x x x2
yp = K1 + K2 ln|x| = +x + x ln|x| + + x ln|x| = +x
4 2 2 4
x2
yG = C1 + C2 ln|x| + +x
4
54
8.3. PRINCIPIO DE SUPERPOSICION
Nota En el ejemplo anterior podramos haber supuesto que la solucion particular fuese
x2
Ax2 + Bx + C y se habra obtenido que yp = 4
+ x + C el cual C podra ser absorbido por C1 .
PERO, ESTE RESULTADO TOTALMENTE CASUAL NO SUCEDE SIEMPRE,
POR LO TANTO, NO UN ES METODO ACEPTABLE, PUES NOS LLEVA A SOLUCIONES
ERRADAS, cuando se esta resolviendo EDL de coeficientes variables.
Encuentre la solucion general de las EDL
1. y y = sec2 x tgx e4x (5x 2)

7. y 3y 4y =
x3
2. y + y = ctgx
e2x

3. y + 4y = secx tgx 8. y 4y + 4y =
(1 + x)
ex
4. y 2y + y = 9. y + 2y + y = ex ln|x|
(1 x)2
5. y y = sen2 x 10. y + 4y = sec(2x)
(2x 1)ex
6. y y = 11. y + 5y + 9y + 5y = 2ex secx
x2
8.3. Principio de superposicion
Teorema 8.1 Sea yg1 una solucion de la EDL
a(x)y + b(x)y + c(x)y = h1 (x) (41)
y sea yg2 una solucion de la EDL
a(x)y + b(x)y + c(x)y = h2 (x), (42)
entonces se puede demostrar que yg = yg1 + yg1 es una solucion de la EDL 43, es decir, la
solucion de 43 se obtiene superponiendo las soluciones de las ecuaciones 41 y 42
a(x)y + b(x)y + c(x)y = h(x) (43)
donde
h(x) = h1 (x) + h2 (x)
55
8.4. METODO DE VP PARA EDL DE ORDEN SUPERIOR
En pocas palabras, el principio de superposicion nos dice que si podemos partir la funcion
h(x) en una suma de dos (o mas) expresiones mas sencillas hk (x), entonces podemos restringir
nuestra atencion a resolver las ecuaciones no homogeneas
a(x)y + b(x)y + c(x)y = hk (x), con k = 1, 2, ...., m,
pues la solucion es simplemente la suma de las soluciones de estas ecuaciones. Esta importante
relacion se denomina Principio de superposicion, el cual se puede extender facilmente a ED
de mayor orden.
Ejemplo 8.8
8.4. Metodo de VP para EDL de Orden Superior
El MVP, puede extenderse facilmente a ED de orden superior. Sea la ED normal definida

en un intervalo I R
y (n) + an1 (x)y (n1) + ........ + a0 (x)y = h(x) (44)
y se supone que la solucion homogenea es de la forma
yh (x) = C1 y1 (x) + C2 y2 (x) + ....... + Cn yn (x) (45)
y se busca una solucion particular de la forma
yp (x) = K1 (x)y1 (x) + K2 (x)y2 (x) + ....... + Kn (x)yn (x) (46)
usando (46) en (44) se obtiene omitiendo los (x)
K1 y1 + K2 y2 + ....... + Kn yn = 0
K1 y1 + K2 y2 + ....... + Kn yn = 0
K1 y1 + K2 y2 + ....... + Kn yn = 0
................ = 0
(n2) (n2) (n2)
K1 y1 + K2 y2 + ....... + Kn yn = 0
(n1) (n1) (n1)
K1 y1 + K2 y2 + ....... + Kn yn = h(x)
resolviendo el sistema por Cramer se obtiene
Vi
Ki (x) = h(x)
W [y1 , y2, ....., yn ]
56
8.4. METODO DE VP PARA EDL DE ORDEN SUPERIOR
donde
0

0

0

Vi =
..

0

1
igual que en el caso de ED de segundo orden se obtiene
Z x
V1 (t)y1 (x) + V2 (t)y2 (x) + ... + Vn (t)yn (x)
yp (x) = h(t) dt
x0 W [y1(t), y2 (t), ....., yn (t)]
donde x0 es cualquier punto en I, y como
Z x
yp (x) = F (x, t) h(t) dt
x0
donde
V1 (t)y1 (x) + V2 (t)y2 (x) + ... + Vn (t)yn (x)
F (x, t) =
W [y1 (t), y2(t), ....., yn (t)]
o usando determinantes
y1 (t) ... yn (t)
y1 (t) ... yn (t)
1
F (x, t) = .. .. ..
W (n2) (n2)
y1 (t) ... yn (t)
y1 (x) ... yn (x)
y el wronskiano W es
y1 (t) ... yn (t)
y1 (t) ... yn (t)
W = .. .. ..
(n2) (n2)
y1 (t) ... yn (t)
(n1) (n1)
y1 (t) ... yn (t)
y + y = tgx
yh (x) = C1 + C2 senx + C3 cosx
yp (x) = K1 (x) + K2 (x)senx + K3 (x)cosx
57
8.5. ECUACION DE EULER
K1 + K2 senx + K3 cosx = 0
K2 cosx K3 senx = 0
K2 senx K3 cosx = tgx
sumando la primera y tercera ecuacion se obtiene
K1 = tgx K1 = ln|cosx|
quedando el sistema
K2 cosx K3 senx = 0
K2 senx K3 cosx = tgx
del cual se obtiene
K2 = senx tgx K2 = ln|secx + tgx| + senx
K3 = cosx tgx K3 = cosx
yp (x) = ln|cosx| senx ln|secx + tgx| + sen2 x + cosx cosx

yp (x) = ln|cosx| senx ln|secx + tgx| + 1
yG (x) = C1 + C2 senx + C3 cosx ln|cosx| senx ln|secx + tgx|
el 1 de la solucion particular lo absorbe C1 de la solucion homogenea.
Ejemplo 8.10
8.5. Ecuacion de Euler
Una EDL homogenea de orden n de forma
xn y (n) + an1 xn1 y (n1) + ..... + a1 xy + a0 y = h(x) (47)
con
a0 , a1 , a2 , .......an1
son constantes reales, donde, el numero del exponente de la variable x es igual al orden de la
derivada a la que multiplica, se conoce como Ecuacion de Euler , ella esta definida en todo
el eje x, pero es normal en intervalos que no contengan el punto x = 0. Si h(x) = 0 entonces
se dice que es una ecuacion de Euler homogenea. La particularidad de esta ecuacion es que
es una de las pocas con coeficientes variables en las que se tiene un metodo analtico capaz
58
de encontrar sus soluciones homogenea y particular. Existen variadas formas de presentar la

solucion, en este caso, usaremos el siguiente:
Considerese la ecuacion de Euler homogenea de orden 2
a2 x2 y + a1 xy + a0 y = 0
Esta ED (47) con h(x) = 0 se puede convertir en una EDL de CC en la variable z al efectuar
el cambio de variable z = lnx, considerando la regla de la cadena, se tiene que
y

z

x
dy dy dz dy 1 1
y (x) = = = = y (z) xy (x) = y (z)
dx dz dx dz x x
As, con el resultado anterior hemos logrado eliminar el coeficiente variable del termino con
primera derivada en la ecuacion de Euler. Usando la notacion de operador lineal, se obtiene que
xD = Dz . En la expresion anterior D denota derivada en variable y(x), y Dz denota derivada
en variable y(z). Ahora, de forma analoga al procedimiento anterior, analizaremos que ocurre
con el termino de la segunda derivada en la ecuacion de Euler. Tenemos que:
dy dy dz
y (x) = =
dx dz dx
dy

d dx d dy dz
dz dx
d dy
dz dz dy d dx dz
y (x) = = = +
dx dx dx dx dz dx
dy dz

d dz dx d2 y dz dz dy d2 z 1 1
y (x) = = 2 + 2 = y (z) 2 y (z) 2
dx dz dx dx dz dx x x
2
x y (x) = y (z) y (z)
Observar que usando la notacion de operador lineal, se tiene que:
x2 D2 = Dz (Dz 1)
As, gracias al cambio de variable z = ln(x) y el uso de la regla de la cadena, la ecuacion de Euler
de coeficientes variables se transforma en una ecuacion de coeficientes constantes (se eliminan
las x), por lo que podemos buscar una solucion homogenea con los metodos ya conocidos.Asi
la ED
a2 x2 y + a1 xy + a0 y = h(x)
59
quedara

a2 x2 D2 + a1 xD + a0 (y(x)) = h(x) a2 x2 Dz (Dz 1) + a1 xDz + a0 (y(z)) = h(ez )
Es util recordar que

xD = Dz
x2 D2 = Dz (Dz 1)
x3 D3 = Dz (Dz 1)(Dz 2)
x4 D4 = Dz (Dz 1)(Dz 2)(Dz 3)
Asi, dada la ecuacion

a2 x2 y + a1 xy + a0 y = 0
b es a2 D(D 1) + a1 D + a0
Entonces el operador diferencial lineal de coeficientes constantes L
y la solucion esta determinada por las races 1 y 2 de la ecuacion cuadratica
a2 m(m 1) + a1 m + a0 = 0, conocida como la ecuacion indicial. As,
i) Si 1 y 2 son races reales y distintas, entonces
yh (x) = C1 |x|1 + C2 |x|2 ;
ii) Si = 1 = 2 son races reales e iguales, entonces
yh (x) = |x| (C1 + C2 ln|x|) ;
iii) Si 1 = a + bi, y 2 = a bi y b > 0 son races complejas conjugadas, entonces
yh (x) = |x|a (C1 sen(b ln|x|) + C2 cos(b ln|x|)) .

x2 y + 2xy 6y = 0
Solucion La ecuacion indicial asociada es
m(m 1) + 2m 6 = 0
cuyas races son reales y distintas: {2, 3}
yh (x) = C1 |x|2 + C2 |x|3
60

xy + y = 0
Solucion multiplicando por x queda
x2 y + xy = 0
La ecuacion indicial asociada es
m(m 1) + m = 0 m2 = 0
cuyas races son reales e iguales: {0, 0}
yh (x) = C1 + C2 ln|x|

x2 y + y + 4y = 0
m(m 1) + m + 4 = 0
cuyas races son complejas: {0 + 2i, 0 2i}
yh (x) = |x|0 (C1 sen(2 ln|x|) + C2 cos(2 ln|x|)) = C1 sen(2 ln|x|) + C2 cos(2 ln|x|).

x2 y + 3xy + y = 0
m(m 1) + 3m + 1 = 0 m2 + 2m + 1 = 0 (m + 1)2 = 0
cuyas races e iguales: {1, 1}
yh (x) = |x|1 (C1 + C2 ln|x|)) .
1. x2 y 2xy + 2y = 0 3. x2 y + 10xy + 20y = 0
2. x2 y + 10xy + 8y = 0 4. x2 y + xy y = 0
61
5. x2 y + 2xy 2y = 0 9. 2x2 y + xy y = 0
6. x2 y + xy + 9y = 0 10. 4x2 y 4xy + 3y = 0
7. 4x2 y 8xy + 9y = 0 11. 4x2 y 8xy + 8y = 0
8. x2 y 3xy + 7y = 0 12. x2 y + 2xy 12y = 0
EDL de Euler de mayor orden

x3 y 2x2 y 17xy 7y = 0
m(m 1)(m 2) 2m(m 1) 17m 7 = m3 5m2 13m 7 = 0 (m + 1)2 (m 7) = 0
tiene las siguientes races reales: {1, 1, 7}
yh (x) = C1 |x|1 + C2 |x|1 ln|x| + C3 x7 .
Nota: cada vez que hay una raz repetida, la segunda solucion va multiplicada por el factor
ln|x|, en la ED de Euler. Asi, por ejemplo si la races de una ED de Euler de orden 6 cuyas
races son {2, 2, 4, 4, 4, 4} la solucion homogenea asociada es
yh (x) = C1 x2 + C2 x2 ln|x| + C3 x4 + C4 x4 ln|x| + C5 x4 ln2 |x| + C6 x4 ln3 |x|.
1. x3 y + 4x2 y 2y = 0 3. x3 y + 4x2 y 8xy + 8y = 0
2. x3 y 3x2 y + 6xy 6y = 0 4. x4 y (4) + 6x3 y + 7x2 y + xy y = 0
Ejemplo 8.16 Encuentre la solucion general de la siguiente ecuacion de Euler no homogenea:
x2 y + 10xy + 8y = x2
m(m 1) + 10m + 8 = m2 + 9m + 8 = 0
tiene las siguientes races reales: {8, 1}
yh (x) = C1 |x|8 + C2 |x|1 .
62
para buscar la solucion particular se puede usar el MVP o transformando la ED de Euler a una
ED de CC
(D + 1)(D + 8)(y) = e2z
como
(D 2)(e2z ) = 0 yp = Ae2z

D2 + 9D + 8 (Ae2z ) = D2 4D + 4 (Ae2z ) + (13D + 4) (Ae2z )

D2 + 9D + 8 (Ae2z ) = D2 4D + 4 (Ae2z ) + (13D + 4) (Ae2z )
| {z }
=0

D2 + 9D + 8 (Ae2z ) = (13D + 4) (Ae2z ) = 26Ae2z + 4Ae2z = 30Ae2z e2z
1 1 1
A= yp (z) = e2z yp (x) = x2
30 30 30
1
yG (x) = C1 |x|8 + C2 |x|1 + x2 .
30
Solucion por MVP:
yh (x) = C1 x8 + C2 x1
yp (x) = K1 (x)x8 + K2 (x)x1

K1 (x)x8 + K2 (x)x1 = 0
x2
8K1 (x)x9 K2 (x)x2 = x2
x8 x1
W = W [y1(x), y2 (x)] = 9 2
= 7x10
8x x
1 1
K1 = x9 K2 = x2
7 7
1 1
K1 = x10 K2 = x3
70 21
1 1 1
yp (x) = K1 (x)x8 + K2 (x)x1 = x10 x8 + x3 x1 = x2
70 21 30
1
yG (x) = C1 |x|8 + C2 |x|1 + x2
30
Ejemplo 8.17 Encuentre la solucion general de la siguiente ecuacion no homogenea:
x4 y + x3 y x2 y = 1
Solucion
x4 y + x3 y x2 y = 0 x2 y + xy y = 0
63
m(m 1) + m 1 = m2 1 = 0
yh (x) = C1 x + C2 |x|1 .
para buscar la solucion particular se debe usar
yh (x) = C1 x + C2 x1
yp (x) = K1 (x)x + K2 (x)x1

K1 (x)x + K2 (x)x1 = 0
1
K1 (x) K2 (x)x2 = x4
x x1 2
W = W [y1 (x), y2 (x)] = =
1 x 2 x
1 1
K1 = x4K2 = x2
2 2
1 1
K1 = 3 K2 =
6x 2x
1
yp (x) = x2
3
1 1 1
yG (x) = C1 x + C2 +
x 3 x2
Ejemplo 8.18 Encuentre la solucion general de la siguiente ecuacion no homogenea:
x2 y + 10xy + 20y = 4lnx x
Solucion Luego, con el cambio de variable z = ln(x), la ecuacion de coeficientes variables se

transforma en una de coeficientes constantes

D2 + 9D + 20 (y(z)) = 4z ez
m2 + 9m + 20 = (m + 4)(m + 5) = 0
yh (z) = C1 e4z + C2 e5z .
64
para buscar la solucion particular se puede (por MCI) suponer que
yp (z) = A + Bz + Cez
aplicando superposicion
yp1 (z) = A + Bz
(D2 (A + Bz) +9D + 20)(A + Bz) = 9B + 20A + 20Bz 4z

| {z }
=0
9
20A + 9B = 0 A = 100 9 1
1
yp1 (x) = + z
20B = 4 B = 5
100 5
yp2 (z) = Cez

D2 + 9D + 20 (Cez ) = D2 2D + 1 (Cez ) + (11D + 19) (Cez )
| {z }
=0

D + 9D + 20 (Ce ) = (11D + 19) (Cez ) = 11Cez + 19Cez ez
2 z
1 1
C= yp2 (z) = ez
30 30
9 1 1
yp (z) = yp1 (z) + yp2 (z) = + z + ez
100 5 30
9 1 1
yG (z) = C1 e4z + C2 e5z + z + ez
100 5 30
9 1 1
yG (x) = C1 |x|4 + C2 |x|5 + ln|x| + x
100 5 30
1. x2 y 2xy + 2y = ln2 x + ln(x2 ) 8. x2 y + 7xy + 5y = x
2. x2 y + xy y = x2 9. x2 y + 3xy 3y = 5x2
3. x2 y + xy 9y = x3 + 1 10. x2 y 2y = lnx
4. x2 y + xy 9y = sen(lnx3 )
11. 4x2 y 4xy + 3y = sen(lnx)
5. x2 y + 4xy + 2y = 2lnx 1 1 1
12. y + y 2 y = 2
6. x3 y + 4x2 y + xy + y = x x x x + x3
2 2 lnx
7. x2 y + xy + 4y = 1 + cos(2lnx) 13. y y + 2 y =
x x x
65
REFERENCIAS REFERENCIAS
Referencias
[1] Ecuaciones Diferenciales. Kreider-Kuller-Ostberg. Fondo Educativo Interamericano S.A,

1973.
[2] Partial Differential Equations in Action: From Modelling to Theory. Sandro Salsa. Springer
Verlag, 2008.
[3] Ordinary and Partial Differential Equations. R. P. Agarwal & D. ORegan. Universitext,
Springer Verlag, 2009.
66

Ecuacionesdiferenciales LHV 0

Cargado por

Copyright:

Formatos disponibles

Ecuacionesdiferenciales LHV 0

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ecuacionesdiferenciales LHV 0

Cargado por

Copyright:

Formatos disponibles

Universidad Tecnica Federico Santa Mara

donde la incognita y = y(x) es una funcion de la variable independiente x.

Donde las funciones ai (x) R, i = 1, . . . , n son llamados los coeficientes de la EDO. Si

2. Ecuaciones Diferenciales de Primer Orden

2.1. Integracion Directa

Supongamos ahora que queremos encontrar la solucion de 3 definida sobre un intervalo I y

2.2. Variables Separables

Definicion 2.2 Estas EDOs se escriben como

donde k = eC > 0. Como vemos, en este caso la solucion es implcita.

2.2.1. Ecuaciones Diferenciales Homogeneas de Primer Orden

f (x, y) = k f (x, y), R \ {0}.

Notemos que si contamos con dos funciones f, g homogeneas de grado k, g 6= 0. Entonces,

Si introducimos esto en 8 se obtiene

Ejemplo 2.7 Consideremos la EDO

Volviendo a las variables originales

Volviendo a las variables originales

2.2.2. Ejercicios Propuestos

2. (2x + 3y)dx + (x 2y)dy = 0 8. xdt (t2 a2 )dx = 0

5. (1 + y)dt (1 t)dy = 0 11. d + tg d = 0

14. 3ex tgy dx + (1 ex )sec2 y dy = 0 16. tgx sen2 y dx + cos2 x ctgy dy = 0

2.2.3. Ecuaciones Reducibles a homogeneas.

Si se tiene que la ecuacion posee la forma

a) a/b 6= /. Es decir, rectas son no paralelas. Se busca el punto de interseccion (x0 , y0 )

b) a/b = /. Las rectas son paralelas, pero esto implica que

Ejemplo 2.11 Consideramos

2.3. Ecuaciones Diferenciales Exactas

dF = M(x, y)dx + N(x, y)dy.

df = M(x, y)dx + N(x, y)dy = 0, (x, y) R. (12)

Ejemplo 2.12 La ecuacion

Ejemplo 2.13 La ecuacion (x2 + y)dx + (x 2y)dy = 0 es exacta pues

(3x2 y + 2xy)dx + (x3 + x2 + 2y)dy = 0

Sera exacta? si pues

M(x, y) = 3x2 y + 2xy, N(x, y) = x3 + x2 + 2y

Ejemplo 2.15 Resolver

Ejemplo 2.16 Resolver la EDO siguiente

La EDO es exacta pues:

2.3.1. Ejercicios Propuestos

2.3.2. Factor de Integracion

Supongamos que la ecuacion 12 no es exacta. Para poder escribirla en dicha forma, se

integrante) de modo que la ecuacion

(M) dx + (N) dy = 0, (14)

resulte ser exacta. Al realizar un calculo sencillo se llega a la ecuacion

Luego el factor integrante depende de y. Este factor confirma la EDO

2.3.3. Ejercicios Propuestos

6. y(1 + xy)dx = xdy

10. (x2 + y 2 + x)dx + ydy = 0 14. xydx + (x2 + 2y 2 + 2)dy = 0

Ejemplo 2.19 Resolver la ecuacion

(x3 + xy 2 y)dx + (x2 y + y 3 + x)dy = 0

si sabe que el factor integrante es = (x2 + y 2 ), podemos ver que

2.4. Ecuaciones Diferenciales Lineales

Una ecuacion de la forma

y + p(x) y = q(x), (16)

donde p(x) = a0 (x)/a1 (x) y q(x) = h(x)/a1 (x).

Teorema 2.2. Principio de Linealidad. La solucion general de 16 esta dada por

y(x) = yh (x) + yp (x), (17)

donde yh = yh (x) es la solucion del problema homogeneo

llamada solucion homogenea e yp = yp (x) es una solucion particular de

Ya que yh es solucion del problema homogeneo. Entonces, yp es la solucion particular y se tiene