Problema 4 Suponga que un virus zombi ataca a la humanidad.

Según la información disponible,

el virus ha infectado una poblacion aislada de 10 millones de habitantes, transformando a Z0 de
ellos en zombies. Al parecer, la cantidad de zombies Z(t) se incrementa solo por transmición directa
del virus a sujetos sanos. Por otro lado, se ha desarrollado una vacuna que revierte completamente
los efectos del virus. Bajo la hipótesis de que la suma de sujetos sanos y zombies es constante, se
tiene que un modelo razonable es

Z 0 (t) = αZ (10 − Z) − β (10 − Z) , Z(0) = Z0

donde α > 0 mide la habilidad de los zombies para reproducirse, β ≥ 0 mide la efectividad en el
uso de la vacuna, y todo está medido en millones. En lo que sigue, considere α > 0 como fijo y
β ≥ 0 como un parámetro.

(a) Encuentre los puntos crı́ticos de la ecuación.

(b) Clasifique cada uno de los puntos crı́ticos encontrados en función del valor de β.

(c) Dado distintos valores de β ≥ 0 determine los valores de Z0 (con 0 ≤ Z0 ≤ 10) para los cuales
la población de zombies tiende asintóticamente a 10.

Solución.
(a) Los puntos crı́ticos son soluciones de 0 = (αZ − β) (10 − Z) = f (Z). Luego, son Z1 = 10
y Z2 = β/α (evidentemente, es posible que Z1 = Z2 ). (1,5 puntos)
(b) f (Z) es una función cuadrática con coeficiente lider negativo. Si sus dos raices coinciden,
entonces f ≤ 0. En caso contrario, en la primera raiz f 0 es positivo, mientras que en la segunda es
negativo. Por lo tanto, si β = α10, Z1 = Z2 = 10 es un nodo (0,5 puntos). Si β > α10 entonces
Z1 = 10 es fuente, y Z2 = β/α es sumidero (1 punto). Si β < α10 entonces Z1 = 10 es sumidero,
y Z2 = β/α es fuente (1 punto).
(c) Es más fácil ver lo que ocurre haciendo diagramas de fase:

caso β < α10 : ← β/α → 10 ←

caso β = α10 : ← β/α = 10 ←
caso β > α10 : ← 10 → β/α ←

o un diagrama de bifurcación (hecho aquı́ con α = 1, pero el dibujo es, cualitativamente, el mismo
para todo α > 0):

Si Z0 = 10 entonces Z(t) = 10 y, trivialmente, se tiene lo deseado independiente del valor

de β ≥ 0 (0,5 puntos). El único otro caso posible es cuando β < α10 y Z0 ∈ (β/α, 10]. (1,5
puntos)
(Observación: no se pide que el estudiante haga diagrama de fase o de bifurcación.)

6
PONTIFICIA UNIVERSIDAD CATÓLICA DE CHILE
FACULTAD DE MATEMÁTICAS
DEPARTAMENTO DE MATEMÁTICA
TEMPORADA ACADÉMICA DE VERANO 2019

INTERROGACIÓN 2
MAT1640

Problema 1 La funcion y1 (x) = x es solución del problema homogéneo

x2 y 00 − x2 + 2x y 0 + (x + 2) y = 0.

Encuentre la solución general de

x2 y 00 − x2 + 2x y 0 + (x + 2) y = x3

en el intervalo (0, ∞).

Solución. Primero observemos que, en el intervalo indicado, el problema es equivalente a

   
00 2 0 1 2
y − 1+ y + + y = x.
x x x2

Por reducción de orden

Z R Z x 2
e (1+ x )dx
2
e x
y2 = y1 2
dx = y1 2
dx = xex
y1 x

es otra solución del problema homogéneo linealmente independiente de y1 . (2 puntos)

Usamos variación de parámetros para encontrar una solución particular yp = u1 y1 + u2 y2 .
Tenemos  
 x xe x 
W (y1 , y2 ) =  x
 2 x x x 2 x
x  = x e + xe − xe = x e (1 punto)
1 xe + e
por lo que
x2 e x x2
u01 = − = −1 , u 0
2 = = e−x (1 punto)
x2 ex x2 e x
y tomamos
u1 = −x , u2 = −e−x . (1 punto)
Entonces, la solución general del problema es

y(x) = c1 x + c2 xex − x2 − x. (1 punto)

1
Problema 2 Considere un sistema masa/resorte compuesto de una masa de 1 kilo y un resorte
con constante del resorte k = 2. Suponga que el medio en el que se encuentra el sistema ofrece una
fuerza de amortiguación (opuesta al movimiento) igual a dos veces la velocidad, y que se aplica
una fuerza externa f (t) = eαt sin(t), donde α ∈ R es un parámetro.

(a) Encuentra la solución general del problema cuando α = 1.

(b) Encuentra la solución general del problema cuando α = −1.

Por 1 punto bonus para el problema, responda la siguiente pregunta: ¿Para qué valores de α la
solución general es no acotada? Justifique.

Solución. (a) El sistema es descrito por la EDO

x00 + 2x0 + 2x = et sin(t).

La solución general del problema homogéneo es

xc = e−t (c1 cos(t) + c2 sin(t)) . (1 punto)

Buscamos una solución particular de la forma xp = et (A sin(t) + B cos(t)). Calculamos

x0p = et (A sin(t) + B cos(t)) + et (A cos(t) − B sin(t)) = et ((A − B) sin(t) + (A + B) cos(t))

x00p = et ((A − B) sin(t) + (A + B) cos(t)) + et ((A − B) cos(t) − (A + B) sin(t))
= et (−2B sin(t) + 2A cos(t))

x00p + 2x0 + 2x = et sin(t) (−2B + 2A − 2B + 2A) + et cos(t) (2A + 2A + 2B + 2B)

= et sin(t) (4A − 4B) + et cos(t) (4B + 4A) .

Luego, necesitamos
4A − 4B = 1 , 4B + 4A = 0 (1 punto)
de donde A = 1/8, B = −1/8. Entonces, la solución general es
1
x(t) = e−t (c1 cos(t) + c2 sin(t)) + et (sin(t) − cos(t)) . (1 punto)
8
(b) Como ya tenemos la solución del problema homogéneo, solo debemos encontrar una solución
particular. Como en este caso f (t) es solución del problema homogéneo, buscamos

xp = te−t (A sin(t) + B cos(t)) (1 punto)

2
Calculamos
 
x0p = e−t (A sin(t) + B cos(t)) + t −e−t (A sin(t) + B cos(t)) + e−t (A cos(t) − B sin(t))
= e−t (A sin(t) + B cos(t)) + te−t [sin(t) (−A − B) + cos(t) (A − B)]

x00p = −e−t (A sin(t) + B cos(t)) + e−t (A cos(t) − B sin(t))

+ e−t [sin(t) (−A − B) + cos(t) (A − B)] − te−t [sin(t) (−A − B) + cos(t) (A − B)]
+ te−t [cos(t) (−A − B) − sin(t) (A − B)]
= e−t [sin(t) (−2A − 2B) + cos(t) (2A − 2B)] + te−t [2B sin(t) − 2A cos(t)]

x00p + 2x0 + 2x = e−t sin(t) (t(2B − 2A − 2B + 2A) − 2A − 2B + 2A)

