Tema 1

ECUACIONES DIFERENCIALES

1.1. Introducción
Gran parte de los sistemas que nos rodean están sometidos al cambio, por tanto,
es un hecho cotidiano para todos nosotros. Las Matemáticas son muy útiles para
investigar, entre otros, fenómenos como el movimiento de los planetas, la desintegra-
ción de sustancias radiactivas, la velocidad de las reacciones quı́micas y los patrones
meteorológicos. Por otro lado, los biólogos investigan en campos tales como la conta-
minación o la dinámica de poblaciones. Incluso en áreas, aparentemente alejadas de
la Matemáticas, como las Ciencias Polı́ticas o la Medicina, es frecuente que recurran
a los modelos matemáticos, en los cuales la clave está en el cambio.
Muchos de estos modelos se expresan a través de una ecuación diferencial. Si y = f (t)
es una función que relaciona las variables t e y, entonces su derivada
dy
y′ =,
dt
nos indica la tasa de cambio o velocidad de cambio de la variable y con respecto
de la variable t.
Cuando estudiamos un problema del mundo real necesitamos usualmente desarrollar
un marco matemático. Sabemos que el proceso por el que se crea y evoluciona este
marco es la construcción de un modelo matemático, siendo algunos de ellos muy
precisos, especialmente los de la Fı́sica. Sin embargo, otros lo son menos, concreta-
mente los que tratan de problemas de Biologı́a o Ciencias Sociales. No obstante, en
los últimos años los enunciados de estas materias se han vuelto lo suficientemente
precisos como para poder expresarlos matemáticamente.
Un ejemplo de creación de un modelo continuo lo tenemos en la predicción del
tiempo. En teorı́a, si pudiésemos programar en un ordenador todas las hipótesis co-
rrectas, ası́ como los enunciados matemáticos apropiados sobre las formas en que las

3
4 Tema 1 Ecuaciones diferenciales

condiciones climáticas operan, tendrı́amos un buen modelo para predecir el tiempo

mundial. En el modelo del clima global, un sistema de ecuaciones calcula los cambios
que dependen del tiempo, siendo las variables el viento, la temperatura y la hume-
dad, tanto en la atmósfera como en la tierra. El modelo1 puede predecir también las
alteraciones de la temperatura en la superficie de los océanos.
Por todo lo comentado anteriormente, hemos puesto de manifiesto que en los mode-
los matemáticos del mundo real tienen gran importancia el estudio de las ecuaciones
diferenciales. En cualquier lugar donde se lleve a cabo un proceso que cambie conti-
nuamente en relación al tiempo (rapidez de variación de una variable con respecto
a otra), suele ser apropiado el uso de las ecuaciones diferenciales.

EJERCICIO 1

Escribir una ecuación diferencial que describa la situación dada.

1 La cantidad de bacterias en un cultivo crece, en cada momento, a un
ritmo que es proporcional al número de bacterias presentes.

2 Cuando los factores ambientales imponen un lı́mite superior sobre su

tamaño, la población crece a un ritmo que es conjuntamente proporcional
a su tamaño actual y a la diferencia entres u lı́mite superior y su tamaño
actual.

3 La razón a la que las personas oyen hablar sobre un nuevo aumento de

precios es proporcional al número de personas en la ciudad que no han
oı́do hablar al respecto.

4 El ritmo con el que se propaga una epidemia en una comunidad es conjun-

tamente proporcional a la cantidad de residentes que han sido infectados
y al número de residentes propensos a la enfermedad que no han sido
infectados.

5 Si es cierto que en una economı́a estable la velocidad de disminución del

número de personas y, con un salario de por lo menos x euros, es directa-
mente proporcional al número de personas e inversamente proporcional
a su salario, obténgase la ley de Pareto, es decir la expresión de y en
función de x.

1.2. ¿Qué es una ecuación diferencial?

Aunque no sepamos que es una ecuación diferencial, sin embargo estamos fami-
liarizados con el problema de resolver ecuaciones y sistemas de ecuaciones algebrai-
cas. Además, sabemos lo que se entiende por solución de la ecuación, aunque en
1
En el Centro Nacional de Investigación Atmosférica de EEUU tienen un superordenador con
el nombre de CRAY que puede ejecutar un modelo parecido.
1.2 ¿Qué es una ecuación diferencial? 5

ecuaciones polinómicas de grado elevado o en ecuaciones donde aparecen funciones

trascendentes no podamos encontrar su valor exacto.
De manera general, F (x1 , x2 , · · · , xn ) = 0 , siendo F una función vectorial de IRn en
IRm , representa a un sistema de m ecuaciones en las variables x1 , x2 , · · · , xn .
Si utilizamos el lenguaje del cálculo diferencial podemos escribir ecuaciones donde
aparezcan una función y = y(t), definida sobre un cierto intervalo I ⊂ IR, la variable
t, y las derivadas de diferentes órdenes de y. Por ejemplo:
y ′ = 6t + 5 y ′ = 6y
y ′ + 3y + t = 0 (y ′′ )2 + 2ty + sen t = 0 .
Llamemos la atención sobre el hecho de que ya hemos tenido ocasión de estudiar
este tipo de situaciones, concretamente cuando se realizó el estudio de las integrales
indefinidas. En efecto, dada la ecuación y ′ (t) = sen t la idea básica era encontrar
una función y(t) = − cos t + C que cumpla la igual anterior.
Los siguientes ejemplos tratan de mostrar como las ecuaciones diferenciales aparecen
al modelar situaciones muy simples.

EJEMPLO 1.1
Un zoológico planea llevar un león marino a otra ciudad. El animal irá cubierto du-
rante el viaje con una manta mojada. En cualquier tiempo t, la manta perderá hume-
dad debido a la evaporación, a una razón proporcional a la cantidad y(t) de agua
presente en la manta. Inicialmente, la sábana contendrá 40 litros de agua de mar.
Estamos interesados en encontrar una ecuación diferencial que describa este proble-
ma.
Al ser y(t) la cantidad de agua en la manta en el tiempo t, del enunciado deducimos
que la razón de cambio de y(t) (su derivada y ′ (t)), es proporcional a y(t). Enton-
ces y ′ (t) = ky(t), donde la constante de proporcionalidad k es negativa, ya que la
cantidad de agua disminuye con el tiempo. Por tanto, nuestro modelo será
y ′ (t) = ky(t) , k ≤ 0, y(0) = 40 .

EJEMPLO 1.2
La tabla siguiente:

Horas 0 1 3 5 7 9 11 13 15 17 19
Conc.(mg/l) 12.0 10.0 7.0 5.0 3.5 2.5 2.0 1.5 1.0 0.7 0.5

muestra la concentración de teofilina, una droga común para combatir el asma, en

el torrente sanguı́neo, como una función del tiempo después de la aplicación de una
dosis inicial.
Si representamos la concentración de teofilina en función del tiempo nos aparece una
gráfica que disminuye de manera exponencial (Figura 1.1 izquierda)
6 Tema 1 Ecuaciones diferenciales

Figura 2.1. Izquierda: escala normal. Derecha: escala logarı́tmica

Si tomamos logaritmos neperianos (Figura 2.1 derecha) de los valores de la con-

centración, podemos ajustar esta nueva nube de puntos por una recta. Este pro-
ceso lo llevamos a cabo con el programa Mathematicar y su solución es la recta
2.45337 − 0.164264t, que corta al eje de ordenadas en el punto (0, 2.45337) y su pen-
diente es −0.164264. Por lo tanto, si la solución del modelo es del tipo exponencial
y(t) = Cekt , entonces ln y = ln C + kt. En consecuencia,

ln C = 2.45338 ⇒ C = e2.45338 = 11.6275 ; k = −0.164265

Figura 2.2. Izquierda: ajuste lineal. Derecha: ajuste exponencial

y(t) = 11.6275e−0.164264 t

Pasemos ahora a precisar algunos de los conceptos sugeridos.

Una ecuación diferencial es aquella en la que aparece una función desconocida y
una o más de sus derivadas. Cuando la función desconocida depende de dos o más
variables, entonces las derivadas que aparecen en la ecuación diferencial serán deri-
vadas parciales, y en este caso diremos que se trata de una ecuación en derivadas
parciales. Si la función depende sólo de una variable independiente, entonces la
ecuación recibe el nombre de ecuación diferencial ordinaria (E.D.O.). En este
curso estudiaremos algunos tipos de ecuaciones diferenciales ordinarias de orden n
que representaremos por
dy d2 y dn y
F (t, y, , 2 , · · · , n ) = F (t, y, y ′ , · · · , y n) ) = 0 ,
dt dt dt
1.2 ¿Qué es una ecuación diferencial? 7

donde F es una expresión matemática en la que aparecen la variable t, una función

desconocida y, y las derivadas de y hasta el orden n.

EJEMPLO 1.3
Las siguientes ecuaciones son ecuaciones diferenciales ordinarias.
−2y ′′ + 3y ′ − y = et

dy
y ′ (t) = = ay − by 2
dt

d2 y dy
−2 + t + 6y = 0
dt2 dt
Las ecuaciones
∂ 2 u(x, t) 1 ∂ 2 u(x, t)
2
= 2
∂x c ∂t2

∂2u ∂2u ∂2u 1 ∂u

2
+ 2 + 2 = , u = u(x, y, z, t) ,
∂x ∂y ∂z k ∂t
son ejemplos de ecuaciones en derivadas parciales.

El orden de una ecuación diferencial es el que corresponde a la derivada de mayor

orden que aparece en la ecuación. De esta manera, y ′ = ay − by 2 es una ecuación
diferencial ordinaria de primer orden, mientras que
∂ 2 u(x, t) 1 ∂ 2 u(x, t)
= ,
∂x2 c2 ∂t2
es una ecuación en derivadas parciales de segundo orden.

EJEMPLO 1.4
Clasificar cada una de las siguientes ecuaciones diferenciales como ordinarias o en
derivadas parciales. Determinar el orden y la linealidad o no linealidad en cada caso.
(a) y ′ + t2 y = tet (b) y ′′′ + 4y ′′ − 5y ′ + 3y = sen t
∂2u ∂2u
(c) + 2 =0 (d) t2 dy + y 2 dt = 0
∂x2 ∂y
( 2 )5
dy d y ∂4u ∂2u ∂2u
(e) +3 + 5y = 0 (f ) + + 2 +u=0
dt dt2 ∂x2 ∂y 2 ∂x2 ∂y
( )3 √ 2
dr d r
(g) y ′′ + y sen t = 0 (h) = +1
ds ds2
d2 y d2 q dq q
(i) + t sen y = 0 (j) L +R + =0
dt2 √ ( )2
dt 2 dt c
d2 ρ 4 dρ
(k) 2
= ρ+
dθ dθ
8 Tema 1 Ecuaciones diferenciales

