Analysis for computer scientists foundations methods and algorithms Second Edition Oberguggenberger 2024 scribd download
Analysis for computer scientists foundations methods and algorithms Second Edition Oberguggenberger 2024 scribd download
com
https://textbookfull.com/product/analysis-for-computer-
scientists-foundations-methods-and-algorithms-second-
edition-oberguggenberger/
OR CLICK BUTTON
DOWNLOAD NOW
https://textbookfull.com/product/analysis-for-computer-scientists-
foundations-methods-and-algorithms-michael-oberguggenberger/
textboxfull.com
https://textbookfull.com/product/numerical-methods-for-engineers-and-
scientists-using-matlab-second-edition-esfandiari/
textboxfull.com
https://textbookfull.com/product/machine-learning-refined-foundations-
algorithms-and-applications-second-edition-borhani/
textboxfull.com
https://textbookfull.com/product/numerical-algorithms-methods-for-
computer-vision-machine-learning-and-graphics-1st-edition-justin-
solomon/
textboxfull.com
Machine Learning Refined: Foundations, Algorithms, and
Applications Second Edition Jeremy Watt
https://textbookfull.com/product/machine-learning-refined-foundations-
algorithms-and-applications-second-edition-jeremy-watt/
textboxfull.com
https://textbookfull.com/product/quantum-computing-for-computer-
scientists-mannucci/
textboxfull.com
https://textbookfull.com/product/introduction-to-experimental-
biophysics-second-edition-biological-methods-for-physical-scientists-
jay-l-nadeau/
textboxfull.com
https://textbookfull.com/product/foundations-of-radio-for-scientists-
and-technologists-christopher-john-coleman/
textboxfull.com
https://textbookfull.com/product/essential-algorithms-a-practical-
approach-to-computer-algorithms-using-python-and-c-second-edition-rod-
stephens/
textboxfull.com
Undergraduate Topics in Computer Science
Michael Oberguggenberger
Alexander Ostermann
Analysis for
Computer
Scientists
Foundations, Methods, and
Algorithms
Second Edition
Undergraduate Topics in Computer
Science
Series editor
Ian Mackie
Advisory Board
Samson Abramsky, University of Oxford, Oxford, UK
Chris Hankin, Imperial College London, London, UK
Mike Hinchey, University of Limerick, Limerick, Ireland
Dexter C. Kozen, Cornell University, Ithaca, USA
Andrew Pitts, University of Cambridge, Cambridge, UK
Hanne Riis Nielson, Technical University of Denmark, Kongens Lyngby, Denmark
Steven S. Skiena, Stony Brook University, Stony Brook, USA
Iain Stewart, University of Durham, Durham, UK
Undergraduate Topics in Computer Science (UTiCS) delivers high-quality
instructional content for undergraduates studying in all areas of computing and
information science. From core foundational and theoretical material to final-year
topics and applications, UTiCS books take a fresh, concise, and modern approach
and are ideal for self-study or for a one- or two-semester course. The texts are all
authored by established experts in their fields, reviewed by an international advisory
board, and contain numerous examples and problems. Many include fully worked
solutions.
123
Michael Oberguggenberger Alexander Ostermann
University of Innsbruck University of Innsbruck
Innsbruck Innsbruck
Austria Austria
This Springer imprint is published by the registered company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Preface to the Second Edition
We are happy that Springer Verlag asked us to prepare the second edition of our
textbook Analysis for Computer Scientists. We are still convinced that the algo-
rithmic approach developed in the first edition is an appropriate concept for pre-
senting the subject of analysis. Accordingly, there was no need to make larger
changes.
However, we took the opportunity to add and update some material. In partic-
ular, we added hyperbolic functions and gave some more details on curves and
surfaces in space. Two new sections have been added: One on second-order dif-
ferential equations and one on the pendulum equation. Moreover, the exercise
sections have been extended considerably. Statistical data have been updated where
appropriate.
Due to the essential importance of the MATLAB programs for our concept, we
have decided to provide these programs additionally in Python for the users’
convenience.
We thank the editors of Springer, especially Simon Rees and Wayne Wheeler,
for their support during the preparation of the second edition.
v
Preface to the First Edition
vii
viii Preface to the First Edition
indispensable tool for modelling and numerical analysis.) For the sake of read-
ability, proofs are only detailed in the main text if they introduce essential ideas and
contribute to the understanding of the concepts. To continue with the example
above, the integral representation of the remainder term of the Taylor expansion is
derived by integration by parts. In contrast, Lagrange’s form of the remainder term,
which requires the mean value theorem of integration, is only mentioned. Never-
theless we have put effort into ensuring a self-contained presentation. We assign a
high value to geometric intuition, which is reflected in a large number of
illustrations.
Due to the terse presentation it was possible to cover the whole spectrum from
foundations to interesting applications of analysis (again selected from the view-
point of computer science), such as fractals, L-systems, curves and surfaces, linear
regression, differential equations and dynamical systems. These topics give suffi-
cient opportunity to enter various aspects of mathematical modelling.
The present book is a translation of the original German version that appeared in
2005 (with the second edition in 2009). We have kept the structure of the German
text, but took the opportunity to improve the presentation at various places.
The contents of the book are as follows. Chapters 1–8, 10–12 and 14–17 are
devoted to the basic concepts of analysis, and Chapters 9, 13 and 18–21 are ded-
icated to important applications and more advanced topics. The Appendices A and
B collect some tools from vector and matrix algebra, and Appendix C supplies
further details which were deliberately omitted in the main text. The employed
software, which is an integral part of our concept, is summarised in Appendix D.
Each chapter is preceded by a brief introduction for orientation. The text is enriched
by computer experiments which should encourage the reader to actively acquire the
subject matter. Finally, every chapter has exercises, half of which are to be solved
with the help of computer programs. The book can be used from the first semester
on as the main textbook for a course, as a complementary text or for self-study.
We thank Elisabeth Bradley for her help in the translation of the text. Further, we
thank the editors of Springer, especially Simon Rees and Wayne Wheeler, for their
support and advice during the preparation of the English text.
1 Numbers. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1 The Real Numbers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Order Relation and Arithmetic on R . . . . . . . . . . . . . . . . . . . . . . 5
1.3 Machine Numbers. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4 Rounding . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.5 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2 Real-Valued Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.1 Basic Notions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2 Some Elementary Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3 Trigonometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.1 Trigonometric Functions at the Triangle . . . . . . . . . . . . . . . . . . . 27
3.2 Extension of the Trigonometric Functions to R . . . . . . . . . . . . . 31
3.3 Cyclometric Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.4 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4 Complex Numbers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.1 The Notion of Complex Numbers. . . . . . . . . . . . . . . . . . . . . . . . 39
4.2 The Complex Exponential Function . . . . . . . . . . . . . . . . . . . . . . 42
4.3 Mapping Properties of Complex Functions . . . . . . . . . . . . . . . . . 44
4.4 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
5 Sequences and Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
5.1 The Notion of an Infinite Sequence . . . . . . . . . . . . . . . . . . . . . . 49
5.2 The Completeness of the Set of Real Numbers . . . . . . . . . . . . . 55
5.3 Infinite Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
5.4 Supplement: Accumulation Points of Sequences . . . . . . . . . . . . . 62
5.5 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
6 Limits and Continuity of Functions . . . . . . . . . . . . . . . . . . . . . . . . . . 69
6.1 The Notion of Continuity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
6.2 Trigonometric Limits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
ix
x Contents
14 Curves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185
14.1 Parametrised Curves in the Plane . . . . . . . . . . . . . . . . . . . . . . . . 185
14.2 Arc Length and Curvature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 193
14.3 Plane Curves in Polar Coordinates . . . . . . . . . . . . . . . . . . . . . . . 200
14.4 Parametrised Space Curves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 202
14.5 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 204
15 Scalar-Valued Functions of Two Variables . . . . . . . . . . . . . . . . . . . . 209
15.1 Graph and Partial Mappings . . . . . . . . . . . . . . . . . . . . . . . . . . . . 209
15.2 Continuity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 211
15.3 Partial Derivatives. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 212
15.4 The Fréchet Derivative . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 216
15.5 Directional Derivative and Gradient . . . . . . . . . . . . . . . . . . . . . . 221
15.6 The Taylor Formula in Two Variables . . . . . . . . . . . . . . . . . . . . 223
15.7 Local Maxima and Minima. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 224
15.8 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 228
16 Vector-Valued Functions of Two Variables . . . . . . . . . . . . . . . . . . . . 231
16.1 Vector Fields and the Jacobian . . . . . . . . . . . . . . . . . . . . . . . . . . 231
16.2 Newton’s Method in Two Variables . . . . . . . . . . . . . . . . . . . . . . 233
16.3 Parametric Surfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 236
16.4 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 238
17 Integration of Functions of Two Variables . . . . . . . . . . . . . . . . . . . . 241
17.1 Double Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 241
17.2 Applications of the Double Integral . . . . . . . . . . . . . . . . . . . . . . 247
17.3 The Transformation Formula . . . . . . . . . . . . . . . . . . . . . . . . . . . 249
17.4 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 253
18 Linear Regression . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 255
18.1 Simple Linear Regression . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 255
18.2 Rudiments of the Analysis of Variance . . . . . . . . . . . . . . . . . . . 261
18.3 Multiple Linear Regression . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 265
18.4 Model Fitting and Variable Selection . . . . . . . . . . . . . . . . . . . . . 267
18.5 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 271
19 Differential Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 275
19.1 Initial Value Problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 275
19.2 First-Order Linear Differential Equations . . . . . . . . . . . . . . . . . . 278
19.3 Existence and Uniqueness of the Solution . . . . . . . . . . . . . . . . . 283
19.4 Method of Power Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 286
19.5 Qualitative Theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 288
19.6 Second-Order Problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 290
19.7 Exercises. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 294
xii Contents
The commonly known rational numbers (fractions) are not sufficient for a rigor-
ous foundation of mathematical analysis. The historical development shows that for
issues concerning analysis, the rational numbers have to be extended to the real num-
bers. For clarity we introduce the real numbers as decimal numbers with an infinite
number of decimal places. We illustrate exemplarily how the rules of calculation and
the order relation extend from the rational to the real numbers in a natural way.
A further section is dedicated to floating point numbers, which are implemented
in most programming languages as approximations to the real numbers. In particular,
we will discuss optimal rounding and in connection with this the relative machine
accuracy.
Two rational numbers nk and m are equal if and only if km = n. Further an integer
k ∈ Z can be identified with the fraction k1 ∈ Q. Consequently, the inclusions N ⊂
Z ⊂ Q are true.
Definition 1.1 Two sets M and N have the same cardinality if there exists a bijective
mapping between these sets. A set M is called countably infinite if it has the same
cardinality as N.
The sets N, Z and Q have the same cardinality and in this sense are equally large.
All three sets have an infinite number of elements which can be enumerated. Each
enumeration represents a bijective mapping to N. The countability of Z can be seen
from the representation Z = {0, 1, −1, 2, −2, 3, −3, . . .}. To prove the countability
of Q, Cantor’s2 diagonal method is being used:
1
1 → 2
1
3
1 → 4
1 ...
1
2
2
2
3
2
4
2 ...
↓
1
3
2
3
3
3
4
3 ...
1
4
2
4
3
4
4
4 ...
.. .. .. ..
. . . .
The enumeration is carried out in direction of the arrows, where each rational number
is only counted at its first appearance. In this way the countability of all positive
rational number (and therefore all rational numbers) is proven.
To visualise the rational numbers we use a line, which can be pictured as an
infinitely long ruler, on which an arbitrary point is labelled as zero. The integers are
marked equidistantly starting from zero. Likewise each rational number is allocated
a specific place on the real line according to its size, see Fig. 1.1.
However, the real line also contains points which do not correspond to rational
numbers. (We say that Q is not complete.) For instance, the length of the diagonal d
in the unit square (see Fig. 1.2) can be measured
√ with a ruler. Yet, the Pythagoreans
already knew that d 2 = 2, but that d = 2 is not a rational number.
−2 −1 − 12 0 1
3
1
2 1 a 2
1 We will rarely use the term mapping in such generality. The special case of real-valued functions,
which is important for us, will be discussed thoroughly in Chap. 2.
2 G. Cantor, 1845–1918.
1.1 The Real Numbers 3
√
Proposition 1.2 2∈
/ Q.
√
Proof
√ This statement is proven indirectly. Assume√ that 2 were rational. Then
2 can be represented as a reduced fraction 2 = nk ∈ Q. Squaring this equa-
tion gives k 2 = 2n 2 and thus k 2 would be an even number. This is only possible if
k itself is an even number, so k = 2l. If we substitute this into the above we obtain
4l 2 = 2n 2 which simplifies to 2l 2 = n 2 . Consequently n would also be even which
is in contradiction to the initial assumption that the fraction nk was reduced.
√
As it is generally known, 2 is the unique positive root of the polynomial x 2 − 2.
The naive supposition that all non-rational numbers are roots of polynomials with
integer coefficients turns out to be incorrect. There are other non-rational numbers
(so-called transcendental numbers) which cannot be represented in this way. For
example, the ratio of a circle’s circumference to its diameter
π = 3.141592653589793... ∈
/Q
is transcendental, but can be represented on the real line as half the circumference
of the circle with radius 1 (e.g. through unwinding).
In the following we will take up a pragmatic point of view and construct the
missing numbers as decimals.
Definition 1.3 A finite decimal number x with l decimal places has the form
x = ± d0 .d1 d2 d3 . . . dl
Each nonzero decimal with a finite number of decimal places can be written as a
periodic decimal (with an infinite number of decimal places). To this end one dimin-
ishes the last nonzero digit by one and then fills the remaining infinitely many decimal
places with the digit 9. For example, the fraction − 50 17
= −0.34 = −0.3399999...
becomes periodic after the third decimal place. In this way Q can be considered as
the set of all decimals which turn periodic from a certain number of decimal places
onwards.
Definition 1.6 The set of real numbers R consists of all decimals of the form
± d0 .d1 d2 d3 ...
with d0 ∈ N0 and digits di ∈ {0, ..., 9}, i.e. decimals with an infinite number of
decimal places. The set R \ Q is called the set of irrational numbers.
are irrational. There are much more irrational than rational numbers, as is shown by
the following proposition.
Proposition 1.7 The set R is not countable and has therefore higher cardinality
than Q.
Proof This statement is proven indirectly. Assume the real numbers between 0 and
1 to be countable and tabulate them:
Then x = 0.d1 d2 d3 d4 ... is not included in the above list which is a contradiction to
the initial assumption of countability.
1.1 The Real Numbers 5
30
1;24,51,10
5
42;25,3
Fig. 1.3 Babylonian cuneiform inscription YBC 7289 (Yale Babylonian Collection, with authori-
sation) from 1900 before our time with a translation of the inscription
√ according to [1]. It represents
a square with side length 30 and diagonals 42; 25, 35. The ratio is 2 ≈ 1; 24, 51, 10
√ 24 51 10
2 ≈ 1; 24, 51, 10 = 1 + + 2 + 3 = 1.41421296... ,
60 60 60
see Fig. 1.3. The somewhat unfamiliar notation is due to the fact that the Babylonians
worked in the sexagesimal system with base 60.
In the following we write real numbers (uniquely) as decimals with an infinite number
of decimal places, for example, we write 0.2999... instead of 0.3.
Definition 1.8 (Order relation) Let a = a0 .a1 a2 ... and b = b0 .b1 b2 ... be non-
negative real numbers in decimal form, i.e. a0 , b0 ∈ N0 .
(a) One says that a is less than or equal to b (and writes a ≤ b), if a = b or if there
is an index j ∈ N0 such that a j < b j and ai = bi for i = 0, . . . , j − 1.
(b) Furthermore one stipulates that always −a ≤ b and sets −a ≤ −b whenever
b ≤ a.
The relation ≤ obviously has the following properties. For all a, b, c ∈ R it holds
that
a ≤ a (reflexivity),
a ≤ b and b ≤ c ⇒ a ≤ c (transitivity),
a ≤ b and b ≤ a ⇒ a = b (antisymmetry).
In case of a ≤ b and a = b one writes a < b and calls a less than b. Furthermore
one defines a ≥ b, if b ≤ a (in words: a greater than or equal to b), and a > b, if
b < a (in words: a greater than b).
Addition and multiplication can be carried over from Q to R in a similar way.
Graphically one uses the fact that each real number corresponds to a segment on
the real line. One thus defines the addition of real numbers as the addition of the
respective segments.
A rigorous and at the same time algorithmic definition of the addition starts from
the observation that real numbers can be approximated by rational numbers to any
degree of accuracy. Let a = a0 .a1 a2 ... and b = b0 .b1 b2 ... be two non-negative real
numbers. By cutting them off after k decimal places we obtain two rational approxi-
mations a (k) = a0 .a1 a2 ...ak ≈ a and b(k) = b0 .b1 b2 ...bk ≈ b. Then a (k) + b(k) is a
monotonically increasing sequence of approximations to the yet to be defined num-
ber a + b. This allows one to define a + b as supremum of these approximations.
To justify this approach rigorously we refer to Chap. 5. The multiplication of real
numbers is defined in the same way. It turns out that the real numbers with addition
and multiplication (R, +, ·) are a field. Therefore the usual rules of calculation apply,
e.g., the distributive law
(a + b)c = ac + bc.
The following proposition recapitulates some of the important rules for ≤. The
statements can easily be verified with the help of the real line.
a≤b ⇒ a + c ≤ b + c,
a ≤ b and c ≥ 0 ⇒ ac ≤ bc,
a ≤ b and c ≤ 0 ⇒ ac ≥ bc.
Note that a < b does not imply a 2 < b2 . For example −2 < 1, but nonetheless
4 > 1. However, for a, b ≥ 0 it always holds that a < b ⇔ a 2 < b2 .
a b c d e f
Fig. 1.4 The intervals (a, b), [c, d] and (e, f ] on the real line
∀a ∈ R : −∞ < a < ∞.
[a, ∞) = {x ∈ R ; x ≥ a}
(−∞, b) = {x ∈ R ; x < b}
and furthermore (−∞, ∞) = R. Note that −∞ and ∞ are only symbols and not
real numbers.
Example 1.12 Find all x ∈ R satisfying −3x − 2 ≤ 5 < −3x + 4. In this example
we have the following two inequalities
7
−3x ≤ 7 ⇔ x ≥− .
3
This is the first constraint for x. The second inequality states
1
3x < −1 ⇔ x <−
3
and poses a second constraint for x. The solution to the original problem must fulfil
both constraints. Therefore the solution set is
7 1 7 1
S = x ∈ R; − ≤ x < − = − ,− .
3 3 3 3
8 1 Numbers
(x − 1)2 = x 2 − 2x + 1 ≥ 4.
x −1≥2 or x − 1 ≤ −2.
The real numbers can be realised only partially on a computer. In exact arithmetic,
√ for example in maple, real numbers are treated as symbolic expressions, e.g.
like
2 = RootOf(_Zˆ2-2). With the help of the command evalf they can be
evaluated, exact to many decimal places.
The floating point numbers that are usually employed in programming languages
as substitutes for the real numbers have a fixed relative accuracy, e.g. double precision
with 52 bit mantissa. The arithmetic rules of R are not valid for these machine
numbers, e.g.
1 + 10−20 = 1
in double precision. Floating point numbers are standardised by the Institute of
Electrical and Electronics Engineers IEEE 754-1985 and by the International Elec-
trotechnical Commission IEC 559:1989. In the following we give a short outline of
these machine numbers. Further information can be found in [20].
One distinguishes between single and double format. The single format (single
precision) requires 32-bit storage space
V e M
1 8 23
V e M
1 11 52
Here, V ∈ {0, 1} denotes the sign, emin ≤ e ≤ emax is the exponent (a signed integer)
and M is the mantissa of length p
M = d1 2−1 + d2 2−2 + . . . + d p 2− p ∼
= d1 d2 . . . d p , d j ∈ {0, 1}.
1.3 Machine Numbers 9
p
x = (−1) 2 V e
d j 2− j .
j=1
To simplify matters we will only describe the key features of floating point numbers.
For the subtleties of the IEEE-IEC standard, we refer to [20].
In our representation the following range applies for the exponents:
emin emax
single precision −125 128
double precision −1021 1024
With M = Mmax and e = emax one obtains the largest floating point number
xmax = 1 − 2− p 2emax ,
whereas M = Mmin and e = emin gives the smallest positive (normalised) floating
point number
xmin = 2emin −1 .
The floating point numbers are not evenly distributed on the real line, but their relative
density is nearly constant, see Fig. 1.5.
In the IEEE standard the following approximate values apply:
xmin xmax
single precision 1.18 · 10−38 3.40 · 1038
double precision 2.23 · 10−308 1.80 · 10308
±INF ... ±∞
NaN ... not a number, e.g. for zero divided by zero.
In general, one can continue calculating with these symbols without program termi-
nation.
10 1 Numbers
1.4 Rounding
u= 0 e b1 . . . b p
and
v=u+ 0 e 00 . . . 01 = u + 0 e − ( p − 1) 10 . . . 00
1
|rd(x) − x| ≤ (v − u) = 2e− p−1
2
holds for the optimal rounding rd(x) of x. With this estimate one can determine the
relative error of the rounding. Due to a1 ≤ 2 it holds that
Proposition 1.15 Let x ∈ R with xmin ≤ |x| ≤ xmax . Then there exists ε ∈ R with
Proof We define
rd(x) − x
ε= .
x
According to the calculation above, we have |ε| ≤ eps.
1= 0 1 100 . . . 00 , z= 0 1 000 . . . 01 = 2 · 2− p .
Thus z = 2 eps. The number z can be determined experimentally and therefore eps
as well. (Note that the number z is called eps in MATLAB.)
1.4 Rounding 11
1.5 Exercises
√
1. Show that 3 is irrational.
2. Prove the triangle inequality
|a + b| ≤ |a| + |b|
for all a, b ∈ R.
Hint. Distinguish the cases where a and b have either the same or different signs.
3. Sketch the following subsets of the real line:
More generally, sketch the set Ur (a) = {x : |x − a| < r } (for a ∈ R, r > 0).
Convince yourself that Ur (a) is the set of points of distance less than r to the
point a.
4. Solve the following inequalities by hand as well as with maple (using solve).
State the solution set in interval notation.
1
(a) 4x 2 ≤ 8x + 1, (b) > 3 + x,
3−x
1+x
(c) 2 − x 2 ≥ x 2, (d) > 1,
1−x
(e) x 2 < 6 + x, (f) |x| − x ≥ 1,
(g) |1 − x 2 | ≤ 2x + 2, (h) 4x 2 − 13x + 4 < 1.
20
8(x − 2) ≥ + 3(x − 7).
x +1
(a) x = y; (b) y < x; (c) y > x; (d) y > |x|; (e) |y| > |x|.
Hint. Consult Sects. A.1 and A.6 for basic plane geometry.
12 1 Numbers
7. Compute the binary representation of the floating point number x = 0.1 in single
precision IEEE arithmetic.
8. Experimentally determine the relative machine accuracy eps.
Hint. Write a computer program in your programming language of choice which
calculates the smallest machine number z such that 1 + z > 1.
Real-Valued Functions
2
The notion of a function is the mathematical way of formalising the idea that one
or more independent quantities are assigned to one or more dependent quantities.
Functions in general and their investigation are at the core of analysis. They help to
model dependencies of variable quantities, from simple planar graphs, curves and
surfaces in space to solutions of differential equations or the algorithmic construction
of fractals. One the one hand, this chapter serves to introduce the basic concepts. On
the other hand, the most important examples of real-valued, elementary functions
are discussed in an informal way. These include the power functions, the exponential
functions and their inverses. Trigonometric functions will be discussed in Chap. 3,
complex-valued functions in Chap. 4.
Experiment 2.1 Study the mapping y = x 2 with the help of MATLAB. First choose
the region D in which the x-values should vary, for instance D = {x ∈ R : −1 ≤
x ≤ 1}. The command
x = −1 : 0.01 : 1;
produces a list of x-values, the row vector
Using
y = x.ˆ2;
a row vector of the same length of corresponding y-values is generated. Finally
plot(x,y) plots the points (x1 , y1 ), . . . , (xn , yn ) in the coordinate plane and con-
nects them with line segments. The result can be seen in Fig. 2.1.
In the general mathematical framework we do not just want to assign finite lists
of values. In many areas of mathematics functions defined on arbitrary sets are
needed. For the general set-theoretic notion of a function we refer to the literature,
e.g. [3, Chap. 0.2]. This section is dedicated to real-valued functions, which are
central in analysis.
Definition 2.2 A real-valued function f with domain D and range R is a rule which
assigns to every x ∈ D a real number y ∈ R.
f : D → R : x → y = f (x). 0
−1 0 1
The graph of the function f is the set Fig. 2.1 A function
Γ ( f ) = {(x, y) ∈ D × R ; y = f (x)}.
In the case of D ⊂ R the graph can also be represented as a subset of the coordinate
plane. The set of the actually assumed values is called image of f or proper range:
x1 = x2 ⇒ f (x1 ) = f (x2 ).
2.1 Basic Notions 15
1.5
y
1
0.5 (x, x2 )
Γ (f )
0
D=R x
−0.5
−1 0 1
∀y ∈ B ∃x ∈ D : y = f (x).
f : D → B : y = f (x),
f −1 : B → D : x = f −1 (y),
1.5
y = x2 y = x3
1.5 1
1 0.5
0.5 0
x
−0.5
0
x −1 injective
−0.5 not injective
−1.5
−1 0 1 −1 0 1
1.5
y=x 4 y = 2x3 − x
1.5 1
1 0.5
0
0.5 x
x −0.5
0
−1 surjective
−0.5 not surjective on B = R
−1.5
−1 0 1 −1 0 1
√
x = f −1 (y) = y
x
which maps each y ∈ B to the unique x ∈ D with y = f (x) is called the inverse
function of the function f .
Once one has found the inverse function f −1 , it is usually written with variables
y = f −1 (x). This corresponds to flipping the graph of y = f (x) about the diagonal
y = x, as is shown in Fig. 2.6.
Experiment 2.7 The term inverse function is clearly illustrated by the MATLAB plot
command. The graph of the inverse function can easily be plotted by interchanging
the variables, which exactly corresponds to flipping the lists y ↔ x. For example,
2.1 Basic Notions 17
y = f (x)
x = 0:0.01:1;
y = x.ˆ2;
plot(x,y)
hold on
plot(y,x)
How the formatting, the dashed diagonal and the labelling are obtained can be learned
from the M-file mat02_1.m.
The elementary functions are the powers and roots, exponential functions and loga-
rithms, trigonometric functions and their inverse functions, as well as all functions
which are obtained by combining these. We are going to discuss the most important
basic types which have historically proven to be of importance for applications. The
trigonometric functions will be dealt with in Chap. 3.
While it was the first aim of the founder to provide dwellings for the
factory worker which should have adequate accommodation and large
gardens, it was not intended that at Bournville provision should be
made alone for the poorer working class. It might be pointed out that
one of the most prominent ideals in the scheme of the Garden Suburb
Trust, already referred to, is “that all classes may live in kindly
neighbourliness,” and the amalgamation of the factory-worker and the
brain-worker in the same district is catered for as being expressly
desirable. At Bournville there has always been a demand for houses
both on the part of the skilled artisan and others, and this demand
has been provided for from the first. Rents in the village range from
4s. 6d. a week, rates not included, to 12s. a week; and there are also
a few houses of a still larger class at higher rentals. Nor are the
houses let to Messrs. Cadbury’s own workpeople exclusively, as the
following figures will show—figures based on a private census taken
during 1901, and here quoted from a booklet issued by the Village
Trust:—
Proportion of Householders working in—
Bournville per
41·2 cent.
Villages within a mile of Bournville 18·6 „
Birmingham 40·2 „
Occupations of Householders—
Employed at indoor work in per
factories 50·7 cent.
Clerks and travellers 13·3 „
Mechanics, carpenters, bricklayers,
and various occupations not
admitting of exact classification 36·0 „
The village is four miles from Birmingham, and is easily accessible
by cycle, rail, or electric car. The last come within easy distance of the
village, workmen’s fares being 2d. return.
Under the founder’s first scheme the land was let upon leases of
999 years, subject to a ground rent varying from ½d. to 1d. per yard
(600 square yards at ½d. and 1d. = £1 5s. and £2 10s. respectively).
Arrangements were made to find capital on mortgages granted at the
rate of Three per cent. to those who paid less than half the cost of the
house and Two and a-half per cent. to those who paid more. Although
a stipulation was made that no one person should be allowed to build
more than four houses, it was found necessary to revise the
arrangement in order to prevent speculation. In 1900, therefore, the
estate was handed over to a Trust on behalf of the nation, the whole
income to be directed to solving the housing problem. The houses
now built are let to tenants at moderate weekly rentals, which include
the annual ground rent, equal to about 1d. per yard (according to its
value), and which should yield Four per cent. net. The revenue of
house and ground rents is employed, after provision has been made
for the maintenance and repair of present property, in the
development of the village itself, and in the laying out and
development of other villages elsewhere, the Trust being empowered
under the deed of foundation to acquire land in any part of Great
Britain. Subsequently to the formation of the Trust, additional land
adjacent to Bournville has been added to the founder’s gift, and
included in the village, which now extends over 458 acres. Already
upwards of 100 acres of land have been laid out for building. There
are now about 450 houses in the hands of the Trust, which number,
added to the 143 sold under the first scheme, makes a total of nearly
600. With the income of the Trust, building is being steadily
proceeded with, and there is a continual demand for houses.
PLATE VI.
SYCAMORE ROAD,
BOURNVILLE.
PLATE VII.
THE SCHOOLS,
BOURNVILLE
SEE PAGE 13.
PLATE VIII.
CARVED STONE PANELS
FOR SCHOOLS.
PLATE IX.
RUSKIN HALL,
BOURNVILLE.
PLATE X.
MEETING HOUSE,
BOURNVILLE.
PLATE XI.
MEETING HOUSE,
BOURNVILLE.
Plate xii gives the plan, with elevation, of a block of eight cottages,
the accommodation of which is the least, and the dimensions the
lowest, that should be provided for homes with one living room.
The accommodation is as follows:—
Ground Floor.
Living Room, 12 ft. 4 in. × 13 ft. Kitchen, 8 ft. × 12 ft. 6 ins. (with “Cabinet”
Bath, and boiler with patent steam exhaust). Larder under stairs.
Bedroom Floor.
First Bedroom, 9 ft. 2 ins. × 13 ft., and recess. Second Bedroom, 8 ft. 4 ins. ×
11 ft. 2 ins. Third Bedroom, 7 ft. 6 ins. × 8 ft.
Total cost, £135 per cottage.
Laying out of garden, £7 10s. extra.
Cubical contents, 64,800 ft. at 4d. per foot cube = £1,080 per
block, or £135 per cottage.
There has been considerable discussion of late with regard to the
building of cheap cottages suitable for labourers and the poorer
artisans, both in the country and elsewhere. Experiments have been
made in which the building materials employed have been other than
brick, the object being a reduction in cost. The bye-laws which do not
at present sanction the erection of cottages in some of these
materials will, it is hoped, before long be altered. Meanwhile, what is
wanted in most districts is the cheap dwelling in brick.
PLATE XIII.
BLOCK OF EIGHT COTTAGES.
SEE PAGE 16.
PLATE XV.
BLOCK OF FOUR COTTAGES.
SEE PAGE 25.
Total cost, including all extras and builder’s profit, £872 per block,
or £218 per cottage. Laying out of gardens, £10 per cottage.
Cubical contents, 48,295 ft. at 4⅓d. per foot cube, £872, or £218
per cottage.*
FRONT ELEVATION
GROUND PLAN
BEDROOM PLAN
PLATE XVI.
PAIR OF COTTAGES.
SEE PAGE 28.
PLATE XVII.
PAIR OF COTTAGES.
SEE PAGE 28.
Total cost of cottages built to this plan, including all extras, £230
per cottage.
Laying out of gardens, £10 per cottage.
Cubical contents, 22,000 ft. at 5d. per foot cube, £460, or £230 per
cottage.
PLATE XVIII.
PAIR OF COTTAGES.
FRONT ELEVATION
GROUND PLAN
BEDROOM PLAN
PLATE XVIII.
PAIR OF COTTAGES.
SEE PAGE 28.
The view here given shows three pairs of cottages built to the plan
shown on Plate xvi., and illustrates how a variety of elevation may be
gained by adding bays, dormers, etc., and by using differing materials.
D E S C R I P T I O N O F P LAT E X X .
Welcome to our website – the ideal destination for book lovers and
knowledge seekers. With a mission to inspire endlessly, we offer a
vast collection of books, ranging from classic literary works to
specialized publications, self-development books, and children's
literature. Each book is a new journey of discovery, expanding
knowledge and enriching the soul of the reade
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
textbookfull.com