Ecuaciones

Ecuaciones Diferenciales Ordinarias.
Verónica Bolaños Granados
11 de noviembre de 2021
2
Índice general
1. Ec Dif Ord de Primer Orden 5

1.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2. Ecuaciones Diferenciales Ordinarias Lineales de Primer Orden . . . . 6
1.2.1. Ejemplos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3. Ecuaciones diferenciales lineales de primer orden de variables separables 13
1.4. Ecuaciones Exactas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.4.1. Factor Integrante . . . . . . . . . . . . . . . . . . . . . . . . . 24
2. Ec Dif Ord de Orden n 29

2.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.2. Ecuación diferencial lineal homogénea con coeficientes constantes de
segundo orden . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.3. Ecuación diferencial lineal homogénea con coeficientes constantes de
segundo orden con condiciones iniciales . . . . . . . . . . . . . . . . . 31
2.4. Dependencia e independencia lineal . . . . . . . . . . . . . . . . . . . 36
2.5. Una formula para el Wronskiano . . . . . . . . . . . . . . . . . . . . . 39
2.6. La ecuación diferencial no homogénea de segundo orden . . . . . . . 40
2.7. Ecuación diferencial homogenea de orden n . . . . . . . . . . . . . . . 42
2.8. Ecuaciones de orden n con condiciones iniciales . . . . . . . . . . . . 44
2.9. Ecuaciones cuyas constantes son reales . . . . . . . . . . . . . . . . . 50
2.10. La ecuación no homogenea de orden n. . . . . . . . . . . . . . . . . . 51
2.11. Método de coeficientes indeterminados para resolver ecuaciones no
homogéneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.11.1. Álgebra de los operadores con coeficientes constantes . . . . . 54
2.11.2. Sobre el Método . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3. Ecuaciones Lineales con Coeficientes Variables 59

3.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
3
4 Índice general
3.2. Problemas con valores iniciales para ecuaciones homogeneas . . . . . 60

3.3. Soluciones de Ecuaciones Homogéneas. . . . . . . . . . . . . . . . . . 62
3.4. Wronskiano e Independencia Lineal . . . . . . . . . . . . . . . . . . . 64
3.5. Reducción del orden de una ecuación homogénea . . . . . . . . . . . . 69
3.6. Ecuación no Homogenea . . . . . . . . . . . . . . . . . . . . . . . . . 71
3.7. Solución con Series de Potencias . . . . . . . . . . . . . . . . . . . . . 72
Capı́tulo 1
Ecuación Diferencial Ordinaria de

Primer Orden.
Terminologı́a y Definiciones
1.1. Introducción
En cálculo se aprende que la derivada dy/dx o y 0 de la función y = φ(x) es otra
2 2
función de x. Por ejemplo si y = ex , entonces y 0 = dy/dx = 2xex . Al reemplazar
2
ex por y se obtiene
y 0 = 2xy. (1.1)
El problema que tenemos que resolver en el curso será dada una ecuación como
1.1, ¿hay algún método para llegar a la función desconocida y = φ(x)?
Veamos primero la definición de ecuación diferencial.
Definición 1 (Ecuación Diferencial) Una ecuación que contiene las derivadas de

una o más variables dependientes con respecto a una o más variables independientes
es una ecuación diferencial.
Como ejemplo tenemos la ecuaciones:

p
y 000 = 1 + (y 00 )2 .
y 00 =2y 3
0
ey =1
5
6 1.2. Ecuaciones Diferenciales Ordinarias Lineales de Primer Orden
Las ecuaciones diferenciales se clasifican de acuerdo a su tipo , orden y linea-

lidad.
Clasificación según el tipo
Si una ecuación sólo tiene derivadas ordinarias de una o más variables depen-
dientes con respecto a una sola variable independiente entonces se dice que es una
ecuación diferencial ordinaria. Por ejemplo
y 0 = ex − 10y y y 00 − y 0 + 6y = 0.
Una ecuación que contiene las derivadas parciales de una o más variables de-
pendientes, respecto de dos o mas variables independientes se llama ecuación en
derivadas parciales. Por ejemplo
∂u ∂v ∂ 2u ∂ 2u ∂u
=− y 2
= 2
−2 .
∂y ∂x ∂x ∂t ∂t
Clasificación según el orden El orden de una ecuación diferencial es el de la
derivada de mayor orden de la ecuación. Por ejemplo
segundo orden ↓ ↓ primer orden

d2 y dy 3

dx2
+5 dx
− 4y = ex
Clasificación según linealidad se dice que una ecuación diferencial de la forma
y (n) = f (x, y, y 0 , . . . , y (n−1) )

es lineal cuando se puede escribir de la forma
an (x)y (n) + an−1 (x)y (n−1) + . . . + a1 (x)y 0 + a0 (x)y = g(x).
1.2. Ecuaciones Diferenciales Ordinarias Lineales

de Primer Orden
Definición 2 Supongamos que f es una función de variable compleja definida para
toda x ∈ R en un intervalo I, y para todos y ∈ S ⊆ C el valor de f en (x, y),
se representa por f (x, y). Una ecuación diferencial ordinaria de primer orden se
representa
y 0 = f (x, y). (1.2)
Capı́tulo 1. Ec Dif Ord de Primer Orden 7
Consideremos una ecuacion diferencial lineal de primer orden.
y 0 + a(x)y = b(x) (1.3)

donde a, b son funciones definidas en I ∈ R. Escribiendo lo anterior en la forma
0
y = f (x, y) tenemos
f (x, y) = −a(x)y + b(x). (1.4)

Cuando b(x) = 0 en 1.3 la ecuación se llama homogénea, mientras que si no es
identicamente cero la ecuación se llama no homogénea.
Tenemos el siguiente teorema.
Teorema 1 Considere la ecuación diferencial ordinaria lineal homogénea
y 0 + ay = 0. (1.5)
donde a es una constante compleja. Si c ∈ C cualquiera, la función definida por
φ(x) = ce−ax
es una solución a esta ecuación, y además toda solución tiene esta forma.
Demostración
Sea c una constante arbitraria, entonces la función definida por 1.6 es solución
de 1.5, ya que
φ0 (x) + aφ(x) = −ace−ax + ace−ax = 0.

Reciprocamente, sea φ una solución de 1.5, entonces φ0 + aφ = 0, esto implica
eax (φ0 + aφ) = 0,
(eax φ)0 = 0.
Por lo que existe una constante c ∈ C tal que eax φ = c o
φ(x) = ce−ax . (1.6)

Por lo que se ha demostrado que toda solución de 1.5 tiene la forma de 1.6.

Consideremos ahora cuando b(x) 6= 0. Tenemos el siguiente teorema.
Teorema 2 Considere la ecuación diferencial ordinaria lineal no homogénea
y 0 + ay = b(x). (1.7)
donde a es una constante y b es una función continua en un intervalo I. Si x0 ∈ I,
y c es una constante cualquiera, entonces la función φ definida por
Z x
−ax
φ(x) = e eat b(t)dt + ce−ax (1.8)
x0
es una solución de 1.7, y cualquier otra solución tiene esta forma.
Demostración
Sea c una constante arbitraria sea φ dada por la ecuación 1.8. Derivando tenemos
Z x
0 −ax
φ (x) = −ae eat b(t)dt + e−ax [eax b(x)] − ace−ax
x0
Sustituyendo en 1.7 tenemos

Rx Rx
−ae−ax x0 eat b(t)dt + e−ax [eax b(x)] − ace−ax + ae−ax x0 eat b(t)dt + cae−ax
= e−ax+ax b(x)
= b(x)
Y Recı́procamente, como en el teorema anterior sea φ una solución de 1.7, entonces
eax (φ0 + aφ) = eax b,
o
(eax φ)0 = eax b.
Sea B una función tal que B 0 (x) = eax b(x), como por ejemplo:
Z x
B(x) = eat b(t)dt,
x0
donde x0 ∈ I es un punto fijo. (Tal función existe debido a la continuidad de la
exponencial y de b)
Como eax φ es otra función cuya derivada es eax b, se tiene que
eax φ(x) = B(x) + c

para alguna constante c, por consiguiente
φ(x) = e−ax B(x) + ce−ax (1.9)

1.2.1. Ejemplos
Ejemplos 1 Considere la ecuación
Ly 0 + Ry = E sen ωx,
donde L, R, E, ω son constantes positivas.
(a) Calcule la solución φ que satisfaga la condición φ(0) = 0.
(b) Demuestre que la solución puede escribirse de la forma:
EωL E
φ(x) = e−(R/L)x + √ sen(ωx − α),
R2 2
+ω L 2
R + ω 2 L2
2
donde α es el ángulo que satisface las siguientes ecuaciones:
R ωL
cos α = , sen α = 2 .
R2 2
+ω L2 R + ω 2 L2
Solución
(a) Por ser las constantes positivas, tenemos que la ecuación es equivalente a
R E sen ωx
y0 + y= ,
L L
de donde tenemos lo siguiente.

R Rx R R
φ(x) = e− L x x0 e L t EL sen ωtdt + ce− L x
R Rx R R
= EL e− L x x0 e L t EL sen ωtdt + ce− L x
por otra parte

R R
e L t EL sen ωtdt = E
R R
L
e L t sen ωtdt
R R
resolviendo por partes e L t sen ωtdt tenemos
R R
u = eLt du = R
L
e L t dt
dv = sen ωtdt v = −1
ω
cos ωt
−1 R R
e L t cos ωt − −1 cos ωt R
R
ω ωR L
e L t dt =
−1 R R R
ω
e L t cos ωt + ωL e L t cos ωtdt =
R R
u = eLt du = RL
e L t dt
1
dv = cos ωtdt v = ω sen ωt
h i
−1 R R 1 R 1 R
ωt R
R
ω
eLt cos ωt + ωL ω
eLt sen ωt − ω
sen L
e L t dt
Pasando la integral al primer miembro de la igualdad y simplificando

h iR R
R2 R eLt R

1+ ω 2 L2
e L t sen ωtdt = ω
− cos ωt + ωL
sen ωt
De donde obtenemos que

R
R ω 2 L2 eLt R
R
e L t sen ωtdt = ω 2 L2 +R2 ω
− cos ωt + ωL
sen ωt
Por lo que
E −R
h
ωL2 R ix R
φ(x) = L
e Lx ω 2 L2 +R2
eLt − cos ωt +
sen ωt + ce− L x =
R
ωL
t=x0
EωL R −R

ω 2 L2 +R2
− cos ωx + ωL sen ωx + Ce x L
la C es la agrupación de la evaluación de xo con la c
con la condición φ(0) = 0 tenemos
0 = φ(0) = − ω2EωL
L2 +R2
+C ⇒C = EωL
ω 2 L2 +R2
R

Por lo tanto φ(x) = EωL
ω 2 L2 +R2
e− L x − cos ωx + R
ωL
sen ωx

Por último consideremos la ecuación diferencial lineal no homogenea de primer
orden con coeficientes variables. Tenemos es siguiente teorema.
Teorema 3 Supongamos que a y b son dos funciones continuas en un intervalo

I ∈ R. Sea A una función tal que A0 = a. Entonces la funcion dada por
Z x
−A(x)
ψ(x) = e eA(t) b(t)dt,
x0
donde xo ∈ I es una solución de la ecuación
y 0 + a(x)y = b(x) (1.10)
en I. La función φ1 dada por

φ1 (x) = e−A(x)
es una solución de la ecuación homogénea
y 0 + a(x)y = 0.
Si c es una constante cualquiera, entonces φ = ψ + cφ1 es una solución de 1.10,

y toda solución de 1.10 tiene esta forma.
Demostración
Excluiremos los puntos singulares (a(x)=0) de la función a. Ası́ supongamos que
a(x) 6= 0 ∀x ∈ I. Supongamos que φ es una solución de 1.10. Notemos que la función
dada por eA φ cumple con
eA (φ0 + aφ) = (eA φ)0 . (1.11)
Si A es una función cuya derivada es a, como por ejemplo
Z x
A(x) = a(t)dt,
x0
donde x0 ∈ I es un punto fijo,dado que
(eA φ)0 = eA φ0 + aeA φ = eA (φ0 + aφ).
Observemos que se tendrá φ0 + aφ = b si y sólo si
eA φ = B + c, (1.12)
donde c es una constante, y B es una función cuya derivada es eA b. Podemos

escoger B dada por Z x
B(x) = eA(t) b(t)dt.
x0
Para demostrar la necesidad y suficiencia tenemos

⇒) (necesidad) Sustituyendo en (1.11) tenemos
(eA φ)0 = eA (φ0 + aφ) = eA b,
la última igualdad es por la hipótesis.

⇐) (suficiencia) Para probar que las dos funciones son iguales lo probaremos de
forma puntual, para ello sea x ∈ I arbitrario, utilizando (1.11) tenemos
(eA (φ0 + aφ))(x) = (eA φ)0 (x) = (eA b)(x)
la última igualdad es por la hipótesis. entonces
eA(x) (φ0 + aφ)(x) = eA(x) b(x)
como eα > 0 tenemos (φ0 +aφ)(x) = b(x) y como x fue arbitrario tenemos φ0 +aφ = b
Por otro lado 1.12 es válida si y sólo si
φ(x) = eA(x) B(x) + ce−A(x) , (1.13)

Para demostrar la necesidad y suficiencia tenemos
⇒) Tenemos
Z x Z x Z x
0
A
e φ= A
(e φ) (t)dt = A
(e b)(t)dt = B 0 (t)dt = B|xx0 = B + c,
x0 x0 x0
la segunda igualdad es por la hipótesis, la primera y la penúltima por el TFC.

⇐)
(eA φ)0 = (B + c)0 = dB + dc = eA b.
la primera igualdad es por la hipótesis y la última por ser c constante.
Ası́ se ha demostrado que toda solución de 1.10 tiene la forma de 1.12, y recı́pro-
camente, que si c es una constante cualquiera la función definida por 1.12 es una
solución de 1.10.

La solución de 1.10, es una solución particular de la ecuación y la solución general
de la homogénea.
EJEMPLOS
1.3. Ecuaciones diferenciales lineales de primer or-

den de variables separables
Se dice que una ecuación de primer orden del tipo 1.2 es de variables separables
si puede escribirse de la forma:
g(x)
f (x, y) = ,
h(y)
donde g, h son funciones de una sola variable. En este caso se puede escribir la
ecuación ası́
dy
h(y) = g(x), (1.14)
dx
o bien:
h(y)dy = g(x)dx,
lo que nos da el origen del término.
Estudiaremos la Ecuación 1.14 para el caso en que g y h sean funciones continuas
de variable real definidas para x y para y reales respectivamente.
Teorema 4 Sean g y h, dos funciones continuas de variable real definidas en los

