1509 04333v2 PDF

See
discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/281838644
An Introduction to Business Mathematics
Book September 2015

DOI: 10.13140/RG.2.1.3063.0889
CITATIONS READS
0 59,500
1 author:
Henk van Elst

Karlshochschule International University
46 PUBLICATIONS 1,050 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
statsVideos View project
Non-knowledge and uncertainty View project
All content following this page was uploaded by Henk van Elst on 17 September 2015.
The user has requested enhancement of the downloaded file. All in-text references underlined in blue are added to the original document
and are linked to publications on ResearchGate, letting you access and read them immediately.
A N I NTRODUCTION TO
B USINESS M ATHEMATICS
arXiv:1509.04333v2 [q-fin.GN] 16 Sep 2015
Lecture notes for the Bachelor degree programmes

IB/IMC/IMA/ITM/IEVM/ACM/IEM/IMM
at Karlshochschule International University
Module
0.1.1 IMQM: Introduction to Management and its Quantitative Methods
H ENK VAN E LST
August 30, 2015
Fakultat I: Betriebswirtschaft und Management

Karlshochschule
International University
Karlstrae 3638
76133 Karlsruhe
Germany
Email: [email protected]
EPrint: arXiv:1509.04333v2 [q-fin.GN]

c 20092015 Karlshochschule International University and Henk van Elst
Abstract
These lecture notes provide a self-contained introduction to the mathematical methods required in a Bachelor
degree programme in Business, Economics, or Management. In particular, the topics covered comprise real-
valued vector and matrix algebra, systems of linear algebraic equations, Leontiefs stationary inputoutput
matrix model, linear programming, elementary financial mathematics, as well as differential and integral
calculus of real-valued functions of one real variable. A special focus is set on applications in quantitative
economical modelling.
Cite as: arXiv:1509.04333v2 [q-fin.GN]
These lecture notes were typeset in LATEX 2 .

Contents
Abstract
Qualification objectives of the module (excerpt) 1
Introduction 2
1 Vector algebra in Euclidian space Rn 5

1.1 Basic concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2 Dimension and basis of Rn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3 Euclidian scalar product . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 Matrices 11
2.1 Matrices as linear mappings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Basic concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.3 Matrix multiplication . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3 Systems of linear algebraic equations 17

3.1 Basic concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.2 Gauian elimination . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.3 Rank of a matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.4 Solution criteria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.5 Inverse of a regular (n n)-matrix . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.6 Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
4 Leontiefs inputoutput matrix model 25

4.1 General considerations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.2 Inputoutput matrix and resource consumption matrix . . . . . . . . . . . . . . . . 26
4.2.1 Inputoutput matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
4.2.2 Resource consumption matrix . . . . . . . . . . . . . . . . . . . . . . . . 28
4.3 Stationary linear flows of goods . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
4.3.1 Flows of goods: endogenous INPUT to total OUTPUT . . . . . . . . . . . 29
4.3.2 Flows of goods: exogenous INPUT to total OUTPUT . . . . . . . . . . . . 30
4.4 Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
CONTENTS
5 Linear programming 33
5.1 Exposition of a quantitative problem . . . . . . . . . . . . . . . . . . . . . . . . . 33
5.2 Graphical solution method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
5.3 Dantzigs simplex algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
6 Elementary financial mathematics 39

6.1 Arithmetical and geometrical sequences and series . . . . . . . . . . . . . . . . . 39
6.1.1 Arithmetical sequence and series . . . . . . . . . . . . . . . . . . . . . . . 39
6.1.2 Geometrical sequence and series . . . . . . . . . . . . . . . . . . . . . . . 40
6.2 Interest and compound interest . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
6.3 Redemption payments in constant annuities . . . . . . . . . . . . . . . . . . . . . 43
6.4 Pension calculations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
6.5 Linear and declining-balance depreciation methods . . . . . . . . . . . . . . . . . 48
6.5.1 Linear depreciation method . . . . . . . . . . . . . . . . . . . . . . . . . 48
6.5.2 Declining-balance depreciation method . . . . . . . . . . . . . . . . . . . 48
6.6 Summarising formula . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
7 Differential calculus of real-valued functions 51

7.1 Real-valued functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
7.1.1 Polynomials of degree n . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
7.1.2 Rational functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
7.1.3 Power-law functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
7.1.4 Exponential functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
7.1.5 Logarithmic functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
7.1.6 Concatenations of real-valued functions . . . . . . . . . . . . . . . . . . . 54
7.2 Derivation of differentiable real-valued functions . . . . . . . . . . . . . . . . . . 55
7.3 Common functions in economic theory . . . . . . . . . . . . . . . . . . . . . . . . 57
7.4 Curve sketching . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
7.5 Analytic investigations of economic functions . . . . . . . . . . . . . . . . . . . . 59
7.5.1 Total cost functions according to Turgot and von Thunen . . . . . . . . . . 59
7.5.2 Profit functions in the diminishing returns picture . . . . . . . . . . . . . . 62
7.5.3 Extremal values of rational economic functions . . . . . . . . . . . . . . . 64
7.6 Elasticities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
8 Integral calculus of real-valued functions 71

8.1 Indefinite integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
8.2 Definite integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
8.3 Applications in economic theory . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
A Glossary of technical terms (GB D) 75
Bibliography 81
Qualification objectives of the module
(excerpt)
The qualification objectives shall be reached by an integrative approach.

A broad instructive range is aspired. The students shall acquire a 360 degree orientation concerning
the task- and personnel-related tasks and roles of a manager and supporting tools and methods and
be able to describe the coherence in an integrative way. The knowledge concerning the tasks and
the understanding of methods and tools shall be strengthened by a constructivist approach based
on case studies and exercises.
Students who have successfully participated in this module will be able to
...,
solve problems in Linear Algebra and Analysis and apply such mathematical methods to
quantitative problems in management.
apply and challenge the knowledge critically on current issues and selected case studies.
1
2 CONTENTS
Introduction
These lecture notes contain the entire material of the quantitative methods part of the first semester
module 0.1.1 IMQM: Introduction to Management and its Quantitative Methods at Karl-
shochschule International University. The aim is to provide a selection of tried-and-tested math-
ematical tools that proved efficient in actual practical problems of Economics and Management.
These tools constitute the foundation for a systematic treatment of the typical kinds of quantita-
tive problems one is confronted with in a Bachelor degree programme. Nevertheless, they provide
a sufficient amount of points of contact with a quantitatively oriented subsequent Master degree
programme in Economics, Management, or the Social Sciences.
The prerequisites for a proper understanding of these lecture notes are modest, as they do not go
much beyond the basic A-levels standards in Mathematics. Besides the four fundamental arith-
metical operations of addition, subtraction, multiplication and division of real numbers, you should
be familiar, e.g., with manipulating fractions, dealing with powers of real numbers, the binomial
formulae, determining the point of intersection for two straight lines in the Euclidian plane, solv-
ing a quadratic algebraic equation, and the rules of differentiation of real-valued functions of one
variable.
It might be useful for the reader to have available a modern graphic display calculator (GDC) for
dealing with some of the calculations that necessarily arise along the way, when confronted with
specific quantitative problems. Some current models used in public schools and in undergraduate
studies are, amongst others,
Texas Instruments TI84 plus,
Casio CFX9850GB PLUS.
However, the reader is strongly encouraged to think about resorting, as an alternative, to a spread-
sheet programme such as EXCEL or OpenOffice to handle the calculations one encounters in
ones quantitative work.
The central theme of these lecture notes is the acquisition and application of a number of effec-
tive mathematical methods in a business oriented environment. In particular, we hereby focus on
quantitative processes of the sort
INPUT OUTPUT ,
for which different kinds of functional relationships between some numerical INPUT quantities
and some numerical OUTPUT quantities are being considered. Of special interest in this context
3
4 CONTENTS
will be ratios of the structure
OUTPUT
.
INPUT
In this respect, it is a general objective in Economics to look for ways to optimise the value of
such ratios (in favour of some economic agent), either by seeking to increase the OUTPUT when
the INPUT is confined to be fixed, or by seeking to decrease the INPUT when the OUTPUT is
confined to be fixed. Consequently, most of the subsequent considerations in these lecture notes
will therefore deal with issues of optimisation of given functional relationships between some
variables, which manifest themselves either in minimisation or in maximisation procedures.
The structure of these lecture notes is the following. Part I presents selected mathematical methods
from Linear Algebra, which are discussed in Chs. 1 to 5. Applications of these methods focus on
the quantitative aspects of flows of goods in simple economic models, as well as on problems in
linear programming. In Part II, which is limited to Ch. 6, we turn to discuss elementary aspects
of Financial Mathematics. Fundamental principles of Analysis, comprising differential and in-
tegral calculus for real-valued functions of one real variable, and their application to quantitative
economic problems, are reviewed in Part III; this extends across Chs. 7 and 8.
We emphasise the fact that there are no explicit examples nor exercises included in these lecture
notes. These are reserved exclusively for the lectures given throughout term time.
Recommended textbooks accompanying the lectures are the works by Asano (2013) [2], Dowl-
ing (2009) [11], Dowling (1990) [10], Bauer et al (2008) [3], Bosch (2003) [6], and Hulsmann
et al (2005) [16]. Some standard references of Applied Mathematics are, e.g., Bronstein et al
(2005) [7] and Arens et al (2008) [1]. Should the reader feel inspired by the aesthetics, beauty,
ellegance and efficiency of the mathematical methods presented, and, hence, would like to know
more about their background and relevance, as well as being introduced to further mathematical
techniques of interest, she/he is recommended to take a look at the brilliant books by Penrose
(2004) [21], Singh (1997) [23], Gleick(1987) [13] and Smith (2007) [24]. Note that most of the
textbooks and monographs mentioned in this Introduction are available from the library at Karl-
shochschule International University.
Finally, we draw the readers attention to the fact that the *.pdf version of these lecture notes con-
tains interactive features such as fully hyperlinked references to original publications at the web-
sites dx.doi.org and jstor.org, as well as active links to biographical information on sci-
entists that have been influential in the historical development of Mathematics, hosted by the web-
sites The MacTutor History of Mathematics archive (www-history.mcs.st-and.ac.uk)
and en.wikipedia.org.
Chapter 1
Vector algebra in Euclidian space Rn
Let us begin our elementary considerations of vector algebra with the introduction of a special
class of mathematical objects. These will be useful at a later stage, when we turn to formulate
certain problems of a quantitative nature in a compact and elegant way. Besides introducing these
mathematical objects, we also need to define which kinds of mathematical operations they can be
subjected to, and what computational rules we have to take care of.
1.1 Basic concepts

Given be a set V of mathematical objects a which, for now, we want to consider merely as a
collection of n arbitrary real numbers a1 , . . . , ai , . . . , an . In explicit terms,

a1

.
.

.

V = a = ai | ai R, i = 1, . . . , n . (1.1)
.

..

an
Formally the n real numbers considered can either be assembled in an ordered pattern as a column
or a row. We define
Def.: Real-valued column vector with n components

a1
..
.

a := ai , ai R, i = 1, . . . , n , (1.2)
.
..
an
Notation: a Rn1 ,
and
5
6 CHAPTER 1. VECTOR ALGEBRA IN EUCLIDIAN SPACE RN
Def.: Real-valued row vector with n components
aT := (a1 , . . . , ai , . . . , an ) , ai R, i = 1, . . . , n , (1.3)
Notation: aT R1n .
Correspondingly, we define the n-component objects

0
..
.
0T := (0, . . . , 0, . . . , 0)

0 := 0 and (1.4)
.
..
0
to constitute related zero vectors.

Next we define for like objects in the set V , i.e., either for n-component column vectors or for
n-component row vectors, two simple computational operations. These are
Def.: Addition of vectors

a1 b1 a1 + b1
.. .. ...
. .

a + b := ai + bi = ai + bi , ai , bi R , (1.5)
. .. ..
.. . .
an bn an + bn
and
Def.: Rescaling of vectors

a1
..
.

a := ai , , ai R . (1.6)
.
..
an
The rescaling of a vector a with an arbitrary non-zero real number has the following effects:
|| > 1 stretching of the length of a
0 < || < 1 shrinking of the length of a
< 0 directional reversal of a.

1.2. DIMENSION AND BASIS OF RN 7
The notion of the length of a vector a will be made precise shortly.
The addition and the rescaling of n-component vectors satisfy the following addition and multipli-
cation laws:
Computational rules for addition and rescaling of vectors
For vectors a, b, c Rn :
1. a + b = b + a (commutative addition)
2. a + (b + c) = (a + b) + c (associative addition)
3. a + 0 = a (addition identity element)
4. For every a, b Rn , there exists exactly one x Rn such that a + x = b

(invertibility of addition)
5. ()a = (a) with R (associative rescaling)
6. 1a = a (rescaling identity element)
7. (a + b) = a + b;
( + )a = a + a with , R (distributive rescaling).
In conclusion of this section, we remark that every set of mathematical objects V constructed in
line with Eq. (1.1), with an addition and a rescaling defined according to Eqs. (1.5) and (1.6), and
satisfying the laws stated above, constitutes a linear vector space over Euclidian space Rn .1
1.2 Dimension and basis of Rn

Let there be given m n-component vectors2 a1 , . . . , ai , . . . , am Rn , as well as m real numbers
1 , . . . , i , . . . , m R. The new n-component vector b resulting from the addition of arbitrarily
rescaled versions of these m vectors according to
m
X
b = 1 a1 + . . . + i ai + . . . + m am =: i ai R n (1.7)
i=1
is referred to as a linear combination of the m vectors ai , i = 1, . . . , m.

Def.: A set of m vectors a1 , . . . , ai , . . . , am Rn is called linearly independent when the condi-
tion m
!
X
0 = 1 a1 + . . . + i ai + . . . + m am = i ai , (1.8)
i=1
1
This is named after the ancient greek mathematician Euclid of Alexandria (about 325 BC265 BC).
2
A slightly s horter notation for n-component column vectors a Rn1 is given by a Rn ; likewise aT Rn
for n-component row vectors aT R1n .
i.e., the problem of forming the zero vector 0 Rn from a linear combination of the m vectors
a1 , . . . , ai , . . . , am Rn , can only be solved trivially, namely by 0 = 1 = . . . = i = . . . = m .
When, however, this condition can be solved non-trivially, with some i 6= 0, then the set of m
vectors a1 , . . . , ai , . . . , am Rn is called linearly dependent.
In Euclidian space Rn , there is a maximum number n (!) of vectors which can be linearly indepen-
dent. This maximum number is referred to as the dimension of Euclidian space Rn . Every set
of n linearly independent vectors in Euclidian space Rn constitutes a possible basis of Euclidian
space Rn . If the set {a1 , . . . , ai , . . . , an } constitutes a basis of Rn , then every other vector b Rn
can be expressed in terms of these basis vectors by
n
X
b = 1 a1 + . . . + i ai + . . . + n an = i ai . (1.9)
i=1
The rescaling factors i R of the ai Rn are called the components of vector b with respect
to the basis {a1 , . . . , ai , . . . , an }.
Remark: The n unit vectors

1 0 0
0 1 0

e1 := .. , e2 := .. , ... , en := .. , (1.10)
. . .
0 0 1
constitute the so-called canonical basis of Euclidian space Rn . With respect to this basis, all
vectors b Rn can be represented as a linear combinationen

b1
b2 Xn

b = .. = b1 e1 + b2 e2 + + bn en = bi ei . (1.11)
. i=1
bn
1.3 Euclidian scalar product

Finally, to conclude this section, we introduce a third mathematical operation defined for vectors
on Rn .
Def.: For an n-component row vector aT R1n and an n-component column vector b Rn1 ,
the Euclidian scalar product

b1
..
. Xn
T
a b := (a1 , . . . , ai , . . . an ) bi = a1 b1 + . . . + ai bi . . . + an bn =: ai bi (1.12)
.
.. i=1
bn
1.3. EUCLIDIAN SCALAR PRODUCT 9
defines a mapping f : R1n Rn1 R from the product set of n-component row and column
vectors to the set of real numbers. Note that, in contrast to the addition and the rescaling of n-
component vectors, the outcome of forming a Euclidian scalar product between two n-component
vectors is a single real number.
In the context of the Euclidian scalar product, two non-zero vectors a, b Rn (wit a 6= 0 =
6 b)
T T
are referred to as mutually orthogonal when they exhibit the property that 0 = a b = b a.
Computational rules for Euclidian scalar product of vectors
For vectors a, b, c Rn :
1. (a + b)T c = aT c + bT c (distributive scalar product)
2. aT b = bT a (commutative scalar product)
3. (aT ) b = (aT b) with R (homogeneous scalar product)
4. aT a > 0 for all a 6= 0 (positive definite scalar product).
Now we turn to introduce the notion of the length of an n-component vector.

Def.: The length of a vector a Rn is defined via the Euclidian scalar product as
v
u n
q uX
|a| := aT a = a21 + . . . + a2i + . . . + a2n =: t a2i . (1.13)
i=1
Technically one refers to the non-negative real number |a| as the absolute value or the Euclidian
norm of the vector a Rn . The length of a Rn has the following properties:
|a| 0, and |a| = 0 a = 0;
|a| = |||a| for R;
|a + b| |a| + |b| (triangle inequality).
Every non-zero vector a Rn , i.e., |a| > 0, can be rescaled by the reciprocal of its length. This
procedure defines the
Def.: Normalisation of a vector a Rn ;
a
a := |a| = 1 . (1.14)
|a|
By this method one generates a vector of length 1, i.e., a unit vector a. To denote unit vectors we
will employ the hat symbol.
Lastly, also by means of the Euclidian scalar product, we introduce the angle enclosed between
two non-zero vectors.
Def.: Angle enclosed between a, b 6= 0 Rn
aT b
cos[(a, b)] = = aT b (a, b) = cos1 (aT b) . (1.15)
|a| |b|
Remark: The inverse cosine function3 cos1 (. . .) is available on every standard GDC or spread-
sheet.
3
The notion of on inverse function will be discussed later in Ch. 7.
Chapter 2
Matrices
In this chapter, we introduce a second class of mathematical objects that are more general than
vectors. For these objects, we will also define certain mathematical operations, and a set of com-
putational rules that apply in this context.
2.1 Matrices as linear mappings

