Soluciones del examen de feb05 de Ecuaciones Diferenciales I (C)

 0
x =y x (0) = 1
0
1. Hallar la solución del sistema y = 4x + z con y(0) = 1 , y estudiar su estabilidad.
 0
z = −z + 4e−2t z (0) = −4

z = Ce−t − 4e−2t

3a desacoplada: z0 = −z+4e−2t , zp = Ae−2t → −2A = − A+4 → → z = −4e−2t
z (0) = −4
0 −2t
x0 = y 00 − 4x = −4e−2t , x = c e2t + c e−2t + x , x = Ate−2t , x p = A(1 − 2t)e

, x 1 2 p p → A=1
y0 = 4x −4e−2t x 00p = A(4t − 4)e−2t

− − x (0) = c1 + c2 = 1
2t 2t
x = c1 e + c2 e + te , 2t → x = (1 + t)e−2t → y = x 0 = (−1 − 2t)e−2t .
x 0 (0) = 2c1 − 2c2 = −1

sX − 1 = Y


 ↓
2 (s+2)(s+3)
Laplace: sY + 1 = 4X + Z s X − s + 1 = 4X − s+4 2 , (s + 2)(s − 2) X = s s++s2−6 =
2
s +2 ,
 sZ + 4 = − Z + 4 → Z = −4s−4 = − 4 ↑ → z = −4e−2t

s +2 (s−2)(s+3) s +2

X= s +3
( s +2)2
= 1
s +2 + 1
( s +2)2
→ x = (1 + t)e−2t , y = x 0 ↑

 
Matrices 3x3: −λ 1 0   
−1 1 1

4 −λ 2
0  = −(λ + 1)(λ − 4) = 0 , λ = −1, −2, 2 → P = 1 −2 2,
(lo más largo)  1 0 −1 − λ 3 0 0

 

e− t 0 0
 
−16e−t

Z t 4et−s
   
0 0 4 1+t
− 1 1 Jt
P = 12 6 −3 3, e = 0 e − 2t 1 − 2t
0 , x = 12 P −3e + 3 P1 − 2t
 3e  ds = e − 2t −1 − 2t
6 3 1 0 0 e2t −e2t 0 e2t−4s −4

Como ∃λ = 2 (con Reλ > 0 ), todas las soluciones del sistema y esta en particular son inestables.
(Que la solución calculada x(t) → 0 cuando t → ∞ no tiene nada que ver).

4. Sea t2 (1 + t) x 00 + t(3 + 2t) x 0 + x = 0 .

a] Hallar, trabajando por series en t = 0 , una solución no trivial que no contenga el ln t .
b] Probar que todas sus soluciones están acotadas cuando t → ∞ .
∞
a] t = 0 singular regular con r (r −1)+3r +1 = 0 , r = −1 doble. Sin ln t tenemos x1 = ∑ ck tk−1 →
n =0
∞
t k −1 tk t k −1 tk t k −1
 
∑ (k − 1)(k − 2)ck + (k − 1)(k − 2)ck + 3(k − 1)ck + 2( k − 1) c k + c k =0→
k =0

t−1 : c0 − 3c0 + c0 = 0 → c0 indeterminado; t0 : 2c0 − 2c0 + c1 = 0 , c1 = 0 ;

1
k2 ck +(k − 1)(k − 2)ck−1 = 0 , ck en función del anterior → x1 = t .

3+2t 1
− 3t dt → x2 = 1t t +1
= 1 + lnt t . x1 y x2 acotadas cuando t → ∞ .
R R  R
b] − t (1+ t )
dt = t +1 t dt

(Calculando la x2 de Frobenios llegaríamos, con más trabajo, a esa misma solución).

Para responder a b] no se necesita la solución, basta estudiar el punto del infinito:
t = 1s → (1 + 1s )(s2 ẍ + 2s ẋ ) + (3 + 2s )(−s ẋ ) + x = 0 → s(1 + s) ẍ − s ẋ + x = 0 ,
s = 0 es singular regular con r = 1, 0 → x1 = s ∑ · y x2 = ∑ · + dx1 ln s
ambas acotadas cuando s → 0+ ( t → ∞ ), pues s ln s → 0 cuando s → 0+ .
 2. Resover x 00 + x = 2(cos t)−3 con x (0) = 0 , x 0 (0) = 0 .

λ2 + 1 = 0 , λ = ±i → x = c1 cos t + c2 sen t + x p .
 
 cos t sen t 2
|W|(t) = − sen t cos t  = 1 → x p = s c22 − c cs3 = 2sc − cc2 = cos1 t − 2 cos t
R R


1 x (0) = c1 + 1 = 0 1 sen2 t
→ x = c1 cos t + c2 sen t + → → x= − cos t =
cos t x 0 (0) = c2 = 0 cos t cos t

dy
3a. Sea la ecuación [e] (1 + x3 y) + ( x4 + x3 y) dx = 0 . a] Hallar la solución general de [e]
sabiendo que tiene un factor integrante que sólo depende de x . b] Precisar cuántas
soluciones de [e] cumplen y(1) = 1 y dar su expresión lo más simplificada posible.
dy
a] (1 + x3 y) g( x ) + ( x4 + x3 y) g( x ) dx es exacta si x3 f ( x ) = (4x3 + 3x2 y) g( x ) + ( x4 + x3 y) g0 ( x )
−3x3 −3x2 y
R
es decir, si g0 ( x ) = x4 + x3 y
g( x ) = − 3x g( x ) → g( x ) = e− (3/x )dx = 1
x3
factor integrante.

1 dy H = − 2x1 2 + xy + p(y) 1
x3
+ y + ( x + y) dx exacta → y2 → H = y2 + 2xy − x2
= C solución general.
H = xy + 2 + q( x )
−1− x 3 y
b] Solución única, pues f ≡ x4 + x3 y y f y continuas en entorno de (1, 1) . La calculamos:
√ ( y = −2x − 1x
y(1) = 1 → 1+ 2−1 = 2 = C → y = − x ± x2 + 2 + x −2 = − x ± ( x + 1x ) → y = 1x
no lo cumple)

x 0 = x4 + x3 y

3b. Dibujar el mapa de fases de [S] . (Utilizar los cálculos realizados en 3a.).
y 0 = −1 − x 3 y

Puntos críticos: x3 ( x + y) = 0 → x = 0, −1 = 0 imposible, ó

1

  
−1
 3
4x + 3x2 y x3
 y
y = − x, x4 = 1 → −1 , 1 . AL: .
−3x2 y − x3
           
1
−1
: ±2 → 11 , −13 y −11 : ±2 → −13 , 11 .
1
Ambos puntos son sillas (del lineal y del no lineal).
Vertical: y = − x , y = 0 x
. Algunos valores del campo:
Horizontal: y = −1/x3 –1 1

4
   
1
v( x, 0) = x , v(±1, y) = (1 ± y) −1 .
−1
–1
En 3a hemos hallado la expresión de las órbitas.
Las separatrices son la curvas que pasan por
(1, −1) y (−1, 1) . Para ambos es C = −2 →
√
y = − x ± x2 − 2 + x −2 = − x ± ( x − 1x ) = − 1x , −2x + 1x .
Y tenemos además las otras dos curvas halladas en 3a.
 Soluciones del examen de sept05 de Ecuaciones Diferenciales I (C)
1−2ty3
1. Sea la ecuación y0 = 3t2 y2 .
a] Hallar su solución general por dos caminos diferentes.
b] Precisar cuántas soluciones cumplen y(1) = 1 y escribirlas explícitamente.

z = y3 R t2 1/3
a] Bernouilli: 3y2 y0 = − 2t y3 + t12 → z0 = − 2z 1 C 1
t + t2 → z = t2 + t2 t2
→ y = [ Ct+2 t ] ;
U = t2 y3 − t + p ( y )
Exacta: 2ty3 − 1 + 3t2 y2 y0 = 0 , My = Nt = 6ty2 , → t2 y3 − t = C ↑ .
U = t2 y3 + q ( t )

b] Como f y f y son continuas en un entorno de (1, 1) , hay solución única con ese dato.
1
Imponiéndolo en la solución recuadrada: 1 = [C + 1]1/3 → C = 0 → y = t1/3 .

