ECUACIONES DIFERENCIALES

1 Introducción a las ecuaciones diferenciales

1.1 Definiciones y clasificación
Definición. Una ecuación que relaciona una función desconocida con una o
más de sus derivadas es una ecuación diferencial.

Ejemplo 1. La ecuación diferencial

dx
= x2 + t2
dt
involucra tanto la función desconocida x(t) como su primera derivada x′ (t) =
dx
dt .

La ecuación diferencial

d2 y dy
+ 3 + 7y = 0
dx2 dx
incluye la función desconocida y de la variable independiente x y sus dos
primeras derivadas.
Clasificación. Las ecuaciones diferenciales se clasifican según el tipo, el
orden y la linealidad.

1
 Según el tipo las ecuaciones diferenciales pueden ser ordinarias o parciales.

Si una ecuación diferencial contiene sólo derivadas de una o más variables

dependientes respecto a una sola variable independiente se dice que es una
ecuación diferencial ordinaria (EDO). Por ejemplo,

dy d2 y dy dx dy
dx + 5y = ex , dx2 + 3 dx + 7y = 0 y dt + dt = 2x + y

son ecuaciones diferenciales ordinarias. Una ecuación que involucra derivadas

parciales de una o más variables dependientes respecto a dos o más variables
independientes se llama ecuación diferencial parcial (EDP). Por ejemplo,

∂2u ∂2u ∂2u ∂2u

∂x2 + ∂y 2 = 0, ∂x2 = ∂t2 − 2 ∂u
∂t y ∂u
∂y = ∂v
∂x

son ecuaciones diferenciales parciales.

El orden de una ecuación diferencial (ya sea EDO o EDP) es el orden de

la mayor derivada en la ecuación. Por ejemplo,

d2 y dy 3
+ 3( ) + 7y = 0
dx2 dx
es una ecuación diferencial ordinaria de segundo orden.

Simbólicamente podemos expresar una ecuación diferencial ordinaria de

2
n-ésimo orden con una variable dependiente por la forma general

F (x, y, y ′ , y ′′ , ..., y (n) ) = 0

, donde F es una función con valores reales de n+2 variables: x, y, y ′ , y ′′ , ..., y (n) .

La ecuación diferencial

dn y
n
= f (x, y, y ′ , y ′′ , ..., y (n−1) ),
dx
donde f es una función continua con valores reales, se conoce como la forma
normal de la ecuación F (x, y, y ′ , y ′′ , ..., y (n) ) = 0

Una ecuación diferencial de n-ésimo orden F (x, y, y ′ , y ′′ , ..., y (n) ) = 0 se

dice que es lineal si F es lineal en y, y ′ , y ′′ , ..., y (n) . Esto significa que una
EDO de n-ésimo orden es lineal cuando la ecuación F (x, y, y ′ , y ′′ , ..., y (n) ) = 0
se puede escribir de la forma

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

Según esta ecuación podemos ver que las dos propiedades características
de una ecuación diferencial lineal son las siguientes:

1. la variable dependiente y y todas sus derivadas y, y ′ , y ′′ , ..., y (n) son de

primer grado, es decir, la potencia de cada término que contiene a y es
1.
2. los coeficientes a0 , a1 , ..., an de y, y ′ , y ′′ , ..., y (n) dependen a lo más de la
variable independiente x.

Las ecuaciones

3
 d3 y
(y − x)dx + 4xdy = 0, y ′′ − 2y ′ + y = 0 y dx3
dy
+ x dx + 5y = ex

son ecuaciones diferenciales lineales. Si una ecuación diferencial no es lin-

eal se dice que es no lineal. Funciones no lineales de la variable dependiente
′
o de sus derivadas, tales como sin y o ey no pueden aparecer en una ecuación
lineal. Por tanto

d2 y d4 y
(1 − y)y ′ + 2y = ex , dx2 + sin y = 0 y dx4 + y2 = 0

son ejemplos de ecuaciones diferenciales ordinarias no lineales.

Definición. Cualquier función ϕ, definida en un intervalo I y que tiene

al menos n derivadas continuas en I, las cuales cuando se sustituyen en una
ecuación diferencial ordinaria de n-ésimo orden reducen la ecuación a una
identidad, se dice que es una solución de la ecuación en el intervalo.

En otras palabras, una solución de una ecuación diferencial ordinaria de

n-ésimo orden F (x, y, y ′ , y ′′ , ..., y (n) ) = 0 es una función ϕ que posee al menos
n derivadas para las que

F (x, ϕ(x), ϕ′ (x), ..., ϕ(n) (x)) = 0

para toda x en I.

Ejemplo 2. Verifique que la función indicada es una solución de la

ecuación diferencial dada en el intervalo (−∞, ∞).

4
 b) y ′′ − 2y ′ + y = 0; y = xex
dy 1 1 4
a) dx = xy 2 ; y = 16 x

En el ejemplo 2, cada ecuación diferencial tiene la solución constante

y = 0, (−∞, ∞). Una solución de una ecuación diferencial que es igual
a cero en un intervalo I se dice que es solución trivial.

Definición. Se dice que una relación G(x, y) = 0 es una solución im-

plícita de una ecuación diferencial ordinaria F (x, y, y ′ , y ′′ , ..., y (n) ) = 0 en un
intervalo I, si existe una función ϕ que satisface tanto la relación como la
ecuación diferencial en I.

Ejemplo 3. Verifique que la relación x2 + y 2 = 25 es una solución

implícita de la ecuación diferencial
dy x
=−
dx y