Consider given a collection of m n arbitrary real numbers a11 , a12 . . . , aij , . . . , amn , which we
arrange systematically in a particular kind of array.
Def.: A real-valued (m n)-matrix is formally defined to constitute an array of real numbers

according to

a11 a12 . . . a1j . . . a1n
a21 a22 . . . a2j . . . a2n
. .. .. ..
. .. ..
. . . . . .
A := , (2.1)
ai1 ai2 . . . aij . . . ain
. .. .. .. .. ..
.. . . . . .
am1 am2 . . . amj . . . amn
where aij R, i = 1, . . . , m; j = 1, . . . , n.
Notation: A Rmn .
Characteristic features of this array of real numbers are:
m denotes the number of rows of A, n the number of columns of A.
aij represents the elements of A; aij is located at the point of intersection of the ith row and
the jth column of A.
elements of the ith row constitute the row vector (ai1 , ai2 , . . . , aij , . . . , ain ), elements of the
11
12 CHAPTER 2. MATRICES

a1j
a2j
.
.
.
jth column the column vector .
aij
.
..
amj
Formally column vectors need to be viewed as (n 1)-matrices, row vectors as (1 n)-matrices.
An (m n)-zero matrix, denoted by 0, has all its elements equal to zero, i.e.,

0 0 ... 0
0 0 ... 0

0 := .. .. . . .. . (2.2)
. . . .
0 0 ... 0
Matrices which have an equal number of rows and columns, i.e. m = n, are referred to as
quadratic matrices. In particular, the (n n)-unit matrix (or identity matrix)

1 0 ... 0 ... 0
0 1 ... 0 ... 0
. .
. . . . . .. . . . ..
. . . .
1 := (2.3)
0 0 . . . 1 . . . 0
. . .
.. .. . . ... . . . ...

0 0 ... 0 ... 1
holds a special status in the family of (n n)-matrices.

Now we make explicit in what sense we will comprehend (m n)-matrices as mathematical
objects.
Def.: A real-valued matrix A Rmn defines by the computational operation

a11 a12 . . . a1j . . . a1n x1
a21 a22 . . . a2j . . . a2n x2
. .. .. ..
. .. .. .
. . . . . . ..
Ax :=
ai1 ai2 . . . aij . . . ain xj
. .. .. .. .. .. .
.. . . . . . ..
am1 am2 . . . amj . . . amn xn

a11 x1 + a12 x2 + . . . + a1j xj + . . . + a1n xn y1
a21 x1 + a22 x2 + . . . + a2j xj + . . . + a2n xn y2
.. .
.
. .
:= =: =y (2.4)
ai1 x1 + ai2 x2 + . . . + aij xj + . . . + ain xn yi

.
..

.
..

am1 x1 + am2 x2 + . . . + amj xj + . . . + amn xn ym
2.2. BASIC CONCEPTS 13
a mapping A : Rn1 Rm1 , i.e. a mapping from the set of real-valued n-component column
vectors (here: x) to the set of real-valued m-component column vectors (here: y).
In loose analogy to the photographic process, x can be viewed as representing an object, A a
camera, and y the resultant image.
Since for real-valued vectors x1 , x2 Rn1 and real numbers R, mappings defined by real-
valued matrices A Rmn exhibit the two special properties
A(x1 + x2 ) = (Ax1 ) + (Ax2 )

(2.5)
A(x1 ) = (Ax1 ) ,
they constitute linear mappings.1

We now turn to discuss the most important mathematical operations defined for (m n)-matrices,
as well as the computational rules that obtain.
2.2 Basic concepts

Def.: Transpose of a matrix
For A Rmn , we define the process of transposing A by
AT : aTij := aji , (2.6)
where i = 1, . . . , m und j = 1, . . . , n. Note that it holds that AT Rnm .

When transposing an (mn)-matrix, one simply has to exchange the matrix rows with its columns
(and vice versa): the elements of the first row become the elements of the first column, etc. It
follows that, in particular,
(AT )T = A (2.7)
applies.
Two special cases may occur for quadratic matrices (where m = n):
When AT = A, one refers to A as a symmetric matrix.
When AT = A, one refers to A as an antisymmetric matrix.
Def.: Addition of matrices

For A, B Rmn , the sum is given by
A + B =: C : aij + bij =: cij , (2.8)
with i = 1, . . . , m and j = 1, . . . , n.
Note that an addition can be performed meaningfully only for matrices of the same format.
1
It is important to note at this point that many advanced mathematical models designed to describe quantitative
aspects of some natural and economic phenomena do not satisfy the conditions (2.5), as they employ non-linear
mappings for this purpose. However, in such contexts, linear mappings often provide useful first approximations.
Def.: Rescaling of matrices
For A Rmn and R\{0}, let
A =: C : aij =: cij , (2.9)
where i = 1, . . . , m and j = 1, . . . , n.
When rescaling a matrix, all its elements simply have to be multiplied by the same non-zero real
number .
Computational rules for addition and rescaling of matrices
For matrices A, B, C Rmn :
1. A + B = B + A (commutative addition)
2. A + (B + C) = (A + B) + C (associative addition)
3. A + 0 = A (addition identity element)
4. For every A and B, there exists exactly one Z such that A + Z = B.

(invertibility of addition)
5. ()A = (A) with , R\{0} (associative rescaling)
6. 1A = A (rescaling identity element)
7. (A + B) = A + B;
( + )A = A + A with , R\{0} (distributive rescaling)
8. (A + B)T = AT + BT (transposition rule 1)
9. (A)T = AT with R\{0}. (transposition rule 2)
Next we introduce a particularly useful mathematical operation for matrices.
2.3 Matrix multiplication

Def.: For a real-valued (m n)-matrix A and a real-valued (n r)-matrix B, a matrix multipli-
cation is defined by
AB =: C
Pn (2.10)
ai1 b1j + . . . + aik bkj + . . . + ain bnj =: k=1 aik bkj =: cij ,
with i = 1, . . . , m and j = 1, . . . , r, thus yielding as an outcome a real-valued (m r)-matrix C.

2.3. MATRIX MULTIPLICATION 15
The element of C at the intersection of the ith row and the jth column is determined by the
computational rule
cij = Euclidian scalar product of ith row vector of A and jth column vector of B . (2.11)
It is important to realise that the definition of a matrix multiplication just provided depends in
an essential way on the fact that matrix A on the left in the product needs to have as many (!)
columns as matrix B on the right rows. Otherwise, a matrix multiplication cannot be defined in a
meaningful way.
GDC: For matrices [A] and [B] edited beforehand, of matching formats, their matrix multiplication
can be evaluated in mode MATRIX NAMES by [A] [B].
Computational rules for matrix multiplication

For A, B, C real-valued matrices of correspondingly matching formats we have:
1. AB = 0 is possible with A 6= 0, B 6= 0. (zero divisor)
2. A(BC) = (AB)C (associative matrix multiplication)
1 A=A
1 = |{z}
3. A |{z} (multiplicative identity element)
Rnn Rmm
4. (A + B)C = AC + BC
C(A + B) = CA + CB (distributive matrix multiplication)
5. A(B) = (A)B = (AB) with R (homogeneous matrix multiplication)
6. (AB)T = BT AT (transposition rule).

Chapter 3
Systems of linear algebraic equations
In this chapter, we turn to address a particular field of application of the notions of matrices and
vectors, or of linear mappings in general.
3.1 Basic concepts

Let us begin with a system of m N linear algebraic equations, wherein every single equation can
be understood to constitute a constraint on the range of values of n N variables x1 , . . . , xn R.
The objective is to determine all possible values of x1 , . . . , xn R which satisfy these constraints
simultaneously. Problems of this kind, namely systems of linear algebraic equations, are often
represented in the form
Representation 1:
a11 x1 + . . . + a1j xj + . . . + a1n xn = b1
..
.
ai1 x1 + . . . + aij xj + . . . + ain xn = bi (3.1)
..
.
am1 x1 + . . . + amj xj + . . . + amn xn = bm .
Depending on how the natural numbers m and n relate to one another, systems of linear algebraic
equations can be classified as follows:
m < n: fewer equations than variables; the linear system is under-determined,
m = n: same number of equations as variables; the linear system is well-determined,
m > n: more equations than variables; the linear system is over-determined.
A more compact representation of a linear system of format (m n) is given by
17
18 CHAPTER 3. SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS
Representation 2:

a11 . . . a1j . . . a1n x1 b1
.. ..
.
.. ..
.
.. .. ..
. . . . .

Ax = ai1 . . . aij . . . ain xj = bi = b . (3.2)
. .. .. .. .. . ..
.. . . . . .. .
am1 . . . amj . . . amn xn bm
The mathematical objects employed in this variant of a linear system are as follows: A takes the
central role of the coefficient matrix of the linear system, of format (m n), x is its variable
vector, of format (n 1), and, lastly, b is its image vector, of format (m 1).
When dealing with systems of linear algebraic equations in the form of Representation 2, i.e.
Ax = b, the main question to be answered is:
Question: For given coefficient matrix A and image vector b, can we find a variable vector x
that A maps onto b?
In a sense this describes the inversion of the photographic process we had previously referred to:
we have given the camera and we already know the image, but we have yet to find a matching
object. Remarkably, to address this issue, we can fall back on a simple algorithmic method due to
the German mathematician and astronomer Carl Friedrich Gau (17771855).
3.2 Gauian elimination

The algorithmic solution technique developed by Gau is based on the insight that the solution set
of a linear system of m algebraic equations for n real-valued variables, i.e.
Ax = b , (3.3)
remains unchanged under the following algebraic equivalence transformations of the linear sys-
tem:
1. changing the order amongst the equations,
2. multiplication of any equation by a non-zero real number c 6= 0,
3. addition of a multiple of one equation to another equation,
4. changing the order amongst the equations.
Specifically, this implies that we may manipulate a given linear system by means of these four
different kinds of equivalence transformations without ever changing its identity. In concrete cases,
however, one should not apply these equivalence transformations at random but rather follow a
target oriented strategy. This is what Gauian elimination can provide.
3.3. RANK OF A MATRIX 19
Target: To cast the augmented coefficient matrix (A|b), i.e., the array
a11 . . . a1j . . . a1n b1

.. .. .. .. .. ..
. . . . . .
ai1 . . . aij . . . ain bi , (3.4)
.. .. .. .. .. ..
. . . . . .
am1 . . . amj . . . amn bm
when possible, into upper triangular form
1 . . . a1j . . . a1n b1
.. . . . .. .. ..
. . .. . . .
0 . . . aij . . . ain bi , (3.5)
.. . . . .. .. ..
. . .. . . .
0 . . . 0 . . . amn bm
by means of the four kinds of equivalence transformations such that the resultant simpler final
linear system may easily be solved using backward substitution.
Three exclusive cases of possible solution content for a given system of linear algebraic equations
do exist. The linear system may possess either
1. no solution at all, or
2. a unique solution, or
3. multiple solutions.
Remark: Linear systems that are under-determined, i.e., when m < n, can never be solved
uniquely due to the fact that in such a case there not exist enough equations to constrain the values
of all of the n variables.
GDC: For a stored augmented coefficient matrix [A] of format (m n + 1), associated with a given
(m n) linear system, select mode MATRIX MATH and then call the function rref([A]). It is
possible that backward substitution needs to be employed to obtain the final solution.
For completeness, we want to turn briefly to the issue of solvability of a system of linear algebraic
equations. To this end, we need to introduce the notion of the rank of a matrix.
3.3 Rank of a matrix

Def.: A real-valued matrix A Rmn possesses the rank
rank(A) = r , r min{m, n} (3.6)
if and only if r is the maximum number of row resp. column vectors of A which are linearly
independent. Clearly, r can only be as large as the smaller of the numbers m and n that determine
the format of A.
For quadratic matrices A Rnn , there is available a more elegant measure to determine its
rank. This (in the present case real-valued) measure is referred to as the determinant of matrix A,
det(A), and is defined as follows.
Def.:
(i) When A R22 , its determinant is given by

a11 a12
det(A) := := a11 a22 a12 a21 , (3.7)
a21 a22
i.e. the difference between the products of As on-diagonal elements and As off-diagonal
elements.
(ii) When A R33 , the definition of As determinant is more complex. In that case it is
given by

a11 a12 a13
det(A) := a21 a22 a23
a31 a32 a33
:= a11 (a22 a33 a32 a23 ) + a21 (a32 a13 a12 a33 ) + a31 (a12 a23 a22 a13 ) (3.8)
.
Observe, term by term, the cyclic permutation of the first index of the elements aij according
to the rule 1 2 3 1.
(iii) Finally, for the (slightly involved) definition of the determinant of a higher-dimensional
matrix A Rnn , please refer to the literature; e.g. Bronstein et al (2005) [7, p 267].
To determine the rank of a given quadtratic matrix A Rnn , one now installs the following
criteria: rank(A) = r = n, if det(A) 6= 0, and rank(A) = r < n, if det(A) = 0. In the first case,
A is referred to as regular, in the second as singular. For quadratic matrices A that are singular,
rank(A) = r (with r < n) is given by the number r of rows (or columns) of the largest possible
non-zero subdeterminant of A.
GDC: For a stored quadratic matrix [A], select mode MATRIX MATH and obtain its determinant
by calling the function det([A]).
3.4 Criteria for solving systems of linear algebraic equations

Making use of the concept of the rank of a real-valued matrix A Rmn , we can now summarise
the solution content of a specific system of linear algebraic equations of format (m n) in a table.
For given linear system
Ax = b ,
with coefficient matrix A Rmn , variable vector x Rn1 and image vector b Rm1 , there
exist(s)
3.5. INVERSE OF A REGULAR (N N)-MATRIX 21
b 6= 0 b=0
1. rank(A) 6= rank(A|b) no solution
2. rank(A) = rank(A|b) = r
(a) r=n a unique x=0

solution
(b) r<n multiple multiple

solutions: solutions:
n r free n r free
parameters parameters
(A|b) here denotes the augmented coefficient matrix.
Next we discuss a particularly useful property of regular quadratic matrices.
3.5 Inverse of a regular (n n)-matrix

Def.: Let a real-valued quadratic matrix A Rnn be regular, i.e., det(A) R\{0}. Then there
exists an inverse matrix A1 to A defined by the characterising properties
A1 A = AA1 = 1 . (3.9)
Here 1 denotes the (n n)-unit matrix [cf. Eq. (2.3)].
When a computational device is not at hand, the inverse matrix A1 of a regular quadratic matrix
A can be obtained by solving the matrix-valued linear system
!
AX = 1 (3.10)
for the unknown matrix X by means of simultaneous Gauian elimination.
GDC: For a stored quadratic matrix [A], its inverse matrix can be simply obtained as [A]1 , where
the x1 function key needs to be used.
Computational rules for the inverse operation
For A, B Rnn , with det(A) 6= 0 6= det(B), it holds that
1. (A1 )1 = A
2. (AB)1 = B1 A1
3. (AT )1 = (A1 )T
1 1
4. (A)1 = A .

The special interest in applications in the concept of inverse matrices arises for the following
reason. Consider given a well-determined linear system
Ax = b ,
with regular quadratic coefficient matrix A Rnn , i.e., det(A) 6= 0. Then, for A, there exists an
inverse matrix A1 . Matrix-multiplying both sides of the equation above from the left (!) by the
inverse A1 , results in
A1 (Ax) = (A1 A)x = 1x = x = 1

|A{z b} . (3.11)
| {z }
left-hand side right-hand side
In this case, the unique solution (!) x = A1b of the linear system arises simply from matrix
multiplication of the image vector b by the inverse matrix of A. (Of course, it might actually
require a bit of computational work to determine A1.)
3.6 Outlook
There are a number of exciting advanced topics in Linear Algebra. Amongst them one finds
the concept of the characteristic eigenvalues and associated eigenvectors of quadratic matrices,
which has particularly high relevance in practical applications. The question to be answered here
is the following: for given real-valued quadratic matrix A Rnn , do there exist real numbers
n R and real-valued vectors v n Rn1 which satisfy the condition
!
Av n = n v n ? (3.12)
Put differently: for which vectors v n Rn1 does their mapping by a quadratic matrix A Rnn
amount to simple rescalings by real numbers n R?
By re-arranging, Eq. (3.12) can be recast into the form
!
0 = (A n 1) v n , (3.13)
3.6. OUTLOOK 23
with 1 an (n n)-unit matrix [cf. Eq. (2.3)] and 0 an n-component zero vector. This condition
corresponds to a homogeneous system of linear algebraic equations of format (n n). Non-trivial
solutions v n 6= 0 to this system exist provided that the so-called characteristic equation
!
0 = det (A n 1) , (3.14)
a polynomial of degree n (cf. Sec. 7.1.1), allows for real-valued roots n R. Note that symmetric
quadratic matrices (cf. Sec. 2.2) possess exclusively real-valued eigenvalues n . When these
eigenvalues turn out to be all different, then the associated eigenvectors v n prove to be mutually
orthogonal.
Knowledge of the spectrum of eigenvalues n R and associated eigenvectors v n Rn1
of a real-valued matrix A Rnn provides the basis of a transformation of A to its diagonal
form An , thus yielding a diagonal matrix which features the eigenvalues n as its on-diagonal
elements; cf. Leon (2009) [19].
Amongst other examples, the concept of eigenvalues and eigenvectors of quadratic real-valued
matrices plays a special role in Statistics, in the context of exploratory principal component
analyses of multivariate data sets, where the objective is to identify dominant intrinsic structures;
cf. Hair et al (2010) [14, Ch. 3] and Ref. [12, App. A].
Chapter 4
Leontiefs stationary inputoutput matrix

