Hoja 4 1

Problemas de E.D.O. (2o de Matemáticas)

1. Dado el problema (
y 0 = y 2/3
y(0) = 0
hállense al menos tres soluciones. Indicación: Combı́nense las dos soluciones que se
pueden obtener de forma sencilla.

2. Calcúlense todos los valores α ∈ [0, ∞) para los que en el problema

(
y 0 = |y|α
y(0) = 0

haya existencia y unicidad (para α = 0 escrı́base |y|α = 1).

1
¡3. Decı́dase
¢ si cada una de las siguientes implicaciones es cierta para funciones y, z ∈
C [0, ∞) , dando un contraejemplo o una demostración:
a) y ≥ z ⇒ y 0 ≥ z 0 .
b) y 0 ≥ z 0 ⇒ y ≥ z.
c) y(0) = z(0), y 0 ≥ z 0 ⇒ y ≥ z.

4. Estúdiese la existencia y unicidad para el problema

(
y 0 = |y| + x
y(0) = y0

y hállense explı́citamente las soluciones cuando y0 = 1 y cuando y0 = −1, indicando a

qué espacio C n (R) pertenecen.

5. Estúdiese si para cada x0 , y0 la solución de

 2
y 0 = xy + y
x2 + y 2 + 2

y(x0 ) = y0

se puede definir en toda la recta real.

6. Considérese (
y0 = y4 + r
y(0) = 0
donde r es una constante positiva. Hállese un entorno de cero (intentando que sea lo
mayor posible) en el que se pueda asegurar existencia y unicidad. Pruébese que si la
solución existe en [0, r−3/4 ] entonces y(r−3/4 ) ≥ r1/4 y utilı́cese este hecho junto con
Hoja 4 2

y 0 ≥ y 4 para encontrar otro entorno en el que se pueda asegurar que no existe solución
regular.

7. Sea y la solución de (
y 0 = y + sen(xy)
y(0) = 1
y sea z la solución de esta ecuación aproximando sen(xy) por xy.
a) Hállese una cota superior para máx |z(x) − y(x)| cuando x ∈ [0, 00 1].
b) Usando el apartado anterior, cálculese una aproximación para y(00 1).
c) ¿Qué cota superior se podrı́a dar para máx |z(x) − y(x)| si x ∈ [−00 1, 0]?

8. Sean los problemas

( y 2
( z
y0 = + e−y z0 =
1 + x2 1 + x2
y(0) = 10 y(0) = 10

Demuéstrese que 0 ≤ y(x) − z(x) ≤ e−100 (ex − 1) para x ∈ [0, 1].

9. Hállense una fórmula para yn cuando se aplica el método de Euler a

(
y0 = y + 1
y(0) = 0

escogiendo xn = nh con n = 0, 1, 2, . . . Si se toma h = 00 001, ¿qué error se obtiene al

aproximar y(1)?

10. Interprétese la ecuación de movimiento correspondiente al enunciado del proble-

ma 6 de la Hoja 1 como la aplicación del método de Euler a cierto problema.

11. En Mecánica, para θ0 moderadamente pequeño se suele aproximar la solución de

la ecuación del péndulo (
y 00 = − sen y
y(0) = θ0 , y 0 (0) = 0
por θ0 cos x, que en realidad resuelve y 00 = −y. Acótese el error en esta aproximación
cuando θ0 = 00 2 y x ∈ [0, 00 1]. Indicación: Escrı́base primero el problema como un
sistema de primer orden.

12. Considérese la solución del problema

(
y 0 = f (y)
y(0) = λ

como función de x y λ, esto es, y = y(x, λ). Suponiendo la regularidad que sea precisa,
escrı́base una ecuación de primer orden cuya solución sea w = ∂y/∂λ.
Hoja 4 1

Problemas de E.D.O. Curso 2008-2009

1. Para cada una de las sucesiones {fn }∞ n=1 siguientes, determı́nese el lı́mite pun-
tual de la sucesión (si existe) en el conjunto o conjuntos indicados, e indı́quese si la
convergencia es uniforme.

a) fn (x) = exp(−n x2 ) , sobre [−1, 1] .

b) fn (x) = x1/n , sobre [0, 1] .
xn
c) fn (x) = en [0, 1 − ε], en [1 − ε, 1 + ε], y en [1 + ε, ∞).
1 + xn
(
0 si x ≤ n
d ) fn (x) = en cada [a, b] y en R .
x − n si x ≥ n
nx
e) fn (x) = en [−1, 1] y en [1, ∞).
1 + n2 x2
f ) fn (x) = x−n ex en (1, ∞).

2
2. Sea la sucesión de funciones {fn }∞ 2
n=1 en [0, 1], dada por fn (x) = n x e
−nx

a) Estúdiese la convergencia puntual e uniforme de {fn }∞ n=1 .

b) Compruébese que a pesar de que converge
R1 puntualmente a una función integrable
y que cada fn es integrable, se tiene lı́m 0 fn = ∞.

∞
3. Encuéntrese una sucesión de funciones Rcontinuas
R ∞{fn }n=1 que converja uniforme-
∞
mente a f en [0, ∞) para las que existan lı́m 0 fn y 0 f pero no coincidan.

4. Sea fn (x) = cos2n (πx) en R.

a) Estúdiese a qué función converge puntualmente la sucesión {fn }∞
n=1 y si la con-
vergencia es uniforme.
b) Descrı́base la función lı́mk→∞ lı́mn→∞ fn (k!x).

5. Sean {fn }∞ ∞ 2
n=1 y {gn }n=1 dos sucesiones de funciones dadas por fn (x) = x + 1/n
−1
y gn (x) = (nx) .
a) Demuéstrese que ambas convergen uniformemente en [1, ∞) y sin embargo la
sucesión de término general fn gn no lo hace.
b) Demuéstrese que a pesar de que {fn }∞ n=1 converge uniformemente en R a una
función f , {fn }n=1 no converge uniformemente en R a f 2 .
2 ∞

6. Sea la sucesión de término general fn (x) = x/(1+n x2 ). Compruébese que converge

uniformemente a cierta f en R y que se se verifica lı́mn→∞ fn0 (x) = f 0 (x) para cualquier
x 6= 0 pero no para x = 0.

7. Encuéntrese una sucesión de funciones derivables en (−1, 1) que converja uni-

formemente a f (x) = |x|.
Hoja 4 2

8. Dado el problema (
y 0 = y 2/3
y(0) = 0
hállense al menos tres soluciones. Indicación: Combı́nense las dos soluciones que se
pueden obtener de forma sencilla.

9. Calcúlense todos los valores α ∈ [0, ∞) para los que en el problema

(
y 0 = |y|α
y(0) = 0

haya existencia y unicidad (para α = 0 escrı́base |y|α = 1).

10. Decı́dase
¡ ¢ si cada una de las siguientes implicaciones es cierta para funciones
1
y, z ∈ C [0, ∞) , dando un contraejemplo o una demostración:
a) y ≥ z ⇒ y 0 ≥ z 0 .
b) y 0 ≥ z 0 ⇒ y ≥ z.
c) y(0) = z(0), y 0 ≥ z 0 ⇒ y ≥ z.

11. Estúdiese la existencia y unicidad para el problema

(
y 0 = |y| + x
y(0) = y0

y hállense explı́citamente las soluciones cuando y0 = 1 y cuando y0 = −1, indicando a

qué espacio C n (R) pertenecen.

12. Estúdiese si para cada x0 , y0 la solución de

 2
y 0 = xy + y
x2 + y 2 + 2

y(x0 ) = y0

se puede definir en toda la recta real.

13. Considérese (
y0 = y4 + r
y(0) = 0
donde r es una constante positiva. Hállese un entorno de cero (intentando que sea lo
mayor posible) en el que se pueda asegurar existencia y unicidad. Pruébese que si la
solución existe en [0, r−3/4 ] entonces y(r−3/4 ) ≥ r1/4 y utilı́cese este hecho junto con
y 0 ≥ y 4 para encontrar otro entorno en el que se pueda asegurar que no existe solución
regular.
Hoja 4 3

14. Sea y la solución de (

y 0 = y + sen(xy)
y(0) = 1
y sea z la solución de esta ecuación aproximando sen(xy) por xy.
a) Hállese una cota superior para máx |z(x) − y(x)| cuando x ∈ [0, 00 1].
b) Usando el apartado anterior, cálculese una aproximación para y(00 1).
c) ¿Qué cota superior se podrı́a dar para máx |z(x) − y(x)| si x ∈ [−00 1, 0]?

15. Sean los problemas

( y 2
( z
y0 = + e−y z0 =
1 + x2 1 + x2
y(0) = 10 z(0) = 10

Demuéstrese que 0 ≤ y(x) − z(x) ≤ e−100 (ex − 1) para x ∈ [0, 1].

*16. En este problema se probará que π 6∈ RQ considerando la sucesión de funciones

1
fn (x) = a2n xn (1 − x)n /n! y las integrales In = 0 fn (x) sen(πx) dx.
a) Empleando la convergencia uniforme, dedúzcase que para cualquier a > 0 se
cumple lı́m In = 0.
b) Pruébese que todas las derivadas de a−2n fn (x) en x = 0 y en x = 1 son números
enteros.
c) Suponiendo π = a/b, a, b ∈ Z+ , integrando por partes y empleando el apartado
anterior, demuéstrese que πIn ∈ Z.
d) Empléese que una sucesión de enteros (estrictamente) positivos no puede tener
lı́mite nulo, para llegar a una contradicción.
Hoja 4 1

Problemas de E.D.O. Curso 2009-2010

1. Para cada una de las sucesiones {fn }∞ n=1 siguientes, determı́nese el lı́mite pun-
tual de la sucesión (si existe) en el conjunto o conjuntos indicados, e indı́quese si la
convergencia es uniforme.

a) fn (x) = exp(−n x2 ) , sobre [−1, 1] .

b) fn (x) = x1/n , sobre [0, 1] .
xn
c) fn (x) = en [0, 1 − ε], en [1 − ε, 1 + ε], y en [1 + ε, ∞).
1 + xn
(
0 si x ≤ n
d ) fn (x) = en cada [a, b] y en R .
x − n si x ≥ n
nx
e) fn (x) = en [−1, 1] y en [1, ∞).
1 + n2 x2
f ) fn (x) = x−n ex en (1, ∞).

