Business Maths

A N I NTRODUCTION TO
B USINESS M ATHEMATICS
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
Fakultät I: Betriebswirtschaft und Management

Karlshochschule
International University
Karlstraße 36–38
76133 Karlsruhe
Germany
E–mail: [email protected]
E–Print: arXiv:1509.04333v2 [q-fin.GN]
c 2009–2015 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, Leontief’s stationary input–output matrix model, linear
programming, elementary financial mathematics, as well as differential and integral cal-
culus 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 7

1.1 Basic concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2 Dimension and basis of Rn . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3 Euclidian scalar product . . . . . . . . . . . . . . . . . . . . . . . . . 10
2 Matrices 13
2.1 Matrices as linear mappings . . . . . . . . . . . . . . . . . . . . . . . 13
2.2 Basic concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.3 Matrix multiplication . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3 Systems of linear algebraic equations 19

3.1 Basic concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.2 Gaußian elimination . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.3 Rank of a matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.4 Solution criteria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.5 Inverse of a regular (n × n)-matrix . . . . . . . . . . . . . . . . . . . . 23
3.6 Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4 Leontief’s input–output matrix model 27

4.1 General considerations . . . . . . . . . . . . . . . . . . . . . . . . . . 27
4.2 Input–output matrix and resource consumption matrix . . . . . . . 28
4.2.1 Input–output matrix . . . . . . . . . . . . . . . . . . . . . . . . 29
4.2.2 Resource consumption matrix . . . . . . . . . . . . . . . . . . 31
4.3 Stationary linear flows of goods . . . . . . . . . . . . . . . . . . . . . 31
4.3.1 Flows of goods: endogenous INPUT to total OUTPUT . . . . . 31
4.3.2 Flows of goods: exogenous INPUT to total OUTPUT . . . . . . 32
4.4 Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
CONTENTS
5 Linear programming 35
5.1 Exposition of a quantitative problem . . . . . . . . . . . . . . . . . . 35
5.2 Graphical solution method . . . . . . . . . . . . . . . . . . . . . . . . 37
5.3 Dantzig’s simplex algorithm . . . . . . . . . . . . . . . . . . . . . . . 38
6 Elementary financial mathematics 41

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

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

8.1 Indefinite integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
8.2 Definite integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
8.3 Applications in economic theory . . . . . . . . . . . . . . . . . . . . . 76
A Glossary of technical terms (GB – D) 79
Bibliography 85
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 man-
ager 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 Karlshochschule International University. The aim is
to provide a selection of tried-and-tested mathematical tools that proved efficient
in actual practical problems of Economics and Management. These tools con-
stitute the foundation for a systematic treatment of the typical kinds of quantita-
tive problems one is confronted with in a Bachelor degree programme. Neverthe-
less, 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 arithmetical 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 for-
mulae, determining the point of intersection for two straight lines in the Euclid-
ian plane, solving 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 cur-
rent models used in public schools and in undergraduate studies are, amongst
others,
• Texas Instruments TI–84 plus,
• Casio CFX–9850GB PLUS.
However, the reader is strongly encouraged to think about resorting, as an alter-
native, to a spreadsheet programme such as EXCEL or OpenOffice to handle
the calculations one encounters in one’s quantitative work.
The central theme of these lecture notes is the acquisition and application of a
number of effective mathematical methods in a business oriented environment.
In particular, we hereby focus on quantitative processes of the sort
INPUT → OUTPUT ,
3
4 CONTENTS
for which different kinds of functional relationships between some numerical
INPUT quantities and some numerical OUTPUT quantities are being consid-
ered. Of special interest in this context 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 de-
crease 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 vari-
ables, 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 program-
ming. In Part II, which is limited to Ch. 6, we turn to discuss elementary aspects
of Financial Mathematics. Fundamental principles of Analysis, comprising dif-
ferential and integral 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], Dowling (2009) [11], Dowling (1990) [10], Bauer et al (2008) [3], Bosch
(2003) [6], and Hülsmann 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 effi-
ciency of the mathematical methods presented, and, hence, would like to know
more about their background and relevance, as well as being introduced to fur-
ther 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 men-
tioned in this Introduction are available from the library at Karlshochschule In-
ternational University.
Finally, we draw the reader’s attention to the fact that the *.pdf version of these
lecture notes contains interactive features such as fully hyperlinked references
to original publications at the websites dx.doi.org and jstor.org, as well as
active links to biographical information on scientists that have been influential
CONTENTS 5
in the historical development of Mathematics, hosted by the websites The Mac-
Tutor History of Mathematics archive (www-history.mcs.st-and.ac.uk) and
en.wikipedia.org.
6 CONTENTS
Chapter 1
Vector algebra in Euclidian space Rn
Let us begin our elementary considerations of vector algebra with the introduc-
tion 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 sub-
jected 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 ∈ Rn×1 ,
7
8 CHAPTER 1. VECTOR ALGEBRA IN EUCLIDIAN SPACE RN
and
Def.: Real-valued row vector with n components
aT := (a1 , . . . , ai , . . . , an ) , ai ∈ R, i = 1, . . . , n , (1.3)
Notation: aT ∈ R1×n .
Correspondingly, we define the n-component objects
0
 
 ..
 .