intervalos a ≤ x ≤ b, c ≤ y ≤ d respectivamente, y consideremos la ecuación
h(y)y 0 = g(x) (1.15)
Si H y G son dos funciones cualesquiera, tales que G0 = g, H 0 = h, y c es

cualquier constante tal que la relación
H(y) = G(x) + c
define para x, en cierto intervalo I, contenido en a ≤ x ≤ b, a una función derivable

φ, de variable real, entonces φ será una solución de 1.15 en I. Reciprocamente, si φ
es una solución de 1.15 en I esta satisfacerá la relación:
H(y) = G(x) + c
en el intervalo I, para cierta constante c.

14 1.3. Ecuaciones diferenciales lineales de primer orden de variables separables
Demostración
Si φ es una solución de variable real para la Ec. 1.15, que está definida en cierto
intervalo I, el cual contiene un punto x0 , entonces se cumple
h(φ(x))φ0 (x) = g(x)
para toda x ∈ I, y por consiguiente:

Z x Z x
0
h(φ(t))φ (t)dt = g(t)dt (1.16)
x0 x0
para toda x en I. Haciendo u = φ(t) en la integral que aparece en el miembro

izquierdo de 1.16, se puede escribir dicha ecuación de la forma
Z φ(x) Z x
h(u)du = g(t)dt
φ(x0 ) x0
Recı́procamente, supongamos que x y y están relacionados por la fórmula

Z y Z x
h(u)du = g(t)dt (1.17)
y0 x0
y que esta define en I una función implı́cita φ de x, que además es derivable. Entonces
esta función satisface la igualdad:
Z φ(x) Z x
h(u)du = g(t)dt
y0 x0
para toda x en I, ası́ que derivando obtenemos:
h(φ(x))φ0 (x) = g(x)
lo que demuestra que φ es una solución de (1.16) definida en I.

En la práctica, el método usual de trabajar con la Ecuación 1.16 es escribiendola
de la siguiente manera:
h(y)dy = g(x)dx
(esto es, separando las variables), y luego se integra para obtener:
Z Z
h(y)dy = g(x)dx + c,
donde c es una constante, y las integrales se interpretan como antiderivadas. Ası́:

Z Z
H(y) = h(y)dy, G(x) = g(x)dx,
representan dos funciones cualesquiera H, G, tales que
H 0 = h, G0 = g
Entonces toda función derivable φ que esté definida implı́citamente por la relación
H(y) = G(x) + c (1.18)

será una solución de (1.15).
Ejemplos 2 1. El ejemplo más sencillo es el caso en el cual h(y) = 1. Entonces
y 0 = g(x), y toda solución φ tiene la forma
φ(y) = G(x) + c (1.19)
donde G es una función cualquiera, definida en a ≤ x ≤ b, tal que G0 = g,

y c es una constante. Más todavı́a, si c es una constante cualquiera, (1.18)
define una solución de y 0 = g(x). Ası́ se han determinado todas las soluciones
de y 0 = g(x) en a ≤ x ≤ b.
2. Otro caso ocurre cuando g(x) = 1, ya que tenemos
1
y0 = (1.20)
h(y)
a sea:
h(y)dy = dx.
Ası́, si H 0 = h, cualquier función derivable definida implı́citamente por la re-
lación
H(y) = x + c, (1.21)
donde c es una constante, será una solución de (1.20).
3. Consideremos la ecuación
y0 = y2. (1.22)
Aquı́ h()y = 1/y 2 , la cual observamos que no es continua en y = 0. Tenemos
dy
= dx
y2
16 1.3. Ecuaciones diferenciales lineales de primer orden de variables separables
y la relación (1.21) se transforma en
1 −1
− = x + c, o y= .
y x+c
Sı́ c es una constante cualquiera, la función dada por

−1
φ(x) = . (1.23)
x+c
es una solución de (1.22), siempre y cuando x 6= c.
4. sea
y 0 = 3y 2/3 . (1.24)
Esto da lugar a :
dy
= 3dx
y 3/2
si y 6= 0, y en consecuencia,
y 1/3 = x + c, o y = (x + c)3
donde c es una constante. La función dada por
φ(x) = (x + c)3 . (1.25)
será una solución de (1.24) para toda constante c. Observemos que la función
identicamente cero es una solución de (1.24), y dicha función no puede obte-
nerse a partir de (1.25).
Es importante recalcar que el método de separación de variables, puede no dar

todas las soluciones de una ecuación. Por ejemplo, de la ec.(1.22) es claro que la
función ψ, la cual es idénticamente nula para toda x es una solución de (1.22). Sin
embargo, la φ de (1.23) no es igual a dicha solución para ningún valor de la constante
c. Frecuentemente, la atención cuidadosa frente a las posibilidades de una división
entre cero, advierte acerca de la existencia de otras soluciones adicionales.
El ejemplo (1.24) ilustra una dificultad adicional que se presenta al trabajar
con ecuaciones no lineales, esto es puede haber varias soluciones que satisfagan una
condición inicial dada. Ası́ las dos funciones φ y ψ dadas por
φ(x) = x3 , ψ(x) = 0, (−∞ < x < ∞),

son soluciones de (1.24) que pasan por el origen. En realidad, la situación es bastante
más difı́cil de lo que parece ya que existe una infinidad de funciones, las cuales son
soluciones de (1.24) que pasan por el origen. Para ver esto, supongamos que k es
cualquier entero positivo, y definamos φk de la siguiente manera
φk (x) = 0, (−∞ < x ≤ k),
φk (x) = (x − k)3 , (k < x < ∞).

Entonces no es difı́cil darse cuenta de que φk es una solución de (1.24) para toda x
real, y es claro que φk (0) = 0.
1.4. Ecuaciones Exactas

Daremos primero la definición de ecuación exacta.
Definición 3 Supongamos que la ecuación diferencial y 0 = f (x, y) se escribe de la

siguiente manera:
−M (x, y)
y0 = ,
N (x, y)
o, en forma equivalente:
M (x, y) + N (x, y)y 0 = 0, (1.26)
donde M, N , son funciones de variable real, definidas para x, y, reales en cierto
rectángulo R. Se dice que la Ec.(1.26) es exacta en R, si existe una función F que
contenga sus derivadas parciales continuas tales que
∂F ∂F
= M, =N (1.27)
∂x ∂y
en R
Teorema 5 Supongamos que la ecuación
M (x, y) + N (x, y)y 0 = 0
es exacta en un rectángulo R, y que F es una función de variable real que satisface

las siguientes condiciones:
∂F ∂F
= M, =N
∂x ∂y
18 1.4. Ecuaciones Exactas
en R. Toda función derivable φ que esté definida implı́citamente por la relación:
F (x, φ(y)) = c, (c = constante),
es una solución de la ecuación y recı́procamente toda solución de la ecuación cuya

gráfica esté en R puede obtenerse de esta manera.
Demostración
Si (1.26) es exacta en R y F es una función que satisface (1.27), entonces (1.26)
se transforma en
∂F ∂F
(x, y) + (x, y)y 0 = 0
∂x ∂y
Si φ es cualquier solución definida en cierto intervaloI, entonces:
∂F ∂F
(x, φ(x)) + (x, φ(x))φ0 (x) = 0 (1.28)
∂x ∂y
para toda x ∈ I. Si φ(x) = F (x, φ(x)) entonces la Ec.(1.28) significa precisamente
que φ0 (x) = 0, y en consecuencia
F (x, φ(x)) = c,
donde c es cierta constante. Ası́, la solución φ debe ser una función que esté dada
implı́citamente por la relación
F (x, y) = c (1.29)
Desarrollando este mismo argumento en orden inverso, vemos que si φ es una fun-
ción derivable en cierto intervalo I, que está definida implı́citamente por la relación
(1.29),
F (x, φ(x)) = c
para toda x ∈ I, y derivando se obtiene (1.28). Ası́, φ es una solución de (1.26).

El problema de resolver una ecuación exacta, queda reducido ahora al problema
de hallar una función F que satisfaga (1.27). Si (1.26) es exacta y la escribimos de
la siguiente manera:
∂F ∂F
M (x, y)dx + N (x, y)dy = (x, y)dx + (x, y)dy = 0
∂x ∂y
reconocemos en el miembro del centro de la ecuación a la diferencial dF de la función
F . Esta es la razón por la que a esta tipo de ecuaciones se les denomina “exactas”,
el miembro central es una diferencial exacta de la función F .
Algunas veces puede determinarse por inspección una F . Por ejemplo, la ecuación
x
y0 = − (1.30)
y
se escribe de la forma
xdx + ydy = 0,
se ve claramente que el miembro izquierdo es la diferencial de (x2 + y 2 )/2. Ası́,
cualquier función diferenciable que esté definida por la relación
x2 + y 2 = c
es una solución de (1.30). Notemos que la Ec. (1.30) no está definida para y = 0.
El ejemplo anterior también es un ejemplo de una ecuación de variables separa-
bles. De hecho, toda ecuación de este tipo es un caso especial de una ecuación exacta,
por lo cual si escribimos la ecuación de la siguiente manera:
g(x)dx = h(y)dy,
es claro que existe una F dada por
F (x, y) = G(x) − H(y),
donde G0 = g, H 0 = h.
El siguiente teorema nos muestra como reconocer cuando una ecuación es exacta.
Teorema 6 Sean N, M dos funciones de variable real, las cuales tienen primeras
derivadas parciales continuas en cierto rectángulo R : |x − x0 | ≤ a, |y − y0 | ≤ b.
Entonces la ecuación
M (x, y) + N (x, y)y 0 = 0
es exacta en R, si y sólo si
∂M ∂N
= (1.31)
∂y ∂x
en R
Demostración
Supongamos que la Ec.(1.31) es exacta, y que F es una función con segundas
derivadas continuas tales que
∂F ∂F
= M, = N.
∂x ∂y
Entonces:
∂ 2F ∂M ∂ 2F ∂N
= , = .
∂y∂x ∂y ∂x∂y ∂x
y dado que para dicha función
∂ 2F ∂ 2F
= ,
∂y∂x ∂x∂y
debemos tener
∂M ∂N
=
∂y ∂x
Recı́procamente supongamos ahora que se comple (1.31) en R. Necesitamos hallar
una función F que satisfaga las siguientes condiciones:
∂F ∂F
= M, = N.
∂x ∂y
Para hacer esto, observemos que si tuvieramos tal función, entonces:
F (x, y) − F (x0 , y0 )=F

R x(x, y) − F (x0 , y)
R y+∂FF (x0 , y) − F (x0 , y0 )
= x0 ∂F (s, y)ds + (x , t)dt
R x ∂x R yy0 ∂y 0
= x0 M (s, y)ds + y0 N (x0 , t)dt.
En forma similar, debemos tener:
F (x, y) − F (x0 , y0 )=F

R y(x, y) − F (x, yR0 ) + F (x, y0 ) − F (x0 , y0 )
x
= y0 ∂y (x, t)dt + x0 ∂F
∂F
(s, y0 )ds
Ry R x ∂x
= y0 N (x, t)dt + x0 M (s, y0 )ds.
Ahora definimos F mediante la fórmula:
Z x Z y
F (x, y) = M (s, y)ds + N (x0 , t)dt. (1.32)
x0 y0
Esta definición implica que F (x0 , y0 ) = 0 y que

∂F
(x, y) = M (x, y)
∂x
para toda (x, y) en R. De acuerdo con (1.4) podemos suponer que F también está
dada por Z y Z x
F (x, y) = N (x, t)dt + M (s, y0 )ds. (1.33)
y0 x0
Esto de hecho es cierto, y es una consecuencia de la hipótesis establecida en (1.31),