model
We now turn to discuss some specific applications of Linear Algebra in economic theory. To
begin with, let us consider quantitative aspects of the exchange of goods between a certain number
of economic agents. We here aim at a simplified abstract description of real economic processes.
4.1 General considerations

The quantitative model to be described in the following is due to the Russian
economist Wassily Wassilyovich Leontief (19051999), cf. Leontief (1936) [20],
for which, besides other important contributions, he was awarded the 1973
Sveriges Riksbank Prize in Economic Sciences in Memory of Alfred Nobel.
Suppose given an economic system consisting of n N interdependent economic agents ex-
changing between them the goods they produce. For simplicity we want to assume that every one
of these economic agents represents the production of a single good only. Presently we intend
to monitor the flow of goods in this simple economic system during a specified reference period
of time. The total numbers of the n goods leaving the production sector of this model constitute
the OUTPUT quantities. The INPUT quantities to the production sector are twofold. On the
one hand, there are exogenous INPUT quantities which we take to be given by m N different
kinds of external resources needed in differing proportions to produce the n goods. On the other
hand, due to their mutual interdependence, some of the economic agents require goods made by
their neighbours to be able to produce their own goods; these constitute the endogenous INPUT
quantities of the system. Likewise, the production sectors total OUTPUT during the chosen ref-
erence period of the n goods can be viewed to flow through one of two separate channels: (i) the
exogenous channel linking the production sector to external consumers representing an open
market, and (ii) the endogenous channel linking the economic agents to their neighbours (thus
respresenting their interdependencies). It is supposed that momentum is injected into the economic
system, triggering the flow of goods between the different actors, by the prospect of increasing
25
26 CHAPTER 4. LEONTIEFS INPUTOUTPUT MATRIX MODEL
the value of the INPUT quantities, in line with the notion of the economic value chain.
Leontiefs model is based on the following three elementary
Assumptions:
1. For all goods involved the functional relationship between INPUT and OUTPUT quantities
be of a linear nature [cf. Eq. (2.5)].
2. The proportions of INPUT quantities to OUTPUT quantities be constant over the refer-
ence period of time considered; the flows of goods are thus considered to be stationary.
3. Economic equilibrium obtains during the reference period of time: the numbers of goods
then supplied equal the numbers of goods then demanded.
The mathematical formulation of Leontiefs quantitative model employs the following

Vector- and matrix-valued quantities:
1. q total output vector Rn1 , components qi 0 units (dim: units)
2. y final demand vector Rn1 , components yi 0 units (dim: units)
3. P inputoutput matrix Rnn , components Pij 0 (dim: 1)
4. (1 P) technology matrix Rnn , regular, hence, invertible (dim: 1)
5. (1 P)1 total demand matrix Rnn (dim: 1)
6. v resource vector Rm1 , components vi 0 units (dim: units)
7. R resource consumption matrix Rmn , components Rij 0, (dim: 1)
where 1 denotes the (n n)-unit matrix [cf. Eq. (2.3)]. Note that the components of all the
vectors involved, as well as of the inputoutput matrix and of the resource consumption matrix,
can assume non-negative values (!) only.
4.2 Inputoutput matrix and resource consumption matrix

We now turn to provide the definition of the two central matrix-valued quantities in Leontiefs
model. We will also highlight their main characteristic features.
4.2.1 Inputoutput matrix

Suppose the reference period of time has ended for the economic system in question, i.e. the
stationary flows of goods have stopped eventually. We now want to take stock of the numbers of
goods that have been delivered by each of the n economic agents in the system. Say that during
the reference period considered, agent 1 delivered of their good the number n11 to themselves, the
number n12 to agent 2, the number n13 to agent 3, and so on, and, lastly, the number n1n to agent n.
4.2. INPUTOUTPUT MATRIX AND RESOURCE CONSUMPTION MATRIX 27
The number delivered by agent 1 to external consumers shall be denoted by y1 . Since in this model
a good produced cannot all of a sudden disappear again, and since by Assumption 3 above the
number of goods supplied must be equal to the number of goods demanded, we find that for the
total output of agent 1 it holds that q1 := n11 + . . . + n1j + . . . + n1n + y1 . Analogous relations
hold for the total output q2 , q3 , . . . , qn of each of the remaining n 1 agents. We thus obtain the
intermediate result
q1 = n11 + . . . + n1j + . . . + n1n + y1 > 0 (4.1)

..
.
qi = ni1 + . . . + nij + . . . + nin + yi > 0 (4.2)
..
.
qn = nn1 + . . . + nnj + . . . + nnn + yn > 0 . (4.3)
This simple system of balance equations can be summarised in terms of a standard inputoutput
table as follows:
[Values in units] agent 1 agent j agent n external consumers : total output

agent 1 n11 ... n1j ... n1n y1 q1
.. .. .. .. .. .. .. ..
. . . . . . . .
agent i ni1 ... nij ... nin yi qi
.. .. .. .. .. .. .. ..
. . . . . . . .
agent n nn1 ... nnj ... nnn yn qn
The first column of this table lists all the n different sources of flows of goods (or suppliers of
goods), while its first row shows the n + 1 different sinks of flows of goods (or consumers of
goods). The last column contains the total output of each of the n agents in the reference period
of time.
Next we compute for each of the n agents the respective values of the non-negative ratios
INPUT (in units) of agent i for agent j (during reference period)

Pij := , (4.4)
OUTPUT (in units) of agent j (during reference period)
or, employing a compact and economical index notation,1
nij
Pij := , (4.5)
qj
1
Note that the normalisation quantities in these ratios Pij are given by the total output qj of the receiving agent j
and not by the total output qi of the supplying agent i. In the latter case the Pij would represent percentages of the
total output qi .
with i, j = 1, . . . , n. These n n = n2 different ratios may be naturally viewed as the elements of
a quadratic matrix P of format (n n). In general, this matrix is given by
n11 n1j
n11 +...+n1j +...+n1n +y1
. . . nj1 +...+njj +...+n jn +yj
. . . nn1 +...+nnjn1n
+...+nnn +yn
.. .. .. .. ..
. . . . .
ni1 nij nin

P = n11 +...+n1j +...+n1n +y1 . . . nj1 +...+njj +...+njn +yj . . . nn1 +...+nnj +...+nnn +yn

,
.. . .. .
.. . .. .
..
.
nn1 nnj nnn
n11 +...+n1j +...+n1n +y1
. . . nj1 +...+njj +...+njn +yj . . . nn1 +...+nnj +...+nnn +yn
(4.6)
and is referred to as Leontiefs inputoutput matrix of the stationary economic system under
investigation.
For the very simple case with just n = 3 producing agents, the inputoutput matrix reduces to
n11 n12 n13
n11 +n12 +n13 +y1 n21 +n22 +n23 +y2 n31 +n32 +n33 +y3
n21 n22 n23
P= n11 +n12 +n13 +y1 n21 +n22 +n23 +y2 n31 +n32 +n33 +y3
.
n31 n32 n33
n11 +n12 +n13 +y1 n21 +n22 +n23 +y2 n31 +n32 +n33 +y3
It is important to realise that for an actual economic system the inputoutput matrix P can be
determined only once the reference period of time chosen has come to an end.
The utility of Leontiefs stationary inputoutput matrix model is in its application for the purpose
of forecasting. This is done on the basis of an extrapolation, namely by assuming that an input
output matrix Preference period obtained from data taken during a specific reference period also is valid
(to an acceptable degree of accuracy) during a subsequent period, i.e.,
Psubsequent period Preference period , (4.7)
or, in component form,

nij nij
Pij |subsequent period = Pij |reference period = . (4.8)
qj subsequent period qj reference period
In this way it becomes possible to compute for a given (idealised) economic system approximate
numbers of INPUT quantities required during a near future production period from the known
numbers of OUTPUT quantities of the most recent production period. Long-term empirical ex-
perience has shown that this method generally leads to useful results to a reasonable approximation.
All of these calculations are grounded on linear relationships describing the quantitative aspects of
stationary flows of goods, as we will soon elucidate.
4.2.2 Resource consumption matrix

The second matrix-valued quantity central to Leontiefs stationary model is the resource con-
sumption matrix R. This may be interpreted as providing a recipe for the amounts of the m dif-
ferent kinds of external resources (the exogenous INPUT quantities) that are needed in the pro-
duction of the n goods (the OUTPUT quantities). Its elements are defined as the ratios
Rij := amounts (in units) required of resource i for the production of one unit of good j ,
(4.9)
4.3. STATIONARY LINEAR FLOWS OF GOODS 29
with i = 1, . . . , m und j = 1, . . . , n. The rows of matrix R thus contain information concerning the
m resources, the columns information concerning the n goods. Note that in general the (m n)
resource consumption matrix R is not (!) a quadratic matrix and, therefore, in general not
invertible.
4.3 Stationary linear flows of goods

4.3.1 Flows of goods: endogenous INPUT to total OUTPUT
We now turn to a quantitative description of the stationary flows of goods that are associated
with the total output q during a specific period of time considered. According to Leontiefs
Assumption 1, there exists a linear functional relationship between the endogenous vector-valued
INPUT quantity q y and the vector-valued OUTPUT quantity q. This may be represented in
terms of a matrix-valued relationship as
n
X
q y = Pq qi yi = Pij qj , (4.10)
j=1
with i = 1, . . . , n, in which the inputoutput matrix P takes the role of mediating a mapping
between either of these vector-valued quantities. According to Assumption 2, the elements of the
inputoutput matrix P remain constant for the period of time considered, i.e. the corresponding
flows of goods are assumed to be stationary.
Relation (4.10) may also be motivated from an alternative perspective that takes the physical sci-
ences as a guidline. Namely, the total numbers q of the n goods produced during the period of
time considered which, by Assumption 3, are equal to the numbers supplied of the n goods satisfy
a conservation law: whatever has been produced of the n goods during the period of time con-
sidered cannot get lost in this period. In quantitative terms this simple relationship may be cast
into the form
q = y + Pq .
|{z} |{z} |{z}
total output final demand (exogenous) deliveries to production sector (endogenous)
For computational purposes this central stationary flow of goods relation (4.10) may be rearranged
as is convenient. In this context it is helpful to make use of the matrix identity q = 1q, where 1
denotes the (n n)-unit matrix [cf. Eq. (2.3)].
Examples:
(i) given/known: P, q
Then it applies that
n
X
y = (1 P)q yi = (ij Pij )qj , (4.11)
j=1
with i = 1, . . . , n; (1 P) represents the invertible technology matrix of the economic

system regarded.
(ii) given/known: P, y
Then it holds that

n
X
1
q = (1 P) y qi = (ij Pij )1 yj , (4.12)
j=1
with i = 1, . . . , n; (1 P)1 here denotes the total demand matrix, i.e., the inverse of the
technology matrix.
4.3.2 Flows of goods: exogenous INPUT to total OUTPUT

Likewise, by Assumption 1, a linear functional relationship is supposed to exist between the ex-
ogenous vector-valued INPUT quantity v and the vector-valued OUTPUT quantity q. In matrix
language this can be expressed by
n
X
v = Rq vi = Rij qj , (4.13)
j=1
with i = 1, . . . , m. By Assumption 2, the elements of the resource consumption matrix R remain

constant during the period of time considered, i.e., the corresponding resource flows are supposed
to be stationary.
By combination of Eqs. (4.13) and (4.12), it is possible to compute the numbers v of resources
required (during the period of time considered) for the production of the n goods for given final
demand y. It applies that
n X
X n
1
v = Rq = R(1 P) y vi = Rij (jk Pjk )1 yk , (4.14)
j=1 k=1
with i = 1, . . . , m.
GDC: For problems with n 5, and known matrices P and R, Eqs. (4.11), (4.12) and (4.14) can
be immediately used to calculate the quantities q from given quantities y, or vice versa.
4.4 Outlook
Leontiefs inputoutput matrix model may be extended in a straightforward fashion to include
more advanced considerations of economic theory. Supposing a closed though not necessarily
stationary economic system G comprising n interdependent economic agents producing n differ-
ent goods, one may assign monetary values to the INPUT quantity v as well as to the OUTPUT
quantities q and y of the system. Besides the numbers of goods produced and the associated
flows of goods one may monitor with respect to G for a given period of time, one can in addition
analyse in time and space the amount of money coupled to the different goods, and the corre-
sponding flows of money. However, contrary to the number of goods, in general there does not
4.4. OUTLOOK 31
exist a conservation law for the amount of money with respect to G. This may render the analysis
of flows of money more difficult, because, in the sense of an increase in value, money can either
be generated inside G during the period of time considered or it can likewise be annihilated; it is
not just limited to either flowing into respectively flowing out of G. Central to considerations of
this kind is a balance equation for the amount of money contained in G during a given period of
time, which is an additive quantity. Such balance equations constitute familiar tools in Physics (cf.
Herrmann (2003) [15, p 7ff]). Its structure in the present case is given by2

rate of change in time
of the amount of money = flux of money rate of generation of money
+ .
into G [in CU/TU] in G [in CU/TU]
in G [in CU/TU]
Note that, with respect to G, both fluxes of money and rates of generation of money can in prin-
ciple possess either sign, positive or negative. To deal with these quantitative issues properly,
one requires the technical tools of the differential and integral calculus which we will discuss
at an elementary level in Chs. 7 and 8. We make contact here with the interdisciplinary science
of Econophysics (cf., e.g., Bouchaud and Potters (2003) [5]), a very interesting and challenging
subject which, however, is beyond the scope of these lecture notes.
Leontiefs inputoutput matrix model, and its possible extension as outlined here, provide the
quantitative basis for considerations of economical ratios of the kind
OUTPUT [in units]
,
INPUT [in units]
as mentioned in the Introduction. In addition, dimensionless (scale-invariant) ratios of the form
REVENUE [in CU]
,
COSTS [in CU]
referred to as economic efficiency, can be computed for and compared between different economic
systems and their underlying production sectors. In Ch. 7 we will briefly reconsider this issue.
2
Here the symbols CU and TU denote currency units and time units, respectively.
Chapter 5
Linear programming
On the backdrop of the economic principle, we discuss in this chapter a special class of quanti-
tative problems that frequently arise in specific practical applications in Business and Manage-
ment. Generally one distinguishes between two variants of the economic principle: either (i) to
draw maximum utility from limited resources, or (ii) to reach a specific target with minimum ef-
fort (costs). With regard to the ratio (OUTPUT)/(INPUT) put into focus in the Introduction, the
issue is to find an optimal value for this ratio under given boundary conditions. This aim can be
realised either (i) by increasing the (positive) value of the numerator for fixed (positive) value of
the denominator, or (ii) by decreasing the (positive) value of the denominator for fixed (positive)
value of the numerator. The class of quantitative problems to be looked at in some detail in this
chapter typically relate to boundary conditions according to case (i).
5.1 Exposition of a quantitative problem
To be maximised is a (non-negative) real-valued quantity z, which depends in a linear functional

fashion on a fixed number of n (non-negative) real-valued variables x1 , . . . , xn . We suppose that
the n variables x1 , . . . , xn in turn are constrained by a fixed number m of algebraic conditions,
which also are assumed to depend on x1 , . . . , xn in a linear fashion. These m constraints, or
restrictions, shall have the character of imposing upper limits on m different kinds of resources.
Def.: Consider a matrix A Rmn , a vector b Rm1 , two vectors c, x Rn1 , and a constant
d R. A quantitative problem of the form

max z = cT x + d |Ax b, x 0 , (5.1)
33
34 CHAPTER 5. LINEAR PROGRAMMING
or, expressed in terms of a component notation,
max z(x1 , . . . , xn ) = c1 x1 + . . . + cn xn + d (5.2)

a11 x1 + . . . + a1n xn b1 (5.3)
..
.
am1 x1 + . . . + amn xn bm (5.4)
x1 0 (5.5)
..
.
xn 0 , (5.6)
is referred to as a standard maximum problem of linear programming with n real-valued vari-

ables. The different quantities and relations appearing in this definition are called
z(x1 , . . . , xn ) linear objective function, the dependent variable,
x1 , . . . , xn n independent variables,
Ax b m restrictions,
x 0 n non-negativity constraints.
Remark: In an analogous fashion one may also formulate a standard minimum problem of
linear programming, which can be cast into the form

min z = cT x + d |Ax b, x 0 .
In this case, the components of the vector b need to be interpreted as lower limits on certain
capacities.
For given linear objective function z(x1 , . . . , xn ), the set of points x = (x1 , . . . , xn )T satisfying
the condition
z(x1 , . . . , xn ) = C = constant R , (5.7)
for fixed value of C, is referred to as an isoquant of z. Isoquants of linear objective functions