en el intervalo (−5, 5)

FAMILIA DE SOLUCIONES. Cuando resolvemos una ecuación difer-

encial de primer orden F (x, y, y ′ ) = 0, generalmente obtenemos una solución
que contiene una sola constante arbitraria o parámetro c. Una solución que
contiene una constante arbitraria representa un conjunto G(x, y, c) = 0 de
soluciones llamado familia de soluciones uniparamétrica. Cuando resolvemos
una ecuación diferencial de orden n F (x, y, y ′ , y ′′ , ..., y (n) ) = 0, buscamos una

5
familia de soluciones n-paramétrica G(x, y, c1 , c2 , ..., cn ) = 0. Esto significa
que una sola ecuación diferencial puede tener un número infinito de solu-
ciones correspondiendo a un número ilimitado de elecciones de los parámet-
ros. Una solución de una ecuación diferencial que esté libre de la elección
de parámetros se llama solución particular. Por ejemplo, la familia uni-
paramétrica y = cx − xcosx es una solución explícita de la ecuación lineal
de primer orden xy ′ − y = x2 sin x en el intervalo (−∞, ∞). La solución
y = −xcosx es una solución particular correspodiente a c = 0

Algunas veces una ecuación diferencial tiene una solución que no es miem-
bro de una familia de soluciones de la ecuación, es decir, una solución que
no se puede obtener usando un parámetro específico de la familia de solu-
1 4
ciones. Esta solución extra se llama solución singular. Por ejemplo, y = 16 x
1
y y = 0 son soluciones de la ecuación diferencial dx = xy 2 en (−∞, ∞).
dy

Esta ecuación tiene la familia de soluciones uniparamétrica y = ( 14 x2 + c)2 .

1 4
Cuando c = 0, se obtiene y = 16 x , pero y = 0 no se obtiene para ningún
valor de c, por tanto y = 0 es una solución singular.

SISTEMAS DE ECUACIONES DIFERENCIALES. Un sistema

de ecuaciones diferenciales ordinarias tiene dos o más ecuaciones que implican
derivadas de dos o más funciones incógnitas de una sola variable indepen-
diente. por ejemplo, si x y y son las variables dependientes y t denota la
variable independiente, entonces un sistema de dos ecuaciones diferenciales
de primer orden está dado por

dx
= f (t, x, y)
dt
dy
= g(t, x, y)
dt
Una solución de este sistema es un par de funciones derivables ϕ1 (t) y ϕ2 (t),

6
defindas en un intervalo común I, que satisface cada ecuación del sistema en
el intervalo.

1.2 Problemas con valores inciales

Con frecuencia nos interesan problemas en los que buscamos una solución
y(x) de una ecuación diferencial tal que y(x) satisface condiciones prescritas,
es decir, condiciones impuestas sobre una y(x) desconocida o sus derivadas.
En algún intervalo I que contiene a x0 , el problema
dn y
Resolver : n
= f (x, y, y ′ , y ′′ , ..., y (n−1) )
dx
Sujeto a y(x0 ) = y0 , y ′ (x0 ) = y1 , ..., y (n−1) (x0 ) = yn−1
donde y0 , y1 , ..., yn−1 son constantes reales arbitrarias dadas, se llama
problema con valores iniciales (PVI). Los valores de y(x) y de sus
primeras n−1 derivadas en un sólo punto x0 , y(x0 ) = y0 , y ′ (x0 ) = y1 , ..., y n−1 (x0 ) =
yn−1 , se llaman condiciones iniciales.

Un PVI de primer y segundo orden se enuncian respectivamente de la

forma siguiente:

dy
Resolver : = f (x, y)
dx
Sujeto a y(x0 ) = y0 (1)
y

d2 y
Resolver : 2
= f (x, y, y ′ )
dx
Sujeto a y(x0 ) = y0 , y ′ (x0 ) = y1 (2)

7
 Ejemplo 4. La función y = cex es una familia uniparamétrica de solu-
ciones de la ecuación de primer orden y ′ = y. Todas las soluciones en esta
familia están definida en el intervalo (−∞, ∞). Determine la solución par-
ticular que satisface la condición inicial y(0) = 3.

Ejemplo 5. c1 cos4t+c2 sin 4t es una familia de soluciones de dos parámet-

ros de x′′ + 16x = 0. Determine una solución del problema con valores
iniciales
π π
x′′ + 16x = 0, x( ) = −2, x′ ( ) = 1
2 2

Existencia y unicidad
Teorema. Sea R Una región rectangular en el plano xy definida por
a ≤ x ≤ b, c ≤ y ≤ d que contiene al punto (x0 , y0 ) en su interior. Si f (x, y)
∂y son continuas en R, entonces existe algún intervalo I0 : (x0 − h, x0 + h),
y ∂f
h > 0, contenido en [a, b], y una función única y(x), definida en I0 , que es
una solución del problema con valores iniciales (2).

Ejemplo. Consideremos el PVI

dy 1
= xy 2 , y(2) = 1
dx
Determine si tiene solución única.

8
1.3 Variables separables
Definición. Una ecuación diferencial de primer orden de la forma
dy
= g(x)h(y)
dx
se dice que separable o que tiene variables separables.

por ejemplo, las ecuaciones

dy dy
= y 2 xe3x+4y y = y + sin x
dx dx
son respectivamente, separable y no separable.

Ejemplo. Resolver:

1. (1 + x)dy − ydx = 0
2. dy
dx = − xy , y(4) = −3

3. dy
dx = y2 − 4