Las soluciones son:

(a) y ′ + t2 y = tet
Ecuación diferencial ordinaria lineal de primer orden.
(b) y ′′′ + 4y ′′ − 5y ′ + 3y = sen t
Ecuación diferencial ordinaria lineal de tercer orden.
∂2u ∂2u
(c) + 2 =0
∂x2 ∂y
Ecuación en derivadas parciales de segundo orden.
(d) t2 dy + y 2 dt = 0
Ecuación diferencial ordinaria de primer orden no lineal.
( 2 )5
dy d y
(e) +3 + 5y = 0
dt dt2

Ecuación diferencial ordinaria de segundo orden no lineal.

∂4u ∂2u ∂2u
(f) + + 2 +u=0
∂x2 ∂y 2 ∂x2 ∂y
Ecuación diferencial en derivadas parciales de cuarto orden.
(g) y ′′ + y sen t = 0
Ecuación diferencial ordinaria lineal de segundo orden.
( )3 √ 2
dr d r
(h) = +1
ds ds2

Ecuación diferencial ordinaria no lineal de segundo orden.

d2 y
(i) + t sen y = 0
dt2
Ecuación diferencial ordinaria no lineal de segundo orden.
d2 q dq q
(j) L +R + =0
dt2 dt c
Ecuación diferencial ordinaria lineal de segundo orden.
√ ( )2
d2 ρ 4 dρ
(k) 2
= ρ+
dθ dθ

Ecuación diferencial ordinaria no lineal de segundo orden

1.3 Solución de una ecuación diferencial 9

1.3. Solución de una ecuación diferencial

Antes de desarrollar esta sección consideremos la ecuación x2 − 4x + 3 = 0. Cuan-
do nos planteamos el problema de encontrar soluciones de esta ecuación estamos
suponiendo que existe un conjunto X donde la variable x puede tomar valores. En
general, la ecuación no es válida para todo valor x ∈ X y el problema de resolver la
ecuación consiste en encontrar S ⊂ X tal que x2 − 4x + 3 = 0. Entonces S será el
conjunto de soluciones, que en nuestro caso es {1, 3}, y por tanto decimos que 1 y
3 son soluciones.
DEFINICIÓN 1.3.1 Una solución de la ecuación diferencial
F (t, y, y ′ , · · · , y n) ) = 0,
es cualquier función y = φ(t), definida en un cierto intervalo I ⊂ IR, con derivada
de orden n en ese intervalo y tal que
F (t, φ(t), φ′ (t), · · · , φ(t)n) ) = 0 , ∀t ∈ I .
El proceso de determinar todas las funciones que son soluciones de una ecuación
diferencial se llama resolver la ecuación diferencial. Por ejemplo, la integración es
un tipo muy simple de resolución de ecuaciones diferenciales.
A diferencia de las ecuaciones algebraicas, las ecuaciones diferenciales tienen por so-
lución una función. Además, una ecuación diferencial tiene generalmente un número
infinito de soluciones que recibe el nombre de solución general. Algunas ecuacio-
nes diferenciales tienen soluciones que no pueden obtenerse de la solución general y
en este caso reciben el nombre de soluciones singulares.
En ocasiones, se desea encontrar una solución particular que satisfaga ciertas condi-
ciones adicionales llamadas condiciones iniciales. Las condiciones iniciales especifican
los valores de una solución y de cierto número de sus derivadas en un valor concreto
de la variable t (con frecuencia es t = 0, pero puede ser cualquier otro). El problema
de determinar una solución de una ecuación diferencial que satisfaga ciertas condi-
ciones iniciales se llama un problema de valores iniciales o de Cauchy.

EJEMPLO 1.5
La ecuación diferencial (y ′ (t))2 + 1 = 0 no tiene solución real, ya que no existe un
número real que elevado al cuadrado y sumado con uno valga cero.

La ecuación t2 + y 2 − 4 = 0 define en forma implı́cita una solución de la ecuación

diferencial t + yy ′ = 0 en el intervalo −2 < t < 2. En efecto, si derivamos en forma
implı́cita la expresión t2 + y 2 − 4 = 0 obtenemos,
2t + 2yy ′ = 0 ⇒ t + yy ′ = 0 .
√
Si despejamos en la solución el valor de y observamos que y = ± 4 − t2 sólo está de-
finida en −2 < t < 2.
10 Tema 1 Ecuaciones diferenciales