luego por (1.33) se tiene
∂F
(x, y) = N (x, y)
∂y
para todo (x, y) en R, con lo que queda determinada F .
Para demostrar que (1.33) es válida, donde F es la función dada por (1.32),
consideremos la siguiente diferencia
Z y Z x
F (x, y) − N (x, t)dt − M (s, y0 )ds
y0 x0
Rx Ry
= x0 [M (s, y) − M (s, y0 )]ds − y0 [N (x, t) − N (x0 , t)dt
R x hR y ∂M i R y hR x ∂N i
= x0 y0 ∂y (s, t)dt ds + y0 x0 ∂x (s, t)ds dt
R x hR y i
= x0 y0 ∂M ∂y
− ∂N
∂x
(s, t) dsdt
la cual es nula en virtud de (1.31).

Ası́ tenemos una manera de encontrar la función F .
1. Intégrese M respecto de x, remplazando la constante de integración por una

función f (y).
2. Derı́vese la función F obtenida en el paso 1 con respecto a y e iguálese el

resultado con N . Con esto podemos calcular f (y) integrando f 0 (y), sin usar
una constante de integración.
3. Sustitúyase el valor de f (y) en la F obtenida en el paso 1 y finalmente, escrı́base

la solución de la ecuación de la forma F (x, y) = c con c una constante arbitraria.
Se puede hacer lo mismo integrando N respecto a y y siguendo el mismo proceso

ahora en función de la variable x.
Ejemplos 3 1. Consideremos la ecuación
3x2 − 2xy
y0 = , (1.34)
x2 − 2y
la cual puede escribirse de la siguiente manera:
(3x2 − 2xy)dx + (x2 − 2y)dy = 0.

Aquı́
M (x, y) = 3x2 − 2xy, N (x, y) = x2 − 2y,
y mediante cálculo directo se demuestra que
∂M ∂N
= = −2x,
∂y ∂x
tenemos
∂F ∂F
= M, = N.
∂x ∂y
Ası́, F satisface la siguiente condición:
∂F
(x, y) = 3x2 − 2xy,
∂x
lo cual implica que para cada y fija,
F (x, y) = x3 − x2 y + f (y) (1.35)
donde f es independiente de x. Ahora ∂F/∂y = N quiere decir que
−x2 + f 0 (y) = 2y − x2 ,
es decir
f 0 (y) = 2y.
Ası́, una selección para f está dada por f (y) = y 2 , de manera que sustituyendo
en (1.35), obtenemos
F (x, y) = x3 − x2 y + y 2 .
Cualquier función derivable φ que esté definida implı́citamente por una rela-
ción:
x3 − x2 y + y 2 = c, (1.36)
donde c es una constante, será una solución de (1.34), y todas las soluciones
de (1.34) pueden obtenerse de esta manera. Frecuentemente, las soluciones se
identifican con las relaciones (1.36).
−2xy
2. Consideremos la ecuación y 0 = 1+x2
, la cual puede escribirse como 2xydx +
(1 + x2 )dy, aquı́
M (x, y) = 2xy, N (x, y) = 1 + x2 ,

∂M ∂N
para saber que es exacta debemos comprobar que 2x = ∂y
= ∂x
= 2x entonces
F (x, y) = x2 y + f (y)
derivando respecto de y se tiene
x2 + f 0 (y) = 1 + x2 .
entonces f (y) = y luego la solución es
x2 y + y = c.
2x−3e3x y
3. Consideremos la ecuación y 0 = e3x
, la cual puede escribirse como (3e3x y −
2x)dx + e3x dy, aquı́
M (x, y) = 3e3x y − 2x, N (x, y) = e3x ,

∂M ∂N
para saber que es exacta debemos comprobar que 3e3x = ∂y
= ∂x
= 3e3x
entonces
F (x, y) = e3x y − x2 + f (y)
e3x + f 0 (y) = e3x
entonces f (y) = c luego la solución es
e3x y − x2 = c.
4. Consideremos la ecuación
y − 2y 2 x
y0 =
x
2
la cual puede escribirse como (2y x−y)dx+xdy = 0, es claro que esta ecuación
no es exacta. Al multiplicarla por y12 tenemos

1 x
2x − dx + 2 dy = 0
y y
Con esta nueva ecuación tenemos
1 x
M (x, y) = 2x − , N (x, y) = ,
y y2
1 ∂M ∂N 1
para saber que es exacta debemos comprobar que y2
= ∂y
= ∂x
= y2
entonces
x
F (x, y) = x2 − + f (y)
y

x 0 x
+ f (y) = .
y2 y2
entonces f (y) = c luego la solución es

x
x2 − = c.
y
1.4.1. Factor Integrante

Si una ecuación diferencial de primer orden no es exacta, se puede considerar la
posibilidad de convertirla en exacta mediante la multiplicación por un factor apro-
piado llamado factor integrante.
Método de inspección Deseamos obtener un método para encontrar un factor
integrante para una ecuación dada M dx+N dy = 0 la cual no es exacta. Supongamos
que P es dicho factor integrante. Entonces P M dx + P N dy = 0 es exacta y cumple
∂(P M ) ∂(P N )
= (1.37)
∂y ∂x
Se tienen los siguientes casos.

Caso 1 Cuando P es solamente función de x. De (1.37) tenemos
∂M ∂N ∂P
P =P +N (1.38)
∂y ∂x ∂x
Despejando tenemos
dP 1 ∂M ∂N
= − dx. (1.39)
P N ∂y ∂x
Como el coeficiente de dx en el segundo miembro de la Ec. (1.39) es una función
solamente de x, digamos f (x), se puede resolver la ecuación diferencial (1.39) por
separación de variables, tenemos
R
f (x)dx
P =e . (1.40)
Para aplicar este caso de factor integrante, es preciso calcular primero el lado derecho
de (1.39) y si resulta solamente función de x, se utilizala (1.40) para obtener P
Caso 2 Cuando P es solamente función de y. De manera similar tenemos
R
g(y)dy
P =e . (1.41)
donde
1 ∂N ∂M
g(y) = − dx. (1.42)
M ∂x ∂y
Caso 3 Cuando P es función de xy. De (1.37) se obtiene
∂ ∂
(P (x, y)M (x, y)) = (P (x, y)N (x, y))
∂y ∂x
Hagamos xy = u. Tenemos
∂M ∂P (u) ∂N ∂P (u)
P (u) +M = P (u) +N (1.43)
∂y ∂y ∂x ∂x
Como ∂P∂y(u) = P 0 (u) ∂u

∂y
= P 0 (u)x, y análogamente ∂P (u)
∂x
= P 0 (u)y la ecuación 1.43
se transforma en
∂M ∂N
P (u) + M xP 0 (u) = P (u) + N yP 0 (u),
∂y ∂x
de la cual se despeja
∂M
dP ∂y
− ∂N
∂x
= du. (1.44)
P Ny − Mx
Como el coeficiente de du en la ecuación 1.44 es función solamente de u = xy, di-
gamos h(u), se puede resolver la ecuación diferencial 1.44 por separacón de variables,
obteniéndose
R
h(u)du
P =e . (1.45)
y
Ejemplos 4 1. Consideremos la ecuación diferencial y 0 = − 3+3x−y . La ecuación
anterior puede verse como ydx + (3 + 3x − y)dy = 0. No es exacta. Calculando
el lado derecho de la ecuación 1.39 tenemos lo siguiente:

1 ∂M ∂N
N ∂y
− ∂x
1
= 3+3x−y
(1 − 3)
−2
= 3+3x−y
,
que no es función solo de x, luego el caso 1 no es aplicable. Calculemos ahora el
lado derecho de la ecuación 1.42, y ası́ obtenemos g(y) = y2 que sı́ está solo en
función de y. Entonces se trata del caso 2. Usando la ecuación 1.41 se obtiene
el factor integrante de la siguiente forma.
R
P = eR g(y)dy
2
= e y dy
= e2 ln y
2
= eln y
= y2,
Multiplicando nuestra ecuación original por este factor integrante, resulta la

ecuación exacta
y 3 dx + (3y 2 + 3xy 2 − y 3 )dy = 0.
Resolviendo por el método de ecuaciones exactas, tenemos
y4
y3x + y3 − = c.
4
2. Considere la ecuación y 0 = x − y. Esta ecuación la podemos representar de la

forma (x − y)dx − dy = 0. Esta ecuación no es exacta, por lo que calcularemos
el lado derecho de 1.41, con lo que obtenemos

1 ∂M ∂N
N ∂y
− ∂x
1
= −1
(−1 − 0)
= 1,
R
Ası́ obtenemos f (x) = 1, por lo que P = e f (x)dx = ex . Al multilicar la ecuación
original por el factor integrante obtenemos (ex x − ex y)dx − ex dy = 0 que es
exacta y cuya solución está dada por −ex y + xex − ex = c.
3 2
3. Consideremos la ecuación y 0 = − xy3 +xy +y
+x2 y+x
. Esta ecuación la podemos represen-
3 2 3 2
tar de la forma (y + xy + y)dx + (x + x y + x)dy = 0. No es exacta, notemos
que los casos 1 y 2 no son aplicables. Ahora calculemos el lado derecho de 1.44
obtenemos ∂M ∂N
∂y
− ∂x
N y−M x
(3y 2 +2xy+1)−(3x2 +2xy+1)

= (x3 +x2 y+x)y−(y 3 +xy 2 +y)x
3(y 2 −x2 )
= x3 y−y 3 x
x −y 2 2
= −3 xy(x 2 −y 2 )
−3
= xy
−3
Ası́, h(u) = u
que nos pone en el caso 3. Ası́ el factor integrante está dado
por
−3du
R
P =e u
du
R
= e−3 u
= e−3 ln u
−3
= eln u
= u−3 ,
P = x31y3 . Multiplicando nuestra ecuación original por el factor integrante te-
nemos la siguiente ecuaciṕon exacta.
( x31y3 )((y 3 + xy 2 + y)dx + (x3 + x2 y + x)dy)
= ( x13 + 1
x2 y
+ 1
x3 y 2
)dx + ( y13 + 1
xy 2
+ 1
x2 y 3
)dy
Siguiendo con el método tenemos lo siguiente.
a) Integrando M con respecto a x tenemos:

Z
1 1 1 −3 −2 −3
F (x, y) = 3
+ 2 + 3 2 dx + f (y) = 4 + 3 + 4 2 + f (y).
x xy xy x xy xy
b) Derivando F con respecto a y tenemos:
∂F (x, y) 2 6
= 3 2 + 4 3 + f 0 (y).
∂y xy xy
c) Igualando a N , tenemos
2 6 1 1 1
+ + f 0 (y) = + 2+ 2 3
x3 y 2 x4 y 3 y 3 xy xy
ası́
0 1 1 6 1 1 2
f (y) = 3 1+ 2 − 4 + 2 −
y x x y x x3
d) Calculando f tenemos

−3 1 6 −2 1 2
f (y) = 4 1+ 2 − 4 + 3 −
y x x y x x3
e) Sustituyendo f en F tenemos

−3 −2 −3 −3 1 6 −2 1 2
F (x, y) = 4 + 3 + 4 2 + 4 1+ 2 − 4 + 3 −
x xy xy y x x y x x3
Capı́tulo 2
Ec Dif Ord de Orden n
2.1. Introducción
La ecuación diferencial lineal de orden n con coeficientes constantes es una ecua-
ción de la forma
ao y (n) + a1 y (n−1) + . . . + an y = b(x),
donde a0 6= 0, a1 , . . . , an son constantes complejas y b una función de variable com-
pleja en un intervalo I. Al dividir la ecuación entre a0 obtenemos una ecuación de la
misma forma pero con coeficiente a0 = 1. Ası́ que siempre supondremos que a0 = 1,
por lo que tenemos la ecuación
y (n) + a1 y (n−1) + . . . + an y = b(x). (2.1)

Representaremos el lado izquierdo de 2.1 por L(y). Tenemos la siguiente defini-
ción.
Definición 4 Consideremos la ecuación dada por 2.1 L(y) = b(x). Si b(x) = 0

∀x ∈ I la ecuación será L(y) = 0 se llama ecuación diferencial lineal ho-
mogénea con coeficientes constantes. Si b(x) 6= 0, L(y) = b(x) se llama ecua-
ción diferencial lineal no homogénea con coeficientes constantes.
Nota: L es un operador diferencial que opera sobre funciones que tienen n deri-
vadas en I, y transforma dichas funciones φ en funciones L(φ) cuyo valor en x está
dado por
L(φ)(x) = φ(n) (x) + a1 φ(n−1) (x) + . . . + an φ(x).
29
30
2.2. Ecuación diferencial lineal homogénea con coeficientes constantes de segundo orden
Ası́
L(φ) = φ(n) + a1 φ(n−1) + . . . + an φ.
Por consiguiente una solución de L(y) = b(x) es una función φ que tiene n
derivadas en I, tales que L(φ) = b.
Si b es continua en I es posible obtener todas las soluciones de L(Y ) = b(x).
Tenemos la siguiente definición.
Definición 5 Consideremos la ecuación diferencial
L(y) = y (n) + a1 y (n−1) + . . . + an y = 0,

hagamos p(r) un polinomio dado por
p(r) = rn + a1 rn−1 + . . . + an r,
p(r) se conoce como el polinomio caracteristico de L.
Note que p(r) se obtiene a partir de L(y) sustituyendo y (k) por rk .

Para resolver 2.1 empezaremos primero con el caso de la ecuación de segundo
orden.
2.2. Ecuación diferencial lineal homogénea con co-

eficientes constantes de segundo orden
Estudiaremos la ecuación
L(y) = y 00 + a1 y 0 + a2 y = 0. (2.2)
donde a1 y a2 son constantes. Note que dada la ecuación 2.2 el polinomio ca-
racterı́stico asociado a ella p(r) tiene sus dos raı́ces complejas (esto lo garantiza el
teorema fundamental del álgebra.)
Teorema 7 Consideremos la ecuación 2.2, con a1 y a2 constantes. Si r1 y r2 son dos

raı́ces diferentes del polinomio caracterı́stico asociado a 2.2, entonces las funciones
φ1 y φ2 definidas por
φ1 (X) = er1 x , φ2 (X) = er2 x (2.3)

son soluciones de L(y) = 0.
Capı́tulo 2. Ec Dif Ord de Orden n 31
Si r es una raiz doble del polinomio caracterı́stico asociado a 2.2, entonces las
funciones φ1 y φ2 definidas por
φ1 (X) = er1 x , φ2 (X) = xer2 x (2.4)

Demostración
Consideremos la ecuación 2.2 y su polinomio caracterı́stico asociado p(r) Sean r
raiz de p(r). Tenemos
L(erx ) = (r2 + a1 r + a2 )erx , (2.5)
entonces L(erx ) = 0 ((erx ) = p(r)erx = 0erx ).
Consideremos el caso ahora de r una raiz doble. Tenemos entonces que p(r) = 0
y además p0 (r) = 0 y
L(xerx ) = [p0 (r) + xp(r)]erx = 0
Con lo que se tiene lo deseado.

NOTA: Si φ es una combinación lineal de φ1 y φ2 entonces φ también seŕa
solución. En efecto, sea φ = c1 φ1 + c2 φ2 tenemos
L(φ)=(c1 φ1 + c2 φ2 )00 + a1 (c1 φ1 + c2 φ2 )0 + a2 (c1 φ1 + c2 φ2 )

=c1 φ001 + c2 φ002 + a1 c1 φ01 + a1 c2 φ02 + a2 c1 φ1 + a2 c2 φ2
=c1 L(φ1 ) + c2 L(φ2 )
=0.
EJEMPLOS
2.3. Ecuación diferencial lineal homogénea con co-

eficientes constantes de segundo orden con
condiciones iniciales
Una ecuación con condiciones iniciales del tipo L(y) = 0 es el problema de
encontrar una función φ que satisfaga las condiciones:
φ(x0 ) = α, φ0 (x0 ) = β, (2.6)

donde x0 es algún número real y α y β son dos constantes dadas. Se representa
L(y) = 0, y(x0 ) = α, y 0 (x0 ) = β, (2.7)

32 con condiciones iniciales
Teorema 8 (Teorema de existencia.) Para cualquier x0 ∈ R y para las constan-

tes α y β existe una solución φ de la ecuación con condiciones iniciales 2.7 ∀x ∈ R
Demostración
Se demostrá que existen constante únicas c1 , c2 , tales que φ = c1 φ1 +c2 φ2 satisface
el sistema 2.6, donde φ01 , φ02 son las soluciones dadas por 2.3 o 2.4. Entonces se deben
satisfacer las siguientes condiciones:
c1 φ1 (x0 ) + c2 φ2 (x0 ) = α
c1 φ01 (x0 ) + c2 φ02 (x0 ) = β

la solución c1 , c2 es única si y sólo si el determinante cumple la condición
∆ = φ1 (x0 )φ02 (x0 ) − φ01 (x0 )φ2 (x0 ) 6= 0.
En el caso r1 6= r2 , tenemos 2.3 y
∆ = r2 er1 x0 er2 x0 − r1 er1 x0 er2 x0 = (r2 − r1 )e(r1 +r2 )x0 6= 0.
Si r1 es raiz doble entonces tenemos 2.4 y
∆ = er1 x0 (er1 x0 + x0 r1 er1 x0 ) − r1 x0 er1 x0 er1 x0 = e2r1 x0 6= 0.
Se tiene lo pedido.

Necesitamos demostrar el siguiente teorema para demostrar la unicidad de la
solución.
Teorema 9 Sea φ una solución cualquiera de
L(y) = y 00 + a1 y 0 + a2 y = 0
en un intervalo I tal que x0 ∈ I. Entonces para toda x ∈ I:
||φ(x0 )||e−k|x−x0 | ≤ ||φ(x)|| ≤ ||φ(x0 )||ek|x−x0 | (2.8)