4. (e2y − y)cosx dx
dy
= ey sin 2x, y(0) = 0.

= e−x ,
dy 2
5. dx y(3) = 5.

1.4 Ecuaciones Lineales

Definición. Una ecuación diferencial de primer orden de la forma
dy
a1 (x) + a0 (x)y = g(x) (1)
dx

9
se dice que es una ecuación lineal en la variable dependiente y.

Forma estándar. Al dividir ambos lados de la ecuación (1) entre a1 (x)

se obtiene la forma estándar de la ecuación lineal:
dy
+ P (x)y = f (x) (2)
dx
Método para resolver ecuaciones lineales
1. Escriba la ecuación en la forma estándar.
dy
+ P (x)y = f (x)
dx
2. Calcula el factor integrante µ(x) mediante la fórmula
∫
P (x)dx
µ(x) = e

3. Multiplique la ecuación en la forma estándar por µ(x) para obtener

dy
µ(x) + µ(x)P (x)y = µ(x)f (x)
dx
d
(µ(x)y) = µ(x)f (x)
dx
4. Integre la última ecuación luego despeje y
Ejemplo. Resolver:

1. dy
dx − 3y = 6
dy
2. x dx − 4y = x6 ex
dy
3. dx + y = f (x), y(0) = 0 donde


 1 si 0 ≤ x ≤ 1

f (x) =

 0 si x > 1


10
1.5 Ecuaciones Exactas

Definición. El diferencial dz, también conocida como diferencial total, se

define como

∂z ∂z
dz = fx (x, y)dx + fy (x, y)dy = dx + dy
∂x ∂y

Definición. Una expresión diferencial M (x, y)dx + N (x, y)dy es una difer-
encial exacta en una región R del plano xy si ésta corresponde a la diferencial
de alguna función f (x, y) definida en R. Una ecuación diferencial de primer
orden de la forma
M (x, y)dx + N (x, y)dy = 0
se dice que es una ecuación exacta si la expresión del lado izquierdo es una
diferencial exacta.

Ejemplo. x2 y 3 dx + x3 y 2 dy = 0 es una ecuación exacta ya que el lado

izquierdo es una diferencial exacta:

1
d( x3 y 3 ) =
3
Teorema. Sean M (x, y) y N (x, y) continuas y que tienen primeras derivadas
parciales continuas en una región R definida por a ≤ x ≤ b, c ≤ x ≤ d. En-
tonces una condición necesaria y suficiente para que M (x, y)dx + N (x, y)dy
sea una diferencial exacta es
∂M ∂N
=
∂y ∂x
.

11
Ejemplo. Resuelve las siguientes ecuaciones diferenciales

(a) 2xydx + (x2 − 1)dy = 0

(b) (e2y − ycoxy)dx + (2xe2y − xcosxy + 2y)dy = 0

Factores integrantes. Si la ecuación diferencial M (x, y)dx+N (x, y)dy = 0

no es exacta, algunas veces es posible encontrar un factor integrante µ(x, y)
de manera que, después de multiplicar el lado izquierdo de la ecuación, la
expresión µM (x, y)dx + µN (x, y)dy es una diferencial exacta.

En el caso en que µ es una función de una variable, se podria obtener de la

siguiente forma:
My −Nx
(a) Si N depende sólo de x entonces
∫ My −Nx
dx
µ(x) = e N

Nx −My
(b) Si M depende sólo de y entonces
∫ Nx −My
dy
µ(y) = e M

Ejemplo. Resolver xydx + (2x2 + 3y 2 − 20)dy = 0

1.6 Soluciones por sustitución

Ecuaciones Homogéneas. Si una función f tiene la propiedad f (λx, λy) =

λn f (x, y) para algún número real n, entonces se dice que es una función ho-
mogénea de grado n. Por ejemplo f (x, y) = x3 + y 3 es homogénea de grado

12
3 mientras que f (x, y) = x3 + y 3 + 1 no es homogénea.

Una ecuación diferencial de la forma M (x, y)dx + N (x, y)dy = 0 es ho-

mogénea si ambas funciones coeficientes M y N son homogéneas de igual
grado. En tal caso cualquiera de las sustituciones y = vx o x = uy, donde u
y v son las nuevas variables dependientes, reducirán una ecuación homogénea
a una ecuación diferencial de primer orden separable.

Ejemplo. Resolver (x2 + y 2 )dx + (x2 − xy)dy = 0

Ecuación de Bernoulli. La Ecuación diferencial

dy
+ p(x)y = f (x)y n
dx
donde n es cualquier número real, se llama ecuación de Bernoulli. Para
n = 0 o n = 1 la ecuación es lineal. Para n ̸= 0 y n ̸= 1 la sustitución
u = y1−n la reduce a una ecuación lineal.

dy
Ejemplo. Resolver x dx + y = x2 y 2 .

Reducción a variables separables. Una ecuación diferencial de la forma

dy
= f (Ax + By + C)
dx
se puede reducir siempre a una ecuación con variables separables por medio
de la sustitución u = Ax + By + C, B ̸= 0.

13
Ejemplo. Resolver dy
dx = (−2x + y)2 − 7, y(0) = 0

1.7 Ecuaciones diferenciales como modelos matemáticos

La descripción matemática de un sistema de fenómenos se llama modelo