2
2. Sea la sucesión de funciones {fn }∞ 2
n=1 en [0, 1], dada por fn (x) = n x e
−nx

a) Estúdiese la convergencia puntual e uniforme de {fn }∞ n=1 .

b) Compruébese que a pesar de que converge
R1 puntualmente a una función integrable
y que cada fn es integrable, se tiene lı́m 0 fn = ∞.

∞
3. Encuéntrese una sucesión de funciones Rcontinuas
R ∞{fn }n=1 que converja uniforme-
∞
mente a f en [0, ∞) para las que existan lı́m 0 fn y 0 f pero no coincidan.

4. Sea fn (x) = cos2n (πx) en R.

a) Estúdiese a qué función converge puntualmente la sucesión {fn }∞
n=1 y si la con-
vergencia es uniforme.
b) Descrı́base la función lı́mk→∞ lı́mn→∞ fn (k!x).

5. Sean {fn }∞ ∞ 2
n=1 y {gn }n=1 dos sucesiones de funciones dadas por fn (x) = x + 1/n
−1
y gn (x) = (nx) .
a) Demuéstrese que ambas convergen uniformemente en [1, ∞) y sin embargo la
sucesión de término general fn gn no lo hace.
b) Demuéstrese que a pesar de que {fn }∞ n=1 converge uniformemente en R a una
función f , {fn }n=1 no converge uniformemente en R a f 2 .
2 ∞

6. Sea la sucesión de término general fn (x) = x/(1+n x2 ). Compruébese que converge

uniformemente a cierta f en R y que se se verifica lı́mn→∞ fn0 (x) = f 0 (x) para cualquier
x 6= 0 pero no para x = 0.

7. Encuéntrese una sucesión de funciones derivables en (−1, 1) que converja uni-

formemente a f (x) = |x|.
Hoja 4 2

8. Dado el problema (
y 0 = y 2/3
y(0) = 0
hállense al menos tres soluciones. Indicación: Combı́nense las dos soluciones que se
pueden obtener de forma sencilla.

9. Calcúlense todos los valores α ∈ [0, ∞) para los que en el problema

(
y 0 = |y|α
y(0) = 0

haya existencia y unicidad (para α = 0 escrı́base |y|α = 1).

10. Decı́dase
¡ ¢ si cada una de las siguientes implicaciones es cierta para funciones
1
y, z ∈ C [0, ∞) , dando un contraejemplo o una demostración:
a) y ≥ z ⇒ y 0 ≥ z 0 .
b) y 0 ≥ z 0 ⇒ y ≥ z.
c) y(0) = z(0), y 0 ≥ z 0 ⇒ y ≥ z.

11. Estúdiese la existencia y unicidad para el problema

(
y 0 = |y| + x
y(0) = y0

y hállense explı́citamente las soluciones cuando y0 = 1 y cuando y0 = −1, indicando a

qué espacio C n (R) pertenecen.

12. Estúdiese si para cada x0 , y0 la solución de

 2
y 0 = xy + y
x2 + y 2 + 2

y(x0 ) = y0

se puede definir en toda la recta real.

13. Considérese (
y0 = y4 + r
y(0) = 0
donde r es una constante positiva. Hállese un entorno de cero (intentando que sea lo
mayor posible) en el que se pueda asegurar existencia y unicidad. Pruébese que si la
solución existe en [0, r−3/4 ] entonces y(r−3/4 ) ≥ r1/4 y utilı́cese este hecho junto con
y 0 ≥ y 4 para encontrar otro entorno en el que se pueda asegurar que no existe solución
regular.
Hoja 4 3

14. Sea y la solución de (

y 0 = y + sen(xy)
y(0) = 1
y sea z la solución de esta ecuación aproximando sen(xy) por xy.
a) Hállese una cota superior para máx |z(x) − y(x)| cuando x ∈ [0, 00 1].
b) Usando el apartado anterior, cálculese una aproximación para y(00 1).
c) ¿Qué cota superior se podrı́a dar para máx |z(x) − y(x)| si x ∈ [−00 1, 0]?

15. Sean los problemas

( y 2
( z
y0 = + e−y z0 =
1 + x2 1 + x2
y(0) = 10 z(0) = 10

Demuéstrese que 0 ≤ y(x) − z(x) ≤ e−100 (ex − 1) para x ∈ [0, 1].

*16. En este problema se probará que π 6∈ RQ considerando la sucesión de funciones

1
fn (x) = a2n xn (1 − x)n /n! y las integrales In = 0 fn (x) sen(πx) dx.
a) Empleando la convergencia uniforme, dedúzcase que para cualquier a > 0 se
cumple lı́m In = 0.
b) Pruébese que todas las derivadas de a−2n fn (x) en x = 0 y en x = 1 son números
enteros.
c) Suponiendo π = a/b, a, b ∈ Z+ , integrando por partes y empleando el apartado
anterior, demuéstrese que πIn ∈ Z.
d) Empléese que una sucesión de enteros (estrictamente) positivos no puede tener
lı́mite nulo, para llegar a una contradicción.
Hoja 5 1

Problemas de E.D.O. (2o de Matemáticas)

1. Dados los sistemas lineales

( (
x0 = 2x + 3y − 7 x0 = −x − 2y + 3
y 0 = −x − 2y + 4 y 0 = 2x − y − 6

a) Determı́nese la naturaleza de sus puntos crı́ticos y sus propiedades de estabilidad.

b) En el caso en el que se obtiene un punto de silla, ¿cómo se podrı́an determinar las
direcciones de sus ejes?

2. En cada uno de estos sistemas determı́nese la naturaleza del punto crı́tico (0, 0) y
sus propiedades de estabilidad:
( (
x0 = x + y − 2 sen(xy) x0 = − sen(x + y) + 1 − e−x2
a) b)
y 0 = −2 x + y + 3 y 2 y 0 = −2 x + 4 y + y sen(x + y)

3. Estúdiense los puntos crı́ticos del siguiente sistema autónomo y esbócese una
posibilidad coherente con estos datos para las trayectorias del sistema en el primer
cuadrante. (
x0 = x (60 − 4 x − 3 y)
y 0 = y (42 − 3 x − 2 y)

4. Considérese el sistema
(
x0 = y + x (x2 + y 2 )
y 0 = −x + y (x2 + y 2 )

a) Estúdiese que tipo de punto crı́tico es (0, 0) en el sistema linealizado.

b) Resuélvase el sistema no lineal empleando coordenadas polares y decı́dase la esta-
bilidad de dicho punto crı́tico.

5. Se llama ecuación de van der Pol al problema

(
x0 = y
y 0 = −x − µy(x2 − 1)

donde µ 6= 2 es una constante. Estúdiese la naturaleza del punto crı́tico (0, 0) en función
de µ. Nota: en el caso en el que la linealización no es válida, se puede resolver el sistema.

6. Transfórmese la ecuación del péndulo x00 + sen x = 0 en un sistema autónomo con

el procedimiento habitual (escribiendo y = x0 ).
a) Calcúlense todos los puntos crı́ticos.
b) Hállese la ecuación de las trayectorias
Hoja 5 2

c) Decı́dase la estabilidad y carácter de todos los puntos crı́ticos.

d) ¿Qué función de Liapunov se podrı́a emplear para probar la estabilidad en (0, 0)?

7. Hállense los puntos crı́ticos de los siguientes sistemas y estúdiese su estabilidad.

( (
x0 = −y − x3 x0 = ex+y + y
a) b)
y0 = x − y3 y0 = y − x y

8. Pruébese que el sistema:

(
x0 = −y + x(1 − x2 − y 2 )
y 0 = x + y(1 − x2 − y 2 )

tiene (al menos) una solución periódica rodeando al origen. Estúdiese la estabilidad de
ese punto. Indicación: Utilı́cense coordenadas polares.

9. Determı́nese una función de Liapunov para

(
x0 = −2y − x3
y 0 = x/2 − 4y 3

10. Hállese una función de Liapunov de la forma α(2x + y)2 + β(x + y)2 que pruebe
la estabilidad en el origen del sistema autónomo:
(
x0 = −3x − 2y − (2x + y)3 + (x + y)3
y 0 = 5x + 3y − 2(x + y)3 + (2x + y)3

11. Estúdiese si los siguientes sistemas autónomos tienen alguna solución periódica.
( ( 2 2
x0 = y x0 = y + x3 (2 − ex +y )
a) b) 2 2
y 0 = −y − x2 y + x5 y 0 = −x + y 3 (2 − ex +y )

*12. Pruébese que si (0, 0) es asintóticamente estable para x0 = F (x, y), y 0 = G(x, y);
entonces es inestable para x0 = −F (x, y), y 0 = −G(x, y). Utilı́cese este hecho para
enuciar un teorema como el de Liapunov cuya conclusión sea la inestabilidad. Indicación:
Considérese la inversión temporal t 7→ −t.

*13. Escribiendo la ecuación x00 + x = ²x2 , x(0) = 1, x0 (0) = 0 (con ² pequeño) como
un sistema autónomo, pruébese que tiene una solución periódica. Ésta deberı́a parecerse
a cos t. El método de Lindstedt permite desarrollar la solución como
5 2
x(t) = cos(ωt) + ²x1 (ωt) + ²2 x2 (ωt) + . . . con ω =1− ² + ...
12
donde xj son funciones que se calculan de forma que al sustituir x(t) en la ecuación, los
términos de cualquier grado en ² desaparezcan. Confiando en este argumento hasta grado
dos, inferir que el centro de la oscilación de la solución es (3 − 2²)²/6 aproximadamente.
Hoja 5 1

Problemas de E.D.O. (Curso 2008-09)

1. Dados los sistemas lineales

( (
x0 = 2x + 3y − 7 x0 = −x − 2y + 3
y 0 = −x − 2y + 4 y 0 = 2x − y − 6

a) Determı́nese la naturaleza de sus puntos crı́ticos y sus propiedades de estabilidad.

b) En el caso en el que se obtiene un punto de silla, determinar las direcciones de sus
ejes.

2. Estudiar los puntos crı́ticos de los problemas

µ 0¶ µ ¶
x 4x + 2y + 2x2 − 3y 2
1. =
y0 4x − 3y + 7xy
µ 0¶ µ ¶
x x(4 − 2x − y)
2. =
y0 y(3 − x − y)
µ 0¶ µ ¶
x −y − x3
3. =
y0 x − y3
µ 0¶ µ ¶
x −y + x3
4. =
y0 x + y3
µ 0 ¶ µ x+y ¶
x e +y
5. =
y0 y − xy

3. Describir el plano de fases para los sistemas:

µ 0¶ µ 2¶
x x
1. 0 =
y y
µ 0¶ µ ¶
x 2xy
2. =
y0 y 2 − x2

4. Discutir según los valores de µ la estabilidad del (0, 0)para el sistema

µ 0¶ µ ¶
x y
=
y0 −x − µ(x2 − 1)y

(Ecuación de Van der Pol).

5. Mediante el método directo de Liapunov estúdiese el tipo de estabilidad del (0, 0)

en los siguientes sistemas:
Hoja 5 2

 
µ 0¶ 1 3
x  −x + y − 4 x 
1. =  1 1 
y0 − x − y − 2y 3
3 2
µ 0¶ µ ¶
x y − 3x3
2. =
y0 −x − 7y 3
µ 0¶ µ ¶
x −xy 4
3. =
y0 yx4
µ 0¶ µ ¶
x x − xy 4
4. =
y0 y − y 3 x2
µ 0¶ µ ¶
x −2y − x3
5. =
y0 3x − 4y 3

6. En cada uno de estos sistemas determı́nese la naturaleza del punto crı́tico (0, 0) y
sus propiedades de estabilidad:
( (
x0 = x + y − 2 sen(xy) x0 = − sen(x + y) + 1 − e−x2
a) b)
y 0 = −2 x + y + 3 y 2 y 0 = −2 x + 4 y + y sen(x + y)

7. Estúdiense los puntos crı́ticos del siguiente sistema autónomo y esbócese una
posibilidad coherente con estos datos para las trayectorias del sistema en el primer
cuadrante. (
x0 = x (60 − 4 x − 3 y)
y 0 = y (42 − 3 x − 2 y)

8. Considérese el sistema
(
x0 = y + x (x2 + y 2 )
y 0 = −x + y (x2 + y 2 )

a) Estúdiese que tipo de punto crı́tico es (0, 0) en el sistema linealizado.

b) Resuélvase el sistema no lineal empleando coordenadas polares y decı́dase la esta-
bilidad de dicho punto crı́tico.
c) Realizar un análisis análogo (sin calcular explı́citamente las soluciones) en el caso
del sistema µ 0¶ µ ¶
x 3x − y − x exp(x2 + y 2 )
=
y0 x + 3y − y exp(x2 + y 2 )

9. Transfórmese la ecuación del péndulo x00 + sen x = 0 en un sistema autónomo con

el procedimiento habitual (escribiendo y = x0 ).
a) Calcúlense todos los puntos crı́ticos.
b) Hállese la ecuación de las trayectorias
c) Decı́dase la estabilidad y carácter de todos los puntos crı́ticos.
Hoja 5 3

d) ¿Qué función de Liapunov se podrı́a emplear para probar la estabilidad en (0, 0)?

10. Pruébese que el sistema:

(
x0 = −y + x(1 − x2 − y 2 )
y 0 = x + y(1 − x2 − y 2 )

tiene (al menos) una solución periódica rodeando al origen. Estúdiese la estabilidad de
ese punto. Indicación: Utilı́cense coordenadas polares.

11. Determı́nese una función de Liapunov para

(
x0 = −2y − x3
y 0 = x/2 − 4y 3

12. Hállese una función de Liapunov de la forma α(2x + y)2 + β(x + y)2 que pruebe
la estabilidad en el origen del sistema autónomo:
(
x0 = −3x − 2y − (2x + y)3 + (x + y)3
y 0 = 5x + 3y − 2(x + y)3 + (2x + y)3

13. Estúdiese si los siguientes sistemas autónomos tienen alguna solución periódica.
( ( 2 2
x0 = y x0 = y + x3 (2 − ex +y )
a) b) 2 2
y 0 = −y − x2 y + x5 y 0 = −x + y 3 (2 − ex +y )

*14. Pruébese que si (0, 0) es asintóticamente estable para x0 = F (x, y), y 0 = G(x, y);
entonces es inestable para x0 = −F (x, y), y 0 = −G(x, y). Utilı́cese este hecho para
enuciar un teorema como el de Liapunov cuya conclusión sea la inestabilidad. Indicación:
Considérese la inversión temporal t 7→ −t.
Hoja 5 1

Problemas de E.D.O. (Curso 2009-2010)

1. Dados los sistemas lineales

( (
x0 = 2x + 3y − 7 x0 = −x − 2y + 3
y 0 = −x − 2y + 4 y 0 = 2x − y − 6

a) Determı́nese la naturaleza de sus puntos crı́ticos y sus propiedades de estabilidad.

b) En el caso en el que se obtiene un punto de silla, determinar las direcciones de sus
ejes.

2. Estudiar los puntos crı́ticos de los problemas

µ 0¶ µ ¶
x 4x + 2y + 2x2 − 3y 2
1. =
y0 4x − 3y + 7xy
µ 0¶ µ ¶
x x(4 − 2x − y)
2. =
y0 y(3 − x − y)
µ 0¶ µ ¶
x −y − x3
3. =
y0 x − y3
µ 0¶ µ ¶
x −y + x3
4. =
y0 x + y3
µ 0 ¶ µ x+y ¶
x e +y
5. =
y0 y − xy

3. Describir el plano de fases para los sistemas:

µ 0¶ µ 2¶
x x
1. 0 =
y y
µ 0¶ µ ¶
x 2xy
2. =
y0 y 2 − x2

4. Discutir según los valores de µ la estabilidad del (0, 0)para el sistema

µ 0¶ µ ¶
x y
=
y0 −x − µ(x2 − 1)y

(Ecuación de Van der Pol).

5. Mediante el método directo de Liapunov estúdiese el tipo de estabilidad del (0, 0)

en los siguientes sistemas:
Hoja 5 2

 
µ 0¶ 1 3
x  −x + y − 4 x 
1. =  1 1 
y0 − x − y − 2y 3
3 2
µ 0¶ µ ¶
x y − 3x3
2. =
y0 −x − 7y 3
µ 0¶ µ ¶
x −xy 4
3. =
y0 yx4
µ 0¶ µ ¶
x x − xy 4
4. =
y0 y − y 3 x2
µ 0¶ µ ¶
x −2y − x3
5. =
y0 3x − 4y 3

6. En cada uno de estos sistemas determı́nese la naturaleza del punto crı́tico (0, 0) y
sus propiedades de estabilidad:
( (
x0 = x + y − 2 sen(xy) x0 = − sen(x + y) + 1 − e−x2
a) b)
y 0 = −2 x + y + 3 y 2 y 0 = −2 x + 4 y + y sen(x + y)

7. Estúdiense los puntos crı́ticos del siguiente sistema autónomo y esbócese una
posibilidad coherente con estos datos para las trayectorias del sistema en el primer
cuadrante. (
x0 = x (60 − 4 x − 3 y)
y 0 = y (42 − 3 x − 2 y)

8. Considérese el sistema
(
x0 = y + x (x2 + y 2 )
y 0 = −x + y (x2 + y 2 )

a) Estúdiese que tipo de punto crı́tico es (0, 0) en el sistema linealizado.

b) Resuélvase el sistema no lineal empleando coordenadas polares y decı́dase la esta-
bilidad de dicho punto crı́tico.
c) Realizar un análisis análogo (sin calcular explı́citamente las soluciones) en el caso
del sistema µ 0¶ µ ¶
x 3x − y − x exp(x2 + y 2 )
=
y0 x + 3y − y exp(x2 + y 2 )

9. Transfórmese la ecuación del péndulo x00 + sen x = 0 en un sistema autónomo con

el procedimiento habitual (escribiendo y = x0 ).
a) Calcúlense todos los puntos crı́ticos.
b) Hállese la ecuación de las trayectorias
c) Decı́dase la estabilidad y carácter de todos los puntos crı́ticos.
Hoja 5 3

d) ¿Qué función de Liapunov se podrı́a emplear para probar la estabilidad en (0, 0)?

10. Pruébese que el sistema:

(
x0 = −y + x(1 − x2 − y 2 )
y 0 = x + y(1 − x2 − y 2 )

tiene (al menos) una solución periódica rodeando al origen. Estúdiese la estabilidad de
ese punto. Indicación: Utilı́cense coordenadas polares.

11. Determı́nese una función de Liapunov para

(
x0 = −2y − x3
y 0 = x/2 − 4y 3

12. Hállese una función de Liapunov de la forma α(2x + y)2 + β(x + y)2 que pruebe
la estabilidad en el origen del sistema autónomo:
(
x0 = −3x − 2y − (2x + y)3 + (x + y)3
y 0 = 5x + 3y − 2(x + y)3 + (2x + y)3

13. Estúdiese si los siguientes sistemas autónomos tienen alguna solución periódica.
( ( 2 2
x0 = y x0 = y + x3 (2 − ex +y )
a) b) 2 2
y 0 = −y − x2 y + x5 y 0 = −x + y 3 (2 − ex +y )

*14. Pruébese que si (0, 0) es asintóticamente estable para x0 = F (x, y), y 0 = G(x, y);
entonces es inestable para x0 = −F (x, y), y 0 = −G(x, y). Utilı́cese este hecho para
enuciar un teorema como el de Liapunov cuya conclusión sea la inestabilidad. Indicación:
Considérese la inversión temporal t 7→ −t.
 lOMoARcPSD|2192898

Hoja1 - HOJA 1

Ecuaciones Diferenciales (Universidad Autónoma de Madrid)

Distributing prohibited | Downloaded by Maria Fernandez ([email protected])

 lOMoARcPSD|2192898

Problemas E.D.O. (2 de Matemáticas). Curso 2016-2017.

Los recomendados estan marcados con un (*)

1. (*) a) Comprobar que, para cada valor de la constante C, la función

Z x
2 2
y = ex (C + e−t dt)
0

es una solución de la ecuación diferencial

y ′ = 2xy + 1.

b) Sean y1 (x), y2 (x) dos funciones de la familia anterior, correspondientes a dos valores distintos
de la constante C. Hallar la ecuación diferencial que satisface u(x) = y1 (x) − y2 (x).

2. (*) a) Hallar los valores de m para los cuales y = emx es solución de la ecuación diferencial