of n = 2 independent variables constitute straight lines, of n = 3 independent variables Euclid-
ian planes, of n = 4 independent variables Euclidian 3-spaces (or hyperplanes), and of n 5
independent variables Euclidian (n 1)-spaces (or hyperplanes).
In the simplest cases of linear programming, the linear objective function z depends on n = 2
variables x1 and x2 only. An illustrative and efficient method of solving problems of this kind will
be looked at in the following section.
5.2. GRAPHICAL SOLUTION METHOD 35
5.2 Graphical method for solving problems with two indepen-
dent variables
The systematic graphical solution method of standard maximum problems of linear programming
with n = 2 independent variables comprises the following steps:
1. Derivation of the linear objective function
z(x1 , x2 ) = c1 x1 + c2 x2 + d
in dependence on the variables x1 and x2 .
2. Identification in the x1 , x2 plane of the feasible region D of z which is determined by the

m restrictions imposed on x1 and x2 . Specifically, D constitutes the domain of z (cf. Ch. 7).
3. Plotting in the x1 , x2 plane of the projection of the isoquant of the linear objective function z
which intersects the origin (0 = x1 = x2 ). When c2 6= 0, this projection is described by the
equation
x2 = (c1 /c2 )x1 .
4. Erecting in the origin of the x1 , x2 plane the direction of optimisation for z which is deter-
mined by the constant z-gradient
z
T x1 c1
(z) = z = .
x2
c2
5. Parallel displacement in the x1 , x2 plane of the projection of the (0, 0)-isoquant of z along
the direction of optimisation (z)T across the feasible region D out to a distance where the
projected isoquant just about touches D.
6. Determination of the optimal solution (x1O , x2O ) as the point resp. set of points of inter-
section between the displaced projection of the (0, 0)-isoquant of z and the far boundary
of D.
7. Computation of the optimal value of the linear objective function zO = z(x1O , x2O ) from
the optimal solution (x1O , x2O ).
8. Specification of potential remaining resources by substitution of the optimal solution

(x1O , x2O ) into the m restrictions.
In general one finds that for a linear objective function z with n = 2 independent variables x1
and x2 , the feasible region D, when non-empty and bounded, constitutes an area in the x1 , x2
plane with straight edges and a certain number of vertices. In these cases, the optimal values of
the linear objective function z are always to be found either at the vertices or on the edges of the
feasible region D. When D is an empty set, then there exists no solution to the corresponding
linear programming problem. When D is unbounded, again there may not exist a solution to the
linear programming problem, but this then depends on the specific circumstances that apply.
Remark: To solve a standard minimum problem of linear programming with n = 2 indepen-
dent variables by means of the graphical method, one needs to parallelly displace in the x1 , x2
plane the projection of the (0, 0)-isoquant of z along the direction of optimisation (z)T until
contact is made with the feasible region D for the first time. The optimal solution is then given by
the point resp. set of points of intersection between the displaced projection of the (0, 0)-isoquant
of z and the near boundary of D.
5.3 Dantzigs simplex algorithm

The main disadvantage of the graphical solution method is its limitation to problems with
only n = 2 independent variables. In actual practice, however, one is often concerned with
linear programming problems that depend on more than two independent variables. To
deal with these more complex problems in a systematic fashion, the US-American mathemati-
cian George Bernard Dantzig (19142005) has devised during the 1940ies an efficient algorithm
which can be programmed on a computer in a fairly straightforward fashion; cf. Dantzig
(1949,1955) [8, 9].
In mathematics, simplex is an alternative name used to refer to a convex polyhedron, i.e., a body of
finite (hyper-)volume in two or more dimensions bounded by linear (hyper-)surfaces which inter-
sect in linear edges and vertices. In general the feasible regions of linear programming problems
constitute such simplexes. Since the optimal solutions for the independent variables of linear
programming problems, when they exist, are always to be found at a vertex or along an edge of
simplex feasible regions, Dantzig developed his so-called simplex algorithm such that it system-
atically scans the edges and vertices of a feasible region to identify the optimal solution (when it
exists) in as few steps as possible.
The starting point be a standard maximum problem of linear programming with n indepen-
dent variables in the form of relations (5.2)(5.6). First, by introducing m non-negative slack
variables s1 , . . . , sm , one transforms the m linear restrictions (inequalities) into an equivalent set
of m linear equations. In this way, potential differences between the left-hand and the right-hand
sides of the m inequalities are represented by the slack variables. In combination with the defining
equation of the linear objective function z, one thus is confronted with a system of 1 + m linear
algebraic equations for the 1 + n + m variables z, x1 , . . . , xn , s1 , . . . , sm , given by
Maximum problem of linear programming in canonical form
z c1 x1 c2 x2 . . . cn xn = d (5.8)
a11 x1 + a12 x2 + . . . + a1n xn + s1 = b1 (5.9)
a21 x1 + a22 x2 + . . . + a2n xn + s2 = b2 (5.10)
..
.
am1 x1 + am2 x2 + . . . + amn xn + sm = bm . (5.11)
As discussed previously in Ch. 3, a system of linear algebraic equations of format (1 + m)

(1 + n + m) is under-determined and so, at most, allows for multiple solutions. The general
5.3. DANTZIGS SIMPLEX ALGORITHM 37
(1 + n + m)-dimensional solution vector
xL = (zL , x1,L , . . . , xn,L , s1,L , . . . , sm,L )T (5.12)
thus contains n variables the values of which can be chosen arbitrarily. It is very important to be
aware of this fact. It implies that, given the linear system is solvable in the first place, one has a
choice amongst different solutions, and so one can pick the solution which proves optimal for the
given problem at hand. Dantzigs simplex algorithm constitues a tool for determining such an
optimal solution in a systematic way.
Let us begin by transferring the coefficients and right-hand sides (RHS) of the under-determined
linear system introduced above into a particular kind of simplex tableau.
Initial simplex tableau
z x1 x2 . . . xn s1 s2 . . . sm RHS
1 c1 c2 . . . cn 0 0 ... 0 d
0 a11 a12 . . . a1n 1 0 ... 0 b1
0 a21 a22 . . . a2n 0 1 ... 0 b2 (5.13)
.. .. .. . . .. .. .. .. . ..
. . . . . . . . .. .
0 am1 am2 . . . amn 0 0 ... 1 bm
In such a simplex tableau one distinguishes so-called basis variables from non-basis variables.
Basis variables are those that contain in their respective columns in the number tableau a (1 + m)-
component canonical unit vector [cf. Eq. (1.10)]; in total the simplex tableau contains 1 + m of
these. Non-basis variables are the remaining ones that do not contain a canonical basis vector
in their respective columns; there exist n of this kind. The complete basis can thus be perceived
as spanning a (1 + m)-dimensional Euclidian space R1+m . Initially, always z and the m slack
variables s1 , . . . , sm constitute the basis variables, while the n independent variables x1 , . . . , xn
classify as non-basis variables [cf. the initial tableau (5.13)]. The corresponding so-called (first)
basis solution has the general appearance
xB1 = (zB1 , x1,B1 , . . . , xn,B1 , s1,B1 , . . . , sm,B1 )T = (d, 0, . . . , 0, b1 , . . . , bm )T ,
since, for simplicity, each of the n arbitrarily specifiable non-basis variables may be assigned the
special value zero. In this respect basis solutions will always be special solutions (as opposed
to general ones) of the under-determined system (5.8)(5.11) the maximum problem of linear
programming in canonical form.
Central aim of the simplex algorithm is to bring as many of the n independent vari-
ables x1 , . . . , xn as possible into the (1 + m)-dimensional basis, at the expense of one of the
m slack variables s1 , . . . , sm , one at a time, in order to construct successively more favourable
special vector-valued solutions to the optimisation problem at hand. Ultimately, the simplex algo-
rithm needs to be viewed as a special variant of Gauian elimination as discussed in Ch. 3, with a
set of systematic instructions concerning allowable equivalence transformations of the underlying
under-determined linear system (5.8)(5.11), resp. the initial simplex tableau (5.13). This set of
systematic algebraic simplex operations can be summarised as follows:
Simplex operations
S1: Does the current simplex tableau show cj 0 for all j {1, . . . , n}? If so, then the
corresponding basis solution is optimal. END. Otherwise goto S2.
S2: Choose a pivot column index j {1, . . . , n} such that cj := min{cj |j

{1, . . . , n}} < 0.
S3: Is there a row index i {1, . . . , m} such that ai j > 0? If not, the objective function z is
unbounded from above. END. Otherwise goto S4.
S4: Choose a pivot row index i such that ai j > 0 and bi /ai j := min{bi /ai j |ai j > 0, i
{1, . . . , m}}. Perform a pivot operation with the pivot element ai j . Goto S1.
When the final simplex tableau has been arrived at, one again assigns the non-basis variables the
value zero. The values of the final basis variables corresponding to the optimal solution of the
given linear programming problem are then to be determined from the final simplex tableau by
backward substitution, beginning at the bottom row. Note that slack variables with positive values
belonging to the basis variables in the optimal solution provide immediate information on existing
remaining capacities in the problem at hand.
Chapter 6
Elementary financial mathematics
In this chapter we want to provide a brief introduction into some basic concepts of financial math-
ematics. As we will try to emphasise, many applications of these concepts (that have immediate
practical relevance) are founded on only two simple and easily accessible mathematical structures:
the so-called arithmetical and geometrical real-valued sequences and their associated finite series.
6.1 Arithmetical and geometrical sequences and series

6.1.1 Arithmetical sequence and series
An arithmetical sequence of n N real numbers an R,
(an )nN ,
is defined by the property that the difference d between neighbouring elements in the sequence be
constant, i.e., for n > 1
an an1 =: d = constant 6= 0 , (6.1)
with an , an1 , d R. Given this recursive formation rule, one may infer the explicit representa-
tion of an arithmetical sequence as
an = a1 + (n 1)d with n N . (6.2)
Note that any arithmetical sequence is uniquely determined by the two free parameters a1 and d,
the starting value of the sequence and the constant difference between neighbours in the sequence,
respectively. Equation (6.2) shows that the elements an in a non-trivial arithmetical sequence
exhibit either linear growth or linear decay with n.
When one calculates for an arithmetical sequence of n + 1 real numbers the arithmetical mean
of the immediate neighbours of any particular element an (with n 2), one finds that
1 1
(an1 + an+1 ) = (a1 + (n 2)d + a1 + nd) = a1 + (n 1)d = an . (6.3)
2 2
39
40 CHAPTER 6. ELEMENTARY FINANCIAL MATHEMATICS
Summation of the first n elements of an arbitrary arithmetical sequence of real numbers leads to
a finite arithmetical series,
n n
X X d
Sn := a1 + a2 + . . . + an = ak = [a1 + (k 1)d] = na1 + (n 1)n . (6.4)
k=1 k=1
2
In the last algebraic step use was made of the Gauian identity1 (cf., e.g., Bosch (2003) [6, p 21])
n1
X 1
k (n 1)n . (6.5)
k=1
2
6.1.2 Geometrical sequence and series