donde
1
||φ(x0 )|| = [|φ(x0 )|2 + |φ0 (x0 )|2 ] 2 , k = 1 + |a1 | + |a2 |
Demostración
Sea u = ||φ(x0 )||2 , luego
u = φφ + φ0 φ0 ,
entonces
u0 = φ0 φ + φφ0 + φ00 φ0 + φ0 φ00 ,
por consiguiente tenemos
|u0 (x)| ≤ 2|φ(x)||φ0 (x)| + 2|φ0 (x)||φ00 (x)|. (2.9)
Dado que φ es solución, tenemos:
φ00 = −a1 φ0 − a2 φ,
luego por la desigualdad del triángulo tenemos
|φ00 (x)| ≤ |a1 ||φ0 (x)| + |a2 ||φ(x)|. (2.10)
Sustitutendo 2.10 en 2.9, obtenemos:
|u0 (x)| ≤ 2|φ(x)||φ0 (x)| + 2|φ0 (x)||φ00 (x)|

≤ 2|φ(x)||φ0 (x)| + 2|φ0 (x)|(|a1 ||φ0 (x)| + |a2 ||φ(x)|)
≤ 2(1 + |a2 |)|φ(x)||φ0 (x)| + 2|a1 ||φ0 (x)|2 .
Ahora notemos que
0 ≤ (|φ(x)| − |φ0 (x)|)2 = |φ(x)|2 + |φ0 (x)|2 − 2|φ(x)||φ0 (x)|,

ası́
2|φ(x)||φ0 (x)| ≤ |φ(x)|2 + |φ0 (x)|2
por lo que obtenemos:
|u0 (x)| ≤ (1 + |a2 |)(|φ(x)|2 + |φ0 (x)|2 ) + 2(|a1 ||φ0 (x)|2 )

≤ (1 + |a2 |)|φ(x)|2 + (1 + 2|a1 | + |a2 |)|φ0 (x)|2
≤ 2(1 + |a1 | + |a2 |)(|φ(x)|2 + |φ0 (x)|2 ),
o
|u0 (x)| ≤ 2ku(x).
Esto es equivalente a
− 2ku(x) ≤ u0 (x) ≤ 2ku(x). (2.11)

34 con condiciones iniciales
Considerando solamente la parte derecha de la desigualdad tenemos
u0 − 2ku ≤ 0.
Al multiplicar por e−2kx nos da
e−2kx (u0 − 2ku) = (e−2kx u)0 ≤ 0.

Consideremos x > x0 e integramos ambos lados de la desigualdad desde x0 a x
tenemos
e−2kx u(x) − e−2kx0 u(x0 ) ≤ 0,
o
u(x) ≤ u(x0 )e2k(x−x0 )
obteniendo la raiz cuadrada equivale a
||φ(x)|| ≤ ||φ(x0 )||ek(x−x0 ) , (x > x0 ),
Considerando ahora la parte izquierda de la desigualdad tenemos
u0 + 2ku ≥ 0.
Al multiplicar por e2kx obtenemos
e2kx (u0 + 2ku) = (e2kx u)0 ≥ 0.

Consideremos ahora cuando x > x0 integrando desde x0 a x tenemos
e2kx u(x) − e2kx0 u(x0 ) ≥ 0,
o
u(x) ≥ u(x0 )e−2k(x−x0 )
que equivale a
||φ(x)|| ≥ ||φ(x0 )||e−k(x−x0 ) , (x > x0 ),
por lo que tenemos
||φ(x0 )||e−k(x−x0 ) ≤ ||φ(x)|| ≤ ||φ(x0 )||ek(x−x0 ) (x > x0 ),

que es lo que se pedı́a.

Teorema 10 (Teorema de Unicidad) Sean α y β dos constantes cualquiera fijas,

y sea x0 ∈ R. En cualquier intervalo I tal que xo ∈ I existe a lo más una solución φ
de la ecuación diferencial con condiciones iniciales
L(y) = 0, y(x0 ) = α, y 0 (x0 ) = β.
Demostración
Supongamos que φ y φ1 son dos soluciones. Sea χ = φ − φ1 . Entonces se tiene
L(χ) = L(φ) − L(φ1 ) = 0, y χ(x0 ) = 0, χ0 (x0 ) = 0.
Ası́ ||χ(x0 )|| = 0, y aplicando las desigualdades de 2.8, vemos que
||χ(x)|| = 0 ∀ x ∈ I.
Por lo que χ(x) = 0 para toda x ∈ I, es decir φ = φ1 .

Ya estamos en posición de probar el siguiente teorema.
Teorema 11 Sean φ1 y φ2 , las dos soluciones de la ecuación L(y) = 0, dadas por

el Teorema 7. Si c1 y c2 son dos constantes cualquiera, la función φ = c1 φ1 + c2 φ2 es
una solución de L(y) = 0 en R.
Recı́procamente, si φ es una solución cualquiera de L(y) = 0 en R, existen dos
constantes únicas tales que
φ = c1 φ1 + c2 φ2 .
Demostración
La primera parte del teorema es consecuencia directa de que L es una transfor-
mación lineal.
Recı́procamente sea φ una solución y sea x0 ∈ R, sea φ(x0 ) = α y φ0 (x0 ) = β.
En el Teorema 8 se probo que existe una solución ψ de L(y) = 0, que satisface las
condiciones ψ(x0 ) = α, ψ 0 (x0 ) = β, y que tiene la forma :
ψ = c1 φ1 + c2 φ2 ,
donde c1 , c2 están unı́vocamente determinadas por α, β. Por el Teorema 10 φ = ψ.

36 2.4. Dependencia e independencia lineal
2.4. Dependencia e independencia lineal

Definición 6 Dadas dos funciones φ1 , φ2 definidas en un intervalo I, se dice que
son Linealmente Dependientes en I, si existen dos constantes c1 , c2 no ambas
nulas tales que se cumpla:
c1 φ1 + c2 φ2 = 0
para toda x ∈ I.
Se dice que son Linealmente Independientes en I, si no son linealmente
dependientes, es decir, si las únicas dos constantes c1 , c2 tales que
c1 φ1 + c2 φ2 = 0
para toda x ∈ I. Son c1 = 0, c2 = 0
Las funciones definidas por 2.3 son l.i. en cualquier intervalo i. Para probarlo
supongamos que
c1 er1 x + c2 er2 x = 0 (2.12)
para toda x ∈ I. Multiplicando por e−r1 x , obtenemos
c1 + c2 e(r2 −r1 )x = 0
derivando tenemos
c2 (r2 − r1 )e(r2 −r1 )x = 0
por lo que c2 = 0. Luego de 2.12 se tiene que c1 = 0.
Ahora, las funciones definidas por 2.4 son l.i. en cualquier intervalo i. Para pro-
barlo supongamos que
c1 er1 x + c2 xer1 x = 0
para toda x ∈ I. Multiplicando por e−r1 x , obtenemos
c1 + c2 x = 0
derivando tenemos c2 = 0. Luego se tiene que c1 = 0.
Definición 7 Sean φ1 y φ2 dos funciones, soluciones de L(y) = 0. Se define el

Wronskiano de φ1 , φ2 como

φ1 φ 2

0
0
φ1 φ2
Es una función, y su valor correspondiente a x, se representa por W (φ1 , φ2 )(x).
El siguiente Teorema nos da una forma sencilla de saber si dos funciones son l.i.
Teorema 12 Dadas φ1 , φ2 dos soluciones de L(y) = 0 son l.i.en I, si y sólo si
W (φ1 , φ2 )(x) 6= 0
para toda x ∈ I.
Demostración
Primero supongamos que
W (φ1 , φ2 ) 6= 0 ∀x ∈ I,
y consideremos constantes c1 , c2 , tales que
c1 φ1 + c2 φ2 = 0 (2.13)
para toda x ∈ I.Entonces tenemos
c1 φ01 + c2 φ02 = 0 (2.14)
para toda x ∈ I. Para x fija, las ecuaciones 2.13, 2.14 son ecuaciones lineales ho-
mogeneas que se satisfacen para c1 , c2 . El determinante del sistema es precisamente
W (φ1 , φ2 )(x), el cual no es nulo, por lo cual c1 = c2 = 0 es la únoica solución del
sistema dado por 2.13, 2.14 por lo que φ1 , φ2 son l.i en I.
Recı́procamente supongamos que φ1 , φ2 son l.i en I. Supongamos que existe un
x0 ∈ I tal que W (φ1 , φ2 )(x0 ) = 0. Esto implica que el sistema de dos ecuaciones:
c1 φ1 + c2 φ2 = 0 (2.15)
c1 φ01 + c2 φ02 = 0, (2.16)
tiene una solución c1 , c2 , con al menos uno distinto de cero. Sea c1 , c2 una solución
de este tipo, y consideremos la función ψ = c1 φ1 + c2 φ2 . Ahora L(ψ) = 0 por 2.16,
tenemos que ψ(x0 ) = 0 y ψ 0 (x0 ) = 0.
Luego por el Teorema 10 tenemos que ψ(x) = 0para toda x ∈ I, luego
c1 φ1 + c2 φ2 = 0
para toda x ∈ I, lo que contradice el hecho de que φ1 , φ2 son l.i. Por lo que la
hipótesis de que exista x0 ∈ I tal que W (φ1 , φ2 )(x0 ) = 0 es falsa, con lo que se ha
demostrado que
W (φ1 , φ2 )(x0 ) 6= 0

38 2.4. Dependencia e independencia lineal
Teorema 13 Dadas φ1 , φ2 dos soluciones de L(y) = 0, y sea x0 ∈ I. Entonces φ1 , φ2

son l.i.en I, si y sólo si
W (φ1 , φ2 )(x0 ) 6= 0
Demostración
Supongamos que φ1 , φ2 son l.i. en I entonces se cumple la condición para toda
x ∈ I, en particular para algún x0 ∈ I.
Recı́procamente, supongamos que W (φ1 , φ2 )(x0 ) 6= 0 y supongamos que c1 , c2 son
dos constantes tales que
c1 φ1 + c2 φ2 = 0
para toda x ∈ I, entonces tenemos
c1 φ1 + c2 φ2 = 0
c1 φ01 + c2 φ02 = 0,

Teorema 14 Dadas φ1 , φ2 dos soluciones cualesquiera de L(y) = 0 l.i. en I. Toda
solución de L(y) = 0 puede escribirse unı́vocamente ası́:
φ = c1 φ1 + c2 φ2
donde c1 , c2 son constantes.
Demostración
Sea x0 ∈ I. Sean φ(x0 ) = α y φ0 (x0 ) = β y consideremos las ecuaciones
c1 φ1 + c2 φ2 = αc1 φ01 + c2 φ02 = β,
para las constantes c1 , c2 . El determinate del sistema es W (φ1 , φ2 )(x0 ) 6= 0 (ya que
φ1 , φ2 son l.i. en I), por lo que existe constantes únicas c1 , c2 que satisfacen estas
ecuaciones. Sea ψ = c1 φ1 + c2 φ2 , tiene
ψ(x0 ) = φ(x0 ), ψ 0 (x0 ) = φ(x0 ), y L(ψ) = 0
luego por el Teorema 10 concluimos que ψ = φ en I.

La importancia del Teorema anterior radica en que solamente necesitamos conocer
dos soluciones cualesquiera de L(y) = 0, linealmente independientes para obtener
todas las soluciones de la ecuación. Por ejemplo, la ecuación
y 00 + y = 0
tiene las soluciones eix , e−ix , las cuales son linealmete independientes, pero también
tiene las dos soluciones linealmente independientes cosx, senx.
2.5. Una formula para el Wronskiano

En esta sección veremos una formula para calcular el Wronskiano de dos solucio-
nes de L(y) = 0 en cualquier punto, a partir de un punto x0 ∈ I dado.
Teorema 15 Si φ1 , φ2 son dos soluciones de L(Y ) = 0 en un intervalo I que con-

tenga x0 , entonces
W (φ1 , φ2 )(x) = e−a1 (x−x0 ) W (φ1 , φ2 )(x0 ) (2.17)
Demostración
Tenemos
φ001 + a1 φ01 + a2 φ1 = 0
φ002 + a1 φ02 + a2 φ2 = 0,
multiplicando la primera ecuación por −φ2 , la segunda por φ1 , y sumado obtenemos:
(φ1 φ002 − φ001 φ2 ) + a1 (φ1 φ02 − φ01 φ2 ) = 0.
Sea W = W (φ1 , φ2 ), observemos que
W = φ1 φ02 − φ01 φ2 , y W 0 = φ1 φ002 − φ001 φ2 .
ası́, W satisface la ecuación de primer orden
W 0 + a1 W = 0.
Por consiguiente
W (x) = ce−a1 x
para alguna constante c. Tomando x = x0 tenemos que
c = ea1 x0 W (x0 ),
luego
W (x) = e−a1 (x−x0 ) W (x0 )

40 2.6. La ecuación diferencial no homogénea de segundo orden
2.6. La ecuación diferencial no homogénea de se-

gundo orden
Hallaremos las soluciones de la ecuación
L(y) = y 00 + a1 y 0 + a2 y = b(x),
donde b es una funcion continua en el intervalo I.

Tenemos el siguiente Teorema.
Teorema 16 Sea b una función continua en el intervalo I. toda solución ψ de la

ecuación L(y) = b(x) en I se puede escribir de la forma:
ψ = ψp + c1 φ1 + c2 φ2
donde ψp es una solución particular, φ1 , φ2 son dos soluciones de L(y) = 0, l.i.y c1 , c2

son constantes. Recı́procamente toda ψ de ese tipo es una solución de la ecuación
L(y) = b(x).
Demostración
Encontremos primero una solución particular de la ecuación diferencial L(y) =
b(x). Para ello utilizaremos el método de variación de parametros el cual consiste en
suponer que existen dos funciones µ1 , µ2 tales que µ1 φ1 + µ2 φ2 es una solución en I.
Entonces tenemos
(µ1 φ1 + µ2 φ2 )00 + (µ1 φ1 + µ2 φ2 )0 + a2 (µ1 φ1 + µ2 φ2 )