2y ′′′ + y ′′ − 5y ′ + 2y = 0.

b) A partir de combinaciones lineales de las soluciones encontradas en el apartado anterior,

hallar una solución que verifique las condiciones

y(0) = 0, y ′ (0) = 1, y ′′ (0) = −1.

3. (*)Comprobar que para cada valor de la constante C > 0, la identidad

Z x
sent
Cx − y dt = 0
0 t
define y como función de x. Hallar la ecuación diferencial que satisfacen todas las funciones
de la familia. Hallar y(0) e y ′ (0) para cada una de ellas.
4. (*) a) Utilizando las isoclinas, esbozar las soluciones de y ′ = y 2 − 1.
b) Resolver explı́citamente la ecuación y ′ = y 2 − 1 , y(0) = 0 y comparar el resultado con lo
obtenido por el método de isoclinas empleado anteriormente.

5. (*) Trazando algunas isoclinas, esbozar las soluciones de las siguientes ecuaciones:
p
a) y ′ = sen(y + x), b) x′ = t2 + x 2

6. (*) a) Describir geométricamente la familia de curvas

x2 + y 2 + 2Cx = 0

y calcular una ecuación diferencial que satisfacen todas ellas.

b) Determinar la ecuación diferencial de la familia de curvas

y(x) = x cos(x + C).

Distributing prohibited | Downloaded by Maria Fernandez ([email protected])

 lOMoARcPSD|2192898

Hoja 1 2

7. (*) Dada la ecuación diferencial

y ′ = (1 + x)y + 1 − 3x + x2 ,

calcular los primeros términos del desarrollo de Taylor de la solución que satisface y(0) = 0.

8. (*) Esbozar las siguientes familias uniparamétricas de curvas y hallar sus trayectorias or-
togonales:
a) xy = C.
b) y = Cex .
c) y = Cxn , donde n es un entero positivo. Explicar qué sucede con las trayectorias ortogonales
cuando aumentamos el valor del entero n.

9. (*) Hallar la familia de curvas ortogonales a la familia de parábolas

y 2 − Cx = C 2 /4.

10. (*) Hallar la familia de curvas ortogonales a la familia de circunferencias definida por

x2 + (y − C)2 = C 2 .

Interpretar el resultado geométricamente.

11. (*) Hallar las trayectorias ortogonales de las siguientes familias uniparamétricas de curvas
expresadas en coordenadas polares:

a) r = C(1 + cos θ). b) r = 2Csenθ.

12. (*) Hallar las curvas que satisfacen las condiciones geométricas siguientes:
a) El segmento de la tangente limitado por los ejes coordenados tiene como punto medio al
punto de tangencia.
b) La proyección sobre el eje OX de la parte de la tangente entre (x, y) y el eje OX tiene
longitud 1.
c) El ángulo entre el radio polar y la tangente es constante.
d) La curva pasa por (0, 0) y está contenida en el primer cuadrante, de modo que el área bajo
la curva desde (0, 0) hasta (x, y) es un tercio del área del rectángulo que tiene a esos puntos
como vértices opuestos.

13. (*) Para cada una de las ecuaciones diferenciales siguientes, hallar la solución particular
que satisface la condición inicial dada:
• y ′ = e3x−2y , con y(0) = 0.

• e−y + (1 + x2 )y ′ = 0, con y(0) = 0.

• xyy ′ = (x + 1)(y + 1), con y(1) = 0.

14. (*) Según la Ley de enfriamiento de Newton la tasa de variación de la temperatura en un

cuerpo es proporcional a la diferencia de temperatura con el ambiente. Si una barra de hierro

Distributing prohibited | Downloaded by Maria Fernandez ([email protected])

 lOMoARcPSD|2192898

Hoja 1 3

a 100C se enfrı́a a 90C en 5 segundos cuando se deja a una temperatura ambiente de 20C,
¿cuánto tardará en estar a 30C?

15. En una sala que está a una temperatura de 20C nos sirven dos tazas de café, a una
temperatura de 40C. Disponemos de leche frı́a, a una temperatura de 10C. En una de las tazas
echamos una cantidad de leche frı́a igual a la de café contenido en su interior, y esperamos cinco
minutos. En la otra taza esperamos cinco minutos, y después agregamos la leche. Determinar
cual de las dos tazas tiene el café con leche más caliente.

16. (*) Una bola de naftalina tiene inicialmente un radio de 1 cm. Al cabo de un mes su radio
se ha reducido a 0,5 cm. Suponiendo que la naftalina se evapora a un ritmo proporcional a la
superficie de contacto con el aire, hallar la evolución del radio de la bola en función del tiempo.

17. Un dı́a comenzó a nevar por la mañana y siguió cayendo la nieve de forma constante todo
el dı́a. A las 12 del mediodı́a una quitanieves comenzó a limpiar una carretera, con velocidad
inversamente proporcional al espesor de la nieve depositada. Sabiendo que a las 2 de la tarde
habı́a limpiado 2 km., y que a las 4 de la tarde habı́a limpiado 1 km. más, determinar a qué
hora comenzó a nevar.

18. Cuatro hormigas situadas en las esquinas de una mesa cuadrada de lado 1 comienzan
a andar simultaneamente a la misma velocidad, cada una en la dirección de su vecina más
próxima en la dirección contraria a las agujas del reloj. Tomando coordenadas polares con
origen en el centro de la mesa y eje polar a lo largo de una diagonal, hallar la trayectoria de la
hormiga que parte del eje polar.

19. (*) Una población de bacterias que sigue la Ley de Malthus (la tasa de variación es
proporcional al número de individuos) se duplica al cabo de 24 horas. ¿Cuánto tardará en
triplicarse?

20. (*) Supongamos que una población sigue el modelo p′ = bp2 − ap con a, b > 0.
a) Representar en un diagrama de fases, clasificando los puntos crı́ticos del sistema.
b) Demostrar que si p(t0 ) < a/b, entonces la población tiende a extinguirse.

21. (*)Dada la ecuación y ′ = cos y, sin encontrar las soluciones explı́citamente, estudiar su
comportamiento cualitativo (en particular, clasificar los puntos crı́ticos según su estabilidad).

Distributing prohibited | Downloaded by Maria Fernandez ([email protected])

 EDO
NOTA
CURSO ACADÉMICO 2008-2009

Nombre y Apellidos
Grupo

Problema 1 Dada la ecuación diferencial ordinaria

dy
= e−y ,
dx
a) Esbozar razonadamente el campo de pendientes utilizando las isoclinas, y esbozar las soluciones.
b) Resolver explı́citamente (separando variables) la ecuación con dato inicial y(0) = 0, dibujarla aproximada-
mente y comparar el resultado con lo obtenido por el método de isoclinas empleado anteriormente.
 3
Problema 2 Dada la familia de curvas definida por y = (x + C) , con C ∈ R, hallar la familia de trayectorias
ortogonales. Esbozar razonadamente las dos familias.
 EDO
CURSO ACADÉMICO 2008-2009
NOTA 10
Nombre y Apellidos
Grupo

Problema 1 Resolver la siguiente ecuación diferencial hallando el factor integrante oportuno.

· ¸
x x (y + 1)ey
e
|{z} dx − e tan(y) + dy = 0.
cos(y)
M (x,y) | {z }
N (x,y)

Solución. La ecuación dada no es exacta:

· ¸
∂M ∂N 1 ∂N ∂M ex tan(y) sin(y)
= 0 6= = −ex tan(y) ⇒ − =− =−
∂y ∂x M ∂x ∂y ex cos(y)
∂µ ∂µ
Busco un factor integrante µ(x, y) que dependa solo de y: µ(x, y) = µ(y), entonces ∂x =0y ∂y = µ0 (y):
· ¸ · ¸
∂µ M ∂µ N ∂µ ∂N ∂M µ0 (y) 1 ∂N ∂M sin(y)
= ⇐⇒ M =µ − ⇐⇒ = − =−
∂y ∂x ∂y ∂x ∂y µ(y) M ∂x ∂y cos(y)

Integrando en dy encuentro la expresión del factor integrante

£ ¤ £ ¤
log µ(y) = log cos(y) + C ⇒ µ(y) = cos(y)

dado que la constante puedo elegirla como quiero, he puesto C = 0. Ahora es fácil averiguar que la ecuación
f dx + N
M e dy = µ M dx + µ N dy = ex cos(y) dx − [ex sin(y) + (1 + y)ey ] dy = 0

es exacta:
f
∂M e
∂N
= −ex sin(y) = = −ex sin(y) .
∂y ∂x
Entonces existe un potencial P (x, y) tal que ∂P fy
=M ∂P e , y puedo encontrarlo integrando:
=N
∂x ∂y
Z Z Z
∂P f dx + c(y) = cos(y)
P (x, y) = dx + c(y) = M ex dx + c(y) = ex cos(y) + c(y)
∂x
e puedo encontrar c(y)
usando la expresión de N

e = ∂P = −ex sin(y) + c0 (y)

−ex sin(y) − (1 + y)ey = N ⇒ c0 (y) = −(1 + y)ey ⇒ c(y) = −yey + K
∂y
donde K es una constante arbitraria. El potencial será

P (x, y) = ex cos(y) − yey + K

y las soluciones del problema serán dadas para la ecuación implicita P (x, y) = K, es decir la solución y(x) es
definida implicitamente en la ecuación

P (x, y) = K ⇐⇒ ex cos(y) − yey + K = K ⇐⇒ ex cos(y) − yey = K − K = C ,

donde C ∈ R es una constante arbitraria que define la familia unı́parametrica de soluciones. Se podrı́a explicı́tar
en términos de x: x = log [(C + yey )/ cos(y)] , pero solo para los y tales que cos(y) 6= 0. ¤
Problema 2 Dada la siguiente ecuación diferencial ordinaria del segundo orden
£ 0 ¤2
00 0 y (x)
y (x) = 2 y (x) − , con y(x) > 0, y 0 (x) > 0, ∀x ∈ R
y(x)