Si derivamos la función {
−t4 si t<0
y=
t4 si t≥0
podemos comprobar que es solución de la ecuación diferencial ty ′ − 4y = 0 en el
intervalo −∞ < t < ∞.

EJEMPLO 1.6

Estudiar si la función y = 1/t es una solución de la ecuación y′ = −y2 en el intervalo

(0, +∞).
La función y = 1/t es derivable en el intervalo (0, +∞) y su derivada viene dada por
y ′ = −1/t2 . Por lo que resulta inmediato que la función y = 1/t satisface la ecuación
diferencial y ′ = −y 2 .

1.3.1. Existencia y unicidad de soluciones

Una vez que sabemos lo que se entiende por ecuación diferencial y solución de la
misma, podemos preguntarnos:

¿Toda ecuación diferencial tiene solución?

En el caso de que ésta exista, ¿cuántas tiene?, ¿quiénes son?

Antes de responder a estas preguntas, veamos el ejemplo siguiente:

EJEMPLO 1.7

La ecuación diferencial (y ′ (t))2 + (y(t))2 + 1 = 0 no tiene solución ya que (y ′ (t))2 +

(y(t))2 ≥ 0 para cualquier pareja de valores reales que tomen las funciones y ′ (t) e
y ′ (t).

Es inmediato comprobar que

y(t) = t3 + C , ∀t ∈ IR ,

es solución de la ecuación diferencial y ′ (t) = 3t2 , para cualquier valor de la constante

C. Por tanto, existe un número infinito de soluciones.

En cuanto a la ecuación y ′′ (t) = 0, cualquier función cuya gráfica sea una recta
será solución. También en este caso existe un número infinito de soluciones.
1.3 Solución de una ecuación diferencial 11

Es bastante corriente que si una ecuación diferencial tiene solución, tenga infinitas
soluciones. En efecto, en el proceso de resolver la ecuación diferencial tenemos que
hacer al menos una integral y en consecuencia nos aparecerá una constante que, al
tomar diferentes valores, nos definirá una gama infinita de soluciones.
A partir de este momento, y salvo que no lo indiquemos, nos centraremos en las
ecuaciones diferenciales de primer orden F (t, y, y ′ ) = 0, donde supondremos que
podemos expresarlas como y ′ = f (t, y).
Consideremos el problema de valores iniciales (P.V.I.):

y ′ = f (t, y) , y(t0 ) = y0 , (1.1)

estamos interesados en saber si dicho problema tiene solución y en caso afirmativo

si ésta es única.

EJEMPLO 1.8

Es fácil ver que la ecuación diferencial ty ′ + y = 1 admite como solución general

y = c/t + 1 , c ∈ IR, en cualquier intervalo que no contenga al cero. En efecto,
derivando la función y(t) se tiene
c c c
y′ = − ⇒ ty ′ + y = − + = 1 .
t2 t t
Si queremos determinar la solución que pasa por el punto (1, 2) tenemos que imponer
la condición y(1) = 2. El valor de c que cumple con el requisito anterior es c = 1,
con lo cual la solución particular pedida es
1
y= + 1. (1.2)
t
En consecuencia, la función (1.2) es una solución del problema de valores iniciales

ty ′ + y = 1 , y(1) = 2 , (1.3)

en el intervalo (0, +∞). Puesto que en

c
y= + 1, c ∈ IR ,
t
están todas las soluciones de la ecuación diferencial ty ′ + y = 1, entonces el problema
(1.3) tiene solución única.
En cambio, no es posible encontrar una solución que pase por el punto (0, 2). Por
tanto, en este caso diremos que el problema de valores iniciales

ty ′ + y = 1 , y(0) = 2 , (1.4)

no tiene solución.
12 Tema 1 Ecuaciones diferenciales

Es inmediato comprobar que el problema de valores iniciales

(y ′ )2 = 4y , y(0) = 1 , (1.5)

tiene dos soluciones: (a) y = (t − 1)2 , (b) y = (t + 1)2 .

TEOREMA 1.3.2 (Teorema de Cauchy-Peano) Sea (t0 , y0 ) ∈ IR2 y suponga-

mos que existe un rectángulo cerrado

R = {(t, y) ∈ IR2 : |t − t0 | ≤ a , |y − y0 | ≤ b}

en el que la función f es continua. Entonces el problema de valores iniciales

y ′ = f (t, y) , y(t0 ) = y0 ,

tiene al menos una solución definida en el intervalo (t0 − δ, t0 + δ), donde

( )
b
δ = mı́n a, , M = máx 2 |f (t, y)| . (1.6)
M (t,y)∈IR

Hemos visto en el teorema de Cauchy-Peano que la continuidad de la función f (t, y)

en una región R garantiza que por cada punto de R pasa una solución de la ecuación
diferencial y ′ = f (t, y) ¿Será también cierto que la continuidad de la función f (t, y)
obliga a que por cada punto de R pase una única solución? El siguiente ejemplo nos
dará la respuesta a esta pregunta.

EJEMPLO 1.9
2
Supongamos la ecuación diferencial y ′ = f (t, y) = y 3 , que podemos escribirla

d
y′y− 3 = 1
2 1
⇒ (3y 3 ) = 1 .
dt
Integrando
( )3
1 t
3y 3 = t + c ⇒ y= +k , k = cte.
3
El problema de valores iniciales
2
′
y = y3 , y(0) = 0 ,

no tiene solución única, ya que y = t3 /27, e y = 0 son dos soluciones del mismo.
1.3 Solución de una ecuación diferencial 13

Este ejemplo muestra una ecuación diferencial con una función f (t, y) = y 2/3 con-
tinua en un rectángulo R que contiene al punto (0, 0), y sin embargo no tiene una
única solución. Si queremos conseguir este último objetivo será necesario exigir a la
función f nuevas condiciones.
TEOREMA 1.3.3 (Teorema de Picard) Sea (t0 , y0 ) ∈ IR2 y supongamos que
existe un rectángulo cerrado

R = {(t, y) ∈ IR2 : |t − t0 | ≤ a , |y − y0 | ≤ b}
∂f
en el que las funciones f y son continuas. Entonces el problema de valores
∂y
iniciales
y ′ = f (t, y) , y(t0 ) = y0 ,
tiene solución única definida en el intervalo (t0 − δ, t0 + δ), donde δ está dado por
(1.6).

OBSERVACIÓN 1.3.4

Los Teoremas 1.3.2 y 1.3.3 nos dan condiciones suficientes pero no necesarias
para garantizar la existencia y unicidad de soluciones para un problema de
valores iniciales.

La solución de un problema de valores iniciales puede existir en un intervalo

mayor que el mencionado en los Teoremas 1.3.2 y 1.3.3.

En los Teoremas 1.3.2 y 1.3.3 hemos considerado rectángulos R cerrados y

acotados. Pueden enunciarse teoremas análogos utilizando rectángulos abiertos

R = {(t, y) ∈ IR2 : |t − t0 | < a , |y − y0 | < b} ,

o bien rectángulos del tipo

R = {(t, y) ∈ IR2 : t0 ≤ t < t0 + a , |y − y0 | < b} .

∂f
En estos casos tenemos que añadir la hipótesis de que las funciones f y
∂y
estén acotadas.

EJEMPLO 1.10

En el problema de valores iniciales

ty ′ = 2y , y(0) = 0 , (1.7)

las funciones
2y ∂f (t, y) 2
f (t, y) = , = ,
t ∂y t
14 Tema 1 Ecuaciones diferenciales

no están definidas en los puntos de la recta t = 0. Por tanto, no es posible encontrar

un rectángulo R que contenga al punto (0, 0) en el cual la función f sea continua.
No podemos aplicar el Teorema 1.3.2 y, en consecuencia, no podemos asegurar nada
sobre la existencia de solución del problema de valores iniciales (1.7). Sin embargo,
es fácil ver que las funciones y = c t2 con c ∈ IR son soluciones del problema (1.7) en
el intervalo −∞ < t < ∞. El problema de valores iniciales tiene infinitas soluciones.

1.4. Análisis geométrico de y ′ = f (t, y)

Recordemos que estamos considerando ecuaciones diferenciales ordinarias de primer
orden F (t, y, y ′ ) = 0, donde F es una función de tres variables, y que será posible
expresarla
y ′ (t) = f (t, y) . (1.8)

1.4.1. Campo de direcciones

Las soluciones de (1.8) son funciones y las podemos representar gráficamente como
una curva en el plano Oty. Supongamos que D sea el dominio de la función f , y
(t0 , y0 ) ∈ D siendo y(t) una solución de (1.8) de tal manera que su gráfica pasa por
el punto (t0 , y0 ), por tanto y(t0 ) = y0 . En consecuencia, la ecuación (1.8) expresa
que f (t0 , y0 ) es el valor de la pendiente de la tangente a la gráfica de y(t) en (t0 , y0 ).
De esta manera, para cada uno de los puntos del dominio D podemos dibujar un
pequeño segmento con la dirección que f (t, y) determina. Un subconjunto del plano
O t y en el que para cada punto se ha definido una dirección se conoce con el nombre
de campo de direcciones.
Resumiendo, lo que hemos hecho al plantear la ecuación (1.8) es definir un campo
direccional y el problema de encontrar sus soluciones es el de encontrar aquellas
curvas con la propiedad de ser tangentes a cada punto del campo de direcciones.

EJEMPLO 1.11

Para dibujar el campo de direcciones y poder trazar algunas de las soluciones de la

ecuación diferencial y ′ = y 2 , veamos qué información podemos extraer de nuestra
ecuación diferencial.

(a) Es evidente que para cualquier valor de y su derivada y ′ es positiva. Por tanto,
las curvas solución son crecientes
(b) Para estudiar la concavidad de las curvas solución necesitamos su segunda
derivada y ′′ = 2yy ′ = 2y 3 . En consecuencia, si y > 0, las curvas solución son
convexas, mientras que si y < 0 son cóncavas.
1.4 Análisis geométrico de y ′ = f (t, y) 15

(c) Campo de direcciones. Nuestra ecuación diferencial define un campo de di-

recciones en todo el plano Oty cuyas direcciones son constantes a lo largo de
rectas paralelas al eje de abscisas t.
Podemos construirlo (véase Figura 2.3) con ayuda del Mathematicar .
<< Graphics‘PlotField‘
PlotVectorField[{1, y2 }, {t, −3, 3}, {y, −3, 3}]
Como la dirección que define el campo de la ecuación y ′ = y 2 , en cada pun-
to del plano depende sólo de la coordenada y, entonces para cualquier y0 los
puntos de la forma (t, y0 ) con t ∈ IR, se encuentran rodeados de un campo
direccional idéntico. En consecuencia, las soluciones pueden obtenerse una de
otra haciendo traslaciones en la dirección del eje t

Figura 2.3. Campo de direcciones de y ′ = y 2 .

(d) Este hecho puede comprobarse si encontramos la solución explı́cita de la ecua-

ción diferencial. Es inmediato comprobar que y(t) = −1/(t + c). Observemos
que esta familia de soluciones no contiene la solución y = 0 para cualquier c
finita.

Para este ejemplo ha sido muy fácil encontrar la solución de la ecuación diferencial,
pero esto no es lo más frecuente. Por tanto, en gran parte de los casos será nece-
sario hacer un estudio geométrico para conocer, al menos, el comportamiento de
las soluciones. Tengamos en cuenta que en muchos de los modelos que analizaremos
estaremos interesados no en la solución concreta del problema, sino en su compor-
tamiento a “largo plazo”.
16 Tema 1 Ecuaciones diferenciales

1.5. Teorı́a cualitativa de EDO autónomas

1.5.1. Introducción.
A finales del 1600 I. Newton y G. Leibnitz descubrieron el Cálculo y pusieron las
bases para el estudios de los Sistemas Dinámicos. En un principio y hasta momen-
tos recientes se ha intentado encontrar de forma exacta la solución de la ecuación
diferencial que modeliza a una determinada situación. Sin embargo, existen modelos
aparentemente sencillos donde esto no es posible, por ejemplo el problema propuesto
a finales del siglo XIX por Poincaré2 conocido con el nombre de los tres cuerpos.
Los matemáticos probaron que para este problema de atracción gravitatoria no era
posible dar su solución explı́cita.
Por tanto, el desarrollo histórico de las ecuaciones diferenciales ha seguido dos cami-
nos diferentes. El primero, se caracteriza por una búsqueda de soluciones explı́citas,
bien sea en fórmulas exactas (lo que rara vez es posible) o bien en términos de series
de potencias. En el segundo, se abandona toda intención de resolver las ecuaciones
diferenciales en sentido tradicional y se intenta obtener información cualitativa sobre
el comportamiento general de las soluciones.
En esta sección realizaremos un estudio geométrico para obtener información so-
bre el comportamiento de las soluciones de ciertas ecuaciones diferenciales llamadas
autónomas. En las próximas secciones estudiaremos la forma de resolver ciertos ti-
pos de ecuaciones diferenciales. En general, resolver una ecuación diferencial es un
problema difı́cil, sin embargo, en muchas ocasiones es posible dar información sobre
las soluciones sin necesidad de calcularlas.

1.5.2. Ecuaciones diferenciales autónomas

Ahora, nos centraremos en el problema de aprender cuanto sea posible sobre las
caracterı́sticas esenciales de las soluciones de ecuaciones diferenciales de la forma
y ′ = g(y) por análisis directo de la propia ecuación. Este tipo de ecuaciones dife-
renciales recibe el nombre de autónomas pues el segundo miembro de la ecuación
es “independiente del tiempo”, en el sentido de no aparecer t. Además, si y(t) es
solución de una ecuación autónoma también lo es la función y(t + c), para cualquier
constante c.

DEFINICIÓN 1.5.1 Los puntos c ∈ IR tales que y(t) = c es solución de la ecua-

ción diferencial se llaman puntos de equilibrio.

Si suponemos que el comportamiento dinámico de un sistema biológico está mode-

lado matemáticamente por las curvas solución de una ecuación diferencial autónoma
2
A.H. Poincaré (1854 - 1912) se le consideró como el matemático más grande de su época.
Fundó la dinámica topológica y la topologı́a. En sus trabajos sobre la mecánica celeste elaboró la
teorı́a de los desarrollos asintóticos, la cual, en la actualidad es una de las herramientas más
poderosa del matemático aplicado
1.5 Teorı́a cualitativa de EDO autónomas 17

y estamos interesados por el comportamiento a largo plazo de las trayectorias (es

decir, de las curvas solución), son de especial interés los estados de equilibrio, que son
aquellos estados y(t) que no cambian con el tiempo. Matemáticamente esto significa
que y(t) = c es una solución de la ecuación y ′ = g(y).

EJEMPLO 1.12

Dada la ecuación diferencial y ′ = 7.5 − 24.25y + 22.25y 2 − 8y 3 + y 4 . Para encontrar

los puntos de equilibrio resolvemos la ecuación y ′ = 0 y obtenemos y = 0.5, 2, 2.5, 3.
Por tanto, las funciones y = 0.5, y(t) = 2, y = 2.5, y(t) = 3 son soluciones cons-
tantes. Por el Teorema 1.3.2 y Teorema 1.3.3 sabemos que la solución es única. En
consecuencia, ninguna otra solución puede tomar los valores 0.5, 2, 2.5 o 3. De este
modo, el plano Oty quedará dividido en regiones horizontales de tal manera que una
solución que comience en una región no podrá salir de ella.

Figura 2.4. Lı́nea fase de y ′ = (y − 0.5)(y − 2)(y − 2.5)(y − 3).

Si la condición inicial y0 es menor de 0.5 tendremos que y ′ (t0 ) es positiva y la

solución será creciente. Si y está entre 0.5 y 2 o entre 2.5 y 3, entonces la derivada
será negativa y la función decrecerá . Finalmente, si una solución comienza entre 2 y
2.5 o se encuentra por encima de 3 será creciente. En general se cumple la siguiente
propiedad.

RESULTADO 1.5.2 Si g es una función con derivada continua en todo IR y con-

sideramos la ecuación diferencial y ′ = g(y). Entonces:
Para cada una de las raı́ces de g(y) = 0, existe una solución constante de la
ecuación diferencial. Si g(c) = 0 entonces y = c es una solución.

Las soluciones constantes dividen al plano Oty en franjas horizontales. Cual-

quier otra solución no constante estará contenida en una franja y será estric-
tamente creciente o estrictamente decreciente.
18 Tema 1 Ecuaciones diferenciales

Cada solución no constante es asintótica a una solución constante, o bien,

crece o decrece sin lı́mite.

En nuestro ejemplo, observamos que si la condición inicial está próxima al 0.5 ó 2.5,
entonces se tiene que la solución del problema de valores iniciales tiende a 0.5 ó 2.5
cuando t tiende a infinito. Por el contrario, si la condición inicial está próxima al 3
pero sin serlo, entonces la solución del problema de valores iniciales crece sin lı́mite
o decrece hacia 2.5. De alguna manera las soluciones constantes 0.5 y 2.5 atraen a
las soluciones mientras que las soluciones constantes 2 y 3 las repelen.
Las ideas anteriores conducen a los conceptos de estabilidad e inestabilidad. Ası́,
las soluciones y(t) = 0.5 e y = 2.5 son estables mientras que y(t) = 2 o y(t) = 3
tienen un comportamiento inestable.
Intuitivamente, desde un punto de vista fı́sico solo interesan los puntos de equilibrio
que son “estables ”. Un péndulo en la posición vertical superior está en equilibrio,
pero es muy improbable que eso ocurra. Además, la menor perturbación altera-
rá completamente el comportamiento del péndulo. Tal equilibrio es inestable. En
cambio, la posición de reposo inferior es estable; si la perturbamos ligeramente, el
péndulo oscilará a su alrededor y (a causa del rozamiento) se aproximará gradual-
mente a ella de nuevo. De aquı́ nace la idea intuitiva de fuente y sumidero.

DEFINICIÓN 1.5.3 Decimos que un punto de equilibrio y0 es:

Un sumidero si cualquier solución con condición inicial “suficientemente cercana.a

y0 es asintótica a y0 cuando t aumenta.

Una fuente, cuando todas las soluciones que comienzan cerca de y0 se alejan
de y0 cuando t aumenta.

Un nodo si no es fuente o sumidero.

En nuestro caso, el eje de ordenadas recibe el nombre de lı́nea fase, siendo los
puntos 0.5 y 2.5 sumideros y los puntos 2 y 3 fuentes.
Por lo comentado anteriormente, si c es un punto de equilibrio y g ′ (c) < 0 entonces
el cambio de signo es de positivo a negativo y las condiciones iniciales justo por
debajo de c dan lugar a funciones crecientes hacia c y las por encima de c funciones
decrecientes a la solución constante. En el caso en que g ′ (c) = 0 no podemos asegurar
nada y es necesario ver si se produce cambio de signo. Si no se produce cambio de
signo tendremos que las soluciones por encima y por debajo de la constante son
ambas crecientes o decrecientes, es decir, por un lado se alejarán de la solución
constante y por otro se acercarán.
1.5 Teorı́a cualitativa de EDO autónomas 19

RESULTADO 1.5.4 En general, se cumple:

Si g(a) = 0 y g ′ (a) < 0, entonces a es un estado de equilibrio estable para la
ecuación diferencial autónoma y ′ = g(y).
Si g(a) = 0 y g ′ (a) > 0, entonces implica que a es un estado de equilibrio
inestable para la ecuación diferencial autónoma y ′ = g(y).
Si a es un estado de equilibrio para y ′ = g(y) en el cual g ′ (a) = 0, debemos
estudiar la situación con más cuidado. Podemos encontrar ejemplos donde a
sea estable o inestable.

EJEMPLO 1.13
En el estudio de los efectos de la selección natural sobre una población aparece la
siguiente ecuación diferencial,

y ′ (t) = 0.01y 2 (t)(1 − y(t)) (1.9)

donde y(t) representa a la frecuencia con que se presenta cierto gen a, ¿contrá quien
va la presión selectiva?
Para conocer el comportamiento a largo plazo del modelo bastará con realizar un
estudio cualitativo de la ecuación diferencial autónoma (1.9)) y para ello será nece-
sario encontrar y clasificar sus puntos de equilibrio.
Los puntos de equilibrio son las soluciones y(t) constantes, por tanto aquellas fun-
ciones donde y ′ (t) = 0, es decir y(t) = 1, y(t) = 0.
Las soluciones constantes dividen a la región {(t, y) ∈ IR2 /t ≥ 0 , y ≥ 0} en dos
franjas (Figura 2.5 colores verde y amarillo). Para valores iniciales de y(t) pertene-
cientes a la primera región 0 < y(t) < 1, la derivada es positiva y en consecuencia las
soluciones y(t) son crecientes. Sin embargo, en la segunda región 1 < y(t) (aunque
sin sentido biológico) la derivada y ′ (t) es negativa lo que indica que las funciones
soluciones y(t) son decrecientes. Estos resultados nos permiten decir que el punto
de equilibrio y(t) = 0 es inestable, mientras que y(t) = 1 es asintóticamente estable
(sumidero). A largo plazo, y para cualquier valor inicial 0 < y(0) < 1 las soluciones
y(t) → 1