matemático. La formulación de un modelo matemático de un sistema se
inicia con
(a) identificación de las variables que ocasionan el cambio del sistema.
(b) Después se establece un conjunto de suposiciones razonables o hipóte-
sis, acerca del sistema que estamos tratando de describir. Esas
hipótesis también incluyen todas las leyes empíricas que se pueden
aplicar al sistema.
Con frecuencia, el modelo matemático de un sistema físico inducirá la vari-
able tiempo t. Una solución del modelo expresa el estado del sistema; en
otras palabras, los valores de la variable dependiente (o variables) para los
valores adecuados de t que describen el sistema en el pasado, presente y fu-
turo.

DINAMICA POBLACIONAL. Uno de los primeros intentos para mod-

elar el crecimiento de la población humana por medio de las matemáticas
fue realizado en 1798 por el economista inglés Thomás Malthus. La idea
detrás del modelo de Malthus es la suposición de que la razón de cambio de
la población de un país en un cierto tiempo es proporcional a la población
total del país en ese tiempo. En términos matemáticos, si P (t) denota la
población al tiempo t, entonces esta suposición se puede expresar como
dP
= kP
dt
14
donde k es una constante de proporcionalidad. este modelo falla si se con-
sideran muchos otros factores que pueden influir en el crecimiento o decrec-
imiento, como por ejemplo, la inmigración y la emigración.

Ejemplo 1. Inicialmente un cultivo tiene un número P0 de bacterias. En

t = 1 h se determina que el número de bacterias es 32 P0 . Si la razón de crec-
imiento es proporcional al número de bacterias P (t) presentes en el tiempo t,
determine el tiempo necesario para que se triplique el número de bacterias.

DECAIMIENTO RADIACTIVO. El núcleo de un átomo está formado

por combinaciones de protones y neutrones. Muchas de estas combinaciones
son inestables, esto es, los átomos se desintegran o se convierten en átomos de
otras sustancias. Se dice que estos átomos son radiactivos. Por ejemplo, con
el tiempo, el radio Ra − 226 , intensamente radiactivo, se transforma en el
radiactivo gas Radón, Rn − 222. Para modelar el fenómeno del decaimiento
radiactivo, se supone que la razón dAdt con la que los núcleos de una sustancia
se desintegran es proporcional al número de núcleos A(t) de la sustancia que
queda al tiempo t:
dA
= kA.
dt

En física la vida media es una medida de la estabilidad de una sustancia

radiactiva. La vida media es el tiempo que tarda en desintegrarse o trans-
mutarse en otro elemento la mitad de los átomos en una muestra inicial A0 .
Mientras mayor sea la vida media de una sustancia, más estable es la sustan-
cia. Por ejemplo, la vida media del radio altamente radiactivo Ra − 226 es
de aproximadamente 1700 años. El isótopo más común del uranio, U − 238,

15
tiene una vida media de 4 500 000 000 años. En aproximadamente 4.5 miles
de millones de años la mitad de una cantidad de U − 238 se transmuta en
plomo 206.

Ejemplo 2. Un reactor de cría convierte uranio 238 relativamente estable

en isótopo plutonio 239. Después de 15 años se ha determinado que 0.043%
de la cantidad inicial A0 de plutonio se ha desintegrado. Determine la vida
media de este isótopo, si la razón de desintegración es proporcional a la can-
tidad que queda.

DATADO CON CARBONO. Alrededor de 1950, el químico Willard

Libby inventó un método que utiliza el carbono radiactivo para determi-
nar las edades aproximadas de fósiles. La teoría del datado con carbono, se
basa en que el isótopo carbono 14 se produce en la atmósfera por acción de la
radiación cósmica sobre el nitrógeno. La razón de la cantidad de C − 14 con
el carbono ordinario en la atmósfera parece ser costante y, en consecuencia,
la cantidad proporcional del isótopo presente en todos los organismos vivos
es igual que la de la atmósfera. Cuando muere un organismo cesa la obsor-
ción del C − 14 ya sea por respiración o por alimentación. Así al comparar
la cantidad proporcional del C − 14 presente, por ejemplo en un fósil con la
razón constante que hay en la atmósfera, es posible obtener una estimación
razonable de la edad del fósil. Este método se basa en que se sabe que la
vida media del C − 14 es de aproximadamente 5600 años.

Ejemplo 3. Se encuentra que un hueso fosilizado contiene la milésima parte

de la cantidad de C − 14 encontrada en la materia viva. Determine la edad
del fósil.

16
LEY DE NEWTON DE CALENTAMIENTO/ENFRIAMIENTO.
De acuerdo con la ley empírica de calentamiento/enfriamiento, la rapidez con
la que cambia la temperatura de un cuerpo es proporcional a la diferencia en-
tre la temperatura del cuerpo y la del medio que le rodea. Si T (t) representa
la temperatura del cuerpo al instante t, Tm es la temperatura del ambiente
que le rodea y dT
dt es la rapidez con que cambia la temperatura del cuerpo,
entonces la expresión matemática de esta ley es
dT
= k(T − Tm )
dt
donde k es una constante de proporcionalidad.

Ejemplo 4. Al sacar un pastel del horno, su temperatura es 300◦ F . Tres

minutos después su temperatura es de 200◦ F . ¿Cuánto tiempo le tomará al
pastel enfriarse hasta una temperatura de 80◦ F si la temperatura del ambi-
ente es 70◦ F ?.

CIRCUITOS EN SERIE. Para un circuito en serie que sólo contiene un

resistor y un inductor la segunda ley de Kirchhoff establece que la suma de la
caida de voltaje a través del inductor (L dI dt ) más la caída de voltaje a través
del resistor (IR) es igual al voltaje aplicado(E(t)) al circuito. Por tanto obte-
nenos la ecuación diferencial para la corriente I(t),

