Índice

1. Ecuaciones diferenciales lineales de orden superior 2

2. Ecuación lineal homogénea con coeficientes constantes 5

3. Método de Coeficientes Indeterminados 9

4. Método de variación de los parámetros 10

1
1. Ecuaciones diferenciales lineales de orden superior
Recordemos que una ecuación diferencial lineal de orden n es una expresión de la forma

an (x)y (n) + an−1 (x)y (n−1) + . . . + a1 (x)y ′ + a0 (x)y = b(x), (1)

donde an (x), . . . , a0 (x) y b(x) son funciones reales de variable real definidas sobre un intervalo
I . En el caso de que n = 1 tenemos la ecuación lineal de orden uno. Por otro lado, siempre
que an (x) sea distinto de cero, la ecuación anterior suele escribirse de la forma

y (n) + pn−1 (x)y (n−1) + ... + p1 (x)y ′ + p0 (x)y = q(x), (2)

ai (x) b(x)
donde pi = para i = 0, 1, . . . , n − 1 y q(x) = . La ecuacion (2) diremos que es
an (x) an (x)
homogenea si q(x) = 0 para todo x ∈ I . En caso contrario diremos que es no homogenea. Un
ejemplo de ecuacion homogenea, de tercer grado es

y ′′′ + xy ′′ + x2 y = 0,
mientras que sería no homogénea la ecuación,

y ′′′ + xy ′′ + x2 y = log x,

Teorema 1
Sean y1 , y2 , . . . , yn soluciones de la ecuación (2), con q(x) = 0, donde x pertenece a un
intervalo I , entonces la combinación lineal

y = c1 y1 + c2 y2 + . . . + cn yn ,

en el que c1 , c2 , . . . , cn son números reales, también es una solución para todo x ∈ I .

Corolario 1
Un multiplo constante y = c1 y1 (x) de una solución y1 (x) de una ecuación diferencial lineal
homogénea también es una solución.

Corolario 2
Una ecuación diferencial lineal homogénea siempre tiene la solución trivial y = 0.

2
1.1 Independencia lineal
Dadas las funciones f1 (x), f2 (x), . . . , fn (x), si existen constantes c1 , c2 , . . . , cn , no todas
cero, tales que:

c1 f1 (x) + c2 f2 (x) + . . . + cn fn (x) = 0 (3)

donde x ∈ I , entonces, diremos que las funciones f1 (x), f2 (x), . . . , fn (x) son linealmente
dependientes (LD), si tal relación no existe, diremos que las funciones son linealmente inde-
pendientes (LI). Esto es, las funciones f1 (x), f2 (x), . . . , fn (x) son linealmente independientes
cuando la ecuación (3) implica que c1 = c2 = . . . = cn = 0.

Observación: Si las funciones de un conjunto son linealmente dependientes, entonces al

menos una de ellas es una combinación lineal de las otras; si ellas son linealmente indepen-
dientes, ninguna es una combinación lineal de las otras.

1.2 El wronskiano
A continuación obtendremos una condición suficiente para que n funciones sean lineal-
mente independientes donde x ∈ I .
Supongamos que cada una de las funciones f1 (x), f2 (x), . . . , fn (x) es derivable al menos
(n − 1) veces, para toda x ∈ I , entonces a partir de la ecuación (3) por derivación sucesiva
obtenemos que:

c1 f1′ (x) + c2 f2′ (x) + ... + cn fn′ (x) = 0,

c1 f1′′ (x) + c2 f2′′ (x) + ... + cn fn′′ (x) = 0,

.. .. ..
. . .

(n−1) (n−1) ..
c1 f1 (x) + c2 f2 (x) + . + cn fn(n−1) (x) = 0.

Considerando como un sistema de ecuaciones en c1 , c2 , . . . , cn , las n ecuaciones lineales

anteriores no tendrán solución, excepto aquella que tenga todas las ci iguales a cero, si el
determinante del sistema no es nulo. Esto es, si

 
 f1′ (x) f2′ (x) ... fn′ (x) 
 
 
 
 f1′′ (x) f2′′ (x) ... fn (x)  ̸= 0
′′
 (4)
 .. .. .. .. 
 . . . . 
 (n−1) 
 f1 (x)
(n−1)
f2 (x) ...
(n−1)
fn (x) 