= µ1 L(φ1 ) + µ2 L(φ2 ) + (µ001 φ1 + µ002 φ2 ) + 2(µ01 φ01 + µ02 φ02 ) + a1 (µ01 φ1 + µ02 φ2 )
(µ001 φ1 + µ002 φ2 ) + 2(µ01 φ01 + µ02 φ02 ) + a1 (µ01 φ1 + µ02 φ2 ) = b
observemos que si
µ01 φ1 + µ02 φ2 = 0 (2.18)
entonces
0 = (µ01 φ1 + µ02 φ2 )0 = (µ001 φ1 + µ002 φ2 ) + (µ01 φ01 + µ02 φ02 ),
por consiguiente
µ01 φ01 + µ02 φ02 = b (2.19)
Ası́ se tiene que si podemos hallar dos funciones µ1 , µ2 tales que satisfagan las
condiciones 2.18 y 2.19, entonces la función µ1 φ1 + µ2 φ2 satisfacerá L(y) = b(x).
Las ecuaciones 2.18 y 2.19 son dos ecuaciones lineales en µ01 , µ02 , cuyo determinante
es precisamente W (φ1 , φ2 ), dado que hicimos la hipótesis de que son l.i. entonces el
determinate nunca se anulará, por lo que existe una solución única de µ01 , µ02 , que
está dada por:
−φ2 b φ1 b
µ01 = , µ02 = .
W (φ1 , φ2 ) W (φ1 , φ2 )
para obtener µ1 , µ2 , basta integrar (se puede integrar ya que ). Tomando x0 ∈ I
podemos considerar µ1 , µ2 como las siguientes funciones:
Z x Z x
φ2 (t)b(t) φ1 (t)b(t)
µ1 = − dt, µ2 = dt
x0 W (φ1 , φ2 )(t) x0 W (φ1 , φ2 )(t)
Por consiguiente la solución particular ψp = µ1 φ1 + µ2 φ2 será

x
(φ1 (t)φ2 (x) − φ1 (x)φ2 (t)) b(t)
Z
ψp = dt. (2.20)
x0 W (φ1 , φ2 )(t)
Ahora sea que ψp nuestra una solución particiular de esta ecuación dada por 2.20,
y que ψ es una solución. Entonces
L(ψ − ψp ) = L(ψ) − L(ψp ) = b − b = 0
en I. Esto demuestra que ψ − ψp es una solución de la ecuacón homogénea L(y) = 0.

Por consiguiente si φ1 , φ2 son soluciones de L(y) = 0. l.i., existen las constantes
únicas c1 , c2 tales que
ψ − ψp = c1 φ1 + c2 φ2
Es decir cualquier solución ψ de L(y) = b(x) puede escribirse de la forma:
ψ = ψp + c1 φ1 + c2 φ2

EJEMPLOS
1. Sea y 00 − y 0 − 2y = e−x
2. Sea 4y 00 − y = ex
42 2.7. Ecuación diferencial homogenea de orden n
2.7. Ecuación diferencial homogenea de orden n

En esta sección generalizaremos los resultados obtenidos en la seccion anterior
para la ecuación diferencial de orden 2 ahora para la de orden n. Supongamos que
L(y) esta dada por
L(y) = y (n) + a1 y (n−1) + . . . + an y,

donde a1 , a2 , . . . , an son constantes. Resolveremos primero el problema L(y) = 0.
Teorema 17 Sea L(y) = 0 una ecuación diferencial ordinaria homogénea con coefi-
cientes constantes de orden n, sea p(r) el polinomio caracterı́stico y sean r1 , r2 , . . . , rs
las raices con multiplicidades mi i = 1, . . . , s, entonces las n funciones
er1 x , xer1 x , . . . , xm1 −1 er1 x ;
er2 x , xer2 x , . . . , xm2 −1 er2 x ; · · ·

ers x , xers x , . . . , xms −1 ers x
Demostración
Si r es una raiz del polinomio, claramente L(er1 x ) = 0, por lo que erx es solución.
Si r1 es raiz de multiplicidad m1 del polinomio, entonces
p(r1 ) = 0, p0 (r1 ) = 0, · · · , pm−1 (r1 ) = 0
Si derivamos la ecuación L(erx ) = p(r)erx con respecto a r obtenemos
∂k
k
rx ∂ rx
L(e ) = L e = L(xk erx )
∂rk ∂rk

(k) (k−1) k(k − 1) (k−2) 2
= p (r) + kp (r)x + p x + · · · + p(r)x erx .
k
2!
Ası́ para k = 0, 1, 2, · · · , m − 1 tenemos que xk erx es una solución de L(y) = 0.

Teorema 18 Sea L(y) = 0 una ecuación diferencial ordinaria homogénea con co-
eficientes constantes de orden n,dada en el teorema 17 las n soluciones de L(y) = 0
son l.i.en cualquier intervalo.
Demostración
Suponga que tenemos cij con i = 1, · · · , s, j = 0, · · · , mi − 1 tal que
s m
X Xi −1
cij xj = 0 (2.21)
i=0 j=0
en I. Sumando sobre j para i fijo, ponemos

m
X i −1
Pi (x) = cij xj
j=0
como el polinimio caracterı́stico de eri x en 2.21. Entonces tenemos

P1 (x)er1 x + P2 (x)er2 x + · · · + Ps (x)ers x = 0 (2.22)
en I.Asumamos que no todas las constantescij son cero. Entonces tendremos al menos
un polinomio Pi que sea distinto de cero. Podemos asumir que es Ps , entonces 2.22
implica
P1 (x) + P2 (x)e(r2 −r1 )x + · · · + Ps (x)e(rs −r1 )x = 0 (2.23)
en I. Derivando suficiente número de veces P1 lo podemos reducir a cero. En este
proceso el grado de los polinomios que multiplican a e(ri −r1 )x no cambian, como
tampoco cambia el hecho de que sean o no identicamente cero. Obtenemos una
expresión de la forma
Q2 (x)e(r2 −r1 )x + Q3 (x)e(r3 −r1 )x + · · · + Qs (x)e(rs −r1 )x = 0
o
Q2 (x)er2 x + Q3 (x)er3 x + · · · + Qs (x)ers x = 0
en I, donde cada Qi es polinomio, con gradPi = gradQi , y Qs no nulo. Continuando
el proceso llegaremos a
Rs (x)ers x = 0 (2.24)
en I, donde Rs es polinomio, con gradPs = gradQs , y Qs no nulo. Pero implica que
Rs (x) = 0 ∀x ∈ I, lo cual es una contradiccón, por lo que la hipotésis de que al
menos uno de los cij sea distinto de cero. Por lo que las soluciones son l.i.

Recordemos que φ1 , · · · , φm son soluciones de L(y) = 0 en un intervalo, entonces
φ = c1 φ1 + · · · + cm φm
También es solución de L(y) = 0 en el intervalo.
44 2.8. Ecuaciones de orden n con condiciones iniciales
Ejemplos 5
y 000 − 3y 0 + 2y = 0
Tiene polinimio caracteristico
p(r) = r3 − 3r + 2
cuyas raices son 1, 1, −2 Ası́ tenemos tres soluciones linealmente independientes
ex , xex , e−2x
y cualquier solución tiene la forma
φ(x) = (c1 + xc2 )ex + c3 e−2x
dode c1 , c2 , c3 son cualquier constante.
2.8. Ecuaciones de orden n con condiciones inicia-

les
El problema de condiciones iniciales para L(y) = 0, es encontrar una solución φ
que satisfaga
φ(x0 ) = α1 , φ0 (x0 ) = α2 , · · · , φ(n−1) (x0 ) = αn ,
donde x0 es un número real y α1 , · · · , αn son constantes dadas. Se representa ası́:
L(y) = 0, , y(x0 ) = α1 , y (n−1) (x0 ) = αn .
Existe solamente una solución para dicho problema, la demostración se basa en

la estimación de la taza de crecimiento de φ, para empezar tenemos la siguiente
definición.
Definición 8 Dada una ecuación L(y) = 0 y una solución φ que tenga las derivadas
φ0 , φ00 , · · · , φn . Definimos ||φ(x)|| como
1/2
||φ(x)|| = |φ(x)|2 + · · · + |φ(n−1) (x)|2 .
L(y) = y (n) + a1 y (n−1) + · · · + an y = 0

en un intervalo I tal que x0 ∈ I. Entonces para toda x ∈ I:
||φ(x0 )||e−k|x−x0 | ≤ ||φ(x)|| ≤ ||φ(x0 )||ek|x−x0 | (2.25)

donde
k = 1 + |a1 | + |a2 | + · · · + |an |.
Demostración
Sea u = ||φ(x0 )||2 , luego
u = φφ + φ0 φ0 + · · · + φ(n−1) φ(n−1) .
Entonces
u0 = φ0 φ + φ00 φ0 + · · · + φ(n) φ(n−1) + φφ0 + φ00 φ0 + · · · + φ(n−1) φ(n)
por consiguiente tenemos
|u0 (x)| ≤ 2|φ(x)||φ0 (x)| + 2|φ0 (x)||φ00 (x)| + · · · + 2|φ(n−1) (x)||φ(n) (x)|. (2.26)
Dado que φ es solución, tenemos:
φ(n) = − a1 φ(n−1) + · · · + an φ ,

luego
|φ(n) (x)| ≤ |a1 ||φ(n−1) (x)| + · · · + |an ||φ(x)|. (2.27)
Sustitutendo 2.27 en 2.26, obtenemos:
|u0 (x)| ≤ 2|φ(x)||φ0 (x)| + 2|φ0 (x)||φ00 (x)| + · · ·

+2|φ(n−2) (x)||φ(n−1) (x)| + 2|a1 ||φ(n−1) (x)|2 +
2|a2 ||φ(n−2) (x)||φ(n−1) (x)| + · · · + 2|an ||φ(x)||φ(n−1) (x)|.
Aplicando la propiedad
2|b||c| ≤ |b|2 + |c|2 ,

obtenemos:
|u0 (x)| ≤ (1 + |an |)|φ(x)|2 +

(2 + |an−1 |)|φ0 (x)|2 + · · · + (2 + |a2 |)|φ(n−2) (x)|2
+(1 + 2|a1 | + |a2 | + · · · + |an |)|φ(n−1) (x)|2
o
|u0 (x)| ≤ 2ku(x).
Esto es equivalente a
− 2ku(x) ≤ u0 (x) ≤ 2ku(x). (2.28)

Considerando solamente la parte derecha de la desigualdad tenemos
u0 − 2ku ≤ 0.
Entonces tenemos una ecuación diferencial de primer orden en u, la que al multiplicar
por e−2kx nos da
e−2kx (u0 − 2ku) = (e−2kx u)0 ≤ 0.

e−2kx u(x) − e−2kx0 u(x0 ) ≤ 0,
o
u(x) ≤ u(x0 )e2k(x−x0 )
elevando al cuadrado equivale a
||φ(x)|| ≤ ||φ(x0 )||ek(x−x0 ) , (x > x0 ),
Considerando ahora la parte izquierda de la desigualdad tenemos
u0 − 2ku ≥ 0.
Entonces tenemos una ecuación diferencial de primer orden en u, la que al multiplicar
por e−2kx obtenemos
e−2kx (u0 − 2ku) = (e−2kx u)0 ≥ 0.

e−2kx u(x) − e−2kx0 u(x0 ) ≥ 0,

o
u(x) ≥ u(x0 )e−2k(x−x0 )
que equivale a
||φ(x)|| ≥ ||φ(x0 )||ek(x−x0 ) , (x > x0 ),
por lo que tenemos
||φ(x0 )||e−k(x−x0 ) ≤ ||φ(x)|| ≤ ||φ(x0 )||ek(x−x0 ) (x > x0 ),

Considerando 2.28 cuando x < x0 e integrando de x a x0 tenemos
||φ(x0 )||e−k(x−x0 ) ≤ ||φ(x)|| ≤ ||φ(x0 )||ek(x−x0 ) (x0 > x),

que es lo que se pedı́a.

Teorema 20 (Teorema de Unicidad) Sean α1 , · · · , αn n constantes cualquiesra
fijas, y sea x0 ∈ R. En cualquier intervalo I tal que xo ∈ I existe a lo más una
solución φ de la ecuación diferencial con condiciones iniciales
L(y) = 0, φ(x0 ) = α1 , · · · φ(n − 1)(x0 ) = αn .
Demostración
Supongamos que φ1 y φ2 son dos soluciones. Sea χ = φ1 − φ2 . Entonces se tiene
L(χ) = L(φ1 ) − L(φ2 ) = 0, y χ(x0 ) = 0, χ0 (x0 ) = 0.
Ası́ ||χ(x0 )|| = 0, y aplicando las desigualdades de 2.10, vemos que
||χ(x)|| = 0 ∀ x ∈ I.
Por lo que χ(x) = 0 para toda x ∈ I, es decir φ1 = φ2 .

Definición 9 Dadas n funciones φ1 , · · · , φn definidas en un intervalo I, se dice que
son Linealmente dependientes en I, si existen n constantes c1 , · · · , cn no todas
nulas tales que se cumpla:
c1 φ1 + · · · · · · + cn φn = 0
para toda x ∈ I.
Se dice que son Linealmente independientes en I, si no son linealmente
dependientes, es decir, si las únicas n constantes c1 , · · · , cn tales que
c1 φ1 + · · · + cn φn = 0
para toda x ∈ I. Son c1 = 0, · · · , cn = 0
Definición 10 Sean φ1 , · · · , φn funciones, soluciones de L(y) = 0. Se define el

Wronskiano de φ1 , · · · , φn como

φ1 φ2
· · · φn
φ0 φ0 · · · φ0n
1 2
.. ..
n−1. n−1
. · · ·
φ1 φ2 · · · φnn−1
Es una función, y su valor correspondiente a x, se representa por W (φ1 , · · · , φn )(x).
Teorema 21 Dadas φ1 , · · · , φn n soluciones de L(y) = 0 son l.i.en I, si y sólo si
W (φ1 , · · · , φn )(x) 6= 0
para toda x ∈ I.
Demostración
W (φ1 , · · · , φn ) 6== 0 ∀x ∈ I,
y consideremos constantes c1 , · · · , cn , tales que
c1 φ1 + · · · + cn φn = 0 (2.29)
para toda x ∈ I.Entonces tenemos
c1 φ01 + · · · + cn φ0n = 0,
..
. (2.30)
(n−1)
c1 φ1 + · · · + cn φ(n−1)
n =0
para toda x ∈ I. Para x fija, las ecuaciones 2.29, 2.30 son ecuaciones lineales
homogeneas que se satisfacen para c1 , · · · , cn . El determinante del sistema es preci-
samente W (φ1 , · · · , φn )(x), el cual no es nulo, por lo cual c1 = · · · = cn = 0 es la
única solución del sistema dado por 2.29, 2.30 por lo que φ1 , · · · , φn son l.i en I.
Recı́procamente supongamos que φ1 , · · · φn son l.i en I. Supongamos que existe un
x0 ∈ I tal que W (φ1 , · · · , φn )(x0 ) = 0. Esto implica que el sistema de n ecuaciones:
c1 φ1 + · · · + cn φn = 0
c1 φ01 + · · · + cn φ0n = 0,
.. (2.31)
.
(n−1)
c1 φ1 + · · · + cn φ(n−1)
n =0
tiene una solución c1 , · · · , cn , con al menos uno distinto de cero. Sea c1 , · · · , cn

