C APÍTULO 1

Métodos elementales de
resolución de ecuaciones
diferenciales ordinarias

1. La población P(t) de un suburbio de una gran ciudad en un instante

cualquiera se rige por
−1 − 10−7 P )
 dP
dt = P ( 10 ,
P(0) = 5000
en donde t se mide en meses. ¿Cuál es el valor límite de la población?
¿En qué momento será la población igual a la mitad de su valor lími-
te?

Solución : Calculamos en primer lugar el tamaño de la población,

P(t), resolviendo el problema de valores iniciales. La ecuación di-
ferencial tiene sus variables separadas:
P0
= 1,
P(10−1 − 10−7 P)
donde hemos denotado P0 = dP dt . Integrando los dos miembros de
esta identidad entre 0 y t obtenemos
Z P(t)
7 dQ
10 = t,
5000 Q(106 − Q)
donde hemos efectuado el cambio de variable Q = P(t). Teniendo
en cuenta ahora que
 
1 −6 1 1
= 10 + ,
Q(106 − Q) Q 106 − Q

1
2

concluimos tras una serie de cálculos simples que la única solución

de nuestro problema es
t
106 e 10
P(t) = t .
199 + e 10

El valor límite de la población es por tanto

lı́m P(t) = 106 ,

t→∞

como se desprende de una simple aplicación de la regla de L’Hôpital.

Para responder a la segunda cuestión tenemos que encontrar el valor
6
t0 para el que P(t0 ) = 102 . Basta entonces con resolver la ecuación
t0
6 e 10 106 t0
10 t0 = ⇔ e 10 = 199 .
199 + e 10 2

Tomando logaritmos neperianos en ambos miembros concluimos que

t0 = 10 log(199) meses ≈ 4,41 años .

2. Resuelve las siguientes ecuaciones diferenciales:

(a) x0 = et − 2t
t2 −1

(b) ( x2 + 9) y0 + xy = 0
dy
(c) dx = 2xe− y

(d) x0 = 1+t
t2 x2

(e) x0 = et+x

Solución : (a) La ecuación tiene sus variables separadas. Integrando

obtenemos

x(t) = et − log(|t2 − 1|) + C , C ∈ R.

2
Métodos elementales 3

6
1·10

800000

600000

400000

200000

20 40 60 80 100 120 140

Figura 1.1: Representación gráfica de la solución del Ejercicio 1 en el intervalo

[0, 150].

3
4

(b) Separando las variables obtenemos

y0 x
=− 2
y x +9

e integrando con respecto a x llegamos a

C
y( x) = √ .
x2 + 9

dy
(c) Separando las variables resulta e y dx = 2x, de donde se obtiene la
solución general

y( x) = log( x2 + C ) , C ∈ R : x2 + C > 0 ,

sin más que integrar ambos miembros con respecto a la variable x.

Obsérvese que, dado cualquier dato inicial y( x0 ) = y0 , la solución
sólo existe si
x2 > −C = x20 − e y0 .

(d) Separando las variables obtenemos

1+t
x2 x0 = .
t2
Integrando entonces con respecto a t en ambos miembros de la ecua-
ción encontramos que la solución general de la misma viene dada
por
   1
1 3
x(t) = 3 log(|t|) − +C , C ∈ R.
t

(e) Separando las variables resulta e−x x0 = et , de donde obtenemos

la solución general

x(t) = − log(C − et ) , C > et ,

integrando la ecuación con respecto a la variable t. Obsérvese que,

dado cualquier dato inicial x(t0 ) = x0 , la solución sólo existe si

t < log(C ) con C = et0 + e−x0 .

4
Métodos elementales 5

3. Un reactor transforma plutonio 239 en uranio 238 que es relativa-

mente estable para uso industrial. Después de 15 años se determina
que el 0.0043 por ciento de la cantidad inicial A0 de plutonio se ha
desintegrado. Determina la semivida1 de este isótopo si la rapidez
de desintegración es proporcional a la cantidad restante.

Solución : Llamemos x(t) a la cantidad de plutonio 239 que queda

en el instante t, con lo que x0 (t) indicará la velocidad o rapidez de
desintegración del mismo. Como la velocidad de desintegración es
proporcional a la cantidad de isótopo restante, la ley diferencial que
rige el proceso de desintegración es

x0 = λx sujeta a la condición inicial x(0) = A0 ,

cuya única solución viene dada por x(t) = A0 eλt . Para tener com-
pletamente determinada la solución necesitamos conocer el valor de
la constante de desintegración λ, el cual puede encontrarse a través
de la relación (establecida en el enunciado del problema)
 0,0043  99,9957
x(15) = 1 − A0 = A0 ,
100 100
por lo que ha de ser
99,9957 1  99,9957 
A0 e15λ = A0 ⇔ λ = log .
100 15 100

Finalmente, la semivida de este isótopo es el valor t0 para el que se

cumple la condición x(t0 ) = A20 , por lo que

A0
 99,9957


1
log t0
A0 e 15 100 =
2
15 log(2)
⇔ t0 = = 241790 años .
log(100) − log(99,9957)

4. Dadas dos funciones f , g derivables, sabemos que la identidad

( f g)0 = f 0 g0
1 Tiempo necesario para que la cantidad inicial de átomos se reduzca a la mitad

5
6

3 +2x
es falsa en general. Si fijamos f ( x) = e x , determina las funciones
g que verifican dicha identidad.

Solución : Por un lado

3 +2x 3 +2x
( f g)0 ( x) = (3x2 + 2) e x g( x) + e x g0 ( x)

mientras que, por otro lado,

3 +2x
f 0 ( x) g0 ( x) = (3x2 + 2) e x g0 ( x) .

Entonces ha de cumplirse
3 +2x 3 +2x
(3x2 + 1) e x g0 ( x) = (3x2 + 2) e x g( x)

o, equivalentemente,

g0 ( x) 3x2 + 2
= 2 .
g( x) 3x + 1
Resolviendo esta ecuación diferencial en variables separadas obtene-
mos √
x+ √1 arctan( 3x)
g( x) = C e 3 , C ∈ R.

5. Resuelve las siguientes ecuaciones diferenciales:

(a) 3x + y − 2 + y0 ( x − 1) = 0

(b) (t2 x2 − 1) x0 + 2tx3 = 0, haciendo x = zα

(c) x + ( x − t) x0 = 0

(d) 2t + 3x + ( x + 2) x0 = 0

Solución : (a) La ecuación puede ser reescrita en forma canónica (con

la derivada despejada) de la siguiente forma:
3x + y − 2
y0 = .
1−x

6
Métodos elementales 7

Veamos cómo podemos reducir esta ecuación diferencial a una ho-

mogénea. Consideramos las rectas de ecuaciones

3x + y − 2 = 0 , 1−x = 0,

de las que resulta el punto de corte ( x = 1, y = −1). A continuación

se procede vía el siguiente cambio de variable:

X = x−1, Y = y+1.

Entonces se tiene que

Y + 3X Y
Y 0 = y0 = − = +3, (1.1)
X X
que es una ecuación diferencial homogénea. Haciendo ahora el cam-
Y
bio de función incógnita u = X y usando (1.1) obtenemos

Y 0 = u + Xu0 = u + 3 ,

de donde se deduce que

3
u0 = .
X
Por tanto
u( X ) = 3 log(| X |) + C , C ∈ R.
Deshaciendo los cambios de variable efectuados para recuperar las
variables originales llegamos a

y( x) = 3( x − 1) log(| x − 1|) + C ( x − 1) − 1 , C ∈ R.

(b) Obviamente x ≡ 0 es solución. Busquemos todas las demás so-

luciones. Efectuando el cambio de función incógnita x = zα en la
ecuación obtenemos

α (t2 z2α − 1) zα −1 z0 + 2tz3α = 0 ,

de donde resulta

α (t2 z2α − 1) z0 + 2tz2α +1 = 0

o, equivalentemente,

tz2α +1 ( z/t)2α +1
   
0 2 2
z =− =−
α t2 z2α − 1 α ( z/t)2α − (1/t2α +2 )

7
8

que es una ecuación homogénea. Haciendo el cambio de variable

u = zt obtenemos la siguiente ecuación equivalente:
 2α +2 2α
u −1

t 1
u0 = .
u t
Integrando los dos miembros de esta ecuación con respecto a la va-
riable t llegamos a la siguiente expresión implícita para u:
1 2α +2 2α
t u − log(|u|) = C + log(|t|) , C ∈ R.
2α
Deshaciendo finalmente los dos cambios de variable efectuados2 ob-
tenemos
x2 t2 − 2 log(| x|) = C , C ∈ R .

(c) La ecuación admite la solución trivial. Busquemos las restantes

soluciones. Resolviendo la ecuación con respecto a la derivada obte-
nemos x
x
x0 = = t x,
t−x 1− t
x
que es homogénea. El consabido cambio de variable u = t nos con-
duce a  2 
0 1 u
u = ,
t 1−u
cuyas soluciones satisfacen la relación implícita
1
− − log(|u|) = log(|t|) + C , C ∈ R.
u
Deshaciendo el cambio de variable obtenemos finalmente

x log(| x|) + t + Cx = 0 , C ∈ R.

(d) Resolviendo la ecuación con respecto a la derivada obtenemos

3x + 2t
x0 = − ,
x+2
que adopta la forma de una ecuación diferencial reducible a homo-
génea. Consideramos las rectas de ecuaciones

3x + 2t = 0 , x+2 = 0,
2u z
= t y x = zα

8
Métodos elementales 9

de las que resulta el punto de corte (t = 3, x = −2). Efectuamos el

siguiente cambio de variable:

T = t−3, X = x+2.

Entonces se tiene que

3X + 2T 3X + 2
X 0 = x0 = − = − TX , (1.2)
X T

X
que es homogénea. Haciendo ahora u = T y usando (1.2) obtenemos

3u + 2
X 0 = u + Tu0 = − ,
u
de donde se deduce que

u2 + 3u + 2
 
0 1
u =− .
T u

Por tanto, la siguiente relación implícita

 2
u( T ) + 2
= C|T | , C > 0,
u( T ) + 1

es satisfecha. Deshaciendo los cambios de variable efectuados para

recuperar las variables originales llegamos primero a

( X + 2T )2
= CT 2
X+T
y después a
( x + 2t − 4)2
= C (t − 3)2 .
x+t−1


6. Encuentra la curva plana que pasa por (1, 1) y verifica que dado un
punto cualquiera ( x, y), al considerar el corte de la recta tangente a
la curva en ese punto con el eje de ordenadas y el corte de la recta
normal con el eje de abscisas, su distancia al origen es la misma.

9
10

Solución : Las ecuaciones de la recta tangente y normal en el punto

(t, x) son, respectivamente,
 1
( T − t) x0 = X − x , ( T − t) − 0 = X − x ,
x
donde las variables están representadas con letras mayúsculas (T pa-
ra las ordenadas y X para las abscisas).
La condición inicial que nos proporciona el problema es x(1) = 1.
Sustituyendo en las ecuaciones anteriores los datos del problema ob-
tenemos:  1
−tx0 = b − x , ( a − t) − 0 = − x .
x
Además, ha de verificarse
| x − tx0 | = | xx0 + t| .
Esta última ecuación nos conduce a la resolución de los siguientes
problemas de valores iniciales:

x = xx−
 0 t
+t ,
x(1) = 1
 0
x = tx−+xt
,
x(1) = 1
correspondientes ambos a ecuaciones homogéneas. El primero de
ellos puede reescribirse de la siguiente forma
x
−1
(
x0 = xt +1
t .
x(1) = 1
Haciendo el cambio de variable u = xt se llega a la siguiente ecuación
diferencial con variables separadas

1 u2 + 1
 
0
u =− ,
t u+1
con dato inicial asociado u(1) = x(1) = 1, cuya única solución3
satisface la siguiente relación implícita:
x π
2 arctan + log( x2 + t2 ) = + log(2) .
t 2
3 Es inmediato verificar las condiciones del teorema de existencia y unicidad de solu-
ciones

10
Métodos elementales 11

El tratamiento del segundo problema es análogo. La ecuación dife-

rencial que se obtiene ahora tras hacer el cambio de variable anterior
es
1 u2 + 1
 
0
u = ,
t 1−u
con dato inicial asociado u(1) = x(1) = 1. La única solución de este
problema de valores iniciales satisface
x π
2 arctan − log( x2 + t2 ) = − log(2) .
t 2

7. Encuentra las soluciones de las siguientes ecuaciones diferenciales

buscando (si es el caso) factores integrantes de la forma que se indica:

(a) sen(tx) + tx cos(tx) + t2 cos(tx) x0 = 0

sen2 ( x)
 
sen(2x)
(b) y + x + y − y2
y0 = 0

(c) (3xy2 − 4y) + (3x − 4x2 y) y0 = 0 con µ ( x, y) = xn ym

(Febrero 1993)

(d) xy0 ( y − 1) − y = 0 con µ ( x, y) = µ ( y)

(e) (t + 1)2 + (1 + t2 ) x0 = 0 con µ (t, x) = µ (t + x)

(Febrero 1996)

(f) (1 + xy + y2 ) + (1 + xy + x2 ) y0 = 0 con µ ( x, y) = µ ( xy)

(g) ( x + y2 ) + 2( y2 + y + x − 1) y0 = 0 con µ ( x, y) = µ (e ax+by )

(h) 2xyy0 = x2 + y2 + 1 con µ ( x, y) = µ ( y2 − x2 )

Solución : (a) Es exacta con

P(t, x) = sen(tx) + tx cos(tx) , Q(t, x) = t2 cos(tx) .

11
12

15

10

-4 -2 2 4

-5

-10

-15

Figura 1.2: Representación gráfica de las soluciones del Ejercicio 6.

12
Métodos elementales 13

En efecto,
∂P ∂Q
= 2t cos(tx) − t2 x sen(tx) = .
∂x ∂t
Entonces
∂F
= P ⇒ F (t, x) = t sen(tx) + H ( x) .
∂t
Por otro lado
∂F
= Q ⇒ H0 ≡ 0 ⇒ H ≡ k ∈ R .
∂x
Por tanto, la solución general responde a la siguiente relación implí-
cita:
t sen(tx) = C ∈ R .

(b) Es exacta con

sen(2x) sen2 ( x)
P( x, y) = +x, Q( x, y) = y − .
y y2
En efecto,
∂P 2 sen( x) cos( x) ∂Q
=− = .
∂y y 2 ∂x
Entonces
∂F cos(2x) x2
= P ⇒ F ( x, y) = − + + H ( y) .
∂x 2y 2
Por otro lado
∂F 1 y2 1
= Q ⇒ H 0 ( y) = y − 2 ⇒ H ( y) = + .
∂y 2y 2 2y
Por tanto, la solución general responde a la siguiente relación implí-
cita:
1 
y2 + x2 + 1 − cos(2x) = C ∈ R .
y
(c) La ecuación no es exacta, ya que para las funciones

P( x, y) = y(3xy − 4) , Q( x, y) = x(3 − 4xy)

no se satisface la condición de exactitud. En efecto,

∂P ∂Q
= 2(3xy − 2) 6= 3 − 8xy = .
∂y ∂x

13
14

Buscamos n, m ∈ Z tales que la función µ ( x, y) = xn ym sea un factor

integrante. Para ello imponemos exactitud sobre las funciones

P? ( x, y) = y(3xy − 4) xn ym , Q? ( x, y) = x(3 − 4xy) xn ym .

Entonces resulta

2(3xy − 2) xn ym + y(3xy − 4) xn (mym−1 )

= (3 − 8xy) xn ym + x(3 − 4xy) ym (nxn−1 )

o, equivalentemente,

xn ym [(3m + 4n) xy − (3n + 4m) − (7 − 14xy)] = 0 . (1.3)

Por tanto, para que µ ( x, y) = xn ym sea un factor integrante para

nuestra ecuación se ha de cumplir

3m + 4n = −14 , 3n + 4m = −7 ,

es decir, m = 2 y n = −5. Luego µ = x−5 y2 . Buscamos ahora F ( x, y)

tal que, en primer lugar,

∂F
= y3 x−5 (3xy − 4) ⇒ F ( x, y) = − x−3 y4 + x−4 y3 + H ( y) .
∂x
Por otro lado
∂F
= x−4 y2 (3 − 4xy) = −4x−3 y3 + 3x−4 y2 + H 0 ( y) ⇒ H ≡ k ∈ R .
∂y
Consecuentemente, la solución responde al siguiente esquema im-
plícito:
x−4 y3 (1 − xy) = C ∈ R .

(d) Sean P( x, y) = − y y Q( x, y) = x( y − 1). En este caso la ecuación

no es exacta, ya que

∂P ∂Q
≡ −1 , = y−1.
∂y ∂x

Buscamos un factor integrante de la forma µ ( x, y) = µ ( y). Para ello

multiplicamos la ecuación por µ ( y) e imponemos la condición de
exactitud a las funciones

P? ( x, y) = − yµ ( y) , Q? ( x, y) = x( y − 1)µ ( y) .

14
Métodos elementales 15

Entonces ha de cumplirse
− yµ 0 ( y) − µ ( y) = ( y − 1)µ ( y) ⇒ µ ( y) = e− y .
Buscamos finalmente F ( x, y) tal que
∂F ∂F
= − ye− y , = x( y − 1)e− y .
∂x ∂y
Se tiene que
F ( x, y) = − xye− y + H ( y) y x( y − 1)e− y = − xe− y + xye− y + H 0 ( y) ,
de lo cual resulta finalmente que H ≡ k ∈ R, por lo que la solución
de la ecuación de partida responde a la siguiente relación implícita:
xye− y = C ∈ R .

(e) La ecuación no es exacta, ya que si

P(t, x) = (t + 1)2 , Q(t, x) = 1 + t2 ,

se tiene que
∂P ∂Q
= 0 6= 2t = .
∂x ∂t
Buscamos un factor integrante de la forma µ (t, x) = µ (t + x). La
condición suficiente y necesaria que ha de cumplirse para que en
nuestro caso µ (t + x) sea factor integrante es

(t + 1)2µ 0 (t + x) = 2tµ (t + x) + (1 + t2 )µ 0 (t + x) .
Resolviendo esta ecuación obtenemos µ (t + x) = et+x . Buscamos en-
tonces F (t, x) tal que
∂F
= (t + 1)2 et+x ⇒ F (t, x) = (1 + t2 )et+x + H ( x) .
∂t
Finalmente, la función H ( x) queda determinada sin más que impo-
ner la condición
∂F
= (1 + t2 )et+x + H 0 ( x) = Q(t, x)µ (t + x) = (1 + t2 )et+x ,
∂x
que nos conduce a H ≡ k ∈ R. Por tanto, la solución x(t) viene dada
por la ley implícita

(1 + t2 ) et+ x = C ∈ R .

15
16

(f) La ecuación no es exacta. En efecto, si llamamos

P( x, y) = 1 + xy + y2 , Q( x, y) = 1 + xy + x2 ,

se tiene que
∂P ∂Q
= x + 2y 6= y + 2x = .
∂y ∂x
Buscamos un factor integrante de la forma µ ( x, y) = µ ( xy). La con-
dición suficiente y necesaria que ha de cumplirse para que en nuestro
caso µ ( xy) sea factor integrante es

( x + 2y)µ ( xy) + x(1 + xy + y2 )µ 0 ( xy)

= ( y + 2x)µ ( xy) + y(1 + xy + x2 )µ 0 ( xy) .

Resolviendo esta ecuación obtenemos µ ( xy) = e xy . Buscamos enton-

ces F ( x, y) tal que

∂F
= (1 + xy + y2 )e xy ⇒ F ( x, y) = ( x + y)e xy + H ( y) .
∂x
Finalmente, la función H ( y) queda determinada al imponer la con-
dición
∂F
= ( x2 + xy + 1)e xy + H 0 ( y) = Q( x, y)µ ( xy) = (1 + xy + x2 )e xy ,
∂y

que nos conduce a H ≡ k ∈ R. Por tanto, la solución y( x) viene dada

por la ley implícita
( x + y) e xy = C ∈ R .

(g) La ecuación no es exacta, ya que si llamamos

P( x, y) = x + y2 , Q( x, y) = 2( y2 + y + x − 1) ,

se tiene que
∂P ∂Q
= 2y 6= 2 = .
∂y ∂x
Buscamos un factor integrante de la forma µ = µ (e ax+by ), el cual
habrá de verificar la siguiente condición suficiente y necesaria:

2yµ (e ax+by ) + b( x + y2 )e ax+byµ 0 (e ax+by )

= 2µ (e ax+by ) + 2a( y2 + y + x − 1)e ax+byµ 0 (e ax+by ) ,

16
Métodos elementales 17

que nos conduce como única opción a elegir un factor integrante de

la forma específica µ ( x, y) = e x+2y . Buscamos ahora F ( x, y) tal que
∂F
= ( x + y2 )e x+2y ⇒ F ( x, y) = ( x − 1 + y2 )e x+2y + H ( y) .
∂x
Finalmente, la función H ( y) queda determinada al imponer la con-
dición
∂F
= 2( y + x − 1 + y2 )e x+2y + H 0 ( y)
∂y
= Q( x, y)µ ( x, y) = 2( y2 + y + x − 1)e x+2y ,
que nos conduce nuevamente a H ≡ k ∈ R. Por tanto, la solución
y( x) viene dada por la ley implícita
( x − 1 + y2 ) e x+2y = C ∈ R .

(h) La ecuación no es exacta, ya que si llamamos

P( x, y) = x2 + y2 + 1 , Q( x, y) = −2xy ,
se tiene que
∂P ∂Q
= 2y 6= −2y = .
∂y ∂x
Buscamos un factor integrante de la forma µ = µ ( y2 − x2 ), el cual
habrá de verificar la siguiente condición suficiente y necesaria:

2yµ ( y2 − x2 ) + 2y( x2 + y2 + 1)µ 0 ( y2 − x2 )

= −2yµ ( y2 − x2 ) + 4x2 yµ 0 ( y2 − x2 ) .
Resolviendo esta ecuación en variables separadas obtenemos
1
µ ( y2 − x2 ) = .
( y2 − x2 + 1)2
Buscamos ahora F ( x, y) tal que
∂F x2 + y2 + 1 x
= 2 ⇒ 4 F ( x, y) = 2 + H ( y) .
∂x 2
( y − x + 1) 2 y − x2 + 1
4 Para x2 + y2 +1
calcular una primitiva de ( y2 −x2 +1)2 con respecto a x basta con observar que la
siguiente descomposición es satisfecha:
!
x2 + y2 + 1 1 1 1
= + .
( y2 − x2 + 1)2
p p
2 ( x − 1 + y2 )2 ( x + 1 + y2 )2

17
18

Finalmente, la función H ( y) queda determinada al imponer la con-

dición

∂F 2xy
=− 2 + H 0 ( y)
∂y ( y − x2 + 1)2
2xy
= Q( x, y)µ ( y2 − x2 ) = − ,
( y2 − x2 + 1)2

que nos conduce a H ≡ k ∈ R. Por tanto, la solución general y( x)

tiene las dos siguientes ramas:
p
y( x) = ± x2 + Cx − 1 , C ∈ R .

8. Responde razonadamente a las siguientes cuestiones:

(a) ¿Puede ser la función φ(t) = t2 (definida para t ∈ R) solución

de una ecuación lineal de primer orden homogénea? ¿Y de una
ecuación lineal de primer orden no homogénea?

(b) ¿Pueden ser las funciones φ(t) = et y ψ(t) = e−t (definidas

para t ∈ R) soluciones de una misma ecuación lineal de primer
orden homogénea? ¿y de una ecuación lineal de primer orden
no homogénea?

En caso afirmativo proporciona ejemplos explícitos.

Solución : (a) La función φ(t) = t2 (definida en R) NO puede ser

solución de una ecuación lineal de primer orden homogénea x0 (t) =
a(t) x(t), ya que de serlo habría de ocurrir

2
2t = a(t)t2 ⇒ a(t) = ∀ t 6= 0
t
y, en cualquier caso, el coeficiente a(t) sólo estaría definido en R \ {0}
(no en R).
Sin embargo, la respuesta es SI para el caso no homogéneo x0 (t) =
a(t) x(t) + b(t). En efecto, bastaría con encontrar funciones continuas

18
Métodos elementales 19

a, b : R → R tales que 2t = a(t)t2 + b(t). Por ejemplo, a ≡ 0 y

b(t) = 2t.
(b) Las funciones φ(t) = et y ψ(t) = e−t (definidas en R) NO pueden
ser soluciones de una misma ecuación diferencial lineal de primer
orden homogénea, ya que de serlo tendrían que verificarse simultá-
neamente las dos siguientes condiciones:

et = a(t)et , − e−t = a ( t ) e−t ,

que a la postre se traducen en que simultáneamente a ≡ 1 y a ≡ −1.

Para el caso no homogéneo los análogos de las condiciones anterio-
res son
et = a ( t ) et + b ( t ) , − e−t = a ( t ) e−t + b ( t ) ,
que tras un cálculo sencillo conducen a las siguientes expresiones de
los coeficientes:
et + e−t 2
a(t) = , b(t) = − .
et − e−t et − e−t
Estas expresiones son sólo válidas en R \ {0}, luego φ(t) y ψ(t) NO
pueden ser soluciones de una misma ecuación lineal de primer orden
no homogénea.


9. Calcula los valores de la constante µ ∈ R que hacen que la ecuación

diferencial x0 = (µ + cos2 t) x tenga una solución π-periódica5 no
trivial.

Solución : Se trata de una ecuación diferencial en variables separadas,

por lo que se puede encontrar fácilmente su solución general6 :
1 1
x(t) = Ce(µ + 2 )t+ 4 sen(2t) .
5 Esdecir, x(t) = x(t + π ) ∀ t ∈ R
6 Parallegar a la expresión final de la misma es necesario calcular previamente una
primitiva de cos2 (t), tarea la cual resulta bastante simple si se escribe

1 + cos(2t)
cos2 (t) = .
2

19
20

Al evaluarla en t + π se obtiene
1 1 1 1
x(t + π ) = Ce(µ + 2 )(t+π )+ 4 sen(2t+2π ) = Ce(µ + 2 )(t+π )+ 4 sen(2t) ,

luego la condición de periodicidad equivale a imponer

1
e(µ + 2 )π = 1 ,

es decir, µ = − 12 .


10. Resuelve las siguientes ecuaciones diferenciales lineales usando el

método de variación de las constantes:

(a) x0 − tx = 3t

(b) y0 = 5y + cos( x)
 
0 2t
(c) x = t2 +1 x + t3

(d) x0 − 2x = 4e2t

Solución : (a) Resolviendo en primer lugar la correspondiente ecua-

ción diferencial homogénea, x0 = tx, la cual tiene sus variables sepa-
radas obtenemos
t2
xh (t) = C e 2 , C ∈ R .
Asumimos ahora C = C (t) y reconsideramos la ecuación diferencial
completa. En este caso se tiene que
t2
 t2 
x0 = (C 0 + Ct)e 2 = t Ce 2 + 3t ,

t2
luego C 0 (t) = 3te− 2 . Por tanto
t2
C (t) = −3e− 2 + K , K ∈ R.

(b) La solución general de la ecuación homogénea, y0 = 5y, es

yh = Ce5x , C ∈ R.

20
Métodos elementales 21

Variando la constante y considerando ahora la ecuación completa

obtenemos

(C0 + 5C )e5x = 5Ce5x + cos( x) ⇒ C0 ( x) = cos( x) e−5x ,

luego ha de ser
1 −5x
C ( x) = e (sen( x) − 5 cos( x)) + K , K∈R
26
y, por tanto,
1
y( x) = (sen( x) − 5 cos( x)) + K e5x , K ∈ R.
26

(c) La solución general de la ecuación homogénea, x0 = 2t

t2 +1
x, es

xh = C (t2 + 1) , C ∈ R.

Variando la constante y considerando ahora la ecuación completa

obtenemos
t3
C 0 (t2 + 1) + 2Ct = 2Ct + t3 ⇒ C 0 (t) = ,
t2 + 1
luego ha de ser
1 2 
C (t) = t − log(t2 + 1) + K , K∈R
2
y, por tanto,
1 2 
x(t) = t − log(t + 1) + K (t2 + 1) ,
2
K ∈ R.
2

(d) La solución general de la ecuación homogénea, x0 − 2x = 0, es

xh = Ce2t , C ∈ R.

Si asumimos C = C (t) y recuparamos la ecuación completa, resulta

que

(C0 + 2C − 2C ) e2t = 4 e2t ⇒ C0 ≡ 4 ⇒ C (t) = 4t + K , K ∈ R .

Entonces
x(t) = (4t + K ) e2t , K ∈ R.

21
22

11. Resuelve los siguientes problemas de valores iniciales mediante la

fórmula de variación de las constantes:

(a) x0 + 3x = e−3t , x(1) = 5

(b) x0 − x
t = 1
1+t2
, x(2) = 0

(c) x0 = cosh(t) x + esenh(t) , x(0) = 1

Solución : (a) La solución general de la ecuación homogénea, x0 +

3x = 0, es
xh (t) = C e−3t , C ∈ R .
Variando la constante, x(t) = C (t) e−3t , y volviendo a la ecuación
completa obtenemos

(C0 − 3C ) e−3t + 3C e−3t = e−3t ⇒ C (t) = t + K , K ∈ R.

Por tanto, la solución general de la ecuación completa es

x(t) = (t + K ) e−3t , K ∈ R.

Imponiendo la condición inicial x(1) = 5 concluimos que

(1 + K ) e−3 = 5 ⇒ K = 5e3 − 1 ,

luego
x(t) = (t + 5e3 − 1) e−3t .

(b) Resolviendo la ecuación homogénea (que tiene sus variables se-

paradas) obtenemos xh (t) = Ct con C ∈ R. Sustituyendo entonces la
expresión x(t) = C (t)t en la ecuación completa llegamos a

1 1
C0 t + C − C = 2
⇒ C 0 (t) =
1+t t(1 + t2 )
 
t
⇒ C (t) = log √ +K, K ∈ R.
1 + t2
Entonces la solución general de la ecuación completa es
   
t
x(t) = log √ +K t, K ∈ R.
1 + t2

22
Métodos elementales 23

Imponiendo finalmente la condición inicial obtenemos

  2    2 
2 log √ + K = 0 ⇒ K = − log √ .
5 5
Por consiguiente, la solución buscada es
√ !
5t
x(t) = t log √ .
2 1 + t2

(c) La solución general de la ecuación homogénea es

xh (t) = C esenh(t) , C ∈ R.

Si variamos la constante y ensayamos en la ecuación completa con

x(t) = C (t) esenh(t) obtenemos

(C0 + C cosh(t)) esenh(t) = C cosh(t) esenh(t) + esenh(t)

⇒ C 0 ≡ 1 ⇒ C (t) = t + K , K ∈ R.

Por tanto
x(t) = (t + K ) esenh(t) , K ∈ R.
Al imponer finalmente x(0) = 1 obtenemos K = 1 y, en consecuen-
cia, la única solución de nuestro problema de valores iniciales es

x(t) = (t + 1) esenh(t) .

12. Sean a, b : R → R funciones continuas tales que a(t) ≥ c > 0 para

todo t ∈ R y
lı́m b(t) = 0 .
t→∞

Demuestra que todas las soluciones de la ecuación diferencial

x0 = − a(t) x + b(t)

tienden a cero cuando t → ∞. (Indicación: usa la regla de L’Hôpital

en el segundo término de la fórmula de variación de las constantes).

23
24

Solución : Resolvemos en primer lugar la ecuación homogénea, x0 =

− a(t) x, obteniendo
Rt
− a(s) ds
x(t) = Ce t0
, C ∈ R.

Variando ahora la constante, es decir, asumiendo C = C (t), conclui-

mos que la solución general de la ecuación lineal x0 = − a(t) x + b(t)
es
Z t  R
− tt a(s) ds
Rs
− t a(s) ds
R
t a(τ ) dτ
x(t) = x0 e 0 + b(s) e 0 ds e t0 ,
t0

donde hemos asumido x(t0 ) = x0 con t0 y x0 arbitrarios.

Estudiamos en primer lugar el límite en infinito del primer término
de la solución:

− tt a(s) ds 
n R o n Rt o
− t c ds
0 ≤ lı́m  x0 e  ≤ | x0 | lı́m e
 0 0
t→∞ t→∞
n o
= | x0 | lı́m e−c(t−t0 ) = 0 .
t→∞

Para el segundo término obtenemos

R Rs   
 t b(s) e t0 a(τ ) dτ ds   b(t) e tt0 a(s) ds 
R
t0
lı́m Rt = lı́m
t→∞   t→∞  a(t) e tt0 a(s) ds 
R
a(s) ds
e t0
 
b(t)
= lı́m = 0,
t→∞ a (t )
 
ya que lı́mt→∞ {b(t)} = 0 y  a(1t)  ≤ 1c para todo t ∈ R. Obsérvese
 

que hemos utilizado la regla de L’Hôpital y el teorema fundamen-

tal del cálculo para resolver la eventual indeterminación en el límite
anterior.


13. Resuelve las siguientes ecuaciones diferenciales de Bernoulli:

t3
(a) 3tx0 − 2x = x2

(b) x0 = et x7 + 2x

24
Métodos elementales 25

y log( x)
(c) y0 + x = x y2

Solución : (a) Dividiendo la ecuación por 3t, t 6= 0, obtenemos

0 2 t2 −2
x − x= x . (1.4)
3t 3
Multiplicando ahora la ecuación (1.4) por x2 llegamos a
2 3 t2
x2 x0 − x = ,
3t 3
que es equivalente a la ecuación

0 2 t2
z − z= (1.5)
t 3
3
vía el cambio de variable z = x3 . La solución general de la ecuación
(1.5) es
2 t
 
z(t) = t +C , C ∈ R,
3
luego deshaciendo el cambio de variable obtenemos
1
x(t) = (t3 + Ct2 ) 3 , C ∈ R.

(b) Multiplicando la ecuación por x−7 obtenemos x−7 x0 = et + 2x−6 .

Planteamos el cambio de variable z = 16 x−6 , en cuyo caso la ecuación
anterior se puede reformular del siguiente modo:
z0 + 12z = − et .
La solución general de esta ecuación es
1 t
e + C e−12t , C ∈ R .
z(t) = −
13
Deshaciendo el cambio de variable obtenemos
 1
6 t −6

−12t
x(t) = C e − e , C ∈ R.
13
log( x)
(c) Multiplicando la ecuación por y−2 obtenemos y−2 y0 = 1
− xy . x
1
Planteamos el cambio de variable z = − y , en cuyo caso la ecuación
anterior se puede reformular del siguiente modo:
z log( x)
z0 = + .
x x

25
26

La solución general de esta ecuación es

z( x) = Cx − log( x) − 1 , C ∈ R.

Deshaciendo el cambio de variable obtenemos

1
y( x) = , C ∈ R.
1 + log( x) − Cx

14. Resuelve las siguientes ecuaciones diferenciales de Riccati usando la

solución particular proporcionada:

(a) y0 = y2 − xy + 1 , y p ( x) = x

(b) y0 = x−4 − y2 , y p ( x) = 1
x − 1
x2

Solución : (a) Efectuamos el cambio de función incógnita

1
u( x) = ,
y( x) − x
de modo que
 0
u0 1
− 2 = = y0 − 1 = y2 − xy
u u
1 2 1  1 x
= +x −x +x = 2 +
u u u u
o, equivalentemente,
u0 = − xu − 1
que es una ecuación diferencial lineal de primer orden cuya solución
general es
 
x2 x2
Z
u( x) = C − e 2 dx e− 2 , C ∈ R .

Deshaciendo el cambio de variable concluimos que la solución de la

ecuación de Riccati de partida es
 −1
x2 x2
Z
y( x) = x + e 2 C − e 2 dx , C ∈ R.

26
Métodos elementales 27

(b) Efectuamos el cambio de función incógnita

1
u( x) = ,
y( x) − 1x + 1
x2

de modo que
 0
u0 1 1 2 1 1 2
− 2 = = y0 + 2 − 3 = 4 − y2 + 2 − 3
u u x x x x x
1  1 1 1  2 1 2 1 2 11
= 4− + − + 2− 3 =− 2− 1−
x u x x2 x x u x x u
o, equivalentemente,
2 1
u0 = 1 + 1− u
x x
que es una ecuación diferencial lineal de primer orden cuya solución
general es
 2 1
u( x) = x2 Ce x + , C ∈ R.
2
Deshaciendo el cambio de variable concluimos que la solución de la
ecuación de Riccati de partida es
!
1 1
y( x) = 2 x − 1 + 2
, C ∈ R.
x Ce x + 1 2

15. Se considera la ecuación diferencial de Riccati

1 y
y0 = − − + y2 . (1.6)
x2 x
(a) Busca una solución particular de la forma y = xα con α ∈ R.
(b) Encuentra la solución que cumple y(1) = 2 y calcula, si es po-
sible, el límite cuando x → ∞ de dicha solución.

Solución : (a) La función y = xα resuelve la ecuación diferencial (1.6)

si y solamente si

αxα −1 = − x−2 − xα −1 + x2α ,

27
28

luego ha de ser α = −1.

(b) Para la solución particular encontrada en (a), y0 ( x) = 1x , conside-
ramos el cambio de función incógnita
1
u( x) = . (1.7)
y ( x ) − x−1
Entonces
 0
u0 1 1 y
− 2 = = y0 + 2 = y2 −
u u x x
 1 1 2 1  1 1  11 1
= + − + = + ,
u x x u x u u x
luego resolver el problema de valores iniciales asociado a (1.6) con
dato inicial y(1) = 2 equivale a resolver la siguiente ecuación dife-
rencial lineal de primer orden
u
u0 + = −1
x
con dato inicial u(1) = 1. La (única) solución de este último proble-
ma se puede calcular fácilmente, obteniéndose

1  3 − x2 
u( x) = .
2 x
Por tanto, deshaciendo el cambio de variable (1.7) llegamos a la so-
lución del problema de Riccati de partida:

x2 + 3
y( x) = .
x(3 − x2 )

Finalmente, por la propia definición de solución de un problema de

valores iniciales concluimos que en nuestro caso no√ se puede tomar
lı́m x→∞ y( x) porque y( x) tiene una asíntota en x = 3.


16. Resuelve las siguientes ecuaciones diferenciales usando el método

que convenga en cada caso:

(a) 3et tan( x) + (2 − et ) sec( x)2 x0 = 0

28
Métodos elementales 29

(b) (3t2 x + x3 ) x0 + 2t3 = 0

(c) xy0 + y = y2 log( x)

(d) t cos(t + x) + sen(t + x) + t cos(t + x) x0 = 0

(e) (3tx + x2 ) + (3tx + t2 ) x0 = 0

(f) tx cos(tx) + sen(tx) + (t2 cos(tx) + e x ) x0 = 0

2
(g) 3x + 3et x 3 + tx0 = 0

(h) x0 = 1
x+t
y
(i) y0 = 2y log( y)+ y− x

(j) y0 = 1 + e2x
y
(k) y0 = x (log( y) − log( x))

(l) xy0 + y = 2x

(m) y0 = log( x y )
h i
C
Solución : (a) Variables separadas: x(t) = arctan (2−et )3
, C ∈ R.
3
(b) Homogénea: x0 = − 3t22t
x+ x3
2
= − 3(x/t)+( x/t)3
. La solución general
viene dada por la siguiente relación implícita7 :

( x + t)3 ( x4 + 3t2 x2 + 2t4 ) = C ( x + 2t)3 , C ∈ R.

y log( x)
(c) Bernoulli: y0 + x = x y2 . La solución general es8

1
y( x) = , C ∈ R.
Cx + log( x) + 1
7 El siguiente cálculo es útil:

u du 1 dv
Z Z
4 2
=
u + 3u + 2 2 v2+ 3v + 2
Z 
1 dv dv 1
Z 
= − = log(v + 1) − log(v + 2) .
2 v+1 v+2 2

8 La
R log( x)dx
integral x2
se resuelve por partes

29
30

(d) Exacta con

P(t, x) = t cos(t + x) + sen(t + x) , Q(t, x) = t cos(t + x) .

En efecto,
∂P ∂Q
= cos(t + x) − t sen(t + x) = .
∂x ∂t
La solución general viene dada por9
C
x(t) = arc sen −t, C ∈ R.
t

2 3( x/t)+( x/t)2
(e) Homogénea: x0 = − 3tx +x
3tx+t2
= − 3(x/t)+1 . La solución general
viene dada por la siguiente relación implícita:

tx
= C, C ∈ R.
(t − x)4

(f) Exacta con

P(t, x) = tx cos(tx) + sen(tx) , Q(t, x) = t2 cos(tx) + e x .

En efecto,
∂P ∂Q
= 2t cos(tx) − t2 x sen(tx) = .
∂x ∂t
La solución general viene dada por la siguiente relación implícita10 :

e x + t2 cos(tx) = C , C ∈ R.

t 2
(g) Bernoulli: Para t 6= 0 se tiene que x0 + 3t x = − 3et x 3 . La solución
general es
 C − et 3
x(t) = , C ∈ R.
t
(h) Haciendo el cambio de variable u = x + t obtenemos la siguiente
ecuación diferencial equivalente:

u+1
u0 = ,
u
9
R
t cos(t + x) dt = t sen(t + x) + cos(t + x) por partes
10
R t 1
t cos(tx) dt = x sen(tx) + x2
cos(tx) por partes

30
Métodos elementales 31

que es una ecuación en variables separadas cuya solución general es

u − log(|u + 1|) = C + t , C ∈ R.
Deshaciendo finalmente el cambio de variable obtenemos la siguien-
te expresión implícita para la solución:
x = log(| x + t + 1|) + C , C ∈ R.

(i) Exacta con P( x, y) = − y y Q( x, y) = 2y log( y) + y − x. En efecto,

∂P ∂Q
= −1 = .
∂y ∂x
La solución general viene dada por
y( x) = y( y log(| y|) − x) = C , C ∈ R.

(j) Variables separadas: Integrando en los dos miembros con respecto

a x obtenemos
1
y( x) = x + e2x + C , C ∈ R .
2
 
y y
(k) Homogénea: y0 = x log x . La solución general viene dada por
la siguiente expresión:
y( x) = x eCx+1 , C ∈ R.
y
(l) Lineal de primer orden: Si x 6= 0 se tiene y0 + x = 2. La solución
general de la correspondiente ecuación homogénea es
C
yh ( x) = , C ∈ R.
x
Usando el método de variación de las constantes encontramos que
la solución general de la ecuación completa es
C
y( x) = x + , C ∈ R.
x
y0
(m) Variables separadas: y = log( x), de donde

y( x) = Ce x(log(x)−1) , C ∈ R.
En particular, esta ecuación admite la solución trivial y ≡ 0.


31
32

17. Decide de forma razonada si las siguientes afirmaciones son verda-

deras o falsas:

(a) Una solución del problema de valores iniciales

 0
x = x2 + t2
x(1) = 2

definida en un intervalo abierto que contenga a [1, 2] satisface

x(2) = 1.

(b) La ecuación de Riccati y0 + y + y2 + e x = 0 se transforma en

una ecuación diferencial lineal de orden dos,

z00 + a( x) z0 + b( x) z + c( x) = 0 ,
z0
mediante el cambio de variable y = z.
(Septiembre 2003)

Solución : (a) FALSA. Claramente x0 ≥ 0, por lo que cualquier solu-

ción es creciente y si ha de satisfacer x(1) = 2 no puede ser x(2) = 1,
ya que en ese caso decrecería.
(b) VERDADERA. Derivando el cambio de variable obtenemos
 0 2
0 z00 z
y = − ,
z z

luego la ecuación resultante en la nueva función incógnita z( x) es

z00 + z0 + e x z = 0 ,

que es lineal de segundo orden.

32