3
entonces las funciones f1 (x), f2 (x), . . . , fn (x) son linealmente independientes. El determi-
nante en (4) se llama el wronskiano de las n funciones involucradas.
La anulación del wronskiano es un intervalo, no es una condición necesaria para la inde-
pendencia lineal. El wronskiano puede anularse también cuando las funciones son linealmente
independientes
Si las n funciones involucradas son soluciones de una ecuación diferenciallineal homogénea,
la situación se simplifica como se muestra en el teorema siguiente.

Teorema 2
Si, sobre un intervalo I al que pertenece x, an (x) ̸= 0, an (x), an−1 (x) . . . , a0 (x), son constantes
e y1 , y2 , . . . , yn son soluciones de la ecuación

an y (n) + an−1 y (n−1) + . . . + a1 y ′ + a0 y = 0, (5)

entonces una condición necesaria y suficiente para que las y1 , y2 , . . . , yn sean linealmente
independientes es la no anulación del wronskiano de y1 , y2 , . . . , yn en el intervalo I .

1.3 Solución general de una ecuación diferncial homogénea

(EDH)
Sean y1 , y2 , . . . , yn soluciones linealmente independientes de la ecuación diferencial ho-
mogénea (5), entonces la solución general de la ecuación (5) es

y = c1 y1 + c2 y2 + . . . + cn yn , (6)
en el que c1 , c2 , . . . , cn son números reales arbitrarios.

1.3 Solución general de una ecuación diferncial no homogénea

(EDnoH)
En este apartado vamos a explicar cómo determinar la solución general de la ecuación
diferencial lineal de orden n no homogénea

Teorema 3
Dada una solución yp (x) de la ecuación diferencial (1) en un intervalo I , entonces para cualquier
solución y = ψ(x) de esta ecuación, existen constantes c1 , c2 , . . . , cn tales que

ψ(x) = c1 y1 + c2 y2 + . . . + cn yn + yp (x), x ∈ I, (7)

donde y1 , y2 , . . . , yn son soluciones LI en I de la ecuación (5) correspondiente.

4
Teniendo en cuenta (7), la solución general de la ecuación diferencial lineal no homogénea
(1) se define como:

y(x) = c1 y1 + c2 y2 + . . . + cn yn + yp (x), x ∈ I, (8)

Llamaremos a yp (x) una solución particular de (1). A la solución general de la ecuación ho-
mogénea (5) se le denomina solución complementaria y la denotaremos por yc (x). Es decir,

yc (x) = c1 y1 + c2 y2 + . . . + cn yn , (9)
En consecuencia, podemos escribir la solución general (8) de la ecuación no homogénea (1)
en la forma,

y(x) = yc (x) + yp (x), (10)

1.4 Método de Reducción de Orden

Dada una solución y1 (x) de la ecuación diferencial de segundo orden:

y ′′ + p(x)y ′ + q(x)y = 0, (11)

puede determinarse una segunda solución y2 (x) que sea linealmente independiente con y1 (x).

Teorema 4
Si y1 (x) es solución de la ecuación diferencial 11 entonces una segunda solución y2 (x) de
(11) LI con y1 (x) es
∫ ∫
e− p(x) dx
y2 (x) = y1 (x) dx. (12)
[y1 (x)]2

2. Ecuación lineal homogénea con coeficientes constantes

superior
Consideremos ecuaciones de la forma (5) donde an (x), an−1 (x) . . . , a0 (x) son constantes
reales. Un ejemplo de este tipo de ecuaciones puede ser,

y iv − y ′′′ + y ′′ + y = 0.
Vamos a ver cómo, de una manera bastante sencilla, es posible obtener todas las solu-
ciones de esta ecuación diferencial. Para ello empezaremos estudiando la ecuación diferencial
Homogénea de orden dos.

5
2.1 Ecuación de orden dos
Consideramos una ecuación de orden dos de la forma

a1 y ′′ + a2 y ′ + a3 y = 0,

dado que a1 ̸= 0 puede escribirse como

y ′′ + ay ′ + by = 0, (13)
a2 a3
donde a = yb= .
a1 a1
Para resolver esta ecuacion proponemos una solucion de la forma y(x) = erx , r ∈ R. Una
justificación para proponer esta solución viene de cómo son las soluciones de las ecuaciones
lineales homogéneas de orden uno con coeficientes constantes. Si

y ′ + ay = 0,
su solución es de la forma
y(x) = ce−ax, c ∈ R.

El buscar soluciones de la forma y(x) = erx , con el aumento del orden de la ecuación, no
debe de cambiar sustancialmente la forma de la solución.