[Si se impone en la solución implícita de la exacta, se tiene: C = 0 → t = 0 ó ty3 − 1 = 0 ;

la primera es curva integral que ni es solución ni pasa por (1, 1), y la segunda es la de arriba].

2. Sea x IV + 4x 000 + cx 00 + 4x 0 + x = 16 et .
a] Discutir su estabilidad según los valores de la constante c .
b] Para c = 6 , hallar la solución que satisface x (0) = x 00 (0) = 0 , x 0 (0) = x 000 (0) = 2 .
c] Hallar una solución de la ecuación para cada uno de los valores de c .

a] P(λ) = λ4 + 4λ3 + cλ2 + 4λ + 1 es el polinomio característico.

si c < 0 es I
Del signo de los coeficientes deducimos que [y si c > 0 no sabemos nada].
 
si c ≤ 0 no es AE
4 1 0 0
4 c 4 1
B= 
0 1 4 c
 → 4 > 0 , c > 1 , c > 2 , 1 > 0 → AE ⇔ c > 2 , c < 2 ⇒ I .
0 0 0 1

q4 − 2q2 + 1 = 0
Si c = 2 , λ = ±qi → → λ = ±i → P(λ) = (λ2 + 1)(λ2 + 4λ + 1)
4q(1 − q2 ) = 0
√
→ λ = ±i simples y λ = −2 ± 3 con Reλ < 0 ⇒ EnoA
4
b] λ4 + 4λ3 + 6λ2 + 4λ + 1 = (λ + 1) ; x p = Aet , 16A = 16 →
d.i.
x = ( c 1 + c 2 t + c 3 t 2 + c 4 t 3 ) e− t + et → · · · · · · x = et − e− t .
16
O Laplace: s4 X − 2s2 − 2 + 4s3 X − 8s + 6s2 X − 12 + 4sX + X = s −1 ,
2( s +1)3
(s + 1)4 X = 2s2 + 8s + 14 + s16
−1 = s −1 →X= 2
(s+1)(s−1)
= 1
s −1 − 1 ↑
s +1

c] λ = 1 es autovalor si 10 + c = 0 ⇒ x p = Aet , (10 + c) A = 16 → x p = 1016 t

+c e , si c 6 = −10 .

Y si c = −10 [ P(λ) = (λ + 1)2 (λ2 + 6λ + 1) ] , x p = At2 et , x 0p = A(t2 + 2t)et , x 00p = A(t2 + 4t + 2)et ,
x 000 2 t IV 2 t t t 2 t
p = A ( t + 6t + 6)e , x p = A ( t + 8t + 12)e → 16Ae = 16e → x p = t e
3. Sea (1 + t2 ) x 00 + tx 0 − x = 0 .
a] Hallar hasta orden t6 el desarrollo en serie de la solución con x (0) = 1 , x 0 (0) = 0 .
b] Probar que posee soluciones que tienden hacia 0 cuando t → ∞ .
√
1+ t2
R dt
[Por si a alguien le es útil: 2
√
2
= − t ].
t 1+ t

∞ ∞ ∞ ∞
a] t = 0 regular, x = ∑ ck tk → ∑ [k(k − 1)ck tk−2 + k(k − 1)ck tk ] + ∑ kck tk + ∑ −ck tk = 0 .
k =0 k =2 k =1 k =0
1
t0 : 2c2 − c0 = 0 → c2 = 2 c0 ; t1 : 6c3 + c1 − c1 = 0 → c3 = 0 ;
t2 : 12c4 + 2c2 + 2c2 − c2 = 0 → c4 = − 14 c2 = − 18 c0 ;
−1
tk : (k + 2)(k + 1)ck+2 + (k2 − 1)ck = 0 → ck+2 = − kk+ 2 ck ,
→ c5 = c7 = · · · = 0 , c6 = − 36 c4 = 1
16 c0 , ...
c0 =1, c1 =0
→ x = c0 [1 + 21 t2 − 18 t4 + 16
1 6
t + · · · ] + c1 t → x = 1 + 21 t2 − 18 t4 + 16
1 6
t +··· .
−
R t √ √
1+ t2 = √ 1 √dt
R
O bien, x1 = t , e 2
, x 2 = t = − 1 + t2 → x = c1 t + c2 1 + t2 .
1+ t t2 1+ t2
1/2 (1/2)(−1/2) 4 (1/2)(−1/2)(−3/2) 6
Con los datos iniciales: x = (1 + t2 ) = 1 + 21 t2 + 2! t + 3! t +···
√ −C
b] A partir de la solución general es claro que las soluciones C (t − 1 + t2 ) = √ → 0.
t + 1+ t2 t → ∞
2
O bien, t = 1
s → s4 (1 + s12 ) ẍ + 2s3 (1 + s12 ) ẋ − ss ẋ − x = s2 (1 + s2 ) ẍ + s(1 + 2s2 ) ẋ − x = 0 ,
para la que s = 0 es punto singular regular con r (r − 1)+ r − 1 = 0 , r = ±1
⇒ s ∑ ck sk son soluciones que tienden a 0 cuando s → 0+ ( t → ∞ ).
[La serie de t = 0 converge si |t| < 1 (donde asegura el teorema) y no informa sobre el infinto].

a] Hallar las órbitas y dibujar el mapa de fases de [S].

x 0 = x2 y

4. Sea [S] . [Comprobar que una separatriz pasa por el punto (3, 4)
y0 = x4 − 1 y que hay una órbita que pasa por (−3, 0) y (−1, 4) ].
b] Hallar la solución de [S] que cumple x (1) = y(1) = 0 .
 
1
2xy x2   + , silla, ±2 → ±2
     
±1 0 1
a] Puntos críticos: . A. lineal: = →
0 4x3 0 ±1 ±4 0 – , centro de la AL
0
f ( x, −y) = − x2 y = − f ( x, y) y
 
−1 4
4 ⇒ centro también del no lineal.
g( x, −y) = x − 1 = g( x, y) 0
dy x 4 −1 2x3 2
Órbitas (separable): y dx = x2
→ y2 = 3 + x +C .
[Confirmamos simétría de las órbitas son respecto a y = 0 ].
x
q 
0 = 32 + 2 + C 2x3 2 8 –3 –1 1 3

Separatrices: (1, 0) → → y =± 3 +x−3

 = ±4 .
x =3
Por (−3, 0) y (−1, 4) pasa la misma órbita cerrada la de C = 56


3 .
Vertical: y = 0 , x = 0 . 2x2 ( x − 1)
     
8 9
v( x, 2( x − 1)) = → ,4 , si x = 2, 3 .
Horizontal: x = ±1 . x4 − 1 15 20

y(1)=0
 
0
b] x (1) = y(1) = 0 → órbita x = 0 → y0 = −1 , y = K − t → x = .
1− t
 Soluciones del examen de feb06 de Ecuaciones Diferenciales I (C)
2t(1−y)
1. Sea la ecuación y0 = y+t2 .
a] Hallar la solución (lo más simplificada que se pueda) con: i) y(0) = 1 , ii) y(0) = −1 .
b] Estudiar si existe solución satisfaciendo y(1) = −1 . [2 puntos]

dy U = t2 y − t2 + p ( y )
a] (2ty − 2t) + (y + t2 ) dt exacta: My = Nt = 2t ⇒ y2 → y2 + 2t2 y − 2t2 = C .
U = t2 y + 2 + q(t)
t2 + y z = t2 y C + y2
dt
(Bernouilli) → z0 = 1−
z C 1
y dy = 2(1−y) = t2 ↑ ].
R
[O bien, dy = 2t (1−y) y + 1− y → z = 1− y + 1− y

Tanto para i) como para ii) se tiene: 1 = C → 0 = y2 − 1 + 2t2 y − 2t2 = (y − 1)(y + 1 + 2t2 ) →
La solución de i) es y = 1 (calculable a ojo) y la de ii) es y = −1 − t2
[soluciones únicas, pues f y f y son continuas en un entorno de (0, 1) y (0, −1) ].

