Complete Elements of Matrix Modeling and Computing With MATLAB 1st Edition Robert E. White PDF For All Chapters
Complete Elements of Matrix Modeling and Computing With MATLAB 1st Edition Robert E. White PDF For All Chapters
com
https://ebookname.com/product/elements-of-matrix-modeling-
and-computing-with-matlab-1st-edition-robert-e-white/
OR CLICK BUTTON
DOWNLOAD EBOOK
https://ebookname.com/product/computational-mathematics-models-
methods-and-analysis-with-matlab-and-mpi-1st-edition-robert-e-
white/
https://ebookname.com/product/signals-and-systems-with-matlab-
computing-and-simulink-modeling-4th-edition-steven-t-karris/
https://ebookname.com/product/micromechatronics-modeling-
analysis-and-design-with-matlab-2nd-ed-edition-giurgiutiu/
https://ebookname.com/product/helixtrolysis-cyberology-the-
joycean-machine-1st-edition-louis-armand/
Negotiating the Therapeutic Alliance 1st Edition Jeremy
D. Safran
https://ebookname.com/product/negotiating-the-therapeutic-
alliance-1st-edition-jeremy-d-safran/
https://ebookname.com/product/exhaustivity-contrastivity-and-the-
semantics-of-mandarin-cleft-related-structures-1st-edition-ying-
liu/
https://ebookname.com/product/advanced-nutrition-and-human-
metabolism-5th-edition-sareen-s-gropper/
https://ebookname.com/product/young-europe-what-young-people-
between-the-ages-of-9-and-14-think-about-participation-politics-
and-europe-1st-edition-hannah-olle/
https://ebookname.com/product/final-cut-pro-workflows-the-
independent-studio-handbook-1st-edition-jason-osder/
Food That Really Schmecks Life Writing Edna Staebler
https://ebookname.com/product/food-that-really-schmecks-life-
writing-edna-staebler/
Elements of Matrix
Modeling and
Computing
with
MATLAB ®
Elements of Matrix
Modeling and
Computing
with
MATLAB ®
Robert E. White
This book contains information obtained from authentic and highly regarded sources. Reasonable efforts
have been made to publish reliable data and information, but the author and publisher cannot assume
responsibility for the validity of all materials or the consequences of their use. The authors and publishers
have attempted to trace the copyright holders of all material reproduced in this publication and apologize to
copyright holders if permission to publish in this form has not been obtained. If any copyright material has
not been acknowledged please write and let us know so we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced, transmit-
ted, or utilized in any form by any electronic, mechanical, or other means, now known or hereafter invented,
including photocopying, microfilming, and recording, or in any information storage or retrieval system,
without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.copyright.
com (http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC), 222 Rosewood
Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that provides licenses and
registration for a variety of users. For organizations that have been granted a photocopy license by the CCC,
a separate system of payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are used
only for identification and explanation without intent to infringe.
Visit the Taylor & Francis Web site at
http://www.taylorandfrancis.com
List of Tables xi
Preface xiii
Introduction xv
2 Vectors in Space 47
2.1 Vectors and Dot Product . . . . . . . . . . . . . . . . . . . . . . 47
2.2 Cross and Box Products . . . . . . . . . . . . . . . . . . . . . . . 56
2.3 Lines and Curves in R3 . . . . . . . . . . . . . . . . . . . . . . . 67
2.4 Planes in R3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
2.5 Extensions to Rq . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
3 Ax = d: Unique Solution 95
3.1 Matrix Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
3.2 Matrix Products . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
3.3 Special Cases of Ax = d . . . . . . . . . . . . . . . . . . . . . . . 117
3.4 Row Operations and Gauss Elimination . . . . . . . . . . . . . . 127
3.5 Inverse Matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
3.6 OX Factorization . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
3.7 Determinants and Cramer’s Rule . . . . . . . . . . . . . . . . . . 159
v
vi CONTENTS
Bibliography 397
Index 399
List of Figures
vii
viii LIST OF FIGURES
xi
Preface
xiii
xiv PREFACE
developed from very simple models to more complex models. The reader can
locate sections pretaining to a particular application by using the index.
°R
MATLAB is used to do some of the more complicated computations. Al-
though the primary focus is to develop by-hand calculation skills, most sec-
tions at the end have some MATLAB calculations. The MATLAB m-files used
in the text are listed in the index and are included in the book’s Web site:
http://www4.ncsu.edu/~white. The approach to using computing tools in-
cludes: first, learn the math and by-hand calculations; second, use a computing
tool to confirm the by-hand calculations; third, use the computing tool to do
more complicated calculations and applications.
I hope this book will precipitate discussions concerning the core mathemat-
ical course work that scientists and engineers are required to study. Discrete
models and computing have become more common, and this has increased the
need for additional study of matrix computation, and numerical and linear al-
gebra. The precise topics, skills, theory and appropriate times to teach these
are certainly open for discussion. The matrix algebra topics in this book are
a small subset of most upper level linear algebra courses, which should be en-
hanced and taken by a number of students. This book attempts to make a
bridge from two- and three-variable problems to more realistic problems with
more variables, but it emphasizes skills more than theory.
I thank my colleagues who have contributed to many discussions about the
content of this text. And, many thanks go to my personal friends and Liz White
who have listened to me emote during the last year.
Bob White
MATLAB is a registered trademark of The MathWorks, Inc. For product
information, please contact:
One can view an p×q matrix as a table of objects with p rows and q columns.
The objects are usually real or complex numbers, but they could be characters
or records of information. A simple example is data for the last 12 months of
car sales where there are p = 12 rows and q = 2 columns. The first column
will have the month’s number and the second column will have the number of
cars sold in the corresponding month. By examining the data one would like
to make a prediction about futures sales. This is where the modeling enters.
If the graph of the sales versus months "looks" like a straight line, then the
data may be modeled by a linear function of time | = pw b + f. The slope p b
and intercept f must be chosen so that the computed sales are "close" to the
car sales data. This is done by appropriate manipulations of the two column
vectors and computing a solution of the resulting system of algebraic equations.
Once p b and f have been found, the predicted sales for w larger than 12 can
easily be calculated by evaluating the linear function. The modeling process is
complicated by incorrect sales data, changing prices and other models such as
a parabolic function of time.
This text examines a variety of applications, which have matrix models and
often have algebraic systems that must be solved either by-hand calculations
or using a computing tool. Applications to projectiles, circuits, mixing tanks,
trusses, heat conduction, motion of a mass, curve fitting and image enhancement
will be initially modeled in very simple ways and then revisited so as to make
the model more accurate. This is typical of the modeling process where there is
an application, a model, mathematical method, computations and assessment
of the results. Then this cycle is repeated so as to enhance the application’s
model.
The first two chapters deal with problems in two- and three-dimensional
space where the matrices have no more than three rows or columns. Here
geometric insight can be used to understand the models. In Section 2.5 the
extension to higher dimensions is indicated for vectors and matrices, solution
to larger algebraic systems, more complicated curve fitting, time dependent
problems with systems of dierential equations and image modeling. Chapters
three, four and five have the basic matrix methods that are required to solve
systems in higher dimensions. Chapters six and seven contain time dependent
models and introduce linear systems of dierential equations. The last two
xv
xvi INTRODUCTION
1 }3 }2 }
Chapter 1
1
2 CHAPTER 1. VECTORS IN THE PLANE
The set of rational numbers has a countable but infinite number of elements.
Also, the addition and product of two rational numbers are rational numbers.
{ = ±(={1 · · · {g · · · )10h
±({1 @10 + · · · + {g @10g + · · · )10h =
Real numbers contain the rational numbers, but not all real numbers are
rational. For example, consider { = 31@2 where 3 is a prime number. If { were
a rational number, then { = p@q giving 3 = p2 @q2 and 31 q2 = p2 . The
left side has an odd number of prime factors 3, and the right side has an even
number of prime factors 3. This contradicts the unique factorization property
and, hence, { cannot be a rational number.
This is a floating point number with base equal to 10 where {1 is not equal
to zero, {l are integers between 0 and 9, the exponent h is an integer between
given integers i and j and g is a positive integer called the precision of the
floating point system
Associated with each real number, {, and its floating point approximate
number, i o({), is the floating point error, i o({) {. This error decreases as
the precision, g, increases. Each computer calculation has some floating point
1.1. FLOATING POINT AND COMPLEX NUMBERS 3
error. Moreover, as additional floating point calculations are done, there may
be an accumulation of these floating point errors.
Example 1.1.1. Let { = 1=5378 and i o({) = 0=154 101 where g = 3. The
floating point error is
i o({) { = =0022=
The error will accumulate with any further operations containing i o({), for
example, i o({)2 = =237 101 and
If } 5 C> then real(}) = d is called the real part and imag(}) = e is called the
imaginary part. Two complex numbers are equal if and only if both the real
and imaginary parts are equal.
Complex numbers are often represented graphically where the real part is
located on the horizontal axis and the imaginary part is located on the vertical
axis. The complex number is then either viewed as point in the complex plane
or as directed line segment or arrow from the origin 0 + l0 to the complex point
} = d + el= Figure 1.1.1 illustrates this for } = 2 + 1l= One point of confusion
concerns the starting point of the arrow, which is the origin, but often the arrow
is moved parallel to its original position, as is illustrated by the top dashed arrow
in Figure 1.1.1 that is parallel to }= This is done to depict addition of two arrows
by placing the beginning of the arrow at the end of another arrow.
There are several important operations: conjugate of }, modulus of }, addi-
tion, subtraction, product and division. The conjugate of } is another complex
number formed by the reflection of } about the real axis. The modulus of } is
the length of the vector. The algebraic operations are defined in the natural
way by using l2 = 1=
Definitions 1.1.1. Attributes of complex numbers } = d + el and z = f + gl=
} + 7 + 2l = 3 + 4l=
Addition of complex numbers is defined by adding the real and imaginary parts,
} + 7 + 2l = ({ + l|) + 7 + 2l = ({ + 7) + (| + 2)l and, then, the above two real
1.1. FLOATING POINT AND COMPLEX NUMBERS 5
equations are
{ + 7 = 3 and
| + 2 = 4=
Thus, the complex solution of } + 7 + 2l = 3 + l4 is } = 4 + 2l=
Definitions 1.1.3. Operations with complex numbers } = d+el and z = f+gl=
}±z (d ± f) + (e ± g)l
}z (d + el)(f + gl) = (df eg) + (ef + dg)l and
d + el d + el f gl df + eg ef dg
}@z = = 2 + 2 l=
f + gl f + gl f gl f + g2 f + g2
The following eleven algebraic properties are easily established. The first five
properties are the additive group properties, the next five are the multiplicative
group properties for non-zero complex numbers and the last is the distributive
property.
Theorem 1.1.1 (Field Properties of Complex Numbers) Let }> z> z1 > z2 5 C=
Then
} + z 5 C, } + 0 = }, } + (}) = 0, } + z = z + },
} + (z1 + z2 ) = (} + z1 ) + z2 ;
}z 5 C, }1 = }, }(} 1 ) = 1 when } 6= 0, }z = z},
}(z1 z2 ) = (}z1 )z2 ;
}(z1 + z2 ) = }z1 + }z2 =
Example 1.1.6. This example illustrates how one can solve complex algebraic
equations where we must solve for both the real and imaginary parts of an
unknown complex number } = { + l|= Here we use the above field properties
1 + 4l
= 2+l
}+1
1 + 4l = (2 + l)(} + 1)
= (2 + l)} + (2 + l)=
1 + 4l 5 + 20l
=
(1@5 + (7@5)l) + 1 6 + 7l
5 + 20l 6 7l
=
6 + 7l 6 7l
170 + 85l
=
85
= 2 + l=
The complex numbers are extensions of the real numbers, and many of
the algebraic and absolute value properties extend to the complex numbers.
For example, if d and f are real numbers and |d| is the absolute value, then
|df| = |d| |f|. When d and f are extended to complex numbers and absolute
value is extended to the modulus of a complex number, then one can show
1.1. FLOATING POINT AND COMPLEX NUMBERS 7
|}z| = |}| |z| = In order to see this is true, let } = d + el and z = f + gl and
show |}z|2 = |}|2 |z|2 :
¯ ¯
If } = z, then we could write this as ¯} 2 ¯ = |}|
2
1.1.6 Exercises
1. Consider the real number 1@3.
(a). Find i o(1@3) with base 10 and three digits, g = 3.
2 2
(b). Compute (i o(1@3)) and the error = (i o(1@3)) (1@3)2 =
1@2
2. Let s be prime number. Show s is not a rational number.
3. Let } = 7 l and z = 2 + 3l=
(a). Graph }, z, z and 2z + 3}=
(b). Compute }z and }@z=
4. Let } = 2 + 3l and z = 2 3l=
(a). Graph }, z, z and 2z + 3}=
(b). Compute }z and }@z=
5. For } = cos(2@3) + l sin(2@3) compute 1 + } + } 2 =
6. For } = cos(2@4) + l sin(2@4) compute 1 + } + } 2 + } 3 =
7. Use the properties in Theorem 1.1.1 to solve for the complex number
} = { + |l
} + 7 + 2l = 3 + 4l=
8. Prove the remaining properties in Theorem 1.1.1.
9. Find the real and imaginary parts of
2+l
=
3 7l
10. Find the real and imaginary parts of
2+l
+ 5=
10 + 7l
11. Solve for the complex number } = { + |l
2 + 4l
= 7 + l=
}+l
12. Solve for the complex number } = { + |l
2 + 4l
= 1 + l=
2} l
13. Prove the remaining properties in Theorem 1.1.2.
14. Does the quadratic formula hold when the coe!cients d> e and f are
complex numbers? You may wish to consider some special cases such as d =
1> e = 0 and f = l= One needs to be able to define the square root of a complex
number!
15. Use MATLAB and enter } = 2 + 3l and z = 2 3l=
(a). Graph } and z=
(b). Compute }z and }@z=
(c). Use quiver() to plot the vectors for } and z=
16. Consider the calculations
s in Figure 1.1.2 given by the MATLAB code
complex_prod.m. Let } = 3 + 1l=
(a). Use by-hand calculations to find } 2 , } 3 , their norms and arguments.
(b). Modify complex_prod.m to confirm your by-hand calculations.
10 CHAPTER 1. VECTORS IN THE PLANE
i (}) = d} + e=
The constants d and e are given complex numbers, and } = {+l| is any complex
number where { and | are real numbers. This means
the functions x({> |) and y({> |) have real values. The solution of i (}) = g =
g1 + lg2 is equivalent to equating the real and imaginary parts and reduces to
solving two real equations for { and |
{ | + 3 = 0 and { + | 1 = 7=
The solution of i (}) = 2 + l can be found by equating the real and imaginary
parts x({> |) = 2, y({> |) = 1 and solving for { and |= Or, one can use the
properties in Theorem 1.1.1 to find } as follows:
}+l
= 2+l
}l
}+l = (2 + l)(} l) = (2 + l)} (2 + l)l
}(1 (2 + l)) = (2 + l)l l
}(1 l) = 3l + 1
1 + 3l 1 l 2 + 4l
} = = = 1 + 2l=
1+l 1l 2
A quadratic complex function is
i (}) = d} 2 + e} + f=
The coe!cients d, e and f are given complex numbers, and one would like to
solve i (}) = d} 2 + e} + f = 0= In order to generalize the quadratic formula
to the complex case, the square root of a complex number needs to be defined.
Since the square root function is the inverse of the square function, let us focus
on
i (}) = } 2 =
12 CHAPTER 1. VECTORS IN THE PLANE
In Section 1.1 we used the polar representation } = |}| (cos() + l sin()) and
Theorem 1.1.2 with } = z to obtain
Thus, one must choose x and y such that { = x2 y 2 and | = 2xy= The use of
trigonometric identities confirms
1@2 1@2
x = |}| cos(@2) and y = |}| sin(@2)
are the solution of these two equations and, hence, equation (1.2.2) is correct
for the definition of square root of }=
s
Example 1.2.3. Use equation (1.2.2) to find the square root of } = 1 + 3l.
The modulus of } is 2 and the angle is @3 or 60 degrees so that } = 2(cos(@3)+
l sin(@3))= Equation (1.2.2) gives
s 1@2
(1 + 3l) = 21@2 (cos(@6) + l sin(@6))
s s
= 2( 3@2 + l(1@2))
p s
= 3@2 + l(1@ 2)=
s
Note the solution of the complex quadratic equation } 2 = 1 + 3l is plus or
minus this square root.
Figure 1.2.1 illustrates the linear, square and square root functions of a
complex variable; this was generated by the MATLAB code complex_power.m.
The complex numbers in the domain correspond to the domain semicircle with
radius equal to three. The linear function has values in the complex plane where
the domain semicircle has been doubled in radius and shifted to the right. The
square function has values in the complex plane where the domain semicircle
has been expanded in both the radial direction and in an angular direction to
one complete circle of radius equal to three squared. The smallest curve depicts
the square root function. Now, the domain semicircle has been contracted in
both the radial and angular components to one half of the semicircle and with
radius equal to the square root of three.
1.2. COMPLEX VALUED FUNCTIONS 13
semester of calculus,
{2 {3 {4 {5
h{ = 1+{+ + + + + ··· >
2! 3! 4! 5!
{2 {4
cos({) = 1 + · · · and
2! 4!
{3 {5
sin({) = { + ··· =
3! 5!
Extend the exponential function’s domain from the real numbers to the complex
numbers by replacing { by l|
(l|)2 (l|)3 (l|)4 (l|)5
hl| = 1 + l| + + + + + ···
2! 3! 4! 5!
|2 |3 |4 |5
= 1 + l| l + + l + ···
2! 3! 4! 5!
|2 |4 |3 |5
= [1 + · · · ] + l[| + ···]
2! 4! 3! 5!
= cos(|) + l sin(|)=
} = { + l|
= |}| (cos() + l sin ())
= |}| hl =
Second, this leads to the representation of the cosine and sine functions
Third, consider the product of hl{ and hl| and use the trigonometric identities
for the sum of angles
Fourth, another identity, which follows from the repeated use of hl{ hl| = hl({+|)
with | = {> is de Moivre’s formula
This leads to all q solutions } = (hl2@q )m = hl(2@q)m where m = 0> 1> ===> q 1=
They are solutions because
}q = (hl(2@q)m )q
= hl(2)m
= cos((2)m) + l sin((2)m)
= 1=
Definition 1.2.2. Let m = 0> 1> · · · > q1= Then hl(2@q)m are solutions to } q = 1
and are called the qwk roots of unity.
Example 1.2.5. Let q = 3 and find the three solutions of } 3 = 1= Here
2@q = 2@3 or 120 degrees. The three solutions are
hl(2@3)0 = cos((2@3)0) + l sin((2@3)0) = 1
hl(2@3)1 = cos((2@3)1) + l sin((2@3)1) = 0=5000 + 0=8660l
hl(2@3)2 = cos((2@3)2) + l sin((2@3)2) = 0=5000 0=8660l=
Another special case with q = 12 is illustrated in Figure 1.2.2 where the 12
vectors with radius equal to one are the solutions. The other 12 vectors have
increasing moduli and have staggered angles
(1=05h(l2(1+=05)@q) )m =
The figure was generated by the MATLAB code euler_form.m.
The solution of the more general problem } q = g is straightforward once
the complex number g is written in polar form
g = |g| hl = |g| hl hl2
1@q
} = g1@q = |g| (hl hl2 )1@q
= |g|1@q hl@q hl2@q = (1.2.10)
l2@q m q
Define z h and use the fact (z ) = 1 to obtain all q solutions } =
g1@q zm where m = 0> 1> ===> q1= A summary of the above is given in the following
theorem.
Theorem 1.2.1 (Solutions of } q = g)= Let z hl2@q and g = |g| hl = Then
z satisfies z = z1 > zq = 1 and 1+z+· · ·+zq1 = 0= Moreover, the solutions
1@q l@q m
of } q = g are given by } = |g| h z = g1@q zm where m = 0> 1> · · · > q 1=
Example 1.2.6. Let q = 3 and find the three solutions of } 3 = 1 + l= Here
2@q = 2@3 and g = 1 + l = 21@2 hl@4 = The three solutions are
21@6 hl@12 hl(2@3)0 = 21@6 (cos(@12 + (2@3)0) + l sin(@12 + (2@3)0))
= 1=0842 + 0=2905l
21@6 hl@12 hl(2@3)1 = 21@6 (cos(@12 + (2@3)1) + l sin(@12 + (2@3)1))
= 0=7937 + 0=7937l
21@6 hl@12 hl(2@3)2 = 21@6 (cos(@12 + (2@3)2) + l sin(@12 + (2@3)2))
= 0=2905 1=0842l=
1.2. COMPLEX VALUED FUNCTIONS 17
square and square root functions. An array operation on a row vector produces
another row vector whose components are given by the operation. For example,
if d = [2 3]> then d=ˆ2 = [22 32 ]=
The MATLAB code euler_form.m, see Figure 1.2.2, uses a for-loop to gen-
erate a sequence of vectors associated with the complex exponential function
raised to a variety of powers. The vectors are created by the MATLAB com-
mand quiver(x,y,dx,dy,1), which graphs a vector beginning at the point ({> |)
and ending at ({ + g{> | + g|).
1.2.6 Exercises
1. Let i (}) = (2 + l)} + 7 + l where } = { + l|
(a). Find x({> |) and y({> |) such that i (}) = x + ly=
(b). Find the solution of i (}) = 10 2l=
2. Let i (}) = (3 + l)} + 1 + 2l where } = { + l|
(a). Find x({> |) and y({> |) such that i (}) = x + ly=
(b). Find the solution of i (}) = 1 2l=
3. Let i (}) = (} + 2l)@(} l) where } = { + l|
(a). Find x({> |) and y({> |) such that i (}) = x + ly=
(b). Find the solution of i (}) = 1 2l=
4. Let i (}) = (2} + l)@(} + l) where } = { + l|
(a). Find x({> |) and y({> |) such that i (}) = x + ly=
(b). Find the solution of i (}) = 1 + 2l=
5. Use equation (1.2.2) to find the square root of 1 + l= s
6. Use equation (1.2.2) to find the square root of 1 + 3l=
7. Use the complex quadratic formula to solve d} 2 + e} + f = 0=
(a). Solve the quadratic equation } 2 + 2} l = 0=
(b). Check your answer.
8. Consider the completing square method of developing the quadratic for-
mula.
(a). Prove: if } and z are complex numbers and }z = 0> then either
} = 0 or z = 0=
(b). Use this fact to derive the quadratic formula for complex numbers.
9. Use de Moivre’s formula in equation (1.2.9) to compute (1 + l)10 =
10. Use de Moivre’s formula in equation (1.2.9) to compute (3 + 2l)7 =
11. Use Theorem 1.2.1 to find all three solutions to } 3 = 3 + 4l=
s
12. Use Theorem 1.2.1 to find all four solutions to } 4 = 1 + 3l=
13. Let z hl2@q = Show z satisfies z = z1 > zq = 1 and 1 + z + · · · +
q1
z = 0=
14. Use MATLAB to verify the third equation in Example 1.2.6.
15. Use MATLAB to experiment with dierent size semicircles in the code
complex_power.m.
16. Use MATLAB to experiment with powers q = 3> 6 and 12 in the code
euler_form.m.
1.3. VECTORS IN R2 19
1.3 Vectors in R2
Complex numbers are represented either by ordered pairs of real numbers or as
directed line segments. The first number in the pair is the real part and the
second number is the imaginary part. The directed line segment or arrow starts
at the origin and ends at the complex point. There are a number of important
similar examples such as displacement in the plane of a mass, force acting on
a mass as well as the velocity and acceleration of a mass. As in the complex
numbers, these quantities can be added, subtracted and scaled, and have length
and angles.
In order to move an object, a force must be exerted on the mass, and this
will also have horizontal and vertical components. By adjusting the size of the
components, the direction and magnitude of the force can be determined to
obtain the desired motion. Let i1 be the force in the horizontal direction, i2
be the vertical force and denote the force vector as
$
i = [i1 i2 ] = force vector.
The resulting motion will have speeds in the horizontal and vertical directions.
The velocity vector lists these speeds so that the velocity vector also has a
direction and magnitude. The speed in the horizontal direction is y1 and when
this is coupled with the vertical speed y2 we get a velocity vector
$
y = [y1 y2 ] = velocity vector.
Newton’s law of motion can be applied in both the horizontal and vertical
directions
gy1 gy2
i1 = p and i2 = p =
gw gw
1.3. VECTORS IN R2 21
Figure 1.3.2: f2 = d2 + e2
Another result pertains to the case when the angle is not 90 degrees as is
illustrated in the triangle in Figure 1.3.3 where the triangle with sides d> e and
f is a union of two right triangles. The left right triangle has base equal to
e cos() and height equal to e sin()= The other right triangle has the same
height, but the base is d e cos() (assuming d A e)= Therefore, applying the
Pythagorean theorem to the right triangle on the right side gives
f2 = (e sin())2 + (d e cos())2
= e2 sin2 () + d2 2de cos() + e2 cos2 ()
= e2 (sin2 () + cos2 ()) + d2 2de cos()
= e2 + d2 2de cos()=
f2 = e2 + d2 2de cos()=
Y our fearfulness, (says the Prefect) surprises me. Tell me, I pray,
what idea hast thou of what is called grandeur, dignities, and
high rank in a state?
I am in this world (answered I) like a traveller, who goes on his
way curiously observing the objects, but desiring none, because he is
but a passenger. Moreover, if things are estimated according to the
happiness they procure, I do not think that the highest places should
be much valued; for, I see, they make no man happy, and are a
misfortune to many.
What of riches? added the Prefect.
Pleasure (said I) is like a very rare commodity, which, however,
every one would fain purchase. Among those that succeed, the rich
buy it very dear, it comes cheap to the rest: One may as well be
among the last as the first. Of the few pleasures that exist, the lower
class enjoy as large a share as the highest.
What of wit, genius, talents? says the Prefect.
One half of the world, replied I, study to amuse the other. The first
class is formed of men of talents; whose brains are wound up by
nature higher than ordinary. They are incessantly striving to please:
If they fail, they waste away with grief; if they succeed, it is never
fully, and a single censure creates them more pain than all the
encomiums together give them pleasure. It is, therefore, better to be
of the second class, I mean among those who are amused by the
others.
As far as I see, said the Prefect, the aspect of the great and their
pomp, of the scholar and his extensive genius, of the rich and his vast
possessions, makes little or no impression on thy mind.
I confess, replied I, that no man was ever less dazzled with all this
than myself. Wrapt in a certain coolness of sense, I am guarded
against all strong impressions. I behold with the same eye the
ignorant who know nothing, and the learned who know all, except
truth; the protector who plans, though he knows his weakness, and
the protected who cringes, though he perceives his superiority; the
peasant that is disgusted with the simplicity of his diet, and the rich
sensual, who with thirty niceties, can hardly make a dinner; the
duchess, loaded with diamonds, and the shepherdess decked with
flowers; vanity, which dwells in the cottage as well as in the palace,
and upholds the low as well as the high; care, which sits on the
throne by the king, or follows the philosopher in his retirement. All
the parts on the stage of this world, seem to me one no better than
another: but I do not desire to act any. I would observe all and be
taken up with nothing. Hence it is, that I dreaded the neighbourhood
of these restless flies....
And hence it is precisely, interrupted the Prefect, that thou hadst
nothing to fear from them. Thou admirest nothing; it is sufficient:
The flies can take no hold of thee. The first impression they must
make, is the impression of surprise and admiration; if they make not
that, they miss their aim. But the moment admiration is admitted, a
crowd of passions quickly follow. For, in the object of wonder, great
hurt or great good is expected. Hence Love or Aversion, and all their
attendants; restless Desire which never sleeps; Joy, which embraces
and devours its objects; Melancholy, which, at a distance, and with
weeping eyes, contemplates and calls for what it dreads: Confidence,
which walks with head erect, and often meets a fall; Despair, which is
preceded by fear and followed by madness, and a thousand others. If
thou wilt rest secure from their attacks, cherish thy coolness of sense,
and never lose sight of the grand principle,
Nil Admirari.
CHAP. X.
The Fantastical Tree.
I have one thing more (said the Prefect) to show thee; prepare thy
eyes and thy ears; and be frightened at nothing.
The rivulet, by the side of which we walked to the Fantastical Tree,
receives several streams as it flows along; and, as if it left with regret
so beautiful a residence, after forming a thousand serpentine
windings in the meadow, it glides gently towards its mouth. In that
place, a hole, formed by an opening of the earth, receives and
transmits it through subterraneous channels.
We came to the place where it was broadest. The bottom was of
smooth gravel, and the water not above an inch deep. The Prefect
went in and I followed him.
I had gone but a few paces, when the bottom gave way: I sunk, but
it was only to my waste; and I remained in that posture, without
being able to get to one side or the other. Fear nothing, says the
Prefect, calmly enjoy the last spectacle I have reserved for thee.
I then gave myself up to the efforts of the waters, which carried me
away, and I soon entered into the subterraneous cavities, where they
were lost. At a little distance, the rivulet flowed into another, and
soon after, both ran into a river. I was carried from stream to stream;
I crossed gulphs, lakes, and seas.
As long as a faint light permitted, I contemplated the internal
frame of the earth. It is a labyrinth of immense caverns, deep grottos,
irregular crevices, which have a communication with one another.
The waters that flow in these subterranean places, spread themselves
sometimes into vast basons, and seem to stagnate; sometimes they
run with a rapid stream through narrow straits; and dash against the
rocks with such impetuosity, as to produce the phosporus and flashes
of lightening; very often they fall from the top of the vaults with a
dreadful noise. The dazzled eye sees, as it imagines, the foundations
of the earth shake; one would think, that the whole was turned
upside down, and falling into chaos.
When the glimmering light, which I had enjoyed some time, came
to fail, I found myself buried in profound darkness, which increased
the horror, I had conceived at what I had seen. A hideous noise,
mixed with the murmuring of the streams, with the whistling of the
gulfs, with the roaring of the torrents, threw me into great
perturbation of mind; and my troubled fancy formed to itself a
thousand frightful images.
I went on a good while in this darkness; and I know not how far I
had gone when a faint light struck my eyes. It was not like that which
precedes sun-rising, or follows sun-set; but that melancholy light,
which a town on fire spreads at a distance in the shade of the night. I
was some time before I saw whence it came: At last, I found myself
close to the most terrible of all the sights.
A vast opening exposed to my eyes in an immense cavern, an abyss
of fire. The devouring flame rapidly consumed the combustible
matter with which the arched roofs of the abyss were impregnated. A
thick smoke mixed with fiery sparks, diffused itself to a great
distance. From time to time, the calcined stones fell down by pieces,
and the liquified metals formed flaming streams. Sometimes whole
rocks, rent from the tops of the vaults, gave passage to water, which
poured down in boiling streams. The moment the water touched the
calcined matters and melted minerals, it caused most shocking
detonations: The concavities of the globe resounded, their
foundations were shaken: And I conceived that such was the cause of
those terrible earth-quakes, that have destroyed so many countries,
and swallowed up so many cities.
I was soon in darkness again; for I still went on. Every moment I
should have been destroyed, if the Prefect of Giphantia had not
watched over me. I saw him no more: But his promise was with me:
And the dangers, I had escaped, heartened me against those I had
still to undergo. By degrees I took courage, and became so easy as to
make some reflections.
Alas! said I, through a frightful desart I came into the most
beautiful mansions in the world, and I am now going thence through
gulfs, abysses, and vulcanos. Good and evil closely follow one
another. It is thus, the light of the day and darkness of the night, the
frosts of the winter and the flowers of the spring, the gentle zephyrs
and the raging storms, succeed one another. However, by this
strange concatenation, is formed the enchanting prospect of nature.
Let us not doubt it: The natural world, notwithstanding its disorders,
is the master-piece of infinite wisdom; the moral world, in spite of its
stains, is worthy the admiration of the philosopher: And Babylon,
with all its faults, is the chief city of the world.
At last, after many days of subterraneous navigation, I once more
saw the light; I came out of these terrible vaults, and the last current
landed me upon a maritime coast. The serenity of the air was not
ruffled with the wind; the calm sea shone with the rays of the rising-
sun; and, like a tender wife who stretches out her arms, and sweetly
smiles on a beloved husband, the earth seemed to resume new life at
the return of that glorious orb, from whence springs all its fertility.
By degrees, my troubled senses were calmed: I looked round me, and
found myself in my own country, six hundred furlongs north-west
from Babylon, to which city I address and dedicate this narrative of
my hazardous travels.
FINIS.
1. The Jansenists (so called from Jansenius bishop of Ypres) explained the
Doctrine of Grace after the Calvinistical or rather Methodistical manner, whilst the
Molinists (so named from Molina a Spanish Jesuit) explained it after the Arminian
or rather Semi-pelagian way. The Gallican clergy were divided between these two
Opinions.
The reader may remember, there are three opinions concerning Grace. Says
the Calvinist and Methodist, Grace does ALL. Says the Arminian and Semi-pelagian,
Grace does HALF. Says the Pelagian, Grace does NOTHING.
2. The city of Ombi stood on the eastern side of the Nile, and Tentyra or
Tentyris on the western; both in Thebais part of Upper Egypt. The Tentyrites were
professed enemies of the Crocodiles, whilst the rest of the Egyptians held them in
great veneration, especially the Ombites, who for their sake waged war with the
Tentyrites.
3. Our author in this and the following chapter gives a very lively summary of
the four great monarchies of the world.
I. The Assyrian or Babylonian founded by Nimrod (or Belus I.) soon after the
dispersion at Babel, and which ended with the taking of Babylon (A. C. 538) by
Cyrus who founded II. The Persian empire which ended with the defeat of Darius
Codomannus (A. C. 334) by Alexander the Great who founded III. The Grecian or
Macedonian empire which in about five years was divided among his successors,
and at length (after the battle of Actium and death of Cleopatra) became subject to
IV. The Roman empire under Augustus Cæsar, of which there are still some
remains.
5. After the death of Sardanapalus (who is said to burn himself, his wives and
concubines, his eunuchs and riches, in one of the courts of his palace) the empire
was divided into the Median over which Arbaces reigned at Nineveh, and the
Assyrian over which Belesis reigned at Babylon. These were united under Cyrus
about 210 years after. Belesis (the Baladan of Scripture) is called also Nabonassar.
From the first year of his reign begins the famous Astronomical Æra of
Nabonassar, containing 908 years from February 26 before Christ 747, to the 23d
year of Antoninus Pius in the year of our Lord 161.
7. By a solemn treaty Ptolemy had Egypt, &c. Cassander had Macedonia and
Greece. Lysimachus had Thrace, Bithynia, &c. Seleucus had Syria, &c. Of these, the
kingdom of Egypt (under 14 monarchs including Cleopatra) and of Syria (under 27
kings) subsisted till subdued by the Romans. The rest soon fell to pieces.
8. His Library is said to consist of above 200,000 volumes. Among the rest
was the Septuagint or Greek translation of the Old Testament A. C. 267. done by
Ptolemy’s order. This library was at last destroyed by fire.
9. This man who from a huntsman raised himself to the throne of Lusitania
(now Portugal) defeated the Romans in several battles; so that Cepion the consul
was forced at last to have him murdered by treachery. He was (says Livy) much
lamented and honorably buried.
10. Rome was taken by Alaric king of the Goths in 410. By Genseric the Vandal
in 455. By Odoacer king of the Heruli in 465, and by Totila the Goth in 546, by
whom it was miserably plundered.
11. Attila king of the Huns, (called the scourge of God) after his other
devastations entered Gaul with 500,000 Men and was defeated in the plains of
Chalons in 451, with the loss of 200,000 Huns. After which he wasted Italy and
destroyed Aquileia and other places. Then returning home, he died on his wedding
night. The Huns were the most terrible of all the northern swarms. By the very
terror of their countenances they are said to over-run the Scythians, Alans and
Goths. They were so ignorant as not to know letters.
12. Mahomet was born at Mecca in Arabia, May 5, 570. He is thought by some
to be persuaded that he was really inspired to propagate the belief of one God, and
to overthrow the idolatrous religion of his country. If he retained some absurd
notions, it was (say they) to induce his countrymen to embrace his religion. The
Mahometan æra begins July 16, 622, when he fled from Mecca to Medina. He died
Jan. 17, 631, after having reduced Arabia to his obedience. His religion has since
spread itself over Asia, Africa, and great part of Europe.
13. Soliman, father of the Othman race, came out of Scythia with 50,000 men
in the year 1214, and pushed his conquests to the Euphrates. In attempting to pass
that river he was drowned in 1219. Othman his grandson was declared sultan in
1300. Mahomet II. the seventh emperor of the Turks, put an end to the Eastern
empire by taking Constantinople in 1453. The Turks embraced the religion of
Mahomet.
14. Camayeu, is a stone, whereon are found various figures formed by nature.
It is the name the orientals give the onyx, on which and on agate, these natural
figures are often found. When the figures are perfected by art, it is still called a
camayeu, as is also a painting in one colour, representing basso relievos.