+ e−t cos(t) (t(−2A − 2B + 2A + 2B) − 2B + 2A + 2B)
= e−t sin(t) (−2B) + e−t cos(t) (2A) . (1 punto)

Luego, A = 0 y B = −1/2, y la solución general es

t
x(t) = e−t (c1 cos(t) + c2 sin(t)) − e−t cos(t). (1 punto)
2
Pregunta extra: La solución del problema homogéneo es siempre acotada. Sabemos por el
punto anterior que si α = −1 entonces la solución particular es acotada. Para todo otro valor
de α la solución particular es de la forma xp = eαt (A sin(t) + B cos(t)). Estas soluciones son no
acotadas si y solo si α > 0. Por lo tanto, la solución general es no acotada si y solo si α > 0.

3
PONTIFICIA UNIVERSIDAD CATÓLICA DE CHILE
FACULTAD DE MATEMÁTICAS
DEPARTAMENTO DE MATEMÁTICA
TEMPORADA ACADÉMICA DE VERANO 2020

INTERROGACIÓN 2 – PAUTA
MAT1640

Problema 1. Considere la EDO autónoma

du ueu − 9u
= .
dt eu

Encuentre los puntos crı́ticos, dibuje sus diagramas de fase, y clasifı́quelos. Indique además,
para cada uno de ellos, si es inestable o asintóticamente estable.

Solución Problema 1. Los punto crı́ticos son u1 = 0 [0,5 pts] y u2 = ln(9) [0,5 pts]. Claramente

(−∞, 0) (0, ln(9)) (ln(9), ∞)

u − + +
u
e −9 − − +
u0 + − +

por lo que u1 = 0 es un atractator [1 pts], ası́ntoticamente estable [1 pts], y u2 = ln(9) es un

repelente [1 pts], inestable [1 pts]. Se asigna [1 pts] por cada diagrama de fase correspondiente
(se puede dibujar un solo diagrama con ambos puntos crı́ticos).

1
Problema 2. Resuelva dos, y sólo dos, de los siguientes problemas:

a) Encuentre la solución general de u(6) + u(5) + u(4) + u000 = 0 en R.

b) Encuentre la solución general de x4 u00 + x3 u0 − 4x2 u = 1 en (0, ∞).

c) Encuentre una familia 2-paramétrica de soluciones para u2 u00 = u0 .

Solución Problema 2. (a) La ecuación auxiliar es

0 = m6 +m5 +m4 +m3 = m3 (m3 +m2 +m+1) = m3 (m+1)(m2 +1) = m3 (m+1)(m+i)(m−i) [1,5 pts]

por lo que la solución general es

u(x) = c1 + c2 x + c3 x2 + c4 e−x + c5 cos(x) + c6 sin(x). [1,5 pts]

(b) La EDO es de Cauchy-Euler, por lo que intentamos u = xm , de donde

0 = m(m − 1) + m − 4 = m2 − 4 [0,5 pts]

por lo que la solución del problema homogéneo asociado es

uc (x) = c1 x2 + c2 x−2 . [0,5 pts]

Buscamos solución particular con variación de parámetros. Tomamos u1 = x2 , u2 = x−2 , y

calculamos
 2     2 
 x x −2   0 x −2   x 0 
W =  
−3  = −4x
−1
, W1 =  −4 
−3  = −x
−6
, W2 =  
−4  = x
−2
2x −2x x −2x 2x x

([0,5 pts] por los tres cálculos en conjunto) por lo que

W1 1 1
c01 = = x−5 ⇒ c1 = − x−4
W 4 16
W 2 1 1
c02 = = − x−1 ⇒ c2 = − ln(x)
W 4 4
([0,5 pts] por ambos cálculos) y entonces una solución particular es
1 −2 1 −2 1
up (x) = − x − x ln(x) o simplemente up (x) = − x−2 ln(x)[0,5 pts]
16 4 4
pues x−2 es solución del problema homogéneo (no descontar si no hace esta simplificación). Luego,
la solución general es
1
u(x) = c1 x2 + c2 x−2 − x−2 ln(x).[0,5 pts]
4
0 00 dw 0
(c) Con w = u y u = du u obtenemos

dw 1 1
= 2 ⇒ w = − + c1 [1 pts]
du u u

2
y entonces Z
u
du = x + c2 [1 pts]
uc1 − 1
por lo que obtenemos las soluciones implı́citas
1
(u + ln(c1 u − 1)) = x + c2 . [1 pts]
c1

3
Problema 3.

a) Encuentre una EDO lineal no homogénea de segundo orden con coeficientes constantes tal que
u1 = 1 y u2 = e−x son soluciones del problema homogéneo asociado, y up = 21 x2 − x es solución
particular.

b) Se puede verificar (no lo haga) que u1 = √1 es solución de 4x2 u00 + 8x3 u0 + (4x2 − 3)u = 0.
x
Encuentre la solución general de

4x2 u00 + 8x3 u0 + (4x2 − 3)u = 4x7/2 .

Solución Problema 3. (a) Dado lo pedido, la ecuación auxiliar deber ser

0 = m(m + 1) = m2 + m [0,5 pts]

por lo que el operador diferencial debe ser D2 + D. Como

(D2 + D)up = 1 + x − 1 = x [0,5 pts]

la EDO buscada es
u00 + u0 = x. [0,5 pts]
(b) Sabemos que una segunda solución del problema homogéneo, L.I. de u1 , es
Z R Z
−1/2 e− 2xdx 2 1 2
u2 = x dx = x−1/2 xe−x dx = − x−1/2 e−x . [1 pts]
x−1 2
Para encontrar una solución particular usamos variación de parámetros. Calculamos
   
 x−1/2 − 12 x−1/2 e−x
2
 −x2 1 −2 1 −2 2

W =  1 −3/2 1 −3/2 −x2  x +1− x = e−x [0,5 pts]
1/2 −x2  = e
− x x e +x e 4 4
 2 4

 0 − 1 x−1/2 e−x
2
 1 −x2
W1 =  3/2 1 −3/2 2 −x2 
1/2 −x  = xe
2 [0,5 pts]
x 4
x e + x e 2
 
 x−1/2 0 
W2 =  1 −3/2 
3/2  = x [0,5 pts]
−2x x

de modo que
W1 1 1
c01 = = x ⇒ c1 = x2 [0,5 pts]
W 2 4
W2 2 1 2
c02 = = xex ⇒ c2 = ex [0,5 pts]
W 2
y una solución particular es
1 1
up (x) = x3/2 − x−1/2 [0,5 pts].
4 4

4
Como x−1/2 es solución del problema homogenéo, concluimos que la solución general es
2 1
u(x) = C1 x−1/2 + C2 x−1/2 e−x + x3/2 . [0,5 pts]
4

5
 Problema 4. Considere el caso de un sistema masa-resorte de masa 1 kilogramo que se encuentra
en equilibrio en un medio que genera una fuerza opuesta al movimiento la cual es igual a 10 veces
la velocidad en magnitud. Suponga la constante del resorte k es un parámetro modificable, y se
desea que en los tiempo t = 0 y t = 1 la masa se encuentre en su posición de equilibrio. Encuentre
los valores de k de modo que esto sea posible de lograr con una masa en movimiento. Para cada
uno de estos valores de k indique cuantas veces la masa en movimiento pasa por la posición de
equilibrio en el intervalo de tiempo (0, 1).