(a) Reducirla a una ecuación del primer orden y resolverla. Escribir la solución y(x) en forma esplicita.
(b) Encontrar aquella solución que satisface las condiciones iniciales y(0) = 1 , y 0 (0) = 1. ¿Es unica?

dz dy
Solución. Ponemos z = y 0 ; entonces y 00 = z 0 = dz 0 dz dz
dx = dy dx = y dy = z dy , es decir: hemos puesto y como
variable independiente y hemos hecho depender z de y en lugar que de x. Substituyo en la ecuación

dz(y) z 2 (y) dz(y) 1

z(y) = 2 z(y) − ⇐⇒ + z(y) −2 = 0
dy y dy y |{z}
|{z} b(y)
a(y)

Que es una ecuación lineal del primer orden completa para z(y), cuya solución es:
(la formula se obtiene en varias maneras: factor integrante, variación de la constante...)
R ½ Z y R ¾ R ½ Z y R ¾
− y a(ξ)dξ η
a(ξ)dξ − y 1 dξ η 1
dξ
z(y) = e y0 z(y0 ) − e y0 b(η)dη = e y0 ξ z(y0 ) − e y0 ξ (−2) dη
y0 y0
½ Z y ¾ ½ Z y ¾
− log(y)+log y0 log(η)−log(y0 ) y 0 2
=e z(y0 ) + 2 e dη = z(y0 ) + ηdη
y0 y y0 y0
1© ª C + y2
= y0 z(y0 ) − y02 + y 2 =
y y

Regresando a la variable y, recordando que C = y0 z(y0 ) − y02 = y0 y00 − y02 , (y recuerdo las notaciones habituales
y0 = y(x0 ) , y00 = y 0 (x0 ))
Z x Z x
C + y2 y y0 y(ξ) y 0 (ξ)
y 0 (x) = ⇐⇒ = 1 ⇒ dξ = 1dξ = x − x0
y C + y2 2
x0 C + y (ξ) x0

R 0
resolviendo el integral (recuerdo que ff = log f ) obtengo
(también se podia resolver como una Bernoulli con parametro −1)
µ ¶ q
1 C + y2 2 2 2(x−x0 )
log = x − x 0 ⇐⇒ C + y = (C + y0 )e ⇐⇒ y(x) = (C + y02 )e2(x−x0 ) − C
2 C + y02
q
⇐⇒ y(x) = (y0 y00 − y02 + y02 )e2(x−x0 ) − y0 y00 + y02

finalmente la solución general es dada por

q
y(x) = y0 y00 (e2(x−x0 ) − 1) + y02 ,

y depende de DOS paramétros: y0 = y(x0 ) y y00 = y 0 (x0 ). La solución que satisface las condiciones del punto
(b) se obtiene ponendo x0 = 0, y(x0 ) = y0 = y(0) = 1, y 0 (x0 ) = y00 = y 0 (0) = 1 en la espresión de y
√
y(x) = e2x = ex .

esta solución es la unica que cumple las condiciones iniciales del apartado (b). ¤
 EDO
NOTA
CURSO ACADÉMICO 2008-2009

Nombre y Apellidos
Grupo

Problema 1 Resolver la siguiente ecuación diferencial del segundo orden

1 00
y + 3 y 0 + 25y = cos(2t) .
4
Solución. El polinomio caractéristico es p(λ) = 14 λ2 + 3λ + 25, y sus soluciones son λ± = −6 ± 8i. Por lo tanto
la solución general de la ecuación homogénea es
£ ¤
yC (t) = c1 eλ+ t + c2 eλ− t = e−6t c1 e8i t + c2 e−8i t

donde c1 , c2 ∈ C son constantes arbitrarias, pero tales que la solución sea real. Es conveniente encontrar dos
soluciones linealmente independientes que sean reales, y eso se puede hacer eligiendo

1 e8i t + e−8i t
c1 = c2 = ⇒ y1 (t) = e−6t = e−6t cos(8t)
2 2
y tambien
1 1 e8i t − e−8i t
c1 = , c2 = − ⇒ y2 (t) = e−6t = e−6t sin(8t)
2i 2i 2i
las dos soluciones son independientes dado que su determinante Wronskiano
· ¸ · ¸
y1 (t) y2 (t) e−6t cos(8t) e−6t sin(8t)
W (t) = det = det 6= 0
y10 (t) y20 (t) −8e−6t sin(8t) − 6e−6t cos(8t) 8e−6t cos(8t) − 6e−6t sin(8t)

es distinto de cero en t = 0 (támbien en 8t = π/2) entonces es distinto de cero en todo punto. Por lo tanto y1 y
y2 son independientes y la solución general de la ecuación homogénea es:

yR (t) = k1 e−6t cos(8t) + k2 e−6t sin(8t)

para cada k1 , k2 ∈ R. Ahora tenemos que buscar una solución particular, y la buscamos de la forma

yP (t) = A cos(2t) + B sin(2t)

yP0 (t) = −2A sin(2t) + 2B cos(2t)
yP00 (t) = −4A cos(2t) − 4B sin(2t)

ahora es claro que

1 00
y + 3 yP0 + 25yP = cos(2t) ⇐⇒ (24A + 6B) cos(2t) + (24B − 6A) sin(2t) = cos(2t)
4 P
entonces se llega a la conclusion que yP satisface la ecuación solo cuando A = 4/102 y B = 1/102, por lo tanto
la solución general de la ecuación es
4 1
yG (t) = yR (t) + yP (t) = k1 e−6t cos(8t) + k2 e−6t sin(8t) + cos(2t) + sin(2t) ,
102 102
donde k1 , k2 ∈ R.
Problema 2 Dado el siguiente sistema homogéneo con coeficientes constantes
 
5 −1 0
X0 =  4 1 0  X .
1 1 2

Hallar la solución general. Encontrar aquella solución tal que X(0) = (1, 2, 0). Es unica?
Solución. El polinomio caractéristico es det(A − λI) = (λ − 2)(λ − 3)2 , entonces tenemos un autovalor λ1 = 2
con multiplicidad algebrica y geometrica 1, y un autovalor λ2 = 3 con multiplicidad algebrica 2.
Buscamos una matriz P de cambio de variable que transforme A en su forma de Jordan J:
   
5 −1 0 2 0 0
A= 4 1 0  ←→ J = P −1 A P =  0 3 1 
1 1 2 0 1 3

Calculamos los autovectores: v1 es el autovector correspondiente a λ1 = 2, por lo tanto es solución del sistema
(A − 2I)v = 0, y tiene la forma v1 = (0, 0, z), y elı́go v1 = (0, 0, 1). v2 es un autovector correspondiente a
λ2 = 3, por lo tanto es solución del sistema (A−3I)v = 0, y tiene la forma v1 = (x, 2x, 3x), y elı́go v1 = (1, 2, 3).
Está claro que la dimension del autoéspacio correspondiente a v2 es 1, es decir su multiplicidad geometrica es 1.
Entonces la matriz no es diagonalizable, siendo la multiplicidad algebrica 2 y la geometrica 1. Busco entonces el
autovector generalizado w, que es solución del sistema (A−3I)v = v2 , que tiene la forma w = (x, 2x−1, 3x−4),
y elı́go w = (1, 1, −1). Por lo tanto la matriz P es dada por
   
0 1 1 5 −4 1
P =  0 2 1 , P −1 =  −1 1 0 
1 3 −1 2 1 0

La matriz de Jordan J = D + N ,
   
2 0 0 e2t 0 0
D= 0 3 0 , eD t= 0 e3t 0 ,
0 0 3 0 0 e3t
   
0 0 0 2 1 0 0
t
N = 0 0 1  eN t N2
= I + N t + |{z} + ... =  0 1 t 
2
0 0 0 =0 0 0 1
dado que N 2 = 0. Adémas DN = N D, es decir las dos matrices N y D comútan, entonces
 2t 
e 0 0
eJt = eDt+N t = eDt eN t =  0 e3t t e3t 
0 0 e3t

En fin calculo la matriz exponencial:

 
(1 + 2t)e3t −te3t 0
eAt = P eJt P −1 =  4e3t (1 − 2t)e3t 0 
(6t − 5)e3t + 5e2t (−3t + 4)e3t − 4e2t e2t

y la solución general tendrá la forma X(t) = eAt X(0). Si X(0) = (1, 2, 0), entonces la unica solución tiene la
forma
     
(1 + 2t)e3t −te3t 0 1 e3t
X(t) = eAt X(0) =  4e3t (1 − 2t)e3t 0   2 = 2e3t 
3t 2t 3t 2t 2t
(6t − 5)e + 5e (−3t + 4)e − 4e e 0 3e − 3e2t
3t

.
 EDO
NOTA
CURSO ACADÉMICO 2008-2009

Nombre y Apellidos
Grupo

Problema 1 Sea la sucesión de funciones {fn }∞

n=1 dada por

³ ´ µ ¶
1 1
fn (x) = n 1 + x n sin e−nx .
n

Probar que converge uniformemente en [1, ∞). Probar que converge puntualmente, pero no uniformemente, en
[0, ∞).
Solución. Antes que todo miramos la convergencia puntual en [0, +∞): para cada x ∈ [0, +∞) fijo, la sucesión
numerica fn (x) converge a
 ³ 1
´ ¡ ¢
 n 1 + x n sin n1 e−nx → 0 si x > 0
fn (x) =
 ¡ ¢
n sin n1 → 1 si x = 0
¡ ¢
dado que la exponencial domina si x > 0 y dado que n sin n1 → 1 para n → +∞. El limite es discontinuo en
[0, +∞) por lo tanto la sucesión no puede converger uniformemente en todo [0, +∞): si la sucesión fn de funciones
continuas, fuese uniformente convergente en todo [0, +∞), el limite deberia ser continuo en todo [0, +∞), y ya
sabemos que no lo es.
De otra manera se podia comprobar que el estremo superior

0≤ sup |fn (x)| = 1

x∈[0,+∞)

por lo tanto no puede converger a cero en [0, +∞) y la convergencia no puede ser uniforme.
Analizamos ora si hay convergencia uniforme en [1, +∞): es suficiente provar que

0≤ sup |fn (x)| ≤ cn → 0

