Tema 4

Ecuaciones diferenciales ordinarias.

Muchos problemas ba´sicos de las ciencias experimentales pueden ser modelados usando
ecuaciones donde aparecen involucradas una funcio´n junto con sus derivadas. En este cap´ıtulo nos
centraremos en el estudio de este tipo de ecuaciones y veremos co´mo pueden ser resueltas utilizando
lo ya estudiado en los temas anteriores.
Comenzaremos definiendo lo que entenderemos por ecuacio´n diferencial ordinaria y dando
varios ejemplos cla´sicos de su uso en las ciencias experimentales. Posteriormente estudiaremos co´mo
resolver algunos tipos de e´stas ecuaciones.

4.1 Ecuaciones diferenciales ordinarias.

Dada una funcio´n y = f (x) vamos a estudiar ecuaciones donde aparecen mezcladas la variable x, la
funcio´n y y algunas de sus derivadas y0 (x), y00 (x)... Ejemplos de este tipo de ecuaciones son
x + y0 (x)
xy0 (x) = y(x), y00 (x) = y(x)y0 (x), y000 (x)2 = .
1 + x2

Definicio´n 46 (Ecuacio´n diferencial ordinaria)

Llamaremos ecuacio´n diferencial ordinaria (abreviado EDO) a una ecuacio´n que involucra a
una variable independiente x, una funcio´n y(x) y una o varias derivadas de y(x).
Llamaremos orden de una EDO al orden de la mayor derivada en la ecuacio´n. As´ı, por ejemplo, las
EDOs
2x + yy0 2 d y
3 ✓ ◆4
dy d2y
= 0, — = cos x cos y 0
dx3 + = ,
dx dx2
1,
x
tienen o´rdenes respectivos 1, 3 y 2.
Adema´s, una EDO de orden n diremos que es lineal si se puede escribir de la forma
dny
a dn—1y dy
n(x) + a n—1(x) + .. . + a1(x) + a0(x)y = b(x),
dxn dxn—1 dx

donde las ai, b son funciones que solo dependen de x.

Ejemplos de EDOs lineales son las siguientes

senx y00 y ex dy
1 d3y x4 1 d4y d2y x4y x
0 x x2
+ = , +( — ) = +
1+ 91
dx3 , dx dx4
x2
 + = .
dx2
92 4.1. Ecuaciones diferenciales ordinarias.

Obse´rvese, sin embargo, que las siguientes EDOs no son lineales

✓ ◆2 dy
dy
yy0 = + y = 0, cos x + cos y = x.
dx dx
1,
Por otro lado, necesitamos saber que´ entenderemos como solucio´n de una EDO.
Definicio´n 47 (Solucio´n de una EDO)
Dada una EDO de orden n, llamaremos solucio´n de esta EDO a una funcio´n y = f (x) definida
en un intervalo I y que tiene al menos n derivadas continuas en el intervalo I, de forma que
cuando sustituimos y = f (x) y sus derivadas en la EDO se cumple la ecuacio´n en todo el
intervalo.
Consideremos la EDO y0 + xy = x, entonces es fa´cil comprobar que la funcio´n y = 1 2
2 es solu-
+e
2
—x
cio´n en todo R. De hecho, para cualquier constante c2 R la funcio´n y = 1 + 2 es una solucio´n de la
ce
EDO en todo R. Es decir, esta EDO tiene infinitas soluciones en R.
Para la EDO y00 + y = 0 las funciones y = cos x, y = sen x son soluciones de la ecuacio´n en todo
R. Y, en general, dadas constantes reales c1 , c2 cualesquiera, la funcio´n y = c1 cos x + c2 sen x es
solucio´n en todo R de la EDO.
Observemos que la funcio´n y = 1/x es solucio´n de la EDO y0 + y2 = 0, pero no esta´ definida en
todo R, por lo que podremos decir que es una solucio´n en el intervalo
— ( •, 0) y que es solucio´n en
el intervalo (0, •).
Como hemos visto, el conjunto de soluciones de una EDO puede ser infinito. A veces nos interesa
conocer solamente alguna solucio´n de la EDO y no todas, por lo que se prefijan algunas
condiciones previas sobre la solucio´n o soluciones buscadas. Un caso especialmente interesante es el
siguiente
Definicio´n 48 (Problema de valores iniciales)
Dada una EDO de orden n, un intervalo I y un punto x0 2 I, llamaremos problema de valores
iniciales (abreviado PVI) al problema de encontrar las funciones que cumplen la EDO de orden n
anterior en el intervalo I y que adema´s satisfacen las condiciones iniciales
y(x0) = c0, y0(x0) = c1, . . . , y(n—1)(x) = cn—1,
donde los ci son nu´meros prefijados.

El PVI en R dado por