Solución Problema 4. El problema se modela con el problema de valores de frontera

x00 (t) + 10x0 (t) + kx(t) = 0 , x(0) = x(1) = 0 [1 pts]

y se nos pide encontrar soluciones no triviales en función del valor de k (el cual se puede asumir
positivo). Las soluciones de la ecuación auxiliar son
√
−10 ± 100 − 4k √
m= = −5 ± 25 − k [0,5 pts]
2
de modo que tenemos tres casos:

Caso 1: 25 − k > 0. En tal caso las soluciones son

√ √
x(t) = c1 e(−5+ 25−k)t
+ c2 e(−5− 25−k)t
[0,5 pts]

y de las condiciones de borde

c1 + c2 = 0
√ √
−5+ 25−k
c1 e + c2 e−5− 25−kt
= 0.

Se deduce fácilmente que c1 = c2 = 0 es la única solución, por lo que en este caso no es

posible obtener lo pedido [0,5 pts].

Caso 2: 25 − k = 0. Ahora
x(t) = c1 e−5t + c2 te−5t [0,5 pts]
y de las condiciones de borde

c1 =0
−5 −5
c1 e + c2 e = 0.

Nuevamente se deduce que c1 = c2 = 0 es la única solución, y por lo tanto no es posible

obtener lo pedido [0,5 pts].

Caso 3: 25 − k < 0. En este caso

 √ √ 
x(t) = e−5t c1 cos(t k − 25) + c2 sin(t k − 25) [0,5 pts]

6
 y de las condiciones de borde

c1 = 0
 √ √ 
e−5 c1 cos( k − 25) + c2 sin( k − 25) = 0.
√
Luego c1 = 0, pero la segunda ecuación se cumple para cualquier c2 provisto k − 25 es un
múltiplo de π. En consecuencia, se puede lograr lo pedido si k es alguno de los números

kn = 25 + (nπ)2 , n ∈ N [1 pts]

en cuyo caso el movimiento es descrito por

p
xn (t) = c2 e−5t sin(t kn − 25) = c2 e−5t sin(nπt). [0,5 pts]

De lo anterior observamos que si k = kn entonces la masa (en movimiento) pasa por la

posición de equilibrio n − 1 veces en el intervalo de tiempo (0, 1) [0,5 pts].

Tiempo: 2 horas.

7
 Interrogación 2
MAT 1640 - Ecuaciones Diferenciales

1. El modelo de Schaefer busca representar la evolución de una población de peces como producto
de su crecimiento logístico y su extracción a causa de la pesca. El modelo es el siguiente:
dy  y
=r 1− y − Ey, y ≥ 0,
dt K
donde todas las constantes E, r, K son números positivos. La tasa de extracción de peces
está dada por el término Ey.

a) Determine los valores de E para los cuales existen 2 puntos críticos.

b) Para los valores de E encontrados en la parte anterior, haga un diagrama de fase, cla-
sifique los puntos críticos y grafique las soluciones de equilibrio y otras curvas solución
alrededor de estas.
c) Para una condición inicial y(0) > 0, calcule lı́mt→∞ y(t). Luego, encuentre el valor de la
constante E que maximiza la tasa de extracción asintótica, es decir, la tasa de extracción
límite cuando t → ∞.

Solución.


a) Definiendo la función f (y) = r 1 − Ky y − Ey, los puntos críticos son aquellos valores
y ≥ 0 tales que f (y) = 0. Luego, como
 r 
f (y) = 0 ⇔ y r − E − y = 0
K

la función f (y) se anula para y = 0 y para y = K(1 − Er ). Concluimos que hay 2 puntos
críticos en la región y ≥ 0 si E < r.

Evaluación.
Asignar (0.5 ptos) por determinar correctamente que y = 0 es un punto crítico.
Asignar (0.5 ptos) por determinar correctamente que y = k(1 − E/r) es un punto crítico.
Asignar (1 ptos) por determinar la condición E < r o 0 < E < r. Este puntaje se asignará únicamente si se
determinó a y = k(1 − E/r) entre los puntos críticos. Si el estudiante encuentra otra condición para E porque
equivocó en determinar el punto crítico no cero, pero esta condición es consistente con su punto crítico encontrado
y la restricción y ≥ 0, asignar ((0.5 ptos)).

1
b) Para obtener el diagrama de fase analizamos los signos de f (y) entorno a los puntos
críticos. Tenemos que (recordando que r > 0)
(0, K(1 − Er )) (K(1 − Er ), ∞)
y + +
r
(r − E) − K
y + −
f (y) + −
El diagrama de fase es entonces:

Concluimos que y = 0 es un punto crítico inestable (fuente), mientras que y = K(1 − Er )

es asintóticamente estable (sumidero). El gráfico de las soluciones de equilibrio y otras
alrededor de estas es de la forma:

Evaluación. El siguiente puntaje se asignará si el estudiante determino correctamente los puntos de equilibrio.
Asignar (0.25 ptos) por determinar correctamente el signo de f (y) entre 0 y K(1 − E/r).
Asignar (0.25 ptos) por determinar correctamente el signo de f (y) entre K(1 − E/r) e ∞.
Asignar (0.75 ptos) por esbozar el diagrama de fase correctamente para y ≥ 0. Este diagrama puede ser como el
presentado en la pauta o uno con segmentos que emulan pendientes.
Asignar (0.25 ptos) por clasificar correctamente a y = 0 como un punto inestable o fuente.
Asignar (0.25 ptos) por clasificar correctamente a y = K(1 − E/r) como un punto estable o sumidero.
Asignar (0.75 ptos) por esbozar el gráfico de soluciones de equilibrio y otros alrededor de éstas correctamente para
y ≥ 0.

Evaluación. El siguiente puntaje se asignará si el estudiante equivocó en determinar el punto crítico diferente de cero,
pero para la condición dada para E el punto crítico es positivo.
Asignar (0.2 ptos) por determinar correctamente el signo de f (y) entre 0 y el otro punto crítico.
Asignar (0.2 ptos) por determinar correctamente el signo de f (y) entre el punto crítico no cero e ∞.
Asignar (0.5 ptos) por esbozar el diagrama de fase correctamente para y ≥ 0. Este diagrama puede ser como el
presentado en la pauta o uno con segmentos que emulan pendientes.
Asignar (0.2 ptos) por clasificar correctamente a y = 0 como un punto inestable o fuente.
Asignar (0.2 ptos) por clasificar correctamente a y = K(1 − E/r) como un punto estable o sumidero.
Asignar (0.5 ptos) por esbozar el gráfico de soluciones de equilibrio y otros alrededor de éstas correctamente para
y ≥ 0.

2
c) Para una condición inicial y(0) > 0, de acuerdo al diagrama de fase, tenemos que

E
lı́m y(t) = K(1 − ).
t→∞ r
Como de acuerdo con el enunciado del problema la tasa de extración de peces está dada
por el producto Ey, en el límite t → ∞
E
lı́m Ey(t) = K(1 − )E.
t→∞ r
Si consideramos la función g(E) = K(1 − E
r
)E, para maximizar su valor, debemos
encontrar el valor de E tal que
E K r
g 0 (E) = 0 ⇒ K(1 − )− E =0 ⇒ E= .
r r 2

