Solucionario de Lomeli
Solucionario de Lomeli
Solucionario de Lomeli
1 de septiembre de 2004
Índice General
5 Análisis cualitativo 15
9 Optimización estática 22
1
Nota para el lector
La presente es una versión preliminar de las soluciones del libro Métodos dinámicos
en economía.
Estamos conscientes de que, a pesar del esfuerzo y cuidado puestos en este tra-
bajo, es probable que existan errores involuntarios en las respuestas que aquí pre-
sentamos. Agradeceremos sus comentarios y correcciones al presente documento
a las siguientes direcciones electrónicas:
lomeli@itam.mx
rumbos@itam.mx
http://matematicas.itam.mx
Los autores
2
Capı́tulo 2
Ecuaciones diferenciales lineales
2.2 b = −3, c = 6, x0 = 5.
2.3 α = 3, β = − 19 , A = 1 1
18 , B
= 18 . Por lo tanto la solución para las condiciones
1 1 1
iniciales dadas es x(t) = − + e3t + e−3t .
9 18 18
2.4 α = 0, β = 7. Por lo tanto la solución para las condiciones dadas se puede
escribir como y(v) = 7e−v sin v.
2.8 P(t) = P0 e(α−β)t . Si α > β, lim P(t) = ∞, es decir que P crece indefinidamente.
t→∞
Si α = β, lim P(t) = P0 , es decir que P es siempre constante. Si α < β,
t→∞
lim → ∞P(t) = 0, es decir que P se extingue.
t
E E at E E
2.9 P(t) = + P0 − e . Si P0 = , entonces P(t) = , por lo tanto lim P(t) =
a a a a t→∞
E E
es decir, la población es constante. Si P0 > , entonces lim P(t) = ∞, es
a a t→∞
3
4
E
decir la población es creciente. Si P0 < , entonces lim P(t) = −∞, es decir
a t→∞
la población es decreciente.
T
2.10 Factor de integración µ(t) = e t r(s)ds . Interpretación: Y(t) es la inversión, B(t)
t
YT
es el precio del bono y + δ(s )ds es la cantidad de bonos que se tienen
BT T
en la inversión.
1
2.11 a) r(t) = r0 − .
t+1
r0 (t−T) T+1
b) B(t) = e .
t+1
c) δ(t) = −er0 (T−t) . Si δ < 0 tenemos retiros.
1 r0 (T−t)
d) Z(T) = e − 1 + 1.
r0
T + 1 r0 (t−T) 1 −r0 (t−T)
e) Y(t) = e 1+ e −1 = B(t)Z(t). Simplificando
t + 1 r0
T+1 1 1
Y(t) = 1− er0 (t−T) + .
t+1 r0 r0
1
2.14 a) x(t) = + ce−2 sin t .
2
1
b) x(t) = + ce−t .
2
2
t3
c) x(t) = 5 + e− 3 .
5
1
d) x(t) = − et + e−6t .
7
1 2 −u3
e) y(u) = + e .
3 3
e αr dα
2.15 a) ṗ + pe = . Resolviendo encontramos que
r−α r−α
d d −( αr )t
p (t) = + p0 −
e e
e r−α
r r
.
∞ b
d
−rt d
b) de dt = lim de−rt dt = − lim e−rb − 1 = .
b→∞ 0 r b→∞ r
0
d d αr d
c) lim pe (t) = + p0e − lim e−( r−α )t = = p∗ ya queα, r > 0 y r > α.
t→∞ r r t→∞ r
d α
d) Sustituyendo (2.29) en (2.28) se tiene que p = + (p − pe ) . Por lo tan-
r r
∗α
to p = p (p − p ) . Ahora bien, usando la solución para pe se obtiene
e
r
que
r ∗ α
p(t) = p − pe .
r−α r−α
r α r
Por lo tanto lim p(t) = p∗ − y lim pe = p∗ −
t→∞ r − α r − α t→∞ r − α
α
p∗ = p∗ .
r−α
αr αr
e) p(t) = p − ∗
(p0e − p∗ ) e−( r−α )t , con r > α. Además
r−α
αr
pe (t) = p∗ + (p0e − p∗ ) e−( r−α )t ,
con r > α.
dv dv
2.16 Sea v = ln y, entonces ev = y y y = ev . Sustituyendo ev + P(t)ev =
dt dt
Q(t)ev v. Por lo tanto v − Q(t)v = −P(t).
t3 c
2.17 Sea v = ln y, resolviendo para v se obtiene v(t) = − + . Como y(t) = ev(t)
4 t
t3
entonces y(t) = e− 4 + t .
c
2.20 a) x(t) = et .
√ √
5 23 13 23
b) x(t) = e 4t 3 cos t + √ sin t .
4 23 4
c) x(t) = e−t [cos t + sin t] .
d) x(t) = (1 − 3t)e3t .
√ √
2.21 a) x(t) = k1 e(1+ 2) t
+ k2 e(1− 2) t
− 7.
b) x(t) = c1 cos t − c2 sin t + 1.
√ √
5 23 23
c) x(t) = e 4 t c1 cos t − c2 sin t + 3.
4 4
d) x(t) = A + Be3t − 4t.
1
e) x(t) = c1 e−t + c2 e2t − .
2
1
f) x(t) = c1 e−3t + c2 te−3t + .
9
β
2.22 p(t) = m̄ + e− 2 t (A cos δt − B sin δt) , con β > 0, δ 0.
β
e− 2 t −βA βB
u(t) = ū − − Bδ cos δt + − Aδ sin δt .
γ 2 2
Además lim p(t) = m̄ y lim u(t) = ū, lo que quiere decir que se satisface el
t→∞ t→∞
mismo comportamiento asintótico que en el caso β > 4αγ.
t
2.23 a) x(t) = c1 cos 2t − c2 sin 2t − cos 2t.
4
14
b) x(t) = c1 e−t + c2 e3t − 3t2 + 4t − .
3
4
c) x(t) = c1 e−t + c2 e2t − te−t .
3
1
d) x(t) = c1 e3t + c2 e−t − et − cos t + 2 sin t.
2
1 4
e) x(t) = e−3t (c1 cos 2t − c2 sin 2t) + e−2t cos 2t + sin 2t .
17 17
27
2
3.2 a) y(x) = sin−1 2
.
x +1
b) y(x) − 2 ln |y(x) + 2| = − ln |x + y| − 1.
1 t 1 3t
c) y(x) = − e + e .
2 2
N∗
3.3 a) N(t) = para N ∗ = N0 . Si N ∗ = N0 entonces N(t) =
N∗
1+ N0 −1 e N ∗ kt
∗
N .
b) lim N(t) = N ∗ , es decir que el número de personas que habrá oído el
t→∞
rumor cuando t sea muy grande tenderá al número total de personas
del pueblito.
s
3.4 Sea w = k1−α . Resolviendo se obtiene w(t) = ce−(1−α)(n+δ)t + . Por lo tanto
n+δ
1−α
1
1
−(1−α)(n+δ)t s
la solución para k es de la formak(t) = w 1−α = ce + .
n+δ
7
8
1−α
1
s
Además lim k(t) = = k∗ .
t→∞ n+δ
L̇ L L 1
3.5 a) Sea Υ = K γ L1−γ . Entonces = α − β = α − β γ 1−γ = α − β γ 1−γ =
L Υ K L K L
Lγ Lγ+1
α − β γ . Por lo tanto L̇ = αL − β γ , donde K es constante.
K K
−γ 1 β β
b) Sea w = L . Resolviendo se obtiene w(t) = γ − γ
e−αγt + .
L0 αK αK γ
1
− γ1 1 β −αγt β −γ
Por lo tanto L(t) = w(t) = γ − e + .
L0 αK γ αK γ
1 γ1
β −γ α
c) lim L(t) = γ
. Por lo tanto lim L(t) = K.
t→∞ αK t→∞ β
1
3.6 a) Sea w = y1−n . Su solución está dada por w(t) = keαt − . Por lo tanto
α
1 C
y(t) = αt . Como P = y + L, entonces
ke − 1/α r
1
P(t) = + L.
1
P0 −L + r
αC eαt − r
αC
L, P0 = C − L
b) lim P(t) = .
t→∞ P0 , P0 = C − L
1
3.7 a) x(t) = .
2t − 2 + ce−t
1 1
b) Sea w = 2 , cuya solución es w(x) = x + + ce2x . Por lo tanto y(x) =
y 2
− 21
1
± x + + ce2x .
2
1 x+c x
c) Sea w = , cuya solución es w(x) = . Por lo tanto y(x) = .
y x x+c
1
d) Sea w = 3 , cuya solución es w(x) = x3 2x3 + c . Por lo tanto y(x) =
y
1
1
.
x [2x3 + c] 3
3.8 a) Sea w = x−6 , entonces la tenemos la solución w(t) = 1 + ce6t . Por lo tanto
x(t) = 1.
4 c
b) Sea w = x−4 , entonces la tenemos la solución w(t) = − + 44 . Por lo
43t t
1
tanto x(t) = .
4 47
43t 44 − 4
43t
9
1 c
c) Sea w = y−2 , entonces la tenemos la solución w(t) = + √ . Por lo
√ t t
tanto y(t) = t, con t > 0.
1
b) ṗ = (p − m0 − µt) . Cuya solución, para las condiciones iniciales da-
λ
das, es:
t
p(t) = (m0 + µλ + µt) + (p0 − m0 − λµ) e λ .
e (1 − τ)dα αr
3.13 a) Sea ṗ = pe . Resolviendo se obtiene
r−α r−α
αr
pe (t) = p∗ + (p0e − p∗ ) e−( r−α )t
y
α αr
p(t) = p∗ − (p0e − p∗ ) e−( r−α )t .
r−α
(1 − τ) d
Además lim pe (t) = lim p(t) = ≡ p∗ .
t→∞ t→∞ r
b) Si τ aumenta a τ̄ > τ entonces en el momento del cambio, el cambio en
el precio ṗ pasa de un valor cero a un valor negativo (el precio tiende a
disminuir) correspondiente a la condición pe = p∗ . Después ṗ aumenta
en el tiempo y el sistema procede asintóticamente hacia un nuevo valor
de equilibrio, pe = p̄∗ < p∗ .
(1 − τ) d rt (1 − τ) d
c) p(t) = p0 − e + .
r r
(1 − τ̄) d
d) El nivel del precio diverge, a menos que p0 = , con τ̄ > τ.
r
Capı́tulo 4
Sistemas de ecuaciones
diferenciales lineales
1 1
4.1 a) X(t) = c1 e−2t + c2 e−4t .
1 −1
1 √ 1 √
b) X(t) = c1 √ e(2+ 3)t + c2 √ e(2− 3)t .
3 − 3
1 1
c) X(t) = c1 + c2 e−4t .
0 −4
2 3
d) X(t) = c1 + c2 et .
1 1
e−t cos t e−t sin t
4.2 a) X(t) = c1 + c2 .
−e−t sin t e−t cos t
4 4t
b) X(t) = c1 e3t + c2 e3t .
−2 1 − 2t
cos 2t t sin 2t
c) X(t) = c1 e + c2 et .
sin 2t − cos 2t
−1 −t
d) X(t) = c1 e t + c2 et .
−2 1 − 2t
1 2t 1 3t 1 1
4.3 a) X(t) = c1 e + c2 e − .
1 2 2 1
11
12
√
3
2 cos t
e− 2 t
3
b) X(t) = c1 √ √2 √
− cos 3
2 t+ 3 sin 3
2 t
√
3
2 sin t − 32 t 2
+c2 √ √2 √ e + .
− sin 2 t−
3
3 cos 3
2 t
5
1 1 1 5
c) X(t) = c1 e2t + c2 e3t − et .
1 2 2 4
−1 t 0
4.4 a) X(t) = e− 2 . Por lo tanto lim X(t) = .
10 t→∞ 0
cos βt
b) X(t) = eαt .
sin βt
5 5 cos t + 15 sin t
c) X(t) = 0 + 20 cos t + 10 sin t et .
0 30 cos t − 10 sin t
1 1
4.5 a) X(t) = (2 + w) et + (1 − w) e−2t .
1 −2
b) w = −2.
4.10 a) Como el conjunto {w1 , w2 , w3 , . . . , wn } es l.i. para todo t, por lo tanto las
columnas de Φ(t) son l.i. de modo que existe la inversa Φ−1 (t).
b) Sabemos que cada wi con i = 1, . . . n es solución de la ecuación Ẋ = AX,
por lo que se tiene que ẇi = Awi para i = 1 . . . n. Así
Φ̇(t) = ẇ1 ẇ2 . . . ẇn
= Aw1 Aw2 . . . Awn = A w1 w2 . . . wn = AΦ(t).
t
c) Sea Υ(t) = Φ(t) Φ−1 (s) f (s)ds. Por lo tanto
0
t
Φ−1 (s) f (s)ds = Φ−1 (t)Υ(t).
0
1 4 2
4.20 y(t) = e−2t + et cos t − et sin t. Además lim y(t) no está definido ya que la
5 5 5 t→∞
función oscila.
2v − 2u − 1 −x
4.21 a) y(x) = e
−5
2v + 3u − 1 x 2v − 2u + 4 x
+ e cos x + e sin x.
5 10
b) u = 1 y v = −1. Con esto se tiene y(x) = e−x . Por lo tanto lim y(x) = 0.
x→∞
1 t 1 −t 1
4.22 y(t) = e − e − sin t.
4 4 2
Capı́tulo 5
Análisis cualitativo
0
5.1 a) P∗ = es un punto silla.
0
0
b) P∗ = es un espiral atractora.
0
∗ 0
c) P = es un espiral repulsora.
0
0
d) P∗ = es un nodo repulsor.
0
∗ − 53
e) P = 1
es un punto silla.
3
0
5.2 a) P∗ = es un punto silla.
0
∗ 0
b) P = es degenerado inestable.
0
0
c) P∗ = es degenerado.
0
∗ −3b
d) P = es un punto fijo para cada b. Por lo tanto hay una infini-
b
dad de puntos fijos, cada uno de ellos degenerado.
9
e) P∗ = 4 es un espiral repulsora.
− 92
15
16
0
5.3 a) P∗ = 0 es degenerado.
0
0
∗
b) P = 0 es nodo repulsor.
0
0
c) P∗ = 0 es degenerado.
0
−4
d) P∗ = 1 es nodo repulsor.
1
0 0 2
5.4 a) Existen cuatro puntos fijos: P1∗ = , P2∗ = , P3∗ = ,y
0 6 0
∗ 4
P4 = . El punto P1∗ es un nodo repulsor, P2∗ un nodo atractor, P3∗
−2
x(t)
un punto silla y P4∗ una espiral atractora. Se tiene que lim =
t→∞ y(t)
0
es decir que la población de x se extinguirá y la de y tenderá a 6,
6
no es posible que las dos poblaciones coexistan en el largo plazo.
∗ 0 ∗ 0 ∗ 3
b) Existen cuatro puntos fijos: P1 = , P2 = , P3 = ,y
0 2 0
4
P4∗ = . El punto P1∗ es un nodo repulsor, P2∗ un punto silla, P3∗ un
−2
∗ x(t) 3
nodo atractor y P4 un punto silla. Se tiene que lim =
t→∞ y(t) 0
es decir que la población de x tenderá a 3 y la de y se extinguirá, no es
posible que las dos poblaciones coexistan en el largo plazo.
0 0 1
c) Existen cuatro puntos fijos: P1∗ = , P2∗ = , P3∗ = ,y
0 3 0
10
P4∗ = 3
14
. El punto P1∗ es un nodo repulsor, P2∗ un punto silla, P3∗ un
3
10
x(t)
punto silla y P4∗ un nodo atractor. Se tiene que lim = 3
14
t→∞ y(t) 3
17
k ∗ 1
5.11 a) El punto fijo es P∗ = = y es un punto silla porque para
c∗ 4
5
ese punto se obtiene λ1 = − 12 < 0 y λ2 = 45 > 0.
b) En P∗ la recta tangente que aproxima la variedad estable Ws (P∗ ) es c =
4
k. Similarmente, la recta tangente que aproxima la variedad inestable
5
1 13
Wu (P∗ ) es c = − k + .
2 10
4 ∗
c) c0 ≈ (1.1) = 0.88 > c .
5
5.12 v = −2u.
1
5.13 a) Los puntos fijos son cuatro: P1∗ = (0, 1), P2∗ = (0, −1), P3∗ = ( √ , 0), y
3
1
P4∗ = (− √ , 0). Los puntos P1∗ y P2∗ son los únicos puntos silla. Los
3
puntos P3∗ y P4∗ dan soluciones cíclicas. Para P1∗ se tiene
0 % &
Es (P1∗ ) = gen = (x, y) ∈ R2 | x = 0 ,
1
1 % &
E u
(P1∗ ) = gen = (x, y) ∈ R2 | y = 1 .
0
Para P2∗ se tiene
1 % &
Es (P2∗ ) = gen = (x, y) ∈ R2 | y = −1 ,
0
0 % &
Eu (P1∗ ) = gen = (x, y) ∈ R2 | x = 0 .
1
ẏ dy/dt dy dt dy
= = = .
ẋ dx/dt dt dx dx
Entonces
dy ẏ 1 − 3x2 − y2 1 − 3x2 y2
= = = − .
dx ẋ 2xy 2xy 2xy
Por lo tanto
dy 1 − 3x2 1
= y−1 − y,
dx 2xy 2x
que es una ecuación de Bernoulli con n = −1.
19
c)
% &
Ws (P1∗ ) = Wu (P2∗ ) = (x, y) ∈ R2 | x = 0
y
% &
Wu (P1∗ ) = Ws (P2∗ ) = (x, y) ∈ R2 | x2 + y2 = 1 .
∗ 0 ∗ 3
5.14 a) Hay dos puntos de equilibrio: P1 = y P2 = 4
. El punto P1∗
0 3
es un nodo repulsor y P2∗ es un punto silla.
0 − 1
b) Existen dos puntos de equilibrio: P1∗ = y P2∗ = 3 . El
0 − 13
punto P1∗ es un punto silla y P2∗ es un centro (soluciones cíclicas).
I(p∗ )
k∗ δ
5.15 El punto fijo = f (k∗) es un punto silla.
p∗ r+δ
r−r2 − 4 det J ∗
a) La tasa de convergencia está dada por λ = < 0, don-
2
de
det J ∗ = −δ(r + δ) + f (k∗ )I (p∗ ).
c
c) y(x) = ex + .
x−1
d) Para P∗ se tiene
% &
Ws (P∗ ) = (x, y) ∈ R2 | x = 1
y
% &
Wu (P∗ ) = (x, y) ∈ R2 | y = ex .
Capı́tulo 6
Conceptos básicos de dinámica
discreta
1 t
6.4 a) xt = − (x0 − 2) + 2. Además lim xt = 2, es decir que es asintótica-
2 t→∞
mente estable.
t
3
b) xt = (x0 + 4) − 4. Además lim xt = ∞, es decir que es asintótica-
2 t→∞
mente inestable.
5 5
c) xt = (−1)t x0 − + . Además lim xt no está definido es decir que
2 2 t→∞
diverge.
1 t
d) xt = − x0 . Ademáslim xt = 0, es decir que es asintóticamente esta-
3 t→∞
ble.
20
21
11 11
b) pt = −pt−1 + . El punto fijo es p∗ = el cual es inestable, se tiene
2 4
que lim pt no existe.
t→∞
c) pt = −3pt−1 + 16. El punto fijo es p∗ = 4 el cual es asintóticamente
inestable.
Capı́tulo 9
Optimización estática
22
23
Por lo tanto,
Entonces
De donde
(h ◦ g) (λ x̄1 + (1 − λ) x̄2 )
= h (g (λ x̄1 + (1 − λ) x̄2 ))
h (λg ( x̄1 ) + (1 − λ) g ( x̄2 ))
λh [g ( x̄1 )] + (1 − λ) h [g ( x̄2 )]
= λ (h ◦ g) ( x̄1 ) + (1 − λ) (h ◦ g) ( x̄2 ) .
(λx1 + (1 − λ) x2 ) (λy1 + (1 − λ) y2 ) k.
Así,
(λx1 + (1 − λ) x2 ) (λy1 + (1 − λ) y2 )
= λ2 (x1 y1 ) + (1 − λ)2 (x2 y2 ) + λ (1 − λ) (x2 y1 + x1 y2 )
2 x2 x1
kλ + k (1 − λ) + kλ (1 − λ)
2
+
x1 x2
2 2 x2 x1
= k λ + (1 − λ) + 2λ (1 − λ) + kλ (1 − λ) + −2
x1 x2
2
2 x1 + x22 − 2x1 x2
= k (λ + (1 − λ)) + kλ (1 − λ)
x1 x2
(x2 − x1 )2
= k + kλ (1 − λ)
x1 x2
(x2 − x1 )2
= k 1 + (1 − λ) k.
x1 x2
2
Como f xx = 2 > 0, f yy = 2 > 0 y |H| = f xx f yy − f xy = 0. Entonces H es
positiva semidefinida. Por lo tanto, f es convexa (no estricta).
c) Como f es doblemente diferenciable, analizamos simplemente el signo
de le matriz hessiana H :
f xx f xy 2 0
H= = .
f yx f yy 0 2
Por lo tanto
− αx12 0 ... 0
1
0 − αx22
... 0
H= 2 .
... ... ... ...
αn
0 0 . . . − x2
n
α1 α1 α2 α1 α2 α3
Como |H1 | = − 2 < 0, |H2 | = 2 2 > 0, |H3 | = − 2 2 2 < 0, . . . , (−1)k
x1 x1 x2 x1 x2 x3
|Hk | > 0 con 1 k n. Entonces H es definida negativa. Por lo tanto f es
estrictamente cóncava.
26
ā b̄
9.9 a) Sean ã = y b̃ = . Entonces
f (ā) f b̄
ā 1
f (ã) = f = f ā
f (ā) f ( ā)
1 f (ā)
= f ( ā) =
f ( ā) f (ā)
= 1.
Similarmente f b̃ = 1. Como CS f (1) = { x̄ ∈ X| f ( x̄) = 1} y como
f ( ã) = f b̃ = 1, por lo tanto ã, b̃ ∈ CS f (1) .
λ f b̄
b) Sea µ = . Como f : X → (0, ∞) , entonces f (ā) > 0
(1 − λ) f ( ā) + λ f b̄
y f b̄ > 0. Además, como 0 < λ < 1 se tiene que λ > 0 y (1 − λ) > 0,
por lo tanto µ > 0. Por otra parte, reescribamos µ como
λ f b̄ + (1 − λ) f ( ā) − (1 − λ) f (ā)
µ=
(1 − λ) f (ā) + λ f b̄
(1 − λ) f (ā)
= 1− < 1,
(1 − λ) f (ā) + λ f b̄
(1 − λ) f ( ā)
ya que > 0. Por lo tanto µ < 1. Se concluye que
(1 − λ) f (ā) + λ f b̄
0 < µ < 1.
c) Como ã, b̃ ∈ CS f (1) , 0 < µ < 1 y CS f (1) es convexo, entonces
(1 − µ) ã + µb̃ ∈ CS f (1) .
Por lo tanto, f (1 − µ) ã + µb̃ 1.
d) De la definición de µ se tiene
λ f b̄ (1 − λ) f (ā)
1−µ = 1− = .
(1 − λ) f (ā) + λ f b̄ (1 − λ) f ( ā) + λ f b̄
27
Entonces
1 f (1 − µ) ã + µb̃
(1 − λ) f (ā) ā λ f b̄ b̄
= f +
(1 − λ) f (ā) + λ f b̄ f ( ā) (1 − λ) f ( ā) + λ f b̄ f b̄
1 ' (
= f (1 − λ) ā + λb̄
(1 − λ) f ( ā) + λ f b̄
1
= f (1 − λ) ā + λb̄ .
(1 − λ) f ( ā) + λ f b̄
1
Es decir, 1 f (1 − λ) ā + λb̄ . Por lo tanto
(1 − λ) f (ā) + λ f b̄
(1 − λ) f ( ā) + λ f b̄ f (1 − λ) ā + λb̄ .
Por lo tanto, (−1)k |Hk | > 0. Entonces, H es definida negativa. Por lo tanto, f
n
es estrictamente cóncava. Si 0 < ∑ αk 1, entonces f es cóncava.
k=1
28
1
9.11 a) w (λx1 , . . . , λxn ) = δ1 (λx1 )ρ + · · · + δn (λxn )ρ ρ
ρ ρ 1
= λρ δ1 x1 + · · · + δn xn ρ
ρ ρ 1
= λ δ1 x1 + · · · + δn xn ρ
= λw (x1 , . . . , xn ) .
n
ρ ρ ρ
d) Sea g (x1 , . . . , xn ) = ∑ δk xk = δ1 x1 + · · · + δn xn . Entonces,
k=1
ρ−2
δ1 ρ (ρ − 1) x1 0 ... 0
ρ−2
H= 0 δ2 ρ (ρ − 1) x2 ... 0 .
ρ−2
0 0 ... δn ρ (ρ − 1) xn
con ā, b̄ ∈ X y 0 < λ < 1. De modo que (problema 9.9), nuevamente 0 <
µ < 1. Como ã, b̃ ∈ CI f (1), 0 < µ < 1 y CI f (1) es convexo (ya que f es
cuasiconvexa), entonces f (1 − µ) ã + µb̃ 1. Finalmente, sustituyendo ã, b̃
y µ en esta desigualdad y, viendo que f es homogénea de grado 1, se obtiene
que
(1 − λ) f ( ā) + λ f b̄ f (1 − λ) ā + λb̄ .
Por lo tanto, f es convexa. Ahora, apliquemos este resultado a la función
CES,
ρ ρ 1
w (λx1 , . . . , λxn ) = δ1 x1 + · · · + δn xn ρ .
Es claro que w es homogénea de grado 1 (problema 9.11.a) y positiva. Fal-
ta ver que w sea cuasicóncava cuando ρ > 1, para aplicar el teorema recién
30
9.13 Sea Ω = (a, b) un conjunto convexo y abierto de R y sea y = f (z) una función
cóncava en Ω.
a) Sean a < z1 < z2 < z3 < b con z2 = λz1 + (1 − λ) z3 , 0 < λ < 1. Como f
es cóncava en Ω, entonces
Por lo tanto,
y2 λy1 + (1 − λ) y3 ....... (1) .
Por lo tanto,
(y2 − y1 ) (y − y2 )
Como z3 − z2 > 0 y z2 − z1 > 0, entonces 3 . Por lo
(z2 − z1 ) (z3 − z2 )
tanto
f (z2 ) − f (z1 ) f (z3 ) − f (z2 )
....... (4) .
z2 − z1 z3 − z2
b) Como Ω es convexo y por la densidad de los reales, siempre se pueden
escoger números r, s, t, u tales que a < r < s < x < y < t < u < b, para
cada x ∈ (a, b) .
c) Sean a < r < s < x < t < u < b. Podemos aplicar la ecuación (4) en
cada trío de puntos en r < s < x < y < t < u, obteniendo
f (s) − f (r) f (x) − f (s)
s−r x−s
f (y) − f (x) f (t) − f (y) f (u) − f (t)
....... (5) ,
y−x t−y u−t
con s, r, u y t fijos. Así, se definen las constantes c1 y c2 como:
f (s) − f (r) f (u) − f (t)
c1 = , c2 = ....... (6) .
s−r u−t
De (5) y (6) se obtiene
f (y) − f (x)
c1 c2 ....... (7) .
y−x
d) Por último, como y − x > 0 entonces lim+ c1 (y − x) f (y) − f (x)
y→x
c2 (y − x) . Por lo tanto,
Como
lim [c1 (y − x)] = lim+ [c2 (y − x)] = 0,
y→x + y→x
entonces
lim [ f (y) − f (x)] = 0.
y→x +
Por lo tanto,
lim f (y) = f (x) .
y→x +
Por lo tanto,
lim f (y) = f (x) .
y→x
1
f (y) ∼ T T
= f (x) + (y − x) ∇ f + (y − x) [H f (x)] (y − x) ,
2
donde H f (x) es el hessiano de f evaluado en x.
(y − x) T [H f (x)] (y − x) 0.
(y − x) T [H f (x)] (y − x) < 0.
(y − x) T [H f (x)] (y − x) > 0.
dy gx
= − , gy = 0.
dx gy
Entonces, ∗
dF (x) dy
= f x∗ + f y∗ = 0.
dx dx
Lo que implica ∗
gx
f x∗ + f y∗ − = 0.
gy
De donde f x∗ gy∗ − f y∗ g∗x = 0. Por lo tanto,
) )
) f∗ f y∗ )
) x )
) ∗ ) = 0,
) gx gy∗ )
con g (x∗ , y∗ ) = 0.
C (w, q) = w · x∗ (w, q) .
∂C
= x∗j (w, q) ,
∂w j
para j = 1, . . . , n.
∂C ∂C ∂C n
∂C
w1
∂w1
+ w2
∂w2
+ · · · + wn
∂wn
= ∑ w j ∂w j
j=1
n
= ∑ w j x∗j (w, q)
j=1
= C (w, q) .
Se tiene que
1
Lx = − 3λ1 − λ2 = 0,
3x
1
Ly = − λ1 − λ2 = 0,
3y
Lλ1 = A − (3x + y) 0, λ1 0, λ1 (3x + y − A) = 0,
Lλ2 = 40 − (x + y) 0, λ2 0, λ2 (x + y − 40) = 0.
9.19 Sea
L (x, q, λ; w, p) = pq − w T x − λ [ f (x) − q] .
Entonces por las condiciones de primer orden, se tiene que la solución es del
tipo
x∗ = x∗ (w, p) ,
q∗ = q∗ (w, p) = f (x∗ (w, p)) ,
λ∗ = λ∗ (w, p) .
∂Π ∂L
= = −x∗j (w, p) , j = 1, ...n.
∂w j ∂w j
∂Π ∂L
= = q∗ (w, p) .
∂p ∂p
36
La función de gasto es
E p, Ū = L xh , λh ; p, Ū = pxh p, Ū − λ ∗ 0.
∂E ∂L
= = xhj p, Ū ,
∂p j ∂p j
para j = 1, . . . , n.
9.21 a) El problema
min p T x
,
s.a U (x) = V (p, m)
implica que
Por lo tanto,
xh = xh (p, V (p, m)) .
b) El problema
max U (x)
,
s.a p T x = E p, Ū
implica que
L x; p, E p, Ū = U (x) − λ p T x − E p, Ū .
37
Por lo tanto,
x∗ = x∗ p, E p, Ū .
αm ∗ βm
9.22 a) x∗ (p, m) = , y (p, m) = .
p1 p2
α α β β
b) V (p, m) = ln m .
p1 p2
p1 α p2 β Ū
c) E p, Ū = e .
α β
αp2 β Ū h βp1 α Ū
d) xh p, Ū = e , y p, Ū = e .
βp1 αp2
1 1
∗ p1r−1 ∗ p2r−1
9.23 a) x (p, m) = m r r , y y (p, m) = m r r .
p1r−1 + p2r−1 p1r−1 + p2r−1
r r
r
1−r
1 r r
− 1
r
yh p, Ū = Ū p2r−1 p1r−1 + p2r−1 .
38
9.24 Se tiene x∗ p, E p, Ū = xh p, Ū . Entonces
∂xi∗ ∂xi∗ ∂E p, Ū ∂xih p, Ū
p, E p, Ū + p, E p, Ū = .
∂p j ∂m ∂p j ∂p j
∂E
Por el Lema de Shepard = xhj , y como Ū = V (p, m) , m = E p, Ū y
∂p j
h ∗
x j (p, V (p, m)) = x j (p, m) . Por lo tanto
39
40
1p
b
11.2 Por demostrar que f p = | f | p dt es una norma.
a
d) Si p = 1, entonces
b b
f + g1 = | f + g| dt (| f | + |g|) dt
a a
b b
= | f | dt + |g| dt = f 1 + g1 .
a a
Por lo tanto,
* * *
b b b
| f + g|2 dt | f |2 dt + |g|2 dt.
a a a
( f + g) + h = f + (g + h) .
y J2 [ f ] = f (0) .
c) Dos posibles ejemplos de funcionales no lineales sobre V son:
1
J1 [ f ] = f 2 (t) dt
0
2
1
y J2 [ f ] = f (t) dt .
0
1
11.4 a) x (t) = − t + 20. Por lo tanto, J [x] = −5.
2
b) x (t) = 9t + 10. Por lo tanto, J [x] = −10710.
1912
c) x (t) = t3 + 4t + 1. Por lo tanto, J [x] = .
5
t2 154
d) x (t) = + 4t + 1. Por lo tanto, J [x] = .
4 3
42
16
e) x (t) = t. Por lo tanto, J [x] = .
3
f) Se tiene que
ẍ = 2x − y
.
ÿ = x
Por lo tanto,
11 2
g) x (t) = t y y (t) = t + 2.
10 5
h) Se tiene que
ẍ = y
.
ÿ = x
1 t −t
Por lo tanto, x (t) = y (t) = π π e − e .
e − e− 2
2
11.8 Teniendo el modelo de inversión de la sección 11.7.1 como base, entonces sus-
tituyendo k(t) y su demanda k̇(t) en la condición de transversalidad se tiene
2
0 = lim e−ρT α k1 r1 er1 T + k2 r2 er2 T + A k1 er1 T + k2 er2 T + kρ
t→∞
2
−ρT
− lim e r1 T r2 T
B k1 e + k2 e + kρ
t→∞
+ ' ( ' (,
= lim e(2r1 −ρ)T k21 αr12 − B + e(2r2 −ρ)T k22 αr22 − B
t→∞
+ ,
+ lim e(r1 +r2 −ρ)T 2k1 k2 [αr1 r2 − B] + e(r1 −ρ)T k1 [A − 2Bkρ]
t→∞
+ ' (,
+ lim e(r2 −ρ)T k2 [A − 2Bkρ] + e−ρT Akρ − Bkρ2 .
t→∞
43
Para que este límite converja es necesario que no aparezcan aquí los términos
e(2r1 −ρ)T y e(r1 −ρ)T que son divergentes, y esto se logra pidiendo que k1 = 0.
En este caso,
+ ' 2 ( ' (,
(2r2 −ρ)T 2 (r2 −ρ)T −ρT 2
lim e k2 αr2 − B + e k2 [A − 2Bkρ] + e Akρ − Bkρ = 0,
t→∞
se verifica automáticamente.
t2 √
11.9 x (t) = y T = 2N.
2
1 1
11.10 a) x (t) = c2 e−t + e −ρt
. Pero como c 2 = x 0 − , entonces
1 − ρ 1 − ρ2
2
1 1
x (t) = x0 − e−t + e−ρt .
1−ρ 2 1 − ρ2
b) Se tiene que f x = e−ρt − x y f ẋ = − ẋ, lo que implica que f xx = −1 < 0,
f x ẋ = 0 y f ẋ = −1. Por lo tanto,
−1 0
H= .
0 −1
De donde |H| = 1 > 1. Por lo tanto, f es cóncava en (x, ẋ) , es decir que
se trata de un máximo.
c) Al imponer la condición c1 = 0, y dado que ρ > 0, entonces
1 −t 1 −ρt
lim x (t) = lim x0 − e + e = 0.
t→∞ t→∞ 1 − ρ2 1 − ρ2
Por lo tanto, sí es cierto que
lim x (t) = 0.
t→∞
k̇ = (A − δ) k − c,
ċ = (A − δ − ρ) c.
Resolviendo se tiene,
k 1 (A−δ)t 1
= c1 e + c2 e(A−δ−ρ)t .
c 0 ρ