MÉTODOS DINÁMICOS

Taller I: Ecuaciones Diferenciales

1. Ejercicios de calentamiento

i) Probar que

a) 𝑥 (𝑡 ) = 𝑡 + 𝐴𝑡 −1 , con 𝐴 constante, es solución de

𝑥
𝑥̇ + =2
𝑡
1
b) 𝑥 (𝑡 ) = 𝐴 cos(3𝑡 ) + 𝐵𝑠𝑒𝑛(3𝑡 ) + cos(𝑡), con 𝐴 y 𝐵constantes, es
8
solución de
𝑥̈ + 9𝑥 − cos(𝑡) = 0
1 5
c) 𝑥 (𝑡 ) = 𝑡 3 (1 + 𝑡 ) + 𝑡 −1 + 1, con 𝑥 (0) = , es solución de
4 2
𝑥 𝑥
𝑥̇ = − + 𝑡2 + 𝑡3
1+𝑡 𝑡
ii) La ecuación 𝑥̇ + 𝑏𝑥 + 𝑐 = 0, con 𝑥 (0) = 𝑥0 , tiene como solución a
𝑥 (𝑡 ) = 3𝑒 3𝑡 + 2. Encontrar los coeficientes 𝑏 y 𝑐 y el valor inicial 𝑥0

iii) Considere las siguientes ecuaciones diferenciales lineales:

6
a) 2𝑥̇ = −12𝑥 − 2𝑒 𝑡 , Con 𝑥 (0) =
7
2 2 ( )
b) 𝑥̇ = −𝑡 𝑥 + 5𝑡 , Con 𝑥 0 = 6
c) 𝑥̇ + 𝑐𝑜𝑠(𝑡 )2𝑥 = 𝑐𝑜𝑠(𝑡), Con 𝑥 (0) = 1

Para cada una de ellas:

- Encuentre la solución de la ecuación, 𝑥(𝑡)
- Encuentre lim 𝑥(𝑡). ¿Indica este resultado que el sistema converge o
𝑡→∞
diverge? Explique
- Grafique la solución utilizando las técnicas para graficar funciones que
conoce.

iv) Obtenga la solución general de las siguientes ecuaciones diferenciales

a) (6𝑡𝑥)𝑑𝑡 + (4𝑥 + 9𝑡 2 )𝑑𝑥 = 0

 b) (10 − 6𝑡 + 𝑒 −3𝑥 )𝑑𝑥 − 2𝑑𝑡 = 0
c) 𝑡𝑑𝑥 − 3𝑥𝑑𝑡 = 𝑡 2 𝑑𝑡
𝑑𝑥
d) 𝑡 + (1 + 𝑡 )𝑥 = 𝑒 𝑡
𝑑𝑡
e) (5𝑥 2 − 2𝑡 2 )𝑑𝑥 − 𝑥𝑡𝑑𝑡 = 0
f) 𝑥𝑑𝑡 + 3𝑡 2 𝑒 𝑡 𝑑𝑡 = 𝑑𝑥
𝑥 𝑡 1 2𝑡 2
g) 2 ( 3 + 2
) 𝑑𝑡 = ( 2 + ) 𝑑𝑥
𝑡 𝑥 𝑡 𝑥3
h) 𝑥𝑡 𝑑𝑥 + (2𝑥 2 + 3𝑡 2 − 20)𝑑𝑡 = 0
𝑑𝑥 2𝑥 𝑡
i) 3 + = 2
𝑑𝑡 𝑡+1 𝑥
𝑑𝑃
j) = 𝐾 − 𝑎(𝑃 + 𝑡)
𝑑𝑡
k) 𝑥 𝑑𝑡 − 𝑡 3 𝑑𝑥 = 0
3
1
l) 𝑥̈ = 2
𝑡
m) −𝑥𝑑𝑡 + (𝑡 + √𝑡𝑥)𝑑𝑥 = 0
𝑑𝑥 𝑥𝑡+𝑥
n) =
𝑑𝑡 𝑡+𝑡𝑥
𝑑𝑉 𝑎𝑉
o) =
𝑑𝑃 𝑃−𝑏