Figura 2.5.
20 Tema 1 Ecuaciones diferenciales

EJERCICIO 2

1 En el estudio de los efectos de la selección natural sobre una población

aparecen las siguientes ecuaciones diferenciales,

y ′ (t) = y(t)(1 − y(t))(0.15 − 0.5y(t))

y ′ (t) = 0.05y(t)(1 − y(t))(2y(t) − 1)

donde y(t) representa a la frecuencia con que se presenta cierto gen a.

Trazar las soluciones representativas considerando distintas condiciones
iniciales entre 0 y 1 y discutir posible interpretaciones genéticas para
estas curvas.

2 Obtener y clasificar los puntos de equilibrio de las ecuaciones diferenciales

autónomas.
y ′ (t) = (1 − y)(y + 1)2

y ′ (t) = y(y − 1)(8y − 2)

(y)
y ′ (t) = sen 2

3 La dinámica de una población viene dada por el siguiente modelo

( )( )
dy(t) y(t) y(t)
= 0.25 −1 1−
dt 10 200

donde y(t) representa al número de individuos en el tiempo t.

(a) Encuentra los valores de y(t) para que la población se encuentre en

equilibrio.
(b) Encuentra los valores de y(t) para los que decrece la población.

1.6. Resolución de E.D.O. de primer orden

1.6.1. Ecuaciones diferenciales en variables separables
Una importante clase de ecuaciones diferenciales está formada por aquellas que
pueden expresarse de la forma: y ′ = p(t)q(y), donde p(t) es una función únicamente
de la variable t y q(y) es una función únicamente de la variable y.
1.6 Resolución de E.D.O. de primer orden 21

Si y ′ = p(t)q(y) entonces (si q(y) ̸= 0) dividimos por q(y) e integramos respecto de

t, obteniendo: ∫ ∫
1
dy = p(t)dt .
q(y)

EJEMPLO 1.14

Si deseamos resolver
dy
= y cos t , y(π/2) = 1.
dt
Estamos ante una ecuación diferencial de variables separables. Procediendo tal y
como hemos comentado anteriormente llegamos a
dy
= cos t dt, (y ̸= 0).
y
Calculamos cada una de estas dos integrales
∫ ∫
dy
= cos t dt ⇒ ln |y| = sen t + c , c ∈ IR.
y
O bien

|y| = esen t+c = esen t ec ⇒ y = k esen t , k ∈ IR \ {0} (k = ± ec ), (1.10)

Observemos que hemos podido separar las variables cuando y era distinto de cero.
No obstante, es inmediato comprobar que la función y = 0 también es solución
de la ecuación diferencial. Dicha solución también podemos obtenerla de (1.10), si
admitimos que k pueda tomar el valor 0. En consecuencia, la solución general de la
ecuación diferencial viene dada por

y = k esen t , k ∈ IR. (1.11)

Ahora, si deseamos conocer la solución particular que pasa por el punto (π/2, 1),
sustituimos los valores en (1.11),

y(π/2) = 1 ⇒ 1 = k esen π/2 ⇒ k = 1/e.

La solución del problema de valores iniciales vendrá dada por

y = esen y−1

EJEMPLO 1.15

En ciertas situaciones se plantea determinar la relación entre algún estı́mulo fı́sico y

la reacción correspondiente que se produce en el sujeto. Supongamos que la fuerza
de un estı́mulo es s y que la intensidad de la reacción es una función de s, f (s).
Algunos datos experimentales sugieren que la razón de cambio de la intensidad de
la reacción con respecto al estı́mulo es directamente proporcional a la intensidad de
la reacción e inversamente proporcional a la fuerza del estı́mulo.
22 Tema 1 Ecuaciones diferenciales

De los comentarios anteriores se desprende que f (s) satisface la ecuación diferencial

f (s)
f ′ (s) = k
s
para alguna constante positiva k. Es inmediato comprobar que la solución general
de esta ecuación diferencial de variables separables viene dada por

f (s) = c sk

EJEMPLO 1.16

La tasa de variación de una población de bacterias viene dada por la ecuación dife-
rencial y ′ (t) = (1 − t)y(t), siendo y(t) el número de bacterias en el minuto t. Si
inicialmente el número de bacterias es y0 , ¿cuántas bacterias habrá después de t
minutos?
La ecuación diferencial es de variables separadas
( )
dy t2 t2
= (1 − t)dt ⇒ ln y = t − +C ⇒ y = ket− 2
y 2

Ahora encontramos la solución particular correspondiente al valor y(0) = y0 , es decir

k = y0 . Por tanto
t2
y(t) = y0 et− 2

EJEMPLO 1.17

Resolver las siguientes ecuaciones diferenciales.

dy et ey t
(1) = ; (2) y′ =
dt 2y ey + t2 ey

y ′ + y = y(tet + 1),
2
(3) y(0) = 1

dy et
(a) = Se trata de una ecuación de variables separables,
dt 2y

dy et
= ⇒ 2y dy = et dt,
dt 2y

que se resuelve integrando en ambos términos de la ecuación

∫ ∫
2y dy = et dx ⇒ y 2 = et + c , c ∈ IR

ey t
(b) y′ =
ey + x2 ey
1.6 Resolución de E.D.O. de primer orden 23

Simplificando se reduce a una ecuación diferencial inmediata

ey t ey t t
y′ = y 2 y
= y 2
= ,
e +t e e (1 + t ) 1 + t2
que se resuelve por integración,
∫ ∫
t 1 2t 1
y= 2
dt = 2
= ln(1 + t2 ) + c, c ∈ IR.
1+t 2 1+t 2

y′ + y = y(tet + 1),
2
(c) y(0) = 1

Simplificando la expresión, la ecuación diferencial se reduce a una de variables

separables,
dy dy
y ′ = y t et ⇒
2 2 2
= y t et ⇒ = t et dt, (y ̸= 0).
dt y
Integrando en ambos términos, se obtiene
∫ ∫
dy 2 1 2
= t et dt, ⇒ ln |y| = et + c, c ∈ IR,
y 2
que puede expresarse en forma explı́cita como
1 t2
y = k e2e , k ∈ IR \ {0} (k = ± ec ). (1.12)

La división por y al separar las variables nos lleva a considerar la función