dI
L + RI = E(t)
dt
, donde L y R son constantes conocidas como la inductancia y la resistencia,
respectivamente.

17
Ejemplo 5. Una batería de 12 voltios se conecta a un circuito en serie en
el que el inductor es de 12 henry y la resistencia es de 10 ohms. Determine la
corriente I si la corrien te inicial es cero.

2 Ecuaciones diferenciales de orden superior

2.1 Ecuaciones lineales

PROBLEMAS CON VALORES INICIALES. Para una ecuación difer-

encial lineal, un problema con valores iniciales de n-ésimo orden es
dn y dn−1 y dy
Resuelva an (x) n + an−1 (x) n−1 + ... + a1 (x) + a0 (x)y = g(x) (1)
x dx dx

Sujeto a y(x0 ) = y0 , y ′ (x0 ) = y1 , ..., y (n−1) (x0 ) = yn−1

.

Teorema 1. Existencia de una solución única.

Sean an (x), an−1 (x),..., a1 (x), a0 (x) y g(x) continuas en un intervalo I, y sea
an (x) ̸= 0 para toda x en ese intervalo. Si x = x0 es cualquier punto en ese
intervalo, entonces una solución y(x) del problema con valores inciales (1)
existe en el intervalo y es única.

18
Ejemplo 1. La función y = 3e2x + e−2x − 3x es una solución del problema
con valores iniciales
y ′′ − 4y = 12x, y(0) = 4, y ′ (0) = 1.
La ecuación diferencial es lineal; los coeficientes, así como g(x) = 12x son
continuos y a2 (x) ̸= 0 en cualquier intervalo I que contenga a x = 0. Con-
cluimos del teorema 1 que la función dada es la única solución en I.

2.1.1 Ecuaciones Homogéneas

Una ecuación diferencial lineal de n-ésimo orden de la forma

dn y dn−1 y dy
an (x) n + an−1 (x) n−1 + ... + a1 (x) + a0 (x)y = 0 (2)
dx dx dx
se dice que es homogénea, mientras que una ecuación

dn y dn−1 y dy
an (x) n + an−1 (x) n−1 + ... + a1 (x) + a0 (x)y = g(x, ) (3)
dx dx dx

con g(x) ̸= 0 se dice que es no homogénea. Por ejemplo, 2y ′′ + 3y ′ − 5y = 0

es una ecuación diferencial lineal homogénea de segundo orden, mientras que
x3 y ′′′ + 6y ′ + 10y = ex es una ecuación diferencial lineal de tercer orden no
homogénea.

OPERADORES DIFERENCIALES. En cálculo la derivación se denota

dy
con frecuencia con la letra D, es decir, dx = Dy. El símbolo D se llama

19
operador diferencial porque convierte una función derivable en otra función.

Las derivadas de orden superior se expresan en términos de D de manera

natural:
d dy d2 y 2 dn y n
(
dx dx ) = dx 2 = D(Dy) = D y y en general dxn = D y,
donde y representa una función suficientemente derivable. Las expresiones
polinomiales en las que interviene D, tales como D + 3 y D2 + 3D − 4
también son operadores diferenciales. Se define un operador diferencial de
n-ésimo orden u operador polinomial como
L = an (x)Dn + an−1 (x)Dn−1 + ... + a1 (x)D + a0 (0). (4)

El operador diferencial lineal L es lineal, es decir,

L{ αf (x) + βg(x)} = αL(f (x)) + βL(g(x)),
donde α y β son constantes.

ECUACIONES DIFERENCIALES. Cualquier ecuación diferencial lin-

eal puede expresarse en términos de la notación D.

Ejemplo 2. Exprese en términos de D la ecuación diferencial y ′′ +5y ′ +6y =

5x − 3.

Usando la ecuación (4), las ecuaciones diferenciales lineales de n-ésimo orden

(2) y (3) se puede expresar en forma compacta como
L(y) = 0 y L(y) = g(x)

20
Teorema 2. Principio de superposición; ecuaciones homogéneas
Sean y1 , y2 , ..., yk soluciones de la ecuación heomogénea de n-ésimo orden (2)
en un intervalo I. Entonces la combinación lineal

y = c1 y1 (x) + c2 y2 (x) + ... + ck yk (x),

donde las ci , i = 1, 2, ..., k son constantes arbitrarias, también es solución en

el intervalo.

Corolario 1. Del Teorema anterior se deduce que:

(a) Un múltiplo constante c1 y1 (x) de una solución y1 (x) de una ecuación

diferencial lineal homogénea es también una solución.
(b) Una ecuación diferencial lineal homogénea tiene siempre la solución
trivial y = 0
Ejemplo 2. Las funciones y1 = x2 y y2 = x2 ln x son soluciones de
la ecuación lineal homogénea x3 y ′′′ − 2xy ′ + 4y = 0 en el intervalo
(−∞, ∞). Por el principio de superposición la combinación lineal

c1 x2 + c2 x2 ln x

es también una solución de la ecuación en el intervalo.

Definición. Se dice que un conjunto de funciones f1 (x), f2 (x), ..., fn (x)

es linealmente dependiente en un intervalo I si existen constantes
c1 , c2 , ..., cc no todas cero, tales que

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

para toda x en I. Si el conjunto de funciones no es linealmente

dependiente en I, se dice que es linealmente independiente.

21
Ejemplo 3. El conjunto de funciones f1 (x) = cos2 x, f2 (x) =
sen2 x, f3 (x) = sec2 x, f4 (x) = tan2 x es linealmente dependiente
en (− π2 , π2 )