una solución de este tipo, y consideremos la función ψ = c1 φ1 + · · · + cn φn . Ahora
L(ψ) = 0 por 2.31, tenemos que ψ(x0 ) = 0, · · · , ψ (n−1) (x0 ) = 0.
Luego por el Teorema 10 tenemos que ψ(x) = 0para toda x ∈ I, luego
c1 φ1 + · · · + cn φn = 0
para toda x ∈ I, lo que contradice el hecho de que φ1 , · · · , φn son l.i. Por lo que
la hipótesis de que exista x0 ∈ I tal que W (φ1 , · · · , φn )(x0 ) = 0 es falsa, con lo que
se ha demostrado que
W (φ1 , · · · , φn )(x0 ) 6= 0

Teorema 22 (Teorema de existencia.) Para cualquier x0 ∈ R y para las cons-

tantes α1 , · · · , αn existe una solución φ de la ecuación con condiciones iniciales
L(y) = 0, φ(x0 ) = α1 , · · · φ(n − 1)(x0 ) = αn .
∀x ∈ R.
Demostración
Se demostrá que existen constante únicas c1 , · · · , cn , tales que φ = c1 φ1 + · · · +
cn φn satisface las condiciones iniciales. Entonces se deben satisfacer las siguientes
condiciones:
c1 φ1 (x0 ) + · · · + cn φn (x0 ) = α1
c1 φ01 (x0 ) + · · · + cn φ0n (x0 ) = α2 ,
.. (2.32)
.
(n−1)
c1 φ1 (x0 ) + · · · + cn φ(n−1)
n (x0 ) = αn
50 2.9. Ecuaciones cuyas constantes son reales
la solución c1 , · · · , cn es única si y sólo si el determinante es distinto de cero,

dicho determinante es precisamente el Wronskiano, el cual es diferente de cero por
el Teorema 21 se tiene lo pedido.

Veremos una formula para calcular el Wronskiano de n soluciones de L(y) = 0
en cualquier punto, a partir de un punto x0 ∈ I dado.
Teorema 23 Si φ1 , · · · , φn son n soluciones de L(Y ) = 0 en un intervalo I que

contenga x0 , entonces
W (φ1 , · · · , φn )(x) = e−a1 (x−x0 ) W (φ1 , · · · , φn )(x0 ) (2.33)
Omitiremos la demostración por que se verá como corolario del teorema que se
refiere a n soluciones con coeficientes variables.
Corolario 1 Dadas φ1 , · · · , φn n soluciones de L(y) = 0, y sea x0 ∈ I. Entonces

φ1 , · · · , φn son l.i.en I, si y sólo si
W (φ1 , · · · , φn )(x0 ) 6= 0
La demostración es consecuencia inmediata del Teorema 21 y de la formula 2.33.
y 000 − 3r1 y 00 + 3r12 y 0 − r13 y
cuyo polinomio caracteristico es p(r) = (r − r1 )3
2. y 000 − 4y 0 = 0
2.9. Ecuaciones cuyas constantes son reales

Teorema 24 Supongamos que todas las constantes a1 , a2 , · · · , an de la ecuación
L(y) = y (n) + a1 y (n−1) + · · · + an y = 0
son reales. Entonces existe un conjunto de n soluciones de variable real que son
linealmente independientes y toda solución de variable real, es una combinación lineal
de ellas con coeficientes reales. Si una solución satisface condiciones iniciales reales,
es de variables reales.
Demostración
Recordemos que si tenemos un polinomio con coeficientes reales entonces si tie-
ne una raı́z compleja su conjugado tambien será raı́z, y más aun tendrá la misma
multiplicidad.

La importancia del Teorema anterior radica en que en muchos problemas prácti-
cos, surgen ecuaciones diferenciales con coeficientes reales, y las soluciones reales son
evidentes. Por ejemplo
y (4) + y = 0
surge en el estudio de deflexiones de vigas. El polinomio caracteristico está dado pòr
p(r) = r4 + 1
y sus raices son:

1 1 1 1
√ (1 + i), √ (1 − i), √ (−1 + i), √ (−1 − i).
2 2 2 2
Ası́ toda solución real de la ecuación tiene la forma
√ √ √ √ √ √
φ = ex/ 2 c1 cos(x/ 2) + c2 sen(x/ 2) + e−x/ 2 c3 cos(x/ 2) + c4 sen(x/ 2)
donde c1 , c2 , c3 , c4 son constantes reales.
2.10. La ecuación no homogenea de orden n.

Teorema 25 Sea b continua en un intervalo I, y sean φ1 , φ2 , · · · , φn n soluciones
linealmente independientes de L(y) = 0 en I. Toda solución ψ de L(y) = b(x) puede
escribirse de la forma
ψ = ψp + c1 φ1 + · · · + cn φn
donde ψp es una solución particular de L(y) = b(x), y c1 , · · · , cn son constantes.
Recı́procamente ψ escrita de esta manera es solución.
Demostración
Encontremos primero una solución particular de la ecuación diferencial L(y) =
b(x). Para ello utilizaremos el método de Variación de Parametros el cual consiste
en suponer que existen dos funciones µ1 , · · · , µn tales que µ1 φ1 + · · · + µn φn es una
solución en I. Entonces tenemos
52 2.10. La ecuación no homogenea de orden n.
(µ1 φ1 + · · · + µn φn )(n) + a1 (µ1 φ1 + · · · + µn φn )(n−1) + · · · + an (µ1 φ1 + · · · + µn φn ) = b,
Si µ01 , · · · , µ0n satisfacen
µ01 φ1 + · · · + µ0n φn = 0
µ01 φ01 + · · · + µ0n φ0n = 0
..
. (2.34)
(n−2)
µ01 φ1 + · · · + µ0n φ(n−2)
n =0
(n−1)
µ01 φ1 + · · · + µ0n φ(n−1)
n =b
Vemos que
µ1 φ1 + · · · + µn φn = ψp
µ1 φ01 + · · · + µn φ0n = ψp0
..
. (2.35)
(n−1)
µ1 φ1 + · · · + µn φ(n−2)
n = ψp(n−1)
(n)
µ1 φ1 + · · · + µn φ(n)n + b = ψp
(n)
Por consiguiente:
L(ψp ) = µ1 L(φ1 ) + · · · + µn L(φn ) + b = b
y, por lo tanto ψp es una solución de L(y) = b(x). Todo el problema se reduce a

resolver el sistema de ecuaciones lineales (2.34) para µ01 , · · · , µ0n . El determinante
de los coeficientes es precisamente W (φ1 , · · · , φn ), el cual nunca se anula cuando
φ1 , · · · , φn son l.i.. Por consiguiente existe un conjunto único de funciones µ01 , · · · , µ0n ,
que satisface (2.34) y estan dadas por
Wk (x)b(x)
µ0k (x) = , (k = 1, 2, · · · , n),
W (φ1 , · · · φn )(x)
donde Wk es el determinante que se obtiene de W (φ1 , · · · , φn ) reemplazando la

k-esima columna por 0, 0, · · · , 0, b.

φ1 φ2
...0 ··· φn

φ0 φ0 . ..0 ··· φ0n
1 2
.. ..
n−1. n−1
. ···
φnn−1

φ1 φ2 ...b ···
Si x0 es cualquier punto de I podemos considerar que µk es la función dada por la
siguiente relación:
Z x
Wk (t)b(t)
µk (x) = dt, (k = 1, 2, · · · , n),
x0 W (φ1 , · · · φn )(t)
Por consiguiente la solución particular ψp = µ1 φ1 + · · · + µn φn será

n Z x
X Wk (t)b(t)
ψp = φk (x) dt. (2.36)
k=1 x0 W (φ1 , · · · , φn )(t)
Ahora sea que ψp nuestra una solución particiular de esta ecuación dada por 2.36,
y que ψ es una solución. Entonces
L(ψ − ψp ) = L(ψ) − L(ψp ) = b − b = 0
en I. Esto demuestra que ψ − ψp es una solución de la ecuacón homogénea L(y) = 0.

Por consiguiente si φ1 , · · · , φn son soluciones de L(y) = 0. l.i., existen las constantes
únicas c1 , · · · , cn tales que
ψ − ψp = c1 φ1 + · · · + cn φn
Es decir cualquier solución ψ de L(y) = b(x) puede escribirse de la forma:
ψ = ψp + c1 φ1 + · · · + cn φn

y 000 + y 00 + y 0 + y = 1
con condiciones iniciales ψ(0) = 0, ψ 0 (0) = 1, ψ 00 (0) = 0 cuyo polinomio

caracteristico es p(r) = r3 + r2 + r + 1 y tiene raices i, −i, −1
2. y (4) + 16y = cosx

54 2.11. Método de coeficientes indeterminados para resolver ecuaciones no homogéneas
2.11. Método de coeficientes indeterminados para

resolver ecuaciones no homogéneas
2.11.1. Álgebra de los operadores con coeficientes constantes
Trabajaremos con una ecuación del tipo L(y) = b(x) y donde toda solución ψ
tiene todas las derivadas en el intervalo (−∞, +∞).
Todos los operadores que se van a definir, se supone que están definidos en el con-
junto de todas las funciones φ en el intervalo (−∞, +∞). Sean L y M los operadores
dados por
L(φ) = a0 φ(n) + a1 φ(n−1) + · · · + an φ,
M (φ) = b0 φ(m) + b1 φ(m−1) + · · · + bm φ,
donde a0 , . . . , an , b0 , . . . , bm son constantes y a0 6= 0 6= b0 . Ası́ los polinomios carac-
terı́sticos de L y M son
p(r) = a0 rn + a1 rn−1 + · · · + an r
y
q(r) = b0 rm + b1 rm−1 + · · · + bm r
respectivamente.
Definimos la suma de L + M como el operador dado por
(L + M )(φ) = L(φ) + M (φ),
y el producto M L como el operador dado por
(M L)(φ) = M (L(φ)).
Si α es constante, definimos αL por
(αL)(φ) = αL(φ).
Observemos que L + M , M L y αL son operadores diferenciable lineales con
coeficientes constantes.
Si D es un operador de diferenciación,
D(φ) = φ0 ,
definimos D2 = DD, D3 = DDD y ası́ sucesivamente, definimos D0 (φ) = φ. Si α es
una constante, entonces
α(φ) = (αD0 )(φ) = αφ
Teorema 26 La correspondencia que asocia a cada
L = a0 Dn + a1 Dn−1 + · · · + an D0 ,
es el polinomio caracterı́stico p dado por
p(r) = a0 rn + a1 rn−1 + · · · + an ,
es una correspondencia biunı́voca entre todos los operadores diferenciales lineales

con coeficientes constantes y todo los polinomios. Si L y M están asociados con p, q
respectivamente, entonces L + M está asociado con p + q, M L con pq y αL está
asociado con αp.
Demostración
Se dice que dos operadores lineales son iguales L y M si
L(φ) = M (φ)
∀φ que tenga un número infinito de derivadas en el intervalo (−∞, +∞). Supongamos

que L, M tienen polinomios caracterı́sticos p, q respectivamente. Vemos que si L =
M , entonces
L(erx ) = p(r)erx = M (erx ) = q(r)erx ,
entonces p(r) = q(r)∀r.
Para el resto de la demostración observemos que
(L + M )(erx ) = L(erx ) + M (erx ) = [p(r) + q(r)]erx ,
(M L)(erx ) = M (L(erx )) = M (p(r)erx ) = p(r)M (erx ) = p(r)q(r)erx ,

(αL)(erx ) = α(L(erx )) = αp(r)erx .

El resultado implica que las propiedades algebraicas de los operadores con coefi-
cientes constantes son las mismas que las de los polinomios.
Teorema 27 Consideremos la ecuación con coeficientes constantes
L(y) = P (x)eax , (2.37)
Donde P es el polinomio dado por
P (x) = b0 xm + b1 xm−1 + · · · + bm , (b0 6= 0). (2.38)

Supongamos que a es una raı́z de multiplicidad j, del polinomio caracterı́stico p

de L. Entonces existe una solución única ψ de la ecuación 2.37 que tiene la forma
ψ(x) = xj (c0 xm + c1 xm−1 + · · · + cm )eax ,
donde c0 , c1 , . . . , cm , son constantes que se determinan por el método de los coeficien-

tes indeterminados.
2.11.2. Sobre el Método

El siguiente método es generalmente más rápido para la resolución de ecuaciones
no homogéneas el tipo L(y) = b(x), cuando b es una solución de una ecuación ho-
mogénea M (y) = 0 con coeficientes constantes. Entonces b debe de ser una suma de
tipo polinomial p(x)eax .
Supongamos que L y M tienen coeficientes constantes y que son de orden m y n,
respectivamente. Si ψ es una solución de L(y) = b, y M (b) = 0, entonces
M (L(ψ)) = M (b) = 0
Esto demuestra que ψ es solución de la ecuación homogénea M (L(y)) = 0. Ası́, ψ

puede escribirse como una combinación lineal de m + n soluciones de M (L(y)) = 0.
Sin embargo, no cualquier combinación lineal es solución de L(y) = b(x). Para hallar
cuáles condiciones deben satisfacer las constantes, volvemos a sustituir en L(y) = b.
Esto determina un conjunto de coeficientes.
Ejemplo: L(y) = y 00 − 3y 0 + 2y = x2 .
Como x2 es una solución de M (y) = y 000 = 0, vemos que toda solución ψ de
L(y) = x2 es una solución de
M (L(y)) = y v − 3y iv + 2y 000 = 0
El polinomio caracterı́stico es r3 (r2 −2r+2), cuyas raı́ces son 0, 1 y 2, por consiguiente

ψ debe tener la forma
ψ(x) = c0 + c1 x + c2 x2 + c3 ex + c4 e2x .
Observemos que c3 ex + c4 e2x es una solución de L(y) = 0. Estamos interesados

en calcular una solución particular ψp de L(y) = x2 , podemos suponer que tiene la
forma
ψp (x) = c0 + c1 x + c2 x2 .
Necesitamos determinar las constantes c0 , c1 y c2 de tal manera que L(ψp ) = x2 .

Tenemos
ψp0 (x) = c1 + 2c2 x, ψp00 (x) = 2c2 ,
y
L(ψp ) = (2c2 − 3c1 + 2c0 ) + (−6c2 + 2c1 )x + 2c2 x2 = x2 .
Ası́,
1 3
2c2 = 1, o c2 = y − 6c2 + 2c1 = 0, entoncesc1 =
2 2
7
2c2 − 3c1 + 2c0 = 0, entonces c0 =
4
1
por consiguiente ψp (x) = 4 (7 + 6x + 2x ) es una solución particular de L(y) = x2 .
2
Este método se llama de coeficientes indeterminados, ya que para resolver la