A geometrical sequence of n N real numbers an R,
(an )nN ,
is defined by the property that the quotient q between neighbouring elements in the sequence be
constant, i.e., for n > 1
an
=: q = constant 6= 0 , (6.6)
an1
with an , an1 R and q R\{0, 1}. Given this recursive formation rule, one may infer the
explicit representation of a geometrical sequence as
an = a1 q n1 with n N. (6.7)
Note that any geometrical sequence is uniquely determined by the two free parameters a1 and q,
the starting value of the sequence and the constant quotient between neighbours in the sequence,
respectively. Equation (6.7) shows that the elements an in a non-trivial geometrical sequence
exhibit either exponential growth or exponential decay with n (cf. Sec. 7.1.4).
When one calculates for a geometrical sequence of n + 1 real numbers the geometrical mean of
the immediate neighbours of any particular element an (with n 2), one finds that
p
an1 an+1 = a1 q n2 a1 q n = a1 q n1 = an . (6.8)
Summation of the first n elements of an arbitrary geometrical sequence of real numbers leads to
a finite geometrical series,
n n n1
X X k1
X qn 1
Sn := a1 + a2 + . . . + an = ak = a1 q = a1 q k = a1 . (6.9)
k=1 k=1 k=0
q1
In the last algebraic step use was made of the identity (cf., e.g., Bosch (2003) [6, p 27])
n1
X qn 1
qk for q R\{0, 1} . (6.10)
k=0
q1
n
X
1
Analogously, the modified Gauian identity (2k 1) n2 applies.
k=1
6.2. INTEREST AND COMPOUND INTEREST 41
6.2 Interest and compound interest
Let us consider a first rather simple interest model. Suppose given an initial capital of positive
value K0 > 0 CU paid into a bank account at some initial instant, and a time interval consisting
of n N periods of equal lengths. At the end of each period, the money in this bank account
shall earn a service fee corresponding to an interest rate of p > 0 percent. Introducing the
dimensionless interest factor2
p
q := 1 + >1, (6.11)
100
one finds that by the end of the first interest period a total capital of value (in CU)
p p
K1 = K0 + K0 = K0 1 + = K0 q
100 100
will have accumulated. When the entire time interval of n interest periods has ended, a final
capital worth of (in CU)
recursively: Kn = Kn1 q , nN, (6.12)
will have accumulated, where Kn1 denotes the capital (in CU) accumulated by the end of n 1
interest periods. This recursive representation of the growth of the initial capital K0 due to a total
of n interest payments and the effect of compound interest makes explicit the direct link with the
mathematical structure of a geometrical sequence of real numbers (6.6).
It is a straightforward exercise to show that in this simple interest model the final capital Kn is
related to the initial capital K0 by
explicitly: Kn = K0 q n , nN. (6.13)
Note that this equation links the four non-negative quantities Kn , K0 , q and n to one another.
Hence, knowing the values of three of these quantities, one may solve Eq. (6.13) to obtain the
value of the fourth. For example, solving Eq. (6.13) for K0 yields
Kn
K0 = =: B0 . (6.14)
qn
In this particular variant, K0 is referred to as the present value B0 of the final capital Kn ; this is
obtained from Kn by an n-fold division with the interest factor q.
Further possibilities of re-arranging Eq. (6.13) are:
(i) Solving for the interest factor q: r
n Kn
q= , (6.15)
K0
(ii) Solving for the contract period n:
ln (Kn /K0 )
n= . (6.16)
ln(q)
2
Inverting this defining relation for q leads to p = 100 (q 1).
From now on, n N shall denote the number of full years that have passed in a specific interest
model.
Now we turn to discuss a second, more refined interest model. Let us suppose that an initial
capital K0 > 0 CU earns interest during one full year m N times at the mth part of a nominal
annual interest rate pnom > 0. At the end of the first out of m periods of equal length 1/m, the
initial capital K0 will thus have increased to an amount
pnom pnom
K1/m = K0 + K0 = K0 1 + .
m 100 m 100
By the end of the kth (k m) out of m periods the account balance will have become
pnom k
Kk/m = K0 1 + ;
m 100
pnom
the interest factor 1 + will then have been applied k times to K0 . At the end of the full
m 100
year, K0 in this interest model will have increased to
pnom m
K1 = Km/m = K0 1 + , mN.
m 100
This relation defines an effective interest factor
pnom m
qeff := 1 + , (6.17)
m 100
with associated effective annual interest rate
h pnom m i
peff = 100 1 + 1 , mN, (6.18)
m 100
peff
obtained from re-arranging qeff = 1 + .
100
When, ultimately, n N full years will have passed in the second interest model, the initial
capital K0 will have been transformed into a final capital of value
pnom nm n
Kn = K 0 1 + = K0 qeff , n, m N . (6.19)
m 100
The present value B0 of Kn is thus given by
Kn
B0 = n
= K0 . (6.20)
qeff
Finally, as a third interest model relevant to applications in Finance, we turn to consider the con-
cept of installment savings. For simplicity, let us restrict our discussion to the case when n N
equal installments of constant value E > 0 CU are paid into an account that earns p > 0 percent
annual interest (i.e., q > 1) at the beginning of each of n full years. The initial account balance
6.3. REDEMPTION PAYMENTS IN CONSTANT ANNUITIES 43
be K0 = 0 CU. At the end of a first full year in this interest model, the account balance will have
increased to
p p
K1 = E + E =E 1+ = Eq .
100 100
At the end of two full years one finds, substituting for K1 ,
K2 = (K1 + E)q = (Eq + E)q = E(q 2 + q) = Eq(q + 1) .
At the end of n full years we have, recursively substituting for Kn1 , Kn2 , etc.,
n1
X
n 2 n1
Kn = (Kn1 + E)q = = E(q + . . . + q + q) = Eq(q + . . . + q + 1) = Eq qk .
k=0
Using the identity (6.10), since presently q > 1, the account balance at the end of n full years can
be reduced to the expression
qn 1
Kn = Eq , q R>1 , nN. (6.21)
q1
The present value B0 associated with Kn is obtained by n-fold division of Kn with the interest
factor q:
Eq. 6.21
Kn z}|{ E(q n 1)
B0 := n = = n1 . (6.22)
q q (q 1)
This gives the value of an initial capital B0 which will grow to the same final value Kn after n
annual interest periods with constant interest factor q > 1.
Lastly, re-arranging Eq. (6.21) to solve for the contract period n yields.
ln [1 + (q 1)(Kn /Eq)]
n= . (6.23)
ln(q)
6.3 Redemption payments in constant annuities

The starting point of the next discussion be a mortgage loan of amount R0 > 0 CU that an
economic agent borrowed from a bank at the obligation of annual service payments of p > 0
percent (i.e., q > 1) on the remaining debt. We suppose that the contract between the agent and
the bank fixes the following conditions:
(i) the first redemption payment T1 amount to t > 0 percent of the mortgage R0 ,
(ii) the remaining debt shall be paid back to the bank in constant annuities of value A > 0 CU
at the end of each full year that has passed.
The annuity A is defined as the sum of the variable nth interest payment Zn > 0 CU and the
variable nth redemption payment Tn > 0 CU. In the present model we impose on the annuity
the condition that it be constant across full years,
!
A = Zn + Tn = constant . (6.24)
For n = 1, for example, we thus obtain

p t p+t t !
A = Z1 + T1 = R0 + R0 = R0 = R0 (q 1) + = constant . (6.25)
100 100 100 100
For the first full year of a running mortgage contract, the interest payment, the redemption payment,
and, following the payment of a first annuity, the remaining debt take the values
p
Z1 = R0 = R0 (q 1)
100
T1 = A Z1
substitute for Z1
z}|{ p
R1 = R0 + Z1 A = R0 + R0 A = R0 q A .
100
By the end of a second full year, these become
Z2 = R1 (q 1)
T2 = A Z2
substitute for Z2 substitute for R1
z}|{ z}|{
R2 = R1 + Z2 A = R1 q A = R0 q 2 A(q + 1) .
At this stage, it has become clear according to which patterns the different quantities involved in
the redemption payment model need to be formed. The interest payment for the nth full year in a
mortgage contract of constant anuities amounts to (recursively)
Zn = Rn1 (q 1) , nN, (6.26)
where Rn1 denotes the remaining debt at the end of the previous full year. The redemption
payment for full year n is then given by (recursively)
Tn = A Zn , nN. (6.27)
The remaining debt at the end of the nth full year then is (in CU)
recursively: Rn = Rn1 + Zn A = Rn1 q A , nN. (6.28)
By successive backward substitution for Rn1 , Rn2 , etc., Rn can be re-expressed as

n1
X
n n1 n
Rn = R0 q A(q + . . . + q + 1) = R0 q A qk .
k=0
Now employing the identity (6.10), we finally obtain (since q > 1)
qn 1
explicitly: Rn = R0 q n A , nN. (6.29)
q1
All the formulae we have now derived for computing the values of the quantities {n, Zn , Tn , Rn }
form the basis of a formal redemption payment plan, given by
6.3. REDEMPTION PAYMENTS IN CONSTANT ANNUITIES 45
n Zn [CU] Tn [CU] Rn [CU]
0 R0
1 Z1 T1 R1 ,
2 Z2 T2 R2
.. .. .. ..
. . . .
a standard scheme that banks must make available to their mortgage customers for the purpose of
financial orientation.
Remark: For known values of the free parameters R0 > 0 CU, q > 1 and A > 0 CU, the simple
recursive formulae (6.26), (6.27) and (6.28) can be used to implement a redemption payment plan
in a modern spreadsheet programme such as EXCEL or OpenOffice.
We emphasise the following observation concerning Eq. (6.29): since the constant annuity A con-
tains implicitly a factor (q 1) [cf. Eq. (6.25)], the two competing terms in this relation each grow
exponentially with n. For the redemption payments to eventually terminate, it is thus essential to
fix the free parameter t (for known p > 0 q > 1) in such a way that the second term on the
right-hand side of Eq. (6.29) is given the possibility to catch up with the first as n progresses (the
latter of which has a head start of R0 > 0 CU at n = 0). The necessary condition following from
!
the requirement that Rn Rn1 is thus t > 0.
Equation (6.29) links the five non-negative quantities Rn , R0 , q, n and A to one another. Given
one knows the values of four of these, one can solve for the fifth. For example:
(i) Calculation of the contract period n of a mortgage contract, knowing the mortgage R0 , the
!
interest factor q and the annuity A. Solving the condition Rn = 0 imposed on Rn for n
yields (after a few algebraic steps)

ln 1 + pt
n= ; (6.30)
ln(q)
the contract period is thus independent of the value of the mortgage loan, R0 .
(ii) Evaluation of the annuity A, knowing the contract period n, the mortgage loan R0 , and the
!
interest factor q. Solving the condition Rn = 0 imposed on Rn for A immediately yields
q n (q 1)
A= R0 . (6.31)
qn 1
Now equating the two expressions (6.31) and (6.25) for the annuity A, one finds in addition
that
t q1
= n . (6.32)
100 q 1
6.4 Pension calculations
Quantitative models for pension calculations assume given an initial capital K0 > 0 CU that
was paid into a bank account at a particular moment in time. The issue is to monitor the subse-
quent evolution in discrete time n of the account balance Kn (in CU), which is subjected to two
competing influences: on the one-hand side, the bank account earns interest at an annual interest
rate of p > 0 percent (i.e., q > 1), on the other, it is supposed that throughout one full year a total
of m N pension payments of the constant amount a are made from this bank account, always
at the beginning of each of m intervals of equal duration per year.
Let us begin by evaluating the amount of interest earned per year by the bank account. An impor-
tant point in this respect is the fact that throughout one full year there is a total of m deductions of
value a from the bank account, i.e., in general the account balance does not stay constant through-
out that year but rather decreases in discrete steps. For this reason, the account is credited by the
bank with interest only at the mth part of p > 0 percent for each interval (out of the total of m)
that has passed, with no compound interest effect. Hence, at the end of the first out of m intervals
per year the bank account has earned interest worth of (in CU)
p (q 1)
Z1/m = (K0 a) = (K0 a) .
m 100 m
The interest earned for the kth interval (out of m; k m) is then given by
(q 1)
Zk/m = (K0 ka) .
m
Summation over the contributions of each of the m intervals to the interest earned then yields for
the entire interest earned during the first full year (in CU)
m m
" m
#
X X (q 1) (q 1) X
Z1 = Zk/m = (K0 ka) = mK0 a k .
k=1 k=1
m m k=1
By means of substitution from the identity (6.5), this result can be recast into the equivalent form

1
Z1 = K0 (m + 1)a (q 1) . (6.33)
2
Note that this quantity decreases linearly with the number of deductions m made per year resp. with
the pension payment amount a.
One now finds that the account balance at the end of the first full year that has passed is given by
Eq. (6.33)
z}|{ 1
K1 = K0 ma + Z1 = K0 q m + (m + 1)(q 1) a .
2
At the end of a second full year of the pension payment contract the interest earned is

1
Z2 = K1 (m + 1)a (q 1) ,
2
6.4. PENSION CALCULATIONS 47
while the account balance amounts to
substitute for K1 and Z2
z}|{ 2 1
K2 = K1 ma + Z2 = K0 q m + (m + 1)(q 1) a(q + 1) .
2
At this stage, certain fairly simple patterns for the interest earned during full year n, and the
account balance after n full years, reveal themselves. For Zn we have

1
Zn = Rn1 (m + 1)a (q 1) , (6.34)
2
and for Kn one obtains
substitute for Kn1 and Zn Xn1

z}|{ n 1
Kn = Kn1 ma + Zn = K0 q m + (m + 1)(q 1) a qk .
2 k=0
The latter result can be re-expressed upon substitution from the identity (6.10). Thus, Kn can
finally be given by
n
n 1 q 1
explicitly: Kn = K0 q m + (m + 1)(q 1) a , n, m N . (6.35)
2 q1
In a fashion practically identical to our discussion of the redemption payment model in Sec. 6.3,
the two competing terms on the right-hand side of Eq. (6.35) likewise exhibit exponential growth
with the number n of full years passed. Specifically, it depends on the values of the parameters
K0 > 0 CU, q > 1, a > 0 CU, as well as m 1, whether the second term eventually manages to
catch up with the first as n progresses (the latter of which, in this model, is given a head start of
value K0 > 0 CU at n = 0).
We remark that Eq. (6.35), again, may be algebraically re-arranged at ones convenience (as long
as division by zero is avoided). For example:
(i) The duration n (in full years) of a particular pension contract is obtained from solving the
!
condition Kn = 0 accordingly. Given that [. . .]a K0 (q 1) > 0, one thus finds3

[...]a
ln [...]aK0 (q1)
n= . (6.36)
ln(q)
(ii) The present value B0 of a pension scheme results from the following consideration: for
fixed interest factor q > 1, which initial capital K0 > 0 CU must be paid into a bank account
such that for a duration of n full years one can receive payments of constant amount a at the
beginning of each of m intervals (of equal length) per year? The value of B0 = K0 is again
!
obtained from imposing on Eq. (6.35) the condition Kn = 0 and solving for K0 . This yields

1 qn 1
B0 = K0 = m + (m + 1)(q 1) a n . (6.37)
2 q (q 1)
3
1

To avoid notational overload, the brackets [. . .] here represent the term m + 2 (m + 1)(q 1) .
(iii) The idea of so-called everlasting pension payments of amount aever > 0 CU is based on the
strategy to consume only the annual interest earned by an initial capital K0 > 0 CU residing
in a bank account with interest factor q > 1. Imposing now on Eq. (6.35) the condition
!
Kn = K0 to hold for all values of n, and then solving for a, yields the result
q1
aever = 1 K0 ; (6.38)
m + (m + 1)(q + 1)
2
Note that, naturally, aever is directly proportional to the initial capital K0 !
6.5 Linear and declining-balance depreciation methods

Attempts at the quantitative description of the process of declining material value of industrial
goods, properties or other assets, of initial value K0 > 0 CU, are referred to as depreciation. In-
ternational tax laws generally provide investors with a choice between two particular mathematical
methods of calculating depreciation. We will discuss these options in turn.
6.5.1 Linear depreciation method

When the initial value K0 > 0 CU is supposed to decline to 0 CU in the space of N full years by
equal annual amounts, the remaining value Rn (in CU) at the end of n full years is described by

K0
Rn = K0 n , n = 1, . . . , N . (6.39)
N
Note for the difference of remaining values for years adjacent one obtains Rn Rn1 =
that
K0
=: d < 0. The underlying mathematical structure of the straight line depreciation
N
method is thus an arithmetical sequence of real numbers, with constant negative difference d
between neighbouring elements (cf. Sec. 6.1.1).
6.5.2 Declining-balance depreciation method

The foundation of the second depreciation method to be described here, for an industrial good of
initial value K0 > 0 CU, is the idea that per year the value declines by a certain percentage rate
p > 0 of the value of the good during the previous year. Introducing a dimensionless depreciation
factor by
p
q := 1 <1, (6.40)
100
the remaining value Rn (in CU) after n full years amounts to
recursively: Rn = Rn1 q , R0 K0 , nN. (6.41)
The underlying mathematical structure of the declining balance depreciation method is thus a
geometrical sequence of real numbers, with constant ratio 0 < q < 1 between neighbouring
6.6. SUMMARISING FORMULA 49
elements (cf. Sec. 6.1.2). With increasing n the values of these elements become ever smaller. By
means of successive backward substitution expression (6.41) can be transformed to
explicitly: Rn = K0 q n , 0<q<1, nN. (6.42)
From Eq. (6.42), one may derive results concerning the following questions of a quantitative na-
ture:
(i) Suppose given a depreciation factor q and a projected remaining value Rn for some industrial
good. After which depreciation period n will this value be attained? One finds
ln (Rn /K0 )
n= . (6.43)
ln(q)
(ii) Knowing a projected depreciation period n and corresponding remaining value Rn , at which
percentage rate p > 0 must the depreciation method be operated? This yields
r r !
R n Rn
q= n p = 100 1 n . (6.44)
K0 K0
6.6 Summarising formula

To conclude this chapter, let us summarise the results on elementary financial mathematics that
we derived along the way. Remarkably, these can be condensed in a single formula which contains
the different concepts discussed as special cases. This formula, in which n represents the number
of full years that have passed, is given by (cf. ZehMarschke (2010) [26]):
qn 1
K n = K0 q n + R , q R>0 \{1} , nN. (6.45)
q1
The different special cases contained therein are:

(i) Compound interest for an initial capital K0 > 0 CU: with R = 0 and q > 1, Eq. (6.45)
reduces to Eq. (6.13).
(ii) Installment savings with constant installments E > 0 CU: with K0 = 0 CU, q > 1 and
R = Eq, Eq. (6.45) reduces to Eq. (6.19).
(iii) Redemption payments in constant annuities: with K0 = R0 < 0 CU, q > 1 and
R = A > 0 CU, Eq. (6.45) reduces to the negative (!) of Eq. (6.29). In this dual formulation,
remaining debt Kn = Rn is (meaningfully) expressed as a negative account balance.
(iv) Pension
payments of constant
amount a > 0 CU: with q > 1 and R =
1
m + (m + 1)(q 1) a, Eq. (6.45) transforms to Eq. (6.35).
2
(v) Declining balance depreciation of an asset of initial value K0 > 0 CU: with R = 0 and
0 < q < 1, Eq. (6.45) converts to Eq. (6.42) for the remaining value Kn = Rn .
Chapter 7
Differential calculus of real-valued

functions of one real variable
In Chs. 1 to 5 of these lecture notes, we confined our considerations to functional relationships

between INPUT quantities and OUTPUT quantities of a linear nature. In this chapter now, we
turn to discuss characteristic properties of truly non-linear functional relationships between one
INPUT quantity and one OUTPUT quantity.
7.1 Real-valued functions of one real variable

Let us begin by introducing the concept of a real-valued function of one real variable. This
constitutes a special kind of a mapping1 that needs to satisfy the following simple but very strict
rule:
a mapping f that assigns to every element x from a subset D of the real numbers R (i.e.,
D R) one and only one element y from a second subset W of the real numbers R (i.e., W R).
Def.: A unique mapping f of a subset D R of the real numbers onto a subset W R of the
real numbers,
f: D W , x 7 y = f (x) (7.1)
is referred to as a real-valued function of one real variable.
We now fix some terminology concerning the concept of a real-valued function of one real variable:
D: domain of f ,
W : target space of f ,
x D: independent variable of f , also referred to as the argument of f ,
y W : dependent variable of f ,
1
Cf. our introduction in Ch. 2 of matrices as a particular class of mathematical objects.
51
52 CHAPTER 7. DIFFERENTIAL CALCULUS OF REAL-VALUED FUNCTIONS
f (x): mapping prescription,
graph of f : the set of pairs of values G = {(x, f (x))|x D} R2 .
For later analysis of the mathematical properties of real-valued functions of one real variable, we
need to address a few more technical issues.
Def.: Given a mapping f that is one-to-one and onto, with domain D(f ) R and target space
W (f ) R, not only is every x D(f ) assigned to one and only one y W (f ), but also every
y W (f ) is assigned to one and only one x D(f ). In this case, there exists an associated
mapping f 1 , with D(f 1 ) = W (f ) and W (f 1 ) = D(f ), which is referred to as the inverse
function of f .
Def.: A real-valued function f of one real variable x is continuous at some value x D(f ) when
for x R>0 the condition
lim f (x x) = lim f (x + x) = f (x) (7.2)

x0 x0
obtains, i.e., when at x the left and right limits of the function f coincide and are equal to the value
f (x). A real-valued function f as such is continuous when f is continuous for all x D(f ).
Def.: When a real-valued function f of one real variable x satisfies the condition
f (a) < f (b) for all a, b D(f ) with a < b , (7.3)
then f is called strictly monotonously increasing. When, however, f satisfies the condition
f (a) > f (b) for all a, b D(f ) with a < b , (7.4)
then f is called strictly monotonously decreasing.

Note, in particular, that real-valued functions of one real variable that are strictly monotonous and
continous are always one-to-one and onto and, therefore, are invertible.
In the following, we briefly review five elementary classes of real-valued functions of one real vari-
able that find frequent application in the modelling of quantitative problems in economic theory.
7.1.1 Polynomials of degree n

Polynomials of degree n are real-valued functions of one real variable of the form
y = f (x) = an xn + an1 xn1 + . . . + ai xi + . . . + a2 x2 + a1 x + a0

(7.5)
with ai R, i = 1, . . . , n, n N, an 6= 0 .
Their domain comprises the entire set of real numbers, i.e., D(f ) = R. The extent of their target
space depends specifically on the values of the real constant coefficients ai R. Functions in this
class possess a maximum of n real roots.
7.1. REAL-VALUED FUNCTIONS 53
7.1.2 Rational functions
Rational functions are constructed by forming the ratio of two polynomials of degrees m resp. n,
i.e.,
pm (x) am xm + . . . + a1 x + a0
y = f (x) = =
qn (x) bn xn + . . . + b1 x + b0 (7.6)
with ai , bj R, i = 1, . . . , m, j = 1, . . . , n, m, n N, am , bn 6= 0 .
Their domain is given by D(f ) = R\{x|qn (x) = 0}. When for the degrees m and n of the
polynomials pm (x) and qn (x) we have
(i) m < n, then f is referred to as a proper rational function, or
(ii) m n, then f is referred to as an improper rational function.
In the latter case, application of polynomial division leads to a separation of f into a purely poly-
nomial part and a proper rational part. The roots of f always correspond to those roots of the
numerator polynomial pm (x) for which simultaneously qn (x) 6= 0 applies. The roots of the de-
nominator polynomial qn (x) constitute poles of f . Proper rational functions always tend for very
small (i.e., x ) and for very large (i.e., x +) values of their argument to zero.
7.1.3 Power-law functions

Power-law functions exhibit the specific structure given by
y = f (x) = x with R . (7.7)
We here confine ourselves to cases with domains D(f ) = R>0 , such that for the coresponding
target spaces we have W (f ) = R>0 . Under these conditions, power-law functions are strictly
monotonously increasingwhen > 0, and strictly monotonously decreasing when < 0. Hence,
they are inverted by y = x = x1/ . There do not exist any roots under the conditions stated here.
7.1.4 Exponential functions

Exponential functions have the general form
y = f (x) = ax with a R>0 \{1} . (7.8)
Their domain is D(f ) = R, while their target space is W (f ) = R>0 . They exhibit strict
monotonous increase for a > 1, and strict monotonous decrease for 0 < a < 1. Hence, they
too are invertible. Their y-intercept is generally located at y = 1. For a > 1, exponential functions
are also known as growth functions.
Special case: When the constant (!) base number is chosen to be a = e, where e denotes the
irrational Eulers number (according to the Swiss mathematician Leonhard Euler, 17071783)
defined by the infinite series

X 1 1 1 1 1
e := = + + + + ... ,
k=0
k! 0! 1! 2! 3!
one obtains the natural exponential function
y = f (x) = ex =: exp(x) . (7.9)
In analogy to the definition of e, the relation

X xk x0 x1 x2 x3
ex = exp(x) = = + + + + ...
k=0
k! 0! 1! 2! 3!
applies.
7.1.5 Logarithmic functions

Logarithmic functions, denoted by
y = f (x) = loga (x) with a R>0 \{1} , (7.10)
are defined as inverse functions of the strictly monotonous exponential functions y = f (x) = ax
and vice versa. Correspondingly, D(f ) = R>0 and W (f ) = R apply. Strictly monotonously
increasing behaviour is given when a > 1, strictly monotonously decreasing behaviour when
0 < a < 1. In general, the x-intercept is located at x = 1.
Special case: The natural logarithmic function (lat.: logarithmus naturalis) obtains when the
constant basis number is set to a = e. This yields
y = f (x) = loge (x) := ln(x) . (7.11)
7.1.6 Concatenations of real-valued functions

Real-valued functions from all five categories considered in the previous sections may be combined
arbitrarily (respecting relevant computational rules), either via the four fundamental arithmetical
operations, or via concatenations.
Theorem: Let real-valued functions f and g be continuous on domains D(f ) resp. D(g). Then
the combined real-valued functions
1. sum/difference f g, where (f g)(x) := f (x) g(x) with D(f ) D(g),
2. product f g, where (f g)(x) := f (x)g(x) with D(f ) D(g),

f f f (x)
3. quotient , where (x) := with g(x) 6= 0 and D(f ) D(g)\{x|g(x) = 0},
g g g(x)
4. concatenation f g, where (f g)(x) := f (g(x)) mit {x D(g)|g(x) D(f )},
are also continuous on the respective domains.

7.2. DERIVATION OF DIFFERENTIABLE REAL-VALUED FUNCTIONS 55
7.2 Derivation of differentiable real-valued functions
The central theme of this chapter is the mathematical description of the local variability of con-
tinuous real-valued function of one real variable, f : D R W R. To this end, let
us consider the effect on f of a small change of its argument x. Supposing we affect a change
x x + x, with x R, what are the resultant consequences for f ? We immediately find that
y y + y = f (x + x), with y R, obtains. Hence, a prescribed change of the argument x
by a (small) value x triggers in f a change by the amount y = f (x + x) f (x). It is of
general quantitative interest to compare the sizes of these two changes. This is accomplished by
forming the respective difference quotient
y f (x + x) f (x)
= .
x x
It is then natural, for given f , to investigate the limit behaviour of this difference quotient as the
change x of the argument of f is made successively smaller.
Def.: A continuous real-valued function f of one real variable is called differentiable at
x D(f ), when for arbitrary x R the limit
y f (x + x) f (x)
f (x) := lim = lim (7.12)
x0 x x0 x
exists and is unique. When f is differentiable for all x D(f ), then f as such is referred to as
being differentiable.
The existence of this limit in a point (x, f (x)) for a real-valued function f requires that the latter
exhibits neither jumps nor kinks, i.e., that at (x, f (x)) the function is sufficiently smooth.
The quantity f (x) is referred to as the first derivative of the (differentiable) function f at posi-
tion x. It provides a quantitative measure for the local rate of change of the function f in the
point (x, f (x)). In general one interprets the first derivative f (x) as follows: an increase of the
argument x of a differentiable real-valued function f by 1 (one) unit leads to a change in the value
of f by approximately f (x) 1 units.
Alternative notation for the first derivative of f :
df (x)
f (x) .
dx
The differential calculus was developed in parallel with the integral calculus (see Ch. 7) during the
second half of the 17th Century, independent of one another by the English physicist, mathemat-
iccian, astronomer and philosopher Sir Isaac Newton (16431727) and the German philosopher,
mathematician and physicist Gottfried Wilhelm Leibniz (16461716).
Via the first derivative of a differentiable function f at an argument x0 D(f ), i.e., f (x0 ), one
defines the so-called linearisation of f in a neighbourhood of x0 . The equation describing the
associated tangent to f in the point (x0 , f (x0 )) is given by
y = f (x0 ) + f (x0 )(x x0 ) . (7.13)

GDC: Local values f (x0 ) of first derivatives can be computed for given function f in mode CALC
using the interactive routine dy/dx.
The following rules of differentiation apply for the five families of elementary real-valued functions
discussed in Sec. 7.1, as well as concatenations thereof:
Rules of differentiation
1. (c) = 0 for c = constant R (constants)
2. (x) = 1 (linear function)
3. (xn ) = nxn1 for n N (natural power-law functions)
4. (x ) = x1 for R and x R>0 (general power-law functions)
5. (ax ) = ln(a)ax for a R>0 \{1} (exponential functions)
6. (eax ) = aeax for a R (natural exponential functions)

1
7. (loga (x)) = for a R>0 \{1} and x R>0 (logarithmic functions)
x ln(a)
1
8. (ln(x)) = for x R>0 (natural logarithmic function).
x
For differentiable real-valued functions f and g it holds that:
1. (cf (x)) = cf (x) for c = constant R
2. (f (x) g(x)) = f (x) g (x) (summation rule)
3. (f (x)g(x)) = f (x)g(x) + f (x)g (x) (product rule)

f (x) f (x)g(x) f (x)g (x)
4. = (quotient rule)
g(x) (g(x))2
5. ((f g)(x)) = f (g)|g=g(x) g (x) (chain rule)
f (x)
6. (ln(f (x))) = for f (x) > 0 (logarithmic differentiation)
f (x)

1 1
7. (f (x)) = , if f is one-to-one and onto.
f (y) y=f 1 (x)
(differentiation of inverse functions).
The methods of differential calculus just introduced shall now be employed to describe the local
change behaviour of a few simple examples of functions in economic theory, and also to deter-
mine their local extremal values. The following section provides an overview of such frequently
occurring economic functions.
7.3. COMMON FUNCTIONS IN ECONOMIC THEORY 57
7.3 Common functions in economic theory
1. total cost function K(x) 0 (dim: CU)
argument: level of physical output x 0 (dim: units)
2. marginal cost function K (x) > 0 (dim: CU/unit)

argument: level of physical output x 0 (dim: units)
3. average cost function K(x)/x > 0 (dim: CU/unit)

argument: level of physical output x > 0 (dim: units)
4. unit price function p(x) 0 (dim: CU/unit)

5. total revenue function E(x) := xp(x) 0 (dim: CU)

6. marginal revenue function E (x) = xp (x) + p(x) (dim: CU/unit)

7. profit function G(x) := E(x) K(x) (dim: CU)

8. marginal profit function G (x) := E (x) K (x) = xp (x) + p(x) K (x) (dim:
CU/unit)
9. utility function U(x) (dim: case dependent)

argument: material wealth, opportunity, action x (dim: case dependent)
The fundamental concept of a utility function as a means to capture in quantitative terms
the psychological value (happiness) assigned by an economic agent to a certain amount
of money, or to owning a specific good, was introduced to economic theory in 1738 by
the Swiss mathematician and physicist Daniel Bernoulli FRS (17001782); cf. Bernoulli
(1738) [4]. The utility function is part of the folklore of the theory, and often taken to be
a piecewise differentiable, right-handedly curved (concave) function, i.e., U (x) < 0, on
the grounds of the assumption of diminishing marginal utility (happiness) with increasing
material wealth.
10. economic efficiency W (x) := E(x)/K(x) 0 (dim: 1)

11. demand function N(p) 0, monotonously decreasing (dim: units)

argument: unit price p, (0 p pmax ) (dim: CU/unit)
12. supply function A(p) 0, monotonously increasing (dim: units)

argument: unit price p, (pmin p) (dim: CU/units).
A particularly prominent example of a real-valued economic function of one real variable
constitutes the psychological value function, devised by the IsraeliUS-American experi-
mental psychologists Daniel Kahneman and Amos Tversky (19371996) in the context of
their Prospect Theory (a pillar of Behavioural Economics), which was later awarded a
Sveriges Riksbank Prize in Economic Sciences in Memory of Alfred Nobel in 2002 (cf. Kahne-
man and Tversky (1979) [18, p 279], and Kahneman (2011) [17, p 282f]). A possible representa-
tion of this function is given by the piecewise description

a log10 (1 + x)
for x R0
v(x) = , (7.14)

2a log10 (1 x) for x R<0
with parameter a R>0 . Overcoming a conceptual problem of Bernoullis utility function, here,
in contrast, the argument x quantifies a change in wealth (or welfare) with respect to some given
reference point (rather than a specific value of wealth itself).
7.4 Curve sketching

Before we turn to discuss applications of differential calculus to simple quantitative problems
in economic theory, we briefly summarize the main steps of curve sketching for a real-valued
function of one real variable, also referred to as analysis of the properties of differentiability of a
real-valued function.
1. domain: D(f ) = {x R|f (x) is regular}
2. symmetries: for all x D(f ), is
(i) f (x) = f (x), i.e., is f even, or

(ii) f (x) = f (x), i.e., is f odd, or
(iii) f (x) 6= f (x) 6= f (x), i.e., f exhibits no symmetries?
!
3. roots: identify all xN D(f ) that satisfy the condition f (x) = 0.
4. local extremal values:
(i) local minima of f exist at all xE D(f ), for which the

!
necessary condition f (x) = 0, and the
!
sufficient condition f (x) > 0 are satisfied simultaneously.
(ii) local maxima of f exist at all xE D(f ), for which the
!
!
sufficient condition f (x) < 0 are satisfied simultaneously.
7.5. ANALYTIC INVESTIGATIONS OF ECONOMIC FUNCTIONS 59
5. points of inflection: find all xW D(f ), for which the
!
!
sufficient condition f (x) 6= 0 are satisfied simultaneously.
6. monotonous behaviour:
(i) f is monotonously increasing for all x D(f ) with f (x) > 0
(ii) f is monotonously decreasing for all x D(f ) with f (x) < 0
7. local curvature:
(i) f behaves left-handedly curved for x D(f ) with f (x) > 0
(ii) f behaves right-handedly curved for x D(f ) with f (x) < 0
8. asymptotic behaviour:
asymptotes to f are constituted by
(i) straight lines y = ax + b with the property limx+ [f (x) ax b] = 0 or
limx [f (x) ax b] = 0
(ii) straight lines x = x0 at poles x0
/ D(f )
9. range: W (f ) = {y R|y = f (x)}.
7.5 Analytic investigations of economic functions

7.5.1 Total cost functions according to Turgot and von Thunen
According to the law of diminishing returns, which was introduced to economic theory by the
French economist and statesman Anne Robert Jacques Turgot (17271781) and also by the Ger-
man economist Johann Heinrich von Thunen (17831850), it is meaningful to model non-negative
total cost functions K(x) (in CU) relating to typical production processes, with argument level
of physical output x 0 units, as a mathematical mapping in terms of a special polynomial of
degree 3 [cf. Eq. (7.5)], which is given by
K(x) = a3 x3 + a2 x2 + a1 x + a0
| {z } |{z}
=Kv (x) =Kf
(7.15)
with a3 , a1 > 0, a2 < 0, a0 0, a22 3a3 a1 < 0 .
The model thus contains a total of four free parameters. It is the outcome of a systematic regres-
sion analysis of agricultural quantitativeempirical data with the aim to describe an inherently
non-linear functional relationship between a few economic variables. As such, the functional
relationship for K(x) expressed in Eq. (7.15) was derived from a practical consideration. It is a
reflection of the following observed features:
(i) for levels of physical output x 0 units, the total costs relating to typical production
processes exhibit strictly monotonously increasing behaviour; thus
(ii) for the total costs there do not exist neither roots nor local extremal values;2 however,
(iii) the total costs display exactly one point of inflection.
The continuous curve for K(x) resulting from these considerations exhibits the characteristic shape
of an inverted capital letter S: beginning at a positive value corresponding to fixed costs, the total
costs first increase degressively up to a point of inflection, whereafter they continue to increase,
but in a progressive fashion.
In broad terms, the functional expression given in Eq. (7.15) to model totals costs in dependence
of the level of physical output is the sum of two contributions, the variable costs Kv (x) and the
fixed costs Kf = a0 , viz.
K(x) = Kv (x) + Kf . (7.16)
In economic theory, it is commonplace to partition total cost functions in the diminishing returns
picture into four phases, the boundaries of which are designated by special values of the level of
physical output of a production process:
phase I (interval 0 units x xW ):

the total costs K(x) possess at a level of physical output xW = a2 /(3a3 ) > 0 units a
point of inflection. For values of x smaller than xW , one obtains K (x) < 0 CU/unit2 , i.e.,
K(x) increases in a degressive fashion. For values of x larger than xW , the opposite applies,
K (x) > 0 CU/unit2 , i.e., K(x) increases in a progressive fashion. The marginal costs,
given by
K (x) = 3a3 x2 + 2a2 x + a1 > 0 CU/unit for all x 0 units , (7.17)
attain a minimum at the same level of physical output, xW = a2 /(3a3 ).
phase II (interval xW < x xg1 ):

the variable average costs
Kv (x)
= a3 x2 + a2 x + a1 , x > 0 units (7.18)
x
become minimal at a level of physical output xg1 = a2 /(2a3 ) > 0 units. At this value
of x, equality of variable average costs and marginal costs applies, i.e.,
Kv (x)
= K (x) , (7.19)
x
2
The last condition in Eq. (7.15) ensures a first derivative of K(x) that does not possess any roots; cf. the case of
!
a quadratic algebraic equation 0 = ax2 + bx + c, with discriminant b2 4ac < 0.
which follows by the quotient rule of differentiation from the necessary condition for an
extremum of the variable average costs,

! Kv (x) (K(x) Kf ) x Kv (x) 1
0= = ,
x x2
and the fact that Kf = 0 CU/unit. Taking care of the equality (7.19), one finds for the tangent
to K(x) in the point (xg1 , K(xg1 )) the equation [cf. Eq. (7.13)]
Kv (xg1 ) Kv (xg1 )
T (x) = K(xg1 )+K (xg1 )(xxg1 ) = Kv (xg1 )+Kf + (xxg1 ) = Kf + x.
xg1 xg1
Its intercept with the K-axis is at Kf .
phase III (interval xg1 < x xg2 ):

The average costs
K(x) a0
= a3 x2 + a2 x + a1 + , x > 0 units (7.20)
x x
attain a minimum at a level of physical output xg2 > 0 units, the defining equation of which
!
is given by 0 CU = 2a3 x3g2 + a2 x2g2 a0 . At this value of x, equality of average costs and
marginal costs obtains, viz.
K(x)
= K (x) , (7.21)
x
which follows by the quotient rule of differentiation from the necessary condition for an
extremum of the average costs,

! K(x) K (x) x K(x) 1
0= = .
x x2
Since a quotient can be zero only when its numerator vanishes (and its denominator remains
non-zero), one finds from re-arranging the numerator expression equated to zero the property
K (x) K (x)
=x =1 for x = xg2 . (7.22)
K(x)/x K(x)
The corresponding extremal value pair (xg2 , K(xg2 )) is referred to in economic theory as the
minimum efficient scale (MES). From a business economics perspective, at a level of phys-
ical output x = xg2 the (compared to our remarks in the Introduction inverted) ratio INPUT
K(x)
over OUTPUT, i.e., , becomes most favourable. By respecting the property (7.21),
x
the equation for the tangent to K(x) in this point [cf. Eq. (7.13)] becomes
K(xg2 ) K(xg2 )
T (x) = K(xg2 ) + K (xg2 )(x xg2 ) = K(xg2 ) + (x xg2 ) = x.
xg2 xg2
Its intercept with the K-axis is thus at 0 CU.
phase IV (half-interval x > xg2 ):
In this phase K (x)/K(x)/x > 1 obtains; the costs associated with the production of one
additional unit of a good, approximately the marginal costs K (x), now exceed the aver-
age costs, K(x)/x. This situation is considered unfavourable from a business economics
perspective.
7.5.2 Profit functions in the diminishing returns picture

In this section, we confine our considerations, for reasons of simplicity, to a market sitution with
only a single supplier of a good in demand. The price policy that this single supplier may thus
inact defines a state of monopoly. Moreover, in addition we want to assume that for the market
situation considered economic equilibrium obtains. This manifests itself in equality of supply
and demand, viz.
x(p) = N(p) , (7.23)
wherein x denotes a non-negative supply function (in units) (which is synonymous with the sup-
pliers level of physical output) and N a non-negative demand function (in units), both of which
are taken to depend on the positive unit price p (in CU/unit) of the good in question. The supply
function, and with it the unit price, can, of course, be prescribed by the monopolistic supplier in
an arbitrary fashion. In a specific quantitative economic model, for instance, the demand func-
tion x(p) (recall that by Eq. (7.23) x(p) = N(p) obtains) could be assumed to be either a linear or
a quadratic function of p. In any case, in order for x(p) to realistically describe an actual demand
unit price relationship, it should be chosen as a strictly monotonously decreasing function, and as
such it is invertible. The non-negative demand function x(p) features two characteristic points,
signified by its intercepts with the x- and the p-axes. The prohibitive price pproh is to be deter-
!
mined from the condition x(pproh ) = 0 units; therefore, it constitutes a root of x(p). The saturation
quantity xsat , on the other hand, is defined by xsat := x(0 CU/unit).
The inverse function associated with the strictly monotonously decreasing non-negative demand
function x(p), the unit price function p(x) (in CU/unit), is likewise strictly monotonously
decreasing. Via p(x), one calculates, in dependence on a known amount x of units sup-
plied/demanded (i.e., sold), the total revenue (in CU) made by a monopolist according to (cf.
Sec. 7.3)
E(x) = xp(x) . (7.24)
Under the assumption that the non-negative total costs K(x) (in CU) underlying the production
process of the good in demand can be modelled according to the diminishing returns picture of
Turgot and von Thunen, the profit function (in CU) of the monopolist in dependence on the level
of physical output takes the form
unit price
z}|{
G(x) = E(x) K(x) = x p(x) a3 x3 + a2 x2 + a1 x + a0 . (7.25)
| {z } | {z }
total revenue total costs
The first two derivatives of G(x) with respect to its argument x are given by

G (x) = E (x) K (x) = xp (x) + p(x) 3a3 x2 + 2a2 x + a1 (7.26)
G (x) = E (x) K (x) = xp (x) + 2p (x) [6a3 x + 2a2 ] . (7.27)
Employing the principles of curve sketching set out in Sec. 7.4, the following characteristic values
of G(x) can thus be identified:
break-even point
xS > 0 units, as the unique solution to the conditions
!
G(x) = 0 CU (necessary condition) (7.28)
and
!
G (x) > 0 CU/unit (sufficient condition) , (7.29)
end of the profitable zone

xG > 0 units, as the unique solution to the conditions
!
G(x) = 0 CU (necessary condition) (7.30)
and
!
G (x) < 0 CU/unit (sufficient condition) , (7.31)
maximum profit
xM > 0 CU, as the unique solution to the conditions
!
G (x) = 0 CU/unit (necessary condition) (7.32)
and
!
G (x) < 0 CU/unit2 (sufficient condition) . (7.33)
At this point, we like to draw the readers attention to a special geometric property of the quan-
titative model for profit that we just have outlined: at maximum profit, the total revenue func-
tion E(x) and the total cost function K(x) always possess parallel tangents. This is due to the
fact that by the necessary condition for an extremum to exist, one finds that
! !
0 CU/unit = G (x) = E (x) K (x) E (x) = K (x) . (7.34)
GDC: Roots and local maxima resp. minima can be easily determined for a given stored function
in mode CALC by employing the interactive routines zero and maximum resp. minimum.
To conclude these considerations, we briefly turn to elucidate the technical term Cournots point,
which frequently arises in quantitative discussions in economic theory; this is named after the
French mathematician and economist AntoineAugustin Cournot (18011877). Cournots point
simply labels the profit-optimal combination of the level of physical output and the associated unit
price, (xM , p(xM )), for the unit price function p(x) of a good in a monopolistic market situation.
Note that for this specific combination of optimal values the AmorosoRobinson formula applies,
which was developed by the Italian mathematician and economist Luigi Amoroso (18861965)
and the British economist Joan Violet Robinson (19031983). This states that
K (xM )
p(xM ) = , (7.35)
1 + p (xM )
with K (xM ) the value of the marginal costs at xM , and p (xM ) the value of the elasticity of the
unit price function at xM (see the following Sec. 7.6). Starting from the defining equation of the
total revenue E(x) = xp(x), the AmorosoRobinson formula is derived by evaluating the first
derivative of E(x) at xM , so
Sec. 7.6
p (xM ) z}|{
E (xM ) = p(xM ) + xM p (xM ) = p(xM ) 1 + xM = = p(xM ) [1 + p (xM )] ,
p(xM )
and then re-arranging to solve for p(xM ), using the fact that E (xM ) = K (xM ).
Remark: In a market situation where perfect competition applies, one assumes that the unit
price function has settled to a constant value p(x) = p = constant > 0 CU/unit (and, hence,
p (x) = 0 CU/unit2 obtains).
7.5.3 Extremal values of rational economic functions

Now we want to address the determination of extremal values of economic functions that constitute
ratios in the sense of the construction
OUTPUT
,
INPUT
a topic raised in the Introduction.
Let us consider two examples for determining local maxima of ratios of this kind.
(i) We begin with the average profit in dependence on the level of physical output x 0 units,
G(x)
. (7.36)
x
!
The conditions that determine a local maximum are [G(x)/x] = 0 CU/unit2 and
!
[G(x)/x] < 0 CU/unit3 . Respecting the quotient rule of differentiation (cf. Sec. 7.2),
the first condition yields
G (x)x G(x)
2
= 0 GE/ME2 . (7.37)
x
Since a quotient can only be zero when its numerator vanishes while its denominator remains
non-zero, it immediately follows that
G (x)
G (x)x G(x) = 0 CU x =1. (7.38)
G(x)
7.6. ELASTICITIES 65
The task at hand now is to find a (unique) value of the level of physical output x which
satisfies this last condition, and for which the second derivative of the average profit becomes
negative.
(ii) To compare the performance of two companies over a given period of time in a meaningful
way, it is recommended to adhere only to measures that are dimensionless ratios, and so
independent of scale. An example of such a dimensionless ratio is the measure referred to
as economic efficiency,
E(x)
W (x) = , (7.39)
K(x)
which expresses the total revenue (in CU) of a company for a given period as a multiple of
the total costs (in CU) it had to endure during this period, both as functions of the level of
physical output. In analogy to our discussion in (i), the conditions for the existence of a
! !
local maximum amount to [E(x)/K(x)] = 0 1/unit and [E(x)/K(x)] < 0 1/unit2 .
By the quotient rule of differentiation (see Sec. 7.2), the first condition leads to
E (x)K(x) E(x)K (x)
= 0 1/unit , (7.40)
K 2 (x)
i.e., for K(x) > 0 CU,
E (x)K(x) E(x)K (x) = 0 CU2 /unit . (7.41)
By re-arranging and multiplication with x > 0 unit, this can be cast into the particular form
E (x) K (x)
x =x . (7.42)
E(x) K(x)
The reason for this special kind of representation of the necessary condition for a local
maximum to exist [and also for Eq. (7.38)] will be clarified in the subsequent section. Again,
a value of the level of physical output which satisfies Eq. (7.42) must in addition lead to
a negative second derivative of the economic efficiency in order to satisfy the sufficient
condition for a local maximum to exist.
7.6 Elasticities
Finally, we pick up once more the discussion on quantifying the local variability of differentiable
real-valued functions of one real variable, f : D R W R, though from a slightly different
perspective. For reasons to be elucidated shortly, we confine ourselves to considerations of regimes
of f with positive values of the argument x and also positive values y = f (x) > 0 of the function
itself.
As before in Sec. 7.2, we want to assume a small change of the value of the argument x and
evaluate its resultant effect on the value y = f (x). This yields
xR yR
x x + x = y = f (x) y + y = f (x + x) . (7.43)
We remark in passing that relative changes of non-negative quantities are defined by the quotient
new value old value
old value
under the prerequisite that old value > 0 applies. It follows from this specific construction that
the minimum value a relative change can possibly attain amounts to 1 (corresponding to a
decrease of the old value by 100%).
Related to this consideration we identify the following terms:
a prescribed absolute change of the independent variable x: x ,
the resultant absolute change of the function f : y = f (x + x) f (x) ,

x
the associated relative change of the independent variable x: ,
x
y f (x + x) f (x)
the associated resultant relative change of the function f : = .
y f (x)
x
Now let us compare the order-of-magnitudes of the two relative changes just envisaged, and
x
y
. This is realised by considering the value of their quotient, resultant relative change of f
y
divided by the prescribed relative change of x:
y f (x + x) f (x)
y f (x)
= .
x x
x x
Since we assumed f to be differentiable, it is possible to investigate the behaviour of this quotient
x
of relative changes in the limit of increasingly smaller prescribed relative changes 0
x
x 0 near some x > 0. One thus defines:
Def.: For a differentiable real-valued function f of one real variable x, the dimensionless (i.e.,
units-independent) quantity
y f (x + x) f (x)
y f (x) f (x)
f (x) := lim = lim =x (7.44)
x0 x x0 x f (x)
x x
is referred to as the elasticity of the function f at position x.

The elasticity of f quantifies its resultant relative change in response to a prescribed infinitesimally
small relative change of its argument x, starting from some positive initial value x > 0. As such it
constitutes a measure for the relative local rate of change of a function f in a point (x, f (x)). In
economic theory, in particular, one adheres to the following interpretation of the elasticity f (x):
when the postive argument x of some positive differentiable real-valued function f is increased by
1 %, then in consequence f will change approximately by f (x) 1 %.
In the scientific literature one often finds the elasticity of a positive differentiable function f of a
positive argument x expressed in terms of logarithmic differentiation. That is,
d ln[f (x)]
f (x) := for x > 0 and f (x) > 0 ,
d ln(x)
since by the chain rule of differentiation it holds that
df (x) df (x)
d ln[f (x)] f (x) f (x)
= = x dx = x .
d ln(x) dx f (x) f (x)
x
The logarithmic representation of the elasticity of a differentiable function f immediately explains
why, at the beginning, we confined our considerations to positive differentiable functions of posi-
tive arguments only.3 A brief look at the list of standard economic functions provided in Sec. 7.3
reveals that most of these (though not all) are positive functions of non-negative arguments.
For the elementary classes of real-valued functions of one real variable discussed in Sec. 7.1 one
finds:
Standard elasticities
1. f (x) = xn for n N and x R>0 f (x) = n (natural power-law functions)
2. f (x) = x for R and x R>0 f (x) = (general power-law functions)
3. f (x) = ax for a R>0 \{1} and x R>0 f (x) = ln(a)x (exponential functions)
4. f (x) = eax for a R and x R>0 f (x) = ax (natural exponential functions)
5. f (x) = loga (x) for a R>0 \{1} and x R>0

1
f (x) = (logarithmic functions)
ln(a) loga (x)
1
6. f (x) = ln(x) for x R>0 f (x) = (natural logarithmic function).
ln(x)
In view of these results, we would like to emphasise the fact that for the entire family of gen-
eral power-law functions the elasticity f (x) has a constant value, independent of the value of
the argument x. It is this very property which classifies general power-law functions as scale-
invariant. When scale-invariance obtains, dimensionless ratios, i.e., quotients of variables of the
3
To extend the regime of applicability of the measure f , one may consider working in terms of absolute values |x|
and |f (x)|. Then one has to distinguish between four cases, which need to be looked at separately: (i) x > 0, f (x) > 0,
(ii) x < 0, f (x) > 0, (iii) x < 0, f (x) < 0 and (iv) x > 0, f (x) < 0.
same physical dimension, reduce to constants. In this context, we would like to remark that scale-
invariant (fractal) power-law functions of the form f (x) = Kx , with K > 0 and R<0 \{1},
are frequently employed in Economics and the Social Sciences for modelling uncertainty of eco-
nomic agents in decision-making processes, or for describing probability distributions of rare
event phenomena; see, e.g., Taleb (2007) [25, p 326ff] or Gleick (1987) [13, Chs. 5 and 6]. This
is due, in part, to the curious property that for certain values of the exponent general power-law
probability distributions attain unbounded variance; cf. Ref. [12, Sec. 8.9].
Practical applications in economic theory of the concept of an elasticty as a measure of relative
change of a differentiable real-valued function f of one real variable x are generally based on the
following linear (!) approximation: beginning at x0 > 0, for small prescribed percentage changes
x
of the argument x in the interval 0 % < 5 %, the resultant percentage changes of f amount
x0
approximately to
(percentage change of f ) (elasticity of f at x0 ) (percentage change of x) , (7.45)
or, in terms of a mathematical formula, to

f (x0 + x) f (x0 ) x
f (x0 ) . (7.46)
f (x0 ) x0
We now draw the readers attention to a special kind of terminology developed in economic theory
to describe the relative local change behaviour of economic functions in qualitative terms. For
x D(f ), the relative local change behaviour of a function f is called
inelastic, whenever |f (x)| < 1,
unit elastic, when |f (x)| = 1, and
elastic, whenever |f (x)| > 1.

For example, a total cost function K(x) in the diminishing returns picture exhibits unit elastic
behaviour at the minimum efficient scale x = xg2 where, by Eq. (7.22), K (xg2 ) = 1. Also, at
the local maximum of an average profit function G(x)/x, the property G (x) = 1 applies; cf.
Eq. (7.38).
Next, we review the computational rules one needs to adhere to when calculating elasticities for
combinations of two real-valued functions of one real variable in the sense of Sec. 7.1.6:
Computational rules for elasticities
If f and g are differentiable real-valued functions of one real variable, with elasticities f and g ,
it holds that:
1. product f g: f g (x) = f (x) + g (x),
f
2. quotient , g 6= 0: f /g (x) = f (x) g (x),
g
3. concatenation f g: f g (x) = f (g(x)) g (x),

1
4. inverse function f 1 : f 1 (x) = .
f (y) y=f 1 (x)
To end this chapter, we remark that for a positive differentiable real-valued function f of one
positive real variable x, a second elasticity may be defined according to

d x df (x)
f [f (x)] := x . (7.47)
dx f (x) dx
Of course, by analogy this procedure may be generalised to higher derivatives of f still.

Chapter 8
Integral calculus of real-valued functions of

one real variable
In the final chapter of these lecture notes we give a brief overview of the main definitions and
laws of the integral calculus of real-valued functions of one variable. Subsequently we consider
a simple application of this tool in economic theory.
8.1 Indefinite integrals

Def.: Let f be a continuous real-valued function of one real variable and F a differentiable real-
valued function of the same real variable, with D(f ) = D(F ). Given that f and F are related
according to
F (x) = f (x) for all x D(f ) , (8.1)
then F is referred to as a primitive of f .
Remark: The primitive of a given continuous real-valued function f cannot be unique. By the
rules of differentiation discussed in Sec. 7.2, besides F also F + c, with c R a real-valued
constant, constitutes a primitive of f since (c) = 0.
Def.: If F is a primitive of a continuous real-valued function f of one real variable, then
Z
f (x) dx = F (x) + c , c = constant R , with F (x) = f (x) (8.2)
defines the indefinite integral of the function f . The following names are used to refer to the
different ingredients in this expression:
x the integration variable,
f (x) the integrand,
dx the differential, and, lastly,
71
72 CHAPTER 8. INTEGRAL CALCULUS OF REAL-VALUED FUNCTIONS
c the constant of integration.
For the elementary, continuous real-valued functions of one variable introduced in Sec. 7.1, the
following rules of indefinite integration apply:
Rules of indefinite integration
R
1. dx = x + c with = constant R (constants)
Z
x2
2. x dx = +c (linear functions)
2
Z
xn+1
3. xn dx = + c for n N (natural power-law functions)
n+1
Z
x+1
4. x dx = + c for R\{1} and x R>0 (general power-law functions)
+1
Z
ax
5. ax dx = + c for a R>0 \{1} (exponential functions)
ln(a)
Z
eax
6. eax dx = + c for a R\{0} (natural exponential functions)
a
R
7. x1 dx = ln |x| + c for x R\{0}.
Special methods of integration need to be employed when the integrand consists of a concatanation
of elementary real-valued functions. Here we provide a list with the main tools for this purpose.
For differentiable real-valued functions f and g, it holds that
R R R
1. (f (x) g(x)) dx = f (x) dx g(x) dx
with , = constant R (summation rule)
R R
2. f (x)g (x) dx = f (x)g(x) f (x)g(x) dx (integration by parts)
u=g(x) and du=g (x)dx
R z}|{ R
3. f (g(x))g (x) dx = f (u) du = F (g(x)) + c (substitution method)
Z
f (x)
4. dx = ln |f (x)| + c for f (x) 6= 0 (logarithmic integration).
f (x)
8.2 Definite integrals

Def.: Let f be a real-valued function of one variable which is continuous on an interval [a, b]
D(f ), and let F be a primitive of f . Then the expression
Z b
f (x) dx = F (x)|x=b
x=a = F (b) F (a) (8.3)
a
defines the definite integral of f in the limits of integration a and b.

For definite integrals the following general rules apply:
8.3. APPLICATIONS IN ECONOMIC THEORY 73
Z a
1. f (x) dx = 0 (identical limits of integration)
a
Z a Z b
2. f (x) dx = f (x) dx (interchange of limits of integration)
b a
Z b Z c Z b
3. f (x) dx = f (x) dx + f (x) dx for c [a, b] (split of integration interval).
a a c
Remark: The main qualitative difference between an (i) indefinite integral and a (ii) definite in-
tegral of a continuous real-valued function of one variable reveals intself in the different kinds of
outcome: while (i) yields as a result a real-valued (primitive) function, (ii) simply yields a single
real number.
GDC: For a stored real-valued function, the evaluation
R of a definite integral can be performed in
mode CALC with the pre-programmed function f(x)dx. The corresponding limits of integration
need to be specified interactively.
As indicated in Sec. 7.6, the scale-invariant power-law functions f (x) = x for R and
x R>0 play a special role in practical applications. For x [a, b] R>0 and 6= 1 it holds
that Z b x=b
x+1 1 +1 +1

x dx = = b a . (8.4)
a + 1 x=a + 1
Problematic in this context can be considerations of taking limits of the form a 0 resp. b ,
since for either of the two cases
(i) case < 1: Z b
lim x dx , (8.5)
a0 a
(ii) case > 1: Z b

lim x dx , (8.6)
b a
one ends up with divergent mathematical expressions.
8.3 Applications in economic theory

The starting point shall be a simple market situation for a single product. For this product, on the
one-hand side, there be a demand function N(p) (in units) which is monotonously decreasing on
the price interval [pu , po ]; the limit values pu and po denote the minimum and maximum prices per
unit (in CU/u) acceptable for the product. On the other hand, the market situation be described by
a supply function A(p) (in units) which is monotonously increasing on [pu , po ].
The equilibrium unit price pM (in CU/unit) for this product is defined by assuming a state of
economic equilibrium of the market, quantitatively expressed by the condition
A(pM ) = N(pM ) . (8.7)
74 CHAPTER 8. INTEGRAL CALCULUS OF REAL-VALUED FUNCTIONS
Geometrically, this condition defines common points of intersection for the functions A(p) and
N(p) (when they exist).
GDC: Common points of intersection for stored functions f and g can be easily determined inter-
actively in mode CALC employing the routine intersect.
In a drastically simplified fashion, we now turn to compute the revenue made on the market by the
suppliers of a new product for either of three possible strategies of market entry.
1. Strategy 1: The revenue obtained by the suppliers when the new product is being sold
straight at the equilibrium unit price pM , in an amount N(pM ), is simply given by
U1 = U(pM ) = pM N(pM ) (in CU) . (8.8)
2. Strategy 2: Some consumers would be willing to purchase the product intially also at a unit
price which is higher than pM . If, hence, the suppliers decide to offer the product initially
at a unit price po > pM , and then, in order to generate further demand, to continuously1 (!)
reduce the unit price to the lower pM , the revenue obtained yields the larger value
Z po
U2 = U(pM ) + N(p) dp . (8.9)
pM
Since the amount of money

Z po
K := N(p) dp (in CU) (8.10)
pM
is (theoretically) safed by the consumers when the product is introduced to the market ac-
cording to strategy 1, this amount is referred to in the economic literature as consumer
surplus.
3. Strategy 3: Some suppliers would be willing to introduce the product to the market initially
at a unit price which is lower than pM . If, hence, the suppliers decide to offer the product
initially at a unit price pu < pM , and then to continuously2 (!) raise it to the higher pM , the
revenue obtaines amounts to the smaller value
Z pM
U3 = U(pM ) A(p) dp . (8.11)
pu
Since the suppliers (theoretically) earn the extra amount

Z pM
P := A(p) dp (in CU) (8.12)
pu
when the product is introduced to the market according to strategy 1, this amount is referred
to in the economic literature as producer surplus.
1
This is a strong mathematical assumption aimed at facilitating the actual calculation to follow.
2
See previous footnote.
Appendix A
Glossary of technical terms (GB D)
A
absolute change: absolute Anderung
absolute value: Betrag
account balance: Kontostand
addition: Addition
analysis: Analysis, Untersuchung auf Differenzierbarkeitseigenschaften
arithmetical mean: arithmetischer Mittelwert
arithmetical sequence: arithmetische Zahlenfolge
arithmetical series: arithmetische Reihe
augmented coefficient matrix: erweiterte Koeffizientenmatrix
average costs: Stuckkosten
average profit: Durchschnittsgewinn, Gewinn pro Stuck
B
backward substitution: ruckwertige Substitution
balance equation: Bilanzgleichung
basis: Basis
basis solution: Basislosung
basis variable: Basisvariable
Behavioural Economics: Verhaltensokonomik
boundary condition: Randbedingung
break-even point: Gewinnschwelle
C
chain rule: Kettenregel
characteristic equation: charakteristische Gleichung
coefficient matrix: Koeffizientenmatrix
column: Spalte
column vector: Spaltenvektor
component: Komponente
compound interest: Zinseszins
concatenation: Verschachtelung, Verknupfung
conservation law: Erhaltungssatz
75
76 APPENDIX A. GLOSSARY OF TECHNICAL TERMS (GB D)
constant of integration: Integrationskonstante
constraint: Zwangsbedingung
cost function: Kostenfunktion
consumer surplus: Konsumentenrente
continuity: Stetigkeit
contract period: Laufzeit
Cournots point: Cournotscher Punkt
curve sketching: Kurvendiskussion
D
decision-making: Entscheidungsfindung
declining-balance depreciation method: geometrischdegressive Abschreibung
definite integral: bestimmtes Integral
demand function: Nachfragefunktion
dependent variable: abhangige Variable
depreciation: Abschreibung
depreciation factor: Abschreibungsfaktor
derivative: Ableitung
determinant: Determinante
difference: Differenz
difference quotient: Differenzenquotient
differentiable: differenzierbar
differential: Integrationsdifferenzial
differential calculus: Differenzialrechnung
dimension: Dimension
direction of optimisation: Optimierungsrichtung
divergent: divergent, unbeschrankt
domain: Definitionsbereich
E
economic agent: Wirtschaftstreibende(r) (meistens ein homo oeconomicus)
economic efficiency: Wirtschaftlichkeit
economic equilibrium: okonomisches Gleichgewicht
economic principle: okonomisches Prinzip
economic theory: Wirtschaftstheorie
Econophysics: Okonophysik
eigenvalue: Eigenwert
eigenvector: Eigenvektor
elastic: elastisch
elasticity: Elastizitat
element: Element
end of profitable zone: Gewinngrenze
endogenous: endogen
equilibrium price: Marktpreis
equivalence transformation: Aquivalenztransformation
77
exogenous: exogen
exponential function: Exponentialfunktion
extrapolation: Extrapolation, uber bekannten Gutigkeitsbereich hinaus verallgemeinern
F
feasible region: zulassiger Bereich
final capital: Endkapital
fixed costs: Fixkosten
forecasting: Vorhersagen erstellen
function: Funktion
G
Gauian elimination: Gauscher Algorithmus
GDC: GTR, grafikfahiger Taschenrechner
geometrical mean: geometrischer Mittelwert
geometrical sequence: geometrische Zahlenfolge
geometrical series: geometrische Reihe
growth function: Wachstumsfunktion
H
I
identity: Identitat
image vector: Absolutgliedvektor
indefinite integral: unbestimmtes Integral
independent variable: unabhangige Variable
inelastic: unelastisch
initial capital: Anfangskapital
installment: Ratenzahlung
installment savings: Ratensparen
integral calculus: Integralrechnung
integrand: Integrand
integration variable: Integrationsvariable
interest factor: Aufzinsfaktor
interest rate: Zinsfu
inverse function: Inversfunktion, Umkehrfunktion
inverse matrix: inverse Matrix, Umkehrmatrix
isoquant: Isoquante
J
L
law of diminishing returns: Ertragsgesetz
length: Lange
level of physical output: Ausbringungsmenge
limits of integration: Integrationsgrenzen
linear combination: Linearkombination
linearisation: Linearisierung
linear programming: lineare Optimierung
local rate of change: lokale Anderungsrate
logarithmic function: Logarithmusfunktion
M
mapping: Abbildung
marginal costs: Grenzkosten
maximisation: Maximierung
minimisation: Minimierung
minimum efficient scale: Betriebsoptimum
monetary value: Geldwert
monopoly: Monopol
monotonicity: Monotonie
mortgage loan: Darlehen
N
non-basis variable: Nichtbasisvariable
non-negativity constraints: Nichtnegativitatsbedingungen
non-linear functional relationship: nichtlineare Funktionalbeziehung
O
objective function: Zielfunktion
one-to-one and onto: eineindeutig
optimal solution: optimalen Losung
optimal value: optimaler Wert
optimisation: Optimierung
order-of-magnitude: Groenordnung
orthogonal: orthogonal
over-determined: uberbestimmt
P
parallel displacement: Parallelverschiebung
pension calculations: Rentenrechnung
percentage rate: Prozentsatz
perfect competition: totale Konkurrenz
period: Periode
pivot column index: Pivotspaltenindex
pivot element: Pivotelement
pivot operation: Pivotschritt
pivot row index: Pivotzeilenindex
pole: Polstelle, Singularitat
polynomial division: Polynomdivision
79
polynomial of degree n: Polynom vom Grad n
power-law function: Potenzfunktion
present value: Barwert
primitive: Stammfunktion
principal component analysis: Hauptkomponentenanalyse
producer surplus: Produzentenrente
product rule: Produktregel
profit function: Gewinnfunktion
prohibitive price: Pohibitivpreis
Prospect Theory: Neue Erwartungstheorie
psychological value function: psychologische Wertfunktion
Q
quadratic matrix: quadratische Matrix
quotient: Quotient
quotient rule: Quotientenregel
R
range: Wertespektrum
rank: Rang
rare event: seltenes Ereignis
rational function: gebrochen rationale Funktion
real-valued function: reellwertige Funktion
reference period: Referenzzeitraum
regression analysis: Regressionsanalyse
regular: regular
relative change: relative Anderung
remaining debt: Restschuld
remaining resources: Restkapazitaten
remaining value: Restwert
rescaling: Skalierung
resources: Rohstoffe
resource consumption matrix: Rohstoffverbrauchsmatrix
restrictions: Restriktionen
root: Nullstelle
row: Reihe
row vector: Zeilenvektor
S
saturation quantity: Sattigungsmenge
scale: Skala, Groenordnung
scale-invariant: skaleninvariant
simplex: Simplex, konvexer Polyeder
simplex tableau: Simplextabelle
singular: singular
sink: Senke
slack variable: Schlupfvariable
source: Quelle
stationary: stationar, konstant in der Zeit
straight line depreciation method: lineare Abschreibung
strictly monotonously decreasing: streng monoton fallend
strictly monotonously increasing: streng monoton steigend
summation rule: Summationsregel
supply function: Angebotsfunktion
T
tangent: Tangente
target space: Wertebereich
technology matrix: Technologiematrix
total demand matrix: Gesamtbedarfsmatrix
total revenue: Ertrag
transpose: Transponierte
U
uncertainty: Unsicherheit
under-determined: unterbestimmt
uniqueness: Eindeutigkeit
unit elastic: proportional elastisch
unit matrix: Einheitsmatrix
unit price: Stuckpreis
unit vector: Einheitsvektor
utility function: Nutzenfunktion
V
value chain: Wertschopfungskette
variability: Anderungsverhalten, Variabilitat
variable average costs: variable Stuckkosten
variable costs: variable Kosten
variable vector: Variablenvektor
vector: Vektor
vector algebra: Vektoralgebra
W
well-determined: wohlbestimmt
Z
zero matrix: Nullmatrix
zero vector: Nullvektor
Bibliography
[1] T Arens, F Hettlich, C Karpfinger, U Kockelkorn, K Lichtenegger and H Stachel (2008)

Mathematik 1st Edition (Heidelberg: Spektrum) ISBN13: 9783827417589
[2] A Asano (2013) An Introduction to Mathematics for Economics (Cambridge: Cambridge

University Press) ISBN13: 9780521189460
[3] C Bauer, M Clausen, A Kerber and H MeierReinhold (2008) Mathematik fur Wirtschaftswis-
senschaftler 5th Edition (Stuttgart: SchafferPoeschel) ISBN13: 9783791027487
[4] Bernoulli D (1738) Specimen theoriae novae de mensura sortis Papers of the Imperial
Academy of Sciences in St. Petersburg
English translation: (1954) Exposition of a new theory on the measurement of risk
Econometrica 22 2336
[5] JP Bouchaud and M Potters (2003) Theory of Financial Risk and Derivative Pricing 2nd
Edition (Cambridge: Cambridge University Press) ISBN13: 9780521741866
[6] K Bosch (2003) Mathematik fur Wirtschaftswissenschaftler 14th Edition (Munchen: Olden-
bourg) ISBN10: 3486273094
[7] I N Bronstein, K A Semedjajew, G Musiol and H Muhlig (2005) Taschenbuch der Mathematik
6th Edition (Frankfurt (Main): Harri Deutsch) ISBN10: 3817120060
[8] G B Dantzig (1949) Programming of interdependent activities: II mathematical model

[9] G B Dantzig (1955) Linear programming under uncertainty Management Science 1 197206
[10] E T Dowling (1990) Schaums Outline of Calculus for Business, Economics, and the Social
Sciences (New York: McGrawHill) ISBN13: 9780070176737
[11] E T Dowling (2009) Schaums Outline of Mathematical Methods for Business and Economics
(New York: McGrawHill) ISBN13: 9780071635325
[12] H van Elst (2015) Foundations of descriptive and inferential statistics Preprint
arXiv:1302.2525v3 [stat.AP]
[13] J Gleick (1987) Chaos Making a New Science nth Edition 1998 (London: Vintage) ISBN
13: 9780749386061
81
82 BIBLIOGRAPHY
[14] J F Hair jr, W C Black, B J Babin and R E Anderson (2010) Multivariate Data Analysis 7th
Edition (Upper Saddle River (NJ): Pearson) ISBN13: 9780135153093
[15] F Herrmann (2003) Physik III: Thermodynamik Skripten zur Experimentalphysik Univer-
sitat Karlsruhe (TH) URL (cited on August 30, 2015):
www.physikdidaktik.uni-karlsruhe.de/skripten/thermod.pdf
[16] J Hulsmann, W Gamerith, U LeopoldWildburger and W Steindl (2005) Einfuhrung in die

Wirtschaftsmathematik 4th Edition (Berlin: Springer) ISBN10: 3540244093
[17] D Kahneman (2011) Thinking, Fast and Slow (London: Penguin) ISBN13: 9780141033570
[18] D Kahneman and A Tversky (1979) Prospect theory: an analysis of decision under risk
[19] S J Leon (2009) Linear Algebra with Applications 8th Edition (London: Pearson) ISBN13:
9780136009290
[20] W W Leontief (1936) Quantitative input and output relations in the economic systems of the
United States The Review of Economics and Statistics 18 105125
[21] R Penrose (2004) The Road to Reality A Complete Guide to the Laws of the Universe 1st
Edition (London: Jonathan Cape) ISBN10: 0224044478
[22] H Rinne (2008) Taschenbuch der Statistik 4th Edition (Frankfurt Main): Harri Deutsch)
ISBN13: 9783817118274
[23] S Singh (1997) Fermats Last Theorem (London: Fouth Estate) ISBN13: 9781857025217
[24] L Smith (2007) Chaos A Very Short Introduction (Oxford: Oxford University Press)
ISBN13: 9780192853783
[25] N N Taleb (2007) The Black Swan The Impact of the Highly Improbable (London: Pen-
guin) ISBN13: 9780141034591
[26] A ZehMarschke (2010) Finanzmathematik Zins und Zinseszins, Tilgungs- und Renten-
rechnung, Investitionsrechnung, Abschreibungen private lecture notes
View publication stats