Definición. Wronskiano
Suponga que cada una de las funciones f1 (x), f2 (x), ..., fn (x) tienen
al menos n − 1 derivadas. El determinante
 
 f1 f ... f 
 2 n 
 f′ ′ ′ 
 1 f 2 ... f n 
W (f1 , f2 , ..., fn ) =  .. .. .. .. 
 . . . . 
 (n−1) (n−1) 
f1 f2
(n−1)
... fn 
donde las primas denotan derivadas, se llama Wronskiano de las
funciones.

Teorema 3. Criterio para soluciones linealmente independientes.

Sean y1 , y2 , ..., yn n soluciones de la ecuación diferencial lineal ho-
mogénea de n-ésimo orden (2) en el intervalo I. El conjunto de solu-
ciones es linealmente independiente en (I) si y sólo si W (y1 , y2 , ..., yn ) ̸=
0 para toda x en el intervalo.

Definición. Conjunto fundamental de soluciones.

Cualquier conjunto y1 , y2 , ..., yn de n soluciones linealmente inde-
pendientes de la ecuación diferencial lineal homogénea de n-ésimo
orden (2) en un intervalo I es un conjunto fundamental de solu-
ciones en el intervalo.

22
Teorema 4. Existencia de un conjunto fundamental .
Existe un conjunto fundamental de soluciones para la ecuación
diferencial lineal homogénea de n-ésimo orden (2) en un intervalo
I.

Teorema 5. Solución general; ecuaciones homogéneas.

Sean y1 , y2 , ..., yn un conjunto fundamental de soluciones para la
ecuación diferencial lineal homogénea de n-ésimo orden (2) en un
intervalo I. Entonces la solución general de la ecuación en el inter-
valo es

y = c1 y1 (x) + c2 y2 (x) + ... + cn yn (x),

donde las ci , i = 1, 2, ..., n son constantes arbitrarias.

Ejemplo 4. Las funciones y1 = e3x y y2 = e−3x son soluciones de

la ecuación lineal homogénea y ′′ − 9y = 0 en el intervalo (−∞, ∞).
Verifique que estas soluciones constituyen un conjunto fundamen-
tal de soluciones y determine la solución general de la ecuación.

Ejemplo . La función y = 4 sinh 3x − 5e3x es una solución de la

ecuación diferencial del ejemplo anterior. Aplicando el teorema 5
explique como obtener esta solución a partir de la solución general
y = c1 e3x + c2 e−3x .

23
2.2 Reducción de orden

Consideremos la ecuación diferencial lineal homogénea de segundo

orden
a2 (x)y ′′ + a1 (x)y ′ + a0 (x)y = 0 (1)
y suponga que y1 es una solución no trivial en un intervalo I. Se
busca una segunda solución y2 tal que y1 y y2 sean un conjunto
linealmente independiente en I. Esto significa que yy12 = u(x).
La función u(x) se determina al sustituir y2 (x) = u(x)y1 (x) en
la ecuación diferencial dada. Este método se llama reducción de
orden porque debemos resolver una ecuación diferencial lineal de
primer orden para encontrar a u.

llevando la ecuación (1) a la forma estándar

y ′′ + P (x)y ′ + Q(x)y = f (x)
podemos probar que una segunda solución y2 , a partir de una solu-
ción y1 , tal que y1 y y2 sean linealmente independiente, está dada
por
∫ ∫
− P (x)dx
e
y2 = y1 dx
y12
Ejemplo . La función y1 = x2 es solución de x2 y ′′ − 3xy ′ + 4y = 0.
Encuentre la solución general de la ecuación en el intervalo (0, ∞)

2.3 Ecuaciones lineales homogéneas con coeficientes con-

stantes

Consideremos la ecuación de segundo orden

ay ′′ + by ′ + cy = 0 (1)

24
donde a, b y c son constantes. Se intenta encontrar una solución
de la forma y = emx , entonces después de sustituir y, y ′ y y ′′ la
ecuación (1) se convierte en

emx (am2 + bm + c) = 0.
Debido a que emx ̸= 0, la única forma en que y = emx es solución
de la ecuación (1) es cuando se elige m como una solución de la
ecuación cuadrática
am2 + bm + c = 0.
Esta última ecuación se llama ecuación auxiliar de la ecuación difer-
encial (1).

Hay tres casos a considerar dependiendo de la naturaleza de las

raices m1 y m2 de la ecuación auxiliar. Si las raíces m1 y m2 son:
i. Reales y distintas, entonces y1 = em1 x y y2 = em2 x y como
estas funciones son linealmente independientes en (−∞, ∞)
entonces se deduce que la solución general de la ecuación (1)
en este intervalo es
y = c1 em1 x + c2 em2 x
.
ii. Reales e iguales, entonces y1 = em1 x y usando reducción de
orden podemos probar que una segunda solución linealmente
independiente en (−∞, ∞) es y2 = xem1 x y por tanto la solu-
ción general de la ecuación (1) en este intervalo es
y = c1 em1 x + c2 xem1 x
.

25
iii. Complejas conjugadas de la forma α±βi, donde α y β > 0 son
reales, entonces usando la fórmula de Euler podemos probar
que y1 = eαx cosβx y y2 = eαx senβx constituyen un conjunto
fundamental de soluciones de (1) en (−∞, ∞) y por tanto la
solución genral de la ecuación es

y = eαx (c1 cosβx + c2 senβx)

Ejemplo . Resolver