ecuación hay que determinar los coeficientes de la ecuación M , ası́ que el problema
es de naturaleza algebraı́ca y no es necesario efectuar integraciones.
(a) eax r−a

(b) xk eax (r − a)k+1
(c) sen(ax), cos(ax)a real r 2 + a2
(d) xk sen(ax), xk cos(ax)a real (r2 + a2 )k+1
Consideremos otro ejemplo para ilustrar el método. Sea la ecuación L(y) = Aeax ,
donde p es el polinomio caracterı́stico de L, y a, A son constantes. Supongamos que a
no es raı́z de p. El operador M dado por M (y) = y 0 −ay, cuyo polinomio caracterı́stico
es r−a, anula la función Aeax . El polinomio caracterı́stico de M L es (r−a)p(r) y a es
una raı́z simple de éste. Ası́ toda solución de la ecuación tiene la forma ψ = Beax + φ,
donde L(φ) = 0 y B es una constante. Sustituyendo ψ en la ecuación, obtenemos:
L(ψ) = BL(eax ) + L(φ) = Bp(a)eax = Aeax

A
como p(a) 6= 0, vemos que B = p(a) . Por consiguiente hemos demostrado que, si a no
es una raı́z del polinomio caracterı́stico de L, existe una solución ψ de la ecuación
de la forma
A ax
ψ(x) = e
p(a)
Capı́tulo 3
Ecuaciones Lineales con

Coeficientes Variables
3.1. Introducción
La ecuación diferencial lineal de orden n con coeficientes variables es una ecuación
de la forma
ao (x)y (n) + a1 (x)y (n−1) + . . . + an (x)y = b(x),
donde a0 6= 0, a1 , . . . , an y b son funciones de variable compleja en un intervalo I.
Supondremos que a( x) 6= 0 Al dividir la ecuación entre a0 obtenemos una ecuación
de la misma forma pero con coeficiente a0 igual a la constante 1. Ası́ que siempre
supondremos que a0 = 1, por lo que tenemos la ecuación
y (n) + a1 (x)y (n−1) + . . . + an (x)y = b(x). (3.1)

Representaremos el lado izquierdo de 3.1 por L(y). Tenemos la siguiente defini-
ción.
Definición 11 Consideremos la ecuación dada por 3.1 L(y) = b(x). Si b(x) = 0

∀x ∈ I la ecuación será L(y) = 0 se llama ecuación diferencial lineal ho-
mogénea con coeficientes variables. Si b(x) 6= 0, L(y) = b(x) se llama ecua-
ción diferencial lineal no homogénea con coeficientes variables.
Nota: L es un operador diferencial que opera sobre funciones que tienen n deri-
vadas en I, y transforma dichas funciones φ en funciones L(φ) cuyo valor en x está
dado por
59
60 3.2. Problemas con valores iniciales para ecuaciones homogeneas
L(φ)(x) = φ(n) (x) + a1 (x)φ(n−1) (x) + . . . + an (x)φ(x).

Por consiguiente una solución de L(y) = b(x) es una función φ que tiene n
derivadas en I, tales que L(φ) = b.
Si b es continua en I es posible obtener todas las soluciones de L(Y ) = b(x).
3.2. Problemas con valores iniciales para ecuacio-

nes homogeneas
Aun cuando no siempre es posible expresar una solución de 3.1 en términos de
funciones elementales, sı́ puede demostrarse que dichas soluciones existen. Admitire-
mos el siguiente teorema sin demostración.
Teorema 28 (Teorema de existencia.) Para cualquier x0 ∈ I y para las cons-

tantes α1 , · · · , αn existe una solución φ de la ecuación con condiciones iniciales
L(y) = 0, φ(x0 ) = α1 , · · · φ(n − 1)(x0 ) = αn .
∀x ∈ R.
Se obtienen dos cosas importantes del Teorema:
1. La solución existe en todo intervalo I donde a1 , . . . , an son continuas.
2. Cualquier problema con condiciones iniciales tiene una solución.
Lo anterior puede no ser cierto si se anula el coeficiente de y (n) para algún punto.
Consideremos la siguiente ecuación
xy 0 + y = 0,
cuyos coeficientes son continuos para toda x ∈ R. Esta ecuación junto con la condi-
ción inicial y(1) = 1, tiene l solución φ1 donde
1
φ1 (x) =
x
pero esta solución sólo existe para x > 0. También si φ es una solución cualquiera,
entonces
xφ(x) = c,
Capı́tulo 3. Ecuaciones Lineales con Coeficientes Variables 61
donde c es una cierta constante constante. Por consiguiente en el origen sólo existe
la solución trivial (c = 0), lo cual implica que el único problema con valor inicial
xy 0 + y = 0 y(0) = α1
que tiene solución, es aquel para el cual α1 = 0.
Definición 12 Dada una ecuación L(y) = 0 y una solución φ que tenga las deriva-
das φ0 , φ00 , · · · , φn−1 . Definimos ||φ(x)|| como
1/2
||φ(x)|| = |φ(x)|2 + · · · + |φ(n−1) (x)|2 .
L(y) = y (n) + a1 (x)y (n−1) + · · · + an (x)y = 0

en un intervalo I tal que x0 ∈ I y
|aj (x)| ≤ bj (j = 1, · · · , n),
para b1 , · · · , bn constantes no negativas.Entonces para toda x ∈ I:
||φ(x0 )||e−k|x−x0 | ≤ ||φ(x)|| ≤ ||φ(x0 )||ek|x−x0 | (3.2)

donde
k = 1 + b 1 + b2 + · · · + bn .
Demostración
Como L(φ) = 0 , tenemos:
φ(n) (x) = −a1 (x)φ(n−1) (x) − · · · − an (x)φ(x)
y, por consiguiente:
|φ(n) (x)| ≤ |a1 (x)||φ(n−1) (x)| + · · · + |an (x)||φ(x)|
≤ b1 |φ(n−1) (x)| + · · · + bn |φ(x)|

sustituyendo bj por |aj | se aplica la misma demostración.

Notemos que I debe de ser un intervalo cerrado y acotado, para poder garantizar
la existencia de las cotas.
62 3.3. Soluciones de Ecuaciones Homogéneas.
Teorema 30 Sean α1 , · · · , αn n constantes cualquiesra fijas, y sea x0 ∈ R. En cual-

quier intervalo I tal que xo ∈ I existe a lo más una solución φ de la ecuación
diferencial con condiciones iniciales
L(y) = 0, φ(x0 ) = α1 , · · · φ(n − 1)(x0 ) = αn .
Demostración
Sean φ, ψ dos soluciones de L(y) = 0 definidas en I que satisfacen las condiciones
del teorema, y consideremos χ = φ − ψ. Deseamos demostrar que χ(x) = 0 para
toda x en I. Aun cuando las funciones aj son continuas en I , no necesariamente
son acotadas en I. Por consiguiente no podemos emplear directamente el teorema
29. Sin embargo supongamos que x es cualquier punto en I, diferente de x0 , y sea J
cualquier intervalo cerrado y acotado que contenga a x, x0 y que este en I. En dicho
intervalo, las funciones aj son acotadas, esto es,
|aj (x)| ≤ bj (j = 1, 2, · · · , n),
en el intervalo J, para ciertas constantes bj , que pueden depender de J. Se cumple que

L(χ) = 0 en J, y que ||χ(x)|| = 0, y, por lo tanto φ(x) = ψ(x). Dada la arbitrariedad
de x entonces se cumple lo deseado.

3.3. Soluciones de Ecuaciones Homogéneas.

Al igual que en coeficientes constantes el espacio de soluciones cumple con las
condiciones bajo la tranformación lineal diferencial.
Teorema 31 Existe n soluciones de L(y) = 0, linealmente independientes en I.
Demostración
Sea x0 un punto en I. De acuerdo con el teorema 29 existe una solución φ1 de
L(y) = 0, y que ademas satisface las condiciones
(n−1)
φ1 (x0 ) = 1, φ01 (x0 ) = 0, · · · , φ1 (x0 ) = 0. (3.3)
En general, para cada i = 1, 2, · · · , n existe una solución φi que satisface las condi-
ciones
(i−1) (j−1)
φi (x0 ) = 1, φi (x0 ) = 0, j 6= i. (3.4)
Las soluciones φ1 , · · · , φn son linealmente independientes en I, para demostrarlo,

supongamos que existen ciertas constantes c1 , · · · cn tales que
c1 φ1 (x) + · · · + cn φn (x) = 0 (3.5)
para toda x ∈ I. Derivando, vemos que
c1 φ01 + · · · + cn φ0n = 0,
..
. (3.6)
(n−1)
c1 φ1 + · · · + cn φ(n−1)
n =0
para toda x ∈ I. Para x fija, las ecuaciones 3.4, 3.5, deben ser válidas en x0 .
Sustituyendo x = x0 en 3.4, vemos que de acuerdo con las condiciones iniciales,
c1 1 + 0 + · · · + 0 = 0 es decir c1 = 0. sustituyendo x = x0 en las ecuaciones 3.6
obtenemos c2 , c3 , · · · , cn = 0, por lo que concluimos que las soluciones φ1 , · · · , φn son
linealmente independientes.

Teorema 32 Sean φ1 , · · · , φn las n soluciones de L(y) = 0, definidas en I, que

además satisfacen las condiciones 3.4. Si φ es cualquier solución de L(y) = 0 en I,
entonces existen n constantes c1 , · · · , cn , tales que
φ = c1 φ1 + · · · + cn φn
Demostración
Sean
φ(x0 ) = α1 , φ0 (x0 ) = α2 , · · · , φ(n−1) (x0 ) = αn ,
y consideremos la función
ψ = α1 φ1 + α2 φ2 + · · · + αn φn
Que es una solución de L(y) = 0, y es claro que
ψ(x0 ) = α1 φ1 (x0 ) + α2 φ2 (x0 ) + · · · + αn φn (x0 ) = α1
ya que
φ1 (x0 ) = 1, φ2 (x0 ) = 0, · · · , φn (x0 ) = 0.
Empleando las relaciones de (3.4), vemos que
ψ(x0 ) = α1 , ψ 0 (x0 ) = α2 , · · · , ψ (n−1) (x0 ) = αn

64 3.4. Wronskiano e Independencia Lineal
Por lo tanto, ψ es una solución de L(y) = 0 que tiene las mismas condiciones iniciales
en x0 que φ. De acuerdo con el teorema 30, debemos tener φ = ψ, esto es:
φ = α1 φ1 + α2 φ2 + · · · + αn φn
Con esto queda demostrado el teorema para las constantes:
c1 = α 1 , c 2 = α 2 , · · · , c n = α n .

Si un conjunto de funciones tiene la propiedad de que, si φ1 , φ2 pertenecen al
conjunto, y c1 , c2 son dos constantes cualesquiera, entonces c1 φ1 + c2 φ2 pertenece
tambien al conjunto, por lo que tenemos un espacio lineal de funciones. Acabamos
de ver que el conjunto de soluciones es un espacio lineal de funciones. Ver base,
dimensión.
3.4. Wronskiano e Independencia Lineal

Tenemos la siguiente definición.
Definición 13 Sean φ1 , · · · , φn funciones, soluciones de L(y) = 0. Se define el

Wronskiano de φ1 , · · · , φn como

φ1 φ2
· · · φ n

0 0 0
φ φ
1 2 · · · φn
.. ..
.
n−1 n−1 . · · ·
n−1

φ1 φ2 · · · φn
Es una función, y su valor correspondiente a x, se representa por W (φ1 , · · · , φn )(x).
Teorema 33 Dadas φ1 , · · · , φn n soluciones de L(y) = 0 son l.i.en I, si y sólo si
W (φ1 , · · · , φn )(x) 6= 0
para toda x ∈ I.
Demostración
W (φ1 , · · · , φn ) 6= 0 ∀x ∈ I,
y consideremos constantes c1 , · · · , cn , tales que
c1 φ1 + · · · + cn φn = 0 (3.7)
para toda x ∈ I. Entonces tenemos
c1 φ01 + · · · + cn φ0n = 0,
..
. (3.8)
(n−1)
c1 φ1 + · · · + cn φ(n−1)
n =0
para toda x ∈ I. Para x fija, las ecuaciones 3.7, 3.8 son ecuaciones lineales homo-
geneas que se satisfacen para c1 , · · · , cn . El determinante del sistema es precisamente
W (φ1 , · · · , φn )(x), el cual no es nulo, por lo cual c1 = · · · = cn = 0 es la única
solución del sistema dado por 3.7, 3.8 por lo que φ1 , · · · , φn son l.i en I.
Recı́procamente supongamos que φ1 , · · · φn son l.i en I. Supongamos que existe un
x0 ∈ I tal que W (φ1 , · · · , φn )(x0 ) = 0. Esto implica que el sistema de n ecuaciones:
c1 φ1 + · · · + cn φn = 0
c1 φ01 + · · · + cn φ0n = 0,
.. (3.9)
.
(n−1)
c1 φ1 + · · · + cn φ(n−1)
n =0
tiene una solución c1 , · · · , cn , con al menos uno distinto de cero. Sea c1 , · · · , cn
una solución de este tipo, y consideremos la función ψ = c1 φ1 + · · · + cn φn . Ahora
L(ψ) = 0 por 3.9, tenemos que ψ(x0 ) = 0, · · · , ψ (n−1) (x0 ) = 0.
Luego por el Teorema 30 tenemos que ψ(x) = 0 para toda x ∈ I, luego
c1 φ1 + · · · + cn φn = 0
para toda x ∈ I, lo que contradice el hecho de que φ1 , · · · , φn son l.i. Por lo que
la hipótesis de que exista x0 ∈ I tal que W (φ1 , · · · , φn )(x0 ) = 0 es falsa, con lo que
se ha demostrado que
W (φ1 , · · · , φn )(x0 ) 6= 0

Teorema 34 Sean φ1 , · · · , φn las n soluciones de L(y) = 0, definidas en I, l.i. Si