v) Resolver las siguientes ecuaciones

𝑡3
a) 𝑥̇ = , Con 𝑥 (1) = 1
𝑥3
√𝑡+1
b) 𝑥̇ = , Con 𝑥 (1) = 1
𝑥2
c) 𝑥̇ 𝑥 (𝑡 + 4) + 𝑥 + 2 = 0, Con 𝑥 (−3) = −1
d) 𝑡 2 𝑥̇ − 11𝑡𝑥 = 𝑥 5 , Con 𝑥 (1) = 1
2
1 𝑒𝑥
e) 𝑑𝑡 + 𝑑𝑥 = 0, Con 𝑥 (1) = 0
𝑥𝑡 𝑡 3 −1

vi) Considere las siguientes ecuaciones diferenciales

a) 𝑥̇ − 8𝑥 = 5𝑡 2 , Con 𝑥 (0) = 4
b) 𝑥̈ + 3𝑥̇ − 4𝑥 = 12, Con 𝑥 (0) = 4 y 𝑥̇ (0) = 2
c) 𝑥̈ + 8𝑥̇ + 16𝑥 = 0, Con 𝑥 (0) = 4 y 𝑥̇ (0) = 2
d) −𝑥̈ = 2𝑥̇ + 2𝑥, Con 𝑥 (0) = 1 y 𝑥̇ (0) = 0
e) 𝑥̈ + 2𝑥̇ + 10𝑥 = 4 cos(2𝑡)

- Encuentre la solución de la ecuación, 𝑥(𝑡)

- Encuentre lim 𝑥(𝑡). ¿Indica este resultado un modelo convergente o
𝑡→∞
divergente? Explique

vii) Encontrar la solución particular de las siguientes ecuaciones

diferenciales
 a) 𝑥̈ − 3𝑥̇ + 4𝑥 = −16𝑡 2 + 24𝑡 − 8
b) 𝑥̈ − 3𝑥̇ + 4𝑥 = 2𝑒 2𝑡
c) 𝑥̈ − 3𝑥̇ + 4𝑥 = 2𝑡𝑒 𝑡 − 𝑒 𝑡
d) 𝑥̈ − 3𝑥̇ + 4𝑥 = −16𝑡 2 + 24𝑡 − 8 + 2𝑒 2𝑡 + 2𝑡𝑒 𝑡 − 𝑒 𝑡
e) 𝑥̈ + 𝑥̇ + 3𝑥 = 𝑠𝑒𝑛(𝑡)
f) 𝑥̈ + 6𝑥̇ + 13𝑥 = 𝑒 −2𝑡 𝑐𝑜𝑠(2𝑡)

viii) Oscilaciones

a) Considere la ecuación diferencial (homogénea) de segunda orden,

dada por:
𝑥̈ + 𝛾𝑥̇ + 𝛿𝑥 = 0
Usted conoce lo siguiente:
La solución de dicha ecuación implica oscilaciones; la solución no
converge ni diverge en el tiempo, y, las oscilaciones tienen períodos que
repiten cada 5 años. Se pide encontrar los valores de 𝛾 y de 𝛿.
b) Considere la siguiente ecuación diferencial de segundo orden
𝑥̈ + 𝑏𝑥̇ + 4𝑥 = 𝑒 0.1𝑡
Se sabe que 𝑥 crece al 10% anual en el largo plazo y en el corto plazo
presenta oscilaciones de 4 años. Indique cuál es el valor de 𝑏 y si los ciclos
tienden a crecer o a desaparecer. Explique su respuesta
c) Considere una ecuación diferencial lineal homogénea de segundo
orden, con coeficientes constantes, que genera ciclos de 5 años, los
mismos que no convergen ni divergen. Encuentre la ecuación

ix) Realizar el diagrama de fase para cada una de las siguientes ecuaciones,
identificar los puntos de equilibrio y decir si son estables e inestables

