Mathematical Economics Lecture Notes: Alexander W. Richter
Mathematical Economics Lecture Notes: Alexander W. Richter
Lecture Notes1
Alexander W. Richter2
Research Department
Federal Reserve Bank of Dallas
Department of Economics
Southern Methodist University
August 2017
1
I am especially grateful to Juergen Jung, Mike Treuden, and Nathaniel Throckmorton for their tremen-
dous contributions to these notes. I also thank Michael Kaganovich and Eric Leeper for their guidance and
for giving me the opportunity to develop and teach this course during graduate school. Several classes of In-
diana University and Auburn University graduate students provided suggestions and corrections. Comments
are welcome and all remaining errors are mine. The views expressed in these notes are my own and do not
necessarily reflect the views of the Federal Reserve Bank of Dallas or the Federal Reserve System.
2
Correspondence: Research Department, Federal Reserve Bank of Dallas, 2200 N Pearl St, Dallas, TX
75201, USA. Phone: +1(922)922-5360. E-mail: [email protected] or [email protected].
Contents
Preface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . v
Chapter 1 Mathematical Preliminaries 1
1.1 Single-Variable Calculus . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1 Limits of Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.2 Definition of a derivative and Tangent Lines . . . . . . . . . . . . . . 2
1.1.3 Properties of the Differential . . . . . . . . . . . . . . . . . . . . . . 3
1.1.4 Single-Variable Maximization . . . . . . . . . . . . . . . . . . . . . 5
1.1.5 Intermediate and Mean Value Theorems . . . . . . . . . . . . . . . . 7
1.1.6 Taylor Approximations . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.1.7 Laws of Logarithms . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.8 Infinite Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.2 Multivariate Calculus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2.1 Level Surfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2.2 Projections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2.3 Gradient Vector and its Relationship to the Level Surface . . . . . . . 16
1.2.4 Gradients and Tangent Planes . . . . . . . . . . . . . . . . . . . . . . 17
1.2.5 Chain Rule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.2.6 Second Order Derivatives and Hessians . . . . . . . . . . . . . . . . . 23
1.3 Basic Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
1.3.1 Induction and Examples . . . . . . . . . . . . . . . . . . . . . . . . . 24
1.3.2 Neighborhoods and Open and Closed Sets . . . . . . . . . . . . . . . 26
1.3.3 Convergence and Boundedness . . . . . . . . . . . . . . . . . . . . . 29
1.3.4 Compactness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Chapter 2 Basic Matrix Properties and Operations 33
2.1 Determinants . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.1.1 Minors, Cofactors, and Evaluating Determinants . . . . . . . . . . . . 33
2.1.2 Properties of Determinants . . . . . . . . . . . . . . . . . . . . . . . 33
2.1.3 Singular Matrices and Rank . . . . . . . . . . . . . . . . . . . . . . . 35
2.2 Inverses of Matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.2.1 Computation of Inverses . . . . . . . . . . . . . . . . . . . . . . . . 36
2.2.2 Properties of Inverses . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.3 Quadratic Forms and Definiteness . . . . . . . . . . . . . . . . . . . . . . . 37
2.3.1 Quadratic Forms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.3.2 Definiteness of Quadratic Forms . . . . . . . . . . . . . . . . . . . . 38
2.4 Eigenvalues and Eigenvectors . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.4.1 Properties of Eigenvalues and Eigenvectors . . . . . . . . . . . . . . 45
2.4.2 Definiteness and Eigenvalues . . . . . . . . . . . . . . . . . . . . . . 45
ii
A. W. Richter CONTENTS
iii
List of Figures
iv
A. W. Richter PREFACE
Preface
These notes are intended for a one-semester course in mathematical economics. The goal of this
course is to help prepare students for the mathematical rigor of graduate economics by providing a
balance between theory and application. There are dozens of examples throughout the notes, which
demonstrate many of the theoretical concepts and theorems. Below is a short list of the notation
used thoughout these notes.
v
Chapter 1
Mathematical Preliminaries
(b) limx→5 x2 − 3x + 1 = 11
Solution: Let ε > 0, suppose |x − 5| < δ, and choose δ = min{1, ε/8}. We can write
|f (x) − 11| = |x2 − 3x − 10| = |(x − 5)(x + 2)|.
To make this small, we need a bound on the size of x + 2 when x is “close” to 5. For example,
if we arbitrarily require that |x − 5| < 1, then
|x + 2| = |x − 5 + 7| ≤ |x − 5| + 7 < 8.
To make f (x) within ε units of 11, we shall want to have |x + 2| < 8 and |x − 5| < ε/8.
Thus, under the above definition of δ
|f (x) − 11| = |(x − 5)(x + 2)| < 8δ ≤ ε.
1
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
f
L+ε
L-ε
a-δ a a+δ
(c) limx→−2 x2 + 2x + 7 = 7
Solution: Let ε > 0, suppose |x + 2| < δ, and choose δ = min{1, ε/3}. We can write
|f (x) − 7| = |x2 + 2x| = |x(x + 2)|.
To make this small, we need a bound on the size of x when x is “close” to −2. For example,
if we arbitrarily require that |x + 2| < 1, then
|x| − |2| ≤ |x + 2| < 1,
since |a| − |b| ≤ ||a| − |b|| ≤ |a ± b| by the triangle inequality. Thus, |x| < 3, which implies
|f (x) − 7| = |x(x + 2)| < 3δ ≤ ε.
(d) limx→2 x3 = 8
ε
Solution: Let ε > 0, suppose |x − 2| < δ, and choose δ = min{1, 19 }. Then
|f (x) − 8| = |x3 − 8| = |x − 2||x2 + 2x + 4|
= |x − 2| · |(x − 2)2 + 6(x − 2) + 12|
< δ(δ2 + 6δ + 12)
ε
≤ (1 + 6 + 12)
19
= ε.
2
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
Definition 1.1.3 (Tangent Line). The tangent line to the curve y = f (x) at the point P (a, f (a)) is
the line through P with slope m, provided that the limit exists.
Remark 1.1.1. In general, when a function f of one variable is differentiable at a point a, the
equation of a tangent line to the graph of f at a is:
y = f (a) + f ′ (a)(x − a)
Example 1.1.2. Find an equation of the tangent line to the hyperbola y = 3/x at the point (3, 1)
using the above definition of a derivative.
Solution: The slope is given by
3
f (3 + h) − f (3) 3+h−1
m = lim = lim
h→0 h h→0 h
1 1
= lim − =− .
h→0 3 + h 3
Therefore, an equation of the tangent line at the point (3, 1) is
y = 1 − (x − 3)/3 → x + 3y − 6 = 0.
Example 1.1.3. Find an equation of the tangent line to the curve y = (x − 1)/(x − 2) at the point
(3, 2) using the above definition of a derivative.
Solution: The slope is given by
h+2
f (3 + h) − f (3) h+1−2
m = lim = lim
h→0 h h→0 h
1
= lim − = −1.
h→0 h + 1
y − 2 = −(x − 3) → x + y = 5.
3
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
Example 1.1.4. Using the properties of the differential, differentiate the following, where a, p, q,
and b are constants
1
(a) y = f (x) = (x2 +x+1)5
−5(2x+1)
Solution: f ′ (x) = (x 2 +x+1)6
√
q p
(b) y = f (x) = 1 + 1 + x
p √
Solution: f ′ (x) = 1/(8y x(1 + x))
Then
f (x) f ′ (x)
lim = lim ′
x→a g(x) x→a g (x)
4
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
Solution:
ln x 1/x 3
lim √ = lim −2/3 = lim √ = 0.
x→∞ 3 x x→∞ x /3 x→∞ 3
x
Theorem 1.1.4 (Inverse Function Theorem). Suppose that f is differentiable on an interval I and
f ′ (x) 6= 0 (no local maxima or minima) for all x ∈ I. Then f is injective, f −1 is differentiable on
f (I), and
∂x 1 1
= (f −1 )′ (y) = ′ = ,
∂y f (x) ∂y/∂x
where y = f (x).
Example 1.1.9. Let n ∈ N and y = f (x) = x1/n for x > 0. Then f is the inverse of the
function g(y) = y n . Use Theorem 1.1.4 to verify the familiar derivative formula for f : f ′ (x) =
(1/n)x1/n−1 .
Solution:
1 1 1 1 1
f ′ (x) = = ′ = = = x1/n−1 .
(f −1 )′ (y) g (y) ny n−1 n(x1/n )n−1 n
θ
y = f (θ) = − .
(1 − θ) log(1 − θ)
(c) If f ′ does not change sign at c, then f has no local maximum or minimum at c.
If the sign of f ′ (x) changes from positive to negative (negative to positive) at c, f is increasing
(decreasing) to the left of c and decreasing (increasing) to the right of c. If follows that f has a local
maximum (minimum) at c.
5
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
If f ′′ (c) > 0(< 0) near c, f is concave upward (downward) near c. Thus, the graph of f lies
above (below) its horizonal tangent at c and so f has a local minimum (maximum) at c.
√
Example 1.1.11. The height of a plant after t months is given by h(t) = t − t/2, t ∈ [0, 3]. At
what time is the plant at its highest?
Solution: The first order condition is given by
1 1 set
h′ (t) = √ − = 0 ⇒ t∗ = 1,
2 t 2
where h(1) = 1 −√1/2 = 1/2. Since h′′ (t) = −1/(4t3/2 ) < 0, t∗ = 1 is a local maximum. Also,
h(0) = 0, h(3) = 3 − 3/2 ≈ 0.27. Therefore, t∗ = 1 is an absolute maximum.
Example 1.1.12. A sports club plans to charter a plane. The charge for 60 passengers is $800 each.
For each additional person above 60, all travelers get a discount of $10. The plane can take at most
80 passengers.
(b) Find the number of passengers that maximizes the total airfare paid by the club members.
Solution: T C ′ (x) = 200− 20x and T C ′′ (x) = −20. Thus, x∗ = 10 and airfare expenditures
are maximized with 70 passengers (T C(x∗ ) = $49, 000).
Example 1.1.13. Let C(Q) be the total cost function for a firm producing Q units of some com-
modity. A(Q) = C(Q)/Q is then called the average cost function. If C(Q) is differentiable, prove
that A(Q) has a stationary point (critical point) at Q0 > 0 if and only if the marginal cost and the
average cost functions are equal at Q0 . (C ′ (Q0 ) = A(Q0 ))
Solution: By definition, QA(Q) = C(Q). Differentiating with respect to Q yields
Assume A(Q) has a stationary point at Q0 > 0. Then it is easy to see that A(Q0 ) = C ′ (Q0 ) as
desired. Now assume A(Q0 ) = C ′ (Q0 ). Then Q0 A′ (Q0 ) = 0, which implies that A′ (Q0 ) = 0 as
desired since Q0 > 0.
Example 1.1.14. With reference to the previous example, let C(Q) = aQ3 + bQ2 + cQ + d, where
a > 0, b ≥ 0, c > 0, and d > 0. Prove that A(Q) = C(Q)/Q has a minimum in the interval (0, ∞).
Then let b = 0 and find the minimum point in this case.
Solution: The average cost function is given by
C(Q) d
A(Q) = = aQ2 + bQ + c + .
Q Q
Differentiating with respect to Q yields
d set
A′ (Q) = 2aQ + b − = 0.
Q2
6
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
f
f(b)
N
f(a)
a c b
Given the restrictions on the constants, as Q → 0, A′ (Q) < 0 and as Q → ∞, A′ (Q) > 0. Thus,
there exists a Q∗ ∈ (0, ∞) that satisfies the first order condition. To determine whether this critical
point is a minimum or maximum, differentiate A′ (Q) with respect to Q to obtain
given the restrictions on the parameters. Thus, A(Q) has a minimum in the interval (0, ∞) by the
second derivative test. When b = 0
1/3
d ∗ ∗ d
∗ 2
= 2aQ → Q = .
(Q ) 2a
7
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
f
f(b)
secant
f(c)
f(a)
tangent at c
a c b a c b
(a) Rolle’s Theorem (b) Mean Value Theorem
Theorem 1.1.6 (Mean Value Theorem). Let f be a continuous function on [a, b] that is differentiable
on (a, b). Then there exists at least one point c ∈ (a, b) such that
f (b) − f (a)
f ′ (c) = .
b−a
Example 1.1.16. As a illustration of one use of the Mean Value Theorem (MVT), we will derive
Bernoulli’s inequality: for x > 0
(1 + x)n ≥ 1 + nx ∀n ∈ N.
Solution: Let f (t) = (1 + t)n on the interval [0, x], which is clearly continuous and differen-
tiable. Then, by the MVT, there exists a c ∈ (0, x) such that
Thus, we have
(1 + x)n − 1 = nx(1 + c)n−1 ≥ nx,
since f ′ (c) = n(1 + c)n−1 , 1 + c > 1, and n − 1 ≥ 0.
Example 1.1.17. Use the Mean Value Theorem (MVT) to establish the following inequalities, as-
suming any relevant derivative formulas.
(a) ex > 1 + x for x > 0
Solution: Define f (x) = ex and recall from the MVT that f (b) − f (a) = f ′ (c)(b − a) for
some c ∈ (0, x). Then
ex − 1 = ec (x − 0) = ec x > x
since ec > 1 for c > 0.
8
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
1
√ 1
(b) 8 < 51 − 7 < 7
√
Solution: Define f (x) = x and consider the interval [49, 51]. Then by the MVT, ∃c ∈
(49, 51) such that √ √
51 − 49 1
= f ′ (c) = √ .
51 − 49 2 c
Consequently, we have
1 √
√ = 51 − 7,
c
which implies
1 1 1 1 √ 1 1
= √ < √ < √ = 51 − 7 < √ = .
8 64 51 c 49 7
√ x−24
(c) 1+x<5+ 10 for x > 24
√
Solution: Define f (x) = 1 + x and consider the interval [24, x]. Then by the MVT, ∃c ∈
(24, x) such that
√ x − 24 x − 24
f (x) − f (24) = 1 + x − 5 = f ′ (c)(x − 24) = √ < ,
2 1+c 10
√ √
since c > 24, 1 + c > 25 = 5.
9
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
Remark 1.1.2. For the special case x0 = 0, the Taylor Series is known as a Maclaurin series.
Example 1.1.19. Find an approximation to the following function about the given point
f (x) = (1 + x)5 , a = 0.
Solution: Applying the above formula
5(1 + a)4 20(1 + a)3 60(1 + a)2
f (x) ≈ (1 + a)5 + (x − a) + (x − a)2 + (x − a)3 +
1! 2! 3!
120(1 + a) 4 120 5
(x − a) + (x − a)
4! 5! a=0
= 1 + 5x + 10x2 + 10x3 + 5x4 + x5 .
Example 1.1.20. Find quadratic approximations to the following functions about the given points:
(a) F (K) = AK α , K0 = 1
Solution: F (K) ≈ A[1 + α(K − 1) + α(α − 1)(K − 1)2 /2]
1/2
(b) f (ε) = 1 + 32 ε + 12 ε2 , ε0 = 0
Solution: f (ε) ≈ 1 + 34 ε − 1 2
32 ε
10
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
Natural Logarithms
The mathematical constant, e, is the unique real number such that the value of the derivative (the
slope of the tangent line) of the exponential function f (x) = ax at the point x = 0 is exactly 1. The
logarithm with base e is called the natural logarithm and has a special notation
loge x ≡ ln x.
ln x = y ⇐⇒ ey = x
ln(ex ) = x for every x ∈ R
eln x = x for every x > 0
ln e = 1.
Finally, when y = loga x, we have ay = x. Thus, applying the natural logarithm to both sides of
this equation, we get y ln a = ln x. Thus, the change of base formula is given by
ln x
y = loga x = .
ln a
Example 1.1.21. Express ln a + 12 ln b as a single logarithm.
√
Solution: ln a + 12 ln b = ln a + ln b1/2 = ln(a b)
Example 1.1.22. Find the inverse function of the following: m = f (t) = 24 · 2−t/25 .
m t
Solution: 24 = 2−t/25 ⇒ ln m − ln 24 = − 25 ln 2 ⇒ t = f −1 (m) = ln252 (ln 24 − ln m)
Example 1.1.23. If f (x) = 2x + ln x, find f −1 (2)
Solution: Define y = f (x). Then f −1 (y) = x. Thus, at y = 2, we have 2x+ln x = 2 ⇒ x = 1.
It immediately follows that f −1 (2) = 1
Example 1.1.24. Calculate limx→∞ (1 + 1/x)x .
Solution: Define y = (1 + 1/x)x . Then ln y = x ln(1 + 1/x). We must first evaluate the limit
of the right-hand-side as x → ∞. Using L’Hospital’s Rule, we obtain
ln(1 + 1/x) 1
lim x ln(1 + 1/x) = lim = lim = 1.
x→∞ x→∞ 1/x x→∞ 1 + 1/x
where x̄ is the stationary value of xt . Then show that percent changes are a good approximation for
log deviations.
Solution:
xt ≈ f (0) + f ′ (0)(x̂t − 0)
= x̄e0 + x̄e0 (x̂t − 0)
= x̄[1 + x̂t ]
11
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
1
xt −x̄
x̄
0.8
Error ln(xt /x̄)
0.6
0.4
Better Approx.
0.2
Around 1
0
−0.2
−0.4
−0.6
−0.8
−1
0 0.5 1 1.5 2 2.5
a1 + a2 + a3 + · · · + an + · · · ,
12
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
The logical question is whether it makes sense to talk about the sum of infinitely many terms? It
would be impossible to find a finite sum for the series
1 + 2 + 3 + 4 + ··· + n + ···
because if we start adding the terms we get the cumulative sums 1, 3, 6, 10, . . . and after the nth
term, we get n(n + 1)/2, which becomes very large as n increases. However, if we start to add the
terms of the series
1 1 1 1 1
+ + + + ··· + n + ···
2 4 8 16 2
we get 12 , 34 , 78 , 15 31 n
16 , 32 , . . . , 1 − 1/2 , . . ., which become closer and closer to 1. Using this idea, we
can define a new sequence {sn } of partial sums given by
n
X
sn = ak = a1 + a2 + · · · + an .
k=1
If {sn } converges to a real number s, we say that the series is convergent and we write
∞
X
an = s.
n=1
1 1 1 1
sn = + + + ··· +
1·2 2·3 3·4 n(n + 1)
1 1 1 1 1 1 1 1
= − + − + − + ··· + −
1 2 2 3 3 4 n n+1
1
=1− .
n+1
This is an example of a telescoping series, so called because of the way in which the terms in the
partial sums cancel. Since the sequence of partial sums converges to 1 for n large, we have
∞
X 1
= 1.
n=1
n(n + 1)
is divergent.
Solution: The partial sums are given by
s1 = 1
13
A. W. Richter 1.1. SINGLE-VARIABLE CALCULUS
1
s2 = 1 +
2
1 1 1 1 1 1 2
s4 = 1 + + + >1+ + + =1+
2 3 4 2 4 4 2
1 1 1 1 1 1 1
s8 = 1 + + + + + + +
2 3 4 5 6 7 8
1 1 1 1 1 1 1
>1+ + + + + + +
2 4 4 8 8 8 8
1 1 1 3
=1+ + + =1+
2 2 2 2
Similarly, s16 > 1 + 42 , s32 > 1 + 52 , and in general
n
s2n > 1 + .
2
This shows s2n → ∞ as n → ∞ and so {sn } is divergent. Therefore the harmonic series diverges.
Theorem 1.1.7. If ∞
P
n=1 an is a convergent series, then limn→∞ an = 0.
P∞
Proof. If n=1 an converges, then the sequence of partial sums {sn } must have a finite limit. But
an = sn − sn−1 , so limn→∞ an = limn→∞ sn − limn→∞ sn−1 = 0.
14
A. W. Richter 1.2. MULTIVARIATE CALCULUS
Example 1.1.29. In the previous example, we found that for |r| < 1
∞ ∞
X
n−1
X a
ar = ar n = .
1−r
n=1 n=0
This result is particularly useful for deriving a closed solution for the expected value of a discrete
geometric random variable.
Example 1.1.30. Find the sum of each series
P∞ 1 n P∞ 1 n
1. n=1 3 2. n=3 2
1 n 1/3 1 n 1/8
Solution: ∞
P∞
= 12 1
P
n=1 3 = 1−(1/3) Solution: n=3 2 = 1−(1/2) = 4
1.2.2 Projections
Definition
P1.2.2 (Scalar Product). For x, y ∈ Rn , the scalar product (or dot product) of x and y is
n
x · y := i=1 xi yi .
Definition 1.2.3 (Orthogonal Vectors). For x, y ∈ Rn , x and y are orthogonal if x · y = 0.
√
Definition 1.2.4 (Vector Length). For x ∈ Rn , the length (or norm) of x is kxk := x · x.
Theorem 1.2.1. If θ is the angle between vectors x and y, then
x · y = kxkkyk cos θ.
15
A. W. Richter 1.2. MULTIVARIATE CALCULUS
x
x - ty
ty y
In order to see these formulas more clearly, note that cos θ = ktyk/kxk. Thus,
which is the formula for the scalar projection of x onto y. Moreover, using the fact that ty =
ktyky/kyk, simple algebra yields the vector projection of x onto y.
For t = (x · y)/kyk2 , (x − ty) · y = 0 (See figure 1.5). Thus, we can decompose x into two
parts, one a multiple of y, ty, and the other orthogonal to y, x − ty.
Given a vector y, which vector x with norm c > 0 maximizes x · y? The set of vectors with
kxk = c consists of all those vectors with heads on the “sphere” (could be a higher dimension) with
radius c. To simplify the problem assume kyk = 1. For any x, if the projection of x on y is ty, then
x·y
projy x ≡ ty = y.
kyk2
16
A. W. Richter 1.2. MULTIVARIATE CALCULUS
Level Curve
Let t0 be the parameter value corresponding to P ; that is r(t0 ) = hx1,0 , x2,0 , . . . , xn,0 i. Since C
lies on S, any point (x1 (t), x2 (t), . . . , xn (t)) must satisfy
Since ∇f = hfx1 , fx2 , . . . , fxn i and r′ (t) = hx′1 (t), x′2 (t), . . . , x′n (t)i, at t = t0 the above condition
can be written
Thus, the gradient vector at P , ∇f (x1,0 , x2,0 , . . . , xn,0 ), is perpendicular to the tangent vector r′ (t)
to any curve C on S that passes through P . To illustrate, consider a function f of two variables and
a point P (x0 , y0 ) in its domain. The gradient vector ∇f (x0 , y0 ) gives the direction of the fastest
increase of f and is perpendicular to the level curve f (x, y) = k that passes through P . This makes
intuitive sense since the values of f remain constant as we move along the curve.
(i) f is continuous at a.
(ii) For all unit vectors, u, the directional derivative in the direction of u is fu (a) = ∇f (a) · u.
∂f (a)
(iii) For all i the ith component of ∇f (a) is ∂xi .
17
A. W. Richter 1.2. MULTIVARIATE CALCULUS
(iv) The equation of the tangent plane to the graph of f at (a, f (a)) is given by
(v) The equation of the tangent plane to the level curve corresponding to f (x) ≡ f (a) at the
point a is given by
0 = ∇f (a) · (x − a).
(vi) The marginal rate of substitution of xi for xj along the level curve corresponding to f (x) ≡
f (a) at the point a is the number of units of xj , which must be removed in order to maintain
a constant “output” f when a unit of xi is added and all other “inputs” are unchanged. The
change in input j (i) is vj (vi ) and all other inputs are unchanged, so vk = 0 for k 6= i, j.
More formally we have
∂f ∂f
0= (a)vi + (a)vj or
∂xi ∂xj
∂f (a)
∂xj vj ∂xi
= = − ∂f (a) .
∂xi x=a vi
f (x)≡f (a) ∂xj
(vii) The direction of change of inputs x which most increases output f (x) starting at a is the
direction ∇f (a).
(b) To find the equation of the tangent plane to the graph of u at (1, 1, 1), u(1, 1, 1) = 3, so
x−1
u(x, y, z) − 3 = ∇u(1, 1, 1) · y − 1 = (x − 1) + (y − 1) + 2(z − 1).
z−1
(c) The equation of the tangent plane to the level curve corresponding to u ≡ 3 at (1, 1, 1) is
x−1
0 = ∇u(1, 1, 1) · y − 1 = (x − 1) + (y − 1) + 2(z − 1).
z−1
(d) To find the marginal rate of substitution of x for z along the level curve u ≡ 3 at the point
(1, 1, 1),
∂u ∂u
0= (1, 1, 1)∆x + (1, 1, 1)∆z
∂x ∂z
so ∂u(1,1,1)
∂z ∂x 1
= − ∂u(1,1,1)
=− .
∂x (x,y,z)=(1,1,1)
2
u=3 ∂z
18
A. W. Richter 1.2. MULTIVARIATE CALCULUS
(e) To find the direction of change of inputs which yields the largest increase in output u at
(1, 1, 1), the direction is
1
1 1
∇u(1, 1, 1) = √ 1 ,
k∇u(1, 1, 1)k 6 2
19
A. W. Richter 1.2. MULTIVARIATE CALCULUS
Example 1.2.3. A firm uses capital, K, and labor, L, to produce output, Y , according to Y =
K α Lβ , where α and β are positive constants.
(a) What is the equation of the isoquant corresponding to a level of output equal to one?
Solution: When Y = 1, the corresponding level curve is K α Lβ = 1.
(b) What is the equation of the tangent plane to the production surface at K = L = 1?
Solution:
∂Y ∂Y
Y = Y (1, 1) + (1, 1)(K − 1) + (1, 1)(L − 1)
∂K ∂L
= 1 + α(K − 1) + β(L − 1)
(c) What is the equation of the tangent “plane” to the level curve corresponding to Y = 1 at
K = L = 1?
Solution:
0 = α(K − 1) + β(L − 1) → αK + βL = α + β.
(d) For small changes along the level curve, starting at K = L = 1, how many units of labor are
needed to replace each unit of capital (i.e. What is the MRSL→K )?
Solution: The Equation of the tangent plane to the level curve at (K, L) = (1, 1) is
0 = (∆K, ∆L) · (fK (1, 1), fL (1, 1)) = (∆K, ∆L) · (α, β),
(e) If it were possible to increase K and L slightly in any proportion so that ||(∆K, ∆L)|| = c,
where c is a very small positive number, what change in K and L would lead to the greatest
increase in output?
Solution: Equation of the tangent plane to the production surface at (1, 1) is
Thus, the maximum value of Y is attained when the above dot product is maximized. This
will occur when
c
(∆K, ∆L) = p (α, β).
α + β2
2
20
A. W. Richter 1.2. MULTIVARIATE CALCULUS
Definition 1.2.8 (Chain Rule–General Version). Suppose that u is a differentiable function of the n
variables x1 , . . . , xn and each xj is a differentiable function of the m variables t1 , . . . , tm . Then u
is a function of t1 , . . . , tm and
∂u ∂u ∂x1 ∂u ∂x2 ∂u ∂xn
= + + ··· +
∂ti ∂x1 ∂ti ∂x2 ∂ti ∂xn ∂ti
for each i = 1, 2, . . . , m.
Example 1.2.4. Given z = f (x, y, t) = yex + t2 where x = g(y, t) = ln(y + t) and y = h(t) =
t3 − 9, use the chain rule to find the total effect of a change in t on z (dz/dt) at t = 2.
Solution:
dz ∂z ∂x dy ∂x ∂z dy ∂z
= + + +
dt ∂x ∂y dt ∂t ∂y dt ∂t
1 1
= yex 3t2 + + 3ex t2 + 2t.
y+t y+t
When t = 2, y = −1 and x = 0, so
dz
= −1(12 + 1) + 12 + 4 = 3.
dt t=2
Example 1.2.5. Assume the following functional forms:
z = x2 + xy 3 , x = uv 2 + w3 , y = u + vew .
Use the chain rule to find the indicated partial derivatives at the given point:
∂z ∂z ∂z
, , , when u = 2, v = 1, w = 0.
∂u ∂v ∂w
Solution: First note that u = 2, v = 1, w = 0 implies x = 2, y = 3
∂z ∂z ∂x ∂z ∂y ∂z ∂z ∂x ∂z ∂y
= + = +
∂u ∂x ∂u ∂y ∂u ∂v ∂x ∂v ∂y ∂v
= (2x + y 3 )(v 2 ) + (3xy 2 ) = (2x + y )(2uv) + (3xy 2 )(ew )
3
= 85, = 178,
∂z ∂z ∂x ∂z ∂y
= +
∂w ∂x ∂w ∂y ∂w
= (2x + y 3 )(3w2 ) + (3xy 2 )(vew )
= 54.
for all t, where n is a positive integer and f has continuous second order partial derivatives.
(a) Verify that f (x, y) = x2 y + 2xy 2 + 5y 3 is homogeneous of degree 3.
Solution: Applying the above definition
21
A. W. Richter 1.2. MULTIVARIATE CALCULUS
Solution: To see this clearly, rewrite the above definition in the following way:
∂2f ∂2f
(x) = (x).
∂xi ∂xj ∂xj ∂xi
Example 1.2.7. Consider the general Cobb-Douglas production function Q = kxa y b . Then,
∂Q ∂Q
= akxa−1 y b , = bkxa y b−1 ,
∂x ∂y
∂2Q ∂2Q
= abkxa−1 y b−1 = .
∂x∂y ∂y∂x
We can continue taking higher order derivatives, and Young’s theorem holds for these cases.
For example, if we take an x1 x2 x4 derivative of order three, then the order of differentiation does
not matter for a C 3 -function. We can keep going and define kth order partial derivatives and C k
functions. For C k functions, the order you take the kth partial derivatives does not matter.
∂ 2 f (x, y) ∂ 2 f (x, y) 2
2 ∂ f (x, y)
x2 + 2xy + y = n(n − 1)f (x, y).
∂x2 ∂x∂y ∂y 2
Solution: To obtain the desired result, differentiate the following result with respect to t
22
A. W. Richter 1.2. MULTIVARIATE CALCULUS
f (a, b) = tn f (x, y)
Example 1.2.10. Verify that partial derivative of f (x, y) = x2 y + 2xy 2 + 5y 3 with respect to x is
homogeneous of degree 2.
Solution: The partial derivative is given by
fx (x, y) = 2xy + 2y 2 .
n2
If all second-order partial derivatives of f exist and are continuous functions of (x1 , x2 , . . . , xn ),
we say that f is twice continuously differentiable or C 2 .
23
A. W. Richter 1.3. BASIC ANALYSIS
The Bordered Hessian matrix of f is a square matrix of second-order partial derivatives that
is bordered by first-order partial derivatives. Given the real-valued function f (x1 , x2 , . . . , xn ), the
bordered Hessian matrix of f is the matrix
∂f ∂f ∂f
0 ∂x1 ∂x2 · · · ∂xn
∂f 2
∂ f 2
∂ f 2
∂ f
∂x
1 ∂x21 ∂x1 ∂x2 · · · ∂x1 ∂xn
∂f 2
∂ f 2
∂ f 2
∂ f
H(f ) = ∂x2 ∂x2 ∂x1
∂x22
· · · ∂x2 ∂xn
. .. .. .. ..
.. . . . .
∂f ∂2f ∂2f ∂2f
∂xn ∂xn ∂x1 ∂xn ∂x2 · · · ∂x2 n
The importance of the Hessian and Bordered Hessian matrices will become clear in later sections.
24
A. W. Richter 1.3. BASIC ANALYSIS
Since the format of a proof using mathematical induction always consists of the same two steps
(establishing the basis for induction and verifying the induction step), it is common practice to
reduce some of the formalism by omitting explicit reference to the statement P (n). It is also ac-
ceptable to omit identifying the steps by name.
Since m and k are natural numbers, so is 7m + 4k . Thus 7k+1 − 4k+1 is also a multiple of 3, and
by induction we conclude that 7n − 4n is a multiple of 3 for all n ∈ N.
In the above example, we have added and subtracted the term 7 · 4k . Where did it come from?
We want somehow to use the induction hypothesis 7k − 4k = 3m, so we break 7k+1 apart into 7 · 7k .
We would like to have 7k − 4k = 3m as a factor instead of just 7k , but to do this we must subtract
(and add) the term 7 · 4k .
1
12 + 22 + · · · + k2 + (k + 1)2 = k(k + 1)(2k + 1) + (k + 1)2
6
1
= [2k3 + 3k2 + k] + k2 + 2k + 1
6
1
= [2k3 + 9k2 + 13k + 6]
6
1
= [(k + 1)(k + 2)(2k + 3)].
6
Thus, the above statement holds for n = k + 1 whenever it holds for n = k, and by principle of
mathematical induction, we conclude that the statement is true for all n.
1 1 1 1 k
+ + + ··· + 2 = .
3 15 35 4k − 1 2k + 1
25
A. W. Richter 1.3. BASIC ANALYSIS
1
Adding 4(k+1)2 −1
to both sides, it follows that
1 1 1 1 1 k 1
+ + + ··· + 2 + 2
= +
3 15 35 4k − 1 4(k + 1) − 1 2k + 1 4(k + 1)2 − 1
k 1
= + 2
2k + 1 4k + 8k + 3
k 1
= +
2k + 1 (2k + 1)(2k + 3)
2k2 + 3k + 1
=
(2k + 1)(2k + 3)
k+1
= .
2k + 3
Thus, the above statement holds for n = k + 1 whenever it holds for n = k, and by principle of
mathematical induction, we conclude that the statement is true for all n.
Example 1.3.5.
(a) Let S be the open interval (0, 5) and let x ∈ S. If ε = min{x, 5 − x}, then we claim that
N (x; ε) ⊆ S. Indeed, for all y ∈ N (x; ε) we have |y − x| < ε, so that
−x ≤ −ε < y − x < ε ≤ 5 − x
Thus 0 < y < 5 and y ∈ S. That is, for some arbitrary point in S, there exists a neighborhood
that is completely contained in S. It follows that every point in S is an interior point of S
(S ⊆ int S). Since the inclusion int S ⊆ S always holds, we have S = int S.
The point 0 is not a member of S, but every neighborhood of 0 will contain positive numbers
in S. Thus 0 is a boundary point of S. Similarly, 5 ∈ bd S and, in fact, bd S = {0, 5}. Note
that none of the boundary of S is contained in S. Of course, there is nothing special about the
open interval (0, 5) in this example. Similar comments would apply to any open interval.
26
A. W. Richter 1.3. BASIC ANALYSIS
Boundary
Point
Interior
Point
(b) Let S be the closed interval [0, 5]. The point 0 is still a boundary point of S, since every
neighborhood of x will contain negative numbers not in S. We have int S = (0, 5) and
bd S = {0, 5}. This time S contains all of its boundary points, and the same could be said of
any other closed interval.
(c) Let S be the interval [0, 5). Then again int S = (0, 5) and bd S = {0, 5}. We see that S
contains some of its boundary, but not all of it.
(d) Let S be the interval [2, ∞). Then int S = (2, ∞) and bd S = {2}. Note that there is no
“point” at ∞ to be included as a boundary point at the right end.
Definition 1.3.4 (Open and Closed Sets). Let S ⊆ R. If bd S ⊆ S, then S is said to be closed. If
bd S ⊆ R\S, then S is said to be open.
Theorem 1.3.2.
(a) A set S is open iff S = int S. Thus, S is open iff every point in S is an interior point of S.
Example 1.3.8. Classify each of the following sets, S, as open, closed, neither, or both.
27
A. W. Richter 1.3. BASIC ANALYSIS
x : |x − 5| ≤ 12
1. N 4.
Solution: Not open: int S = ∅ 6= S, Solution: Not open: int S = (4.5, 5.5) 6=
Closed: bd S = S S, Closed: bd S = {4.5, 5.5} ∈ S
2. Q 5. {x : x2 > 0}
Solution: Neither, int S = ∅ 6= S Solution: Open: int S = R\{0} = S,
and bd S = R * S Not Closed: bd S = {0} 6∈ S
3. ∞ 1
6. n1 : n ∈ N
T
n=1 0, n
Solution: Both: (0, 1) ∩ (0, 12 ) ∩ · · · = ∅ Solution: Neither: int S = ∅ 6= S
(int ∅ = ∅ and bd ∅ = ∅) and bd S * S (0 6∈ S).
Example 1.3.9. True/False: If S ⊆ R2 is an open set, then f (x, y) = x + y cannot have a global
maximizer subject to (x, y) ∈ S.
Solution: True. Suppose (x∗ , y ∗ ) is a global maximizer of x + y subject to (x, y) ∈ S. Since
S is open, we can always find ε small enough so that an open ε-neighborhood around (x∗ , y ∗ ) is
entirely contained in S. But then a point such as (x∗ + ε/2, y ∗ + ε/2) would lie in S with a value
of the objective function (x∗ + ε/2) + (y ∗ + ε/2); i.e.,
ε ε
f x∗ + , y ∗ + = x∗ + y ∗ + ε > x∗ + y ∗ = f (x∗ , y ∗ ).
2 2
Thus, (x∗ , y ∗ ) cannot be a global maximizer, and our original assumption must be false. Specifi-
cally, ∄(x, y) ∈ S that is a global maximizer of f (x, y) = x + y.
Theorem 1.3.3.
(b) The intersection of any finite collection of open sets is an open set.
Proof.
(b) First note that this result does not hold for infinite collections of sets. To see why, notice that
T∞each n ∈ N, if we define An = (−1/n, 1/n), then each An is an open set. However,
for
n=1 An = {0}, which is not open. Thus we see that we cannot generalize the above result
to the intersection of an infinite collection of open sets.
Consider the finite case. Define S := ni=1 Gi , where G = {G1 , G2 , . . . , Gn } is an arbitrary
T
collection of open sets. If S = ∅, we are done, since ∅ is open (and closed) (int ∅ = ∅). If
S 6= ∅, let x ∈ S. Then x ∈ Gi for all i = 1, 2, . . . , n. Since each Gi is open, there exist
neighborhoods Ni (x; εi ) of x such that Ni (x; εi ) ⊆ Gi . Let ε = min{ε1 , . . . , εn }. Then
N (x; ε) ⊆ Gi for each i = 1, . . . , n, so N (x; ε) ⊆ S. Thus x is an interior point of S, and S
is open.
28
A. W. Richter 1.3. BASIC ANALYSIS
Corollary 1.3.1.
(a) The intersection of any collection of closed sets is a closed set.
(b) The union of any finite collection of closed sets is a closed set.
Proof.
(a) To prove the above result, define T := ∞
T
i=1 Fi , where F = {F1 , F2 , . . .} is an arbitrary
infinite collection of closed sets. If T = ∅, we are done, since ∅ is closed (and open)
(bd ∅ = ∅ ⊆ ∅). R \ ( ∞
T S∞
i=1 iF ) = i=1 (R \ Fi ) (the complement of the intersection will be
the union of the individual complements, which can be seen using a simple venn diagram).
Thus, we have R \ ( ∞
T
i=1 F i ) equal to the union of open sets, since a set is closed if and only
if its complement is open. Since weThave shown above thatTthe union of any collection of
open sets is open, we have that R \ ( ∞ ∞
i=1 Fi ) is open. Thus, i=1 Fi must be closed.
Sn
(b) To prove the above result, define T := Sn i=1 Fi , where
Tn F = {F1 , F2 , . . . , Fn } is an arbitrary
finite collection of closed sets. R \ ( i=1 Fi ) = i=1 (R \ Fi ) (the complement of the union
will be the intersection of the individualS complements, which again can be seen using a simple
venn diagram). Thus, we have R \ ( ni=1 Fi ) equal to the intersection of open sets. Since we
have shown S above that the intersectionS of any finite collection of open sets is open, we have
that R \ ( ni=1 Fi ) is open. Thus, ni=1 Fi must be closed by definition.
Definition 1.3.6 (Bounded Sequence). A sequence {sn } is said to be bounded if the range {sn :
n ∈ N} is a bounded set. That is if there exists an M ≥ 0 such that |sn | ≤ M for all n ∈ N.
Theorem 1.3.4. Every convergent sequence in Rn is bounded.
Proof. Let sn be a convergent sequence and let sn → s. From the definition of convergence with
ε = 1, we obtain N ∈ R such that |sn − s| < 1 whenever n > N . Thus for n > N the triangle
inequality implies that |sn | < |s| + 1. We know that sn is bounded if the range {sn : n ∈ N} is a
bounded set, that is, if there exists an M ≥ 0 such that |sn | ≤ M for all n ∈ N. Thus, let
29
A. W. Richter 1.3. BASIC ANALYSIS
(so that M could either be |s| + 1 or the largest absolute value among the first N terms) then we
have |sn | ≤ M for all n ∈ N, so sn is bounded.
Theorem 1.3.5. Suppose that {sn } and {tn } are convergent sequences with lim sn = s and
lim tn = t. Then
(a) lim(sn + tn ) = s + t
(c) lim(sn tn ) = st
Proof.
(a) To show that sn + tn → s + t, we need to make the difference |(sn + tn ) − (s + t)| small.
Using the triangle inequality, we have
Now, given any ε > 0, since sn → s, there exists N1 such that n > N1 implies that |sn − s| <
ε ε
2 . Similarly, since tn → t, there exists N2 such that n > N2 implies that |tn − t| < 2 . Thus
if we let N = max{N1 , N2 }, then n > N implies that
ε ε
|(sn − s) + (tn − t)| ≤ |sn − s| + |tn − t| < + = ε.
2 2
(b) To show that ksn → ks for any k ∈ R, we need to make the difference |ksn − ks| small. We
have
Now, given any ε > 0, since sn → s, there exists N such that n > N implies that |sn − s| <
ε/|k|. Thus for n > N
ε
|ksn − ks| = |k||sn − s| < |k| = ε.
|k|
30
A. W. Richter 1.3. BASIC ANALYSIS
We know that every convergent sequence is bounded (Theorem 1.3.4). Thus, sn is bounded,
and there exists M1 > 0 such that |sn | ≤ M1 for all n. Letting M = max{M1 , |t|}, we
obtain the inequality
|sn tn − st| ≤ M |tn − t| + M |sn − s|.
ε
Now, given ant ε > 0, there exists N1 and N2 such that |tn − t| < 2M when n > N1 and
ε
|sn − s| < 2M when n > N2 . Let N = max{N1 , N2 }. Then n > N implies that
(d) Since sn /tn = sn (1/tn ), it suffices from part (c) to show that lim 1/tn = 1/t. That is, given
ε > 0, we must show that
1 1 t − tn
− = <ε
tn t tn t
for all n sufficiently large. To get a lower bound on how small the denominator can be, we
note that since t 6= 0, there exists an N1 such that n > N1 implies that |tn − t| < |t|/2. Thus
for n > N1 we have
|t| |t|
|tn | = |t − (t − tn )| ≥ |t| − |t − tn | > |t| − =
2 2
by the reverse triangle inequality. There also exists an N2 such that n > N2 implies that
|tn − t| < 12 ε|t|2 . Let N = max{N1 , N2 }. Then n > N implies that
1
− 1 = t − tn < 2 |tn − t| < ε.
tn t tn t |t|2
Theorem 1.3.6. A set S ⊆ Rn is closed if and only if every convergent sequence in S converges to
a point in S.
Proof. Suppose that S is closed and that xn is a sequence in S converging to x. Suppose that x 6∈ S.
BecauseT S is closed, Rn \S is open, so that for some ε > 0, there exists N (x; ε) ⊆ Rn \S. That is
N (x; ε) S = ∅. Since lim xn = x, there is an N such that xn ∈ N (x; ε), if n > N . But then
xn 6∈ S, whenever n > N , which is impossible since xn is a sequence in S. Thus it must be the
case that if a set S is closed, then every convergent sequence in S must converge to a point in S.
Suppose that STis not closed. Then Rn \S is not open, implying that there exists an x ∈ Rn \S
such that N (x; ε) S 6= ∅, for every ε > 0. Practically, this means that at least part of the
neighborhood of x lies in S. In particular, for every positive integer n, there exists an xn ∈ S, such
that |xn − x| < 1/n. This gives us a sequence {xn } in S converging to a point x not in S. Thus,
by the contra-positive argument (not closed implies there exists a convergent sequence in S that
converges to a point not in S), we have that if every convergent sequence in S converges to a point
in S, then S must be closed.
31
A. W. Richter 1.3. BASIC ANALYSIS
1.3.4 Compactness
Theorem 1.3.7 (Heine-Borel). A subset S of Rn is compact iff S is closed and bounded.
32
Chapter 2
2.1 Determinants
2.1.1 Minors, Cofactors, and Evaluating Determinants
Definition 2.1.1 (Minor and Cofactor of a Matrix). Let A be an n × n matrix. Let Aij be the
(n − 1) × (n − 1) submatrix obtained by deleting row i and column j from A. Then, the scalar
33
A. W. Richter 2.1. DETERMINANTS
(a) Rows and columns can be interchanged without affecting the value of a determinant. That is
|A| = |AT |.
(b) If two rows (or columns) are interchanged, the sign of the determinant changes. For example
3 4
= − 1 −2 .
1 −2 3 4
(c) If a row (or column) is changed by adding to or subtracting from its elements the correspond-
ing elements of any other row (or column), the determinant remains unaltered. For example
3 4 3 + 1 4 − 2 4 2
∼ = = −10.
1 −2 1 −2 1 −2
(d) If the elements in any row (or column) have a common factor α, then the determinant equals
the determinant of the corresponding matrix in which α = 1, multiplied by α. For example
6 8
= 2 3 4 = 2 × (−10) = −20.
1 −2 1 −2
(e) When at least one row (or column) of a matrix is a linear combination of the other rows (or
columns), the determinant is zero. Conversely, if the determinant is zero, then at least one
row and/or one column are linearly dependent on the other rows and columns, respectively.
For example the
3 2 1
det 1 2 −1
2 −1 3
is zero because the first column is a linear combination of the second and third columns
Similarly, there is a linear dependence between the rows, which is given by the relation
7 4
row 1 = row 2 + row 3.
8 5
(f) The determinant of an upper triangular or lower triangular matrix is the product of the main
diagonal entries. For example:
3 2 1
0 2 −1 = 3 × 2 × 4 = 24.
0 0 4
(g) The determinant of the product of two square matrices is the product of the individual deter-
minants. Fore example
|AB| = |A||B|.
This rule can be generalized to any number of factors. One immediate application is to matrix
powers: |A2 | = |A||A| = |A|2 , and more generally |An | = |A|n for integer n.
34
A. W. Richter 2.2. INVERSES OF MATRICES
Ax = y,
where x and y are n × 1 column vectors. Premultiplying both sides by A−1 we get the inverse
relationship
x = A−1 y.
More generally, consider the matrix equation for multiple (m) right-hand sides:
A X = Y ,
n×n n×m n×m
X = A−1 Y.
35
A. W. Richter 2.2. INVERSES OF MATRICES
Cji
bij = ,
|A|
where bij denotes the entries of A−1 . Cij is defined as (i, j)th cofactor of A or, more precisely,
the determinant of the submatrix of order (n − 1) × (n − 1) obtained by deleting the ith row and
j th column of A, multiplied by (−1)i+j . The n × n matrix whose (i, j)th entry is Cji is called the
adjugate (or classical adjoint) of A and is written adj A.1
This direct inversion procedure is useful only for small matrix orders: 2 or 3. In the examples
below the inversion formulas for second and third order matrices are listed.
where
a a23 a a13 a a13
C11 = 22 , C21 = − 12 , C31 = 12 ,
a32 a33 a32 a33 a22 a23
a21 a23 a11 a13 a11 a13
C12 = − , C22 = , C32 = − ,
a31 a33 a31 a33 a21 a23
a a22 a a12 a a12
C13 = 21 , C23 = − 11 , C33 = 11 ,
a31 a32 a31 a32 a21 a22
and |A| is given by equation 2.2 on page 33.
Example 2.2.3.
2 4 2 1 −4 2
1
A = 3 1 1 , A−1 =− −2 0 4 .
8
1 0 1 −1 4 −10
If the order exceeds 3, the general inversion formula becomes rapidly useless as it displays combi-
natorial complexity.
1
The adjugate has sometimes been called the “adjoint”, but that terminology is ambiguous. Today, “adjoint” of a
matrix normally refers to its conjugate transpose.
36
A. W. Richter 2.3. QUADRATIC FORMS AND DEFINITENESS
(a) The inverse of the transpose is equal to the transpose of the inverse. That is
because
AA−1 = (AA−1 )T = (A−1 )T AT = I.
(b) The inverse of a symmetric matrix is also symmetric. Given the previous rule, (AT )−1 =
A−1 = (A−1 )T , hence A−1 is also symmetric.
(c) The inverse of a matrix product is the reverse product of the inverses of the factors. That is
(AB)−1 = B −1 A−1 .
(d) For a diagonal matrix D in which all diagonal entries are nonzero, D −1 is again a diagonal
matrix with entries 1/dii .
0 0 0 . . . Snn −1
(f) The inverse of an upper (lower) triangular matrix is also an upper (lower) triangular matrix.
Q(x) = xT Ax
37
A. W. Richter 2.3. QUADRATIC FORMS AND DEFINITENESS
Note that
a11 . . . a1n x1
a21 . . . a2n x2
Q(x) = xT Ax = x1 x2 . . . xn .
. . . .. ..
. . . .
an1 . . . ann xn
Pn
Pi=1 a1i xi
n
i=1 a2i xi
= x1 x2 . . . xn
..
Pn .
i=1 ani xi
X
= aij xi xj .
i,j
(e) indefinite if xT Ax > 0 for some x in Rn and < 0 for some other x in Rn .
38
A. W. Richter 2.3. QUADRATIC FORMS AND DEFINITENESS
Note that for any real vector x 6= 0, Q(x) will be positive, since the square of any number is positive,
the coefficients of the squared terms are positive, and the sum of positive numbers is always positive.
Thus, A is positive definite.
Example 2.3.2. Now consider an alternative 3 × 3 matrix given by
−2 1 0
D = 1 −2 0 .
0 0 −2
The general quadratic form is given by
−2 1 0 x1
Q(x) = xT Ax = x1 x2 x3 1 −2 0 x2
0 0 −2 x3
x1
= −2x1 + x2 x1 − 2x2 −2x3 x2
x3
= −2x21 + 2x1 x2 − 2x22 − 2x23
= −2x21 − 2[x22 − x1 x2 ] − 2x23 .
Note that Q(x) will be negative if x1 and x2 are of opposite sign or equal to one another. Now
consider the case where |x1 | > |x2 |. Write Q(x) as
The first, third, and fourth terms are clearly negative. With |x1 | > |x2 |, |2x21 | > |2x1 x2 |, the first
term is more negative than the second term is positive, and hence the whole expression is negative.
Now consider the case where |x1 | < |x2 |. The first, third, and fourth terms are still negative. But,
with |x1 | < |x2 |, |2x22 | > |2x1 x2 | so that the third term is more negative than the second term is
positive, and so the whole expression is negative. Thus, this quadratic form is negative definite for
all real values of x 6= 0.
Remark 2.3.1. A matrix that is positive (negative) definite is positive (negative) semi-definite.
The definiteness of a matrix plays an important role. For example, for a function f (x) of one
variable, the sign of the second derivative f ′′ (x0 ) at a critical point x0 gives a sufficient condition
for determining whether x0 is a maximum, minimum, or neither (proposition 1.1.2). This test gen-
eralizes to more than one variable using the definiteness of the Hessian matrix H (more on this
when we get to optimization). There is a convenient way to test for the definiteness of a matrix, but
before we can formulate this test we first need to define the concept of principal minors of a matrix.
39
A. W. Richter 2.3. QUADRATIC FORMS AND DEFINITENESS
Definition 2.3.4 (Leading Principle Submatrix and Leading Principle Minor). Let A be an n × n
matrix. The kth order principle submatrix of A obtained by deleting the last n − k rows and the last
n − k columns from A is called the kth order leading principle submatrix of A. Its determinant,
denoted ∆k , is called the kth order leading principle minor of A.
there is one third order principle minor: B3 = det(A) = |A|. There are three second order principle
minors:
(1) a11 a12
1. B2 = , formed by deleting column 3 and row 3 from A;
a21 a22
(2) a11 a13
2. B2 = , formed by deleting column 2 and row 2 from A;
a31 a33
(3) a22 a23
3. B2 = , formed by deleting column 1 and row 1 from A.
a32 a33
There are three first order principle minors:
(1)
1. B1 = a11 , formed by deleting the last 2 rows and columns;
(2)
2. B1 = a22 , formed by deleting the first and third columns and rows, and
(3)
3. B1 = a33 , formed by deleting the first 2 rows and columns.
(a) A is positive definite if and only if ∆k > 0 for k = 1, 2, . . . , n (all of its leading principle
minors are strictly positive);
(b) A is negative definite if and only if (−1)k ∆k > 0 for k = 1, 2, . . . , n (every leading principle
minor of odd order is strictly negative and every leading principle minor of even order is
strictly positive);
(c) If some kth order principle minor of A (or some pair of them) is nonzero but does not fit either
of the above two sign patterns, then A is indefinite.
40
A. W. Richter 2.3. QUADRATIC FORMS AND DEFINITENESS
Remark 2.3.2. Given an n × n symmetric matrix, the conditions ∆k ≥ 0 for positive semi-
definiteness and (−1)k ∆k ≥ 0 for negative semi-definiteness are necessary conditions but not
sufficient conditions. To see this, consider the following 2 × 2 symmetric matrix.
0 0
A= .
0 −4
For this matrix, ∆1 = 0 and ∆2 = 0. Thus, looking only at leading principle minors, one could
falsely conclude that this matrix is positive semi-definite. However, we have to check all principle
minors to deduce the correct form of semi-definiteness. If one checks all principle minors then
(1)
B1 = 0 ≤ 0,
(2)
B1 = −4 ≤ 0,
0 0
B2 =
= 0 ≥ 0,
0 −4
which violates the definition of positive semi-definiteness. In fact, this matrix is negative semi-
definite.
(f) ∆1 > 0, ∆2 = 0, ∆3 > 0, ∆4 > 0 → A is not definite, not negative semi-definite but might
be positive semi-definite. However, to establish this, we must to check all 15 principle minors
Bk , for k = 1, 2, 3, 4.
Example 2.3.5. Use the above theorems to check the following matrices for definiteness:
(a)
1 2
A= , ∆1 = 1 > 0, ∆2 = −3 < 0
2 1
Thus A is indefinite.
41
A. W. Richter 2.4. EIGENVALUES AND EIGENVECTORS
(b)
−1 1 (1) (2)
A= , B1 = B1 = −1 < 0, ∆2 = 0
1 −1
Thus A is negative semi-definite.
(c)
1 0 −1
B = 0 3 0 , ∆1 = 1 > 0, ∆2 = 3 > 0, ∆3 = 9 > 0
−1 0 4
Thus, B is positive definite.
(d)
1 3 2 3
3 5 −2 4
B=
2 −2 2 1 ,
∆1 = 1 > 0, ∆2 = −4 < 0, ∆3 = −56 < 0, ∆4 = 142 > 0
3 4 1 2
Thus B is indefinite.
42
A. W. Richter 2.4. EIGENVALUES AND EIGENVECTORS
P (λ) = λn + α1 λn−1 + · · · + αn = 0.
This is known as the characteristic equation of the matrix A. The left-hand side is known as the
characteristic polynomial. We know that a polynomial of degree n has n (generally complex) roots
λ1 , λ2 , . . . , λn . These n numbers are called eigenvalues, eigenroots, or characteristic values of the
matrix A. The following theorem summarizes the results.
Theorem 2.4.1. The number λ is an eigenvalue of the n × n matrix A if and only if λ satisfies the
characteristic equation
|A − λI| = 0.
Axi = λxi .
and we can solve for the associated eigenvalues by setting the above characteristic equation to zero
and solving for λ.
Example 2.4.1. Find the eigenvalues and associated eigenvectors of the matrix
5 7
A= .
−2 −4
Solution:
5−λ 7
A − λI = (2.4)
−2 −4 − λ
so the characteristic equation of A is
set
0 = det(A − λI)
5 − λ 7
=
−2 −4 − λ
= (5 − λ)(−4 − λ) − (−2)(7)
= λ2 − λ − 6
= (λ + 2)(λ − 3).
Thus, the matrix A has two eigenvalues −2 and 3. To distinguish them, we write λ1 = −2 and
λ2 = 3. To find the associated eigenvectors, we must separately substitute each eigenvalue into
(2.4) and then solve the resulting system (A − λI)x = 0.
43
A. W. Richter 2.4. EIGENVALUES AND EIGENVECTORS
T
Case 1: λ1 = −2. With x = x1 x2 , the system (A − λI)x = 0 is
7 7 x1 0
= .
−2 −2 x2 0
Each of the two scalar equations here is a multiple of the equation x1 + x2 = 0, and any nontrivial
T T
solution x = x1 x2 of this equation is a nonzero multiple of 1 −1 . Hence, to within a
T
constant multiple, the only eigenvector associated with λ1 = −2 is x1 = 1 −1 .
T
Case 2: λ1 = 3. With x = x1 x2 , the system (A − λI)x = 0 is
2 7 x1 0
= .
−2 −7 x2 0
Again, we have only a single equation, 2x1 + 7x2 = 0, and any nontrivial solution of this equa-
tion will suffice. The choice x2 = −2 yields x1 = 7 so (to within a constant multiple) the only
T
eigenvector associated with λ2 = 3 is x2 = 7 −2 .
Example 2.4.2. Find the eigenvalues and eigenvectors of the 2 × 2 matrix
−1 3
A= .
2 0
set
Solution: Given A, |A − λI| = λ2 + λ − 6 = 0. It follows that λ1 = −3 and λ2 = 2. We can
then set up the equation Axi = λi xi or (A − λi I)xi = 0 for i ∈ {1, 2}. Specifically, we have
2 3 x1 0 −3 3 x1 0
= and = ,
2 3 x2 0 2 −2 x2 0
which results in
3 1
x1 = and x2 =
−2 1
or a multiple thereof.
Theorem 2.4.2. Let A be a k × k matrix with eigenvalues λ1 , . . . , λk . Then,
(a) λ1 + λ2 + ... + λk = tr(A) and
(b) λ1 · λ2 · · · λk = det (A).
a11 a12
Proof. Consider the case for 2 × 2 matrices. Assume A = , then |A − λI| = λ2 −
a21 a22
(a11 + a22 )λ + (a11 a22 − a12 a21 ), so that the characteristic equation is
p(λ) = λ2 − (a11 + a22 )λ + (a11 a22 − a12 a21 )
= λ2 − tr(A)λ + det(A).
If λ1 and λ2 are the roots of this polynomial, we can rewrite the solution as
p(λ) = β(λ1 − λ)(λ2 − λ)
= βλ2 − β(λ1 + λ2 )λ + βλ1 λ2 ,
for some constant β. Comparing the two expressions, we find that β = 1 and
tr(A) = β(λ1 + λ2 ),
det(A) = βλ1 λ2 .
The theorem naturally extends to higher dimensional cases.
44
A. W. Richter 2.4. EIGENVALUES AND EIGENVECTORS
(b) A square matrix A is singular if and only if 0 is an eigenvalue. However, the zero vector
cannot be an eigenvector.
(c) A matrix is invertible if and only if none of its eigenvalues is equal to zero.
(d) Every eigenvalue has an infinite number of eigenvectors associated with it, as any nonzero
scalar multiple of an eigenvector is also an eigenvector.
(f) The eigenvalues of A and AT are the same (as their characteristic polynomials are the same),
but there is no simple relationship between their eigenvectors.
(g) The eigenvalues of a shifted matrix A − αI are λ − α and the eigenvectors are the same as
those of A since
Ax = λx ⇒ (A − αI)x = (λ − α)x.
(h) The eigenvalues of A−1 are 1/λ and the eigenvectors are the same as those of A since
Ax = λx ⇒ A−1 x = λ−1 x.
(i) The eigenvalues of A2 are λ2 , and the eigenvectors are the same as those of A since
Ax = λx ⇒ A2 x = A(λx) = λ(Ax) = λ2 x.
This property naturally generalizes to high-order powers (i.e. the eigenvalues of Ak are λk ).
(j) Let x be the eigenvector of an n × n matrix A with eigenvalue λ. Then, the eigenvalue of
αk Ak + αk−1 Ak−1 + · · · + α1 A + α0 I
αk λk + αk−1 λk−1 + · · · + α1 λ + α0 ,
(a) A is positive definite if and only if all of its eigenvalues are strictly positive (λi > 0 ∀i)
(b) A is negative definite if and only if all of its eigenvalues are strictly negative (λi < 0 ∀i)
(c) A is positive semi-definite if and only if all of its eigenvalues are nonnegative (λi ≥ 0 ∀i)
(d) A is negative semi-definite if and only if all of its eigenvalues are nonpositive (λi ≤ 0 ∀i)
45
A. W. Richter 2.4. EIGENVALUES AND EIGENVECTORS
1 2
Example 2.4.3. If A = , then ∆1 = 1 > 0, ∆2 = −3 < 0, which implies that A is
2 1
indefinite. Or, using eigenvalues, we can solve
1 2
det − λI = λ2 − 2λ − 3 = (λ − 3)(λ + 1) = 0
2 1
which implies λ1 = 0 and λ2 = −2 and fulfills the condition for negative semi-definiteness.
1 0 −1
1 0
Example 2.4.5. If A = 0 3 0 , then ∆1 = 1 > 0, ∆2 = = 3 > 0 and ∆3 = 9 > 0,
0 3
−1 0 4
which implies A is positive definite. Alternatively, using eigenvalues,
1−λ 0 −1
det 0 3−λ 0 = (3 − λ)(λ2 − 5λ + 3) = 0,
−1 0 4−λ
√
which implies λ1,2 = (5 ± 13)/2 and fulfills the condition for positive definiteness.
46
Chapter 3
In this chapter, we will look at linear algebra more from a transformational viewpoint than system
of equations viewpoint (e.g., function from one vector space to another).
(VS 6) For each pair of elements a, b in F and each element x in V , (ab)x = a(bx).
(VS 7) For each element a in F and each pair of elements x, y in V , a(x + y) = ax + ay.
The elements of the field F are called scalars and the elements of the vector space V are called
vectors. The following are examples of vector spaces:
1. An object of the form (a1 , a2 , . . . , an ), where the entries ai , i = {1, . . . , n} are elements
from a field F , is called an n-tuple with entries from F . The set of all n-tuples with entries
from F is a vector space, denoted F n , under the operations of coordinate-wise addition and
scalar multiplication; that is, if u = (a1 , a2 , . . . , an ) ∈ F n , v = (b1 , b2 , . . . , bn ) ∈ F n and
c ∈ F , then
47
A. W. Richter 3.1. VECTOR SPACES AND SUBSPACES
2. The set of all m×n matrices with entries from a field F is a vector space, denoted Mm×n (F ),
under the following operations of addition and scalar multiplication: For A, B ∈ Mm×n (F )
and c ∈ F
for 1 ≤ i ≤ m and 1 ≤ j ≤ n.
3. Let S be any nonempty set and F be any field, and let F(S, F ) denote the set of functions
from S to F . The set F(S, F ) is a vector space under the operations of addition and scalar
multiplication defined for f, g ∈ F(S, F ) and c ∈ F by
for each s ∈ S. Note that these are the familiar operations of addition and scalar multiplica-
tion for the functions used in algebra and calculus.
Since (VS 1), (VS 2), and (VS 8) fail to hold, S is not a vector space under these operations.
Theorem 3.1.1 (Cancelation Law for Vector Addition). If x, y, and z are elements of a vector space
V such that x + z = y + z, then x = y.
x = x + 0 = x + (z + v) = (x + z) + v
= (y + z) + v = y + (z + v)
=y+0=y
Proof. By (VS 3), there exists 0 ∈ V such that x + 0 = x, for all x ∈ V . Now suppose that
x + z = x, for all x ∈ V . Then x + 0 = x + z, which implies 0 = z by Theorem 3.1.1.
Proof. Suppose x + z = 0 = x + y. Then z = y by Theorem 3.1.1. This means there is only one
additive inverse to x. Thus, we say y = −x.
Theorem 3.1.2. In any vector space V the following statements are true:
(a) 0x = 0 for each x ∈ V .
48
A. W. Richter 3.2. LINEAR COMBINATIONS AND SPANNING CONDITIONS
(a) 0 ∈ W .
This theorem provides a simple method for determining whether or not a given subset of a vector
space is a subspace. The following are examples of subspaces:
for any scalar c. Hence A + B and cA are diagonal matrices for any scalar c and also belong
to W . Therefore the set of diagonal matrices, W , is a subspace of V = Mn×n (F ).
2. The transpose At of an m×n matrix A is the n×m matrix obtained from A by interchanging
the rows with the columns; that is, (At )ij = Aji . A symmetric matrix is a matrix A such that
At = A. Define V = Mn×n (R) and W = {A ∈ Mn×n : At = A}. The zero matrix is equal
to its transpose and hence belongs to W . Let A, B ∈ W and c ∈ R, then
4. The set of matrices in Mn×n (R) having nonnegative entries is not a subspace of Mn×n (R)
because it is not closed under scalar multiplication.
49
A. W. Richter 3.2. LINEAR COMBINATIONS AND SPANNING CONDITIONS
Equating coefficients
1 3 2 1 0 −4
−2 −5 −2 rref 0 1 2
−5 −4 12 ∼ 0
0 0
−3 −9 −6 0 0 0
Thus a = −4 and b = 2, which proves that 2x3 − 2x2 + 12x − 6 is a linear combination of
x3 − 2x2 − 5x − 3 and 3x3 − 5x2 − 4x − 9.
Definition 3.2.2 (Span). Let S be a nonempty subset of a vector space V . The span of S, denoted
span(S), is the set consisting of all linear combinations of the elements of S. For convenience we
define span(∅) = {0}.
In R3 , for instance, the span of the set {(1, 0, 0), (0, 1, 0)} consists of all vectors in R3 that
have the form a(1, 0, 0) + b(0, 1, 0) = (a, b, 0) for some scalars a and b. Thus the span of
{(1, 0, 0), (0, 1, 0)} contains all the points in the xy-plane. In this case, the span of the set is a
subspace of R3 . This fact is true in general.
Theorem 3.2.1. The span of any subset S of a vector space V is a subspace of V . Moreover, any
subspace of V that contains S must also contain the span of S.
Proof. This result is immediate if S = ∅ because span(∅) = {0}, which is a subspace that is
contained in any subspace of V .
If S 6= ∅, then S contains an element z and 0z = 0 is an element of span(S). Let x, y ∈
span(S). Then there exists elements u1 , u2 , . . . , um , v1 , v2 . . . , vn in S and scalars a1 , a2 , . . . , am
and b1 , b2 , . . . , bn such that
x = a1 u1 + a2 u2 + · · · + am um and y = b1 v1 + b2 v2 + · · · + bm vm .
Then
x + y = a1 u1 + a2 u2 + · · · + am um + b1 v1 + b2 v2 + · · · + bm vm
are clearly linear combinations of the elements of S; so x + y and cx are elements of span(S).
Thus, span(S) is a subspace of V .
Now let W denote any subspace of V that contains S. If w ∈ span(S), then w has the form w =
c1 w1 +c2 w2 +· · ·+ck wk for some elements w1 , w2 , . . . , wk ∈ W and some scalars c1 , c2 , . . . , ck .
Since S ⊆ W , we have w1 , w2 , . . . , wk ∈ W in W . Therefore w = c1 w1 + c2 w2 + · · · + ck wk is
an element of W (since W is a subspace of V , it is closed under addition and scalar multiplication).
Since w, an arbitrary element of span(S), belongs to W , it follows that span(S) ⊆ W .
50
A. W. Richter 3.2. LINEAR COMBINATIONS AND SPANNING CONDITIONS
Example 3.2.2. Suppose V = R3 and S = {(0, −2, 2), (1, 3, −1)}. Is (3, 1, 5) ∈ span(S). To
answer this question, try to find constants a and b so that
Equating coefficients
0 1 3 1 0 4
−2 3 1 rref
∼ 0 1 3
2 −1 5 0 0 0
p(x) = a0 + a1 x + · · · + an xn
= a0 (1) + a1 (x) + · · · + an (xn ) ∈ span(S).
4. Let V = R3 and S = {(1, 1, 0), (1, 0, 1), (0, 1, 1)}. Let x ∈ span(S). Then for a, b, c ∈ R
51
A. W. Richter 3.3. LINEAR INDEPENDENCE AND LINEAR DEPENDENCE
Equating coefficients
1 1 0 x1 1 0 0 (x1 + x2 − x3 )/2
1 0 1 x2 rref
∼ 0 1 0 (x1 − x2 + x3 )/2 ,
0 1 1 x3 0 0 1 (x2 + x1 + x3 )/2
which is consistent. Hence x ∈ span(S), R3 ⊆ span(S), and S generates R3 .
3. A set is linearly independent if and only if the only representations of 0 as linear combinations
of its elements are trivial representations.
The condition in 3 provides a useful method for determining if a finite set is linearly independent.
This technique is illustrated in the following example.
Example 3.3.1. Determine whether the following sets are linearly dependent or linearly indepen-
dent.
1. In P2 (R), let S = {3 + x + x2 , 2 − x + 5x2 , 4 − 3x2 }. Consider the equation:
a(3 + x + x2 ) + b(2 − x + 5x2 ) + c(4 − 3x2 ) = (3a + 2b + 4c) + (a − b)x + (a + 5b − 3c)x2 = 0
for all x. Equating coefficients
3 2 4 1 0 0
1 −1 0 rref
∼ 0 1 0 .
1 5 −3 0 0 1
Thus, the only solution is a = b = c = 0, which implies that S is linearly independent in
P2 (R).
52
A. W. Richter 3.4. BASES AND DIMENSION
Equating coefficients
1 −2 1 −2
−3 6 rref 0 0
−2 4 ∼ 0 0 ,
4 −8 0 0
which implies that a = 2b. Therefore, there are infinitely many nontrivial solutions and so S
is linearly dependent in M2×2 (R).
v = a1 u1 + a2 u2 + · · · + an un
Proof. First let β be a basis for V . If v ∈ V , then v ∈ span(β) because span(β) = V . Thus v is a
linear combination of the elements in β. Suppose that
v = a1 u1 + a2 u2 + · · · + an un and v = b1 u1 + b2 u2 + · · · + bn un
are two such representations of v. Subtracting the second equation from the first gives
53
A. W. Richter 3.4. BASES AND DIMENSION
Example 3.4.1.
1. Since span(∅) = {0} and ∅ is linearly independent, ∅ is a basis for the vector space {0}.
2. In F n , {e1 , e2 , . . . , en }, where ej denotes a vector whose jth coordinate is 1 and whose other
coordinates are 0 is a basis for F n and is called the standard basis.
3. In Mm×n (F ), let M ij denote the matrix whose only nonzero entry is a 1 in the ith row and
jth column. Then {M ij : 1 ≤ i ≤, 1 ≤ j ≤ n} is a basis for Mm×n (F ).
4. In Pn (F ) the set {1, x, x2 , . . . , xn } is a basis. We call this basis the standard basis for Pn (F ).
Proof. If S = ∅ or S = {0}, then V = {0} and ∅ is a subset of S that is a basis for V . Otherwise
suppose there exists u1 ∈ S and u1 6= 0. Define S1 = {u1 }. Then S1 ⊆ S and S1 is independent.
If span(S1 ) = V , we are done. If not, then there exists a u2 ∈ S where u2 6∈ span(S1 ). Define
S2 = S1 ∪ {u2 }. Then S2 is independent. If span(S2 ) = V , we are done. If not,. . ., then there
a um ∈ S where um 6∈ span(Sm−1 ). Define Sm = Sm−1 ∪ {um }. Then Sm is independent. If
k > m and uk ∈ S, then u ∈ span(Sm ). Thus span(Sm ) = V and Sm ⊆ S is a basis for V .
Equating coefficients
2 1 1 1 1 0 2 0
−1 −1 1 −2 rref
∼ 0 1 −3 0 ,
4 3 −1 1 0 0 0 1
which is consistent for all x. Hence x ∈ span(S) and S generates R3 . By Theorem 3.4.2,
S \ {(1, 1, −1)} is a basis for R3 .
which is consistent. Hence span(S) = M2×2 (R). By Theorem 3.4.2, S \ {S3 } is a basis for
M2×2 (R).
54
A. W. Richter 3.4. BASES AND DIMENSION
Then AX = 0. Since m > n, we will always have at least one free variable and so some xj 6= 0.
Thus S is dependent.
Definition 3.4.2. A vector space is called finite-dimensional if it has a basis consisting of a finite
number of elements. The unique number of elements in each basis for V is called the dimension of
V and is denoted dim(V ). A vector space that is not finite-dimensional is call infinite-dimensional.
Example 3.4.3.
1. dim(F n ) = n
3. dim(Mm×n ) = mn
4. dim({0}) = 0
5. dim(P ) = ∞
Corollary 3.4.1. Let V be a vector space, dim(V ) = n < ∞, and S ⊆ V . If |S| < n, then
span(S) 6= V .
1
For any set S, |S| corresponds to the number of vectors in S.
55
A. W. Richter 3.4. BASES AND DIMENSION
Proof. Suppose span(S) = V . and |s| < n. By Theorem 3.4.2, there exists β ⊆ S so that β is a
basis for V . Then n = |β| ≤ |s| < n, which is a contradiction.
Remark 3.4.1. Let V be a vector space, dim(V ) = n < ∞, and S ⊆ V . Then the these results
directly follow from earlier results:
1. If S is independent, then |S| ≤ dim(V ) (contrapositive of Theorem 3.4.3).
2. If S generates V (span(S) = V ), then |S| ≥ dim(V ) (contrapositive of Corollary 3.4.1).
Corollary 3.4.2. Let V be a vector space, dim(V ) = n < ∞, S ⊆ V , and |S| = dim(V ). Then
1. If span(S) = V , then S is a basis for V .
2. If S is linearly independent, then S is a basis for V .
Proof.
1. Since span(S) = V , by Theorem 3.4.2 there exists β ⊆ S so β is a basis for V . Thus,
|β| = dim(V ) = |S| and so β = S.
2. On the contrary, suppose span(S) 6= V . Then there exists v ∈ V so v 6∈ span(S). Hence
S ∪ {v} is independent. This, however is a contradiction to Theorem 3.4.3, which says that
if |S ∪ {v}| > dim(V ), then S ∪ {v} is dependent. Therefore span(S) = V and so S is a
basis.
Example 3.4.4. Do the polynomials x3 2x2 + 1, 4x2 x + 3, and 3x2 generate P3 (R)? No, since
|S| < dim(R3 ), span(S) 6= P3 (Corollary 3.4.1).
Example 3.4.5. Determine whether S = {(1, 0, −1), (2, 5, 1), (0, −4, 3)} form a basis for R3 .
Since |S| = 3 = dim(R3 ), it is sufficient to show that S is linearly independent (Corollary 3.4.2).
Equating coefficients
1 2 0 1 0 0
0 5 −4 rref ∼ 0 1 0 .
−1 1 3 0 0 1
Hence S is linearly independent and forms a basis for R3 .
Example 3.4.6. Find a basis for the following subspaces of F 5 :
W = {(a1 , a2 , a3 , a4 , a5 ) ∈ F 5 : a1 − a3 − a4 = 0}.
What is the dimensions of W ? (a1 , a2 , a3 , a4 , a5 ) ∈ W1 if and only if (a1 , a2 , a3 , a4 , a5 ) =
(s, t, r, t, 2s) for some r, s, t ∈ R. Thus, a spanning set for W1 is
β = {(0, 0, 1, 0, 0), (1, 0, 0, 0, 2), (0, 1, 0, 1, 0)}.
Since
0 1 0 1 0 0
0 0 1 0 1 0
rref
1 0 0
∼ 0 0 1
,
0 0 1 0 0 0
0 2 0 0 0 0
the spanning set, β, is linearly independent and thus forms a basis for W . The dim(W ) = 4.
56
A. W. Richter 3.5. LINEAR TRANSFORMATIONS
Definition 3.5.1 (Linear Transformation). Let V and W be vector spaces over F . We call a function
T : V → W a Linear Transformation from V into W if for all x, y ∈ V and c ∈ F we have
Theorem 3.5.1. The following are basic facts about the function T : V → W :
Now, assume T (cx + y) = cT (x) + T (y) for all x, y ∈ V and c ∈ F . Let x, y ∈ V . Then we
obtain T (x + y) = T (1x + y) = 1 · T (x) + T (y) = T (x) + T (y). Next, let x ∈ V and c ∈ F .
Then we obtain T (cx) = T (cx + 0) = c · T (x) + T (0) = c · T (x). This proves T is linear. The
same type of reasoning can be used to show that if T satisfies
n n
!
X X
T ai xi = ai T (xi ),
i=1 i=1
Example 3.5.1. Given A ∈ Mm×n (F ), define LA : Mn×1 → Mm×1 (F ) by LA (x) = Ax. To see
that this is a linear transformation, note that for any x, y ∈ Mn×1 and any c ∈ R
57
A. W. Richter 3.5. LINEAR TRANSFORMATIONS
T
V W
0V R(T)
N(T) 0W
The next theorem provides a method for finding a spanning set for the range of a linear trans-
formation. With this accomplished, a basis for the range is easy to discover.
58
A. W. Richter 3.5. LINEAR TRANSFORMATIONS
Reflecting on the action of a linear transformation, we see intuitively that the larger the nullity,
the smaller the rank. In other words, the more vectors that are carried into 0, the smaller the range.
The same heuristic reasoning tells us that the larger the rank, the smaller the nullity. The balance
between rank and nullity is made precise in the next theorem.
Theorem 3.5.3 (Dimension Theorem). Let V and W be vector spaces, and let T : V → W be
linear. If V is finite-dimensional, then
nullity(T ) + rank(T ) = dim(V ).
Proof. Suppose that dim(V ) = n, dim(N (T )) = k, and {v1 , . . . , vk } is a basis for N (T ). We may
extend {v1 , . . . , vk } to a basis β = {v1 , . . . , vn } for V . We claim that S = {T (vk+1 ), . . . , T (vn )}
is a basis for R(T ).
First we prove that S generates R(T ). Using Theorem 3.5.2 and the fact that T (vi ) = 0 for
1 ≤ i ≤ k, we have
R(T ) = span({T (vk+1 ), . . . , T (vn )})
Now we prove that S is linearly independent. Form
n
X
bi T (vi ) = 0 for bk+1 , . . . , bn ∈ F.
i=k+1
Since β is a basis for V , we have bi = 0 for all i. Hence S is linearly independent. Notice that this
argument also shows that T (vk+1 ), . . . , T (vn ) are distinct, and hence rank(T ) = n − k.
59
A. W. Richter 3.5. LINEAR TRANSFORMATIONS
Theorem 3.5.4. Let V and W be vector spaces, and let T : V → W be linear. Then T is one-for-
one if and only if N (T ) = {0}.
Proof. First note that T is one-to-one if and only if T (x) = T (y) implies x = y. Suppose N (T ) =
{0} and T (x) = T (y). Since T is linear, T (x − y) = T (x) − T (y) = 0. Thus, x − y ∈ N (T ).
By assumption, N (T ) = {0}. Hence, x − y = 0, which implies x = y. Now assume that T is
injective. Let x ∈ N (T ), then T (x) = 0 = T (0). Hence x = 0, since T is injective.
To find the null space, we must find an x ∈ M2×3 (R) such that
Equating coefficients
2 −1 0 rref 1 0 1/4
∼ .
0 2 1 0 1 1/2
Thus, N (T ) = {(−r/4, −r/2, r, s, t, u)} for r, s, t, u ∈ R. Hence, a basis for the null space can be
written
−1/4 −1/2 1 0 0 0 0 0 0 0 0 0
βN = , , ,
0 0 0 1 0 0 0 1 0 0 0 1
and rank(T ) = dim(R(T )) = 2. This verifies the dimension theorem, since dim(M2×3 (R)) =
6 = 2 + 4. Since nullity(T ) 6= 0, T is not one-to-one. Since rank(T ) = 2 6= dim(M2×3 (R)), T is
not onto.
60
A. W. Richter 3.5. LINEAR TRANSFORMATIONS
and define LA : M5×1 → M3×1 by LA (x) = Ax. To find a basis for N (LA ), solve Ax = 0. We
have
1 −1 −1 2 1 1 −1 0 1 2
rref
A = 2 −2 −1 3 3 ∼ 0 0 1 −1 1 .
−1 1 −1 0 −3 0 0 0 0 0
Thus, N (T ) = (r − s − 2t, r, s − t, s, t) for r, s, t ∈ R. Hence, a basis for the null space can be
written
1 −1 −2
1 0 0
βN = 0 , 1 , −1
0 1 0
0 0 1
61
Chapter 4
Proof. Let X and Y be any two convex subsets of Rn . Suppose αx′ and αx′′ are any two points in
αX, naturally with x′ and x′′ in X. Given λ ∈ [0, 1], since X is convex, we know λx′ +(1−λ)x′′ ∈
X. Thus,
λ(αx′ ) + (1 − λ)(αx′′ ) = α[λx′ + (1 − λ)x′′ ] ∈ αX.
Therefore αX is shown to be convex. Now suppose x′ + y′ and x′′ + y′′ are any two points in
X + Y , naturally with x′ and x′′ in X and y′ and y′′ in Y . Given λ ∈ [0, 1], since X and Y are
both convex, we know λx′ + (1 − λ)x′′ ∈ X and λy′ + (1 − λ)y′′ ∈ Y . Thus,
λ(x′ + y′ ) + (1 − λ)(x′′ + y′′ ) = λx′ + (1 − λ)x′′ + λy′ + (1 − λ)y′′ ∈ X + Y
Therefore X + Y is shown to be convex.
62
A. W. Richter 4.2. CONCAVE AND CONVEX FUNCTIONS
(a) Strictly Convex (b) Convex, but not Strictly Convex (c) Not Convex
and is strictly concave if the inequality holds strictly for λ ∈ (0, 1), that is, if
∀x′ , x′′ ∈ X and λ ∈ (0, 1), (1 − λ)f (x′ ) + λf (x′′ ) < f [(1 − λ)x′ + λx′′ ] ≡ f (xλ ).
Remark 4.2.1. Reversing the direction of the inequalities in the theorem, we obtain the definitions
of convexity and strict convexity.
Many introductory calculus texts call convex functions “concave up” and concave functions
“concave down”, as we did in section 1.1.4. Henceforth, we will stick with the more classical
terms: “convex” and “concave”.
Moreover, f is strictly concave if and only if the inequality holds strictly, that is, if and only if
Remark 4.2.2. This theorem says that a function f is concave if and only if the graph of f lies
everywhere on or below any tangent plane. Equivalently, it says that a function is concave if and
only if the slope of the function at some arbitrary point, say x0 < x, is greater than the slope of the
secant line between points x and x0 . Reversing the direction of the inequalities in the theorem, we
obtain a theorem corresponding to convexity and strict convexity.
63
A. W. Richter 4.2. CONCAVE AND CONVEX FUNCTIONS
Proof. Consider any x′ , x′′ ∈ X and λ ∈ [0, 1]. Since f is concave, we have
This establishes that h is concave. If f is strictly concave and g is strictly increasing, then the
inequality is strict and hence h is strictly concave.
Example 4.2.1. Let the domain be R++ . Consider h(x) = e1/x . Let f (x) = 1/x and let g(y) = ey .
Then h(x) = g[f (x)]. Function f is strictly convex and g is (strictly) convex and strictly increasing.
Therefore, by Theorem 4.2.2, h is strictly convex.
Remark 4.2.3. It is important in Theorem 4.2.2 that g is increasing. To see this, let the domain be
2
R++ .√Consider h(x) = e−x , which just the standard normal density except that it is off by a factor
2
of 1/ 2π. Let f (x) = ex and let g(y) = 1/y. Then h(x) = g[f (x)]. Now, f is convex on R++
and g is also convex on R++ . The function h is not, however, convex. While it is strictly convex for
|x| sufficiently large, for x near zero it is strictly concave. This does not contradict Theorem 4.2.2
because g here is decreasing.
64
A. W. Richter 4.3. CONCAVITY, CONVEXITY, AND DEFINITENESS
Proof. Consider any x′ , x′′ ∈ X and λ ∈ [0, 1]. If each fi is concave, we have
Therefore,
n
X n
X
f (xλ ) ≡ f (λx′ + (1 − λ)x′′ ) = αi fi (λx′ + (1 − λ)x′′ ) ≥ αi [λfi (x′ ) + (1 − λ)fi (x′′ )]
i=1 i=1
n
X Xn
=λ αi fi (x′ ) + (1 − λ) αi fi (x′′ ) ≡ λf (x′ ) + (1 − λ)f (x′′ ).
i=1 i=1
This establishes that f is concave. If some fj is strictly concave and αj > 0, then the inequality is
strict.
Remark 4.3.1. Note that if f is strictly concave, the Hessian can either be negative semi-definite
or negative definite. Thus, if you show that the Hessian is not negative definite, but only negative
semi-definite, you cannot conclude that f is not strictly concave.
12x2 + 2y 2
4xy
H(f ) = .
4xy 2x2 + 12y 2
(1) (2)
The principle minors, B1 = 12x2 + 2y 2 , B1 = 2x2 + 12y 2 , and B2 = 24x4 + 132x2 y 2 + 24y 4
are all weakly positive for all values of x and y, so f is a convex function on all Rn .
65
A. W. Richter 4.4. QUASI-CONCAVE AND QUASI-CONVEX FUNCTIONS
Example 4.3.2. A simple utility or production function is f (x, y) = xy. Its Hessian is
0 1
H(f ) = ,
1 0
whose second order principle minor is det H(f ) = −1. Since this second order principle minor is
negative, H(f ) is indefinite and f is neither concave nor convex.
Example 4.3.3. Consider the monotonic transformation of the function f in the previous example
by the function g(z) = z 1/4 : g[f (x, y)] = x1/4 y 1/4 , defined only on the positive quadrant R2+ . The
hessian of g is 3 −7/4 1/4 1 −3/4 −3/4
− x y x y
H(g) = 1 16−3/4 −3/4 16 3 1/4 −7/4 .
16 x y − 16 x y
The first order principle minors are both non-positive and the second order principle minor, x−3/2 y −3/2 /32,
is non-negative. Therefore, H(g) is negative semi-definite on R2+ and G is a concave function on
R2+ .
We say that f is strictly quasi-concave (quasi-convex) if for all x′ and x′′ in X and all λ ∈ (0, 1)
we have
Theorem 4.4.1. Let f be a real-valued function defined on a convex set X ⊆ Rn . Then f is quasi-
concave (quasi-convex) if and only if the upper contour sets (lower contour sets) of f are all convex,
that is, if for any a ∈ R the set
Ua = {x ∈ X : f (x) ≥ a}
(La = {x ∈ X : f (x) ≤ a})
is convex.
Proof. Assume f is quasi-concave. Fix a and let x′ , x′′ ∈ Ua . Then for all λ ∈ [0, 1],
Since x′ , x′′ ∈ Ua , f (x′ ) ≥ a and f (x′′ ) ≥ a, which implies the min{f (x′ ), f (x′′ )} ≥ a. There-
fore, f (λx′ + (1 − λ)x′′ ) ≥ a and thus λx′ + (1 − λ)x′′ ∈ Ua . Hence, the upper contour set is
convex. Now assume the upper contour set is convex. Then for all λ ∈ [0, 1] and for x′ , x′′ ∈ Ua ,
we have λx′ + (1 − λ)x′′ ∈ Ua . This implies that f (λx′ + (1 − λ)x′′ ) ≥ a. Since this result
must hold for any a, it must hold for a = min{f (x′ ), f (x′′ )}. Thus, f is quasi-concave. A similar
argument holds for quasi-convexity.
66
A. W. Richter 4.4. QUASI-CONCAVE AND QUASI-CONVEX FUNCTIONS
Theorem 4.4.2. All concave (convex) functions are quasi-concave (quasi-convex) and all strictly
concave (strictly convex) functions are strictly quasi-concave (strictly quasi-convex).
Proof. Without loss of generality assume f (x′ ) ≥ f (x′′ ). Since f is concave, for 0 ≤ λ ≤ 1
Thus, f is also quasi-concave. If f is strictly concave then the inequalities become strict and hence
f is strictly quasi-concave. The convex case can be proved in a similar fashion.
It is important to note that the converse of the above theorem is not valid in general (see
figure 4.3). The function f defined on X = {x : x ≥ 0} by f (x) = x2 is quasi-concave (UCS is
convex) but not concave on X, actually it is strictly convex on X.
φ[f (λx′ + (1 − λ)x′′ )] ≥ φ[min{f (x′ ), f (x′′ )}] = min{φ[f (x′ )], φ[f (x′′ )]}.
67
A. W. Richter 4.5. QUASI-CONCAVITY, QUASI-CONVEXITY, AND DEFINITENESS
φ[f (λx′ + (1 − λ)x′′ )] ≤ φ[min{f (x′ ), f (x′′ )}] = max{φ[f (x′ )], φ[f (x′′ )]}.
Remark 4.4.1. The sum of quasi-concave functions need not be quasi-concave unlike the sum of
concave functions which is concave. For instance f1 (x) = x3 and f2 (x) = −x are both quasi-
concave, but the sum f3 (x) = f1 (x) + f2 (x) = x3 − x is neither quasi-concave nor convex.
which is concave since for all i, ln xi is concave and the sum of concave functions is concave
(Theorem 4.2.3). Thus, since concavity implies quasi-concavity, ln f (x) is also quasi-concave. The
exponent et is strictly increasing function: R → R, hence f (x) = exp(ln f (x)) is quasi-concave
by Theorem 4.4.3.
(−1)k+1 ∆k ≥ 0 ∀k = 2, . . . , n + 1 ∀x ∈ X
(∆k ≤ 0 ∀k = 2, . . . , n + 1 ∀x ∈ X)
68
A. W. Richter 4.5. QUASI-CONCAVITY, QUASI-CONVEXITY, AND DEFINITENESS
Remark 4.5.1. Be very careful with the direction of the above definitions. In figure 4.4, the function
f (x) = x2 is not quasi-concave since the upper contour set: {x ∈ R|f (x) ≥ 3} = (−∞, −3] ∪
[3, ∞) is not a convex set. However,
0 2x
H(f ) =
2x 2
so (−1)2+1 ∆2 = 4x2 ≥ 0 for all x with strict inequality everywhere except at x = 0. Although
this function fulfills the necessary condition for quasi-concavity, the function is not quasi-concave.
Example 4.5.1. Prove or give a counterexample: If f is a strictly convex function, then f cannot be
quasi-concave.
Solution: False. For f (x) = 1/x defined on R++ , the upper contour set {x : f (x) ≥ c} is
convex for all c; For example, for c = 1, the upper contour set is the interval (0,1], which is clearly
convex. Thus, this function is quasiconcave. However, the function is also strictly convex since
(f ′′ (x) = 2/x3 > 0 for all x > 0).
Example 4.5.2. Consider f (x) = x3 + x. For x ∈ R, the second order condition, fxx = 6x, is not
always nonpositive. Thus, this function is not concave. The bordered Hessian is given by
3x2 + 1
0
H(f ) =
3x2 + 1 6x
69
A. W. Richter 4.5. QUASI-CONCAVITY, QUASI-CONVEXITY, AND DEFINITENESS
Example 4.5.3.
For the region with x > 0 and −x < y < x, define f (x, y) = x2 − y 2 . Is f concave where defined?
Is f quasiconcave where defined?
Solution: Given f (x, y) = x2 − y 2 ,
(a) fx = 2x, fxx = 2, fy = −2y, fxy = 0, fyy = −2. The Hessian is
2 0
H= .
0 −2
Since
2 0
∆1 = 2 > 0 and ∆2 =
= −4 < 0,
0 −2
H is not negative semidefinite so f is not concave.
70
A. W. Richter 4.5. QUASI-CONCAVITY, QUASI-CONVEXITY, AND DEFINITENESS
For all x, y > 0, (−1)2+1 ∆2 = 4x2 y 4 > 0 and (−1)3+1 ∆3 = 16x4 y 4 > 0. Thus, for (x, y) ∈
R2++ , the sufficient conditions for quasi-concavity hold since the inequalities are both strict. It
remains to check whether the function f is quasiconcave on R2+ . Since {(x, y) ∈ R2+ |f (x, y) ≥
0} = R2+ , all upper contour sets are convex, and f is quasi-concave on R2+ . Note f is not quasi-
concave on R2 since
Example 4.5.6. Is f (x, y) = ln(x + y) quasi-concave on the set of strictly positive x and y values?
Solution: The Hessian is given by
" #
1 1
− (x+y)2 − (x+y)2
H= 1 1 .
− (x+y)2 − (x+y)2
(1) (2) 1
For all x, y > 0, B1 = B1 = − (x+y) 2 < 0 and ∆2 = |H| = 0. Thus, H is negative semidefinite
2
on R++ , implying that f is concave, and hence quasi-concave.
√ √
Example 4.5.7. Is f (x, y, z) = x + y + z 2 concave and/or quasi-concave on R3++ ?
Solution: The Hessian is given by
1 −3/2
−4x 0 0
0 − 14 y −3/2 0 ,
0 0 2
which is not negative semi-definite (∆1 < 0, ∆2 > 0, and ∆3 > 0), so the function is not concave.
The bordered Hessian 1 −1/2 1 −1/2
0 2x 2y 2z
1 x−1/2 − 1 x−3/2 0 0
2 4
1 y −1/2 0 − 1 −3/2
y 0
2 4
2z 0 0 2
has determinant (2z 2 x3 − x5/2 − x3/2 y)/8, which is positive for some (x, y, z) ∈ R3++ and negative
for others. Thus the function is not quasi-concave on R3++ .
Thus, f is not negative semidefinite since ∆2 = −4x2 < 0. Hence f is not concave. The bordered
Hessian is given by
2xy x2
0
H 2xy 2y 2x
x2 2x 0
Thus, (−1)3 ∆2 = 4x2 y 2 > 0 for all (x, y) ∈ R2++ and (−1)4 ∆3 = 6x4 y > 0 for all (x, y) ∈ R2++ .
Thus f is quasi-concave.
71
Chapter 5
Optimization
(b) (Sufficient Condition) f has a strict local maximum (minimum) at x∗ if ∇f (x∗ ) = 0 and
H(f (x∗ )) is negative (positive) definite.
(c) If H(f (x∗ )) is indefinite, then x∗ is neither a local maximum nor a local minimum.
Remark 5.1.1. Note that in the univariate case ∇f = 0 is replaced by f ′ = 0 and H(f ) NSD
(PSD, ND, and PD, respectively) is replaced by f ′′ ≤ (≥, <, >) 0.
Example 5.1.1. Let f (x, y) = −3x2 + xy − 2x + y − y 2 + 1. Then
−6x + y − 2 −6 1
∇f = , H(f ) = ,
x + 1 − 2y 1 −2
−6 1
∆1 = −6 < 0, and ∆2 = = 11 > 0,
1 −2
72
A. W. Richter 5.1. UNCONSTRAINED OPTIMIZATION
Here H(f ) depends on the (x, y) at which it is evaluated. ∇f (x, y) = 0 if 3x2 + 2xy = 0 and
x2 + 4y = 0. From the last equation y = −x2 /4. Substituting this value into the previous equation,
(x2 /2)(6 − x) = 0. Thus, (x, y) = (0, 0) or (x, y) = (6, −9).
0 0
H[f (0, 0)] = ,
0 4
(1) (2)
which is PSD (B1 = 0 = 0, B1 = 4, and B2 = 0).
18 12
H[f (6, −9)] =
12 4
Thus, since the Hessian is indefinite at (6, −9), this point is a saddle point for the function. The
point (0, 0) satisfies the necessary conditions for a local minimum, but not the sufficient conditions.
However, f (x, 0) = x3 so f cannot attain either a local maximum or a local minimum at (0, 0).
This function has no local maxima or local minima.
Remark 5.1.2. It is interesting to compare Theorems 5.1.1 and 5.1.2. In order to guarantee that a
critical point x∗ of a C 2 function f is a strict local maximum (minimum), we need to show that
H(f (x∗ )) is negative (positive) definite; showing that H(f (x∗ )) is negative (positive) semi-definite
is not strong enough. However, if we can show that H(f (y)) is negative (positive) semi-definite not
just at x∗ but for all y in a neighborhood about x∗ , then by Theorem 5.1.2, we can conclude that x∗
is a maximum (minimum) of f .
Remark 5.1.3. A global maximum (or minimum) does not necessarily have to be a strict maximum
(or minimum). Moreover a strict maximum (or minimum) does not necessarily have to be a global
maximum (or minimum). To see this consider figure 5.1.
Example 5.1.3. Find all local maxima and minima of f (x, y, z) = x2 + x(z − 2) + 3(y − 1)2 + z 2
Solution: The associated first order conditions are
∂f set
= 2x + z − 2 = 0
∂x
∂f set
= 6(y − 1) = 0
∂y
∂f set
= x + 2z = 0
∂z
73
A. W. Richter 5.2. CONSTRAINED OPTIMIZATION I: EQUALITY CONSTRAINTS
(a) Global not Strict (b) Global and Strict (c) Strict not Global
∆1 = 2 > 0
2 0
∆2 =
= 12 > 0
0 6
∆3 = |H| = 18 > 0
Thus, H is PD, and therefore f is strictly convex. From the previous theorem, we can conclude that
the point (x, y, z) = ( 43 , 1, − 23 ) is the unique global minimizer (and a strict local minimizer).
where g(x) = 0 denotes an m × 1 vector of constraints, m < n. The condition m < n is needed
to ensure a proper degree of freedom. Without this condition, there would not be a way for the
variables to adjust toward the optimum.
74
A. W. Richter 5.2. CONSTRAINED OPTIMIZATION I: EQUALITY CONSTRAINTS
The most intuitive solution method for problem (5.1) involves the elimination of m variables
from the problem by use of the constraint equations, thereby converting the problem into an equiv-
alent unconstrained optimization problem. The actual solution of the constraint equations for m
variables in terms of the remaining n − m can often prove a difficult, if not impossible, task. More-
over, the elimination of variables is seldom applicable to economic problems, as economic theory
rarely allows for the specification of particular functional forms. Nevertheless, the theory underly-
ing the method of elimination of variables can be used to obtain analytically useful characterizations
of solutions to equality constrained problems.
Alternatively, the solution can be obtained using the Lagrangian function defined as
where λ1 , λ2 . . . , λm multiply the constraints and are known as Lagrange multipliers. In order to
solve the problem, we find the critical points of the Lagrangian by solving the equations
∂L ∂L ∂L
= 0; = 0; · · · ; = 0;
∂x1 ∂x2 ∂xn
(5.3)
∂L ∂L ∂L
= 0; = 0; · · · ; = 0,
∂λ1 ∂λ2 ∂λm
which represent n + m equations for the n + m variables x1 , x2 , . . . , xn , λ1 , λ2 . . . , λm . Thus, we
have transformed what was a constrained problem of n variables into an unconstrained problem
∂L
of n + m variables. Note that since λ1 , λ2 . . . , λm simply multiply the constraints, ∂λ i
, for i =
1, 2, . . . , m, is equivalent to each multipliers’ respective constraint. Thus, the system of equations,
(5.3), can be written more compactly as
which are perpendicular to the level sets of f and g. Since the level sets of f and g have the same
slope at x∗ , the gradient vectors ∇f (x) and ∇g(x) must line up at x∗ . Thus they point in the same
direction or opposite directions (see figure 5.2). In either case, the gradients are scalar multiples
of each other. If the corresponding Lagrange multiplier is λ∗ , then ∇f (x∗ ) = λ∗ ∇g(x∗ ) as the
Lagrange formulation, given in (5.4), suggests.
75
A. W. Richter 5.2. CONSTRAINED OPTIMIZATION I: EQUALITY CONSTRAINTS
x′
x′′
∇f (x∗ ) ∇f (x∗ )
x∗ ∇g(x∗ ) x∗
∗
x
∇g(x∗ )
C C C
Remark 5.2.1. In order for this transformation to remain valid, we must place a restriction on the
constraint set known as constraint qualification, which requires:
(a) ∇g(x∗ ) 6= 0 if the problem defined in (5.1) has only one constraint; and
equals m (full rank) if the problem defined in (5.1) has m constraints, m > 1.
Example 5.2.1. This example illustrates why the rank m condition is required for the transformation
given in (5.2). Suppose our problem is to
maximize f (x1 , x2 , x3 ) = x1
subject to g1 (x1 , x2 , x3 ) = (x1 − 1)2 − x3 = −1
g2 (x1 , x2 , x3 ) = (x1 − 1)2 + x3 = 1.
The set of points satisfying both constraints is {(1, y, 1)|y ∈ R}. If the transformation in (5.2) is
valid, (5.4) implies
which is not possible since the gradient vectors are linearly dependent. The problem here is that the
transformation is not valid, since constraint qualification is not satisfied (rank J(x∗ ) = 1 < m).
76
A. W. Richter 5.2. CONSTRAINED OPTIMIZATION I: EQUALITY CONSTRAINTS
Remark 5.2.2. Constraint qualification says that for transformation (5.2) to be valid, no point satis-
fying the constraint set can be a critical point of the constraint set. This means that if the constraint
set is linear, constraint qualification will automatically be satisfied.
(b) (Sufficient Condition) If there exist vectors x∗ ∈ Rn , λ∗ = (λ∗1 , λ∗2 , . . . , λ∗m ) ∈ Rm such that
∇L(x∗ , λ∗ ) = 0
∇gi (x∗ ) · z = 0, i = 1, 2, . . . , m
it follows that
Conveniently, these conditions for a maximum or minimum can be stated in terms of the Hessian
of the Lagrangian function, which turns out to be a bordered Hessian. The following rules work with
the bordered Hessian of a constrained optimization problem of the form:
..
L · · · Lλ1 λm . Lλ1 x1 · · · Lλ1 xn
λ.1 λ1 .. .. .. ..
. .. ..
. . . . . . .
.. ..
0 . B T
L · · · L . L · · · L
λm λ1 λm λm λm x1 λm xn
H=. . . . . . . . . . = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. ..
B .. A L
x1 λ1 · · · L x1 λm . L x1 x1 · · · L x1 xn
. .. .. .. .. .. ..
.. . . . . . .
..
Lxnλ1 · · · Lxn λm . Lxn x1 · · · Lxn xn
..
0 ··· 0 . ∂g1 /∂x1 · · · ∂g1 /∂xn
.. .. .. .. .. .. ..
. . . . . . .
..
0 ··· 0 . ∂gm /∂x1 · · · ∂gm /∂xn
= . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
∂g /∂x · · · ∂g /∂x ...
1 1 m 1 L x1 x1 · · · L x1 xn
.. .. .. .. .. .. ..
. . . . . . .
..
∂g /∂x · · · ∂g /∂x .
1 n m n L xn x1 ··· Lxn xn
77
A. W. Richter 5.2. CONSTRAINED OPTIMIZATION I: EQUALITY CONSTRAINTS
Criterion 5.2.1 (Sufficient Conditions for strict local Maximum with constraints). Let f and con-
straints g1 , g2 , . . . , gm be twice continuously differentiable real-valued functions. If there exist vec-
tors x∗ ∈ Rn , λ∗ = (λ∗1 , λ∗2 , . . . , λ∗m ) ∈ Rm such that
∇L(x∗ , λ∗ ) = 0
and if H is negative definite on the constraint set, which is the case if (−1)k ∆m+k > 0 for k =
m + 1, ..., n, where m is the number of constraints that hold with equality, n is the number of
endogenous variables, and ∆k is the leading principle minor of order k (Note: Lagrange multipliers
do not count as endogenous variables), then f has a strict local maximum at x∗ .
Criterion 5.2.2 (Sufficient Conditions for strict local Minimum with constraints). Let f and con-
straints g1 , g2 , . . . , gm be twice continuously differentiable real-valued functions. If there exist vec-
tors x∗ ∈ Rn , λ∗ = (λ∗1 , λ∗2 , . . . , λ∗m ) ∈ Rm such that
∇L(x∗ , λ∗ ) = 0
and if H is positive definite on the constraint set, which is the case if (−1)m ∆m+k > 0 for k =
m + 1, ..., n, then f has a strict local minimum at x∗ .
Remark 5.2.3. In short, we must check n − m leading principle minors starting with the principle
minor of highest order and working backwards. For example, if a problem contains 5 variables and
3 constraints, it will be necessary to check the signs of two principle minors: ∆7 and ∆8 .
Example 5.2.2 (Minimizing Cost subject to an Output Constraint). Consider a production function
given by
Let the prices of x1 and x2 be 10 and 5 respectively with an output constraint of 55. Then to mini-
mize the cost of producing 55 units of output given these prices, we set up the following Lagrangian
Now plug this into the last first order condition to obtain
78
A. W. Richter 5.2. CONSTRAINED OPTIMIZATION I: EQUALITY CONSTRAINTS
Therefore, we have two potential solutions (x1 , x2 ) = (19, 12) and (x1 , x2 ) = (1, 3). The Lagrange
multiplier λ is obtained by plugging the solutions into the above first order conditions to obtain
5
10 − λ(20 − 2(19)) = 0 → λ=−
9
5
10 − λ(20 − 2(1)) = 0 → λ= .
9
Given that ∇y(19, 12) = (−18, −9) 6= 0 and ∇y(1, 3) = (18, 9) 6= 0, constraint qualification
holds. To check for a maximum or minimum, we set up the bordered Hessian. Consider first the
point (19, 12, −5/9). The bordered Hessian in this case is
∂ 2 L(x∗ ,λ∗ ) ∂ 2 L(x∗ ,λ∗ ) ∂g(x∗ )
∂x21
∂ 2 L(x ∂x1 ∂x2 ∂x1 2λ 0 20 − 2x1
∗ ,λ∗ ) ∂ 2 L(x∗ ,λ∗ ) ∂g(x∗ )
H = ∂x2 ∂x1 ∂x2
= 0 2λ 15 − 2x2
∂x22
∂g(x )
∗ ∂g(x )∗ 20 − 2x1 15 − 2x2 0
∂x1 ∂x2 0
10
−9 0 −18
= 0 − 109 −9 .
−18 −9 0
Since we have only two endogenous variables (n = 2) and one constraint (m = 1), it is sufficient
to check only the principle minor of highest magnitude (∆3 or, more precisely, the determinant of
the bordered Hessian) in order to determine definiteness.
10 − 10 − 10
2 9 −9 4
0
9
|H| = (−1) − + (−1) (−18)
9 −9 0 −18 −9
10
= − (−81) + (−18)(−20)
9
= 450.
Since k = 2, (−1)2 ∆3 = (−1)2 (450) > 0. Therefore, H is negative definite on the constraint set,
and thus this point is a strict local maximum.
Now consider the other point, (1, 3, 5/9). The bordered Hessian is given by
10
2λ 0 20 − 2x1 9 0 18
0 2λ 15 − 2x2 = 0 109 9 .
20 − 2x1 15 − 2x2 0 18 9 0
Again, it is sufficient to check only the determinant of the bordered Hessian in order to determine
definiteness. In this case, det H = −450 and (−1)∆3 = (−1)(−450) > 0. Therefore, H is
positive definite on the constraint set, and thus this point is a strict local minimum. The minimum
cost is obtained by substituting this point into the cost expression (objective function) to obtain
79
A. W. Richter 5.3. CONSTRAINED OPTIMIZATION II: NON-NEGATIVE VARIABLES
Example 5.2.3. Consider the problem of maximizing x2 y 2 z 2 subject to the constraint g(x, y, z) =
x2 + y 2 + z 2 = 3. The Lagrangian function is given by
L = x2 y 2 z 2 + λ(3 − x2 − y 2 − z 2 )
Since n = 3 and m = 1, we have to check the signs of the two leading principle minors of
highest order, ∆3 and ∆4 . After computation, we find ∆3 = 32 and ∆4 = −192. For k = 2,
(−1)2 ∆3 = (−1)2 (32) > 0 and for k = 3, (−1)3 ∆4 = (−1)3 (−192) > 0. Therefore, H
is negative definite on the constraint set, and thus this point is a strict local maximum. By the
properties of determinant, the remaining seven critical points also satisfy the sufficiency condition
conditions and are classified as local maxima. This is an example of a situation where the solution
is globally optimal, but not unique.
If the optimum, x∗ , happens to be that these requirements are not binding, that is, x1 , x2 , . . . , xj
are in fact strictly positive, then the procedure outlined in the preceding section for determining
optimal points remains unaltered. That is, assigning a Lagrange multiplier λi , i = 1, 2, . . . , m to
each constraint, the Lagrangian function can once again be written as
80
A. W. Richter 5.3. CONSTRAINED OPTIMIZATION II: NON-NEGATIVE VARIABLES
and the first order necessary conditions satisfied at the optimum x∗ are
∂L
= 0, j = 1, 2, . . . , n (5.5)
∂xj
∂L
= 0, k = 1, 2, . . . , m, (5.6)
∂λk
so long as constraint qualification is satisfied. However, it is often the case that some components
are positive while others are zero. Thus, an equation like (5.5) should hold for the partial derivative
of the Lagrangian with respect to every component that is strictly positive and an inequality with
respect to every component that is zero. In other words for x1 , x2 , . . . , xn , we should have
∂L
≤ 0, xj ≥ 0, (5.7)
∂xj
with at least one of these holding with equality. The requirement that at least one inequality in (5.7)
should hold as an equality is sometimes stated more compactly as
∂L
xj = 0.
∂xj
The point is that the product is zero only if at least one of the factors is zero. A pair of inequalities
like (5.7), not both of which can be strict, is said to show Complementary Slackness, which we will
denote “CS”. A single inequality, say xj ≥ 0, is binding if it holds as an equality, that is, if xj is at
the extreme limit of its permitted range; the inequality is said to be slack if xj is positive, meaning it
has some room to maneuver before hitting its extreme. Each one of the pair of inequalities in (5.7)
therefore complements the slackness of the other; if one is slack the other is binding.
The intuition is as follows: if x∗j > 0, the constraint is not binding and it is possible to adjust xj
until the marginal benefit of further adjustments is zero (∂L/∂xj = 0), given the other constraints.
If, on the other hand, x∗j = 0, then the constraint is binding and the marginal benefit of increasing
xj is negative (∂L/∂xj ≤ 0). In this case, it is possible that the objective value could be improved
if negative xj values were permitted.
For all (x, y), the Hessian matrix is negative definite, since ∆1 = −4 < 0 and ∆2 = 8 > 0. Hence
f is a concave function, and any (x, y) that satisfies complementary slackness is a constrained global
maximum. It remains to locate such a point.
Our first pass at the problem is to look for a solution with x > 0, y > 0 so that ∇f (x, y) = 0.
Equating the gradient to the zero-vector gives
x − y = −2 and 2x − 3y = −5.
81
A. W. Richter 5.4. CONSTRAINED OPTIMIZATION III: INEQUALITY CONSTRAINTS
These equations are satisfied only when x = −1 and y = 1, which obviously violates the fact that
x > 0. Thus our problem is not yet solved. However, the calculation just made was not a wasted
effort, for we have in fact found the unconstrained maximum. And, since this has x < 0, it is likely
that x = 0 at the constrained maximum. We therefore look for a solution with x = 0 and y > 0, so
that ∂f /∂y = 0. Equating x and ∂f /∂y to zero we see that 10 − 6y = 0, so y = 5/3. Thus the
point (0, 5/3) satisfies the conditions
x = 0, y > 0, ∂f /∂y = 0,
and it remains to show that this point satisfies the remaining condition for a constrained maximum,
namely ∂f /∂x < 0. At (0, 5/3),
so the condition is satisfied. Thus the constrained maximum is attained where x = 0 and y = 5/3;
the constrained maximum value of f is therefore 28/3. This is of course less than the value taken
by f at the unconstrained maximum (−1, 1), which is in fact 10.
where x and y are the quantities of two different goods. We wish to choose x and y to maximize
Π(x, y) subject to the constraints x ≥ 0 and y ≥ 0.
1 1
In this case, ∂Π/∂x = 10 (24 − 6x − 3y) and ∂Π/∂y = 10 (78 − 3x − 8y). The Hessian Matrix
− 35 3
− 10
H= 3
− 10 − 45
has principle minors ∆1 < 0 and ∆2 > 0. Therefore Π(x, y) is negative definite. Hence the profit
function is concave, and the first order conditions give a global constrained maximum. It is not
hard to see that the only values of x and y for which ∂Π/∂x = ∂Π/∂y = 0 are x = −14/13 and
y = 132/13, which clearly violates the constraints. We therefore look for a solution (x, y) such that
x = 0, y > 0, ∂Π/∂x ≤ 0, and ∂Π/∂y = 0. The first, second, and fourth conditions are satisfied
where x = 0 and y = 9.75. Since 9.75 > 8, the third condition is also satisfied. Hence the solution
is x = 0, y = 9.75 and the maximal profit is 30.025.
Unfortunately, the method for finding the constrained maxima in problems with inequality con-
straints is a bit more complex than the method we used for equality constraints. The first order
82
A. W. Richter 5.4. CONSTRAINED OPTIMIZATION III: INEQUALITY CONSTRAINTS
conditions involve both equalities and inequalities and their solution eentaiils the investigations of
a number of cases. To see this more clearly, consider the following problem:
maximize f (x1 , x2 , . . . , xn )
subject to g1 (x1 , x2 , . . . , xn ) ≤ 0
g2 (x1 , x2 , . . . , xn ) ≤ 0
.. (5.8)
.
gm (x1 , x2 , . . . , xn ) ≤ 0
x1 , x2 , . . . , xn ≥ 0.
As an alternative to the procedure outlined in the previous section, we could introduce n new
constraints in addition to the m original ones:
gm+1 (x) = −x1 ≤ 0, ..., gm+n (x) = −xn ≤ 0. (5.9)
Then, if we introduce Lagrange multipliers λ1 , . . . , λm that are associated with the constraints and
µ1 , . . . , µn to go with the non-negativity constraints, our Lagrangian function is of the form
m
X n
X
L(x, λ, µ) = f (x) − λj gj (x) − µi (−xi ), (5.10)
j=1 i=1
where λ = {λ1 , . . . , λm } and µ = {µ1 , . . . , µn }. The necessary conditions for x∗ to solve this
problem are
m
∂f (x∗ ) X ∂gj (x∗ )
− λj + µi = 0, i = 1, . . . , n (i)
∂xi ∂xi
j=1
With the possibility of inequality constraints, there are now two kinds of possible solutions: one
where the constrained optimum lies in the region where gj (x) < 0, in which case constraint j
is slack, and one where the constrained optimum lies on the boundary gj (x) = 0, in which case
constraint j is binding. In the former case, the function gj (x) plays no role. x∗ still corresponds to
the optimum of the Lagrangian given in (5.10), but this time with λj = 0. The latter case, where the
optimum lies on the boundary of each constraint, is analogous to the equality constraint discussed
previously and corresponds to the optimum of the Lagrangian with λj 6= 0 for all j. In this case,
however, the sign of the Lagrange multiplier is crucial, because the objective function f (x) will
only be at a maximum if its gradient is oriented away from the region g(x) < 0 (i.e. ∇g(x) and
∇f (x) point in the same direction). We therefore have ∇f (x) = λ∇g(x) for λj ≥ 0 for all j (if
the constraint was written as g(x) ≥ 0, the gradient vectors would point in opposite directions and
∇f (x) = −λ∇g(x) for λj ≥ 0 for all j).
83
A. W. Richter 5.4. CONSTRAINED OPTIMIZATION III: INEQUALITY CONSTRAINTS
Remark 5.4.1 (Constraint Qualification). In order for the transformation in (5.10) to be valid, the
gradient vectors ∇gj (x∗ ) (j = 1, . . . , m) corresponding to those constraints that are binding at x∗
must be linearly independent. In other words, the corresponding Jacobian matrix must be full rank.
Remark 5.4.2. When solving optimization problems subject to inequality constraints, it is helpful to
map the problem into the standard form given in (5.8). If the problem is one of minimizing f (x), the
equivalent problem of maximizing −f (x) should be solved. Also, all inequality constraints should
be written as gj (x) ≤ 0 (i.e. if the original constraint was rj (x) ≤ bj , then gj (x) = rj (x) − bj ,
while if the original was rj (x) ≥ bj , then gj (x) = bj − rj (x)).
Theorem 5.4.1 (Kuhn-Tucker Necessary Conditions). Suppose that x∗ = (x∗1 , . . . , x∗n ) solves
(5.8). Suppose further that the constraint qualification is satisfied. The there exist unique num-
bers λ∗1 , . . . , λ∗m such that
∂f (x∗ ) Pm ∂gj (x∗ )
(a) ∂xi − j=1 λj ∂xi ≤ 0 (= 0 if x∗j > 0), i = 1, . . . , n
Theorem 5.4.2 (Kuhn-Tucker Sufficient Conditions). Consider problem (5.8) and suppose that x∗
∗ ∗ and (b) in theorem (5.4.1). If the Lagrangian L = f (x∗ ) −
Pm λ1 , ∗. . . λm∗ satisfy conditions (a)
and
∗
j=1 λj gj (x ) is concave, then x is optimal.
Note that in this formulation of the necessary/sufficient conditions we use the ordinary La-
grangian, not the extended Lagrangian used earlier. The exhaustive procedure for finding a solution
using this theorem involves searching among all 2m+n patterns that are possible from the (m + n)
complementary slackness conditions. Fortunately, short-cuts are usually available.
(a) If the feasible set is compact and the objective function is continuous, then the best of the
local solutions is the global solution.
(b) If the feasible set is convex and the objective function is concave, then any point satisfying
the first-order conditions is a global maximizer. If the feasible set is convex and the objective
function is strictly concave, then any point satisfying the first-order conditions is the unique
global maximizer. (Similar conclusions hold for convex objective functions and minimizers.)
(c) If the feasible set is convex and the objective function is quasi-concave, then any point satisfy-
ing the first-order conditions [with ∇f 6= 0] is a global maximizer. If, in addition, the feasible
set is strictly convex or the objective function is strictly quasi-concave, then any point satis-
fying the first-order conditions (with ∇f 6= 0) is the unique global maximizer. [To see why
we need ∇f 6= 0, consider the problem of maximizing f (x, y) = xy subject to x ≥ 0, y ≥ 0,
and x + y ≤ 2. The feasible set is convex and f is quasi-concave on R2+ . The first-order
conditions hold at (0,0) with ∇f (0, 0) = 0, but (0,0) is clearly not even a local maximizer.]
Example 5.4.1. Now let us apply the above rules to the following maximization problem
max xy
subject to x + y ≥ −1
x+y ≤2
84
A. W. Richter 5.4. CONSTRAINED OPTIMIZATION III: INEQUALITY CONSTRAINTS
Case 1: λ = 0 = µ results in s∗1 = (0, 0, 0, 0). Thus, we no longer have any binding constraints,
that is, they drop out of the Lagrangian. Therefore, m = 0 and n = 2 and we must check
the two leading principle minors of highest magnitude. With the constraint dropping out, the
bordered Hessian becomes
0 1
H= .
1 0
We then see that for k = 2, (−1)2 ∆2 = −1 < 0, which violates the rule for negative
definiteness. Thus, this point is not a local maximum.
Case 2: λ = 0, µ > 0 results in s∗2 = (1, 1, 0, 1). Then we have one binding constraint, so that
m = 1 and n = 2. Thus, we must check only the last leading principle minor. The bordered
Hessian is
0 −1 −1
H = −1 0 1 .
−1 1 0
Since (−1)2 ∆3 = det H = 2 > 0, this point is a local maximum.
Case 3: λ > 0, µ = 0 results in s∗3 = (−1/2, −1/2, 1/2, 0). Again, we have one binding con-
straint, so that m = 1 and n = 2, and we must check only the last leading principle minor.
The bordered Hessian is
0 1 1
H = 1 0 1 .
1 1 0
Since (−1)2 ∆3 = det H = 2 > 0, this point is a local maximum.
85
A. W. Richter 5.4. CONSTRAINED OPTIMIZATION III: INEQUALITY CONSTRAINTS
For all x, y ∈ R, the feasible set is not compact. However, for x, y > 0 or x, y < 0, it is, which is
the relevant set since xy < 0 when x and y have opposite signs. Comparing the output for f (1, 1)
and f (−1/2, −1/2) and noting that the feasible set is compact when x and y have the same sign
implies that (x∗ , y ∗ ) = (1, 1) is the unique global maximum.
min x2 + 2y 2 + 3z 2 ,
subject to 3x + 2y + z ≥ 17.
The Lagrangian function is (note the negative sign, so we in fact minimize once we find the maxi-
mizer of the negative of the objective function):
Case 2 (λ > 0): From (iv), 3x + 2y + z − 17 = 0. Solving (i)-(iii) along with this equation, results
in (x, y, z, λ) = (9/2, 3/2, 1/2, 3). Now check the second order conditions. The bordered
Hessian of the Lagrangian is
Lλλ Lλx Lλy Lλz 0 3 2 1
Lxλ Lxx Lxy Lxz 3 −2 0 0
H= Lyλ Lyx Lyy Lyz 2 0 −4 0 .
=
In this case, n = 3 and m = 1, and we must check whether (−1)2 ∆3 > 0 and (−1)3 ∆4 > 0.
Since (−1)2 ∆1+2 = 44 > 0 and (−1)3 ∆1+3 = 272 > 0, we have found a local maximum.
The original objective function is strictly convex (positive definite), and the feasible set is convex
(linear), so (x∗ , y ∗ ) = (9/2, 3/2) is the unique global minimizer.
√
Example 5.4.3. Find all local maximizers for the function f (x, y) = − x + y subject to x ≥ 0
and x2 y ≥ 108. Then find all global maximizers for the problem or show none exist.
Solution: The Lagrangian function is
√
L = − x + y − λ(108 − x2 y)
86
A. W. Richter 5.4. CONSTRAINED OPTIMIZATION III: INEQUALITY CONSTRAINTS
Both λ = 0 and x = 0 are inconsistent with Ly = 0. Thus, λ > 0 and x > 0 is the only possible
case, and the unique potential solution is (6, 3, 1/216). The associated bordered Hessian is
In this case, n − m = 1, so we only have one condition to check: (−1)2 ∆3 = 108 > 0. Thus,
(x∗ , y ∗ ) = (6, 3) is a strict local maximizer and the only local maximizer. In order to assess
globality, first note that both x and y must be strictly positive to be feasible, so the feasible set
turns out to be the subset of R2++ , where g(x, y) = x2 y ≥ 108. Checking the bordered Hessian
for g, we find g is strictly quasi-concave on R2++ . Thus, the feasible set, an upper contour set for
g, is convex. The objective function is quasi-concave: for any c ≤ 0, the set of (x, y) such that
√
− x + y ≥ c is {(x, y) ∈ R2 |0 ≤ x + y ≤ c2 }, which is convex (for c > 0 the set is empty).
With a quasi-concave objective function and a convex feasible set, the unique local maximizer is the
unique global maximizer.
87
Chapter 6
Comparative Statics
In many economic problems we need to know how an optimal solution or an equilibrium solution
changes when a parameter in the problem changes. For example, how does the utility-maximizing
bundle for a competitive consumer change when a price changes, or how does a market equilibrium
price change when a tax on the good changes? These are examples of comparative statics questions.
In each case we are interested in how changes in exogenous variables (the parameters determined
outside the model) affect the endogenous variables (those determined within the model).
For the consumer choice problem, the endogenous variables are the quantities demanded (cho-
sen by the consumer), while the exogenous variables are prices (outside the control of the compet-
itive consumer). For the market example, the endogenous variable is the market equilibrium price
(determined by supply and demand in the market), while the exogenous variable is the tax rate (de-
termined outside the market in some political process). In this section, you will find two extremely
helpful tools for evaluating comparative statics.
Theorem 6.1.1 (Cramer’s Rule). Let A be a nonsingular matrix. Then the unique solution x =
(x1 , . . . , xn ) of the n × n system Ax = y is
det Bi
xi = , for i = 1, . . . , n,
det A
where Bi is the matrix A with the right-hand side y replacing the ith column of A.
88
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
f (x, y) = 0. (6.3)
89
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
∂G ∂G
(x0 , y0 ) + (x0 , y0 ) · y ′ (x0 ) = 0.
∂x ∂y
Theorem 6.2.1 (Implicit Function Theorem-One Exogenous Variable). Let G(x, y) be a C 1 func-
tion on an ε-neighborhood about (x0 , y0 ) in R2 . Suppose that G(x0 , y0 ) = c and consider the
expression
G(x, y) = c.
If (∂G/∂y)(x0 , y0 ) 6= 0, then there exists a C 1 function y = y(x) defined on an interval I about
the point x0 such that:
2xey + 1
′ f1 (x, y)
y (−1, 0) = − =− 2 y = −1.
f2 (x, y) x=−1 x e − 2 x=−1
y=0 y=0
about the point (x0 , y0 ) = (4, 3). Notice that f (4, 3) satisfies (6.5). The first order partials are
∂f ∂f
= 2x − 3y and = −3x + 3y 2 .
∂x ∂y
Since (∂f /∂y)(4, 3) = 15 6= 0, Theorem (6.2.1) tells us that (6.5) does indeed define y as a C 1
function of x around x0 = 4, y0 = 3. Furthermore,
′ 2x − 3y 1
y (4, 3) = − 2 = .
3y − 3x x=4 15
y=3
90
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
Around a given point (x∗1 , . . . , x∗k , y ∗ ), we want to vary x = (x1 , . . . , xk ) and then find a y-value
which corresponds to each such (x1 , . . . , xk ). In this case, we say that equation (6.6) defines y as a
implicit function of (x1 , . . . , xk ). Once again, given G and (x∗ , y ∗ ), we want to know whether this
functional relationship exists and, if it does, how does y change if any of the xi ’s change from x∗i .
Since we are working with a function of several variables (x1 , . . . , xk ), we will hold all but one of
the xi ’s constant and vary one exogenous variable at a time. However, this puts us right back in the
two-variable case that we have been discussing. The natural extension of Theorem (6.2.1) to this
setting is the following.
G(x∗1 , . . . , x∗k , y ∗ ) = c
∂G ∗
and (x , . . . , x∗k , y ∗ ) 6= 0.
∂y 1
F1 (y1 , . . . , ym , x1 , . . . , xn ) = c1
F2 (y1 , . . . , ym , x1 , . . . , xn ) = c2
.. .. (6.8)
. .
Fm (y1 , . . . , ym , x1 , . . . , xn ) = cm ,
where y1 , . . . , ym are endogenous and x1 , . . . , xn are exogenous. Totally differentiating (6.8) the
above system of equations about the point (y∗ , x∗ ), we obtain
where all the partial derivatives are evaluated at the point (y∗ , x∗ ). By the Implicit Function Theo-
rem, the linear system (6.9) can be solved for dy1 , . . . , dym in terms of dx1 , . . . , dxn if and only if
91
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
is nonsingular at (y∗ , x∗ ). Since this system is linear, when the coefficient matrix (6.10) is nonsin-
gular, we can use the inverse of (6.10) to solve the system (6.9) for the dyi ’s in terms of the dxj ’s
and everything else. Thus, in matrix notation we obtain
∂F1 ∂F1 ∂F1 ∂F1
∂y1 · · · ∂y m
dy1 ∂x1 · · · ∂x n
dx1
.. .. .. .. = − .. .. .. .. ,
. . . . . . . .
∂Fm ∂Fm dym ∂Fm ∂Fm dxn
∂y1 ··· ∂ym ∂x1 ··· ∂xn
which implies
∂F1 ∂F1 −1 Pn ∂F1
dy1 ∂y1 ··· ∂ym i=1 ∂xi dxi
.. .. .. .. ..
. = − . . (6.11)
. . .
∂Fm ∂Fm Pn ∂Fm
dym ∂y1 ··· ∂ym i=1 ∂xi dxi
Since the linear approximation (6.9) of the original system (6.8) is a implicit function of the
dyi ’s in terms of the dxj ’s, the nonlinear system (6.8) defines the yi ’s as implicit functions of the
xj ’s in a neighborhood of (y∗ , x∗ ). Furthermore, we can use the linear solution of the dyi ’s in terms
of the dxj ’s, (6.11), to find the derivatives of the yi ’s with respect to the xj ’s at (x∗ , y∗ ). To compute
∂yk /∂xh for some fixed indices h and k, recall that this derivative estimates the effect on yk of a
one unit increase in xh (dxh = 1). So, we set all the dxj ’s equal to zero in (6.9) or (6.11) except
dxh , and then we solve (6.9) or (6.11) for the corresponding dyi ’s. Thus, (6.11)) reduces to
dy1 ∂F1 ∂F1 −1 ∂F1
dxh ∂y1 · · · ∂y m ∂xh
.. .. . . .
. ..
. = − . . . . (6.12)
dym ∂Fm ∂Fm ∂Fm
dx ∂y1 · · · ∂ym ∂xh
h
92
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
y1 = f1 (x1 , . . . , xn )
.. ..
. .
ym = fm (x1 , . . . , xn )
Furthermore, one can compute (∂fk /∂xh )(y∗ , x∗ ) = (∂yk /∂xh )(y∗ , x∗ ) by setting dxh = 1 and
dxj = 0 for j 6= h in (6.9) and solving the resulting system for dyk . This can be accomplished:
(a) by inverting the nonsingular matrix (6.10) to obtain the solution (6.12) or
F1 (x, y, a) ≡ x2 + axy + y 2 − 1 = 0
(6.15)
F2 (x, y, a) ≡ x2 + y 2 − a2 + 3 = 0
93
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
Thus, we can solve system (6.15) for x and y as functions of a near (0, 1, 2). Furthermore, using
Cramer’s rule, we obtain
∂F1 ∂F1
!
∂x ∂a 2x + ay xy
det ∂F2 ∂F2 det
dy ∂x ∂a 2x −2a
=− =− .
da ∂F1 ∂F1 2x + ay ax + 2y
∂x ∂y det
det ∂F2 ∂F2 2x 2y
∂x ∂y
Evaluating at x = 0, y = 1, a = 2, gives
2 0
det
dy 0 −4 8
(0, 1, 2) = − = = 2 > 0.
da 2 2 4
det
0 2
Therefore, if a increases to 2.1, y will increase to 1.2. Let us now use the method of total differen-
tiation to compute the effect on x. Total differentiating the non-linear system (6.15), we obtain
Evaluating at x = 0, y = 1, a = 2:
2 dx + 2 dy = 0 da
0 dx + 2 dy = 4 da.
Clearly, dy = 2 da (as we just computed above) and dx = −dy = −2da. Thus, if a increases to
2.1, x will decrease to −.2.
Y =C +I +G
C = C(Y − T )
(6.16)
I = I(r)
M s = M (Y, r),
where the nonlinear functions x 7→ C(x), r 7→ I(r), and (Y, r) 7→ M (Y, r) satisfy
∂M ∂M
0 < C ′ (x) < 1, I ′ (r) < 0, > 0, < 0. (6.17)
∂Y ∂r
System (6.16) can be reduced to
Y − C(Y − T ) − I(r) = G
M (Y, r) = M s ,
where we have defined Y and r as implicit functions of G, M s , and T . Suppose that the current
(G, M s , T ) is (G∗ , M s∗ , T ∗ ) and that the corresponding (Y, r)-equilibrium is (Y ∗ , r ∗ ). If we vary
94
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
(G, M s , T ) a little, is there a corresponding equilibrium (Y, r) and how does it change? Totally
differentiating system (6.16), we obtain
∂C ∂I ∂C
1− dY − dr = dG − dT
∂Y ∂r ∂T
∂M ∂M
dY + dr = dM s
∂Y ∂r
or, in matrix notation,
!
∂C
1− − ∂I dG − ∂C
∂Y ∂r dY ∂T dT
= (6.18)
∂M
∂Y
∂M
∂r
dr dM s
is negative by (6.17) and therefore nonzero. By Theorem 6.2.3, the system (6.16) does indeed define
Y and r as implicit functions of G, M s , and T around (Y ∗ , r ∗ , G∗ , M s∗ , T ∗ ). Inverting (6.18), we
compute !
∂M ∂I
dG − ∂C
∂Y 1 ∂r ∂r ∂T dT
= .
∂r D − ∂M 1 − ∂C dM s
∂Y ∂Y
Minimize x + 2y + 4z
(6.19)
subject to x2 + y 2 + z 2 = 21
95
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
∂L set
= −2 + 2µy = 0
∂y
∂L set
= −4 + 2µz = 0
∂z
∂L set
= −21 + x2 + y 2 + z 2 = 0.
∂µ
The first three equations can be solved to obtain x = 1/2µ, y = 1/µ, and z = 2/µ. Substituting
these results into the fourth equation yields µ2 = 1/4. Thus, there are two potential solutions:
(1, 2, 4, 1/2) and (−1, −2, −4, −1/2). The associated bordered Hessian is
Lxx Lxy Lxz Lxµ 2µ 0 0 2x
Lyx Lyy Lyz Lyµ 0 2µ 0 2y
Lzx Lzy Lzz Lzµ = 0
.
0 2µ 2z
Lµx Lµy Lµz Lµµ 2x 2y 2z 0
Minimize x + 2y + 4z
subject to x2 + y 2 + z 2 = 21 + w,
how does the optimal choice of x change as w changes from w = 0? The endogenous variables are
x, y, z, and µ. The exogenous variable is w. The first order conditions remain the same except for
the constraint (∂L/∂µ), which becomes
∂L
= −21 − w + x2 + y 2 + z 2 = 0.
∂µ
When w = 0, we already know (−1, −2, −4, −1/2) is the global maximizer (minimizer for the
original problem), since it satisfies the first order conditions. Moreover, the bordered Hessian at this
point,
−1 0 0 −2
0 −1 0 −4
H̄ ≡ 0
,
0 −1 −8
−2 −4 −8 0
remains the same and its determinant: ∆4 = −84 6= 0. Thus, we can apply Theorem 6.2.3 to obtain
dx ∂2L
dw ∂x∂w 0 −1/42
dy ∂2L
−1 0
−1/21
dw = −H̄ −1 ∂y∂w
2 L = −H̄
0 = −2/21 .
dz
∂
dw ∂z∂w
dµ ∂2L −1 1/84
dw ∂µ∂w
Although we are not interested in dµ, µ was one of the endogenous variables in the first order
conditions, so it must be included in this stage, making H̄ a 4 × 4 matrix. Also, note that almost all
terms can be read off from the problem where w = 0.
96
A. W. Richter 6.2. IMPLICIT FUNCTION THEOREM
97
Chapter 7
98
A. W. Richter 7.1. BASIC OPERATIONS
7.1.2 Moduli
p modulus or absolute value of a complex number, z, is defined as the nonnegative real number
The
x2 + y 2 and is denoted by |z|; that is,
p
|z| = x2 + y 2 .
Geometrically, the number |z| is the distance between the point (x, y) and the origin, or the length
of the vector representing z. It reduces to the absolute value in the real number system when y = 0.
Note that while the inequality z1 < z2 is meaningless unless both z1 and z2 are real, the statement
|z1 | < |z2 | means that the point z1 is closer to the origin than the point z2 .
99
A. W. Richter 7.2. EXPONENTIAL FORM
5i
(b) (1−i)(2−i)(3−i)
Solution: Multiply by the conjugate of each of the terms in the denominator
(c) (1 − i)4
Solution: Note that (1 − i)2 = −2i. Thus,
where cos(θ) = α/r and sin(θ) = β/r as illustrated in figure 7.1. In section 7.1.2, we saw that
the real number r is not allowed to be negative and is the length of the radius vector for z; that is
r = |z|. The real number θ represents the angle, measured in radians, that z makes with the positive
real axis. As in calculus, θ has an infinite number of possible values, including negative ones, that
differ by multiples of 2π. Those values can be determined from the equation tan(θ) = β/α, where
the quadrant containing the point corresponding to z must be specified. Each value of θ is called
an argument of z, and the set of all such values is denoted by arg z. The principle value of arg z,
denoted by Arg z, is the unique value Θ such that −π < Θ ≤ π. Note that
Example 7.2.1. The complex number z = −1 − i, which lies in the third quadrant, has principle
argument −3π/4. That is,
3π
Arg(−1 − i) = − .
4
It must be emphasized that, because of the restriction −π < Θ ≤ π of the principle argument Θ, it
is not true that Arg(−1 − i) = 5π/4.
100
A. W. Richter 7.2. EXPONENTIAL FORM
Imaginary
Real
Proof. In order to derive this result, first recall the Maclaurin expansions of the functions ex , cos x,
and sin x given by:
x2 x3
ex = 1 + x + + + ··· ,
2! 3!
x2 x4 x6
cos x = 1 − + − + ··· ,
2! 4! 6!
x3 x5 x7
sin x = x − + − + ··· .
3! 5! 7!
For a complex number z, define each of the functions by the above series, replacing the real variable
x with the complex variable z. We then find that
This result enables us to write the polar form (7.1) more compactly in exponential form as
z = reiθ .
101
A. W. Richter 7.3. COMPLEX EIGENVALUES
Example 7.2.2. The complex number z = −1 − i in the previous example has exponential form
√
−1 − i = 2 exp [i(−3π/4)] .
This expression is, of course, only one of an infinite number of possibilities for the exponential form
of z.
It is geometrically obvious that
z = Reiθ
is a parametric representation of the circle |z| = R centered at the origin with radius R.
Another nice aspect of complex numbers written in polar form is that their powers are easily
computed. For example
= r 2 [cos(2θ) + i sin(2θ)] ,
using the double angle formulas, which state that cos(2a) = cos2 (a) − sin2 (a) and sin(2a) =
2 sin(a) cos(a). Continuing in this manner, we obtain the following result.
Definition 7.2.1 (DeMoivre’s formula). For complex number z = α + iβ with polar representation
r [cos(θ) + i sin(θ)] and any positive integer n,
This result can alternatively be derived in a much simpler manner by noting that (reiθ )n = r n einθ .
√
Example 7.2.3. In order to put ( 3 + i)7 is rectangular form, one need only write
√ √
( 3 + i)7 = (2eiπ/6 )7 = 27 ei7π/6 = (26 eiπ )(2eiπ/6 ) = −64( 3 + i).
102
A. W. Richter 7.3. COMPLEX EIGENVALUES
a+d 1 p
λ1 = + i |(a − d)2 + 4bc|
2 2
a+d 1 p
λ2 = λ1 = − i |(a − d)2 + 4bc|.
2 2
The fact that complex eigenvalues come in complex conjugate pairs and the fact that λ 6= λ, imply
that complex eigenvalues of a 2 × 2 system are always distinct. It is only with 4 × 4 matrices that the
possibility of repeated complex eigenvalues arises. We next show that complex eigenvectors also
come in conjugate pairs.
Suppose the general form an eigenvalue of matrix A is given by λ = α + iβ. Then the corre-
sponding eigenvector, w, is a non-zero solution to
[A − (α + iβ) I] w = 0,
or reformulated
Aw = (α + iβ)w. (7.3)
Now write the complex vector w in its general form: w = u + iv, where u and v are real vectors.
Then (7.3) becomes
A (u + iv) = (α + iβ) (u + iv) . (7.4)
Applying the conjugate to both sides and recalling that for any two complex numbers z1 and z2 it
holds that z1 z2 = z 1 z 2 , we obtain
Then from (7.4) and (7.5) we see that complex eigenvectors also come in conjugate pairs. More
specifically, if u + iv is an eigenvector for α + iβ, then u − iv is an eigenvector for α − iβ. The
following theorem summarizes the discussion thus far.
p(λ) = λ2 − 2λ + 10,
103
A. W. Richter 7.3. COMPLEX EIGENVALUES
Row-reducing the coefficient matrix in the above equation implies −3iw1 + w2 = 0. Normalizing
w1 to 1, we get w2 = 3i. Thus, the first eigenvector is
1 1 0
w1 = = +i .
3i 0 3
Since we have already seen that eigenvectors or complex eigenvalues come in conjugate pairs, it
must be the case that the second eigenvector is given by
1 0
w2 = −i .
0 3
Note that we could have also found the second eigenvector using the second eigenvalue λ2 = 1− 3i.
104
Chapter 8
Lxt = xt−1
Ln xt = xt−n for n = . . . , −2, −1, 0, 1, 2, . . . .
Multiplying a variable xt by Ln thus gives the value of x shifted back n periods. Notice that if
n < 0, the effect of multiplying xt by Ln is to shift x forward in time by n periods.
Consider polynomials in the lag operator given by
∞
X
A(L) = a0 + a1 L + a2 L2 + · · · = aj Lj ,
j=0
where the aj ’s are constant and L0 ≡ 1. Operating on xt with A(L) yields a moving sum of x′ s:
105
A. W. Richter 8.2. FIRST-ORDER DIFFERENCE EQUATIONS
and if |λ| > 1, we get that the above sum is unbounded. In some instances, we will want a solution
infinitely far back in time and where the infinite sum is bounded. Thus, we sometimes require
|λ| < 1. It is useful to realize that we can use an alternative expansion for the geometric polynomial
1/(1 − λL). Formally, if the normal expansion for this infinite geometric series applies, then
(−λL)−1
1 1 1 1
= =− 1+ + + ···
1 − λL 1 − (λL)−1 λL λL (λL)2
2 3
1 −1 1 −2 1
=− L − L − L−3 − · · · , (8.2)
λ λ λ
which is particularly useful when |λ| > 1. In this case, operating on xt gives
2 ∞ j
1 1 1 X 1
xt = − xt+1 − xt+2 − · · · = − xt+j ,
1 − λL λ λ λ
j=1
which shows [1/(1 − λL)]xt to be a geometrically declining weighted sum of future values of x.
Notice that
P for this infinite sum to be finite for a constant time path xt+j = x for all j and t, the
series − ∞ j=1 (1/λ)j must be convergent, which requires that |1/λ| < 1 or, equivalently, |λ| > 1.
where xt is an exogenous (determined outside the model) variable and a is a constant. This is a first-
order difference equation, since yt is dependent on only its first lag yt−1 . Here, we are interested in
finding a solution of yt in terms of current, past, or future values of the exogenous variable xt and
(less importantly) the constant a. Put it different, we want to characterize the endogenous sequence
yt in terms of the exogenous sequence xt .
Using the lag operator defined above, we can rewrite (8.3) as follows:
Operating on both sides of this equation by (1 − λL)−1 , we can obtain a particular solution for
(8.4), denoted by ŷt , as follows:
bxt a
ŷt = + . (8.5)
1 − λL 1 − λ
Note that since a is a constant, a/(1 − λL) = a/(1 − λ) irrespective of the size of |λ| (This can
be verified by considering the expansion given in (8.1) when |λ| < 1 and the expansion given in
(8.2) when |λ| > 1). In order to obtain the general solution, we need to add a term to (8.5). For this
106
A. W. Richter 8.2. FIRST-ORDER DIFFERENCE EQUATIONS
purpose, suppose that ỹ = ŷ + wt is also a solution to (8.4). Then, using the particular solution to
(8.4), we obtain
Therefore, as long as wt = λwt−1 , ỹt is also a solution. Note that we can iterate on this condition
to obtain
where w0 ≡ c, an arbitrary initial value. Hence, the general solution to (8.3) is given by
bxt a
yt = + + λt c
1 − λL 1 − λ
∞
X a
=b λj xt−j + + λt c, (8.6)
1−λ
j=0
where c is an arbitrary constant. Notice that for yt , defined by (8.6), to be finite λj xt−j must be
small for large j. That is, we require
∞
X
lim λj xt−j = 0 for all t.
n→∞
j=n
For the case ofPxt−j = x for all j and t, the above condition requires |λ| < 1. Notice also that the
infinite sum a ∞ j
i=0 λ in (8.6) is also finite only if |λ| < 1, in which case it equals a/(1 − λ) for
a 6= 0 and 0 otherwise. Tentatively, assume that |λ| < 1.
In order to analyze (8.6), rewrite the equation for t ≥ 1 as
t−1
X ∞
X t−1
X ∞
X
j j j
yt = a λ +a λ +b λ xt−j + b λj xt−j + λt c
j=0 j=t j=0 j=t
t−1 ∞
a(1 − λt ) aλt X X
= + +b λj xt−j + bλt λj x0−j + λt c
1−λ 1−λ
j=0 j=0
t−1 ∞
a(1 − λt ) X a X
= +b λj xt−j + λt +b λj x0−j + λ0 c
1−λ 1−λ
j=0 j=0
t−1
a(1 − λt ) X
= +b λj xt−j + λt y0 (using (8.6))
1−λ
j=0
t−1
a t a X
= + λ y0 − +b λj xt−j , t ≥ 1.
1−λ 1−λ
j=0
Consider the special case in which xt = 0 for all t. Under this condition, we obtain
a t a
yt = + λ y0 − . (8.7)
1−λ 1−λ
107
A. W. Richter 8.2. FIRST-ORDER DIFFERENCE EQUATIONS
Notice that if the initial condition y0 = a/(1 − λ), then yt = y0 . In the case, y is constant across
all future time periods and a/(1 − λ) is known as a stationary point. Moreover, it is easy to see that
for |λ| < 1, the second term in (8.7) tends to zero and thus
a
lim yt = .
t→∞ 1−λ
This shows that the system is stable, tending to approach the stationary value as time passes.
The difference equation (8.4) can also be solved using the alternative representation of (1 −
λL)−1 given in (8.2). Using this this result, the general solution is given by
(−λL)−1 (−λL)−1
yt = a + bxt + λt c
1 − (λL)−1 1 − (λL)−1
∞
a X
= −b λ−j xt+j + λt c. (8.8)
1−λ
j=1
The equivalence of the solutions (8.6) and (8.8) will hold whenever
b (λL)−1
xt and bxt
1 − λL 1 − (λL)−1
are both finite. However, it is often the case that one of these two conditions fails to hold. For
example, if the sequence {xt } is bounded, this is sufficient to imply that {[b/(1 − λL)]xt } is a
bounded sequence if |λ| < 1, but not sufficient to imply that
(λL)−1
bxt
1 − (λL)−1
is a convergent sum for all t. Similarly, if |λ| > 1, boundedness of the sequence {xt } is sufficient
to imply that
(λL)−1
bxt
1 − (λL)−1
is a bounded sequence, but fails to guarantee finiteness of b/(1 − λL)xt . In instances where one of
b (λL)−1
xt or bxt
1 − λL 1 − (λL)−1
is always finite and the other is not, we shall take our solution to the first-order difference equation
(8.4) as either (8.6), where the backward sum in xt is finite, or (8.8), where the forward sum in xt is
finite. This procedure assures us that we shall find the unique solution of (8.4) that is finite for all t,
provided that such a solution exists.
If we want to guarantee that the sequence {yt } given by (8.6) or (8.8) is bounded for all t, it is
evident that we must set c = 0. This is necessary since if λ > 1 and c > 0,
lim cλt = ∞,
t→∞
lim cλt = ∞.
t→−∞
108
A. W. Richter 8.3. SECOND-ORDER DIFFERENCE EQUATIONS
Thus, when |λ| < 1, the bounded sequence yt from (8.6) can be obtained by following the backward
representation with the initial condition c = 0 and is given by
a
yt = b 1 + λL + (λL)2 + (λL)3 + · · · xt +
1−λ
∞
X a
=b λj xt−j + .
1−λ
j=0
On the other hand, when |λ| > 1, we need to use the forward representation in order to get the
bounded sequence yt as follows:
(−λL)−1 (−λL)−1
yt = bx t + a
1 − (λL)−1 1 − (λL)−1
∞
X a
= −b λ−j xt+j + ,
1−λ
j=1
again setting c = 0. In general, the convention is to solve stable roots (|λ| < 1) backward and
unstable roots (|λ| > 1) forward.
where xt is again an exogenous sequence of real numbers for t = . . . , −1, 0, 1, . . .. Using the lag
operator, we can write (8.9) as follows:
(1 − φ1 L − φ2 L2 )yt = bxt + a.
1 − φ1 L − φ2 L2 = (1 − λ1 L)(1 − λ2 L)
= 1 − (λ1 + λ2 )L + λ1 λ2 L2 ,
where λ1 λ2 = −φ2 and λ1 + λ2 = φ1 . To see how λ1 and λ2 are related to the roots or zeros of the
polynomial A(z) = 1 − φ1 z − φ2 z 2 , notice that
1 1
(1 − λ1 z)(1 − λ2 z) = λ1 λ2 −z −z .
λ1 λ2
Note that we use a function of a number z (possibly complex) instead of the lag operator L since
it does not really make much to talk about roots or zeros of a polynomial that is a function of a lag
operator. If we set the above equation to zero in order to solve for its roots, it is clear that the equation
is satisfied at the two roots z1 = 1/λ1 and z2 = 1/λ2 . Given the polynomial A(z) = 1−φ1 z−φ2 z 2 ,
the roots 1/λ1 and 1/λ2 are found by solving the characteristic equation
1 − φ1 z − φ2 z 2 = 0
109
A. W. Richter 8.3. SECOND-ORDER DIFFERENCE EQUATIONS
for two values of z. Given that λi = zi−1 for i = 1, 2, multiplying the above equation by z −2 yields
z −2 − z −1 φ1 − φ2 = λ2 − φ1 λ − φ2 = 0.
Particular Solution
If both eigenvalues are distinct as in (8.11), then the above coefficient can be written as
1 1 λ1 λ2
= − . (8.12)
(1 − λ1 L)(1 − λ2 L) λ1 − λ2 1 − λ1 L 1 − λ2 L
Thus, if either a = 0 or the magnitude of both eigenvalues is strictly less than unity (that is, if
|λ1 | < 1 and |λ2 | < 1), (8.11) can be written as
a 1 λ1 λ2
yt = + − bxt + λt1 c1 + λt2 c2
(1 − λ1 )(1 − λ2 ) λ1 − λ2 1 − λ1 L 1 − λ2 L
∞ ∞ ∞ ∞
X j
X j λ1 b X j λ2 b X j
=a λ1 λ2 + λ1 xt−j − λ2 xt−j + λt1 c1 + λt2 c2 (8.13)
λ1 − λ2 λ1 − λ2
j=0 j=0 j=0 j=0
provided that
∞
λji xt−j = 0,
X
lim for all t
n→∞
j=n
for i = 1, 2. Note that this stipulation is needed so that the corresponding geometric sums are finite.
In order to analyze (8.13), let’s first consider the special case where a = 0, so that this equation
holds regardless of the magnitude of λi . Then rewriting (8.13) for t ≥ 1 gives
t−1 t−1 ∞
λ1 b X j λ2 b X j λ1 b X j
yt = λ1 xt−j − λ2 xt−j + λ1 xt−j
λ1 − λ2 λ1 − λ2 λ1 − λ2
j=0 j=0 j=t
110
A. W. Richter 8.3. SECOND-ORDER DIFFERENCE EQUATIONS
∞
λ2 b X j
− λ2 xt−j + λt1 c1 + λt2 c2
λ1 − λ2
j=t
t−1 t−1
λ1 b X j λ2 b X j
= λ1 xt−j − λ2 xt−j + λt1 θ0 + λt2 η0 ,
λ1 − λ2 λ1 − λ2
j=0 j=0
where
∞ ∞
λ1 b λ2 b
λj1 x0−j λj2 x0−j
X X
θ0 ≡ c1 + and η0 ≡ c2 − .
λ1 − λ2 λ1 − λ2
j=0 j=0
If, without loss of generality, |λ1 | < 1 and |λ2 | > 1, then we can write
L−1
1 1 λ1
= +
(1 − λ1 L)(1 − λ2 L) λ1 − λ2 1 − λ1 L 1 − (λ2 L)−1
∞ ∞
λ1 X λ2 X
= (λ1 L)j + (λ2 L)−j .
λ1 − λ2 λ1 − λ2
j=0 j=1
111
A. W. Richter 8.3. SECOND-ORDER DIFFERENCE EQUATIONS
Moreover, using the useful polar representation defined in section 7.2, we have
Returning to the special case where a = 0 and xt = 0, when the eigenvalues are complex, (8.14),
becomes
yt = θ0 (reiw )t + η0 (re−iw )t
= θ0 (r t eiwt ) + η0 (r t e−iwt )
= θ0 r t [cos wt + i sin wt] + η0 r t [cos wt − i sin wt]
= (θ0 + η0 )r t cos wt + i(θ0 − η0 )r t sin wt.
Since yt must be a real number for all t, it follows that θ0 + η0 must be real and θ0 − η0 must be
imaginary. Therefore, θ0 and η0 must be complex conjugates, say θ0 = peiθ and η0 = pe−iθ . Thus,
we can write
where we have made use of the fact that cos is an even function and sin is an odd function (i.e.
for any input x, cos(x) = cos(−x) and sin(−x) = − sin(x)). The path of yt oscillates with a
frequency determined by w. The dampening factor, r t , is determined by the amplitude, r, of the
complex roots. When r < 1, the stationary point of the difference equation, yt = 0, is approached
as t → ∞. Moreover, as long as w 6= 0, the system displays damped oscillations. If r = 1, yt
displays repeated oscillations of unchanging amplitude and the solution is periodic. If r > 1 the
path of yt displays explosive oscillations, unless the initial conditions are say, y0 = 0 and y1 = 0 so
that y starts out at the stationary point for two successive values.
Now if we once again consider the more general case where we are interested in a bounded
sequence {yt } mapped from a bounded sequence {xt }, we need to set both of the constants c1
and c2 to zero, and focus on the associated particular solution. If we note that moduli of complex
eigenvalues are same, then when |λ| < 1, we can write
∞ ∞
1 λ1 X λ2 X
= (λ1 L)j − (λ2 L)j
(1 − λ1 L)(1 − λ2 L) λ1 − λ2 λ1 − λ2
j=0 j=0
∞
1
λj+1 − λj+1
X
= 1 2 Lj
λ1 − λ2
j=0
112
A. W. Richter 8.3. SECOND-ORDER DIFFERENCE EQUATIONS
∞
1 X iw j+1 −iw j+1
j
= iw (re ) − (re ) L
re − re−iw
j=0
∞
1 X
= 2r j+1 i sin[w(j + 1)]Lj
2ri sin w
j=0
∞
X sin[w(j + 1)] j
= rj L .
sin w
j=0
which implies that φ2 is negative. When the above condition is satisfied, the roots are given by
p p
φ1 −(φ21 + 4φ2 ) φ1 −(φ21 + 4φ2 )
λ1 = +i ≡ a + ib, λ2 = −i ≡ a − ib.
2 2 2 2
Once again, recall that in polar form
113
A. W. Richter 8.3. SECOND-ORDER DIFFERENCE EQUATIONS
√
where r ≡ a2 + b2 and tan w = β/α. Thus we have that
s
φ1 2 (φ21 + 4φ2 ) p
r= − = −φ2 .
2 4
For the oscillations to be damped,
√ meaning that in the long-run the difference equation will be
stable, we require that r = −φ2 < 1, which requires that φ2 > −1.
If the roots are real, the difference equation will be stable if both roots are less that one in
magnitude. This requires
p p
φ1 + φ21 + 4φ2 φ1 − φ21 + 4φ2
−1 < < 1 and − 1 < < 1.
2 2
Note that it is sufficient to find conditions such that statement on the left hand side is less than unity
while the condition on the right hand side is greater than minus one. The former condition requires
1
q q
2
φ1 + φ1 + 4φ2 < 1 → φ21 + 4φ2 < 2 − φ1
2
→ φ21 + 4φ2 < 4 + φ21 − 4φ1 → φ1 + φ2 < 1.
The latter condition requires
q
1
q
2
φ1 + 4φ2 − φ1 < 1 → φ21 + 4φ2 < 2 + φ1
2
→ φ21 + 4φ2 < 4 + φ21 + 4φ1 → φ2 − φ1 < 1.
Therefore, when φ2 > −1, φ1 + φ2 < 1, and φ2 − φ1 < 1 hold, the roots, regardless of whether
they are real (φ21 + 4φ2 ≥ 0) or complex (φ21 + 4φ2 < 0), will yield a stable second order difference
equation. The following figure summarizes these results.
114
A. W. Richter 8.4. SYSTEMS OF LINEAR DIFFERENCE EQUATIONS
Particular Solution
If we are interested in a bounded sequence {yt } mapped from a bounded sequence {xt }, then we
need to set both of the constants c1 and c2 to zero, and focus on the associated particular solution.
When λ1 = λ2 ≡ λ and |λ| < 1, we can show that
∞
X λL
(λL)j+1 = .
1 − λL
j=0
zt+1 = Azt ,
115
A. W. Richter 8.4. SYSTEMS OF LINEAR DIFFERENCE EQUATIONS
zt+1 = Azt .
We can go through the following steps: First multiply the above equation by P −1 and use equation
(8.17) to obtain
P −1 zt+1 = P −1 Azt
λ1 0 · · · 0
0 λ2 · · · 0 −1
= . .. P zt .
.. . .
.. . . .
0 0 ··· λk
116
A. W. Richter 8.4. SYSTEMS OF LINEAR DIFFERENCE EQUATIONS
The above form represents k separate difference equations that can be recursively solved to obtain
zt+1 = P D t+1 P −1 z0 .
which, if we multiply both sides of the above result by the projection matrix P , implies
t+1
λ1 c1
λt+1 c2
2
zt+1 = [v1 , v2 , ...vk ] .
.
.
λt+1
k ck
= v1 λt+1 t+1 t+1
1 c1 + v2 λ2 c2 + · · · + vk λk ck ,
where vk ∈ Rk represents the kth eigenvector and the constant ck ≡ Zk,0 = P −1 zk,0 . The
following theorem summarizes this alternative approach.
117
A. W. Richter 8.4. SYSTEMS OF LINEAR DIFFERENCE EQUATIONS
Theorem 8.4.1 (Difference Equation with Real Roots). Let A be a k × k matrix with k distinct
real eigenvalues λ1 , . . . , λk and corresponding eigenvectors v1 , . . . , vk . The general solution of
the system of difference equations is
xt+1 = axt
yt+1 = dyt
xt = at x0 and yt = dt y0 .
When the equations are coupled (b 6= 0 or c 6= 0), the technique for solving the system is to find
a change of variables that decouples these equations. This is precisely the role of eigenvalues and
eigenvectors.
Example 8.4.1. Consider the following coupled system of difference equations
xt+1 = xt + 4yt
yt+1 = 0.5xt
which implies that λ1 = 2 and λ2 = −1 are the the eigenvalues to this system. To find the
corresponding eigenvectors, row-reduce A − 2I and A + I to obtain the following equations
xt = 4yt
2xt = −4yt .
Normalizing yt to 1, we get the following basis vectors that form the relevant projection matrix
4 −2
P = .
1 1
118
A. W. Richter 8.4. SYSTEMS OF LINEAR DIFFERENCE EQUATIONS
Zt+1 = (P −1 AP )Zt
t
λ1 0 λ1 0
= Z = Z
0 λ2 t 0 λt2 0
t
2 0 c1
= .
0 (−1)t c2
which is the same equation that we would have arrived at had we applied Theorem 8.4.1.
zt+1 = Azt ,
where A is a 2 × 2 matrix with complex eigenvalues α ± iβ. Applying the change of variables
z = P Z to the above difference equation yields
P Zt+1 = AP Zt → Zt+1 = P −1 AP Zt .
Xt = k1 (α + iβ)t
Yt = k2 (α − iβ)t ,
where the constant K ≡ [k1 , k2 ]T = Z0 = P −1 z0 could be real or complex. Using the fact that
zt = P Zt , we can transforming transform the variables into their original form to obtain
k1 (α + iβ)t
xt
zt = = [u + iv, u − iv]
yt k2 (α − iβ)t
= k1 (α + iβ)t (u + iv) + k2 (α − iβ)t (u − iv) . (8.21)
119
A. W. Richter 8.4. SYSTEMS OF LINEAR DIFFERENCE EQUATIONS
Notice that this solution takes the same form as (8.20), except with complex eigenvalues and eigen-
vectors replacing the real eigenvectors and eigenvalues.
Since the original problem contained only real numbers, we would like to find a solution that
contains only real numbers. Since every solution of the system is contained in equation (8.21) for
different choices of K = [k1 , k2 ]T , we want to know if we can find parameters k1 and k2 so that
equation (8.21) is real.
Notice that except for the constant factors, the first term in equation (8.21) is the complex
conjugate of the second. Since the sum of any complex number and its conjugate is the real number
2α, we want to choose the first constant, k1 , to be any complex constant c1 + ic2 and let the second
constant, k2 , be its conjugate pair, c1 − ic2 . Then the first and second term in (8.21) turn out to be
complex conjugates and the sum of them will be a real solution given by
120
A. W. Richter 8.4. SYSTEMS OF LINEAR DIFFERENCE EQUATIONS
is given by
√
xt 1 0
= ( 10)t (c1 cos(tθ ∗ ) − c2 sin(tθ ∗ )) − (c2 cos(tθ ∗ ) + c1 sin(tθ ∗ ))
yt 0 3
√ ∗ ∗
c1 cos(tθ ) − c2 sin(tθ )
= ( 10)t .
−3c2 cos(tθ ∗ ) − 3c1 sin(tθ ∗ )
Remark 8.4.2. In higher dimensions, a given matrix can have both real and complex eigenvalues.
The solution of the corresponding system of difference equations is the obvious combination of the
solutions described in Equation 8.4.1 and Theorem 8.4.2.
121
Bibliography
A DDA , J. AND R. C OOPER (2003): Dynamic Economnics, Cambridge, MA: The MIT Press.
BARRO , R. J. AND X. S ALA - I -M ARTIN (2003): Economic Growth, Cambridge, MA: The MIT
Press, 2nd ed.
B EAVIS , B. AND I. M. D OBBS (1990): Optimization and Stability Theory for Economic Analysis,
Cambridge, MA: Cambridge University Press.
B ROWN , J. AND R. C HURCHILL (2008): Complex Variables and Applications, New York, NY:
McGraw-Hill, 8th ed.
DE LAF UENTE , A. (2000): Mathematical Methods and Models for Economists, Cambridge, MA:
Cambridge University Press.
D IXIT, A. K. (1990): Optimization in Economic Theory, New York, NY: Oxford University Press,
2nd ed.
E DWARDS , C. H. AND D. E. P ENNEY (2008): Differential Equations and Linear Algebra, Upper
Saddle River, NJ: Prentice Hall, 3rd ed.
F RIEDBERG , S. H., A. J. I NSEL , AND L. E. S PENCE (2002): Linear Algebra, Upper Saddle River,
NJ: Prentice Hall, 4th ed.
H AMILTON , J. D. (1994): Time Series Analysis, Princeton, NJ: Princeton University Press.
K. S YDSÆTER , K., A. S TROM , AND P. B ERCK (2005): Economists Mathematical Manual, New
York, NY: Springer, 4th ed.
L AY, S. R. (2004): Analysis with an Introduction to Proof, Upper Saddle River, NJ: Prentice Hall,
4rd ed.
122
A. W. Richter BIBLIOGRAPHY
N OVSHEK , W. (1993): Mathematics for Economists, Cambridge, MA: Emerald Group Publishing
Limited.
S IMON , C. P. AND L. B LUME (1994): Mathematics for Economists, New York, NY: W. W. Norton
and Company.
S TEWART, J. (2011): Calculus: Early Transcendentals, Belmont, CA: Thomson Learning, 7th ed.
S YDSÆTER , K. AND P. H AMMOND (2002): Essential Mathematics for Economic Analysis, Upper
Saddle River, NJ: Prentice Hall.
T URKINGTON , D. A. (2007): Mathematical Tools for Economists, Malden, MA: Blackwell Pub-
lishing.
123