Para que y(x) = erx sea solución de la ecuación (13) derivamos dos veces y(x) y sustitu-
imos en la ecuación, es decir, y ′ (x) = rerx y y(x) = r2 erx y la sustituimos en (13), teniendo

r2 erx + arerx + berx = 0 ⇒ erx (r2 + ar + b) = 0.

Como erx ̸= 0 para todo x ∈ R, necesariamente r2 + ar + b = 0, es decir, r debe ser una
raíz del polinomio p(x) = x2 + ax + b, que llamaremos polinomio característico de la ecuación
(13). Entonces, para obtener las soluciones de la ecuación diferencial debemos encontrar las
raíces de p(x). Analizando las soluciones de p(x) distinguimos los siguientes casos.

Caso 1: Raíces reales distintas

Consideremos la ecuación y ′′ + ay ′ + by = 0 y supongamos que las raíces del polinomio
x2 + ax + b son dos números reales r1 , r2 tal que r1 ̸= r2 , entonces las funciones y1 = er1 x
e y2 = er2 x son soluciones de la ecuación diferencial (13). Estas funciones son linealmente
independientes en (∞, ∞) y en consecuencia forman un conjunto fundamental, entonces, la
solución general de la ecuación (13) en ese intervalo es

y = C1 er1 x + C2 er2 x ,

donde C1 y C2 son números reales.

6
Caso 2: Raíces reales iguales
Consideremos la ecuación y ′′ + ay ′ + by = 0 y supongamos que las raíces del polinomio
2
x + ax + b son dos números reales r1 , r2 tal que r1 = r2 , entonces llegamos a una solución
b
y1 = er1 x , teniendo en cuenta la fórmula cuadrática r1 = − y a partir de la reducción de
2a
una ecuación diferencial de segundo orden a una de primer orden, obtenemos la segunda
solución y2 = xer1 x , en consecuencia estas funciones son linealmente independientes en
(∞, −∞) y en consecuencia forman un conjunto fundamental, entonces, la solución general
de la ecuación (13) en ese intervalo es

y = C1 er1 x + C2 xer1 x

Caso 3: Raíces complejos conjugados

Consideremos la ecuación y ′′ + ay ′ + by = 0 y supongamos que las raíces del polinomio
2
x +ax+b, si r1 y r2 son raíces complejas conjugadas podemos escribir r1 = a+ib y r2 = a−ib,
donde a y b > son números reales e i2 = −1, como son dos raíces y distintas, entonces la
solución la podemos escribir como,

y = C1 e(a+ib)x + C2 e(a−ib)x
A continuación obtendremos la forma de la solución general teniendo en cuenta las raíces
complejas conjugadas para ello utilizaremos la fórmula de Euler

eiy = cos(y) + i sen(y),

donde y es un número real. Teniendo en cuenta la fórmula de Euler tenemos

eiby = cos(by) + i sen(by) (14)

y
e−iby = cos(by) − i sen(by) (15)
A continuación si primero sumamos (14) + (15) y luego restamos (14) - (15), obtenemos

eiby + e−iby = 2 cos(bx) y eiby − e−iby = 2i sen(bx).

Como y = C1 e(a+ib)x + C2 e(a−ib)x es una solución de la ecuación (13) para cualquier elec-
ción de las constantes C1 y C2 , en particular si C1 = C2 = 1 y C1 = 1, C2 = −1, obtenemos
las soluciones:

y1 = e(a+ib)x + e(a−ib)x = eax (eiby + e−iby ) = 2eax cos(bx) (16)

y
y1 = e(a+ib)x + e(a−ib)x = eax (eiby + e−iby ) = 2ieax sen(bx). (17)

7
 En consecuencia teniendo en cuenta que “un multiplo constante de una solución de una
ecuación diferencial lineal homogénea también es una solución” los resultados (16) (17) de-
muestran que las funciones reales eax cos(bx) y eax sen(bx) son soluciones de la ecuación (13),
además forman un conjunto fundamental en (∞, −∞), por lo tanto la solución general está
dada por,
y = eax [C1 cos(bx) + C2 sen(bx)]

2.2 Ecuación de orden superior

En general para resolver una ecuación diferencial de orden n como:

an y (n) + an−1 y (n−1) + . . . + a1 y ′ + a0 y = 0, (18)

donde ai , con i = 0, 1 . . . , n, son constantes reales, debemos resolver una ecuación polinomial
de grado n:

an rn + an−1 rn−1 + . . . + a2 r2 + a1 r + a0 = 0. (19)