b] f no es continua en (1, −1) : el TEyU, aplicado a la ecuación inicial, no nos dice nada.
2
dt
= 2tt(1+−yy) . Se tiene que dy
dt 
= 0 y que tanto 1f

Estudiemos la ‘equivalente’: dy (1,−1)
como 1f t son continuas en un entorno de (1, −1) (el denominador es no nulo).
 

Existe una única solución t(y) con t(−1) = 1 y tiene pendiente t0 = 0 [única curva
integral (de pendiente infinita) por (1, −1)] ⇒ No hay solución y(t) con y(1) = −1 .
√ p
[Llevando el dato a la solución general: C = −3 → y = −t2 ± t4 + 2t2 − 3 = −t2 ± (t2 − 1)(t2 + 3) ,
y ninguna de las dos funciones son derivables en t = 1 ].

2. Sea t(1 + t) x 00 + (2 + 3t) x 0 + x = 0 .

a] Hallar el desarrollo en serie de una solución no trivial que esté acotada en t = 0 .
b] Probar que no todas sus soluciones son analíticas en t = −1 .
c] Comprobar a] y b] utilizando el hecho de que x = 1t es solución. [3 puntos]

∞
a] t = 0 singular regular, r (r − 1) + 2r + 0 = 0 → r = 0, −1 → x1 = ∑ ck tk acotada en t = 0 .
k =0
∞ ∞ ∞
∑ [ k ( k − 1) c k t k −1 + k ( k − 1) c k tk ] + ∑ [2kck t k −1 + 3kck tk ] + ∑ ck tk = 0 .
k =2 k =1 k =0

t0 : 2c1 + c0 = 0 → c1 = − 12 c0 ; t1 : 2c2 + 4c2 + 3c1 + c1 = 0 → c2 = − 23 c1 = 31 c0 ;

+1
tk : [(k + 1)k + 2(k + 1)]ck+1 + [k2 + 2k + 1]ck = 0 → ck+1 = − kk+ 2 ck →
k −1
c3 = − 34 c2 = − 14 c0 , c4 = − 45 c3 = 15 c0 , . . . ; ck = − k+k 1 ck−1 = k +1 c k −2 = ···
log(1+t)
→ x1 = 1 − 21 t + 13 t2 − 41 t3 + · · · + (−1)k k+1 1 tk + · · · = t .

b] Si se identifica la serie, es claro que la solución que define no es analítica en t = −1 .

O bien, s = t + 1 → s(s − 1) x 00 + (3s − 1) x 0 + x = 0 → r = 0 doble →
∞ ∞
x1 = ∑ ck sk analítica, pero x2 = s ∑ bk sk + x1 log s no analítica en s = 0 ( t = −1 ).
k =0 k =0

2+3t
R
− [ 2t + 1+
1
R
log(1+t)
c] Utilizando que x1 = 1t : e = e− t] 1 1 dt
R
t (1+ t ) = t2 (1+ t )
, x2 = t 1+ t = t ,
solución acotada en t = 0 con el desarrollo de arriba y claramente no analítica en t = −1 .
3. a] Hallar una solución de la ecuación x 000 + x 00 + 2x 0 + 8x = eat para: i) a = 1 , ii) a = −2 .
b] Precisar la estabilidad de la solución hallada en cada caso. [1.5 puntos]
a] P(λ) = λ3 + λ2 + 2λ + 8 es el polinomio característico.
1 t
ai] P(1) = 12 ⇒ λ = 1 no es autovalor ⇒ x p = Aet , 12A = 1 → x p = 12 e
aii] P(−2) = 0 ⇒ λ = −2 sí es autovalor ⇒ x p = At e−2t , x 0p = A(1 − 2t) e−2t ,
x 00p = A(4t − 4) e−2t , x 000
p = A (12 − 8t ) e
−2t → 10A e−2t = e−2t → x =
p
1
10 t e
−2t

b] La estabilidad sólo depende de la ecuación homogénea, en concreto de las raíces de P(λ) .

√
λ = −2 autovalor → P(λ) = (λ + 2)(λ2 − λ + 4) → λ = −2 y λ = 12 ± i 215 .
Hay autovalores con Reλ > 0 ⇒ todas las soluciones son inestables (para cualquier a ).
 
 1 1 0 
Lo confirma R-H: B = 8 2 1 → 1 > 0 , −6 < 0 ⇒ I. El signo de los coeficientes no decía nada.
0 0 1

4a. Determinar la estabilidad de la solución de x 0 = x + x2 que cumple x (0) = −2 . [1 punto]

x
La ecuación es autónoma y el dibujo basta para precisar la estabilidad.
t
x = 0 , −1 , > 0 ⇒ crecen si x < −1 ó x > 0 . 0
x (1 + x )
soluciones constantes. < 0 ⇒ decrecen si −1 < x < 0 .
–1

La solución con x (0) = −2 y todas las cercanas están definidas

para todo t ≥ 0 y tienden a −1 cuando t → ∞ ⇒ es AE.
–2

x 0 = x + x2

4b. Sea [S] . i) Hallar las órbitas y dibujar el mapa de fases de [S].
y0 = 2x + y    
x (t) −2
ii) Hallar lı́m x (t) y lı́m y(t) , si y(t) es la solución de [S] con x(0) = 3 . [2.5 puntos]
t→∞ t→∞

dy y R dx x
i) Ecuación de las órbitas: dx = x+ x2 + 1+2 x (lineal).
R 1 1

x+x 2 = x − 1+ x dx = log 1+ x
x (C +2 log | x |)
R
→ e a = 1+x x → y = 1Cx 2x
R dx
+ x + 1+ x x = 1+ x .
3
     
0 −1 1 + 2x 0
Puntos críticos: 0 2
y . A. lineal: 2 1
.
    2
0 0
En 0 hay λ = 1 doble asociado a 1 (nodo1tgI).
       
−1 −1 0 0 1 1
2
→
2 1
→ λ = 1,
1
y λ = −1, −1 : silla.
Campo vertical: x = 0,−1 (órbitas). Horizontal: y = −2x . —2 —1
    1
1+ x x
v( x, 0) = x
2
, v( x, 2) = (1 + x )
2
.
 
x
La separatriz estable se deforma: v( x, 1 − x ) = (1 + x )
1
.
log | x | 2x log | x |
Como 1+ x x →−→ −1 [L’Hôpital], esta separatriz es la órbita de C = 0 : y = 1+ x
1
[su valor en x = −2 es: 4 log 2 ≈ 4× 0.69 < 3 ].

ii) lı́m x (t) = −1 , según dice el dibujo de 4a. lı́m y(t) = ∞ , puesto que la órbita asociada
t→∞ t→∞
está por encima de la separatriz [ y(t) , solución de lineal, está definida ∀t ≥ 0 ].
Tanto 4a como 4bii se pueden responder, con bastante más esfuerzo, hallando las soluciones de [S]:
x (0)=−2
x 0 = x + x2 (Bernouilli o separable) → x = 1
Ce−t −1
→ x= 2
→ −1
e− t −2 t → ∞
;
y(0)=3
y0 = y + 4
y = et 3 − 4 log |2 − e−t | → ∞.


e− t −2
→
t→∞
 Soluciones del examen de sept06 de Ecuaciones Diferenciales I (C)
x0 = y −2

x (0) = 0
1. Hallar la solución de con , y precisar su estabilidad. [2.5 puntos]
y0 = 2x − y y (0) = 4

Convirtiendo el sistema en ecuación: x 00 = y0 = 2x − y = 2x − x 0 − 2 , x 00 + x 0 − 2x = −2