0 :=  0 and 0T := (0, . . . , 0, . . . , 0) (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 + b 1
     
 ..   ..   ..
 .   .   .


a + b :=  ai  +  b i  =  ai + b i  , ai , bi ∈ R , (1.5)
     
 .   ..  ..
 ..  
 
.   . 
an bn an + b n
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

1.2. DIMENSION AND BASIS OF RN 9
• λ < 0 — directional reversal of a.
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 ad-
dition and multiplication 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 ∈ Rn (1.7)
i=1
is referred to as a linear combination of the m vectors ai , i = 1, . . . , m.

1
This is named after the ancient greek mathematician Euclid of Alexandria (about 325 BC–
265 BC).
2
A slightly s horter notation for n-component column vectors a ∈ Rn×1 is given by a ∈ Rn ;
likewise aT ∈ Rn for n-component row vectors aT ∈ R1×n .
Def.: A set of m vectors a1 , . . . , ai , . . . , am ∈ Rn is called linearly independent
when the condition
m
!
X
0 = λ1 a1 + . . . + λi ai + . . . + λm am = λi ai , (1.8)
i=1
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 independent. This maximum number is referred to as the dimension
of Euclidian space Rn . Every set of n linearly independent vectors in Euclid-
ian 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 e i . (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 .
1.3. EUCLIDIAN SCALAR PRODUCT 11
Def.: For an n-component row vector aT ∈ R1×n and an n-component column
vector b ∈ Rn×1 , the Euclidian scalar product
b1
 
 .. 
 .  Xn
T
a · b := (a1 , . . . , ai , . . . an )  bi  = a1 b1 + . . . + ai bi . . . + an bn =: ai b i (1.12)
 
 . 
 ..  i=1
bn
defines a mapping f : R1×n × Rn×1 → 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) are referred to as mutually orthogonal when they exhibit the property
that 0 = aT · b = bT · 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)] = · = âT · b̂ ⇒ ϕ(a, b) = cos−1 (âT · b̂) . (1.15)
|a| |b|
Remark: The inverse cosine function3 cos−1 (. . .) is available on every standard

GDC or spreadsheet.
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 mathe-
matical operations, and a set of computational 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 ∈ Rm×n .
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 ),
13
14 CHAPTER 2. MATRICES
a1j
 
 a2j 
 . 
 . 
 . 
elements of the 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 ∈ Rm×n 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 15
a mapping A : Rn×1 → Rm×1 , i.e. a mapping from the set of real-valued n-
component column vectors (here: x) to the set of real-valued m-component col-
umn 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 ∈ Rn×1 and real numbers λ ∈ R, mappings
defined by real-valued matrices A ∈ Rm×n 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 ∈ Rm×n , 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 ∈ Rn×m .

When transposing an (m × n)-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.

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.: Addition of matrices
For A, B ∈ Rm×n , 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.
2.3. MATRIX MULTIPLICATION 17
Def.: Rescaling of matrices
For A ∈ Rm×n 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 ∈ Rm×n :
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 multiplication 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.
The element of C at the intersection of the ith row and the jth column is deter-
mined 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 pro-
vided 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. Other-
wise, a matrix multiplication cannot be defined in a meaningful way.
GDC: For matrices [A] and [B] edited beforehand, of matching formats, their ma-
trix 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)
3. A |{z}
1 = |{z}
1 A=A (multiplicative identity element)
∈Rn×n ∈Rm×m
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 sin-
gle 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,
19
20 CHAPTER 3. SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS
• 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
• 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 Repre-
sentation 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 pre-
viously 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ß (1777–1855).
3.2 Gaußian 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 system:
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.

3.2. GAUSSIAN ELIMINATION 21
Specifically, this implies that we may manipulate a given linear system by means
of these four different kinds of equivalence transformations without ever chang-
ing its identity. In concrete cases, however, one should not apply these equiv-
alence transformations at random but rather follow a target oriented strategy.
This is what Gaußian elimination can provide.
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 . . . ã1j ... ã1n b̃1

.. . . . .. .. ..
. . .. . . .
0 . . . ãij ... ãin b̃i , (3.5)
.. . . . ... .. ..
. . .. . .
0 ... 0 . . . ãmn b̃m
by means of the four kinds of equivalence transformations such that the resul-
tant simpler final linear system may easily be solved using backward substitu-
tion.
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 ∈ Rm×n 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 ∈ Rn×n , 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 ∈ R2×2 , 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 A’s on-diagonal elements and
A’s off-diagonal elements.
(ii) When A ∈ R3×3 , the definition of A’s 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 ele-
ments 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 ∈ Rn×n , please refer to the literature; e.g. Bronstein
et al (2005) [7, p 267].
To determine the rank of a given quadtratic matrix A ∈ Rn×n , 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. SOLUTION CRITERIA 23
3.4 Criteria for solving systems of linear algebraic

equations
Making use of the concept of the rank of a real-valued matrix A ∈ Rm×n , 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 ∈ Rm×n , variable vector x ∈ Rn×1 and image vector
b ∈ Rm×1 , there exist(s)
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 ∈ Rn×n be regular, i.e., det(A) ∈ R\{0}.
Then there exists an inverse matrix A−1 to A defined by the characterising
properties
A−1 A = AA−1 = 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 A−1 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 Gaußian elimination.
GDC: For a stored quadratic matrix [A], its inverse matrix can be simply obtained
as [A]−1 , where the x−1 function key needs to be used.
3.6. OUTLOOK 25
Computational rules for the inverse operation
For A, B ∈ Rn×n , with det(A) 6= 0 6= det(B), it holds that
1. (A−1 )−1 = A
2. (AB)−1 = B−1 A−1
3. (AT )−1 = (A−1 )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 ∈ Rn×n , i.e., det(A) 6= 0. Then, for A,
there exists an inverse matrix A−1 . Matrix-multiplying both sides of the equation
above from the left (!) by the inverse A−1 , results in
A−1 (Ax) = (A−1 A)x = 1x = x = A −1