y = 0 que también resulta ser solución de la ecuación diferencial. Dicha so-
lución se obtiene de (1.12) si admitimos el valor k = 0. La solución general
vendrá definitivamente dada por
1 t2
y = k e2e , k ∈ IR. (1.13)

Para determinar la solución particular que verifica la condición inicial y(0) = 1,

sustituimos los valores t = 0, y = 1 en (1.13),

y(0) = 1 ⇒ 1 = k e1/2 ⇒ k = e−1/2 .

Sustituyendo en (1.13) obtenemos la solución

( 2
)
1
e7 −1
y=e 2
.

1.6.2. Ecuaciones diferenciales exactas.

Una forma de obtener una ecuación diferencial es suponer F (t, y) = C y calcular su
diferencial total. En efecto,
∂F (t, y) ∂F (t, y)
dt + dy = 0 . (1.14)
∂t ∂y
24 Tema 1 Ecuaciones diferenciales

Es frecuente encontrarnos con ecuaciones diferenciales escritas en la forma

M (t, y) dt + N (t, y) dy = 0 ,

y por comparación con (1.14), podemos preguntarnos si existirá una función F (t, y)
tal que
∂F (t, y) ∂F (t, y)
= M (t, y) , = N (t, y) .
∂t ∂y
Es un hecho conocido (Teorema de Schwartz) que si la función F (t, y) es “razo-
nablemente buena”, entonces sus derivadas cruzadas coinciden. En consecuencia,
tenemos una condición necesaria
∂ 2 F (t, y) ∂ 2 F (t, y) ∂M ∂N
= ⇒ = . (1.15)
∂t∂y ∂y∂t ∂y ∂t

Puede demostrarse, que esta condición también es suficiente.

DEFINICIÓN 1.6.1 Diremos que la ecuación diferencial

M (t, y) dt + N (t, y) dy = 0 ,

es exacta, si cumple
∂M ∂N
= .
∂y ∂t
Si la ecuación diferencial es exacta, entonces
∫
∂F (t, y)
= M (t, y) ⇒ F (t, y) = M (t, y)dt + φ(y) .
∂t
Ahora, podemos derivar respecto de la variable y
[∫ ]
∂F (t, y) ∂
= M (t, y)dt + φ′ (y) = N (t, y) .
∂y ∂y

En consecuencia, [∫ ]
′ ∂
φ (y) = N (t, y) − M (t, y)dt .
∂y
Integramos respecto de y para encontrar el valor de φ(y). Finalmente, la solución
de la ecuación diferencial es F (t, y) = c.

EJEMPLO 1.18

Para la ecuación diferencial (6ty + 2y 2 − 5)dt + (3t2 + 4ty − 6)dy = 0 se tiene

M (t, y) = 6ty + 2y 2 − 5, N (t, y) = 3t2 + 4ty − 6,

1.6 Resolución de E.D.O. de primer orden 25

y puesto que
∂M ∂N
= 6t + 4y = ,
∂y ∂t
es exacta. Por tanto, existirá una función F (t, y) tal que

∂F ∂F
= M (t, y), = N (t, y).
∂t ∂y

Aplicando la técnica de resolución expuesta anteriormente

∫ ∫
∂F
= M (t, y) ⇒ F (t, y) = M (t, y) dt = (6ty + 2y 2 − 5)dt
∂t
= 3t2 y + 2ty 2 − 5t + φ(y),
∂F
= N (t, y) ⇒ 3t2 + 4ty + φ′ (y) = 3t2 + 4ty − 6.
∂y

Entonces, ∫
φ′ (y) = −6 ⇒ φ(y) = −6 dy = −6 y.

La función F (t, y) será: F (t, y) = 3t2 +2ty 2 −5t−6y, y la solución general vendrá dada
en forma implı́cita por 3t2 + 2ty 2 − 5t − 6y = c, c ∈ IR .

Factor integrante

A veces podemos encontrarnos con ecuaciones diferenciales

M (t, y)dt + N (t, y)dy = 0 (1.16)

que no son exactas, pero es posible buscar una función µ(t, y) tal que la ecuación

µ(t, y)M (t, y)dt + µ(t, y)N (t, y)dy = 0 ,

si sea exacta. En este caso, la función µ(t, y) recibe el nombre de factor integrante
de la ecuación (1.16). Notemos que un método para encontrar la función µ(t, y) es
resolver la siguiente ecuación diferencial en derivadas parciales

∂(µ M ) ∂(µ N )
= ,
∂t ∂y

problema que como podemos comprender es bastante complejo. Por esta razón lo
que se hace es simplificarlo. Por ejemplo suponer que la función µ depende solo de
t, sólo de y, o bien de ty, de t + y, etc.
26 Tema 1 Ecuaciones diferenciales

EJEMPLO 1.19

La ecuación diferencial
(t + t4 + t4 y 2 )dt + ydy = 0 (1.17)
no es exacta, ya que
∂M ∂N
M (t, y) = t + t4 + t4 y 2 ⇒ = 2t4 y , N (t, y) = y ⇒ = 0.
∂y ∂t

Si multiplicamos la ecuación (1.17) por 1/(t2 + y 2 ) se obtiene la ecuación

( )
t 2 y
+ t dt + 2 dy = 0 . (1.18)
t2 + y 2 t + y2

Esta nueva ecuación diferencial es exacta. En efecto,

( )
t ∂M1 −2ty
M1 (t, y) = 2 2
+t 2
⇒ = 2
t +y ∂y (t + y 2 )2

y ∂N1 −2ty
N1 (t, y) = ⇒ = 2
t2 + y2 ∂t (t + y 2 )2

Por tanto, la función

1
µ(t, y) =
t2 + y2
es un factor integrante de la ecuación diferencial (1.17). Ahora podemos resolver la
ecuación diferencial exacta (1.18). Es decir, existe una función F (t, y) tal que

∂F t ∂F y
= M1 (t, y) = 2 + t2 , = N1 (t, y) = 2 .
∂t t + y2 ∂y t + y2
Operando
∫
∂F y y 1
= 2 ⇒ F (t, y) = dy = ln(t2 + y 2 ) + φ(t) .
∂y t + y2 t2 +y 2 2
Por otro lado
∂F t t t
= 2 + t2 ⇒ + φ′ (t) = 2 + t2 ,
∂t t + y2 t2 +y 2 t + y2
es decir
1
φ′ (t) = t2 ⇒
φ(t) = t3 .
3
La solución general de (1.17) vendrá dada en forma explı́cita por
1 1
ln(t2 + y 2 ) + t3 = c , c ∈ IR .
2 3
1.6 Resolución de E.D.O. de primer orden 27

1.6.3. Ecuaciones diferenciales lineales de primer orden

La teorı́a de ecuaciones diferenciales lineales ha sido objeto de profundos estudios a
lo largo de los últimos 200 años y es un campo muy bien conocido y muy completo.
Por el contrario, no se sabe casi nada de carácter general acerca de las ecuaciones
diferenciales no lineales.

DEFINICIÓN 1.6.2 Una ecuación diferencial lineal de primer orden es una ecua-
ción del tipo
y ′ + p(t)y = q(t) . (1.19)
Su ecuación homogénea asociada es

y ′ + p(t)y = 0 . (1.20)

TEOREMA 1.6.3 El problema de valores iniciales con una ecuación diferencial

lineal de primer orden tiene solución y es única si las funciones p(t) y q(t) son
continuas.

La resolución de la ecuación homogénea (1.20) es fácil pues es una ecuación de

variables separables y su solución es de la forma
∫
− p(t)dt
y = ce .

Para la resolución de la ecuación lineal completa, se utiliza un método llamado va-

riación de las constantes que consiste en tomar la solución general de la ecuación
homogénea e imponerla como solución de la ecuación completa haciendo depender
de t a la constante c de integración.
Existe un segundo método de resolución que consiste en encontrar el factor
integrante ∫
p(t) dt
µ(t) = e .
Multiplicando la ecuación diferencial por µ(t), se obtiene

µ(t)y ′ (t) + µ(t)p(t)y(t) = µ(t)q(t) ,

que puede expresarse como

∫
′
(µ(t)y(t)) = µ(t)q(t) ⇒ µ(t)y(t) = µ(t)q(t)dt + c .

Tan sólo queda despejar el valor de y(t) para encontrar la solución de la ecuación.
28 Tema 1 Ecuaciones diferenciales

EJEMPLO 1.20
dy 1
Para resolver la ecuación diferencial lineal + y = 3t utilizamos el primer método,
dt t
encontrando la solución de la ecuación homogénea
∫ −1 ) 1
φ(t) = e− = e− ln t = eln(t
1
dt
t = .
t
Calculamos ∫ ∫
q(t) 3t
c(t) = dt = dt = t3 ,
φ(t) 1/t
1
c(t)φ(t) = t3 . = t2 .
t
La solución buscada será
1
y(t) = t2 + c .
t
Utilizando el segundo método encontramos el factor integrante,
∫ 1
µ(t) = e t dt = t ,

multiplicando la ecuación diferencial por esta función, obtenemos

1
y ′ t + y = 3t2 ⇒ (y t)′ = 3t2 ⇒ yt = t3 + c ⇒ y(t) = t2 + c .
t

1.7. E.D.O. lineales de segundo orden

Las ecuaciones diferenciales ordinarias podemos clasificarlas en dos grandes bloques:
las lineales y las no lineales. Las más sencillas de estudiar son las del primer tipo
ya que debido a las propiedades de sus soluciones pueden caracterizarse de manera
general y además disponemos de métodos para resolver muchas de ellas.

DEFINICIÓN 1.7.1 Una ecuación diferencial lineal de orden n es una ecuación

del tipo
an (t)y n) + an−1 (t)y n−1) + · · · + a1 (t)y ′ + a0 (t)y = b(t) (1.21)
donde ai (t) i = 0, 1, 2, · · · , n y b(t) son funciones continuas en algún intervalo I y
además an (t) ̸= 0 , ∀t ∈ I.

DEFINICIÓN 1.7.2 La ecuación diferencial

an (t)y n) + an−1 (t)y n−1) + · · · + a1 (t)y ′ + a0 (t)y = 0 (1.22)

se llama ecuación diferencial lineal homogénea asociada a la ecuación (1.21).

1.7 E.D.O. lineales de segundo orden 29

Si las funciones ai (t) , i = 0, 1, 2, · · · , n son funciones constantes, entonces la ecuación

(1.21) se llama ecuación diferencial lineal de orden n con coeficientes constantes.
Nos centraremos en las ecuaciones diferenciales lineales de segundo orden por un
doble motivo. En primer lugar, podemos hacer un desarrollo teórico relativamente
simple y, en segundo lugar, son muy importantes desde el punto de vista práctico.
Una ecuación diferencial lineal de segundo orden es una ecuación del tipo