Evaluación. Este puntaje se asignará tanto si el estudiante determinó correctamente los puntos críticos o a lo más se
equivocó en determinar el punto crítico no cero, pero que la condición sobre E es coherente con y ≥ 0.
Asignar (1 ptos) por obtener la función a optimizar g(E).
Asignar (0.5 ptos) por determinar el punto crítico de g.

3
2. Considere la ecuación diferencial lineal

sen2 (x)y 00 − 3 sen(x) cos(x)y 0 + (1 + 2 cos2 (x))y = 0, x ∈ (−π, π)

para la cual y1 (x) = sen(x) es una solución.

a) Encuentre una segunda solución y2 , linealmente independiente de y1 . Verique esto último

calculando el Wronskiano.
1 
Nota: tenga presente la identidad trigonométrica sen(a) cos(b) = sen(a + b) + sen(a − b)
2

b) Considere el problema de valor inicial en el intervalo (−π, π) y con condiciones iniciales

y(0) = y 0 (0) = 0. Encuentre dos soluciones al PVI y diga por qué no se contradice el
teorema de existencia y unicidad.

Solución.

a) Método 1: Utilizando el método de reducción de orden, buscamos una segunda solución

de la forma y2 (x) = v(x)y1 (x) = v(x) sen(x). Esto quiere decir que y2 debe satisfacer

sen2 (x)y200 − 3 sen(x) cos(x)y20 + (1 + 2 cos2 (x))y2 = 0

Como
y20 = v 0 sen(x) + v cos(x); y200 = v 00 sen(x) + 2v 0 cos(x) − v sen(x),
entonces

sen2 (x)(v 00 sen(x)+2v 0 cos(x)−v sen(x))−3 sen(x) cos(x)(v 0 sen(x)+v cos(x))+(1+2 cos2 (x))v sen(x) = 0

lo cual se simplifica:

sen3 (x)v 00 + (2 sen2 (x) cos(x) − 3 sen2 (x) cos(x))

+ (− sen3 (x) − 3 sen(x) cos2 (x) + sen(x) + 2 cos2 (x) sen(x)) v = 0;
| {z }
− sen3 (x) + sen(x) − sen(x) cos2 (x)
| {z }
2 2
sen(x) cos (x) − sen(x) cos (x)
| {z }
0

y finalmente obtenemos la ecuación de segundo orden reducible

cos(x) 0
v 00 − v = 0.
sen(x)
Tomando u = v 0 y resolviendo la ecuación para u mediante factor integrante, se tiene
que R cos(x)
u = Ce sen(x) dx = Celog(sen(x)) = C sen(x), x ∈ (0, π)
Entonces, Z Z
v= u(x)dx = C sen(x)dx = −C cos(x),

y por lo tanto y2 (x) = −C cos(x) sen(x) = − C2 sen(2x). Como la constante es arbitraria,

podemos simplemente tomar C = −2 y obtener

y2 (x) = sen(2x)

4
sin embargo cualquier función de la forma C sen(2x) con C 6= 0 también es válida como
solución.

Método 2. Usando la fórmula de Abel (para el wronskiano de soluciones), si y1 , y2 son

dos soluciones a la ecuación homogénea

cos(x) 0 (1 + 2 cos2 (x))

y 00 − 3 y + y = 0, en (0, π).
sen(x) sen3 (x)

por lo que deben satisfacer la igualdad

R cos(x)
y1 y20 − y10 y2 = W (y1 , y2 ) = Ce− (−3 sen(x) )dx
.

Usando que y1 (x) = sen(x), la igualdad se transforma en una ecuación de primer orden
para y2 :
sen(x)y20 − cos(x)y2 = C sen3 (x)
Resolvemos la ecuación
cos(x)
y20 − y2 = C sen2 (x)
sen(x)
R cos(x)
usando el factor integrante e− sen(x)
dx
= 1
sen(x)
y obtenemos que
Z
1 1
y2 =C sen2 dx = C cos(x)
sen(x) sen(x)

Por lo tanto, tomando por ejemplo C = 2, se llega a que

y2 (x) = 2 cos(x) sen(x) = sen(2x).

Evaluación. Si se procede por el Método 1.

Asignar (1 ptos) por remplazar correctamente y2 (x) = v(x) sin(x) en la ecuación diferencial y obtener una edo de
segundo orden para v reducible.
Asignar (0.5 ptos) por hacer el cambio de variable u = v 0 .
Asignar (1 ptos) por determinar la solución a la edo para u.
Asignar (0.5 ptos) por determinar v integrando v 0 = u.

Evaluación. Si se procede por el Método 2.

R
Asignar (0.5 ptos) por calcular correctamente
cos(x)
sin(x)
dx.
Asignar (1 ptos) por obtener la ecuación diferencial sen(x)y20 − cos(x)y2 = C sen3 (x), donde C puede aparecer como
constante genérica o con un valor asignado.
Asignar (0.5 ptos) por determinar un factor integrante para resolver la ecuación para y20 .
Asignar (1 ptos) por determinar una solución y2 . Este puntaje se asignará únicamente si hay un desarrollo que
justifica su elección y no por inspección.

5
b) Buscamos ahora resolver el PVI

cos(x) 0 (1 + 2 cos2 (x))

y 00 − 3 y + y = 0, y(0) = y 0 (0) = 0.
sen(x) sen3 (x)

Es claro que una solución está dada por la función constante y(x) = 0. Sin embargo, por
el principio de superposición, podemos encontrar otra solución a la ecuación anterior de
la forma
y(x) = C1 sen(x) + C2 sen(2x).
Imponiendo las condiciones iniciales, escogemos C1 , C2 tales que

y(0) = C1 sen(0) + C2 sen(0) = 0,

lo cual se satisface para todo par de constante C1 , C2 , y por otro lado, pedimos que

y 0 (0) = C1 cos(0) + 2C2 cos(0) = 0,

es decir C1 + 2C2 = 0. Concluimos que para cualquier constante C ∈ R, la función

C
y(x) = C sen(x) − sen(2x),
2
es solución al PVI.
Lo anterior no contradice el teorema de existencia y unicidad ya que la existencia de
multiples soluciones al PVI se explica por lo discontinuidad de los coeficientes de la
ecuación en x = 0.

Evaluación. Si el estudiante determinó correctamente y2 .

Asignar (0.5 ptos) por escribir la solución general de la edo.
Asignar (0.5 ptos) por obtener una ecuación para C1 y C2 como causa de evaluar la solución general en t = 0.
Asignar (0.5 ptos) por obtener una ecuación para C1 y C2 como causa de evaluar la derivada de la solución general
en t = 0.
Asignar (0.5 ptos) por determinar que la solución al PVI depende de una constante.
Asignar (1 ptos) por explicar por qué no se contradice el teorema de existencia y unicidad.

Evaluación. Si el estudiante encontró una y2 que no resuelve la edo.