a) 2y ′′ − 5y ′ − 3y = 0 b) y ′′ − 10y ′ + 25y = 0

c) y ′′ + 4y ′ + 7y = 0

ECUACIONES DE ORDEN SUPERIOR. Para resolver una

ecuación diferencial lineal de n-ésimo orden

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

donde ai , i = 0, 1..., n, son constantes reales, se debe resolver una
ecuación polinomial de n-ésimo grado
an mn + an−1 mn−1 + ... + a1 m + a0 = 0

Si todas las raíces de esta ecuación son reales y distintas entonces

la solución general de (2) es
y = c1 em1 x + c2 em2 x + ... + cn emn x

Si m1 es una raíz de multiplicicada k de la ecuación auxiliar, en-

tonces la solución general de la ecuación (2) debe contener la com-
binación lineal

26
 c1 em1 x + c2 xem1 x + ... + ck xk−1 em1 x

Si m1 = α + βi, β > 0 es una raíz compleja de multiplicidad k

entonces m2 = α − βi también es una raíz de multiplicidad k.
Con ayuda de la fórmula de Euler podemos probar que la solución
general de la ecuacón (2) debe contener la combinación lineal de
las 2k soluciones linealmente independientes:

eαx cosβx, xeαx cosβx, ...., xk−1 eαx cosβx

eαx senβx, xeαx senβx, ...., xk−1 eαx senβx

Ejemplo . Resolver

a) y ′′′ + 3y ′′ − 4y = 0

b) y 4 + 2y ′′ + y = 0

2.4 Ecuaciones no homogéneas

Consideremos la ecuación diferencial lineal no homogénea de orden

n

dn y dn−1 y dy
an (x) + an−1 (x) + ... + a1 (x) + a0 (x)y = g(x, ) (1)
dxn dxn−1 dx

27
Cualquier función yp libre de parámetros arbitrarios, que satisface
(1) se dice que es una solución particular o integral particular de
la ecuación.

Teorema . Solución general: ecuaciones no homogéneas

Sea yp cualquier solución particular de la ecuación diferencial lineal

no homogénea de orden n (1) en un intervalo I, y sea y1 , y2 , ... yn
un conjunto fundamental de soluciones de la ecuación diferencial
lineal homogénea asociada en I. Entonces la solución general de la
ecuación en I es

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

donde las ci , i = 1, 2, ..., n son constantes arbitrarias.

La función yc = c1 y1 (x)+c2 y2 (x)+...+cn yn (x), que forma parte de

la solución general de la ecuación diferencial lineal no homogénea
y que es la solución general ecuación diferencial lineal homogénea
asociada, se llama función complentaria para la ecuación (1). Esto
significa que para resolver una ecuación diferencial lineal no ho-
mogénea primero se resuelve la ecuación homogénea asociada y
luego se encuetra una solución particular de la ecuación no ho-
mogénea. La solución general de la ecuación no homogénea es

y = yc + yp .

28
Método de coeficientes indeterminados

La primera de las dos formas que se consideran para obtener una

solución particular yp de una EDL no homogénea se llama método
de coeficientes indeterminados. La idea fundamental de este
método es una conjetura acerca de la forma de yp motivada por
las clases de funciones que forman la función de entrada g(x). El
método general se limita a EDL como (1) donde

• Los coeficientes ai , i = 0, 1, ..., n son constantes y

• g(x) es una constante k, una función polinomial, una función
exponencial eax , una función senβx o cosβx o sumas finitas y
productos de estas funciones.

Ejemplo . Resolver

a)y ′′ + 4y ′ − 2y = 2x2 − 3x + 6

b)y ′′ − y ′ + y = 2sen3x

c) y ′′ − 5y ′ + 4y = 8ex

Método de variación de parámetros

Consideremos la ecuación diferencial lineal no homogénea de se-

gundo orden en su forma estándar

29
 y ′′ + P (x)y ′ + Q(x)y = f (x).

Sea yc = c1 y1 + c2 y2 la función complementaria. El método de vari-

acíon de parámetros consiste en determinar la solución particular yp
a partir de la función complementaria sustituyendo los parámetros
c1 y c2 por variables, es decir,

yp = u1 y1 + u2 y2
donde u1 y u2 son funciones de x. Podemos probar que
∫ ∫
w1 w2
u1 = dx y u2 = dx
w w
donde w es el wronskiano de y1 y y2 , w1 = −y2 f (x) y w2 = y1 f (x).

Ejemplo . Resolver

a)y ′′ − 4y ′ + 4y = (x + 1)e2x b)4y ′′ + 36y = csc3x

ECUACIONES DE ORDEN SUPERIOR. Consideremos la

ecuación diferencial lineal no homogénea de orden n en su forma
estándar

y (n) + Pn−1 (x)y (n−1) + ... + P1 (x)y ′ + P0 (x)y = f (x) (2)

Si yc = c1 y1 (x)+c2 y2 (x)+...+cn yn (x) es la función complementaria

para (2) entonces una solución particular es

30
 yp = u1 y1 (x) + u2 y2 (x) + ... + un yn (x),

donde los u′k , k = 1, 2, ..., n se determinan por las n ecuaciones

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

y1′ u′1 + y2′ u′2 + ... + yn′ u′n = 0
.. ..
. .
(n−1) ′ (n−1) ′
y1 u1 + y2 u2 + ... + yn(n−1) u′n = f (x)

Usando la regla de Cramer se obtiene