λ=1,−2; x p =1 x (0)=0, x 0 (0)=4−2
→ x = c1 et + c2 e−2t + 1 → x = 1 − e−2t → y = 2 + x 0 = 2 + 2 e−2t .
(
sX = Y − 2s Y = sX + 2s ↓
Por Laplace: →
sY − 4 = 2X − Y (s + 1)(sX + 2s )− 4 = 2X → (s + 2)(s − 1) X = 2 − 2
s →
2( s −1)
X= s(s+2)(s−1)
= 2
s ( s +2)
= 1
s − 1
s +2 →Y= 2
s + 2
s +2 → x = 1 − e−2t , y = 2 + 2 e−2t .
   −2t 
 −λ 1 
   
e 0
1 1 − 1 1 1 −1 Jt
Con matrices:   = 0 , λ = −2, 1 → P = ,P = 3 2 1 ,e = .

 2 −1 − λ  −2 1 0 et
t !
e2s−2t 0

−4e−2t −1 − 3e−2t 1 − e−2t
  Z t       
1 Jt −4 1 J ( t − s ) 2 1 1
x = 3 Pe −3P e ds = 3 P − = 3P = .
4et 2 + 2e−2t
t
4 0 4 −4et−s  4
0

Para evitar la integración se puede buscar una x p constante (que podría no existir):
1 − e−2t
         
a 0 = b−2 1 1 1 d.i.
x p = b → 0 = 2a − b → x = c1 −2 e−2t + c2 1 et + 2 → x = .
2 + 2e−2t

Como existe un autovalor positivo λ = 1 todas las soluciones del sistema son inestables,
la calculada en particular (que ella esté acotada, no tiene nada que ver).

2. Sea x 00 + [2 − 2t] x 0 + [1 − 2t] x = 0 .

a] Hallar el desarrollo en serie hasta t4 de la solución que satisface x (0) = 0 , x 0 (0) = 1 .
b] Hallar esta solución en términos de una integral, sabiendo que x = e−t es otra solución.
c] Utilizar b] para comprobar a]. [2.5 puntos]

∞
a] t = 0 regular → llevamos x = ∑ ck tk a la ecuación, teniendo en cuenta que c0 = 0 y c1 = 1 .
k =0
∞ ∞ ∞
∑ k(k − 1)ck tk−2 + ∑ [2kck tk−1 − 2kck tk ] + ∑ [ck tk − 2ck tk+1 ] = 0 →
k =2 k =1 k =0

t0 : 2c2 + 2c1 + c0 = 2c2 + 2 = 0 → c2 = −1 ; t1 : 6c3 + 4c2 − c1 − 2c0 = 6c3 − 5 = 0 → c3 = 65 ;

t2 : 12c4 + 6c3 − 3c2 − 2c1 = 12c4 + 6 = 0 → c4 = − 12 . Tenemos: x = t − t2 + 56 t3 − 21 t4 + · · · .
−2(k+1)ck+1 +(2k−1)ck +2ck−1
[La regla de recurrencia ck+2 = (k+2)(k+1)
no es útil para tan pocos términos].

%
O bien: x 00 (0) + 2x 0 (0) + x (0) = 0 , x 00 (0) = −2 . Y derivando la ecuación:
↑
x 000 + [2 − 2t] x 00 − [1 + 2t] x 0 − 2x = 0 → x 000 (0)+ 2x 00 (0)− x 0 (0)− 2x (0) = 0 , x 000 (0) = 5
↑
xiv + [2 − 2t] x 000 − [3 + 2t] x 00 − 4x 0 = 0 → xiv (0)+ 2x 000 (0)− 3x 00 (0)− 4x 0 (0) = 0 , xiv (0) = −12
2 t 2
R
b] x1 = e−t , e− a = et −2t → x2 = e−t 0 es ds (con límites para fijar la primitiva).
R
Rt 2 Rt 2 −t d.i. Rt 2
x = ce−t + ke−t 0 es ds , x 0 = −ce−t − ke−t 0 + ket → x = e− t 0 es ds .
Rt
c] x = e−t 0 1 + s2 + 12 s4 + · · · ds = 1 − t + 21 t2 − 16 t3 + · · · t + 13 t3 + 10
1 5
   
t +···
= t − t2 + [ 13 + 12 ] t3 − [ 31 + 16 ] t4 + · · · = t − t2 + 56 t3 − 21 t4 + · · ·
 dv 2u − v
3. Sea (E) = . a] Hallar su solución general y la particular que cumple v(1) = 1 .
du v
b] Dibujar aproximadamente sus curvas integrales. [2.5 puntos]

z2 + z −2
a] Homogénea. z = uv → uz0 = 2− z z dz
R R du
z − z = − z , (z+2)(z−1)
= − u +C,
z A B A(z−1)+ B(z+2)
(z+2)(z−1)
= z +2 + z −1 = (z+2)(z−1)
→ A = 32 , B = 1
3 →

2 log(z + 2) + log(z − 1) = C − 3 log u , (z + 2)2 (z − 1) = C

u3
,

( uv + 2)2 ( uv − 1) = C
u3
, (v + 2u)2 (v − u) = C .
[única solución con ese dato, pues
v (1) = 1 → C = 0 → v = u f , f v continuas en un entorno].
2
b] Las isoclinas son v = mu : K = f (u, mu) = m −1.
Hay rectas solución si K = m → m2 + m − 2 = 0 , m = −2, 1
[estas se veían ya en la solución general].
m 2 −1 − 12 12
0
. Además, sobre u = 0 es K = −1 .
K 0 ∞ −3 −5 1
[Las curvas integrales pintadas son las órbitas del sistema lineal u0 = v ; v0 = 2u − v cuyos autovalores
y autovectores se calcularon en el problema 1: confirma que es una silla con esas rectas separatrices].

4. Clasificar los puntos críticos y dibujar el mapa de fases de [e] x 00 = x3 − x − xx 0 .

[Ayudas: Comprobar que hay órbitas que son parábolas.
Haciendo u = 12 [ x2 − 1] en la ecuación de las órbitas se obtiene la (E) del problema 3]. [2.5 puntos]

x0 = v
      
0 1 −1
. Puntos críticos: , y .
v0 = x3 − x − xv 0 0 0
       
0 1  
±1 0 1 0 0 1
A.L. en es ; en , .
3x2 − 1 − v − x 0 2 ∓1 0 −1 0
   
1 −1
→ λ = −2 , 1 ; → λ = 2 , −1 . Sillas.
0 0  
 1 
En toda ecuación el v asociado a λ es .
λ
centro de la aproximación lineal
 
0
→ λ = ±i :
0 (centro o foco de nuestro sistema).
La ecuación no es exacta, pero como g(− x, v) = − g( x, v) ,
las órbitas son simétricas respecto a x = 0 y sigue siendo
un centro del no lineal.
dv x − x3
La ecuación de las órbitas [o] = − x no es de ningún tipo conocido.
dx v
Si tiene soluciones que son parábolas, por la simetría, deben tener la forma v = Ax2 + B :
(
2A2 + A − 1 = 0 → A = 1/2, −1 &
( Ax2 + B)2Ax = x − x3 − x ( Ax2 + B) →
2AB + B + 1 = 0 B = −1/2, 1.
2
Las parábolas halladas y = 1 − x2 e y = x 2−1 son precisamente las separatrices de los puntos silla.
dv dv du dv  2u dv 2u − v
Como u = 12 [ x2 − 1] →

x= = = − 1 x → (E) = ,
du du dx dx v du v
2
deducimos que las órbitas son (v + x2 − 1) (2v − x2 + 1) = C , y aparecen otra vez las separatrices.
Para completar el dibujo obtenemos información del campo v :
x=0
 
v
Campo vertical (como en toda ecuación) en v = 0 . Horizontal en 2 . v(±1, v) = .
v = x −1 ∓ v
Soluciones del examen de febrero de 2008 de Ecuaciones Diferenciales I (A y C)

1. Sea y0 = −3y + 3ty2/3 . Hallar su solución general y una o dos soluciones (si las hay) con:
i) y(1) = 0 , ii) y(0) = 1 . [2 puntos]

z=y1/3
Es una ecuación de Bernouilli. 13 y−2/3 y0 = −y−1/3 + t −→ z0 = −z + t →
3
z = Ce−t + e−t et tdt = Ce−t + t − 1 (ó z p = At + B ...) → y = (Ce−t + t − 1) .
R

Como f es continua en todo R2 y f y = −3 + 2ty−1/3 lo es en R2 −{y = 0}, hay solución única