a) 𝑥̇ = 𝑥 3 − 15𝑥 2 + 36𝑥
b) 𝑥̇ = 𝑒 𝑥 𝑠𝑒𝑛 𝑥
𝐾
c) 𝑥̇ = 2𝑥 ln ( )
𝑥
d) 𝑥̇ = 𝑥 2 (1 − 𝑥)
𝑥2
e) 𝑥̇ =
2+𝑥 2
f) 𝑥̇ = (2 − 𝑥)𝑠𝑒𝑛 𝑥
g) 𝑥̇ = 𝑒 −𝑥 𝑥(𝑥 − 2)
2. Wronskiano
Suponga que cada una de las funciones 𝑓1 (𝑥), 𝑓2 (𝑥), . . , 𝑓𝑛 (𝑥), tiene al menos
𝑛 − 1 derivadas. El determinante
𝑓1 𝑓2 … 𝑓𝑛
𝑓1′ 𝑓2′ … 𝑓𝑛′
( )
𝑊 𝑓1 , 𝑓2 , … , 𝑓𝑛 = | |
⋮ ⋮ ⋱ ⋮
(𝑛−1) (𝑛−1) (𝑛−1)
𝑓1 𝑓2 ... 𝑓𝑛
se denomina el Wronskiano de las funciones.
Ahora, sean 𝑥1 , 𝑥2 , … , 𝑥𝑛 , 𝑛 soluciones de la ecuación diferencial lineal
homogénea de n-ésimo orden
𝑎𝑛 (𝑡)𝑥 (𝑛) + 𝑎𝑛−1 (𝑡)𝑥 (𝑛−1) +. … … … . . +𝑎1 (𝑡)𝑥̇ + 𝑎0 (𝑡)𝑥 = 0
en el intervalo 𝐼. El conjunto de soluciones es linealmente independiente en
𝐼 si y solo si 𝑊(𝑥1 , 𝑥2 , … , 𝑥𝑛 ) ≠ 0 para toda 𝑡 en el intervalo.
Utilizando el Wronskiano, determine si las siguientes funciones forman un
conjunto de soluciones linealmente independientes de la ecuación diferencial
dada, en el intervalo que se indica.
a) 4𝑥̈ − 4𝑥̇ + 𝑥 = 0; 𝑒 𝑥/2 , 𝑥𝑒 𝑥/2 ; (−∞, ∞)
b) 𝑡 3 𝑥⃛ + 6𝑡 2 𝑥̈ + 4𝑡𝑥̇ − 4𝑥 = 0; 𝑥, 𝑥 −2 , 𝑥 −2 ln(𝑥); (0, ∞)
c) 𝑡 2 𝑥̈ + 𝑡𝑥̇ + 𝑥 = 0; cos(ln 𝑥), 𝑠𝑒𝑛(ln 𝑥); (0, ∞)

3. Deducción de una segunda solución.

Suponga la ecuación diferencia lineal homogénea de segundo orden
𝑎2 (𝑡)𝑥̈ + 𝑎1 (𝑡)𝑥̇ + 𝑎0 (𝑡)𝑥 = 0
Sea 𝑥1 una solución de esta ecuación. Lo que se busca es una segunda solución
de la ecuación 𝑥2 , que sea linealmente independiente de la primera; es decir
que cumpla la condición
𝑥2 (𝑡 )
= 𝑢(𝑡)
𝑥1 (𝑡 )
La idea es hallar un método para encontrar 𝑥2 (𝑡 ) = 𝑢(𝑡)𝑥1 (𝑡 ). Partimos
dividiendo la ecuación por 𝑎2 (𝑡), de donde
𝑥̈ + 𝑎(𝑡)𝑥̇ + 𝑏(𝑡 )𝑥 = 0
𝑎1 (𝑡) 𝑎0 (𝑡)
Con 𝑎(𝑡 ) = y 𝑏 (𝑡 ) = .
𝑎2 (𝑡) 𝑎2 (𝑡)

Definiendo 𝑥2 = 𝑢(𝑡)𝑥1 (𝑡 ), se tiene que

 𝑥̇ 2 = 𝑢𝑥̇ 1 + 𝑢̇ 𝑥1