Si todas las raíces de la ecuación (19) son reales y distintas, la solución general de la
ecuación (18) es
y = C1 er1 x + C2 er2 x + . . . + Cn ern x .
Por otro lado si todas las raíces de la ecuación (19) son reales e iguales, es decir, de
multiplicidad n la solución general de la ecuación (18) es

y = C1 erx + C2 xerx + . . . + Cn xn−1 erx .

Cuando r1 = a+ib es una raíz compleja de multiplicidad k de (19), su conjugado r1 = a−ib

es también raízmultiplicidad k . En este caso, la solución general de la ecuación diferencial (18)
debe contener una combinación lineal de las siguientes 2k soluciones LI

eax cos(bx), xeax cos(bx), x2 eax cos(bx), . . . xk−1 eax cos(bx), k ∈ Z+

eax sen(bx), xeax sen(bx), x2 eax sen(bx), . . . xk−1 eax sen(bx), k ∈ Z+ .

Por tanto la solución general está dada por:

y = eax [C1 cos(bx) + C2 sen(bx)] + xeax [C3 cos(bx) + C4 sen(bx)]+

+x2 eax [C5 cos(bx)+C6 sen(bx)]+. . .+xk−1 eax [Cq cos(bx)+Cq+1 sen(bx)], q = 1, 3, 5, . . . 2k −1

8
3. Método de Coeficientes Indeterminados
Este método nos permite encontrar una solución particular yp (x) para ecuaciones diferen-
ciales lineales no homogéneos, mostraremos para el caso no homogéneo de segundo orden,
en el caso de orden n se sigue un procediminto análogo, sea

ay ′′ + by ′ + cy = g(x), (20)
donde a, b, c son constantes y


 bn xn + bn−1 xn−1 + . . . + b2 x2 + b1 x + b0 función polinomial

eax función exponenial
g(x) =

 sen(bx) función seno

cos(bx) función coseno

con el término no homogéneo g(x) correspondiente a una función polinómica, exponencial,

seno, coseno o suma finita de éstas, c.c. no es aplicable.

Este método consiste en hacer una suposición inicial para la solución particular con co-
eficientes sin especificar. La expresión supuesta se sustituye en la ecuación diferencial y se
intenta determinar los coeficientes. Si no es posible determinarlos significa que no es válida la
solución propuesta y debemos, por tanto, modificar la suposición inicial.

La limitación del método es que sólo suele funcionar bien con las ecuaciones con coefi-
cientes constantes y cuando el término no homogéneo tiene la forma indicada.

A continuación discutiremos algunos casos para hallar una solución particular de (20), de-
pendiendo de la forma de g(x).

Caso I: El término no homogéneo es un polinomio de grado n,

g(x) = an xn + an−1 xn−1 + . . . + a2 x2 + a1 x + a0 .

En este caso se propone como solución particular un polinomio del mismo grado al que le
multiplicamos el término xh :

yp (x) = (An xn + . . . + A1 + A0 )xh .

siendo An , An−1 , . . . , A1 , A0 coeficientes a determinar. El término xh se añade si 0 es una de

las soluciones de la ecuación característica, siendo h su orden de multiplicidad.

Caso II: El término no homogéneo es un producto de una exponencial y un polinomio de

grado n:
g(x) = eαx (an xn + an−1 xn−1 + . . . + a2 x2 + a1 x + a0 ).

9
En este caso se propone como solución particular el producto de un polinomio del mismo
grado multiplicado por la misma exponencial que aparece en g(x) y añadimos el término xh si
α es una de las raíces de la ecuación característica, siendo h su orden de multiplicidad:

yp (x) = eαx (An xn + . . . + A1 + A0 )xh .

con An , An−1 , . . . , A1 , A0 a determinar

Caso III: El término no homogéneo es un producto de una exponencial y una combinación

de coseno y seno del mismo ángulo multiplicados por polinomios:

g(x) = eαx [Pn (x)cos(bx) + Qm (x)sen(bx)].

En este caso la propuesta de solución particular es el producto de la función exponencial
eαx que aparece en g(x) multiplicado por, cos(bx) por un polinomio y sin(bx) por otro polinomio.
Estos polinomios tienen el mismo grado, el mayor grado de los que aparecen en g(x) y con
coeficientes diferentes a determinar. Añadimos el término xh si α + ib es una de las raíces de
la ecuación característica, siendo h su orden de multiplicidad.

g(x) = eαx [PN (x)cos(bx) + QN (x)sen(bx)]xh .