con y(to ) = yo si yo 6= 0 y hay solución, aunque no sabemos si es única o no, cuando yo = 0 .
Por tanto, con la condición ii) habrá una sola solución. Imponiendo el dato:
3
1 = (C − 1)3 → C = 2 → y = (2e−t + t − 1) .
3
Para i) la solución general nos da una solución: 0 = (Ce−1 ) → C = 0 → y = (t − 1)3 .
Pero hay otra solución clara perdida en los cálculos: y ≡ 0 . Así pues, no hay unicidad.

 0
 x = 2−y a] ¿Para qué valores de c es asintóticamente estable?
2. Sea y0 = −2y + cz . [3 puntos]
 0
z = 2x − y b] Resolver si c = −1 con x (0) = 0, y(0) = 0, z(0) = −4 .

 
−λ −1 0 
a] P(λ) = − 0 = λ3 + 2λ2 + cλ + 2c [ = (λ − 1)(λ + 1)(λ + 2) si c = −1 ].

−2 − λ c 
 2 −1 −λ
Por el signo de los coeficientes, no puede ser asintóticamente estable cuando c ≤ 0 .
Para ver lo que sucede cuando c > 0 acudimos a Routh-Hurwitz:
 
2 1 0
B= 2c c 2 → 2 , 0 , 0 → nunca es AE.
0 0 2c
√
[De hecho, con un poco de vista se observa que ∀c : P(λ) = (λ+2)(λ2 +c) → λ = −1, ± −c →
EnoA si c > 0 , I si c < 0 (y si c = 0 también I pues hay un único v.p. asociado a λ = 0 doble)].
b] Para resolverlo por matrices comenzamos hallando los vectores propios:
           
−1 −1 0 1 1 −1 0 1 2 −1 0 1
λ = 1:  0 −3 −1 v = 0 → −1,
 λ = −1: 0 −1 −1 v = 0 → 1 ,
 λ = −2: 0 0 −1 v = 0 → 2,
2 −1 −1 3 2 −1 1 −1 2 −1 2 0

Para ahorrarnos calcular P−1 y las integrales buscamos ( λ = 0 no autovalor):

1 − e−2t
    
a  0 = 2− b a=1 
1
  
1
 
1
 
1 d.i.
x p = b  → 0 = −2b − d → b = 2 . x=   t   − t   − 2t
c 1 −1 e + c 2 1 e + c 3 2 e +  2  → 2 − 2e−2t 
d

0 = 2a − b d = −4 3 −1 0 −4 −4

z = x 00 + 2x 0 − 4
Convirtiendo en ecuación: y = 2 − x 0 → → x 000 + 2x 00 − x 0 − 2 = −2 →
z0 = 2x + x 0 − 2
x = c1 et + c2 e−t + c3 e−2t + 1

→ x = 1 − e−2t , y = 2 − x 0 , z = −2y − 2y0 .
x (0) = 0, x 0 (0) = 2, x 00 (0) = −4
4(s+1)(s−1)

2
 sX = s − Y
 (s3 + 2s2 − s − 2)Y = 4s − 4s = s → Y= 4
s ( s +2) ↓
Laplace: sY = −2Y − Z → Z = −(s + 2)Y & Z = − 4s

sZ + 4 = 2X − Y 2X = 4 +(1 − 2s − s2 )Y ↑

0
z = −4 , Y = A B
s + s +2 = · · · = 2s − s+2 2 → y = 2 − 2e−2t , x = y−2 z = 1 − e−2t .
3. a] Dibujar el mapa de fases de (E) x 00 = 2x3 − 2x . b] ¿Para qué valores de b es periódica la
solución x (t) de (E) que cumple x (0) = 0 , x 0 (0) = b ? c] ¿Qué ecuación de primer orden
satisface la solución x (t) de (E) con x (0) = 0 , x 0 (0) = 1 ? [2.5 puntos]

a] Al ser exacta conviene dibujar la función potencial:

V ( x ) = − (2x3 − 2x )dx = x2 − 12 x4 ,
R

que tiene un mínimo en x = 0 (centro del mapa de fases),

y dos máximos en x = ±1 (sillas). Es V (0) = 0 , V (±1) = 21 .
Las órbitas vienen dadas por:
v2
√
4 2
2 + V ( x ) = C → v = ± x − 2x + 2C .

 Con las técnicas generales de mapas de fases:

√
x0 = v    
0 3
v = 2x − 2x
→ puntos críticos 00 , ±01 , con λ = ±i 2 , λ = ±2 .
El origen es centro del no lineal (simetrías o exactitud) y los
 
vectores propios de las sillas son, como en toda ecuación, λ1 .
La pendiente es vertical si v = 0 (en toda ecuación)
y es horizontal si x = 0, ±1 (rectas por los puntos).
Las separatrices son las órbitas que pasan por (±1, 0) :
√
0 = ± −1 + 2C → C = 12 → v = ±(1 − x2 ) , parábolas que cortan el eje v en ±1 .
b] A la vista de los cálculos anteriores y del mapa de fases está claro que es periódica si |b| < 1 .

c] La órbita asociada a esos datos es v = 1 − x2 = x 0 (la solución de esta autónoma es x = th t ).

4. Sea 3tx 00 +(2 − 6t) x 0 + 2x = 0 . a] Hallar una solución que no sea analítica en t = 0 . b] Hallar
4 términos del desarrollo de una solución no trivial que sea analítica en t = 0 . [2.5 puntos]

(2−6t) 0 2t
t2 x 00 + t 3 x + 3 x=0 → t = 0 es singular regular con λ(λ − 1)+ 23 λ = 0 → λ = 31 , 0 .
∞
a] Es no analítica x1 = t1/3 ∑ ck tk → que llevada a la ecuación nos da:
k =0
∞
3(k + 31 )(k − 32 )ck tk−2/3
+ 2(k + 13 )ck tk−2/3 − 6(k + 31 )ck tk+1/3 + 2ck tk+1/3
 
∑
k =0
∞ 
= ∑ (3k + 1)kck tk−2/3 − 6kck tk+1/3 = 0 → t−2/3 : 0c0 = 0 ; t1/3 : 4c1 = 0 ;

k =0

6( k −1)
tk−2/3 : ck = c
k(3k+1) k−1
→ c2 = c3 = · · · = 0 → x1 = t1/3 .
b] Conocida una solución tan sencilla, podemos abandonar Frobenius:
2)
1+2t+2t2 + 43 t3 +···
R
(2− 3t n n
t1/3 e
dt = t1/3 = −3(1 − t − 52 t2 − 16 t4 +· · · ) = −3 ∑ n!(21−t 3n) .
R R
x2 = t2/3 t4/3
dt
O bien, continuar con las series:
∞ ∞
x2 = ∑ bk tk → ∑ [(3k − 1)kbk tk−1 − 2(3k − 1)bk tk ] = 0 → t0 : b1 = −b0 ; t1 : b2 = 52 b1 = − 25 b0 ;
k =0 k =0

2(3k−4)
tk−1 : bk = k(3k−1) bk−1 → b3 = 12
5
b2 = − 61 b0 ; . . . → x2 = 1 − t − 52 t2 − 61 t3 +· · · .
Soluciones de septiembre 2008 de Ecuaciones Diferenciales I (grupos A y C)
y y2
1. Sea y0 = 1 + t − t2 . a] Hallar su solución general y la(s) solución(es), si existe(n), con:
i) y(1) = 1 , ii) y(1) = −1 . b] Dibujar aproximadamente sus curvas integrales. [2.5 puntos]

y Ct −1 2 Ct −1 2
a] z = t → tz0 = 1 − z2 → 2dz +z
= ln 11− 1+ z
= Ct2 → z = Ct
R
1− x 2 z = 2 ln t + C → 1− z 2 +1 → y = t Ct 2 +1 .
z=y−t 2 u =1/z
[O Riccati: y = t solución −→ z0 = − zt − zt2 −→ u0 = ut + t12 → u = Ct + t t−3 dt = Ct − 2t1
R

2t Ct +1 2
→ z= Ct2 −1
→ y = t + Ct2t
2 −1 = t Ct2 −1 , casi igual].