Asignar (0.25 ptos) por escribir la solución general de la edo.
Asignar (0.25 ptos) por obtener una ecuación para C1 y C2 como causa de evaluar la solución general en t = 0.
Asignar (0.25 ptos) por obtener una ecuación para C1 y C2 como causa de evaluar la derivada de la solución general
en t = 0.

6
3. Dada la ecuación diferencial

x00 (t) − 2x0 (t) + 2x(t) = et sec(t)

a) Determine la solución general a la ecuación homogénea.

b) Use el método de variación de parámetros para determinar una solución particular a esta
ecuación.
c) Escriba la solución general de esta ecuación.

Solución.

a) La ecuación homogénea asociada es x00 (t) − 2x0 (t) + 2x(t) = 0 con ecuación característica
r2 − 2r + 2 = 0, cuyas raíces son r = 1 ± i. Por lo anterior, la solución general de la
ecuación homogénea está dada por

xh (t) = c1 et cos(t) + c2 et sin(t)

Evaluación.
Asignar (0.5 ptos) por determinar correctamente la ecuación característica.
Asignar (0.25 ptos) por cada raíz correctamente calculada de la ecuación característica. Total (0.5 ptos). Este
puntaje se asignará únicamente si se determinó correctamente la ecuación característica.
Asignar (1 pto) por cada una de las soluciones homogéneas linealmente independiente (x1 (t) = et cos(t) y x1 (t) =
et sin(t)). Total (2 ptos).
Si el estudiante determinó mal la ecuación característica pero calcula correctamente las raíces de su ecuación (aunque
no sean las correctas del ejercicio) y determina correctamente (de acuerdo a las raíces obtenidas) las dos soluciones
homogéneas linealmente independientes, asignar un total de (1 ptos) a la parte a).

7
b) Si x1 (t) = et cos(t) y x2 (t) = et sin(t) son las soluciones homogéneas l.i. asociadas a esta
ecuación, el método de variación de parámetros nos dice que una solución particular está
dada por xp (t) = u1 (t)x1 (t) + u2 (t)x2 (t) donde
Z Z
x2 (t)f (t)
u1 (t) = − dt = − tan(t) dt = ln(cos(t)) + k1
W (t) |{z}
0
Z Z
x1 (t)f (t)
u2 (t) = dt = 1 dt = t + k2
W (t) |{z}
0

siendo f (t) = et sec(t) y W (t) = e2t el Wronskiano de x1 y x2 . De este modo, una solución
particular es
xp (t) = ln(cos(t))et cos(t) + t et sin(t)

Evaluación. El siguiente puntaje se asignará únicamente si se calcularon correctamente las soluciones homogéneas en la
parte a).
Asignar (0.5 ptos) por calcular correctamente el Wronskiano: W = e2t o W = −e2t .
Z
x2 (t)f (t)
Asignar (0.5 ptos) por calcular correctamente . Este puntaje se asignará únicamente si calculó correc-
W (t) Z
x2 (t)f (t)
tamente el Wronskiano. Si el estudiante calculó mal el Wronskiano pero calcula correctamente para el
W (t)
Wronskiano obtenido, asignar (0.25 ptos)
Z
x1 (t)f (t)
Asignar (0.5 ptos) por calcular correctamente . Este puntaje se asignará únicamente si calculó correc-
W (t) Z
x1 (t)f (t)
tamente el Wronskiano. Si el estudiante calculó mal el Wronskiano pero calcula correctamente para el
W (t)
Wronskiano obtenido, asignar (0.25 ptos)
Asignar (0.5 ptos) por determinar correctamente la solución particular de la edo. Este puntaje se asignará Z única-
x2 (t)f (t)
mente si su respuesta es consecuencia de haber calculado correctamente el Wronskiano y las integrales
Z W (t)
x1 (t)f (t)
y . Si el estudiante equivocó en el cálculo de a lo más el Wronskiano y determina una expresión (libre
W (t)
de integrales) para la solución particular, asignar (0.25 ptos)

Evaluación. El siguiente puntaje se asignará a aquellos estudiantes que se equivocaron al calcular alguna (una o ambas)
de las soluciones homogéneas de la parte a).
Asignar (0.25 ptos) por calcular correctamente el Wronskiano.
Z
x2 (t)f (t)
Asignar (0.25 ptos) por calcular correctamente . Este puntaje se asignará únicamente si calculó co-
W (t)
rrectamente el Wronskiano y logra calcular la integral.
Z
x1 (t)f (t)
Asignar (0.25 ptos) por calcular correctamente . Este puntaje se asignará únicamente si calculó co-
W (t)
rrectamente el Wronskiano y logra calcular la integral.
Asignar (0.25 ptos) por determinar correctamente la solución particular de la edo. Este puntaje se asignará
Z única-
x2 (t)f (t)
mente si su respuesta es consecuencia de haber calculado correctamente el Wronskiano y las integrales
Z W (t)
x1 (t)f (t)
y .
W (t)

8
c) Finalmente, la solución general de la ecuación está dada por

x(t) = c1 et cos(t) + c2 et sin(t) + ln(cos(t))et cos(t) + t et sin(t)

Evaluación. El siguiente puntaje se asignará únicamente si se calcularon correctamente las soluciones homogéneas y
particular en las partes a) y b).
Asignar (1 ptos) por escribir la solución general de la ecuación diferencial.

Evaluación. El siguiente puntaje se asignará únicamente a aquellos estudiantes que en la parte a) obtienen la solución
homogénea habiendo cometido como único error la ecuación característica y/o en la parte b) como único error el cálculo
del Wronskiano.
Asignar (0.5 ptos) por escribir la solución general de la ecuación diferencial.

9
4. Resuelva la ecuación diferencial

x000 (t) − x0 (t) = −1 + 2 sin(t)

a) Determine la solución general a la ecuación homogénea.

b) Use el método de coeficientes indeterminados para determinar una solución particular a
esta ecuación.
Sugerencia. Se recomienda encontrar una solución particular para −1 y otra para 2 sin(t) y luego superposición de
soluciones.

c) Escriba la solución general de esta ecuación.

Solución.

a) La ecuación homogénea asociada es x000 (t) − x0 (t) = 0 con ecuación característica r3 − r =

0, cuyas raíces son r = 0, r = −1 y r = 1. Por lo anterior, la solución general de la
ecuación homogénea está dada por

xh (t) = c1 + c2 e−t + c3 et

Evaluación.
Asignar (0.5 ptos) por determinar correctamente la ecuación característica.
Asignar (0.25 ptos) por cada raíz correctamente calculada de la ecuación característica. Total (0.75 ptos). Este
puntaje se asignará únicamente si se determinó correctamente la ecuación característica.
Asignar (0.5 pto) por cada una de las soluciones homogéneas linealmente independiente (x1 (t) = 1, x2 (t) = et y
x3 (t) = e−t ). Total (1.5 ptos).
Si el estudiante determinó mal la ecuación característica pero calcula correctamente las raíces de su ecuación (aunque
no sean las correctas del ejercicio) y determina correctamente (de acuerdo a las raíces obtenidas) las tres soluciones
homogéneas linealmente independientes, asignar un total de (1 ptos) a la parte a).