(
2yy0 = ex
y(0)= 1
x/2
tiene como solucio´n a la funcio´n y = e . De hecho puede comprobarse que es la u´nica solucio´n al
PVI. De manera ana´loga el PVI en R que tiene la misma EDO anterior pero diferente condicio´n
( inicial
2yy0 = ex
y(0)= —1
tiene como u´nica solucio´n a la funcio´n y = —ex/2 .
La EDO de segundo orden y00 — y = x tiene por soluciones a las funciones y = —x + c1ex + c2e—x para
cualesquiera constantes c1, c2 2 R. As´ı, el PVI en R dado por
( 00
y —y=x
y(1)= 0, y 0 (1)= 0
tiene como u´nica solucio´n y = —x + ex—1 .
Hay que tener en cuenta que aunque en los ejemplos anteriores los PVIs solo tienen una solucio
´n, podr´ıa ocurrir en otros casos que haya ma´s de una.
Tema 4 93

4.2 Modelos matema´ticos inspirados en problemas de las ciencias

experi- mentales.

En muchas ocasiones es deseable describir en te´rminos matema´ticos el comportamiento de

algunos sistemas o feno´menos de la vida real, tanto de tipo f´ısico, qu´ımico, sociolo´gico, econo
´mico... La des- cripcio´n matema´tica de un sistema de feno´menos se llama modelo matema´tico y se
construye con ciertos objetivos, como el de entender que´ ocurrira´ en el futuro o que´ paso´ en el
pasado. Por ejemplo, podemos desear entender los mecanismos de cierto ecosistema al estudiar el
crecimiento de la poblacio´n animal en ese sistema, o podemos desear datar fo´siles y analizar el
decaimiento de una sustancia radiactiva ya sea en el fo´sil o en el estrato en que e´ste fue descubierto.
Para la realizacio´n de un modelo matema´tico sobre un sistema, primero hay que identificar todas
las variables que ocasionan que el sistema cambie y posteriormente establecer de que´ manera estas
variables afectan al sistema.
Es claro que cuantas ma´s variables que afectan al sistema se an˜adan mejor resolucio´n
tendra´ el modelo que se obtenga. Sin embargo, el modelo sera´ cada vez ma´s complejo cuantas ma´s
variables entren en juego. Por ello, a veces un modelo de baja resolucio´n (es decir, con pocas
variables) es suficiente para determinar de forma aproximada la solucio´n de nuestro problema.
A continuacio´n estudiaremos algunos modelos matema´ticos cla´sicos en diferentes a´reas de las
Cien- cias Experimentales.

4.2.1 Dina´mica poblacional.

El economista Thomas Malthus hizo uno de los primeros intentos para modelar el crecimiento de
una poblacio´n. As´ı, supuso que la razo´n de la poblacio´n en un cierto tiempo t es proporcional a la
poblacio´n total en ese tiempo. Es decir, si llamamos P(t) a la poblacio´n en el tiempo t entonces se
obtiene que
dP
= k P,
dt
donde k es una constante que depende de la poblacio´n estudiada.
Aunque este modelo es de baja resolucio´n, au´n as´ı sigue siendo u´til para el estudio de algunas
pobla- ciones a tiempo corto, como, por ejemplo, poblaciones de cultivos de bacterias.

4.2.2 Decaimiento radiactivo.

El nu´cleo de algunos a´tomos esta´n formados por combinaciones de protones y neutrones
inestables, es decir, los a´tomos se desintegran o se convierten en a´tomos de otras sustancias. En estos
casos se dice que los nu´cleos son radiactivos. Por ejemplo, el radio Ra-226 intensamente radiactivo se
transforma en gas rado´n Rn-222.
Para modelar el feno´meno de decaimiento radiactivo, dada una cantidad A(t) de una sustancia
en el tiempo t, se supone que la razo´n con la que los nu´cleos se desintegran es proporcional a la
cantidad
existente, esto es,
dA
= k A,
dt
donde k es una constante que depende de la sustancia estudiada.

4.2.3 Ley de enfriamiento/calentamiento de Newton.

La ley emp´ırica de enfriamiento/calentamiento de Newton establece que la rapidez con la que cambia
la temperatura de un cuerpo es proporcional a la diferencia entre la temperatura del cuerpo y la del
medio que lo rodea (la temperatura ambiente).
94 4.2. Modelos matema´ticos inspirados en problemas de las ciencias experimentales.

Si denotamos por T (t) a la temperatura del cuerpo en el tiempo t y Ta a la temperatura del ambiente,
entonces, la ley de enfriamiento/calentamiento de Newton determina que
dT
= k (T (t) — Ta),
dt
donde k es la constante de proporcionalidad.

4.2.4 Reacciones qu´ımicas.

En algunas reacciones qu´ımicas la rapidez con la que las mole´culas de una sustancia A se
descom- ponen es proporcional a la cantidad de sustancia que queda sin reaccionar. As´ı, si llamamos
c(t) a la cantidad de la sustancia A en el tiempo t entonces dA/dt = k A, donde k es una constante
(dependiente de la sustancia A).
Un ejemplo de este tipo, llamada reaccio´n de primer orden, esta´ dada por la conversio´n del
cloruro de terbutilo (CH3)3CCl en alcohol t-but´ılico (CH3)3COH:

(CH3)3CCl + NaOH —! (CH3)3COH +

NaCl.

Solo la concentracio´n de cloruro de terbutilo controla la rapidez de la reaccio´n.

Sin embargo, en la reaccio´n

CH3Cl + NaOH —! CH3OH + NaCl,

la razo´n con la que avanza la reaccio´n es proporcional al producto de las concentraciones de cloruro
de metilo CH3Cl y de hidro´xido de sodio NaOH que quedan.
Para describir en general una reaccio´n de este tipo, conocida como reaccio´n de segundo orden,
su- pongamos que se combina una mole´cula de una sustancia A con una mole´cula de una sustancia B
para formar una mole´cula de una sustancia C (ma´s otras sustancias). Si c(t) denota la cantidad de la
sustancia C que se ha formado en el tiempo t y si a0, b0 son, respectivamente, las cantidades de las
sustancias A y B en el momento inicial t = 0, entonces las cantidades de A y B en tiempo—t son a0 c(t)
— y b0
c(t), respectivamente. As´ı, la razo´n de formacio´n de C esta´ dada por
dc
= k(a0 — c)(b0 — c),
dt
donde k es una constante de proporcionalidad.

4.2.5 Mezclas.
Supongamos que en un tanque mezclador tenemos 300 litros de agua con sal. Por otro lado, un
grifo de entrada introduce 3 litros de agua por minuto, que tiene una concentracio´n de 2 gramos por
litro. La solucio´n bien mezclada sale del tanque con la misma rapidez de 3 litros por minuto.
De esta manera, si denotamos por A(t) a la cantidad de sal en el tanque al tiempo t, entonces la razo´n
con la que A(t) cambia viene dada por
dA
= (razo´n de entrada de sal)—(razo´n de salida de sal).
dt
La razo´n de entrada de sal es de (2g/l) (3l/min) = 6g/min. Por otra parte, ya que la cantidad de
agua en el tanque es constantemente de 300 litros, la razo´n de salida de sal300es de ( A(t ) g/l)
(3l/min) =
A(t)
100 g/min.
De esta manera,
dA A
=6— .
dt 100
 Tem 9

4.2.6 Cuerpos en ca´ıda y resistencia del aire.

Aunque a un cuerpo en ca´ıda en el vac´ıo solo le afecta la fuerza de la gravedad, en general,
cuando no se encuentra en el vac´ıo, el aire ejerce una resistencia al movimiento en la ca´ıda de un cuerpo.
As´ı, el peso del cuerpo es una fuerza que actu´a en la direccio´n de ca´ıda, mientras que la
resistencia del aire actu´a en direccio´n opuesta. La fuerza Fr dada por la resistencia del aire se
denomina amortiguamiento viscoso y es proporcional a la velocidad del cuerpo, esto es, Fr = kv, donde
k es una constante que depende del cuerpo y v su velocidad.
As´ı, la segunda ley del movimiento de Newton determina que la suma de las fuerzas ejercidas sobre
el cuerpo es igual a su masa m por su aceleracio´n a = dtd v, es decir,
dv
m = mg —
kv, dt
donde g es la aceleracio´n de la gravedad.
O equivalentemente, si denotamos por s(t) al espacio recorrido en el tiempo t entonces
d2s ds
m = mg— k .
dt2 dt

4.3 EDOs de primer orden.

A lo largo de esta seccio´n estudiaremos algunas clases de ecuaciones diferenciales ordinarias de

primer orden y co´mo encontrar sus soluciones.
Debemos tener en cuenta que en algunas ocasiones la solucio´n y(x) del problema nos aparecera´ de
forma impl´ıcita, esto es, por ejemplo podemos obtener soluciones del tipo

y2 + 2xy + x2 = 0, y + cos y = x3 — 1.

En el primer caso podremos despejar la funcio´n y en te´rminos de x con lo que obtendremos su

forma expl´ıcita y = —x, pero en el segundo caso no se podra´ obtener la forma expl´ıcita de y.

4.3.1 EDOs separables.

Un caso especial de EDOs de primer orden son aquellas en las que la ecuacio´n diferencial
puede reescribirse de manera que un lado de la igualdad solo dependa de la funcio´n y y el otro lado
solo de la variable x. Ma´s concretamente,

Definicio´n 49 (EDO separable)

Diremos que una EDO de primer orden es separable o que tiene variables separables si se puede escribir de la fo

dy
g(y) dx = h(x).

Por ejemplo, la ecuacio´n y0 = ex—2y cos x es una EDO separable ya que se puede reescribir como e2y d y =
dx
ex cos x. Sin embargo, la EDO y0 = y2 — sen x no es separable.
Debe observarse que una ecuacio´n del tipo y0 = f (y)g(x) puede reescribirse como una EDO
separable
de la forma f 1 dxdy = g(x), pero en este proceso hay que tener en cuenta que f (y) no puede ser cero
porque no podemos
(y) dividir por cero. Adema´s, los nu´meros c para los cuales f (c) = 0 cumplen que
la funcio´n constante y(x) = c es una solucio´n de la EDO, ya que sustituyendo en la EDO original
cumplen la ecuacio´n.
9 4.3. EDOs de primer

1
As´ı,
dy por ejemplo, la EDO y0 = (1 — y2)x2 es separable ya que la podemos reescribir de la forma
2 = x2. Aqu´ı hay que tener en cuenta que las funciones y(x)= —1, y(x)= 1 son dos soluciones de
1—y dx
la EDO original; e´stos son los valores donde 1 — 2 = 0.
y

¿Co´mo se resuelven las EDOs separables?

Dada la EDO g(y) dy = h(x), basta escribirla como g(y)dy = h(x)dx e integrar a ambos lados de la
dx
igualdad, esto es,
Z g(y)dy = Z h(x)dx.

De esta manera se obtendra´n todas las soluciones de la EDO de manera impl´ıcita.

EJEMPLOS

Ejemplo 1 : Dada la EDO y2 y0 = x — 2, podemos encontrar sus soluciones reescribiendo la ecuacio

´n como y2dy = (x— 2)dx e integrando a ambos lados, con lo que obtenemos

y3 + c1 = x2 — 2x + c2.

De
p esta forma, las soluciones vienen dadas por y3 = x2 — 2x + c, o equivalentemente, y =
3 2
x — 2x + c.
y
Ejemplo 2 : La EDO y0 = 2 puede resolverse integrando sobre la igualdad 1 dy = 1
2 dx, salvo cuan-
1+x y 1+x
do y = 0. La funcio´n y(x) = 0 es solucio´n de la EDO y el resto de soluciones pueden
calcularse integrando la igualdad anterior que nos dara´

ln |y| = arctgx + c.

Tomando exponenciales sobre la igualdad obtenemos y = c0earctgx, donde c0 =±ec. Por

tanto, todas las soluciones vienen dadas por y = c0 earctgx , para c0 cualquier nu´mero
real, incluyendo el caso c0 = 0 que nos dara´ la solucio´n previa y(x) = 0.

Ejemplo 3 : Consideremos el PVI dado por la EDO ex+y y0 = x con condicio´n inicial y(0) = 0. La
ecuacio´n puede reescribirse como

ey dy = xe—x dx

y, as´ı, si integramos se obtiene ey = —(1 + x)e—x + c. Ahora, la condicio´n inicial

y(0) = 0 nos dice que 1 = —1 + c, de donde c = 2. Esto es, la solucio´n del PVI
viene dada por ey = —(1 + x)e—x + 2, o equivalentemente y = ln(2 — (1 + x)e—x).

Ejemplo 4 : El PVI con ecuacio´n y0 = —2xy2 y condicio´n inicial y(3) =10 1 puede ser resuelto reescri-
biendo la EDO como
dy
= —2xdx
y2
e integrando. Obse´rvese que estamos dividiendo por y2 y podr´ıa ocurrir que y(x) = 0,
que es solucio´n de la EDO pero no es la solucio´n buscada del PVI ya que no cumple la
condi- cio´n inicial. As´ı, tenemos que 1 = x2 + c, esto es, y = 1 2 . Ahora, de la
condicio´n inicial
y c+x
1
obtenemos que 10 = 3 + c, por lo que y =
2
1+x 2

.
⇤
Tem 9

4.3.2 E
DOs lineales de primer orden.
Otras EDOs de primer orden que pueden ser resueltas de una forma sencilla son aquellas de tipo
lineal, ma´s concretamente,

Definicio´n 50 (EDO lineal de primer orden)

Una EDO de primer orden es lineal si se escribe de la forma

a1(x)y0 + a0(x)y = b(x),

donde a1, a0, b son funciones que solo dependen de la variable x.

Por ejemplo, las ecuaciones y0 + x2y = ex cos x, o, exy0 + y = x3 son EDOs lineales de primer orden. Sin
embargo, las ecuaciones y0 + y2 = sen x, y0 + cos y = 0, o, yy0 + xy = ex no son EDOs lineales.
Una EDO lineal se dice que es homoge´nea cuando la funcio´n b(x) de la definicio´n anterior es
cero. As´ı, ejemplos de EDOs lineales de primer orden homoge´neas son sen(x)y0 + cos(x)y = 0, ex y0
+ y = 0, o, (x2 — 1)y0 + xy = 0.

¿Co´mo se resuelven las EDOs lineales de primer orden?

Para resolver la ecuacio´n a1 (x)y0 + a0 (x)y = b(x), primero hemos de resolver la ecuacio´n
ho- moge´nea asociada a1 (x)y0 + a0 (x)y = 0. Esta u´ltima EDO es separable y puede resolverse
como se explico´ en la seccio´n anterior. La solucio´n de la EDO homoge´nea es de la forma
y(x) = c0 y1 (x) para una cierta constante c0.
Si b(x) no es cero, las soluciones de la ecuacio´n a1 (x)y0 + a0 (x)y = b(x), vendra´n dadas de la
forma y(x) = u(x)y1 (x), para una cierta funcio´n u(x) a determinar. Es decir, se sustituira´ la
constante c0 dada en las soluciones de la EDO homoge´nea por la funcio´n u(x).
Entonces se sustituye la funcio´n y(x) = u(x)y1 (x) sobre la ecuacio´n a1 (x)y0 + a0 (x)y = b(x) y
esto nos dara´ una ecuacio´n de donde podremos despejar u0 (x). Finalmente, integrando la
expresio´n de u0(x) calculamos u(x) y, por tanto, las soluciones de la EDO.

EJEMPLOS

Ejemplo 1 : La EDO (1 + x2 )y0 + x y = 0 es lineal homoge´nea, por lo que se puede resolver

simplemente teniendo en cuenta que es de tipo separable. La ecuacio´n se puede escribir
como
dy x
=— ,
y 1 + x2

siempre que y no se anule. De hecho, y(x) = 0 es solucio´n del problema. Si se integra la

ecuacio´n anterior obtenemos ln |y| = — 1 ln(1 + x2 )+ c. O equivalentemente y = p c0 ,
para
2 2
1+x
c0 un nu´mero real cualquiera, que incluye la solucio´n y(x) = 0 tomando c0 =
0.

Ejemplo 2 : Para resolver la EDO de primer orden lineal arctg(x)y0 — 1 2 2+xy = arctg3(x) primero se
debe resolver la EDO homoge´nea arctg(x)y0 — +x 2 ´
1 2 y = 0. E sta la podemos reescribir
(para
y 6= 0) como 2dx
= .
dy
y arctg(x)(1 + x2)

Con lo que integrando obtenemos y = c0arctg2(x).

9 4.3. EDOs de primer

Por tanto, las soluciones de la EDO arctg(x)y0 — 1+x2 2 y = arctg3(x) las obtendremos de la
forma y = u(x)arctg2 (x) donde, para calcular u(x), sustituimos en la ecuacio´n, con lo que

2
arctg3(x)= arctg(x)(u(x)arctg2(x))0 — (u(x)arctg2(x)) = u0(x)arctg3(x).
1

Tenemos entonces que u 0 (x)= 1 o bien u(x)= x + c, por lo que las soluciones de la EDO
inicial son y = (x + c)arctg2(x).

Ejemplo 3 : Resolvamos ahora el PVI dado por la EDO x y0 —x2 x2

y = e 2 con condicio´n inicial y(1) = —1.
Primero resolvemos la EDO homoge´nea asociada x y0 — x2 y = 0 que puede ser escrita como
dy
= xdx
y
x2
para y 6= 0. Si integramos y simplificamos, las soluciones vienen dadas por y = c0e 2 .
x2 x2
As´ı, las soluciones de la EDO x y 0— x2 y = e 2 sera´n del tipo y = u(x)e 2 , para una cierta
funcio´n u(x) a determinar. Sustituyendo en la ecuacio´n se tiene:

x2 x2 x2 x2
e 2 = x y0 — x2 y = x(u(x)e2 )0 — x2(u(x)e2 ) = xu0(x)e 2 .

Por tanto, u 0 (x)= 1/x, es decir, u = ln |x| + c. De esta manera, todas las soluciones son
x2
y = (c + ln |x|)e 2 .
La condicio´n inicial y(1) = —1 nos da
⇣ —1 = (c + ln
⌘ 1)e2 2 . Por tanto, c = —
1 p1 , es decir, la
x e
funcio´n que resuelve el PVI es y = — p1 + ln x e .
2
e

2
Ejemplo 4 : Para resolver el PVI con ecuacio´n xy0 — y =1 x y condicio´n inicial y(—1) = 4p primero
resolvemos la EDO lineal homoge´nea xy0 — y 2= 0. Esta ecuacio´n la escribimos como
dy dx
=
y x

para y 6= 0. Si integramos y simplificamos, las soluciones vienen dadas por y = c0 x.

2
x
Buscamos entonces las soluciones de la ecuacio´n xy0 — y =1 de la forma y = u(x) x.
Sustituyendo obtenemos 2

x2 0 0 2

0
= xy — y = x(u(x) x) — u(x) x = u (x)x .
1 + x2
1
Por tanto, u 0 (x)= 1 +x2 , y as´ı u(x)= c + arctg x.
Entonces
x
2 las soluciones de la EDO lineal xy0 — y = son y = (c +arctg x)x. Finalmente,
la condicio´n inicial y c1 p
nos da p

11 c p . As´ı, c 0 y
1+x2
(— ) = 4 4 = ( + arctg(— ))(— ) = —( — 4 ) =
la solucio´n del PVI esta´ dada por la funcio´n y = x arctg x.

⇤
 Tem 9

4.3.3 E
DOs exactas.
Consideremos una funcio´n F(s,t) que depende de dos variables. Llamaremos derivada parcial de
F respecto de s a la derivada, si existe, de la funcio´n F respecto de la variable s cuando se considera
∂ F
a t como una constante, y la denotaremos ∂s . Ana´logamente se define la derivada parcial de F
respecto de
t, ∂∂t F . 2
´
Por ejemplo, dada la funcion F(s,t)= s + cos(s + t)+ arctgt sus derivadas parciales vienen dadas
por
∂F ∂F 1
= 2s— sen(s + t), = —sen(s + t)+ .
∂s ∂t 1 + t2
Definicio´n 51 (EDO exacta)
Diremos que una EDO de primer orden es exacta si se escribe de la forma

f (x, y) dy + g(x, y)= 0,

dx
y se cumple que ∂ f = ∂ g .
∂x∂y

2y
As´ı, por ejemplo, las ecuaciones x+y 2 y0 + x+y1 2 = 0, xexy y0 + yexy + 2x = 0 son exactas ya que
✓ ◆ ✓ ◆
∂ 2y — 2 y ∂ 1 ∂ (xexy) ∂ (yexy + 2x)
= = ,
xy
=
∂x (x + y2)2 x+ = (1 + xy) ∂y
2 ∂x
e .
x+ ∂y y
y2
Las ecuaciones 3x2 yy0 + sen(x + y )= 0, ex y2 y0 x3 + 2xy = 0 no son EDOs exactas ya que
—
no cumplen la igualdad necesaria de las derivadas parciales.
¿Co´mo se resuelven las EDOs exactas?
Si
dx
la EDO f (x, y) dy + g(x, y)= 0 (o equivalentemente f (x, y)dy + g(x, y)dx = 0) es exacta entonces
las soluciones vendra´n dadas de manera impl´ıcita de la forma F(x, y) = 0, donde la funcio´n de dos
variables F(x, y) se puede calcular de una cualquiera de las dos siguientes formas:
1. Se toma F(x, y)= F1(x)+ Z f (x, y)dy,

donde para calcular F1(x) se debe tener en cuenta que

∂x
∂ F = g(x, y).

2. Se toma F(x, y)= F2(y)+ Z g(x, y)dx,

donde para calcular F2(y) se debe tener en cuenta que

∂y
∂ F = f (x, y).

EJEMPLOS

Ejemplo 1 : Es fa´cil ver que la EDO (x — 2y)y0 + 2x + y = 0 es exacta ya que

∂ ( x — 2y) ∂ (2 x + y )
∂x =1= ∂y .

Las soluciones vendra´n dadas por F(x, y) = 0, donde tomamos por ejemplo
Z
F(x, y)= F1(x)+ (x— 2y)dy = F1(x)+ xy— y2.
1 4.3. EDOs de primer

Y la funcio´n F1 (x) se calcula si se tiene en cuenta que

∂F
2x + y = = F10 (x) + y.
∂x

Por tanto, F10 (x) = 2x o lo que es lo mismo F1 (x) = x2 + c.

Como conclusio´n, las soluciones de nuestra EDO vienen dadas de forma impl´ıcita como
x2 + x y — y2 + c = 0. Despejando la funcio´n y en la igualdad anterior se obtiene
⇣ p ⌘ ⇣ p ⌘
1 2
1 2
y= x + 4c + 5x , y= x— 4c + 5x
2 2
.
Ejemplo 2 : La EDO 2x2 yy0 + 2xy2 + 1
x = 0 es exacta ya que
✓ ◆
∂ (2x2 y) ∂ 2 1
.
= 4xy = 2xy +
∂x ∂y x

Sus soluciones son de la forma F(x, y)= 0, donde tomamos por ejemplo
Z ✓ ◆
1
F(x, y)= F2(y)+ 2xy2 + x dx = F2(y)+ x2y2 + ln |x|.

Aqu´ı la funcio´n F2 (y) se calcula tomando

∂F
2x2 y = = F20 (y) + 2x2 y.
∂y

Por tanto, F20 (y) = 0, con lo que F2 (y) = c.

De esta manera las soluciones de nuestra EDO vienen dadas de forma impl´ıcita como x2y2 +
ln |x| + c = 0.

Ejemplo 3 : Consideremos el PVI con ecuacio´n (1—x2 )y0 + 3x2—2x y 1 = 0 y condicio´n inicial y(0) =
1. Ya que la EDO es exacta porque —

∂ (1 x2)
— ∂ (3x2 — 2xy—
= —2x =
∂x 1)
∂y ,
podemos resolverla tomando
Z
F(x, y)= F1(x)+ (1 — x2)dy = F1(x)+ (1 — x2)y.

Y la funcio´n F1 (x) se calcula teniendo en cuenta que

∂F
3x2 — 2x y — 1 = = F10 (x) — 2x y.
∂x

Por tanto, F10 (x) = 3x2 — 1 o lo que es lo mismo F1 (x) = x3 — x + c. Y las soluciones de
la EDO son las dadas por x3 — x + c + (1 — x2)y = 0.
Usando la condicio´n inicial y(0) = 1 se tiene que c + 1 = 0, por lo que la solucio´n del
PVI es y = x—x+ 1 1 2 .
 Tem 1

Ejemplo 4 : El PVI, con ecuacio´n 3(1 + x2 )y2 y0 + 2x y3 —+x1 2x 4 = 0 y condicio´n inicial y(0) = 2,
tiene
una EDO exacta ya que ✓ ◆
∂ (3(1 + x2)y2) 2 ∂ 3 2 x
= 6xy =
∂x ∂y 2xy — .
1 + x4
As´ı, la ecuacio´n se puede resolver
tomando ✓ ◆
Z
2x
F(x, y)= F2(y)+ 2xy3 — dx = F2(y)+ x2y3 — arctg(x2),
1
donde F2(y) se calcula tomando
∂F
3(1 + x2 )y2 = = F20 (y) + 3x2 y2 .
∂y

Por tanto, F20 (y) = 3y2 , con lo que F2 (y) = y3 + c. Y las soluciones de la EDO esta´n
dadas por las funciones y3 + c + x2y3 — arctg(x2)= 0.
Si usamos la condicio´n inicial y(0) = 2 obtenemos 8 + c = 0, por lo que la solucio´n del
PVI q
3
es y = 8+arctg(x2)
1+x2
.

4.4 EDOs lineales de segundo orden.

En esta seccio´n trataremos las ecuaciones lineales de segundo orden que son de la forma

a2 y00 + a1 y0 + a0 y = b(x),

donde a0 , a1 , a2 son constantes reales, con a2 0, y b(x) es una funcio´n que depende de la variable x.
Comenzaremos con las ecuaciones lineales de segundo orden homoge´neas, esto es, con aquellas
en las que b(x)= 0.
¿Co´mo se resuelven las EDOs lineales de segundo orden homoge´neas?
Una EDO del tipo a2 y00 + a1 y0 + a0 y = 0 tiene por soluciones aquellas del tipo y(x)= c1 y1(x)+ c2 y2
Se considera la ecuacio´n de segundo grado a2 m2 + a1 m + a0 = 0 y tomamos las dos soluciones
m1, m2 de esta ecuacio´n. Ahora se distinguen casos segu´n los valores de estas soluciones:

1. Si m12, m son dos nu´meros reales diferentes (esto ocurre cuando21 a —24a0 a > 0) entonces

y1(x)= em1x,y2(x)= em2x.

2. Si m12 = m (esto ocurre cuando21a — 4a

2 0 a = 0) entonces

y1(x)= em1x, y2(x)= xem1x.

2
3. Si m1 = a + ib , 2m = a — ib son dos nu´meros complejos (esto ocurre cuando
1 a 2—0 4a a < 0)
con a, b reales, entonces
y1(x)= eax cos(bx),
y2(x)= eaxsen(bx).
1 4.4. EDOs lineales de segundo

EJEMPLOS

Ejemplo 1 : Para resolver la EDO y00 — 3y0 + 2y = 0 consideramos la ecuacio´n de segundo grado
m2 — 3m + 2 = 0.
Sus soluciones son
p
=3 ± (—3)2 — 4 · 2
m
2,
es decir, m1 = 1, m2 = 2. Ya que las soluciones son dos nu´meros reales diferentes,
entonces las soluciones de la EDO son
y = c1ex + c2e2x,
para cualesquiera constantes c1, c2.
Ejemplo 2 : La EDO 4y00 — 4y0 + y = 0 tiene por ecuacio´n de segundo grado asociada
4m2 — 4m + 1 = 0.

Las soluciones de esta ecuacio´n son

p
m= 4± (—4)2 — 4 · 4
2·
,
4 las dos ra´ıces del polinomio de segundo
es decir, las dos soluciones m1 = m2 = 12. Ya que
grado son iguales, entonces las soluciones de la EDO son
y = c 1x 1
e + c xe x,
2 2
1 2

para cualesquiera constantes c1, c2.

Ejemplo 3 : La EDO de segundo orden y00 — 4y0 + 13y = 0 tiene por ecuacio´n de segundo grado asociada
m2 — 4m + 13 = 0.
Sus soluciones son
p p
4± (—4)2 — 4 · 2 ± —9,
m= =
13 2
es decir, las dos soluciones son los nu´meros complejos conjugados m1 = 2 + 3i, m2 =
—2 3i.
Por tanto, las soluciones de la EDO son
y = c1 e2x cos(3x)+ c2 e2x sen(3x),
para cualesquiera constantes c1, c2.
Ejemplo 4 : Resolvamos el PVI ⇢
y00 4y0 + 5y = 0
— y(0)= 3, y 0 (0)= 1

Primero resolvemos la EDO que tiene asociado el polinomio m2 — 4m + 5 = 0. Las solu-

ciones de esta igualdad son los nu´meros complejos conjugados m1 = 2 + i, m2 = 2 — i.
Por tanto, las soluciones de la EDO son y = c1e2x cos x + c2e2xsen x.
Si ahora usamos las condiciones iniciales obtenemos 3 = y(0) = c1 y, adema´s, ya
que y0 = e2x((2c1 + c2) cos x— (c1 — 2c2)sen x) entonces 1 = y0(0)= 2c1 + c2, es decir, c2
= —5. As´ı, la solucio´n del PVI es y = e2x (3 cos x — 5sen x).
Tem 1

Cuando la EDO es del

tipo a2 y00 + a1 y0 + a0 y = b(x),
donde b(x) no es cero, su solucio´n se calcula de forma diferente. En esta caso se utiliza un me´todo
parecido al que se usaba para las EDOs lineales de primer orden.

¿Co´mo se resuelven las EDOs lineales de segundo orden generales?

Para resolver una EDO del tipo a2 y00 + a1 y0 + a0 y = b(x) primero se resuelve la ecuacio´n ho- moge
Y las soluciones de la EDO original sera´n de la forma y(x) = u1(x) y1(x) + u2(x) y2(x) donde las funciones

(
u1 (x) y1 (x)+u0(x) y2 (x)= 0
b(x)
u00
1 (x) y (x)+u022y
(x)(x)= a2

las derivadas u0i(x) e integrar.

EJEMPLOS

Ejemplo 1 : Resolvamos la EDO y00 + y0 — 2y = 3ex. Para ello, primero consideramos la EDO ho- moge
´nea asociada y00 + y0 — 2y = 0.
Su polinomio de segundo grado asociado es

m2 + m— 2 = 0,

que tiene por ra´ıces m1 = 1, m2 =—2. As´ı, las soluciones de la EDO homoge´nea son y
= c1ex + c2e—2x.
Para las soluciones de la EDO no homoge´nea tomamos y = u1 (x)ex + u2 (x)e—2x , donde
⇢ 0
u 1(x) ex + u02(x) e—2x = 0
u01 (x) ex — 2 u20 (x) e—2x = 3ex

Resolviendo el sistema anterior u01 (x) = 1, u20 (x) = —e3x , con lo que u1 (x) = x + c1 ,
u2 (x) =
—31 e3x + c2. Por tanto, las soluciones de la EDO son
1
y = (x + c1)ex + (— e3x + c2)e—2x = (x + d1)ex + c2 e—2x,
3
para ciertas constantes d1, c2.

Ejemplo 2 : Ahora resolvamos la EDO 9y00 + 6y0 + y = sen x. Su EDO homoge´nea asociada es 9y00
+ 6y0 + y = 0, que tiene por ecuacio´n de segundo grado

9m2 + 6m + 1 = 0.

Las soluciones son m1 = m2 = — 13 . Por tanto, las soluciones de la EDO homoge´nea son
1 1
y = c1 e— 3 x + c2 xe— 3 x.
1 4.4. EDOs lineales de segundo

Consideramos entonces las soluciones de la EDO inicial de la forma y = u1(x) — 1x

e (x) xe— 1 x, donde las funciones u (x) se calculan usando que 3 +
2 i
u 3

1 1
0
(x) e—3 x +1 u0 (x) xe—3 x = 0 sen x
: u1 0 (x) e— x 2— 1 (x— 3)u0 (x) e— x =
1

1 3

— 3 u1 3 3
1 1x
02 10 19
x
Si resolvemos el sistema anterior obtenemos que u1(x)= — 9 xe
3 sen x, u2(x)=
3
9 e sen x,
con lo que
1 1x
u (x)= c + e (3(5x— 3) cos x— (5x + 12)sen x) , u
(x)= c — 1 e 1 (3 cos x— sen x) .
3 x
3
1 1
150 2 2
30
Por tanto, las soluciones de la EDO son
1 1 1
y = c e— 3 x + c xe—3 x — (3 cos x + 4sen x).
1 2
50
Ejemplo 3 : Consideremos la EDO y00 + 2y0 + 2y = e—x . Para resolver la EDO homoge´nea asociada
y00 + 2y0 + 2y = 0 consideramos la ecuacio´n

m2 + 2m + 2 = 0,

que tiene por soluciones los nu´meros complejos conjugados m1 = —1 + i, m2 = —1 — i.

As´ı, las soluciones de la EDO homoge´nea son

y = c1 e—x cos x + c2 e—x sen x.

Las soluciones de la EDO no homoge´nea son de la forma y = u1 (x) e—x cos x +u2 (x) e—x sen x,
donde ⇢ 0
u1 (x) e—x cos x + u02 (x) e—x sen x = 0
—u01 (x) e—x (cos x + sen x) + u20 (x) e—x (cos x — sen x) = e—x
De esta forma, u01 (x) = —sen x, u20 (x) = cos x, con lo que u1 (x) = c1 + cos x, u2 (x) = c2 +
sen x.
Las soluciones de la EDO son

y = e—x(1 + c1 cos x + c2sen x).

Ejemplo 4 : Para resolver el PVI

⇢
y00 — 6y0 + 9y = 2e3x
y(0)= 1, y 0 (0)= —2

Primero resolvemos la EDO homoge´nea y00 — 6y0 + 9y = 0 que tiene asociado el polinomio
m2 — 6m + 9 = 0. Las soluciones de esta igualdad son m1 = m2 = 3.
Por tanto, las soluciones de la EDO homoge´nea son y = c1 e3x + c2 xe3x .
Para resolver y00 — 6y0 + 9y = 2e3x tomamos soluciones de la forma y = u1(x)e3x +
u2(x)xe3x donde (
u01(x) e3x + u02(x) xe3x = 0
3 u01(x) e3x + (3x + 1)u02(x) e3x = 2e3x

As´ı, u01 (x) = —2x, u20 (x) = 2, con lo que, u1 (x) = —x2 + c1 , u2 (x) = 2x + c2 .
Tem 1

Con todo esto, las soluciones de la EDO y00 — 6y0 + 9y = 2e3x son

y = e3x(x2 + c2x + c1).

Finalmente, usando las condiciones iniciales y(0)= 1, y 0 (0)= —2 obtenemos que la solu-
cio´n al PVI viene dada por la funcio´n

y = e3x(x2 — 5x + 1).