1509 04333v2 PDF

Uploaded by

Copyright:

Available Formats

1509 04333v2 PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

1509 04333v2 PDF

Uploaded by

Copyright:

Available Formats

See

An Introduction to Business Mathematics

Book September 2015

Henk van Elst

statsVideos View project

Non-knowledge and uncertainty View project

Lecture notes for the Bachelor degree programmes

0.1.1 IMQM: Introduction to Management and its Quantitative Methods

H ENK VAN E LST

August 30, 2015

Fakultat I: Betriebswirtschaft und Management

EPrint: arXiv:1509.04333v2 [q-fin.GN]

Cite as: arXiv:1509.04333v2 [q-fin.GN]

These lecture notes were typeset in LATEX 2 .

Qualification objectives of the module (excerpt) 1

1 Vector algebra in Euclidian space Rn 5

3 Systems of linear algebraic equations 17

4 Leontiefs inputoutput matrix model 25

6 Elementary financial mathematics 39

7 Differential calculus of real-valued functions 51

8 Integral calculus of real-valued functions 71

A Glossary of technical terms (GB D) 75

The qualification objectives shall be reached by an integrative approach.

Texas Instruments TI84 plus,

Casio CFX9850GB PLUS.

Vector algebra in Euclidian space Rn

1.1 Basic concepts