10
b) Una solución particular a esta ecuación es de la forma xp (t) = u(t)+v(t) donde u000 −u0 =
−1 y v 000 − v 0 = 2 sin(t). Como −1 es parte de la solución homogénea, entonces un
candidato es u(t) = at, que luego de remplazar en la ecuación queda que −a = −1 y
por tanto a = 1. Por otro lado, sin(t) no es parte de la solución homogénea, entonces un
candidato es v(t) = b cos(t) + c sin(t), que luego de remplazar en la ecuación queda

2b sin(t) − 2c cos(t) = 2 sin(t)

de donde se deduce que b = 1 y c = 0. Por lo anterior, una solución particular es

xp (t) = t + cos(t)

Evaluación. El siguiente puntaje se asignará únicamente si se calcularon correctamente las soluciones homogéneas en la
parte a).
Asignar (1.25 ptos) por determinar a t o t + homogénea como parte de la solución particular.
Asignar (1 ptos) por determinar a cos(t) o cos(t) + homogénea como parte de la solución particular.

11
 c) La solución general de esta ecuación está dada por

x(t) = c1 + c2 et + c3 e−t + t + cos(t)

Evaluación. El siguiente puntaje se asignará únicamente si se calcularon correctamente las soluciones homogéneas y
particular en las partes a) y b).
Asignar (1 ptos) por escribir la solución general de la ecuación diferencial.

Evaluación. El siguiente puntaje se asignará únicamente a aquellos estudiantes que en la parte a) obtienen la solución
homogénea habiendo cometido como único error la ecuación característica.
Asignar (0.5 ptos) por escribir la solución general de la ecuación diferencial.

Duración: 120 min

No se permite el uso de calculadora

12
 Pontificia Universidad Católica De Chile
Facultad De Matemáticas
Departamento De Matemática
Primer Semestre 2022

Ecuaciones Diferenciales - MAT1640

Pauta Solución Interrogación 2

1. Resuelva la EDO homogénea

y (8) − 2y (4) + y = 0.

Solución
La ecuación caracterı́stica de la EDO es

r8 − 2r4 + 1 = 0.
⇐⇒ (r4 − 1)2 = 0
⇐⇒ (r2 − 1)2 (r2 + 1)2 = 0
⇐⇒ (r − 1)2 (r + 1)2 (r + i)2 (r + i)2 = 0
Las raı́ces de este polinomio son: r = 1, r = −1, r = i, r = −i todas de multiplicidad 2.
Por lo tanto, la solución de la EDO es de la forma
y(x) = c1 cos x + c2 sin x + c3 x cos x + c4 x sin x + c5 ex + c6 xex + c7 e−x + c8 xe−x .

Puntaje
• Determina ecuación caracterı́stica correcta: 2 puntos.
• Determina raı́ces correctas de ecuación caracterı́stica: 2 puntos.
• Determina solución general correcta y justificada: 2 puntos
• El puntaje de los ı́tem anteriores no se subdivide.
2. Determine el valor del parámetro a ∈ R tal que la solución del siguiente PVI
 00
x = x + ate−t
x(0) = 41 , x0 (0) = 0
es acotada cuando t → +∞.

Solución
La ecuación homogénea asociada tiene solución

xc (t) = c1 et + c2 e−t .

Para buscar una solución particular usamos el método de variación de parámetros:

Sean x1 (t) = et , x2 (t) = e−t . El wronskiano es W (x1 , x2 )(t) = 2. Luego:
Z Z
et −t −t e−t −a −t
xp (t) = − (e )(ate )dt + (et )(ate−t )dt = e (1 + t + 2t2 ).
2 2 8
Por lo tanto la solución general de la EDO es
a
x(t) = xc (t) + xp (t) = c1 et + c2 e−t − e−t (1 + t + 2t2 ).
8
Dadas las condiciones iniciales, obtenemos que
a+1 1
c1 = , c2 =
8 8
por lo que  
a+1 1 a
x(t) = et + e−t − e−t (1 + t + 2t2 ).
8 8 8
Finalmente limt→∞ x(t) resulta una función acotada si y sólo si a = −1 , ya que para
ese valor del parámetro se anula el primer término de la solución.
Puntaje
• Resuelve correctamente la EDO homogénea: 1 punto
• Determina correctamente solución particular: 2 puntos
• Determina correctamente solución general de la EDO: 1 punto
• Encuentre correctamente las constantes c1 y c2 : 2 puntos
• Determina correcta y justificadamente el valor de a: 1 punto
• El puntaje de los ı́tem anteriores no se subdivide.
3. Considere el siguiente PVI que modela el movimiento de un objeto unido a un resorte:
 00
x (t) + 5x0 (t) + 4x(t) = 0
x(0) = x0 (0) = 1
(a) Determine el tiempo t > 0 para el cual la posición x(t) del objeto llega a su valor
máximo.
(b) Si en el instante t se aplica una fuerza externa f (t) = e−t cos(t), ¿cuál será el
comportamiento del objeto cuando t → ∞? Justifique su respuesta.

Solución
(a) El sistema corresponde a un movimiento sobreamortiguado por lo que el punto
crı́tico de la solución que modela la posición corresponde al instante donde el ob-
jeto alcanza su máximo desplazamiento.

La ecuación caracterı́stica es r2 + 5r + 4 = 0 y sus raı́ces r = −1, r = −4. Dadas

las condiciones iniciales, la solución es
5 2
x(t) = e−t − e−4t .
3 3
Para encontrar el punto de máximo desplazamiento, resolvemos x0 (t) = 0:
 
0 5 −t 8 −4t −3t 1 5
x (t) = 0 ⇐⇒ − e + e = 0 ⇐⇒ 8e = 5 ⇐⇒ t = ln
3 3 3 8

1.0

0.8

0.6

0.4

0.2

1 2 3 4

(b) La dinámica del movimiento del objeto expuesto a la fuerza externa es gobernado
por la ecuación no homogénea x00 + 5x0 + 4x = e−t cos(t) la cual, por el método de
coeficientes indeterminados, tiene una solución general de la forma

x(t) = c1 e−t + c2 e−4t + Ae−t cos(t) + Be−t sin(t).

Ası́, cuando t → ∞, la solución x(t) → 0, independiente de los valores de c1 , c2 , A
y B, es decir, el objeto tiende a volver a la posición de origen.
Puntaje
• Determina la solución x(t) para la ecuación homogénea: 1 punto.
• Encuentra el valor de t buscado: 1 punto.
• Justifica correctamente que el valor encontrado de t es un punto de máximo: 1
punto
• Determina la forma que tiene la solución de la ecuación no homogénea o la encuentra
explı́citamente: 2 puntos.
• Calcula correctamente el lı́mite y determina el comportamiento del objeto: 1 punto.
• El puntaje de los ı́tem anteriores no se subdivide.
 MAT1640, semestre 2021/2
Pauta Interrogación 2
Profesores Alejandro Ramirez & Christian Sadel

Problema 1, Versión 1:
Encuentre la solución general de la ecuación
y 00 − 2y 0 + y = e3t

Solución La ecuación se puede escribir en la forma

(D − 1)2 y = e3t .
Luego tenemos que buscar las soluciones de

(D − 3)(D − 1)2 y = 0.
Pero la solución general de esta ecuación es