x∈[1,+∞)

para una sucesión cn positiva y independiente de x. Para todo x ≥ 1,

³ ´ ¯¯ µ 1 ¶¯¯
1
|fn (x)| = n 1 + x n ¯¯sin ¯ e−nx ≤(i) 2nx n1 e−nx ≤(ii) 2nx n1 c1 = 2c 1 1 ≤(iii) 2c → 1
n ¯ (nx)5 n4 x5− n n4

donde:
(i) he usado el hecho que x ≥ 1, ası́ que 1 + x1/n ≤ 2x1/n , y que | sin(1/n)| ≤ 1.
(ii) he usado el hecho que e−nx ≤ c/(nx)6 para una constante c > 0.
1
(iii) he usado el hecho que x ≥ 1, ası́ que 1/x5− n ≤ 1.
Elijo cn = 2c/n4 y eso me garantiza la convergencia uniforme en todo [1, +∞).
 Problema 2 Se considera el problema:

 y 0 = 1 y 4 + 1 y + 1 |y| 14 = f (t, y)
3 61+t 3

y(0) = 1.

a) ¿ Qué puedes decir sobre existencia y unicidad de las soluciones ?

1 y
b) Demostrar que y(t) ≥ (1 + t) 6 para todo t ≥ 0. (Indicación: comparar y 0 con 16 1+t .)
c) Supongamos que y(t) es una solución definida en un intervalo [0, T ). Demostrar que T ≤ 1.
(Indicación: comparar y 0 con 13 y 4 .)
d) Demostrar que la solución está definida en el intervalo [0, 13 ). (Indicación: comparar y 0 con y 4 .)
Solución. (a) La función f (t, y) es continua para (t, y) ∈ [0, +∞) × R, por lo tanto el Teorema de Peano garantiza
la existencia local en un intervalo maximal [0, T ) ⊆ [0, +∞). La función f (t, y) no es Lipschitz para todo (t, y) ∈
[0, +∞) × R, dado que
∂f 4 1 1 1 y
(t, y) = y 3 + +
∂y 3 6 1 + t 12 |y| 74

y es claro que ∂f ∂y (t, 0) = ∞. Pero el dato inicial es y(0) = 1, y la función f (t, y) es Lipschitz en la banda
(t, y) ∈ [0, +∞) × [1, Y ], para todo Y ≥ 1, de hecho
¯ ¯ ¯ ¯
¯ ∂f ¯ ¯4 3 1¯
|f (t, y1 ) − f (t, y2 )| ≤ sup ¯ ¯ ¯ ¯
¯ (t, y)¯ |y1 − y2 | ≤ ¯ 3 Y + 3 ¯ |y1 − y2 | = L|y1 − y2 |.
(t,y)∈[0,+∞)×[1,Y ] ∂y

lo cual asegura la unicidad local.

Entonces concluimos que hay existencia y unicidad local en [0, T ), pero no tenemos informaciones acerca de T > 0.
1 y
(b) Dado que f (t, y) ≥ entonces tenemos que
61+t
 
 y 0 = 1 y 4 + 1 y + 1 |y| 14 = f (t, y)  y0 = 1 1 y
3 61+t 3 se compara con 61+t (1)
 
y(0) = 1. y(0) = 1.

y sabemos que y(t) ≥ y(t) para todo t ≥ 0, dado que y(0) = y(0) = 1. Integrando la ecuación para y encuentro que
1 1
y(t) = (1 + t) 6 ≥ 1 y puedo concluir que y(t) ≥ y(t) = (1 + t) 6 ≥ 1, para todo t ≥ 0.
1 4
(c) Dado que f (t, y) ≥ y entonces tenemos que
3
 
 y 0 = 1 y 4 + 1 y + 1 |y| 41 = f (t, y)  y0 = 1 y4
3 61+t 3 se compara con 3 (2)
 
y(0) = 1. y(0) = 1.

y sabemos que y(t) ≥ y(t) para todo t ≥ 0, dado que y(0) = y(0) = 1. Integrando la ecuación para y , encuentro
1
que y(t) = (1 − t)− 3 y descubro que y(t) → +∞ cuando t → 1, por lo tanto y(t) (que es mayor que y(t)) támbien
explota cuando t → 1. Eso quiere decir que, siendo y(t) definida en todo [0, T ), entonces necesariamente T ≤ 1.