y que
𝑥̈ 2 = 𝑢𝑥̈ 1 + 2𝑥̇ 1 𝑢̇ + 𝑥1 𝑢̈
Reemplazando en la ecuación diferencial, tenemos
𝑢𝑥̈ 1 + 2𝑥̇ 1 𝑢̇ + 𝑥1 𝑢̈ + 𝑎(𝑢𝑥̇ 1 + 𝑢̇ 𝑥1 ) + 𝑏(𝑢𝑥1 )
Donde por simplicidad se ha obviado el símbolo (𝑡 ). Reordenando términos
𝑢̈ 𝑥1 + 𝑢̇ (2𝑥̇ 1 + 𝑎𝑥1 ) + 𝑢(𝑥̈ 1 + 𝑎𝑥̇ 1 + 𝑏𝑥1 ) = 0
Dado que el tercer término es igual a cero, entonces tenemos que
𝑢̈ 𝑥1 + 𝑢̇ (2𝑥̇ 1 + 𝑎𝑥1 ) = 0
Que es una ecuación de segundo orden en 𝑢. Definamos 𝑤 = 𝑢̇ , entonces la
ecuación previa se puede escribir como
𝑤̇𝑥1 + 𝑤 (2𝑥̇ 1 + 𝑎𝑥1 ) = 0
Esta última ecuación, de primer orden en 𝑤, es separable.
𝑑𝑤 𝑥̇ 1
𝑥1 = − (2 + 𝑎) 𝑑𝑡
𝑤 𝑥1
Planteando la integral a cada término
1 𝑥̇ 1
∫ 𝑑𝑤 = −2 ∫ 𝑑𝑡 − ∫ 𝑎𝑑𝑡
𝑤 𝑥1
E integrando, obtenemos

ln 𝑤 = −2 ln 𝑥1 − ∫ 𝑎𝑑𝑡 + 𝑘

Donde 𝑘 es una constante. Reagrupando términos

ln(𝑤 𝑥12 ) = − ∫ 𝑎𝑑𝑡 + 𝑘

Obteniendo el antilogaritmo

𝑤𝑥12 = 𝐾1 𝑒 − ∫ 𝑎𝑑𝑡
Re expresando en términos de 𝑢

𝑒 − ∫ 𝑎𝑑𝑡
𝑢(𝑡) = 𝐾1 ∫ + 𝐾2
𝑥12
Sin pérdida de generalidad, hacemos que 𝐾1 = 1 y 𝐾2 = 0
 𝑒 − ∫ 𝑎𝑑𝑡
𝑢(𝑡) = ∫
𝑥12
Finalmente, obtenemos la segunda solución como

𝑒 − ∫ 𝑎𝑑𝑡
𝑥2 (𝑡 ) = 𝑥1 (𝑡 ) ∫
𝑥12
Se pide:
Usando la expresión previa, encuentre una segunda solución, dada la primera
𝑥1 (𝑡 )
a) 𝑥̈ + 2𝑥̇ + 𝑥 = 0; 𝑥1 = 𝑥𝑒 −𝑥
b) 𝑡 2 𝑥̈ − 7𝑡𝑥̇ + 16𝑥 = 0; 𝑥1 = 𝑥 4
c) 𝑡𝑥̈ + 𝑥̇ = 0; 𝑥1 = ln 𝑥
d) Derive la segunda solución linealmente independiente, en el caso de
la ecuación de segundo orden homogénea y autónoma con raíces
repetidas (caso II)

4. Aplicaciones

i) Este primer ejercicio aplica los conceptos de ecuaciones diferenciales a