a2 (t)y ′′ + a1 (t)y ′ + a0 (t)y = b(t) (1.23)

donde a2 (t), a1 (t), a0 (t) y b(t) son funciones continuas en algún intervalo I y además
a2 (t) ̸= 0 , ∀t ∈ I.
Lo usual es escribir la ecuación (1.23) en su forma canónica

y ′′ + p(t)y ′ + q(t)y = g(t) . (1.24)

Empezaremos su estudio analizando la ecuación diferencial lineal homogénea aso-

ciada a (1.24), dada por
y ′′ + p(t)y ′ + q(t)y = 0 . (1.25)

TEOREMA 1.7.3 Sean p(t) y q(t) dos funciones continuas en algún intervalo I.
Entonces, para cualquier t ∈ I, el problema de valores iniciales

y ′′ + p(t)y ′ + q(t)y = 0 , y(t0 ) = y0 , y ′ (t0 ) = y0′

tiene una única solución definida en el intervalo I, cualesquiera que sean los valores
y0 , y0′ ∈ IR.

RESULTADO 1.7.4 Si y1 (t) e y2 (t) son dos soluciones de (1.25), entonces cual-
quier combinación lineal de ambas,

y(t) = c1 y1 (t) + c2 y2 (t) , c1 , c2 ∈ IR

es también solución de (1.25).

EJEMPLO 1.21

La ecuación diferencial lineal de segundo orden y ′′ + 4y = 0 tiene por soluciones

y1 (t) = cos 2t, y2 (t) = sen 2t. Por tanto, si hacemos uso del Resultado 1.7.4 la
función
y(t) = c1 sen 2t + c2 cos 2t , c1 , c2 ∈ IR
será también solución de y ′′ + 4y = 0.
30 Tema 1 Ecuaciones diferenciales

Si consideramos la ecuación diferencial de segundo orden no lineal

ty ′′ + 2yy ′ = 0 , (1.26)

es inmediato comprobar que las funciones

t
y1 (t) = 1 , y2 (t) =
1+t
son soluciones de (1.26). Sin embargo la función
c2 t
y(t) = c1 y1 (t) + c2 y2 (t) = c1 + , c1 , c2 ∈ IR
1+t
no es solución de (1.26). En efecto,
2c2 (c1 + c2 − 1) + 2c1 c2
ty ′ + 2yy ′ = ,
(t + 1)3
no es idénticamente nula para cualquier valor de c1 y c2 .

RESULTADO 1.7.5 Sean y1 (t) e y2 (t) dos soluciones de la ecuación diferencial

y ′′ + p(t)y ′ + q(t)y = 0 , (1.27)

definidas en el intervalo I, tales que

 
 y1 (t0 ) y2 (t0 ) 
 ′ ̸ 0
= (1.28)
 y1 (t0 ) y2′ (t0 ) 

para algún t0 ∈ I. Entonces cualquier solución de (1.27) es combinación lineal de

y1 (t) e y2 (t). Es decir,

y(t) = c1 y1 (t) + c2 y2 (t) , c1 , c2 ∈ IR

es la solución general de (1.27).

DEFINICIÓN 1.7.6 Dadas dos funciones y1 , y2 ∈ C 1 (I), se define el Wrons-

kiano de y1 e y2 como la función
 
 y1 (t) y2 (t) 
W [y1 , y2 ](t) =  ′  , t∈I. (1.29)
y1 (t) y2′ (t) 

Observemos que el resultado anterior lo que hace es reducir el problema de resolver

la ecuación diferencial (1.27) a encontrar dos soluciones particulares y1 (t) e y2 (t) que
cumplan con la condición
W [y1 , y2 ](t0 ) ̸= 0 ,
para algún t0 ∈ I. Dos funciones y1 (t) e y2 (t) con estas caracterı́sticas se dicen que
forman un conjunto fundamental de soluciones de la ecuación diferencial
(1.27).
1.7 E.D.O. lineales de segundo orden 31

EJEMPLO 1.22

Las funciones y1 (t) = e−2t e y2 (t) = e−4t son soluciones de la ecuación diferencial

y ′′ + 6y ′ + 8y = 0 (1.30)

en el intervalo −∞ < t < ∞. Además

   
 y1 (t) y2 (t)   e−2t e−4t 

W [y1 , y2 ](t) =  ′ =  = −2e−6t ̸= 0 , ∀t ∈ (−∞, ∞) .
y1 (t) y2′ (t)   −2e−2t − 4e −4t 

Por lo tanto, forman un conjunto fundamental de soluciones de (1.30) en el intervalo

(−∞, ∞). La solución general de (1.30) será

y(t) = c1 e−2t + c2 e−4t , c1 , c2 ∈ IR .

1.7.1. Método de reducción del orden

Si conocemos una solución particular de la ecuación lineal homogénea

y ′′ + p(t)y ′ + q(t)y = 0 (1.31)

podemos encontrar otra solución de (1.31) aplicando el método de reducción del

orden.
Sea y1 (t) una solución particular de (1.31) hacemos el cambio de variable y =
z(t)y1 (t) y derivamos

y ′ = z ′ y1 + zy1′ , y ′′ = z ′′ y1 + 2z ′ y1′ + zy1′′ .

Si sustituimos estos valores en (1.31) y tenemos en cuenta que y1 es una solución

particular de (1.31), la ecuación diferencial inicial se transforma en esta otra

y1 (t)z ′′ + (2y1′ (t) + p(t)y1 (t))z ′ = 0 .

Ahora el cambio v = z ′ reduce la ecuación anterior a la ecuación lineal homogénea

de primer orden
y1 (t)v ′ + (2y1′ (t) + p(t)y1 (t))v = 0 ,
que podemos resolver por separación de variables
( ∫ )
c
v(t) = 2 exp − p(t)dt , c ∈ IR ,
y1 (t)
y como sólo necesitamos una solución podemos tomar c = 1. Entonces
∫
′
z = v ⇒ z = v(t)dt ,
32 Tema 1 Ecuaciones diferenciales

y la nueva solución de (1.31) será

∫
y2 (t) = zy1 (t) = y1 (t) v(t)dt .

Puede probarse que estas dos soluciones forman un conjunto fundamental de solu-
ciones de la ecuación (1.31). La solución general de (1.31) podemos escribirla como

y(t) = c1 y1 (t) + c2 y2 (t) .

EJEMPLO 1.23

La función y1 (t) = e2t es una solución particular de la ecuación diferencial lineal

homogénea
y ′′ − 4y ′ + 4y = 0 .

Podemos encontrar una segunda solución utilizando el método de reducción del grado
de la ecuación diferencial.
Sea y(t) = z(t)y1 (t) = z(t)e2t , derivando

y ′ (t) = z ′ (t)e2t + 2z(t)e2t , y ′′ (t) = z ′′ (t)e2t + 4z ′ (t)e2t + 4z(t)e2t ,

sustituimos estas derivadas en la ecuación diferencial lineal homogénea inicial y

simplificamos
z ′′ (t) = 0 .

A continuación procedemos a rebajar el orden, para ello llamamos v(t) := z ′ (t) y

resolvemos la ecuación diferencial que aparece
∫
′
v (t) = 0 ⇒ v(t) = c = 1 ⇒ z(t) = v(t)dt = t ,

y la segunda de las soluciones buscada será y2 (t) = u(t)y1 (t) = te2t .

Estas dos soluciones forman un conjunto fundamental de soluciones de la ecuación
diferencial inicial. En efecto
   
 y1 (t) y2 (t)   e2t te 2t 
W [y1 , y2 ](t) =  ′ =
 
 = e4t ̸= 0 , ∀t ∈ (−∞, ∞) .
′
y1 (t) y2 (t) 2e 2t e + 2te 
2t 2t

La solución general de la ecuación diferencial inicial es

y(t) = c1 e2t + c2 te2t , c1 , c2 ∈ IR .

1.7 E.D.O. lineales de segundo orden 33

1.7.2. EDO lineal de segundo orden completa

La solución general de la ecuación diferencial lineal

y ′′ + p(t)y ′ + q(t)y = g(t) , (1.32)

la obtendremos a partir de las soluciones de su ecuación lineal homogénea asociada

y ′′ + p(t)y ′ + q(t)y = g(t) ,

y una solución particular de (1.32).

RESULTADO 1.7.7 Sea yp (t) una solución particular de la ecuación diferencial
lineal
y ′′ + p(t)y ′ + q(t)y = g(t) (1.33)
e {y1 (t), y2 (t)} un conjunto fundamental de soluciones de su ecuación diferencial
lineal homogénea asociada

y ′′ + p(t)y ′ + q(t)y = 0 . (1.34)

Entonces
y(t) = yp (t) + c1 y1 (t) + c2 y2 (t) , c1 , c2 ∈ IR .
será la solución general de (1.33)

1.7.3. Método de variación de parámetros

Como hemos visto en el Resultado 1.7.7, para poder encontrar la solución general
de (1.33) necesitamos conocer una solución particular. El método de variación de
parámetros nos proporciona un procedimiento para calcular dicha solución particu-
lar.
Supongamos que {y1 (t), y2 (t)} es un conjunto fundamental de soluciones de la ecua-
ción diferencial lineal homogénea asociada (1.33), entonces su solución general yh (t)
viene dada por
yh (t) = c1 y1 (t) + c2 y2 (t) , c1 , c2 ∈ IR .
El objetivo es encontrar una solución particular de (1.33) que sea de la forma

yp (t) = c1 (t)y1 (t) + c2 (t)y2 (t) , (1.35)

donde c1 (t) y c2 (t) son dos funciones a determinar. La duda que surge de forma
natural es saber si es posible encontrar dos funciones c1 (t) y c2 (t) tales que

yp (t) = c1 (t)y1 (t) + c2 (t)y2 (t),

sean una solución particular de la ecuación diferencial (1.33). Observemos que lo

que hemos realizado ha sido en la solución yh (t) = c1 y1 (t) + c2 y2 (t), reemplazar las
constantes por los parámetros variables c1 (t), c2 (t).
34 Tema 1 Ecuaciones diferenciales

Derivando
yp′ = c1 y1′ + c′1 y1 + c2 y2′ + c′2 y2 .
Si además exigimos que c1 (t) y c2 (t) sean funciones tales que

c′1 y1 + c′2 y2 = 0 ⇒ yp′ = c1 y1′ + c2 y2′ . (1.36)

Volviendo a derivar
yp′′ = c1 y1′′ + c′1 y1′ + c2 y2′′ + c′2 y2′ ,
y sustituyendo estos valores en (1.33)