(d) Dado que f (t, y) ≤ y 4 , siendo siempre y(t) ≥ 1 gracias al apartado (b), entonces tenemos que
 ( 0
 y 0 = 1 y 4 + 1 y + 1 |y| 14 = f (t, y) y = y4
3 61+t 3 se compara con (3)
 y(0) = 1.
y(0) = 1.

y sabemos que y(t) ≤ y(t) para todo t ≥ 0, dado que y(0) = y(0) = 1. Integrando la ecuación para y , encuentro
1
que y(t) = (1 − 3t)− 3 y descubro que y(t) → +∞ cuando t → 1/3, por lo tanto y(t) (que es menor que y(t))
esta acotada cuando t < 1/3. Eso quiere decir que, siendo y(t) definida en todo [0, T ), entonces estoy seguro que
T ≥ 1/3.

Concluyendo, hemos garantizado que existe una unica solución y(t) ≥ 1 en un intervalo [0, T ), donde 1/3 ≤ T ≤ 1.
 EDO
NOTA
CURSO ACADÉMICO 2008-2009

Nombre y Apellidos
Grupo

Problema 1 Dado el sistema

 ¡ ¢
 x0 (t) = y(t) ex2 (t)−1 − x(t) ex(t)−y(t) = F x(t), y(t)
 y 0 (t) = x2 (t) y(t) − y(t) = G¡x(t), y(t)¢

Determinar la naturaleza de los (tres) puntos crı́ticos utilizando el método de la linealización.

Dibujar las trayectorias ( en el plano de fase (x, y) ) de los sistemas linealizados, y también del sistema completo.
¿ Qué sabes decir sobre la estabilidad de estos puntos crı́ticos, para los sistemas lineales y para el completo ?
Solución. Los puntos crı́ticos son los puntos (x, y) tales que F (x, y) = G(x, y) = 0. Es facı́l averiguar que los
unı́cos puntos crı́ticos son: P0 (0, 0), P1 = (1, 1), P2 = (−1, −1). La ecuación diferencial de las trayectorias en el
plano de fases:
dy G(x, y) x2 y − y
= = 2
dx F (x, y) y ex −1 − x ex−y
las soluciones representan, en el plano de fases (x, y), las trayectorias del sistema“lejos” de los puntos crı́ticos. Se
puede integrar lejos de los puntos crı́ticos? es exacta? (ejercicio de la parte 2).
Antes que continuar es practico calcular unas derivadas parciales:
2
∂F
∂x (x, y) = ∂x F (x, y) = 2xyex −1
− (x + 1)ex−y , ∂x F (P0 ) = −1, ∂x F (P1 ) = 0, ∂x F (P2 ) = 2,
2
∂F
∂y (x, y) = ∂y F (x, y) = ex −1
+ xex−y , ∂y F (P0 ) = e−1 , ∂y F (P1 ) = 2, ∂y F (P2 ) = 0,
∂G
∂x (x, y) = ∂x G(x, y) = 2xy, ∂x G(P0 ) = 0, ∂x G(P1 ) = 2, ∂x G(P2 ) = 2,
∂G
∂y (x, y) = ∂y G(x, y) = x2 − 1, ∂y G(P0 ) = −1, ∂y G(P1 ) = 0, ∂y G(P2 ) = 0,

Para intender lo que pasa en los puntos crı́ticos, se tiene que linealizar al rededor de ellos: pongo Pi = (xi , yi ) y
( 0 ¡ ¢
x = F x(t), y(t) = ∂x F (Pi )(x − xi ) + ∂y F (Pi )(y − yi ) + Err[F ](x, y)
¡ ¢
y 0 = G x(t), y(t) = ∂x G(Pi )(x − xi ) + ∂y G(Pi )(y − yi ) + Err[G](x, y)

donde he definido Err[H](x, y) := F (x, y) − ∂x H(Pi )(x − xi ) − ∂y H(Pi )(y − yi ).

La traslación (x − xi ) 7→ X y (y − yi ) 7→ Y nos da el sistema lineal:
( 0
X (t) = ∂x F (Pi )X + ∂y F (Pi )Y

Y 0 (t) = ∂x G(Pi )X + ∂y G(Pi )Y

cuyo comportamiento alrededor de (X, Y ) = (0, 0), nos dice el comportamento del sistema no-lineal al rededor del
punto critico Pi , siempre cuando los errores son pequeños, es decir cuando
Err(F )(x, y) Err(G)(x, y)
p → 0, y p → 0, para (x, y) → (xi , yi )
(x − xi )2 + (y − yi )2 (x − xi )2 + (y − yi )2

Dado que las funciones F y G son regulares (al menos de clase C 2 ), entonces el desarrollo de Taylor de orden 1 (es
decir, la linealización) tiene un error pequeño:
³p ´
Err[H](x, y) = o (x − xi )2 + (y − yi )2 para H = F, G.
En el caso de P0 es facı́l averiguarlo directamente:
¯ ¯ ¯ ¯ ¯ ³ ´ ¯
¯ Err(F )(x, y) ¯ ¯ y ex2 −1 − x ex−y + x − e−1 y ¯ ¯ x (ex − 1) + e−1 ex2 − 1 y ¯
¯ ¯ ¯ ¯ ¯ ¯
0 ≤ lı́m ¯ p ¯= lı́m ¯ p ¯= lı́m ¯ p ¯
(x,y)→P0 ¯ 2
x +y 2 ¯ (x,y)→(0,0) ¯ 2
x +y 2 ¯ (x,y)→(0,0) ¯ 2
x +y 2 ¯
¯ ¯
x2 + e−2 x2 |y| (x2 + y 2 )(1 + e−2 y)
= lı́m p ≤ lı́m p =0
(x,y)→(0,0) x2 + y 2 (x,y)→(0,0) x2 + y 2
2
dado que (ex − 1) ∼ x y (ex − 1) ∼ e−1 x2 cuando x → 0. El ¡ error de G ¢es analogo. En los otros puntos hay que
tener cuidado si se quiere hacerlo en las variables (X, Y ) = x − xi , y − yi → (0, 0).
Ahora analizamos separadamente los sistemas lineales correspondientes a P0 , P1 , P2 .
• P0 = (0, 0). Analizamos alrededor de (X, Y ) = (0, 0) el comportamiento del sistema lineal
½ 0
X = −X + e−1 Y
un autovalor doble λ = −1 y un solo autovector v1 = (1, 0)
Y 0 = −Y

Estamos en el caso “Jordan”, y la solución general es dada con la matriz exponencial Z(t) = eAt Z(0), donde
→ ¡
− ¢
Z = X(t), Y (t) y A es la matriz del sistema:
·µ ¶ µ ¶ ¸ µ −t ¶µ ¶ µ −t ¶
−1 0 0 1 e 0 1 e−1 t e e−1 te−t
Z(t) = exp t+ e−1 t Z(0) = Z(0) = Z(0)
0 −1 0 0 0 e−t 0 1 0 e−t

la solución general se puede escribir también como

µ ¶ µ ¶ · µ ¶ µ −1 ¶¸
X(t) c1 e−t + c2 e−1 te−t −t 1 e t
= =e c1 + c2
Y (t) c2 e−t 0 1

con v1 = (1, 0) y v2 = (e−1 t, 1). Por lo tanto el punto critico (X, Y ) = (0, 0) es un nodo impropio estable (lineas
rectas que se van hace al origen) y es un atractor para el sistema lineal.
Usando el Teorema de linealización, cuyas hypótesis están verificadas, podemos garantizar que el punto critico
nolineal P0 = (0, 0) o bien es un nodo impropio estable o bien nodo estable para el sistema nolineal, de todas formas
se queda un atractor. (Dibujar el lineal y las dos posibilidades del nolineal)
• P1 = (1, 1). Analizamos alrededor de (X, Y ) = (0, 0) el comportamiento del sistema lineal
½ 0
X = 2Y
dos autovalores λ1 = 2, λ2 = −1 y dos autovectores v1 = (1, 1) v2 = (1, −1)
Y 0 = 2X

La solución general es dada por

Z(t) = c1 e2t v1 + c2 e−2t v2
Por lo tanto el punto critico (X, Y ) = (0, 0) es un punto de silla con ejes principales v1 = (1, 1) (repulsor) y
v2 = (1, −1) (atractor) y en general no es ni atractor ni repulsor, por lo tanto es inestable. Las unicas soluciones
que lo alcanzarán son las que corresponden a c1 = 0, es decir Z(t) = c2 e−2t v2 .
Usando el Teorema de linealización, cuyas hypótesis están verificadas, podemos garantizar que el punto critico
nolineal P0 = (1, 1) es un punto de silla para el sistema nolineal, por lo tanto no es ni atractor ni repulsor, tambien
para el sistema no-lineal. (Dibujar el lineal y el nolineal)
• P2 = (−1, −1). Analizamos alrededor de (X, Y ) = (0, 0) el comportamiento del sistema lineal
½ 0
X = 2X
dos autovalores λ1 = 0, λ2 = 2 y dos autovectores v1 = (0, 1) v2 = (1, 1)
Y 0 = 2X

La solución general es dada por

Z(t) = c1 e0t v1 + c2 e2t v2 = c1 v1 + c2 e2t v2

Por lo tanto el punto critico (X, Y ) = (0, 0) es un nodo impropio inestable, dado que repele a lo largo de cualquier
dirección menos v1 . De hecho cuando c2 = 0 la solución es constante.
Usando el Teorema de linealización, cuyas hypótesis están verificadas, podemos garantizar que el punto critico
nolineal P0 = (−1, −1) o bien es un nodo inestable o bien en un foco inestable para el sistema nolineal, de todas
formas se queda un repulsor, y repele tambien a lo largo de la dirección v1 = (0, 1). (Dibujar el lineal y las dos
posibilidades del nolineal)
• Dibujar el sistema nolineal global: ojo que sea coerente el dibujo.
Grupo de ma~
nana 1

Resumen de E.D.O. (2o de Matemáticas)

En estas páginas se enumeran y repasan los contenidos del primer examen. Obvia-
mente un resumen de la teorı́a, no es un sustituto de ella ni una razón para no realizar
problemas. Los ejemplos están tomados de las listas propuestas, indicándose el número
de hoja y ejercicio, ası́ 2.3 se refiere al ejercicio 3 de la hoja 2.

Capı́tulo 1
Después de la motivación fı́sica y algunas generalidades, se estudian diferentes méto-
dos para resolver EDO de primer orden. Los principales, con una breve descripción de
la idea subyacente, son:

Variables separables: Transformar la ecuación en una del tipo f (y)y 0 = g(x), que
se integra trivialmente (Ej. 1.5).

Ecuaciones homogéneas: Sustituyendo y por xz las ecuaciones homogéneas se con-

vierten en ecuaciones de variables separables (Ej. 1.9).

Ecuaciones exactas y factor integrante: Las ecuaciones de la forma ∂F ∂x

+ ∂F
∂y
y0 = 0
equivalen a F (x, y) = cte. Si no fueran de esta forma (exactas), quizá sı́ lo sean
tras multiplicar por cierta función (factor integrante), (Ej. 1.10b, 1.11).

Ecuaciones lineales: Se pueden ver como un caso especial de factores integrantes

(Ej. 1.15).

Cambio de variable: A veces un cambio de variable permite reducir una ecuación

a otra más sencilla (Ej. 1.17, 1.18).

Además de estos métodos, se vio cómo razonamientos geométricos (isoclinas) podı́an dar
información cualitativa de interés sobre las soluciones (Ej. 1.1); que algunas ecuaciones
de orden dos, o incluso mayor, se reducı́an al caso de orden uno (Ej. 1.16), y que con
una solución es posible a veces calcular otra (Ej. 3.15-16).

A través de sencillos modelos fı́sicos o geométricos, multitud de “problemas con

palabras” se reducen a problemas de EDO. En la asignatura Modelización I (y en parte
en Geometrı́a II) se profundizará sobre ello, pero ahora hay que saber hacer algunos
ejemplos básicos (Ej. 1.4, 1.6, 1.20-21, 1.23).

Capı́tulo 2
Tiene dos partes bien diferenciadas, de la segunda parte no se comenta aquı́ lo ex-
plicado tras la clase del 23/11/2005, que no es materia del primer examen.
Grupo de ma~
nana 2

2.A Convergencia puntual y uniforme

Digamos que {fn }∞ n=1 es una sucesión de funciones fn = fn (x) definidas en un con-
junto A (tı́picamente un intervalo).

Convergencia puntual: fn −→ f puntualmente, significa que para cualquier punto

x de A la sucesión de números fn (x) tiende al número f (x). Es una definición natu-
ral pero no muy conveniente para hacer análisis porque no respeta la continuidad
(Ej. 2.5) ni permite “tomar lı́mites” bajo integrales (Ej. 2.2) o derivadas.

Convergencia uniforme: fn −→ f uniformemente, significa intuitivamente que la

distancia entre las fráficas de fn y f tiende a cero. Su gran ventaja es que respe-
ta la continuidad, permite tomar lı́mites bajo integrales (propias) y, con ciertas
condiciones (Ej. 2.9), bajo derivadas.

La convergencia puntual o uniforme de series de funciones se reduce a la de sucesiones,

porque una serie no es más que la sucesión de sus sumas parciales, aunque hay algunas
técnicas que adquieren formas especı́ficas (por ejemplo el test de Weierstrass). En el
caso particular de series de potencias (digamos centradas en cero, para simplificar),
recordando el Cálculo I se puede ir más allá y probar que existe un radio de convergencia
R, que es computable con el criterio de la raı́z o del cociente (admitimos R = 0, ∞),
tal que la serie converge en el intervalo |x| < R y diverge en |x| > R. El caso x = ±R
hay que tratarlo directamente (Ej. 2.7 a,d,e,f). La convergencia es uniforme en cualquier
[a, b] incluido en |x| < R.
La convergencia de series de potencias no es ajena a las EDO, en realidad algunas
de ellas se resuelven sustituyendo formalmente una serie de potencias y calculando sus
coeficientes (Ej. 3.28, 3.29).

Desde el punto de vista topológico, se puede introducir en las funciones continuas

C([a, b]) la estructura de espacio métrico con la distancia d(f, g) = kf − gk∞ . La con-
vergencia en este espacio se traduce en la convergencia uniforme. La convergencia de
las sucesiones de Cauchy implica que se cumple el teorema del punto fijo (esto es útil
estudiando la existencia y unicidad).

La convergencia uniforme es un concepto delicado e importante. Se aconseja hacer

muchos ejercicios. Todos los de la hoja 2 son adecuados, excluyendo, si es necesario, 2.15
y 2.16 que son más avanzados, y 2.13 que es un anticipo de un tema que se tratará con
más detalle en EDAF.

2.B Ecuaciones lineales

Tres puntos iniciales importantes son:

Las EDO lineales de cualquier orden se pueden “reducir” a orden uno si se admiten
vectores. Esto es, equivalen a sistemas lineales Y~ 0 = A(x)Y~ + ~b(x) con A(x) ∈
Mn×n , ~b(x) ∈ Rn , lo cual es una gran ventaja teórica.
Grupo de ma~
nana 3

Bajo la hipótesis de continuidad de A(x) y ~b(x), los sistemas lineales tienen solución
única cuando se fija una condición inicial Y~ (x0 ) = Y~0 .

En el caso homogéneo (~b = ~0) la solución general (sin imponer un valor inicial) de
un sistema lineal da lugar a un espacio vectorial de dimensión n. Esto es, a partir
de n soluciones independientes (wronskiano 6= 0) se generan todas (Ej. 3.5-6). En
el caso no homogéneo el espacio es afı́n: basta sumar una solución particular al
espacio de soluciones de la homogénea.

Para resolver una EDO lineal con coeficientes constantes y (n) + an−1 y (n−1) + · · · + a0 y =
b(x) se siguen los siguientes pasos (Ej. 3.8, 3.10-14, 3.17-18):

Si r es raı́z del polinomio caracterı́stico xn + an−1 xn−1 + · · · + a0 entonces erx es

solución de la homogénea, y si r tiene multiciplidad k, xj erx también lo es con
j < k. Con ello se consiguen n soluciones independientes y por tanto la solución
general de la homogénea. En el caso de raı́ces complejas se toman partes reales e
imaginarias para presentar la solución en forma real.

Se halla una solución particular y se suma a la solución general de la homogénea.

Las soluciones particulares de la no homogénea se pueden hallar con dos métodos

cuya idea se describe aquı́:

Coeficientes indeterminados: Ajustando combinaciones lineales de funciones senci-

llas. Sólo es útil si ~b es de cierto tipo.

Variación de las constantes: Se supone que la solución particular buscada es cierta

perturbación de la solución de la homogénea.

Los sistemas homogéneos con coeficientes constantes tienen un tratamiento paralelo

al de las ecuaciones. Hay soluciones asociadas a cada raı́z (autovalor) λ del polinomio
caracterı́stico. Concretamente son de la forma eλx~v donde ~v es autovector (A~v = λ~v ). Si
la matriz no es diagonalizable, se usa la base de Jordan y hay que introducir un factor x
en el caso 2 × 2. El caso complejo no aporta nuevas dificultades. Tomando partes reales
e imaginarias se pueden hallar soluciones en forma real.

Los métodos de coeficientes indeterminados y de variación de las constantes también

se aplican en el caso de sistemas. El segundo se puede escribir como Φ(x)C ~ 0 (x) = ~b(x)
donde Φ es una matriz fundamental (las columnas son soluciones independientes de la
homogénea, su wronskiano es no nulo).

La introducción de exponenciales de matrices permite un atractivo tratamiento sinté-

tico de la teorı́a de sistemas lineales e independiente de la dimensión.

Si es necesario, uno deberı́a repasar algún tema de Álgebra Lineal antes de resolver
ejercicios de sistemas (Ej. 3.19-3.24).
17/11/05 1

Problemas de la hoja 3 de E.D.O. resueltos el

17/11/05 y algunos más

1. a) W 0 = (y1 y20 − y10 y2 )0 = y200 y1 − y100 y2 . Utilizando la ecuación

W 0 = (−P y20 − Qy2 )y1 − (−P y10 − Qy1 )y2 = P (y1 y20 − y10 y2 ) = −P w.
R
b) W 0 /W = −P ⇒ W = cte · e− P
que se anula sólo si la constante es nula.

2. Porque al despejar y 00 los coeficientes no son continuos en un entorno de cero.

3. a) Escribiendo y1 = y, y2 = y 0 , Y~ = (y1 , y2 )t , el sistema es