con PN (x) y QN (x) polinomios de grado N = máx{n, m} y con coeficientes a determinar para
cada uno de ellos.

Caso IV: Si el término no homogéneo es una suma de los casos anteriores, por el principio
de superposición, la propuesta de solución particular es una suma de las correspondientes
propuestas.

4. Método de variación de los parámetros

El método de variación de los parámetros es más general que el método de los coeficientes
indeterminados porque se puede aplicar también a ecuaciones lineales con coeficientes vari-
ables y para cualquiera que sea la forma del término no homogéneo g(x), mostraremos el
procediminto para el caso no homogéneo de segundo orden, el caso de orden n se sigue un
procediminto análogo.
Consideremos la ecuación lineal no homogénea de orden dos expresada en forma:

y ′′ + p(x)y ′ + q(x)y = g(x). (21)

El método consiste en buscar una solución de la forma

yp (x) = u1 (x)y1 (x) + u2 (x)y2 (x), (22)

donde y1 y y2 son dos soluciones LI de la ecuación diferencial homogénea asociada, y u1 , u2
son dos funciones a determinar de modo que (22) sea una solución de (21) y satisfagan una
condición arbitraria, pero seleccionada de tal forma que se simplifiquen los cálculos.

10
 Derivando (22) se tiene que

yp (x) = u′1 (x)y1 (x) + u1 (x)y1′ (x) + u′2 (x)y2 (x) + u2 (x)y2′ (x)

= [u1 (x)y1′ (x) + u2 (x)y2′ (x)] + [u′1 (x)y1 (x) + u′2 (x)y2 (x)]
Podemos simplificar esta expresión, imponiendo a u1 , u2 la condición de que

u′1 y1 + u′2 y2 = 0.

En este caso
yp′ = u1 (x)y1′ (x) + u2 (x)y2′ (x)
Por lo tanto

yp (x)′′ = u′1 (x)y1′ (x) + u1 (x)y1′′ (x) + u′2 (x)y2′ (x) + u2 (x)y2′′ (x)
Sustituyendo las expresiones de yp , yp (x)′ , yp (x)′′ en (21), y usando el hecho de que y1 y
y2 son soluciones de la ecuación diferencial homogénea, resulta

u′1 (x)y1′ (x) + u1 (x)y1′′ (x) + u′2 (x)y2′ (x) + u2 (x)y2′′ (x) + p(x)(u1 (x)y1′ (x) + u2 (x)y2′ (x)) +

+q(u1 (x)y1 (x) + u2 (x)y2 (x)) = g(x)

u′1 (x)y1′ (x) + u1 (x)y1′′ (x) + u′2 (x)y2′ (x) + u2 (x)y2′′ (x) + p(x)u1 (x)y1′ (x) + pu2 (x)y2′ (x) +

+q(x)u1 (x)y1 (x) + q(x)u2 (x)y2 (x) = g(x)

z }| { z }| {
u′1 (x)y1′ (x) + u1 (x)y1′′ (x) +u′2 (x)y2′ (x) ′′ ′ ′
+ u2 (x)y2 (x) + p(x)u1 (x)y1 (x) + p(x)u2 (x)y2 (x) +
| {z } | {z }
z }| {
+ q(x)u1 (x)y1 (x) + q(x)u2 (x)y2 (x) = g(x)
| {z }

u1 [y1′′ + p(x)y1′ + q(x)y1 ] + u2 [y2′′ + p(x)y2′ + q(x)y2 ] + u′1 (x)y1′ (x) + u′2 (x)y2′ (x) = g(x)
| {z } | {z }
cero cero

u′1 (x)y1′ (x) + u′2 (x)y2′ (x) = g(x)

11
 Así, buscamos una solución particular de la forma (22), con u1 y u2 funciones que satis-
facen las ecuaciones

u′1 y1 + u′2 y2 = 0,
(23)
u′1 (x)y1′ (x) + u′2 (x)y2′ (x) = g(x).

Es fácil resolver el sistema de ecuaciones (23) para las incógnitas u′1 y u′2 , empleando la
regla de Cramer. Obtenemos

y2 g(x) y1 g(x)
u′1 = − , u′2 = . (24)
W (x) W (x)

donde W (x) denota al wronskiano W (y1 , y2 )(x). Finalmente, integrando las expresiones
(24) resulta:

∫ ∫
y2 g(x) y1 g(x)
u1 = − dx, u2 = dx. (25)
W (x) W (x)

Sustituyendo (4.85) en (4.81) se obtiene la solución particular yp .

12