wk
u′k = , k = 1, 2, ..., n
w
donde w es el Wronskiano de y1 , y2 , ..., yn y wk es el determinante
que se obtiene al reeplazar la k-ésima columna del Wronskiano por
la columna que consta de (0, 0, ..., f (x)).

2.5 Ecuación de Cauchy-Euler

Una ecuación diferencial de la forma

n n−1
nd y n−1 d y dy
an x + an−1 x + ... + a1 x + a0 y = g(x)
dxn dxn−1 dx

donde los coeficientes an , an−1 , ..., a0 son constantes, se conoce como

ecuación de Cauchy-Euler. La característica observable de este tipo
de ecuación es que el grado k = n, n − 1, ..., 1, 0 de los coeficientes

31
 dk y
monomiales xk coincide con el orden k de la derivación dxk
.

Consideremos la ecuación homogénea de segundo orden

2
2d y dy
ax + bx + cy = 0
x2 dx
La solución de ecuaciones de orden superior se deduce de manera
análoga. También, podemos resolver la ecuación no homogénea
ax2 y ′′ + bxy ′ + cy = g(x) por variación de parámetros, una vez que
se ha determinado la función complementaria yc .

NOTA El coeficiente ax2 de y ′′ es cero en x = 0. por lo que,

para garantizar que los resultados fundamentales del teorema de
existencia y unicidad que hemos enunciados antes sean aplicables
a la ecuación de Cauchy-Euler, centramos nuestra atención en en-
contrar soluciones generales definidas en el intervalo (0, ∞). Las
soluciones en el intervalo (−∞, 0) se obtienen al sustituir t = −x
en la ecuacion diferencial.

Método de solución. Se prueba con una solución de la forma

y = xm , donde m es un valor que se debe determinar.

Cuando sustituimos y = xm en la ecuación

2
2d y dy
ax + bx + cy = 0
x2 dx

32
se obtiene la ecuación auxiliar am2 + bm + c = 0. Semejante a lo
que sucede con la ecuaciones lineales con coeficientes constantes,
aquí hay tres casos distintos a considerar dependiendo de la natu-
raleza de las raíces de la ecuación auxiliar.

Si las raíces m1 y m2 son:

i. Reales y distintas, entonces y1 = xm1 y y2 = xm2 y como estas
funciones son linealmente independientes en (0, ∞) entonces se
deduce que la solución general de la ecuación en este intervalo
es
y = c1 xm1 + c2 xm2
.
ii. Reales e iguales, entonces y1 = xm1 y usando reducción de
orden podemos probar que una segunda solución linealmente
independiente en (0, ∞) es y2 = xm1 ln x y por tanto la solución
general de la ecuación (1) en este intervalo es
y = c1 xm1 + c2 xm1 ln x
.
iii. Complejas conjugadas de la forma α ± βi, donde α y β >
0 son reales, entonces usando la fórmula de Euler podemos
probar que y1 = xα cosβ ln x y y2 = xα senβ ln x constituyen
un conjunto fundamental de soluciones en (0, ∞) y por tanto
la solución general de la ecuación es

y = xα (c1 cos(β ln x) + c2 sen(β ln x))

Ejemplo. Resolver

33
 2
2 − 2x dx − 4y = 0
d y dy
a) x2 dx
2
d y dy
b) 4x2 dx 2 + 8x dx + y = 0

c) 4x2 y ′′ + 17y = 0

d) x2 y ′′ − 3xy ′ + 3y = 2x4 ex

REDUCCION A COEFICIENTES CONSTANTES. Cuando

las raíces de las ecuaciones auxiliares para ay ′′ + by ′ + cy = 0 y
ax2 y ′′ +bxy ′ +cy = 0 son reales y distintas, las soluciones generales
respectivas son:

y = c1 em1 x + c2 em2 x y y = c1 xm1 + c2 xm2 , x > 0

Usando la identidad eln x = x, x > 0 la segunada solución puede
expresarse en la misma forma que la primera

y = c1 em1 ln x + c2 em2 ln x = c1 em1 t + c2 em2 t ,

donde t = ln x. Este último resultado ilustra el hecho de que
cualquier ecuación de Cauchy-Euler se puede escribir como una
ecuación diferencial lineal con coeficientes constantes sustituyendo
x = et .

Ejemplo. Resolver x2 y ′′ − xy ′ + y = ln x.

34
2.6 Ecuaciones diferenciales no lineales

Las ecuaciones diferenciales no lineales de segundo orden F (x, y ′ , y ′′ ) =

0, donde falta la variable dependiente y, y F (y, y ′ , y ′′ ) = 0,donde
falta la variable independiente x, a veces se resuelven usando méto-
dos de primer orden

En una ecuación de la forma F (x, y ′ , y ′′ ) = 0 si hacemos la susti-

tución u = y ′ entonces la ecuación diferencial se convierte en
F (x, u, u′ ) = 0. Si podemos resolver esta última ecuación para
u, podemos encontrar a y por integración.

Ejemplo. Resolver y ′′ = 2x(y ′ )2 .

Si la ecuación tiene la forma F (y, y ′ , y ′′ ) = 0, hacemos u = y ′ , pero

debido a que falta la variable independiente x, esta sustitución se
usa para convertir la ecuación diferencial en una en la que variable
independiente es y y la variable dependiente es u. Entonces uti-
lizando la regla de la cadena para calcular la segunda derivada de
y:

du du dy du
y ′′ = = =u .
dx dy dx dy
En este caso la ecuación de primer orden que debemos resolver es

35
 du
F (y, u, u ) = 0.
dy
Ejemplo. Resolver yy ′′ = (y ′ )2

36