f y f y continuas en un entorno de (1, 1) y de (1, −1) ⇒ solución única para cada dato.
C −1 C −1
y (1) = C +1 = 1 imposible (es la y = t perdida). y(1) = C +1 = −1 → C = 0 → y = − t .
b] Isoclinas (como en toda homogénea) y = mt → K = 1 + m − m2 .
Rectas solución si K = m → m = ±1 . m = 0 → K = 1 .
√
m = 2 → K = −1 . 1.6
Horizontal: K = 0 → m = 1±2 5 ≈ –0.6
m = −2 → K = −3 .
t = 0 es curva integral vertical ( t = 0 → K = ∞ ).
2y2 t−2y
t2 +ty−y2
2y3 −ty2 −2t2 y+t3
y
y00 = − t2 + t3 + t2 t 2 = t4
= (y−t)(y+t4t)(2y−t)
y00 = 0 → rectas solución e y = 2t genuinos puntos de inflexión.

2. Sea x 000 + 3x 00 + 4x 0 + bx = 2 . a] Discutir su estabilidad según los valores de b .

b] Hallar una solución particular para todo valor de la constante b .
c] Hallar si b = 2 la solución con x (0) = 0 , x 0 (0) = 1 , x 00 (0) = 0 . [3 puntos]

a] P(λ) = λ3 + 3λ2 + 4λ + b [ = (λ + 1)(λ2 + 2λ + 2) si b = 2 ].

Por el signo de los coeficientes, no puede ser asintóticamente estable cuando b = 0 y es
inestable si b < 0 . Para ver lo que sucede cuando b > 0 acudimos a Routh-Hurwitz:
 
3 1 0
B= b 4 3 → 3 , 12 − b > 0 , b > 0 → AE si 0 < b < 12 ( e I si b > 12 ).
0 0 b
√
Para b = 0 : P(λ) = λ(λ3 + 3λ + 4) → λ = 0 (simple) y λ = −3±2i → EnoA 7

Para b = 12 : P(λ) = (λ3 + 3)(λ2 + 4) → λ = ±2i (simples) y λ = −3 → EnoA

b] Si b 6= 0 , λ = 0 no es autovalor y hay x p = A → x p = 2b .
Si b = 0 , λ = 0 es autovalor y la x p = At → x p = 2t .
c] La solución general para b = 2 es: x = c1 e−t + e−t (c2 cos t + c3 sen t) + 1 →

x 0 = e− [−c1 +(c3 − c2 )c −(c2 + c3 )s] d.i. c1 + c2 + 1 = 0 c =c =0

→ − c1 + c3 − c2 = 1 → 1 3 → x = 1 − e−t cos t .
x 00 = e− [c1 − 2c3 c + 2c2 s] c1 − 2c3 = 0
c2 = −1

2 s2 +3s+2 s +2 A Bs+C
Laplace: s3 X − s + 3s2 X − 3 + 4sX + 2X = s → X= s(s+1)(s2 +2s+2)
= s(s2 +2s+2)
= s + s2 +2s+2
↑
→ A(s2 + 2s + 2)+( Bs + C )s = s + 2 → X = 1s − (s+s1+)12 +1 → x = 1 − e−t L−1 s

s2 +1
 a] Hallar la expresión de sus órbitas y dibujar su mapa de fases.
x 0 = x + 2y − 3

3. Sea [S] . b] ¿Para qué valores de a tiende a −∞ cuando t → ∞ la y(t)
y0 = 4x − y − 3
de la solución de [S] con x (7) = 0 , y(7) = a ? [2.5 puntos]

a] El sistema (además de claramente lineal) es exacto:

Hy = x + 2y − 3
→ → H = y2 + xy − 2x2 − 3y + 3x = C .
Hx = y − 4x + 3
O con un poco de vista, (y − x )(y + 2x − 3) = C , órbitas.

x = 3 − 2y &
 
1
El punto crítico: ha de ser silla o centro.
9 − 9y = 0 → 1
     
1 2
4 −1
→ λ = ±3 silla, con v = 11 y −12 , respectivamente.
v es vertical si x = 3 − 2y y horizontal si y = 4x − 3 .
Para un sistema lineal no se necesita más para hacer un buen dibujo.
b] A la vista del mapa de fases, las soluciones (definidas ∀t por ser lineal con los coeficientes
continuos) tales que su y(t) → −∞ son las situadas a la izquierda de la separatriz estable:
t→∞
 
y = 3 − 2x recta de pendiente −2 por (1, 1) , o a partir de las órbitas .
Como esta separatriz corta x = 0 en y = 3 , los a pedidos son los a<3 .
[Para ver la utilidad del mapa de fases vamos a llegar al mismo resultado resolviendo el sistema:
y 0 + y +3 y00 +y0 0
x= 4 → 4 = y +4y+3 + 2y − 3 → y00 − 9y = −9 → y = c1 e3t + c2 e−3t + 1 ; de los datos:
y (7) = a c e21 + c e−21 + 1 = a
0 → 1 21 2 −21 → 6c1 e21 = 2a − 3 → y = ( 3a − 1)e3(t−7) + · · · ].
y (7) = − a − 3 3c1 e − 3c2 e = −a −3

4. Sea [e] 2t2 x 00 + t(1 + t2 ) x 0 − x = 0 . a] Hallar los 3 primeros términos no nulos del desarrollo
de una solución de [e] que no sea analítica en t = 0 . b] ¿Cuántas soluciones de [e] satisfacen:
i) x (0) = x 0 (0) = 3 ; ii) x (3) = x 0 (3) = 0 ? [2 puntos]

a] t = 0 es singular regular con λ(λ − 1)+ 21 λ + 12 = 0 → 2λ2 − λ − 1 = 0 → λ = 1, − 12 .

∞ ∞
Es analítica x1 = t ∑ ck tk y no lo es x2 = t−1/2 ∑ bk tk , que llevada a la ecuación nos da:
k =0 k =0
∞
∑ 2(k − 12 )(k − 32 )bk tk−1/2 + (k − 12 )bk tk−1/2 + (k − 21 )bk tk+3/2 − 6bk tk−1/2
 
k =0
∞ 
= ∑ (2k2 − 3k)bk tk−1/2 + (k − 21 )bk tk+3/2 = 0 →

k =0

t−1/2 : 0b0 = 0 , b0 indeterminado; t1/2 : −b1 = 0 ;

−5
tk−1/2 : k(2k − 3)bk + (k − 25 )bk−2 = 0 → bk = − 2k2k b
(2k−3) k−2
→
(−1)
b3 = b5 = · · · = 0 y además: b2 = − 4·1 b0 = 41 b0 , b4 = − 83·5 b2 = − 160
3
b0 , . . .

x2 = t−1/2 1 + 41 t2 − 3 4
 
Por tanto 160 t +··· .

b] El TEyU asegura solución única (la trivial x ≡ 0 ) para ii) y no dice nada para i)
por ser a(t) y b(t) discontinuas en t = 0 . Como la solución general de [e] es
x = c0 [t + · · · ] + b0 [t−1/2 + · · · ] , ninguna solución cumple i).
Ecuaciones Diferenciales I (C - piloto) Soluciones de Febrero09

dy (y−t)2
1. Sea dt = + 1 . a] Hallar sus rectas solución y resolverla como: i) homogénea, ii) Riccati.
t2
b] Discutir cuántas soluciones hay con: i) y(0) = 2, ii) y(1) = 2, iii) y(2) = 2, hallándolas si existen.

a] En una homogénea (de isoclinas y = mt ) las rectas solución salen de:

ƒ (m) = (m−1)2 +1 = m → m = 1, 2 → y=t e y = 2t .
[Con esto y el crecimiento, sería muy fácil hacer el dibujo de las soluciones].
dz dz dz z−2 y−2t
R R R
i) tz 0 +z = z 2 −2z+2 , (z−1)(z−2)
= z−2
− z−1 = ln z−1 = ln t +C , = Ct
y−t ↓

2 2 2
ii)  = y−t → 0 = t 2 +1−1 = t 2 → 1

t
= 1t −C → y = t + 1−Ct → y = 2t−Ct
1−Ct
.

b] Como ƒ y ƒy = 2(y−t)/ t 2 son continuas en un entorno de cada uno de los dos puntos, el TEyU
asegura solución única para ii) y iii). Son, respectivamente, las rectas y = 2t e y = t .
4−4C
Ojo con las soluciones mentirosas; imponiendo iii): 2 = 1−2C → 2 = 4 ; y = t se ha perdido .
 