un campo muy distinto al de la economía y las finanzas.
Sin embargo, resulta sumamente interesante y ayuda a esclarecer la amplia
variedad de usos de estas herramientas matemáticas. Imagine un recipiente
con agua caliente a 80° de temperatura, en donde se introduce una lámina
metálica que está a 50°. La temperatura del medio ambiente externo que
rodea al recipiente, que se asume constante, es de 0° (Ver Gráfica 1.)
Sea 𝑥(𝑡) la ecuación que describe la evolución de la temperatura del agua en
el tiempo t, y sea 𝑦(𝑡) la ecuación que describe la evolución de la temperatura
de la lámina en el momento t. Las ecuaciones que describen el movimiento
de la temperatura del agua y de la lámina son:
𝑥̇ = −𝑎(𝑥 − 𝑇𝐴 )
𝑦̇ = −0.2(𝑦 − 𝑥)
Donde 𝑎 es una constante y 𝑇𝐴 representa la temperatura del medio ambiente.
Usted posee, además, la siguiente información: 𝑥 (0) = 80; 𝑥 (10) = 29.43;
𝑦(0) = 50. Se pide:
a) Encuentre el valor de la constante 𝑎
b) Encuentre 𝑦(𝑡)
 c) ¿En qué tiempo la lámina de metal adquiere su temperatura
máxima?

ii) Suponga que el área (superficie) de una ciudad es proporcional al

tamaño de su población y que las ciudades son círculos perfectos.
Suponga también que la tasa de crecimiento de la población es de 6% anual.
Se pide que plantee la ecuación diferencial que rige la evolución del radio, r,
de la ciudad.
iii) Considere una piscina con 10000 litros de agua, a la que se le abre el
tapón. ¿En cuánto tiempo tendrá la piscina 4000 litros, si, por hora, se
le va un 5% del agua?

iv) Expresar mediante ecuaciones diferenciales los siguientes fenómenos

y cuando se los requiera, resolver la ecuación

a) La población de una ciudad crece a una velocidad proporcional a la

población existente en cada momento, así como la diferencia entre
600000 y la población
b) En una operación financiera, el aumento del capital depositado en
un banco, es directamente proporcional al tiempo trascurrido desde
que se colocó el capital
c) Suponga una elasticidad precio de la demanda dada por 𝜀 = −𝛽𝑃,
donde 𝛽 es una constante y 𝑃 es el precio del bien. Encuentre la
función de demanda 𝑄(𝑃)
d) Considere la siguiente elasticidad precio de la demanda, dada por:

2𝑃2 − 𝛽𝑃
𝜀=
𝑄
Donde 𝛽 es una constante y 𝑃 es el precio del bien. Encuentre la función
de demanda 𝑄(𝑃), sabiendo además que 𝑄(0) = 500

v) La población de peces en un lago es de 500 toneladas y crece a una tasa

del 4% anual. Si todos los años se capturan 40 toneladas ¿Cuál sería la
población de peces dentro de 5 años?

vi) Usted deposita $ 1000 en una cuenta bancaria a una tasa del 10% anual,
con capitalización continua. Decide retirar un monto fijo anual al final
de los próximos cuatro años. ¿Cuál es el monto máximo que puede
retirar cada año?
vii) Usted tiene 1000 para gastar en los próximos 5 años. Usted cree que lo
mejor que puede hacer es gastar un monto decreciente en el tiempo,
donde la tasa de decrecimiento en el consumo sea de 10% anual.
Suponga una tasa de interés del mercado de 12% anual.

a) Plantee la ecuación diferencial correspondiente que dé cuenta del

monto o cantidad de dinero disponible en t
b) Resuelva la ecuación planteada previamente

viii) Un tío suyo, que tiene 50 años, le comenta a un amigo que tiene USD
$30 000 ahorrados en su cuenta de ahorro, luego de muchos años en
que ha depositado todos los meses 10% de su sueldo, a una tasa del
0.5% mensual. Su sueldo es hoy USD $600 mensuales, y asegura que
todos los meses, desde que empezó a trabajar, le han subido el sueldo
en 0.2% mensual. ¿Le cree a su tío? Explique por qué. ¿A qué edad
empezó a trabajar? Resuelva el problema suponiendo tiempo continuo