yp′′ + p(t)yp′ + q(t)yp = c1 (y1′′ + p(t)y1′ + q(t)y1 )

+ c2 (y2′′ + p(t)y2′ + q(t)y2 )

+ c′1 y1′ + c′2 y2′ = g(t) .

Pero al ser y1 (t) e y2 (t) soluciones de la ecuación lineal homogénea

y1′′ + p(t)y1′ + q(t)y1 = 0 , y2′′ + p(t)y2′ + q(t)y2 = 0 .

Es decir
y1′ c′1 + y2′ c′2 = g(t) . (1.37)
De (1.36) y (1.37) obtenemos el sistema lineal de ecuaciones
{
y1 c′1 + y2 c′2 = 0
y1′ c′1 + y2′ c′2 = g(t) ,

que resolviéndolo, encontramos las soluciones:

W1 W2
c′1 (t) = , c′2 (t) = ,
W W
donde    
 0 y2   y 0 

W1 =   = −y2 g(t) , W2 =  1′  = y1 g(t)
g(t) y2′  y1 g(t) 

y W es el Wronskiano de y1 (t), y2 (t), que como sabemos viene dado por

 
 y1 y2 
W [y1 , y2 ](t) =  ′ .
y1 y2′ 

Resumiendo, para resolver la ecuación diferencial (1.33) procedemos de la manera

siguiente:

(a) Encontramos la función yh (t) = c1 (t)y1 (t) + c2 (t)y2 (t) y posteriormente evalua-
mos el Wronskiano W [y1 , y2 ](t).
1.7 E.D.O. lineales de segundo orden 35

(b) Obtenemos c1 (t), c2 (t) integrando las expresiones

−y2 (t) g(t) y1 (t) g(t)

c′1 (t) = , c′2 (t) = .
W [y1 , y2 ](t) W [y1 , y2 ](t)

(c) Construimos la solución particular

yp (t) = c1 (t)y1 (t) + c2 (t)y2 (t)

EJEMPLO 1.24

Supongamos que queremos resolver la ecuación diferencial lineal completa de segun-

do orden
2t 2
y ′′ − 2
y′ + y = 1 + t2 . (1.38)
1+t 1 + t2
(a) En primer lugar necesitamos encontrar un conjunto fundamental de soluciones
de la ecuación lineal homogénea asociada
2t 2
y ′′ − 2
y′ + y = 0. (1.39)
1+t 1 + t2
Es inmediato comprobar que una solución particular de (1.39) es y1 (t) = t.
Para calcular otra solución particular aplicamos el método de reducción del
grado. Para ello, si realizamos el cambio de variable

y(t) = z(t)y1 (t) = t z(t) ,

se llega a la ecuación diferencial

2
tz ′′ + z′ = 0 .
1 + t2
Llamando v(t) = z ′ (t) la ecuación diferencial anterior se transforma en

2 1 + t2
tv ′ + v=0 ⇒ v(t) = .
1 + t2 t2
Por tanto, ∫ ∫
1 + t2 1
z(t) = v(t)dt = 2
dt = t − .
t t
En consecuencia, la otra solución particular es
( )
1
y2 (t) = t t − = t2 − 1 .
t

La solución general de (1.39) viene dada por

yh (t) = c1 t + c2 (t2 − 1) , c1 , c2 ∈ IR .
36 Tema 1 Ecuaciones diferenciales

(b) Ahora buscamos una solución particular de la forma

yh (t) = c1 (t) t + c2 (t)(t2 − 1) ,

siendo    
 0 y2 (t)   0 t2 − 1 
  
 g(t) y2′ (t)   1+t 2 2t 
′
c1 (t) =  =   = 1 − t2 ,
 y1 (t) y2 (t)   t t2 − 1 
  
 y ′ (t) y2′ (t)   1 2t 
1

integramos
∫
t3
c1 (t) = (1 − t2 )dt = t − .
3
Del mismo modo
   
 y1 (t) 0   t 0 
   
 y ′ (t) g(t)   1 1 + t2 
c′2 (t) =  1 =  = t,
 y1 (t) y2 (t)   t t2 − 1 
   
 y ′ (t) y ′ (t)   1 2t 
1 2

y, por tanto,
∫
t2
c2 (t) = (t)dt = .
2

(c) Por consiguiente

( )
t3 t2 1
yp (t) = t− t + (t2 − 1) = (t4 + 3t2 ) .
3 2 6

(d) Finalmente, la solución general de (1.38) vendrá dada por

1
y(t) = (t4 + 3t2 ) + c1 t + c2 (t2 − 1) , c1 , c2 ∈ IR .
6

1.8. E.D.O. lineales de segundo orden con coefi-

cientes constantes
En esta sección estudiaremos ecuaciones diferenciales lineales de segundo orden del
tipo
y ′′ + a1 y ′ + a2 y = g(t) , (1.40)
donde a1 y a2 son constantes.
Para poder resolver estas ecuaciones procedemos tal y como lo hicimos en la sección
anterior.
1.8 E.D.O. lineales de segundo orden con coeficientes constantes 37

1.8.1. La ecuación lineal homogénea con coeficientes cons-

tantes
Consideremos la ecuación diferencial

y ′′ + a1 y ′ + a2 y = 0 , (1.41)

donde a1 y a2 son constantes.

Sabemos que la ecuación diferencial lineal de primer orden y ′ + ay = 0, siendo a una
constante, tiene por solución

y(t) = ce−at , −∞ < t < ∞ .

En consecuencia, es lógico tratar de determinar si existen soluciones exponenciales en

−∞ < t < ∞, de la ecuación lineal homogénea (1.41). Comprobaremos que todas
sus soluciones son funciones exponenciales o se construyen a partir de funciones
exponenciales.
Consideremos la ecuación diferencial

y ′′ + a1 y ′ + a2 y = 0 . (1.42)

Probamos una solución de la forma y(t) = eλt . Para ello derivamos y sustituimos en
(1.42)
eλt (λ2 + a1 λ + a2 ) = 0 .
Como eλt ̸= 0, ∀t ∈ (−∞, ∞), debe ocurrir que

λ2 + a1 λ + a2 = 0 .

Esta ecuación se conoce con el nombre de ecuación caracterı́stica de la ecuación

diferencial (1.42). Examinemos los diferentes casos que pueden presentarse:

Primer caso. La ecuación caracterı́stica tiene dos raı́ces, λ1 , λ2 , reales y dis-

tintas. Las soluciones

y1 (t) = eλ1 t , y2 (t) = eλ2 t ,

son linealmente independientes en −∞ < t < ∞ y por lo tanto forman un con-

junto fundamental de soluciones de la ecuación diferencial (1.42). La solución
general es
y(t) = c1 eλ1 t + c2 eλ2 t , c1 , c2 ∈ IR

Segundo caso. Cuando λ1 = λ2 , entonces sólo existe una solución exponencial

y(t) = eλt . Podemos encontrar una segunda solución utilizando el método de
reducción del grado de la ecuación diferencial.
38 Tema 1 Ecuaciones diferenciales

Sea y(t) = z(t)eλt , derivando

y ′ (t) = z ′ eλt + zλeλt , y ′′ (t) = z ′′ eλt + 2λz ′ eλt + zλ2 eλt ,

sustituimos estas derivadas en la ecuación diferencial (1.42), y simplificamos

z(λ2 eλt + a1 λeλt + a2 eλt ) + (z ′′ eλt + 2λz ′ eλt + a1 z ′ eλt ) = 0 ,

pero al ser eλt una solución de la ecuación diferencial, podemos simplificar la

expresión anterior y nos queda z ′′ + 2λz ′ + a1 z ′ = 0.
A continuación procedemos a rebajar el orden, para ello llamamos v(t) := z ′ (t)
y resolvemos la ecuación de variables separadas que aparece

v′
v ′ + (2λ + a1 )v = 0 ⇒ = −(2λ + a1 ) .
v
Es decir
∫
∫ − (2λ + a1 )dt
ln |v| = − (2λ + a1 )dt ⇒ v = k1 e = z′ .

Calculando el valor de z(t)

∫
∫ − (2λ + a1 )dt
z(t) = k1 e dt + k2 .

Si k1 = 1, k2 = 0
∫
∫ − (2λ + a1 )dt ∫
z(t) = e dt = e−(2λ+a1 )t dt .

Por otro lado, para que la ecuación λ2 + a1 λ + a2 = 0 tenga una raı́z doble,
tiene que ocurrir que su discriminante se anule. Calculando el valor de la raı́z
a1
λ=− .
2
Es decir, ∫ ∫
a1
−(2(− )+a1 )t
z(t) = e 2 dt = dt = t ,

y la segunda de las soluciones buscada será y2 (t) = z(t)y1 (t) = teλt .

La solución general de la ecuación diferencial inicial es

y(t) = c1 eλt + c2 teλt , c1 , c2 ∈ IR .

1.8 E.D.O. lineales de segundo orden con coeficientes constantes 39

Tercer caso. Si λ1 y λ2 son raı́ces complejas

λ1 = α + iβ , λ2 = α − iβt , α, β ∈ IR+ .

Estamos dentro del primer caso y por tanto

y(t) = k1 eα+iβ)t + k2 e(α−iβ)t .

A continuación, aplicamos la fórmula de Moivre para los números complejos y

simplificamos

y(t) = k1 eαt (cosβt + isenβt) + k2 eαt (cos(−βt) + isen(−βt))

= eαt ((k1 + k2 )cosβt + (k1 i − k2 i)senβt)

= eαt (c1 cosβt + c2 senβt)

EJEMPLO 1.25

La ecuación diferencial lineal homogénea

y ′′ + 6y ′ + 8y = 0 (1.43)

tiene por ecuación caracterı́stica

λ2 + 6λ + 8 = (λ + 2)(λ + 4) = 0 .

Las raı́ces son λ1 = −2 y λ2 = −4. En consecuencia, la solución general de la

ecuación diferencial (1.43) es de la forma

y(t) = c1 e−2t + c2 e−4t , c1 , c2 ∈ IR .

La ecuación diferencial lineal homogénea

y ′′ − 4y ′ + 4y = 0 (1.44)

tiene por ecuación caracterı́stica

λ2 − 4λ + 4 = 0 ,

que admite la solución real doble λ = 2. Por tanto, la solución general de (1.44) es

y(t) = c1 e2t + c2 te2t , c1 , c2 ∈ IR .

40 Tema 1 Ecuaciones diferenciales

La ecuación diferencial lineal homogénea

y ′′ − 4y ′ + 5y = 0 (1.45)

tiene por ecuación caracterı́stica

λ2 − 4λ + 5 = 0 ,

que admite las soluciones complejas conjugadas