dt 2
El teorema no dice nada para i), pero aplicándolo a dy
= (y−t)t 2 +t 2 , deducimos que por (0, 2)
pasa la única curva integral vertical t ≡ 0 ⇒ no existe solución y(t) cumpliendo y(0) = 2 .
[La solución general nos decía que todas satisfacen y(0) = 0 , pero se podía haber perdido alguna].

¨ 0 = −3+4y+cz
a] Para c = −2 , hallar la solución con (0) = 0 , y(0) = 1 , z(0) = 2 .
2. Sea y 0 = −+y− z .
z 0 = −2y b] Discutir la estabilidad del sistema según los valores de c ∈ R.

¨ 0 = −3+4y−2z 
−3−λ

4 −2
y0 = −+y− z = −λ3 −2λ2 −λ = −λ(λ+1)2 .

a] .. Mediante matrices:  −1 1− λ −1
z0 = −2y

 1 −2 −λ 
−3 4 −2 2 −2 4 −2 2 1
         
λ=0 : −1 1 −1v = 0 →  1  ; λ = −1 : −1 2 −1v = 0 → 1 y  0  →
1 −2 0 −1 1 −2 −1 0 −1

2c1 +2c2 +c3 = 0

 
2
 
2
  
1
 0
x= c1  1 + c2 1e−t + c3  0 e−t → c1 + c2 =1 → c1 = 0, c2 = 1, c3 = −2 , x = e 
 −t 
.
−1 0 −1 d.i. −c1 − c3 =2 2e−t
 
0 2 2 1 −1 2 −1 1 0 0
     
Más largo: x = PeJt P−1 1 con P = 1 1 0  → P−1 = 1 −1 1 , eJt =
0 e−t 0 .

2 −1 0 −1 1 2 0 0 0 e−t

z 00 +2y 0 = −3z 0 −2y−2z , z 00 +3z 0 +2z+2(y 0 +y) = 0


O bien, de la 3 ecuación:  = z 0 +2y → →
y 0 = −z 0 −y−z z 0 +z + y 0 +y = 0
z(0)=2,z 0 (0)=−2 z 0 +z=0 y(0)=1
z 00 +z 0 = 0 → z = c1 +c2 e−t −→ z = 2e−t −→ y 0 +y = 0 → y = Ce−t −→ y = e−t
→  = z 0 +2y = −2e−t +2e−t →  = 0
¨ sX = −3X+4Y −2Z (s+1)X = 4−2(s+1)Z &
2
O Laplace: sY −1 = −X+Y −Z (s+1)X = 4−2s+(s2 −s−2)Z → 2s = (s2 +s)Z → Z = s+2
1 s
sZ−2 = X−2Y → Y = 2
X− 2 Z+1 % s 1
→ X = 0 → Y = − s+1 +1 = s+1 . Como arriba.
 
−3−λ 4 c 
= λ3 +2λ2 −(c+1)λ−(c+2) . Del signo de los coeficientes:

b] −  −1 1− λ −1
 1 −2 −λ 
El tercer coeficiente es negativo si c > −1 y el cuarto si c > −2 ⇒ si c > −2 alguno es
negativo y el sistema es I. Si c = −2 no es AE y si c < −2 aún no lo sabemos.
2 1 0
 

Routh-Hurwitz: B = 
−c−2 −c−1 2 

→ si 2 > 0 , −c > 0 y −c−2 > 0 es AE →
0 0 −c−2
AE ⇔ c < −2 e I si c > −2 (ya lo sabíamos).

Si c = −2 (el caso de a]) es EnoA [ λ = −1 doble y λ = 0 simple].

3. Sea t00 + (2t 2 −1)0 − 4αt = 0 .
a] Precisar para qué valores de α hay soluciones que son polinomios que se anulan en t = 0 .
b] Para α = 1 , hallar una solución analítica en t = 0 y determinar si todas las soluciones lo son.
c] Para α = 0 , hallar la solución general sin utilizar series.

a] t 2 00 + t(2t 2 −1)0 −4αt 2  = 0 → t = 0 singular regular con r = 2, 0 . Se anula en t = 0 :

∞ ∞
1 = ck t k+2 → (k +2)(k +1)ck t k+1 + 2(k +2)ck t k+3 − (k +2)ck t k+1 − 4αck t k+3 = 0
X X  
k=0 k=0

t 1 : 2c0 −2c0 = 0 , c0 indeterminado; t 2 : 6c1 −3c1 = 0 , c1 = 0 ;

t 3 : 12c2 +4c0 −4c2 −4αc0 = 0 , c2 = α−1
2
c0 ;
4α−2k
t k+1 : (k +2)kck −(4α−2k)ck−2 = 0 , ck = (k+2)k ck−2 → c3 = c5 = c7 = 0 .
α−k (α−1)(α−2) (α−1)(α−2)(α−3)
c2k = (k+1)k c2k−2 → c4 = α−2
6
c2 = 12
c0 , c6 = α−3 c =
12 4 144
c0 , . . .
6 N , ninguno de los c2k se anula, pero si α = 1, 2, . . . , la serie pasa a ser un polinomio
Si α ∈
P2n de grado 2n , pues el c2n = 0 y también lo son todos los siguientes.
h i
En particular: P1 = t 2 , P2 = t 2 + 21 t 4 , P3 = t 2 + t 4 + 16 t 6 , . . . .

b] 1 = t 2 es solución analítica de t00 + (2t 2 −1)0 − 4t = 0 en t = 0 . Otra solución es:

Z R 1 −2t
2
e−t
Z Zh
2 e t t3 4 t6
i
1 1 t
2 = t = t 2 = t 2 − + − + · · · = − 12 + t4 − 24 + · · · − t 2 log t .
t4 t3 t 3 t 2 6

No todas las soluciones son analíticas en t = 0 . Se podría ver siguiendo con Frobenius:
∞ ∞ ∞
2 = bk t k + d t 2 ln t , 02 = kbk t k−1 + 2d t ln t + d t , 00 = k(k −1)bk t k−2 + 2d ln t + 3d
X X X

k=0 k=1
2 k=2
∞ ∞ ∞
k(k −1)bk t k−1 + 2kbk t k+1 −kbk t k−1 − 4bk t k+1 + d 2t + 2t 3 = 0 →
X X   X  
→
k=2 k=1 k=0

t 0 : b1 = 0 ; t 1 : 2b2 −2b2 −4b0 +2d = 0 → d = 2b0 6= 0 (como antes; y b2 indeterminado).


1 2 2 Las series de Frobenius son las de
c] 0 = y → y 0 = t −2t y → y = Cte−t  = K + Ce−t


→ . 2
1 = 1−e−t y 2 = 1 (analíticas) .


4. Sea (E) 00 = 1 − 2 − (0 )2 . a] Resolver la ecuación de las órbitas de (E).

b] Dibujar el mapa de fases y decidir si es periódica la solución de (E) con (0) = 0 (0) = 1 .

0 = 

d 1−2 z= 2 dz
a] , = −  −→ = −2z + 2−22 lineal de coeficientes constantes.
 0 = 1−2 − 2 d  d

Más corto que integrar por partes: zp = A2 +B+C → z = Ce−2 − 2 ++ 12 =  2 .

b] Puntos críticos:  = 0 →  = ±1 .
(±1,0)
‚ Œ ‚ Œ
0 1 0 1
−→ → λ2 ±2 = 0 →
−2 −2y ∓2 0

−1
 p h   i
0
silla con λ = ± 2 v = λ1 , como en toda ecuación .
 
1
0
centro de la AL. Como g es par en  se conserva.

Campo horizontal sobre la circunferencia 2 + 2 = 1

[y vertical sobre el eje  como en toda ecuación].
La separatriz es la órbita que pasa por (−1, 0) →
Ce2 −1−1+ 21 = 0 →  2 = 3 −2−2
2
e − 2 ++ 12
[El 2o miembro se hace 0 para un  > 1 ( → −∞ ) , y
si  muy negativo se parece a un múltiplo de e−2 ].
Hay una órbita sencilla (circunferencia) para C = 0 : 

 