ix) El Sr. López acaba de ser contratado en una empresa con un sueldo
mensual de $ 350. Su jefe ha prometido subirle el sueldo a una tasa del
7% mensual. Por otra parte, el Sr. López ha decidido partir con un
consumo de $ 300 al mes. Como su sueldo irá creciendo, él ha decidido
que también lo hará su consumo, solo que este último lo hará a una tasa
del 2% mensual. Finalmente, el Sr. López ha decidido depositar lo que
le sobre en una cuenta de ahorro al 6% mensual. Se pide:

a) Plantee sin resolver, en tiempo continuo, la ecuación diferencial que

rige la evolución del saldo en la cuenta de ahorro a través del
tiempo. Justifique su respuesta, explicitando cualquier supuesto
adicional que considere necesario. Nota: La ecuación debe incluir
sólo la variable “saldo en la cuenta de ahorro” y sus derivadas
b) Resuelva la ecuación diferencial planteada previamente.

x) Sea 𝑤 la riqueza y 𝑢(𝑤) la utilidad proporcionada por la misma.

Supondremos que la función de utilidad cumple con las características
usuales, a saber 𝑢´(𝑤) > 0 y 𝑢´´(𝑤) < 0. El coeficiente de aversión
relativa al riesgo constante se puede representar por

𝑢´´(𝑤)𝑤
𝑠=−
𝑢´(𝑤)
 Y se define como la variación porcentual en la utilidad marginal dados
cambios unitarios en la riqueza. Encuentre 𝑢(𝑤)

xi) Sea 𝑝 el precio de un activo, 𝑝𝑒 el precio esperado del mismo, 𝑑 los

dividendos constantes y 𝑟 la tasa de un instrumento alternativo, por los
bonos del tesoro de USA, también constante. Asumamos que se cumple
con la siguiente condición de no arbitraje
𝑝𝑟 = 𝑑 + 𝑝̇ 𝑒
Por otro lado, las expectativas acerca del precio esperado evolucionan de
acuerdo con
𝑝̇ 𝑒 = 𝛼(𝑝 − 𝑝𝑒 )
Donde 𝑟 > 𝛼.
a) Considerando las dos ecuaciones previas, obtener una ecuación
diferencial para 𝑝̇ 𝑒 y resolverla suponiendo que 𝑝𝑒 (0) = 𝑝0𝑒
∞ 𝑑
b) Probar que ∫0 𝑑𝑒 −𝑟𝑡 𝑑𝑡 = , que es simplemente el valor presente
𝑟
de los dividendos futuros del activo descontados a la tasa 𝑟
𝑑 𝑑
c) Probar que se cumple lim 𝑝𝑒 (𝑡) = , donde = 𝑝∗
𝑡→∞ 𝑟 𝑟

xii) Modelo de empleo de Haavelmo. Sea 𝐾 el nivel de capital y 𝐿 el nivel

de empleo. La función de producción está dada por 𝑌 = 𝐾 𝛾 𝐿1−𝛾 . Se
supone además que el crecimiento en la tasa de empleo está dado por

𝐿̇ 𝐿
=𝛼−𝛽
𝐿 𝑌
Donde 𝛼, 𝛽 > 0 y 0 < 𝛾 < 1.

a) Demostrar que bajo las hipótesis anteriores la ecuación para el

empleo puede escribirse como
𝐿1+𝛾
̇𝐿 = 𝛼𝐿 − 𝛽
𝐾𝛾
b) Resolver la ecuación previa si se supone que al inicio 𝐿(0) = 𝐿0
c) Encontrar lim 𝐿(𝑡)
𝑡→∞

xiii) Suponga que la población de un país evoluciona de acuerdo con la

ecuación
𝑃̇ = (𝛼 − 𝛽)𝑃 − 𝐸
Donde 𝛼 > 0 representa la tasa de natalidad, 𝛽 > 0 la tasa de mortalidad y 𝐸
representa un número fijo de personas que emigran cada año.
Resolver para 𝑃(𝑡) dado 𝑃(0) = 𝑃0 , y analizar qué pasa con lim 𝑃(𝑡) para
𝑡→∞
distintos valores de los parámetros