to constitute related zero vectors.

|| > 1 stretching of the length of a

0 < || < 1 shrinking of the length of a

< 0 directional reversal of a.

3. a + 0 = a (addition identity element)

4. For every a, b Rn , there exists exactly one x Rn such that a + x = b

5. ()a = (a) with R (associative rescaling)

6. 1a = a (rescaling identity element)

1.2 Dimension and basis of Rn

is referred to as a linear combination of the m vectors ai , i = 1, . . . , m.

Remark: The n unit vectors

1.3 Euclidian scalar product

1. (a + b)T c = aT c + bT c (distributive scalar product)

2. aT b = bT a (commutative scalar product)

3. (aT ) b = (aT b) with R (homogeneous scalar product)

4. aT a > 0 for all a 6= 0 (positive definite scalar product).

Now we turn to introduce the notion of the length of an n-component vector.

|a| 0, and |a| = 0 a = 0;

|a| = |||a| for R;

|a + b| |a| + |b| (triangle inequality).

2.1 Matrices as linear mappings

Def.: A real-valued (m n)-matrix is formally defined to constitute an array of real numbers

Characteristic features of this array of real numbers are:

m denotes the number of rows of A, n the number of columns of A.

holds a special status in the family of (n n)-matrices.

A(x1 + x2 ) = (Ax1 ) + (Ax2 )

they constitute linear mappings.1

2.2 Basic concepts

AT : aTij := aji , (2.6)

where i = 1, . . . , m und j = 1, . . . , n. Note that it holds that AT Rnm .

Def.: Addition of matrices

A =: C : aij =: cij , (2.9)