| {z b} . (3.11)
| {z }
left-hand side right-hand side
In this case, the unique solution (!) x = A−1 b 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 A−1 .)
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 ∈ Rn×n , do there exist real numbers λn ∈ R
and real-valued vectors v n ∈ Rn×1 which satisfy the condition
!
Av n = λn v n ? (3.12)
Put differently: for which vectors v n ∈ Rn×1 does their mapping by a quadratic
matrix A ∈ Rn×n 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)
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 ∈ Rn×1 of a real-valued matrix A ∈ Rn×n provides the basis of a transformation
of A to its diagonal form Aλn , 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
Leontief’s stationary input–output

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 simpli-
fied 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 (1905–1999), 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 eco-
nomic agents exchanging between them the goods they produce. For simplicity
we want to assume that every one of these economic agents represents the pro-
duction 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 consti-
tute 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 sector’s
total OUTPUT during the chosen reference period of the n goods can be viewed
to flow through one of two separate channels: (i) the exogenous channel linking
27
28 CHAPTER 4. LEONTIEF’S INPUT–OUTPUT MATRIX MODEL
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 dif-
ferent actors, by the prospect of increasing the value of the INPUT quantities,
in line with the notion of the economic value chain.
Leontief’s model is based on the following three elementary
Assumptions:
1. For all goods involved the functional relationship between INPUT and OUT-
PUT quantities be of a linear nature [cf. Eq. (2.5)].
2. The proportions of “INPUT quantities to OUTPUT quantities” be constant
over the reference 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 de-
manded.
The mathematical formulation of Leontief’s quantitative model employs the fol-
lowing
Vector- and matrix-valued quantities:
1. q — total output vector ∈ Rn×1 , components qi ≥ 0 units (dim: units)
2. y — final demand vector ∈ Rn×1 , components yi ≥ 0 units (dim: units)
3. P — input–output matrix ∈ Rn×n , components Pij ≥ 0 (dim: 1)
4. (1 − P) — technology matrix ∈ Rn×n , regular, hence, invertible (dim: 1)
5. (1 − P)−1 — total demand matrix ∈ Rn×n (dim: 1)
6. v — resource vector ∈ Rm×1 , components vi ≥ 0 units (dim: units)
7. R — resource consumption matrix ∈ Rm×n , components Rij ≥ 0, (dim: 1)
where 1 denotes the (n × n)-unit matrix [cf. Eq. (2.3)]. Note that the compo-
nents of all the vectors involved, as well as of the input–output matrix and of the
resource consumption matrix, can assume non-negative values (!) only.
4.2 Input–output matrix and resource consumption

matrix
We now turn to provide the definition of the two central matrix-valued quantities
in Leontief’s model. We will also highlight their main characteristic features.
4.2. INPUT–OUTPUT MATRIX AND RESOURCE CONSUMPTION MATRIX 29
4.2.1 Input–output 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. 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 stan-

dard input–output table as follows:
[Values in units] agent 1 · · · agent j ··· agent n external consumers Σ: total ou

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
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 n1n 
n11 +...+n1j +...+n1n +y1
... nj1 +...+njj +...+njn +yj
... nn1 +...+nnj +...+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 Leontief’s input–output matrix of the stationary economic
system under investigation.
For the very simple case with just n = 3 producing agents, the input–output
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 input–output

matrix P can be determined only once the reference period of time chosen has
come to an end.
The utility of Leontief’s stationary input–output 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
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 .
4.3. STATIONARY LINEAR FLOWS OF GOODS 31
In this way it becomes possible to compute for a given (idealised) economic sys-
tem 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 experience 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 quantita-
tive aspects of stationary flows of goods, as we will soon elucidate.
4.2.2 Resource consumption matrix

The second matrix-valued quantity central to Leontief’s stationary model is the
resource consumption matrix R. This may be interpreted as providing a recipe
for the amounts of the m different kinds of external resources (the exogenous
INPUT quantities) that are needed in the production 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)
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 con-
sidered. According to Leontief’s 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 − y i = Pij qj , (4.10)
j=1
with i = 1, . . . , n, in which the input–output matrix P takes the role of medi-

ating a mapping between either of these vector-valued quantities. According to
Assumption 2, the elements of the input–output 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 sciences 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 considered 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 exogenous 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
4.4. OUTLOOK 33
with i = 1, . . . , m. By Assumption 2, the elements of the resource consump-
tion 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
v = Rq = R(1 − P)−1 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
Leontief’s input–output matrix model may be extended in a straightforward fash-
ion to include more advanced considerations of economic theory. Supposing
a closed though not necessarily stationary economic system G comprising n in-
terdependent economic agents producing n different 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 associ-
ated 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 corresponding flows of money. However, contrary
to the number of goods, in general there does not 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 principle possess either sign, positive or negative. To deal with
2
Here the symbols CU and TU denote “currency units” and “time units,” respectively.
these quantitative issues properly, one requires the technical tools of the dif-
ferential 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.
Leontief’s input–output 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.
Chapter 5
Linear programming
On the backdrop of the economic principle, we discuss in this chapter a special