xiv) La población del suburbio de una gran ciudad está descrita por el
siguiente modelo
𝑃̇ = 𝑃(10−1 − 10−7 𝑃)
Con 𝑃(0) = 5000, donde 𝑡 se expresa en meses. ¿Cuál es el valor límite de la
población? ¿Cuánto tardará la población en alcanzar el 85% de ese valor
límite?
xv) Modelo de crecimiento logístico con migración. El modelo es el
siguiente
𝑃
𝑃̇ = 𝑟𝑃 (1 − ) − 𝐸
𝐶
Con 𝑃(0) = 𝑃0 . Se supone que 𝐸, 𝑟, 𝐶 > 0 y que 𝐶 > 4𝐸/𝑟.
1
a) Se define 𝛼 = √𝑟 2 − 4𝐸𝑟/𝐶 y 𝐿 = (𝐶 + √𝐶 2 − 4𝐸𝐶/𝑟). Sea
2
𝑟
𝑦(𝑡 ) = (𝑃(𝑡 ) − 𝐿). Demostrar que la ecuación diferencial en
𝐶
términos de 𝑦 es de la forma 𝑦̇ + 𝛼𝑦 = −𝑦 2
b) Resolver la ecuación diferencial para 𝑦(𝑡 ) y luego encontrar la
solución correspondiente para 𝑃(𝑡 )
c) Demostrar que tanto 𝐿 como 𝐶 − 𝐿 son puntos de equilibrio

5. Métodos Numéricos

i) Usted acaba de adquirir un bosque que cuenta con 18 mil m 3 de una

especia forestal.
Su problema es determinar cuanta madera (en miles de m3) tendrá en cuatro
años más. Un experto en el tema le indica que, para la especia forestal en
cuestión, el crecimiento (en miles de m3) se rige por la siguiente ecuación
diferencial
1000𝑋 4
𝑋̇ =
𝑒𝑋
Use el método de Euler y el método de Runge-Kutta para determinar el
volumen de madera en cuatro años más. Utilice ℎ = 0.4
ii) Resuelva mediante los métodos de Euler mejorado y Runge-Kutta, las
siguientes ecuaciones diferenciales:

a) 𝑥̇ = [ln 𝑥]2 ; 𝑥 (0) = 9

b) 𝑥̇ + 𝑡(ln 𝑥 ) + 𝑥 2 = 8𝑡 2 ; 𝑥 (0) = 1
c) 𝑥̇ = 𝑒 −𝑥 𝑥(𝑥 − 2); 𝑥 (0) = 1
Utilice, cada caso, un valor de ℎ = 0.25

6. Con ayuda de la computadora

La siguiente ecuación, adaptada de la logística, modela la población ajustada
por la inmigración
𝑃̇ = 𝑃(𝑎 − 𝑏𝑃) + 𝑐𝑒 −𝑘𝑃
Donde 𝑐 y 𝑘 son parámetros, mayores que cero.
a) ¿Cómo se explica la inclusión del término adicional en el modelo?
b) Suponga el caso particular
𝑃̇ = 𝑃(1 − 𝑃) + 0.3𝑒 −𝑃
Encuentre manualmente los equilibrios del sistema, utilizando el
método de Newton-Rhapson. Luego utilice un programa
informático y verifique sus resultados. ¿Cuáles de los equilibrios
encontrados tienen coherencia, desde el punto de vista del
fenómeno modelado?
c) Usando Wolfrand Alpha, resuelva numéricamente la ecuación, para
𝑃(0) = 0.5, 𝑃(0) = 1 y 𝑃(0) = 1.5. Considere un ℎ = 0.25. Luego
grafique en Excel las soluciones encontradas, en el intervalo 𝑡 ∈
(0; 3).
d) Utilizando un programa de graficación, grafique el campo vectorial
correspondiente
e) Finalmente, grafique el diagrama de fase, utilizando un programa de
graficación. Explique como puede utilizar esa gráfica para
determinar si los equilibrios encontrados previamente son
atractores.