µ 0¶ µ ¶ µ ¶
~ 0 y y2 0 1 ~
Y = = = Y.
y 00 −by1 − ay2 −b −a
El polinomio caracterı́stico de la matriz es (−λ)(−a − λ) + b = λ2 + aλ + b.
b) Pol. caráct. ec.= (−1)n Pol. caráct. matriz. Se puede hacer una prueba por induc-
ción (no es inmediata, se deja como ejercicio).

5. a) Digamos x1 (a) = x2 (a) = 0. Sean cuales sean x01 (a), x02 (a); siempre se pueden
escoger µ, λ no simultáneamente nulos tales que λx01 (a)+µx02 (a) = 0. Entonces la función
x = λx01 + µx02 satisface la ecuación con las condiciones x(a) = x0 (a) = 0. Por la unicidad
debe ser la función idénticamente nula.
b) Lo mismo que en el apartado anterior pero ahora x01 (a) = x02 (a) = 0 y se escogen
µ, λ tales que λx1 (a) + µx2 (a) = 0.

6. a) Son linealmente independientes porque W (cos x − sen x, 2 sen x) = 2 6= 0 . Dos

vectores independientes en un espacio de dimensión 2, siempre forman una base.
b) Hallando la solución, se tiene y = cos x + sen x. La expresión y = 1 · (cos x −
sen x) + 1 · (2 sen x) implica que las coordenadas son (1, 1).

8. Las raı́ces del polinomio caracterı́stico λ2 − 4λ + 3 son λ = 1 y λ = 3, por tanto

la solución general es y = Aex + Be3x . Al imponer las condiciones iniciales, A + B = 6,
A + 3B = 10, lo que conduce a la solución y = 4ex + 2e3x .

9. Por el método con el que hallamos la soluciones, está claro que la solución tiende
a cero si y sólo si las raı́ces de λ2 + aλ + b tienen ambas parte real negativa.
Si las raı́ces son complejas, α ± iβ, entonces a = −2α, b = α2 + β 2 , por tanto en este
caso α < 0 equivale a que a y b sean positivos.
Si las raı́ces son reales, α1 y α2 , a = −α1 − α2 y b = α1 α2 , por tanto α1 , α2 < 0
equivale también en este caso a que a, b > 0.

10. El polinomio caracterı́stico es λ3 − 3λ2 + 3λ − 1 = (λ − 1)3 . Como 1 es una raı́z

triple, la solución general es y = Aex + Bxex + Cx2 ex . Haciendo las cuentas, se obtiene
¡ ¢ ¡ ¢
y 0 = A + B + (B + 2C)x + Cx2 ex , y 00 = A + 2B + 2C + (B + 4C)x + Cx2 ex .
17/11/05 2

Al imponer y(0) = 1 se tiene A = 1, de la misma forma, y 0 (0) = 2 lleva a B = 1 e

y 00 (0) = 3 a C = 0. En definitiva, la solución es y = (x + 1)ex .

11. El polinomio es P = λ4 + 4λ3 + 8λ2 + 8λ + 4. Si i − 1 es raı́z, −i − 1 también debe

serlo, por tanto P es divisible por (λ − (i − 1))(λ − (−i − 1)) = λ2 + 2λ + 2. Haciendo
la división, se obtiene P = (λ2 + 2λ + 2)2 . Entonces las raı́ces son −1 ± i, cada una de
ellas doble, y la solución general es
y = Ae−x cos x + Be−x sen x + Cxe−x cos x + Dxe−x sen x.

13. La solución general de la homogénea es yH = A + Bex + Cxex porque las raı́ces

del polinomio caracterı́stico son λ = 0 y λ = 1 (doble). El segundo miembro es x2 · e0·x
y como 0 es raı́z de multiplicidad uno, por el método de coeficientes indeterminados
debemos ensayar soluciones particulares de la forma yp = x(α + βx + γx2 ). Sustituyendo
6γ − 2(2β + 6γx) + (α + 2βx + 3γx2 ) = x2 ,
que igualando los términos de grados 0, 1 y 2, lleva a α = 6, β = 2 y γ = 1/3. La
solución general es pues
1
y = x3 + 2x2 + 6x + A + Bex + Cxex .
3

14. Tras la fórmula trigonométrica cos(x + α) = cos α cos x − sen α sen x, se podrı́a
aplicar el método de los coeficientes indeterminados en su forma habitual, pero es más
cómodo y natural trabajar directamente con cos(x + α) y sen(x + α) (a fin de cuentas,
cos x y sen x se pueden obtener como combinaciones lineales de ellos y forman una base
de las soluciones de la homogénea). Ensayando con yp = Ax cos(x + α) + Bx sen(x + α),
yp00 + yp = 2(−A sen(x + α) + B cos(x + α))
que al igualar a cos(x + α) implica A = 0, B = 1/2, que corresponde a la solución
particular yp = 21 x sen(x + α).

15. Hacemos el cambio y = xu. Derivando y 0 = u + xu0 , y 00 = 2u0 + xu00 que al

sustituir y simplificar conduce a
u00 2 − 4x2
(2 − 4x2 )u0 + (x − x3 )u00 = 0 ⇒ = .
u0 x3 − x
Descomponiendo en fraciones simples
u00 −2 1 1
0
= − −
u x x+1 x−1
que integrando (y eligiendo la constante de integración 0) implica u0 = x−2 (x2 − 1)−1 .
Una nueva descomposición en fracciones simples del segundo miembro e integración con
constante nula, lleva a
Z Z Z
1 1 1 1 1 1 1 ¯x − 1¯
u=− 2
+ − = + log ¯ ¯.
x 2 x−1 2 x+1 x 2 x+1
17/11/05 3

Multiplicando esto por x tenemos por tanto una nueva solución, ası́ pues la solución
general es µ ¶
x ¯x − 1¯
y = Ax + B 1 + log ¯ ¯ .
2 x+1

16. Buscamos una segunda solución de la forma t2 u. Sustituyendo en la ecuación,

después de algunos cálculos, se deduce tu00 + 5u0 = 0. Es fácil obtener de aquı́ la solución
u = t−4 que corresponde a x2 (t) = t−2 . La solución general es entonces x(t) = At2 +Bt−2 .
Las condiciones x(1) = 2, x0 (1) = 0 implican A = B = 1, por tato x(t) = t2 + t−2 es la
solución buscada.

17. La solución de la homogénea es yH = Aex + Be2x . Por el método de variación de

las constantes, para hallar una solución particular resolvemos
(
A0 ex + B 0 e2x = 0
A0 ex + B 0 2e2x = 1/(1 + e−x )

Resolviendo el sistema lineal, A0 = −1/(ex + R 1), B 0 = 1/(e2x + ex ). Se puede escribir

B 0 = e−x − 1/(ex + 1) y ası́ basta con hallar dx/(ex + 1), que se puede calcular con el
cambio de variable ex = t. Con ello,

A = log(ex + 1) − x, B = log(ex + 1) − x − e−x

que corresponde a la solución particular

yp = (log(ex + 1) − x)ex + (log(ex + 1) − x − e−x )e2x .

La solución general es y = yp + yH que operando es

y = (A − 1)ex + Be2x + (log(ex + 1) − x)(ex + e2x ).

Obviamente se puede sustituir A − 1 por A.

18. La solución de la homogénea es yH = A cos x+B sen x. Por el método de variación

de las constantes, consideramos
(
A0 cos x + B 0 sen x = 0
⇒ A0 = −f (x) sen x, B 0 = f (x) cos x.
−A0 sen x + B 0 cos x = f (x)

Por tanto la solución particular obtenida es

Z x Z x
yp = cos x f (s)(− sen s) ds + sen x f (s) cos s ds
0 0

que usando una sencilla fórmula trigonométrica da la función del enunciado.