2 p , v(±1, ) = − 2 . v(0, ) =
 2 + − 12 = 34 corta  = 0 en  = 1+2 3 ≈ 1.4,–0.4 .
  1− 2
q
La órbita por (1, 1) no es cerrada, pues la separatriz corta  = 1 en  = 3e 2 +1 < 1 (también
−4

lo dice el campo). Por tanto, una solución de la que es proyección no es periódica.

Ecuaciones Diferenciales I (C - piloto) Soluciones de septiembre09

dy 2+t−y a] Hallar su solución general (en forma implícita) calculando

1. Sea dt = y−t .
e +1 un factor integrante g(t) y por otro método diferente.
b] ¿Cuántas soluciones hay con y(0) = 0 ? Dar su forma explícita si existen.
[Ayuda: ¿Qué forma tienen todas sus isoclinas? ¿Cuál es la asociada a la pendiente K = 1 ?].

a] g(t)(y−t−2) + g(t)(ey−t +1)y 0 = 0 exacta si g(t) = g0 (t)(ey−t +1) − g(t)ey−t , g0 (t) = g(t) , g(t) = et
Ut = et (y−t −2) → U = et (y−t −1) + p(y)
→ et (y−t −2) + (ey +et )y 0 = 0 → → ey + (y−t −1) et = C .
Uy = ey + et → U = ey +yet + q(t)
↑
= − e e+z−1 ez +1
2−z z R
O bien, z = y−t → z 0 = ez +1
−1 z +1 → ez +z−1
dz = ln(ez +z−1) = C − t → ez +z−1 = Ce−t .

b] ƒ y ƒy continuas en todo R2 ( ƒ cociente de funciones que lo son y su denominador no se anula)

⇒ existe solución única para cualquier dato inicial, y(0) = 0 en particular.
Las isoclinas son las rectas paralelas y = t +b → K = e2−b
b +1 .

La de K = 1 (es decir, aquella en que la pendiente de los

segmentos es igual que la de la isoclina) cumple eb = 1−b ,
lo que sólo sucede si b = 0 . Por tanto, y = t es solución,
que es la única buscada, pues satisface y(0) = 0 .
[También se puede localizar esta solución observando que z = 0 ( y = t ) es solución constante
de la autónoma en z de arriba, o tanteando en la ecuación inicial, o viendo que y = t cumple
ey + (y−t −1) et = 0 , expresión que sale de la solución general al imponer el dato].

2. Sea 000 + 200 + 0 = 4et + 2t . a] Para  = 1 , hallar su solución con (0) = −2 , 0 (0) = 00 (0) = 0 .
b] Discutir la estabilidad de la ecuación según los valores de  .

a] λ3 +2λ2 +λ = λ(λ+1)2 → solución general de la homogénea:  = c1 + c2 e−t + c3 te−t .

p = Aet + Bt 2 +Ct → 4Aet + (4B+C) + 2Bt = 4et + 2t → A = 1 , B = 1 , C = −4 →

Solución general de la no homogénea:  = c1 + c2 e−t + c3 te−t + et + t 2 −4t →
datos
0 = −c2 e−t + c3 (1−t)e−t + et + 2t −4 , 00 = c2 e−t + c3 (t −2)e−t + et + 2 −→
 c1 + c2 + 1 = −2
c1 = 0
−c2 + c3 − 3 = 0 o % →  = et − 3e−t + t 2 −4t .
c3 = 0 , c2 = −3
c2 − 2c3 + 3 = 0
4 4s2 +2s+2 2(2s−1)(s+1) 2(2s−1)
O Laplace: (s3+2s2+s)X + 2(s2+2s+1) = s−1
+ s22 = s2 (s−1)
= s2 (s−1)
, X = − 2s + s3 (s+1)(s−1) .

2(2s−1) A B C D E As2 (s−1)(s+1)+Bs(s−1)(s+1)+C(s−1)(s+1)+Ds3 (s+1)+Es3 (s−1)

s3 (s−1)(s+1)
= s
+ s2
+ s3
+ s−1
+ s+1
= s3 (s−1)(s+1)

s = 1 → D = 1 , s = 1 → E = −3 , s = 0 → C = 2 , s1 → B = −4 , s4 → A+D+E = 0 → A = 2
X = − s42 + 2
s3
+ 1
s−1
− 3
s+1
→  = −4t + t 2 + et − 3te−t .

b] λ(λ2 +2λ+) . Como es un λ = 0 , a simple vista deducimos que no puede ser AE.
p
Y los otros dos autovalores son fácilmente calculables: λ± = −1 ± 1− . De aquí:
Si  = 0 , hay λ = 0 doble, lo que, al tratarse de una ecuación, basta para asegurar que es I.
Si  < 0 , es λ+ > 0 (y λ− < 0 ), con lo que la ecuación también es I.
Si  > 0 , los λ± o son reales negativos o complejos con Reλ < 0 . Por tanto, es EnoA.

[El signo de los coeficientes o RH (poco útiles conociendo los λ ), sólo dan la inestabilidad para  < 0 ;
de aplicar RH al polinomio de grado 2 sí deduciríamos (considerando el λ = 0 ) la EnoA para  > 0 ].
3. Sea t00 − 20 + 4et  = 0 . a] Hallar los 3 primeros términos no nulos del desarrollo en serie de
una solución que se anule en t = 0 . b] ¿Están acotadas en t = 0 todas las soluciones?
∞
a] t 2 00 + t (−2) 0 + 4tet  = 0 → r(r−1)−2r = 0 → r1 = 3 , r2 = 0 → 1 = ck t k+3 se anula en t = 0 .
X

k=0

t 6c0 t + 12c1 t 2 + 20c2 t 3 + 30c2 t 4 + · · · − 2 3c0 t 2 + 4c1 t 3 + 5c2 t 4 + 6c3 t 5 + · · ·

+ 4 + 4t + 2t 2 + · · · c0 t 3 + c1 t 4 + c2 t 5 + · · · = 0 →

t 2 : 0c0 = 0 → c0 indeterminado;
t 3 : 12c1 − 8c1 + 4c0 = 0 → c1 = −c0 ;
t 4 : 20c2 − 10c2 + 4c1 + 4c0 = 0 → c2 = 0 ;
1
t 5 : 30c3 − 12c3 + 4c2 + 4c1 + 2c0 = 0 → c3 = c .
9 0

Ya tenemos 3 no nulos: 1 = t 3 − t 4 + 19 t 6 + · · · .

∞
b] Como tanto 1 como 2 = bk t k + d t 3 (1 − t + · · · ) ln t están acotadas en t = 0 t 3 ln t → 0 ,
X

k=0 t→0
todas las soluciones de la ecuación están acotadas en t = 0 .

a] Dibujar el mapa de fases de [S].

 0
 = y − 
4. Sea [S] .
y 0 = 2y−3 b] Hallar la solución de [S] que cumple (2) = 0 , y(2) = 1 .

=0 y=1     ‚
y−1 
Œ
0 2/ 3
a] ↓ ↓ → 0
y 1
puntos críticos. M = .
y=0  = 2/ 3 −3 2

−1 0
       
0 1 0
0
→ −3 2
→ λ = −1 , 1
; λ=2 , 1
. Silla.
0 2/ 3
   
2/ 3
1
→ −3 2 → λ = 1±  , foco inestable.

Campo horizontal sobre la recta y = 32  .

Vertical sobre  = 1 y sobre el eje y [órbita].
   
v(, 0) = −  1
3
, v(2/ 3, y) = 23 (y−1) 1
3
.
2 1−
   
 −
v(, ) = −
= − .
1

Como para  < 0 la primera componente es mayor

que la segunda, la separatriz estable se deforma ^ .
No sabemos resolver la ecuación de las órbitas.

b] La órbita que pasa por el punto (0, 1) es la  = 0 (única órbita sencilla y calculable).

Sobre ella es 0 = 0 (claro, es (t) ≡ 0 para todo t ), y además:

‚ Œ
0
y 0 = 2y → y = Ke2t −→ y = e2t−4 . La solución pedida es, pues, x= .
y(2)=1 e2t−4