c1 et + c2 tet + c3 e3t .
Los primero dos términos corresponden a la solución general de la e versión
homogénea de la ecuación original. Por esa razón buscamos una solución
particular de la forma Ae3t . Luego

9A − 6A + A = 1,
1
lo que implica A = 4. Por lo tanto la solución general de la ecuación es

1
yg (t) = c1 et + c2 tet + e3t (6 pt)
4

Problema 1, Versión 2:
Encuentre la solución general de la ecuación
y 00 − 2y 0 + y = e4t

Solución La ecuación se puede escribir en la forma

(D − 1)2 y = e4t .
Luego tenemos que buscar las soluciones de

(D − 4)(D − 1)2 y = 0.
Pero la solución general de esta ecuación es
1
2

c1 et + c2 tet + c3 e4t .
Los primero dos términos corresponden a la solución general de la e versión
homogénea de la ecuación original. Por esa razón buscamos una solución
particular de la forma Ae4t . Luego

16A − 8A + A = 1,
1
lo que implica A = 9. Por lo tanto la solución general de la ecuación es
1
yg (t) = c1 et + c2 tet + e4t (6 pt).
9

Problema 2, Versión 1:

Para y = y(t) con t > 0 sea dado la EDO

(?) t2 y 00 + 7t y 0 + 9y = 0
a) Verifiqué, que y(t) = t−3 es una solución de (?).
b) Usando la solución dado en a), obtiene la solución general de (?).

Solución
a) (t−3 )0 = −3t−4 , (t−3 )00 = 12t−5 . En (?)
12t2 t−5 − 21t · t−4 + 9t−3 = 0 ⇒ 0t−3 = 0
ok, tenemos solución. (1 pt)

b) Metodo 1: Variación constante:

y(t) = u(t)t−3 . En (?) se obtiene

t2 (u00 t−3 − 6ut−4 ) + 7t(u0 t−3 ) = 0 ⇔ u00 + t−1 u0 = 0 (1 pt)
Factor integrante: eln(t) = t (1 pt).
0 = tu00 +u0 = (tu0 )0 ⇔ u0 (t) = C1 /t ⇔ u(t) = C1 ln(t)+C2 (2 pt)
y por es
y(u) = u(t)t−3 = C1 t−3 ln(t) + C2 t−3 ; . (1 pt)

Metodo 2: Identidad de Abel:

y1 = t−3 , Wronskiano W 0 + 7t−1 W = 0, factor integrante t7 (0,5 pt),

solución: W = Ct−7 . (1 pt).
Podemos usar C = 1 para obtener alguna solución y2 independiente de y1 .
t−7 = W = y1 y20 − y10 y2 = t−3 y20 + 3t−4 y2
⇔ y20 + 3t−1 y2 = t−4 (1 pt)
 3

Factor integrante t3 . (0,5 pt)

t−1 = t3 y20 + 3t2 y2 = (t3 y2 )0 ⇔ y2 = t−3 ln(t) + ct−3 (1 pt)
se puede usar c = 0. La solución general ce (?) es:
y(t) = c1 y1 + c2 y2 = c1 t−3 + c2 t−3 ln(t) . (1 pt)

Problema 2, Versión 2:

Para y = y(t) con t > 0 sea dado la EDO

(?) t2 y 00 + 5t y 0 + 4y = 0
a) Verifiqué, que y(t) = t−2 es una solución de (?).
b) Usando la solución dado en a), obtiene la solución general de (?).

Solución
a) (t−2 )0 = −2t−3 , (t−2 )00 = 6t−4 . En (?)
6t2 t−4 − 10t · t−3 + 4t−2 = 0 ⇒ 0t−2 = 0
ok, tenemos solución. (1 pt)

b) Metodo 1: Variación constante:

y(t) = u(t)t−2 . En (?) se obtiene

t2 (u00 t−2 − 4u0 t−3 ) + 5t(u0 t−2 ) = 0 ⇔ u00 + t−1 u0 = 0 (1 pt)
Factor integrante: eln(t) = t (1 pt).
0 = tu00 +u0 = (tu0 )0 ⇔ u0 (t) = C1 /t ⇔ u(t) = C1 ln(t)+C2 (2 pt)
y por es
y(u) = u(t)t−2 = C1 t−2 ln(t) + C2 t−2 ; . (1 pt)

Metodo 2: Identidad de Abel:

y1 = t−3 , Wronskiano W 0 + 5t−1 W = 0, factor integrante t5 (0,5 pt),

solución: W = Ct−5 . (1 pt).
Podemos usar C = 1 para obtener alguna solución y2 independiente de y1 .
t−5 = W = y1 y20 − y10 y2 = t−2 y20 + 2t−3 y2
⇔ y20 + 3t−1 y2 = t−3 (1 pt)
Factor integrante t2 . (0,5 pt)
t−1 = t2 y20 + 2ty2 = (t2 y2 )0 ⇔ y2 = t−2 ln(t) + ct−2 (1 pt)
se puede usar c = 0. La solución general de (?) es:
y(t) = c1 y1 + c2 y2 = c1 t−2 + c2 t−2 ln(t) . (1 pt)
4

Problema 3, Versión 1:
A un resorte colgado en el techo con constante de k = 1, 2N/m se pone una
masa de m = 300g. En una eje vertical orientado hacia abajo se pone el
punto y = 0 donde el sistema masa-resorte esta en equilibrio. A la masa
al tiempo t se pone una fuerza externa de 0, 3 sin(3t) N (t en segundos). Al
tiempo t = 0, la masa sea en el lugar de equilibrio (y = 0) con velocidad 0.
Describe el movimiento vertical de la masa por una EDO.
¿Donde esta la masa después de 10 segundos, (cuantos metros arriba o abajo
del punto de equilibrio)? y ¿cuál velocidad tiene?

Solución
Usamos tiempo t en segundos, y en metros, (y” en m/s2 ) fuerzas en Newton.
Fuerza de acceleración: N = kg m/s2 , ası́ necesitamos la masa en kg. m =
0, 3 kg. my 00 = −ky + f , o my 00 + ky = f . Asi, en unidades de Newton:
0, 3y 00 + 1, 2y = 0, 3 sin(3t)
cambiar masa correctamente (1 pt), obtener ecuación correcta: (0,5 pt.)
Ademas tenemos y(0) = 0, é, y 0 (0) = 0. (0,5 pt.) para usar condiciones
correctos
Para solucionarlo, lo ponemos en forma normal:
y 00 + 4y = sin(3t)
Solucion de ecuacion homogenea yc = c1 cos(2t) + c2 sin(2t). (1 pt.) No
aparece termino sin(3t).
Asi, se puede probar Y (t) = c3 cos(3t)+c4 sin(3t) para la solución particular.
sin(3t) = Y 00 + 4Y = −5c3 cos(3t) − 5c4 sin(3t)
Asi, c3 = 0, c4 = − 51 , Y (t) = − 15 sin(3t) (1 pt.) y tenemos
1
y(t) = − sin(3t) + c1 cos(2t) + c2 sin(2t)
5
3 3
0 = y(0) = c1 , 0 = y 0 (0) = − + 2c2 ⇒ c1 = 0, c2 =
5 10
y por eso
1 3
y(t) = − sin(3t) + sin(2t) (1 pt.)
5 10
y(10) = − 51 sin(30) + 10 3
sin(20) ≈ 0, 47 ;
y 0 (10) = − 35 cos(30) + 53 cos(20) ≈ 0, 15
Despues de 10 segundos, a masa esta aproximadamente 0,47 m = 47 cm
abajo del punto de equilibrio y tiene velocidad de aproximadamente 0,15
m/s = 15 cm/s (1 pt.), tb. esta ok dejar las respuestas con formulas, pero
debe tener la correcta unidad
 5