class of quantitative problems that frequently arise in specific practical applica-
tions in Business and Management. 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 effort (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 con-
ditions. 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 vari-
ables 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 ∈ Rm×n , a vector b ∈ Rm×1 , two vectors c, x ∈ Rn×1 , and
a constant d ∈ R. A quantitative problem of the form
max z = cT · x + d |Ax ≤ b, x ≥ 0 ,

(5.1)
35
36 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 variables. 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 ob-

jective functions of n = 2 independent variables constitute straight lines, of n = 3
independent variables Euclidian planes, of n = 4 independent variables Euclid-
ian 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 37
5.2 Graphical method for solving problems with
two independent 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 de-

termined by the m restrictions imposed on x1 and x2 . Specifically, D consti-
tutes 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 determined 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 re-
gion 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 intersection 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 opti-

mal solution (x1O , x2O ) into the m restrictions.
In general one finds that for a linear objective function z with n = 2 indepen-
dent 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 corre-
sponding 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 independent variables by means of the graphical method, one needs to par-
allelly 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 re-
gion 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 Dantzig’s simplex algorithm

The main disadvantage of the graphical solution method is its limitation to prob-
lems 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 mathematician George Bernard Dantzig
(1914–2005) 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 poly-
hedron, i.e., a body of finite (hyper-)volume in two or more dimensions bounded
by linear (hyper-)surfaces which intersect 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 program-
ming 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 algo-
rithm such that it systematically 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 independent variables in the form of relations (5.2)–(5.6). First, by in-
troducing 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
5.3. DANTZIG’S SIMPLEX ALGORITHM 39
Maximum problem of linear programming in canonical form
z − c1 x 1 − c2 x 2 − . . . − cn x n = 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 (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.
Dantzig’s 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 sim-
plex tableau.
Initial simplex tableau
z x1 x2 . . . x n 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 respec-
tive 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 ini-
tial 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 sys-
tem (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 indepen-
dent variables 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 algorithm needs to be
viewed as a special variant of Gaußian 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 correspond-
ing to the optimal solution of the given linear programming problem are then
to be determined from the final simplex tableau by backward substitution, be-
ginning at the bottom row. Note that slack variables with positive values belong-
ing 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 mathematics. 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 )n∈N ,
is defined by the property that the difference d between neighbouring elements

in the sequence be constant, i.e., for n > 1
an − an−1 =: d = constant 6= 0 , (6.1)
with an , an−1 , d ∈ R. Given this recursive formation rule, one may infer the ex-
plicit representation 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 pa-
rameters 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.
41
42 CHAPTER 6. ELEMENTARY FINANCIAL MATHEMATICS
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
(an−1 + an+1 ) = (a1 + (n − 2)d + a1 + nd) = a1 + (n − 1)d = an . (6.3)
2 2
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 Gaußian identity1 (cf., e.g., Bosch
(2003) [6, p 21])
n−1
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 )n∈N ,
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)
an−1
with an , an−1 ∈ 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 n−1 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 ge-
ometrical mean of the immediate neighbours of any particular element an (with
n ≥ 2), one finds that
√ p
an−1 · an+1 = a1 q n−2 · a1 q n = a1 q n−1 = an . (6.8)
n
X
1
Analogously, the modified Gaußian identity (2k − 1) ≡ n2 applies.
k=1
6.2. INTEREST AND COMPOUND INTEREST 43
Summation of the first n elements of an arbitrary geometrical sequence of real
numbers leads to a finite geometrical series,
n n n−1
X X k−1 X qn − 1
Sn := a1 + a2 + . . . + an = ak = a1 q = a1 q k = a1 . (6.9)
k=1 k=1 k=0
q − 1
In the last algebraic step use was made of the identity (cf., e.g., Bosch (2003) [6,
p 27])
n−1
X qn − 1
qk ≡ for q ∈ R\{0, 1} . (6.10)
k=0
q−1
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 in-
stant, 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 cor-
responding 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 = Kn−1 q , n∈N, (6.12)
will have accumulated, where Kn−1 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 com-
pound 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 , n∈N. (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
2
Inverting this defining relation for q leads to p = 100 · (q − 1).
solve Eq. (6.13) to obtain the value of the fourth. For example, solving Eq. (6.13)
for K0 yields
Kn
K0 = n =: B0 . (6.14)
q
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)
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
m · 100
end of the full year, K0 in this interest model will have increased to
pnom m
K1 = Km/m = K0 1 + , m∈N.
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 , m∈N, (6.18)
m · 100
6.2. INTEREST AND COMPOUND INTEREST 45
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 n·m n
K n = K0 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 concept 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 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 Kn−1 , Kn−2 , etc.,
n−1
X
n 2 n−1
Kn = (Kn−1 + 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 , n∈N. (6.21)
q−1
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 = = n−1 . (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 remain-
ing debt take the values
p
Z1 = R0 · = R0 (q − 1)
100
T1 = A − Z1
substitute for Z1
p
R1 = R0 + Z1 − A = R0 + R0 · − A = R0 q − A .
z}|{
100
By the end of a second full year, these become
Z2 = R1 (q − 1)
T2 = A − Z2
substitute for Z2 substitute for R1
R2 = R1 + Z2 − A = R1 q − A = R0 q 2 − A(q + 1) .
z}|{ z}|{
At this stage, it has become clear according to which patterns the different quan-
tities involved in the redemption payment model need to be formed. The inter-
est payment for the nth full year in a mortgage contract of constant anuities
amounts to (recursively)
Zn = Rn−1 (q − 1) , n∈N, (6.26)

6.3. REDEMPTION PAYMENTS IN CONSTANT ANNUITIES 47
where Rn−1 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 , n∈N. (6.27)
The remaining debt at the end of the nth full year then is (in CU)
recursively: Rn = Rn−1 + Zn − A = Rn−1 q − A , n∈N. (6.28)
By successive backward substitution for Rn−1 , Rn−2 , etc., Rn can be re-expressed

as
n−1
X
n n−1 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 , n∈N. (6.29)
q−1
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
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 con-
stant annuity A contains implicitly a factor (q − 1) [cf. Eq. (6.25)], the two com-
peting 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 ≤ Rn−1 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 q−1
= 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 subsequent 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 important 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 throughout 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
6.4. PENSION CALCULATIONS 49
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)
1
K1 = K0 − ma + Z1 = K0 q − m + (m + 1)(q − 1) a .
z}|{
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
while the account balance amounts to
substitute for K1 and Z2
2 1
K2 = K1 − ma + Z2 = K0 q − m + (m + 1)(q − 1) a(q + 1) .
z}|{
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 = Rn−1 − (m + 1)a (q − 1) , (6.34)
2
and for Kn one obtains
substitute for Kn−1 and Zn Xn−1

n 1
Kn = Kn−1 − ma + Zn = K0 q − m + (m + 1)(q − 1) a qk .
z}|{
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 q−1
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 one’s con-
venience (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 [...]a−K0 (q−1)
n= . (6.36)
ln(q)
(ii) The present value B0 of a pension scheme results from the following con-
sideration: 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
qn − 1

1
B0 = K0 = m + (m + 1)(q − 1) a n . (6.37)
2 q (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
q−1
aever = 1 K0 ; (6.38)
m + (m + 1)(q + 1)
2
Note that, naturally, aever is directly proportional to the initial capital K0 !

3 1

To avoid notational overload, the brackets [. . .] here represent the term m + 2 (m + 1)(q − 1) .
6.5. LINEAR AND DECLINING-BALANCE DEPRECIATION METHODS 51
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. International tax laws generally provide investors
with a choice between two particular mathematical methods of calculating de-
preciation. 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 that for the difference

of remaining values for years adjacent one obtains
K0
Rn − Rn−1 = − =: d < 0. The underlying mathematical structure of the
N
straight line depreciation 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 = Rn−1 q , R0 ≡ K0 , n∈N. (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 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, n∈N. (6.42)

From Eq. (6.42), one may derive results concerning the following questions of a
quantitative nature:
(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 !
Rn 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 con-
densed 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. Zeh–Marschke (2010) [26]):
qn − 1
K n = K0 q n + R , q ∈ R>0 \{1} , n∈N. (6.45)
q−1
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 ,
1
Cf. our introduction in Ch. 2 of matrices as a particular class of mathematical objects.
53
54 CHAPTER 7. DIFFERENTIAL CALCULUS OF REAL-VALUED FUNCTIONS
• 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 ,
• 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)

∆x→0 ∆x→0
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 variable that find frequent application in the modelling of quantitative
problems in economic theory.
7.1. REAL-VALUED FUNCTIONS 55
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 + an−1 xn−1 + . . . + 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.2 Rational functions

Rational functions are constructed by forming the ratio of two polynomials of
degrees m resp. n, i.e.,
pm (x) am x m + . . . + a1 x + a0
y = f (x) = =
qn (x) bn x n + . . . + 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 polynomial part and a proper rational part. The roots of f always corre-
spond to those roots of the numerator polynomial pm (x) for which simultaneously
qn (x) 6= 0 applies. The roots of the denominator 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 increasing when α > 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 Euler’s number (according to the Swiss mathematician
Leonhard Euler, 1707–1783) 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) ob-
tains when the constant basis number is set to a = e. This yields
y = f (x) = loge (x) := ln(x) . (7.11)

7.2. DERIVATION OF DIFFERENTIABLE REAL-VALUED FUNCTIONS 57
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 func-

tions
The central theme of this chapter is the mathematical description of the local
variability of continuous 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 differen-
tiable at x ∈ D(f ), when for arbitrary ∆x ∈ R the limit
∆y f (x + ∆x) − f (x)
f 0 (x) := lim = lim (7.12)
∆x→0 ∆x ∆x→0 ∆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 re-
quires that the latter exhibits neither “jumps” nor “kinks,” i.e., that at (x, f (x))
the function is sufficiently “smooth.” The quantity f 0 (x) is referred to as the first
derivative of the (differentiable) function f at position x. It provides a quantita-
tive 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 0 (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 0 (x) · 1 units.
Alternative notation for the first derivative of f :
df (x)
f 0 (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, mathematiccian, astronomer and philosopher Sir Isaac
Newton (1643–1727) and the German philosopher, mathematician and physicist
Gottfried Wilhelm Leibniz (1646–1716).
Via the first derivative of a differentiable function f at an argument x0 ∈ D(f ),
i.e., f 0 (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 0 (x0 )(x − x0 ) . (7.13)
GDC: Local values f 0 (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 = 0 for c = constant ∈ R (constants)
2. (x)0 = 1 (linear function)
3. (xn )0 = nxn−1 for n ∈ N (natural power-law functions)
4. (xα )0 = αxα−1 for α ∈ R and x ∈ R>0 (general power-law functions)
5. (ax )0 = ln(a)ax for a ∈ R>0 \{1} (exponential functions)
6. (eax )0 = aeax for a ∈ R (natural exponential functions)

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

0
f 0 (x)g(x) − f (x)g 0 (x)

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

−1 0 1
7. (f (x)) = 0 , 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 eco-
nomic theory, and also to determine their local extremal values. The following
section provides an overview of such frequently occurring economic functions.
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 0 (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 0 (x) = xp0 (x) + p(x) (dim: CU/unit)

7. profit function G(x) := E(x) − K(x) (dim: CU)
8. marginal profit function G0 (x) := E 0 (x) − K 0 (x) = xp0 (x) + p(x) − K 0 (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 eco-
nomic 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 (1700–1782); 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 00 (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
Israeli–US-American experimental psychologists Daniel Kahneman and Amos
Tversky (1937–1996) in the context of their Prospect Theory (a pillar of Be-
havioural Economics), which was later awarded a Sveriges Riksbank Prize in
Economic Sciences in Memory of Alfred Nobel in 2002 (cf. Kahneman and Tver-
sky (1979) [18, p 279], and Kahneman (2011) [17, p 282f]). A possible represen-
tation of this function is given by the piecewise description

a log10 (1 + x)
 for x ∈ R≥0
v(x) = , (7.14)

−2a log10 (1 − x) for x ∈ R<0

with parameter a ∈ R>0 . Overcoming a conceptual problem of Bernoulli’s 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 61
7.4 Curve sketching

Before we turn to discuss applications of differential calculus to simple quan-
titative 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 0 (x) = 0, and the
!
sufficient condition f 00 (x) > 0 are satisfied simultaneously.
(ii) local maxima of f exist at all xE ∈ D(f ), for which the
!
!
sufficient condition f 00 (x) < 0 are satisfied simultaneously.
5. points of inflection: find all xW ∈ D(f ), for which the

!
!
sufficient condition f 000 (x) 6= 0 are satisfied simultaneously.
6. monotonous behaviour:
(i) f is monotonously increasing for all x ∈ D(f ) with f 0 (x) > 0

(ii) f is monotonously decreasing for all x ∈ D(f ) with f 0 (x) < 0
7. local curvature:
(i) f behaves left-handedly curved for x ∈ D(f ) with f 00 (x) > 0

(ii) f behaves right-handedly curved for x ∈ D(f ) with f 00 (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
Thünen
According to the law of diminishing returns, which was introduced to eco-
nomic theory by the French economist and statesman Anne Robert Jacques
Turgot (1727–1781) and also by the German economist Johann Heinrich von
Thünen (1783–1850), it is meaningful to model non-negative total cost func-
tions 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 regression analysis of agricultural quantitative–empirical 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
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.
7.5. ANALYTIC INVESTIGATIONS OF ECONOMIC FUNCTIONS 63
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 desig-
nated 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 ob-
tains K 00 (x) < 0 CU/unit2 , i.e., K(x) increases in a degressive fashion. For
values of x larger than xW , the opposite applies, K 00 (x) > 0 CU/unit2 , i.e.,
K(x) increases in a progressive fashion. The marginal costs, given by
K 0 (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 0 (x) , (7.19)
x
which follows by the quotient rule of differentiation from the necessary con-
dition for an extremum of the variable average costs,
0
(K(x) − Kf )0 · x − Kv (x) · 1

! Kv (x)
0= = ,
x x2
and the fact that Kf0 = 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 0 (xg1 )(x−xg1 ) = Kv (xg1 )+Kf + (x−xg1 ) = Kf + x.
xg 1 xg 1
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 0 (x) , (7.21)
x
which follows by the quotient rule of differentiation from the necessary con-
dition for an extremum of the average costs,
0
K 0 (x) · x − K(x) · 1

! K(x)
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 0 (x) K 0 (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 eco-

nomic theory as the minimum efficient scale (MES). From a business
economics perspective, at a level of physical output x = xg2 the (compared
to our remarks in the Introduction inverted) ratio “INPUT over OUTPUT,”
K(x)
i.e., , becomes most favourable. By respecting the property (7.21), the
x
equation for the tangent to K(x) in this point [cf. Eq. (7.13)] becomes
K(xg2 ) K(xg2 )
T (x) = K(xg2 ) + K 0 (xg2 )(x − xg2 ) = K(xg2 ) + (x − xg2 ) = x.
xg 2 xg 2
Its intercept with the K-axis is thus at 0 CU.
• phase IV (half-interval x > xg2 ):

In this phase K 0 (x)/K(x)/x > 1 obtains; the costs associated with the
production of one additional unit of a good, approximately the marginal
costs K 0 (x), now exceed the average costs, K(x)/x. This situation is consid-
ered 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 synony-
mous with the supplier’s 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 function 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 or-
der 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 determined 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 depen-
dence on a known amount x of units supplied/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 Thünen, 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
G0 (x) = E 0 (x) − K 0 (x) = xp0 (x) + p(x) − 3a3 x2 + 2a2 x + a1

(7.26)
00 00 00 00 0
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
!
G0 (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
!
G0 (x) < 0 CU/unit (sufficient condition) , (7.31)
• maximum profit
xM > 0 CU, as the unique solution to the conditions
!
G0 (x) = 0 CU/unit (necessary condition) (7.32)
and
!
G00 (x) < 0 CU/unit2 (sufficient condition) . (7.33)
At this point, we like to draw the reader’s attention to a special geometric prop-
erty of the quantitative model for profit that we just have outlined: at maximum
profit, the total revenue function 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 = G0 (x) = E 0 (x) − K 0 (x) ⇔ E 0 (x) = K 0 (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
Cournot’s point, which frequently arises in quantitative discussions in eco-
nomic theory; this is named after the French mathematician and economist
Antoine–Augustin Cournot (1801–1877). Cournot’s 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 monopolis-
tic market situation. Note that for this specific combination of optimal values the
Amoroso–Robinson formula applies, which was developed by the Italian math-
ematician and economist Luigi Amoroso (1886–1965) and the British economist
Joan Violet Robinson (1903–1983). This states that
K 0 (xM )
p(xM ) = , (7.35)
1 + εp (xM )
with K 0 (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 Amoroso–
Robinson formula is derived by evaluating the first derivative of E(x) at xM ,
so
Sec. 7.6
p0 (xM ) z}|{

0 0
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 0 (xM ) = K 0 (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, p0 (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 func-
tions 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 = 0 CU/unit2
!
and [G(x)/x]00 < 0 CU/unit3 . Respecting the quotient rule of differentiation
(cf. Sec. 7.2), the first condition yields
G0 (x)x − G(x)
= 0 GE/ME2 . (7.37)
x2
Since a quotient can only be zero when its numerator vanishes while its
denominator remains non-zero, it immediately follows that
G0 (x)
G0 (x)x − G(x) = 0 CU ⇒ x =1. (7.38)
G(x)
The task at hand now is to find a (unique) value of the level of physical out-
put 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 dis-
cussion in (i), the conditions for the existence of a local maximum amount
! !
to [E(x)/K(x)]0 = 0 × 1/unit and [E(x)/K(x)]00 < 0 × 1/unit2 . By the quotient
rule of differentiation (see Sec. 7.2), the first condition leads to
E 0 (x)K(x) − E(x)K 0 (x)
= 0 × 1/unit , (7.40)
K 2 (x)
i.e., for K(x) > 0 CU,
E 0 (x)K(x) − E(x)K 0 (x) = 0 CU2 /unit . (7.41)
By re-arranging and multiplication with x > 0 unit, this can be cast into the
particular form
E 0 (x) K 0 (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,
7.6. ELASTICITIES 69
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
∆x∈R ∆y∈R
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 de-

fined by the quotient
new value − old value
old value
under the prerequisite that “old value > 0” applies. It follows from this spe-
cific 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
• the associated resultant relative change of the function f : =
y
f (x + ∆x) − f (x)
.
f (x)
Now let us compare the order-of-magnitudes of the two relative changes just
∆x ∆y
envisaged, and . This is realised by considering the value of their quotient,
x y
“resultant relative change of f 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 of relative changes in the limit of increasingly smaller prescribed
∆x
relative changes → 0 ⇒ ∆x → 0 near some x > 0. One thus defines:
x
Def.: For a differentiable real-valued function f of one real variable x, the dimen-
sionless (i.e., units-independent) quantity
∆y f (x + ∆x) − f (x)
y f (x) f 0 (x)
εf (x) := lim = lim =x (7.44)
∆x→0 ∆x ∆x→0 ∆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 pre-
scribed infinitesimally small relative change of its argument x, starting from
some positive initial value x > 0. As such it constitutes a measure for the rel-
ative 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 elastic-
ity ε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 differenti-
ation. 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 0 (x)
= = x dx = x .
d ln(x) dx f (x) f (x)
x
The logarithmic representation of the elasticity of a differentiable function f im-
mediately explains why, at the beginning, we confined our considerations to pos-
itive differentiable functions of positive 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)

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.
7.6. ELASTICITIES 71
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 general 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 clas-
sifies general power-law functions as scale-invariant. When scale-invariance
obtains, dimensionless ratios, i.e., quotients of variables of the 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 Sci-
ences for modelling uncertainty of economic agents in decision-making pro-
cesses, 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 α gen-
eral 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: begin-
ning at x0 > 0, for small prescribed percentage changes of the argument x in
∆x
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 reader’s 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 ex-
hibits 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 func-
tion 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 elas-
ticities ε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 func-
tion 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 0 (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 = 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 0 (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,
73
74 CHAPTER 8. INTEGRAL CALCULUS OF REAL-VALUED FUNCTIONS
• f (x) — the integrand,
• dx — the differential, and, lastly,
• 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)
x2
Z
2. x dx = +c (linear functions)
2
xn+1
Z
3. xn dx = + c for n ∈ N (natural power-law functions)
n+1
xα+1
Z
4. xα dx = + c for α ∈ R\{−1} and x ∈ R>0 (general power-law
α+1
functions)
ax
Z
5. ax dx = + c for a ∈ R>0 \{1} (exponential functions)
ln(a)
eax
Z
6. eax dx = + c for a ∈ R\{0} (natural exponential functions)
a
7. x−1 dx = ln |x| + c for x ∈ R\{0}.

R
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)
2. f (x)g 0 (x) dx = f (x)g(x) − f 0 (x)g(x) dx

R R
(integration by parts)
u=g(x) and du=g 0 (x)dx

f (g(x))g 0 (x) dx
R R
3. = f (u) du = F (g(x)) + c (substitution
z}|{
method)
Z 0
f (x)
4. dx = ln |f (x)| + c for f (x) 6= 0 (logarithmic integration).
f (x)
8.2. DEFINITE INTEGRALS 75
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:
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
a a c
interval).
Remark: The main qualitative difference between an (i) indefinite integral and a
(ii) definite integral 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 of a definiteR 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)
a→0 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 assum-
ing a state of economic equilibrium of the market, quantitatively expressed by
the condition
A(pM ) = N (pM ) . (8.7)
Geometrically, this condition defines common points of intersection for the func-
tions A(p) and N (p) (when they exist).
GDC: Common points of intersection for stored functions f and g can be easily
determined interactively 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 according to strategy 1, this amount is referred to in the
economic literature as consumer surplus.
1
This is a strong mathematical assumption aimed at facilitating the actual calculation to
follow.
8.3. APPLICATIONS IN ECONOMIC THEORY 77
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.
2
See previous footnote.
Appendix A
Glossary of technical terms (GB – D)
A
absolute change: absolute Änderung
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: Stückkosten
average profit: Durchschnittsgewinn, Gewinn pro Stück
B
backward substitution: rückwertige Substitution
balance equation: Bilanzgleichung
basis: Basis
basis solution: Basislösung
basis variable: Basisvariable
Behavioural Economics: Verhaltensökonomik
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, Verknüpfung
conservation law: Erhaltungssatz
79
80 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
Cournot’s point: Cournotscher Punkt
curve sketching: Kurvendiskussion
D
decision-making: Entscheidungsfindung
declining-balance depreciation method: geometrisch–degressive Abschreibung
definite integral: bestimmtes Integral
demand function: Nachfragefunktion
dependent variable: abhängige 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, unbeschr”ankt
domain: Definitionsbereich
E
economic agent: Wirtschaftstreibende(r) (meistens ein homo oeconomicus)
economic efficiency: Wirtschaftlichkeit
economic equilibrium: ökonomisches Gleichgewicht
economic principle: ökonomisches Prinzip
economic theory: Wirtschaftstheorie
Econophysics: Ökonophysik
eigenvalue: Eigenwert
eigenvector: Eigenvektor
elastic: elastisch
elasticity: Elastizität
element: Element
end of profitable zone: Gewinngrenze
endogenous: endogen
equilibrium price: Marktpreis
equivalence transformation: Äquivalenztransformation
81
exogenous: exogen
exponential function: Exponentialfunktion
extrapolation: Extrapolation, über bekannten Gütigkeitsbereich hinaus verallge-
meinern
F
feasible region: zulässiger Bereich
final capital: Endkapital
fixed costs: Fixkosten
forecasting: Vorhersagen erstellen
function: Funktion
G
Gaußian elimination: Gauß’scher Algorithmus
GDC: GTR, grafikfähiger Taschenrechner
geometrical mean: geometrischer Mittelwert
geometrical sequence: geometrische Zahlenfolge
geometrical series: geometrische Reihe
growth function: Wachstumsfunktion
H
I
identity: Identität
image vector: Absolutgliedvektor
indefinite integral: unbestimmtes Integral
independent variable: unabhängige 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: Länge
level of physical output: Ausbringungsmenge
limits of integration: Integrationsgrenzen
linear combination: Linearkombination
linearisation: Linearisierung
linear programming: lineare Optimierung
local rate of change: lokale Änderungsrate
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: Nichtnegativitätsbedingungen
non-linear functional relationship: nichtlineare Funktionalbeziehung
O
objective function: Zielfunktion
one-to-one and onto: eineindeutig
optimal solution: optimalen Lösung
optimal value: optimaler Wert
optimisation: Optimierung
order-of-magnitude: Größenordnung
orthogonal: orthogonal
over-determined: überbestimmt
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, Singularität
polynomial division: Polynomdivision
83
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: regulär
relative change: relative Änderung
remaining debt: Restschuld
remaining resources: Restkapazitäten
remaining value: Restwert
rescaling: Skalierung
resources: Rohstoffe
resource consumption matrix: Rohstoffverbrauchsmatrix
restrictions: Restriktionen
root: Nullstelle
row: Reihe
row vector: Zeilenvektor
S
saturation quantity: Sättigungsmenge
scale: Skala, Größenordnung
scale-invariant: skaleninvariant
simplex: Simplex, konvexer Polyeder
simplex tableau: Simplextabelle
singular: singulär
sink: Senke
slack variable: Schlupfvariable
source: Quelle
stationary: stationär, 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: Stückpreis
unit vector: Einheitsvektor
utility function: Nutzenfunktion
V
value chain: Wertschöpfungskette
variability: Änderungsverhalten, Variabilität
variable average costs: variable Stückkosten
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) ISBN–13:
9783827417589
[2] A Asano (2013) An Introduction to Mathematics for Economics (Cambridge:

Cambridge University Press) ISBN–13: 9780521189460
[3] C Bauer, M Clausen, A Kerber and H Meier–Reinhold (2008) Mathematik für

Wirtschaftswissenschaftler 5th Edition (Stuttgart: Schäffer–Poeschel) ISBN–
13: 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 23–36
[5] J–P Bouchaud and M Potters (2003) Theory of Financial Risk and Deriva-
tive Pricing 2nd Edition (Cambridge: Cambridge University Press) ISBN–13:
9780521741866
[6] K Bosch (2003) Mathematik für Wirtschaftswissenschaftler 14th Edition

(München: Oldenbourg) ISBN–10: 3486273094
[7] I N Bronstein, K A Semedjajew, G Musiol and H Mühlig (2005) Taschen-

buch der Mathematik 6th Edition (Frankfurt (Main): Harri Deutsch) ISBN–10:
3817120060
[8] G B Dantzig (1949) Programming of interdependent activities: II mathemat-

ical model Econometrica 17 200–211
[9] G B Dantzig (1955) Linear programming under uncertainty Management

Science 1 197–206
[10] E T Dowling (1990) Schaum’s Outline of Calculus for Business, Economics,

and the Social Sciences (New York: McGraw–Hill) ISBN–13: 9780070176737
[11] E T Dowling (2009) Schaum’s Outline of Mathematical Methods for Business

and Economics (New York: McGraw–Hill) ISBN–13: 9780071635325
85
86 BIBLIOGRAPHY
[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
[14] J F Hair jr, W C Black, B J Babin and R E Anderson (2010) Multivari-
ate Data Analysis 7th Edition (Upper Saddle River (NJ): Pearson) ISBN–13:
9780135153093
[15] F Herrmann (2003) Physik III: Thermodynamik — Skripten zur Experimental-
physik Universität Karlsruhe (TH) URL (cited on August 30, 2015):
www.physikdidaktik.uni-karlsruhe.de/skripten/thermod.pdf
[16] J Hülsmann, W Gamerith, U Leopold–Wildburger and W Steindl (2005)
Einführung in die Wirtschaftsmathematik 4th Edition (Berlin: Springer) ISBN–
10: 3540244093
[17] D Kahneman (2011) Thinking, Fast and Slow (London: Penguin) ISBN–13:
9780141033570
[18] D Kahneman and A Tversky (1979) Prospect theory: an analysis of decision
under risk Econometrica 47 263–292
[19] S J Leon (2009) Linear Algebra with Applications 8th Edition (London: Pear-
son) ISBN–13: 9780136009290
[20] W W Leontief (1936) Quantitative input and output relations in the eco-
nomic systems of the United States The Review of Economics and Statistics
18 105–125
[21] R Penrose (2004) The Road to Reality — A Complete Guide to the Laws of the
Universe 1st Edition (London: Jonathan Cape) ISBN–10: 0224044478
[22] H Rinne (2008) Taschenbuch der Statistik 4th Edition (Frankfurt Main): Harri
Deutsch) ISBN–13: 9783817118274
[23] S Singh (1997) Fermat’s Last Theorem (London: Fouth Estate) ISBN–13:
9781857025217
[24] L Smith (2007) Chaos — A Very Short Introduction (Oxford: Oxford University
Press) ISBN–13: 9780192853783
[25] N N Taleb (2007) The Black Swan — The Impact of the Highly Improbable
(London: Penguin) ISBN–13: 9780141034591
[26] A Zeh–Marschke (2010) Finanzmathematik — Zins und Zinseszins, Tilgungs-
und Rentenrechnung, Investitionsrechnung, Abschreibungen private lecture
notes

Business Maths

Uploaded by

Copyright:

Available Formats

Business Maths

Uploaded by

Document Information

Original Description:

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

Business Maths

Uploaded by

Copyright:

Available Formats

A N I NTRODUCTION TO

0.1.1 IMQM: Introduction to Management and its Quantitative Methods

H ENK VAN E LST

August 30, 2015

Fakultät I: Betriebswirtschaft und Management

E–Print: arXiv:1509.04333v2 [q-fin.GN]

c 2009–2015 Karlshochschule International University and Henk van Elst

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 7

3 Systems of linear algebraic equations 19

4 Leontief’s input–output matrix model 27

6 Elementary financial mathematics 41

7 Differential calculus of real-valued functions 53

8 Integral calculus of real-valued functions 73

A Glossary of technical terms (GB – D) 79

The qualification objectives shall be reached by an integrative approach.

Vector algebra in Euclidian space Rn

1.1 Basic concepts

to constitute related zero vectors.

• 0 < |λ| < 1 — shrinking of the length of a

The notion of the length of a vector a will be made precise shortly.

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)

In conclusion of this section, we remark that every set of mathematical objects

1.2 Dimension and basis of Rn

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

The rescaling factors βi ∈ R of the ai ∈ Rn are called the components of vector

Remark: The n unit vectors

constitute the so-called canonical basis of Euclidian space Rn . With respect to

1.3 Euclidian scalar product

defines a mapping f : R1×n × Rn×1 → R from the product set of n-component

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| = |λ||a| for λ ∈ R;

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

Remark: The inverse cosine function3 cos−1 (. . .) is available on every standard

In this chapter, we introduce a second class of mathematical objects that are

2.1 Matrices as linear mappings

a11 a12 ... a1j ... a1n

• 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

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 ∈ Rn×m .

• When AT = A, one refers to A as a symmetric matrix.

• When AT = −A, one refers to A as an antisymmetric matrix.

For A, B ∈ Rm×n , the sum is given by

A + B =: C : aij + bij =: cij , (2.8)

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

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