λ1 = 2 + i , λ2 = 2 − i .

En consecuencia, la solución general de (1.45) es

y(t) = c1 e2t cos t + c2 e2t sen t , c1 , c2 ∈ IR .

1.8.2. La ecuación diferencial lineal completa de segundo

orden con coeficientes constantes
Como sabemos por la sección 9.7.2, una vez resuelta la ecuación lineal homogénea
asociada, la resolución de la ecuación diferencial completa
y ′′ + a1 y ′ + a2 y = g(t) (1.46)
se reduce a buscar una solución particular de (1.46), y podemos utilizar el método
de variación de parámetros para encontrarla.

EJEMPLO 1.26 texto

Supongamos que queremos encontrar la solución general de la ecuación diferencial

y ′′ − 2y ′ + y = (t − 1)et . (1.47)

(a) El polinomio caracterı́stico λ2 − 2λ + 1 tiene por raı́ces λ1 = 1, λ2 = 1. Por

tanto, yc (t) = c1 et + c2 tet .
Si y1 (t) = et , y2 (t) = tet , su Wronskiano vale

W [et , tet ] = W = e2t ̸= 0, ∀t ∈ (−∞, ∞)

.
(b) Calculamos
−y2 g(t) −tet (t − 1)et t3 t2
c′1 = = = −t2 + t ⇒ c1 = − +
W e2t 3 2
y
y1 f (t) et (t − 1)et t2
c′2 = = =t−1 ⇒ c2 = − −t
W e2t 2
1.8 E.D.O. lineales de segundo orden con coeficientes constantes 41

(c) Por consiguiente,

( 3 ) ( 2 ) ( 3 )
t t2 t t t2
yp (t) = − + t
e + − t te =
t
− et
3 2 2 6 2
(d) La solución general de (1.47) vendrá dada por
(3 2
)
y(t) = t6 − t2 et + c1 et + c2 tet , c1 , c2 ∈ IR .

1.8.3. Método de los coeficientes indeterminados

Ahora presentaremos un nuevo método para encontrar una solución particular de
la ecuación diferencial lineal completa con coeficientes constantes que no requiere el
cálculo de integrales.

EJEMPLO 1.27
Supongamos que queremos encontrar una solución particular de
y ′′ + 4y = e3t .
El método consiste en conjeturar la solución a la vista de la función g(t) = e3t .
Como en este caso estamos ante una función exponencial probamos con la solución
y = Ae3t . Si sustituimos en la ecuación diferencial llegamos a
9Ae3t + 4Ae3t = Ae3t ⇒ A = 1/13 ,
y la solución particular buscada es y(t) = 1/13 e3t .

Repitamos el método para encontrar una solución particular de la ecuación diferen-

cial
y ′′ + 2y ′ + y = t2 + e−t .
Ahora el segundo miembro g(t) está compuesto por dos tipos de funciones. La pri-
mera de ellas t2 sugiere ensayar con un polinomio de segundo grado At2 + Bt + C.
La segunda es la función exponencial e−t la cual nos indica que debemos buscar una
función del tipo Dt2 e−t , ya que tanto e−t como te−t son soluciones de la ecuación
homogénea. Por tanto, probamos con la función
y(t) = At2 + Bt + C + Dt2 e−t .
Al sustituir en la ecuación diferencial e identificar coeficientes se obtiene un sistema
de ecuaciones lineales que una vez resuelto presenta las soluciones
1
A = 1 , B = −4 , C = 6 , D = .
2
La solución particular buscada es
1
y(t) = t2 − 4t + 6 + t2 e−t .
2
42 Tema 1 Ecuaciones diferenciales

El procedimiento descrito en el ejemplo anterior se denomina método de los coeficien-

tes indeterminados. Se aplica cuando la función g(t) es de algunos tipos particulares.
Como regla general, probamos con una solución particular del mismo tipo que la
función g(t) y con coeficientes indeterminados, multiplicando por t o t2 , si fuese
necesario, para conseguir que ninguno de los términos de la solución ensayada sea
solución de la ecuación lineal homogénea asociada.

1.9. Notas históricas

Si no se tienen ciertos conocimientos de ecuaciones diferenciales y de los métodos
usados para resolverlas, es difı́cil estudiar la historia y el desarrollo de esta impor-
tante rama de las matemáticas. Más aún, la evolución de la teorı́a de las ecuaciones
diferenciales está ı́ntimamente ligada al desarrollo general de las matemáticas, y no
puede separarse de ella.
La teorı́a de las ecuaciones diferenciales se origina en los inicios del cálculo, con Isacc
Newton (1642-1727) y Gottfried Wilhelm Leibnitz (1646-1716) en el siglo XVII. Aún
cuando Newton realizó, relativamente, poco trabajo en la teorı́a de las ecuaciones
diferenciales, su desarrollo del cálculo y la aclaración de los principios básicos de
la mecánica proporcionaron una base para el desarrollo de sus aplicaciones, en el
siglo XVIII, con mayor alcance por parte de Euler. Newton clasificó las ecuaciones
de primer orden de acuerdo con las formas
dy/dx = f (x) ; dy/dx = f (y) ; dy/dx = f (x, y)
Para la última desarrolló un método de solución, usando series infinitas, cuando
f (x, y) es un polinomio en x e y. Era muy sensible a la crı́tica y, como consecuencia
de ello, tardó bastante en publicar muchos de sus descubrimientos.
Leibnitz llegó a los resultados fundamentales del cálculo independientemente, aunque
un poco más tarde que Newton. Nuestra notación moderna para la derivada dy/dx
y el signo de la integral se deben a Leibnitz. Descubrió el método de separación de
las variables, ası́ como procedimientos para resolver las ecuaciones homogéneas de
primer orden y las ecuaciones lineales de primer orden. Mantuvo una abundante
correspondencia con otros matemáticos, especialmente con los hermanos Bernoulli.
En el curso de esta correspondencia se resolvieron muchos problemas de ecuaciones
diferenciales, durante las últimas décadas del siglo XVII.
A Newton y Leibnitz le siguieron los hermanos Jakob Bernoulli (1654-1705) y Johann
Bernoulli (1667-1748) y, el hijo de Johann, Daniel Bernoulli (1700-1782). Justamen-
te, éstos son tres de los ocho miembros de la familia Bernoulli, quienes en su tiempo,
fueron prominentes matemáticos y hombres de ciencia. Con ayuda del cálculo, formu-
laron y resolvieron las ecuaciones diferenciales de muchos problemas de mecánica.
Un problema (1696-1697) al cual contribuyeron ambos hermanos, y el cual pro-
vocó problemas entre ellos, fue el de la braquistócrona3 que conduce a la ecuación
3
Determinación de la curva de descenso más rápido
1.9 Notas históricas 43

no lineal de primer orden

y(1 + (y ′ )2 ) = c
donde c es una constante. Newton también resolvió el problema antes, en 1697. Se
dice, tal vez no sea cierto, que Newton supo del problema al final de la tarde de
un fatigoso dı́a en la Casa de la Moneda, y lo resolvió en la noche, después de la
cena. Publicó la solución en forma anónima, pero J. Bernoulli, al verla, exclamó “...
conozco al león por su zarpa ...”
En 1690 J. Bernoulli publicó la solución de la ecuación diferencial, que en forma
diferencial se escribe
(b2 y 2 − a3 )1/2 dy = a3/2 dx
Actualmente esta ecuación se toma como un simple ejercicio, pero, en aquel tiempo,
encontrar la solución, constituyó un avance trascendental.
A finales del siglo XVII, muchos de los métodos elementales de solución para ecuacio-
nes diferenciales ordinarias de primer orden se conocı́an y, la atención se dirigió hacia
las ecuaciones diferenciales ordinarias de orden superior y hacia las ecuaciones dife-
renciales parciales. Jacob Riccati (1676-1754), matemático italiano, consideró ecua-
ciones de la forma f (y, y ′ , y ′′ ) = 0. Tambien consideró una importante ecuación no
lineal, conocida como ecuación de Ricatti

dy/dx = a0 (x) + a1 (x)y + a2 (x)y 2

aunque no en forma general.

Leonhard Euler uno de los matemáticos más grandes de todos los tiempos, tam-
bién vivió en el siglo XVII. Sus trabajos reunidos llenan más de sesenta volúmenes.
Aunque quedó ciego, durante los últimos diecisiete a+os de su vida, su trabajo no
disminuyó. De particular interés es su trabajo sobre el planteamiento de problemas
de la mecánica y su desarrollo de métodos de solución para estos problemas ma-
temáticos. Refiriéndose al trabajo de Euler en la mecánica, Lagrange dijo que era el
primer gran trabajo en el que se aplica el análisis a la ciencia del movimiento. Euler
también consideró la posibilidad de reducir ecuaciones de segundo orden a ecuacio-
nes de primer orden, mediante un cambio adecuado de variables; creó el concepto de
factor integrante, en 1739 dio un tratamiento general de las ecuaciones diferenciales
lineales ordinarias con coeficientes constantes. Contribuyó al método de las solucio-
nes en series de potencias y dio un procedimiento numérico para resolver ecuaciones
diferenciales. También hizo contribuciones importantes a la teorı́a de las series de
Fourier y creó la primera discusión sistemática del cálculo de variaciones.
En el siglo XVIII, los grandes matemáticos franceses Joseph- Louis Lagrange (1736-
1813) y Pierre-Simon Laplace (1749- 1827) hicieron importantes aportaciones a la
teorı́a de las ecuaciones diferenciales. Posiblemente sea la ecuación de Laplace, la
ecuación diferencial en derivadas parciales más conocida en la fı́sica matemática, la
ecuación del potencial
uxx + uyy = 0
44 Tema 1 Ecuaciones diferenciales

donde los subı́ndices indican derivadas parciales. El trabajo monumental de La-

grange, Mecanique analytique, contiene las ecuaciones generales del movimiento de
un sistema dinámico, conocidas actualmente como las ecuaciones de Lagrange. Pa-
ra Laplace la naturaleza era esencial y las matemáticas, eran su herramienta en el
aprendizaje de sus secretos; para Lagrange las matemáticas eran un arte que justifi-
caba su propio ser. Sin embargo, ambos hombres realizaron avances de gran alcance,
tanto en la teorı́a como en las aplicaciones de las matemáticas.
En los últimos años, algunos matemáticos dedicados al estudio de las ecuaciones dife-
renciales ordinarias y parciales han tratado de elaborar una teorı́a sistemática (pero
general) rigurosa. La finalidad no es tanto crear métodos de solución para ecuaciones
diferenciales particulares, sino desarrollar técnicas apropiadas para el tratamiento
de diferentes clases de ecuaciones .