Problema 3, Versión 2:
A un resorte colgado en el techo con constante de k = 1, 8N/m se pone una
masa de m = 200g. En una eje vertical orientado hacia abajo se pone el
punto y = 0 donde el sistema masa-resorte esta en equilibrio. A la masa
al tiempo t se pone una fuerza externa de 0, 2 sin(3t) N (t en segundos). Al
tiempo t = 0, la masa sea en el lugar de equilibrio (y = 0) con velocidad 0.
Describe el movimiento vertical de la masa por una EDO.
¿Donde esta la masa después de 10 segundos, (cuantos metros arriba o abajo
del punto de equilibrio)? y ¿cuál velocidad tiene?

Solución
Usamos tiempo t en segundos, y en metros, (y” en m/s2 ) fuerzas en Newton.
Fuerza de acceleración: N = kg m/s2 , ası́ necesitamos la masa en kg. m =
0, 3 kg. my 00 = −ky + f , o my 00 + ky = f . Asi, en unidades de Newton:
0, 2y 00 + 1, 8y = 0, 2 sin(2t)
cambiar masa correctamente (1 pt), obtener ecuación correcta: (0,5 pt.)
Ademas tenemos y(0) = 0, é, y 0 (0) = 0. (0,5 pt.) para usar condiciones
correctos
Para solucionarlo, lo ponemos en forma normal:
y 00 + 9y = sin(2t)
Solucion de ecuacion homogenea yc = c1 cos(3t) + c2 sin(3t). (1 pt.) No
aparece termino sin(2t).
Asi, se puede probar Y (t) = c3 cos(2t)+c4 sin(2t) para la solución particular.
sin(2t) = Y 00 + 9Y = 5c3 cos(2t) + 5c4 sin(2t)
Asi, c3 = 0, c4 = 51 , Y (t) = 1
5 sin(2t) (1 pt.) y tenemos
1
y(t) = sin(2t) + c1 cos(3t) + c2 sin(3t)
5
2 2
0 = y(0) = c1 , 0 = y 0 (0) = + 3c2 ⇒ c1 = 0, c2 = −
5 15
y por eso
1 2
y(t) = sin(2t) − sin(3t) (1 pt.)
5 15
y(10) = 51 sin(20) − 152
sin(30) ≈ 0, 31 ;
y 0 (10) = 25 cos(20) − 52 cos(30) ≈ 0, 10
Despues de 10 segundos, a masa esta aproximadamente 0,31 m = 31 cm
abajo del punto de equilibrio y tiene velocidad de aproximadamente 0,10
m/s = 10 cm/s (1 pt.), tb. esta ok dejar las respuestas con formulas, pero
debe tener la correcta unidad
6

Problema 4, Versión 1:
Encuentre la solución general de la ecuación
3
y 00 + 4y =
sin(t)
Puede usar que
Z  
1  1 + sin(t) 

dt = ln  +C
cos(t) cos(t) 
Z  
1  1 + cos(t) 
dt = − ln  +C
sin(t) sin(t) 

Solución Hay dos soluciones de la versión homogénea de la ecuación:

y1 (t) = cos(2t) y2 (t) = sin(2t) (1 pt)

Por otra parte, el Wronskiano de estas soluciones es

W [y1 , y2 ](t) = 2 cos2 (2t) + 2 sin2 (2t) = 2. (1 pt)

Luego, una solución particular de la ecuación es

Z t 3
y2 (s) sin(s) Z t 3
y1 (s) sin(s)
yp (t) = −y1 (t) ds + y2 (t) ds
t0 2 t0 2
Z t Z t
3 sin(2s) 3 cos(2s)
= − cos(2t) ds + sin(2t) ds. (2 pt)
t0 2 sin(s) t0 2 sin(s)
Pero
Z t Z t Z t
3 sin(2s) 3 2 sin(s) cos(s)
ds = ds = 3 cos(s)ds = 3 sin(t) + c
t0 2 sin(s) 2 t0 sin(s) t0
(de hecho c = 3 sin(t0 ) pero el parte −cy1 (t) se absorba al parte de solución
homogénea). Y
Z t Z Z t Z
3 cos(2s) 3 t 1 − 2 sin2 (s) 3 t 1
ds = ds = −3 sin(s)ds + ds
t0 2 sin(s) 2 sin(s) 2 t0 sin(s)
 t0  t0
3  1 + cos(t) 
= 3 cos(t) − ln  + c̃
2 sin(t) 
Luego la solución general de la ecuación es

yg (t) =c1 cos(2t) + c2 sin(2t) − 3 cos(2t) sin(t)

 
3  1 + cos(t) 
+ 3 sin(2t) cos(t) − sin(2t) ln  . (2 pt)
2 sin(t) 
 7

Problema 4, Versión 2:
Encuentre la solución general de la ecuación
3
y 00 + 4y =
cos(t)
Puede usar que
Z  
1  1 + sin(t) 

dt = ln  +C
cos(t) cos(t) 
Z  
1  1 + cos(t) 
dt = − ln  +C
sin(t) sin(t) 

Solución Hay dos soluciones de la versión homogénea de la ecuación:

y1 (t) = cos(2t) y2 (t) = sin(2t). (1 pt)

Por otra parte, el Wronskiano de estas soluciones es

W [y1 , y2 ](t) = 2 cos2 (2t) + 2 sin2 (2t) = 2. (1 pt)

Luego, una solución particular de la ecuación es

Z t 3
y2 (s) cos(s) Z t 3
y1 (s) cos(s)
yp (t) = −y1 (t) ds + y2 (t) ds
t0 2 t0 2
Z t Z t
3 sin(2s) 3 cos(2s)
= − cos(2t) ds + sin(2t) ds. (2 pt)
t0 2 cos(s) t0 2 cos(s)
Pero
Z t Z t Z t
3 sin(2s) 3 2 sin(s) cos(s)
ds = ds = 3 sin(s)ds = −3 cos(t) + c
t0 2 cos(s) 2 t0 cos(s) t0
(de hecho c = 3 cos(t0 ) pero el parte −cy1 (t) se absorba al parte de solución
homogénea). Y
Z t Z Z t Z
3 cos(2s) 3 t 2 cos2 (2) − 1 3 t 1
ds = ds = 3 cos(s)ds − ds
t0 2 cos(s) 2 t0 cos(s) t0 2 t0 cos(s)
 
3  1 + sin(t) 
= 3 sin(t) − ln  + c̃ .
2 cos(t) 
Luego la solución general de la ecuación es

yg (t) =c1 cos(2t) + c2 sin(2t) + 3 cos(2t) cos(t)

 
3  1 + sin(t) 
+ 3 sin(2t) sin(t) − sin(2t) ln  . (2 pt)
2 cos(t)