φ es cualquier solución de L(y) = 0 en I, entonces existen n constantes c1 , · · · , cn ,
tales que
φ = c1 φ1 + · · · + cn φn .
Ası́, cualquier conjunto de n soluciones de L(y) = 0, l.i. en I, es una base para el
espacio de soluciones de L(y) = 0 en I
Demostración
Sean x0 un punto en I, y supongamos que
φ(x0 ) = α1 , φ0 (x0 ) = α2 , · · · , φ(n−1) (x0 ) = αn ,
ya demostramos que existen únicas constantes c1 , c2 , · · · , cn tales que
ψ = c1 φ1 + c2 φ2 + · · · + cn φn
es una solución de L(y) = 0, y es claro que
ψ(x0 ) = α1 , ψ 0 (x0 ) = α2 , · · · , ψn(n−1) (x0 ) = αn
Por el Teorema 30 se tiene entonces lo siguiente: φ = ψ, o bien
φ = c1 φ1 + · · · + cn φn
las condiciones iniciales para ψ son equivalentes a las siguientes ecuaciones para
c1 , · · · , cn :
c1 φ1 + · · · + cn φn = α1
c1 φ01 + · · · + cn φ0n = α2 ,
.. (3.10)
.
(n−1)
c1 φ1 + · · · + cn φ(n−1)
n = αn
Este es un conjunto de n ecuaciones lineales para c1 , · · · , cn . El determinante de

los coeficientes es W (φ1 , · · · , φn )(x0 ) el cual no es nulo ya que φ1 , · · · , φn son l.i.. Por
consiguiente, existe una única solución c1 , · · · , cn de las ecuaciones (3.10).

Veremos una formula para calcular el Wronskiano de n soluciones de L(y) = 0
en cualquier punto, a partir de un punto x0 ∈ I dado.
Teorema 35 Si φ1 , · · · , φn son n soluciones de L(Y ) = 0 en un intervalo I que

contenga x0 , entonces
Rx
− a1 (t)dt
W (φ1 , · · · , φn )(x) = e x0
W (φ1 , · · · , φn )(x0 ) (3.11)
Demostración
Se demostrará primero el caso para n = 2, y luego para el caso general.
Demostración para n = 2. En este caso
W (φ1 , φ2 ) = φ1 φ02 − φ01 φ2 ,
y, por consiguiente,
W 0 (φ1 , φ2 ) = φ1 φ002 − φ001 φ2 .
Como φ1 , φ2 satisface la ecuación y 00 + a1 (x)y 0 + a2 (x)y = 0, obtenemos:
φ001 = −a1 φ01 − a2 φ1
φ002 = −a1 φ02 − a2 φ2

Ası́:
W 0 (φ1 , φ2 )=φ1 (−a1 φ02 − a2 φ2 ) − (−a1 φ01 − a2 φ1 )φ2
=−a1 (φ1 φ02 − φ01 φ2 ) = −a1 W (φ1 , φ2 ).
Vemos que W (φ1 , φ2 ) satisface la ecuación lineal de primer orden
y 0 + a1 (x)y = 0,
y, por consiguiente, Rx
− a1 (t)dt
W (φ1 , φ2 )(x) = ce x0
donde c es una constante. Sustituyendo x = x0 , obtenemos
c = W (φ1 , φ2 )(x0 )
demostrando ası́ lo pedido para n = 2.

Demostración para el caso general n.
Para abreviar, escribiremos W = W (φ1 , . . . , φn ). A partir de la definición de W
como un determinante se sigue que su derivada W 0 es una suma de n determinantes
W 0 = V1 + · · · + Vn
Donde Vk difiere de W solamente en su k-ésimo renglón y donde el k-ésimo renglón

de Vk se obtiene derivando el k-ésimo renglón de W . Ası́:

φ01 · · · φ 0
n
φ 1 · · · φ n
φ1 ··· φn
0 0
00 00

φ1
· · · φ n
φ1
· · · φ n

φ01 ··· φ0n

00 00 00 00
W 0 = φ1
··· φn + φ1
··· φn + · · · +
φ001 ··· φ00n
.. .. .. .. .. ..

. . . . . .
(n−1) (n−1)
(n−1) (n−1)
(n) (n)

φ1 · · · φn φ1 · · · φn φ1 ··· φn
Los primeros n − 1 determinantes V1 , . . . , Vn−1 son iguales a cero por tener dos
renglones idénticos. como φ1 , . . . , φn son soluciones de L(y) = 0 se tiene lo siguiente:
(n) (n−1)
φi = −a1 φi − · · · − an φi , (i = 1, . . . , n),
y, por tanto,

φ1 ··· φn
φ01 φ0n

···
.. ..
W0 = . .

(n−1) (n−1)
φ1 ··· φn
(j) (j)
Pn−1 Pn−1
− · · · − j=0 an−j φn
j=0 an−j φ1

De donde obtenemos lo siguiente.

φ1 ··· φn
φ01 φ0n

···
.. ..

W0 = . .

φ(n−1) (n−1)

1 ··· φn
−a1 φ(n−1) (n−1)

1 · · · −a1 φn
por consiguiente, W satisface la ecuación lineal de primer orden y 0 + a1 (x)y = 0

y ası́
− xx a1 (t)dt
R
W (x) = e 0 W (x0 ).

Corolario 2 Si los coeficientes de L son constantes, entonces
W (φ1 , . . . , φn )(x) = e−a1 (x−x0 ) W (φ1 , . . . , φn )(x0 )

3.5. Reducción del orden de una ecuación homogénea

Supongamos que hemos hallado una solución φ1 de la ecuación:
L(y) = y (n) + a1 (x)y (n−1) + · · · + an (x)y = 0.
Encontraremos soluciones φ de la ecuación L(y) = 0, que tengan la forma φ = uφ1 .

Teorema 36 Sea φ una solución de la ecuación L(y) = 0 en un intervalo I, y

supongamos que φ1 (x) 6= 0 en I. Si v2 , · · · , vn es una base cualquiera en I, para el
espacio de soluciones de la ecuación lineal de orden n − 1 y si
vk = u0k , (k = 2, · · · , n)
entonces φ1 , u2 φ1 , · · · , un φ1 forman una base para el espacio de soluciones de L(y) =

0 en I.
Demostración
Supongamos que φ = uφ1 , es una solución, de la ecuación, entonces debe de
cumplirse lo siguiente:
0 = (uφ1 )(n) + a1 (uφ1 )(n−1) +· · · + an−1 (uφ1 )0 + an (uφ1 )
(n) (n−1)
= u(n) φ1 + · · · + uφ1 +a1 u(n−1) φ1 + · · · + a1 uφ1
+· · ·
+an−1 u0 φ1 + an−1 uφ01
+an uφ1
El coeficiente de u en esta ecuación es precisamente L(φ1 ) = 0. Por consiguiente, si

v = u0 , ésta es una ecuación lineal de orden n − 1 en v,
(n−1) (n−2)
φ1 v (n−1) + · · · + (nφ1 + (n − 1)φ1 + · · · + an−1 φ1 )v = 0 (3.12)
El coeficiente de v (n−1) es φ1 , y, por consiguiente, φ1 (x) 6= 0, en un intervalo I, esta

ecuación tiene n − 1 soluciones v2 , · · · , vn , linealmente independientes en I. Si x0 es
algún punto en I, y si
Z x
uk (x) = vk (t)dt, (k = 2, · · · , n),
x0
entonces se cumple u0k = vk , y las funciones
φ1 , u2 φ1 , · · · , un φ1 (3.13)
70 3.5. Reducción del orden de una ecuación homogénea

Para demostrar que forman una base del espacio de soluciones, es suficiente de-
mostrar que son l.i.. Para ello supongamos que existen constantes c1 , · · · , cn tales
que
c1 φ1 + c2 u2 φ1 + · · · + cn un φ1 = 0.
Como φ1 (x) 6= 0 en I, tenemos
c1 + c2 u2 + · · · + cn un = 0 (3.14)
y derivando obtenemos:
c2 u02 + · · · + cn u0n = 0
es decir
c2 v2 + · · · + cn vn = 0.
Como v1 , · · · , vn son l.i. entonces c2 , · · · , cn = 0 y usando 3.14, se tiene que c1 = 0.
Luego las funciones dadas por 3.13 son una base.

Veremos sólamente el caso para n = 2 del teorema, ya que en este caso la ecuación
para v es lineal de primer orden , y entonces la podemos resolver explı́citamente por
lo visto en el capı́tulo 1. Tenemos
L(y) = y 00 + a1 (x)y 0 + a2 (x)y = 0,
y si φ1 es una solución en I tenemos
L(uφ1 ) = (uφ1 )00 + a1 (uφ1 )0 + a2 (uφ1 )

= u00 φ1 + 2u0 φ01 + uφ001 + a1 u0 φ1 + a1 uφ01 + a2 uφ1
= u00 φ1 + u0 (2φ01 + a1 φ1 ).
Ası́, si v = u0 , y si u es tal que L(uφ1 ) = 0
φ1 v 0 + (2φ01 + a1 φ1 )v = 0 (3.15)
Notemos que (3.15) es una ecuación lineal de primer orden que se puede resolver
mediante un φ1 (x) 6= 0 en I. v satisface
φ21 v 0 + (2φφ01 + a1 φ21 )v = 0, (3.16)

El cual es justo (3.15) multiplicada por φ1 .Ası́
(φ21 v)0 + a1 (φ21 v) = 0,

lo cual implica que Rx

− a1 (t)dt
φ21 (x)v(x) = ce x0
,
donde x0 es un punto en I, y c es una constante, debido a que cualquier solución
multiplicada por una constante es también solución , tenemos que
1 − Rxx a1 (t)dt
v(x) = e 0 ,
φ21 (x)
es también solución de (3.16), y también de (3.15). Por lo que tenemos dos soluciones
independientes de
L(y) = y 00 + a1 y 0 + a2 y = 0 (3.17)
en I que son φ1 y φ2 , donde
Z x
1 − Rxs a1 (t)dt
φ2 (x) = φ1 (x) e 0 ds. (3.18)
x0 φ21 (s)
Teorema 37 Si φ1 es una solución de (3.17) en un intervalo I, y φ1 (x) 6= 0 en

I, una segunda solución φ2 de (3.17) está dado por (3.18). Las funciones φ1 y φ2
forman una base del espacio de solucionesde la ecuación diferencial en I.
Ejemplo.
3.6. Ecuación no Homogenea

Sean a1 , · · · , an , bfunciones continuas en un intervalo I y considere la ecuación
L(y) = y (n) + a1 (x)y (n−1) + · · · + an (x)y = b(x) (3.19)
Usaremos el método de variación de parámetros para resolver este tipo de ecuaciones,
ya que el método no depende de que los coeficientes sean o no constantes.
Teorema 38 Sea b una función continua en un intervalo I, y sean φ1 , · · · , φn una
base para el espacio de soluciones de L(y) = 0 en I. Cada solución ϕ de L(y) = b(x)
se puede escribir como
ϕ = ϕp + c1 φ1 + · · · + cn φn ,
donde ϕp es una solución particularde L(y) = b(x), y c1 , · · · , cn son constantes. Cada
ϕ es una solución de la ecuación no homogenea, Una solución particular está dada
por el método de variación de parámetros.
Ejemplo.
72 3.7. Solución con Series de Potencias
3.7. Solución con Series de Potencias

Ecuaciones

Cargado por

Copyright:

Formatos disponibles

Ecuaciones

Cargado por

Información del documento

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

Ecuaciones

Cargado por

Copyright:

Formatos disponibles

Ecuaciones Diferenciales Ordinarias.

Verónica Bolaños Granados

1. Ec Dif Ord de Primer Orden 5

2. Ec Dif Ord de Orden n 29

3. Ecuaciones Lineales con Coeficientes Variables 59

3.2. Problemas con valores iniciales para ecuaciones homogeneas . . . . . 60

Ecuación Diferencial Ordinaria de

Definición 1 (Ecuación Diferencial) Una ecuación que contiene las derivadas de

Como ejemplo tenemos la ecuaciones:

Las ecuaciones diferenciales se clasifican de acuerdo a su tipo , orden y linea-

segundo orden ↓ ↓ primer orden

y (n) = f (x, y, y 0 , . . . , y (n−1) )

an (x)y (n) + an−1 (x)y (n−1) + . . . + a1 (x)y 0 + a0 (x)y = g(x).

1.2. Ecuaciones Diferenciales Ordinarias Lineales

Consideremos una ecuacion diferencial lineal de primer orden.

y 0 + a(x)y = b(x) (1.3)

f (x, y) = −a(x)y + b(x). (1.4)

Teorema 1 Considere la ecuación diferencial ordinaria lineal homogénea

donde a es una constante compleja. Si c ∈ C cualquiera, la función definida por

φ0 (x) + aφ(x) = −ace−ax + ace−ax = 0.

eax (φ0 + aφ) = 0,

φ(x) = ce−ax . (1.6)

Teorema 2 Considere la ecuación diferencial ordinaria lineal no homogénea

Sustituyendo en 1.7 tenemos

eax φ(x) = B(x) + c

φ(x) = e−ax B(x) + ce−ax (1.9)

Ejemplos 1 Considere la ecuación

donde L, R, E, ω son constantes positivas.

(a) Calcule la solución φ que satisfaga la condición φ(0) = 0.

(b) Demuestre que la solución puede escribirse de la forma:

donde α es el ángulo que satisface las siguientes ecuaciones:

de donde tenemos lo siguiente.

por otra parte

Pasando la integral al primer miembro de la igualdad y simplificando

De donde obtenemos que

la C es la agrupación de la evaluación de xo con la c

con la condición φ(0) = 0 tenemos

Teorema 3 Supongamos que a y b son dos funciones continuas en un intervalo

donde xo ∈ I es una solución de la ecuación

y 0 + a(x)y = b(x) (1.10)

en I. La función φ1 dada por

Si c es una constante cualquiera, entonces φ = ψ + cφ1 es una solución de 1.10,

donde x0 ∈ I es un punto fijo,dado que

(eA φ)0 = eA φ0 + aeA φ = eA (φ0 + aφ).

Observemos que se tendrá φ0 + aφ = b si y sólo si

donde c es una constante, y B es una función cuya derivada es eA b. Podemos

Para demostrar la necesidad y suficiencia tenemos

(eA φ)0 = eA (φ0 + aφ) = eA b,

la última igualdad es por la hipótesis.

(eA (φ0 + aφ))(x) = (eA φ)0 (x) = (eA b)(x)

la última igualdad es por la hipótesis. entonces

eA(x) (φ0 + aφ)(x) = eA(x) b(x)

φ(x) = eA(x) B(x) + ce−A(x) , (1.13)

la segunda igualdad es por la hipótesis, la primera y la penúltima por el TFC.

1.3. Ecuaciones diferenciales lineales de primer or-

Teorema 4 Sean g y h, dos funciones continuas de variable real definidas en los

h(y)y 0 = g(x) (1.15)

Si H y G son dos funciones cualesquiera, tales que G0 = g, H 0 = h, y c es

define para x, en cierto intervalo I, contenido en a ≤ x ≤ b, a una función derivable

en el intervalo I, para cierta constante c.