E The EXPONENTIAL - The Magic Number Of: Growth
E The EXPONENTIAL - The Magic Number Of: Growth
E The EXPONENTIAL - The Magic Number Of: Growth
GROWTH
Keith Tognetti
School of Mathematics and Applied Statistics
University of Wollongong NSW 2522 Australia
12 February 1998
Introduction
This Module is written as a self contained introduction to e, bringing together the main
theorems and important properties of this fundamental constant of natural growth processes.
Not only is an axiomatic treatment given but this is complemented by some theorems that
have been selected for their unexpected beauty.
It assumes only an elementary understanding of integration. Some anecdotes and history
are also included. It will be seen that Eulers gamma constant also comes into the main
theorem and this highlights the intimate interconnection between e, the area under the curve
1
x
, and the the truncated Harmonic series.
The number e is the base of Natural logarithms but it is not the base of Naperian logarithms
as will be seen by glancing at the short historical note on Naperian logarithms in Appendix
A.
1
e owes its importance to mathematics because we come across it in all natural growth
processes.
e like is transcendental and thus it cannot be expressed as the root of a nite polynomial
equation with rational coefcients (go to Appendix Cfor some more about irrationals). If you
read nothing else at least to help you get a feel for arithmetic look at Appendix D, the article
by the Noble Prize winning physicist and sometimes lovable character Richard Feynman.
For a very general, exhaustive and more elementary introduction to e see Maor (Refer-
ences).
Let us begin our description of e by writing down the rst digits
e = 2.7182818284590452353602874...
We immediately note that 1828 is repeated at the beginning. This is rather exceptional as
most of the important constants of mathematics that are irrational have no obvious pattern at
their beginning. However it is almost certain that if we searched through the digits of most of
these numbers we would come across any pattern we liked to dene, as long as the pattern was
not periodic. Thus for example the above list of digits will most certainly occur somewhere
in the expansion of ; although be warned that it might take us a long time to come across
them.
Our rst introduction to the essential nature of e is through a problem in commerce -
Suppose that we have $1 and invest it such that at the end of the rst year we are given $1
interest, thats right 100% interest. Now if instead of the full year we invested our money for
1
10
year it would be only fair that we would get
1
10
of this amount in interest, that is we would
have a total amount of $1 +
1
10
at the end of this rst period.
With simple interest we would get this same amount for each of our ten periods and hence
the total interest would be $1 which is the same as in the rst case. This is the same as saying
that our money has grown from $1 to $2 by an amount of $10
1
10
.
2
Now if we compound our interest the following happens. At the end of the rst
1
10
of
a year we would get $(1 +
1
10
) as before. However at the beginning of the second period
we would invest this amount instead of $1. It follows that at the end of the second period
we would have $(1 +
1
10
)
2
= 1.21 and so on until at the end of the year we would have
$(1 +
1
10
)
10
= $2.5937. Similarly at the end of the year with 100 periods we would have
$(1 +
1
100
)
100
= 2.7048.. and for 1000 periods $2.71692..
From this we nd experimentally that as the number of periods increases the total amount
appears to approach $e. We note that the derivative of (1 +
1
n
)
n
is
(1 +
1
n
)
n
{ln(1 +
1
n
)
1
(1 + n)
}
and that this is positive and decreases to zero as n increases. So this is further evidence that
(1 +
1
n
)
n
converges as n becomes unbounded. However we require the Main Theorem to
prove this conclusively.
You might like to check the following gures for reinvesting on a monthly ($2.63), weekly
($2.69) and daily ($2.71) basis. Even if we reinvest at the end of a half year the amount is
$2.25 which is 12.5% better than simple interest.
This is one of the reasons why banks make such a prot from simple minded customers.
In general if instead of an increase of $1 we had an increase of $r for the year (r is called
the effective annual rate of interest) we can readily show in a similar way to the above that,
with an initial investment of $1, the amount at the end of the year is (1 +
r
n
)
n
and that this
expression approaches e
r
as n increases without bound (see A4, below).
For a rather dramatic example of the absurdity of believing that any economic system can
really sustain exponential growth at even a reasonable" growth rate see Appendix D.
Biological Growth Models
A similar exponential process occurs during the initial growth of many biological organ-
3
isms. This occurs whenever the rate of change is proportional to the size or the number of
the organisms. In this case we nd that for example if y is the length of a root or the number
of cells in an organism then
y = y(0)e
rt
, (1)
where y(0) is the initial value , t is time and r is the intrinsic rate of increase.
This relationship is what is nowgenerally referred to as Malthusian or exponential growth.
However the original growth process referred to by Malthus was actually geometric. Knowing
that
de
x
dx
= e
x
, (from proposition Q4), it is seen that if we differentiate this equation on each
side we get
dy
dt
= y(0)re
rt
= ry. (2)
Hence
r =
dy
dt
y
. (3)
Thus r, which we call the intrinsic rate of increase, is the actual rate of increase
dy
dt
divided
by the total length of, or number of, organisms.
We note that if y were say the number of individuals in a population then
dy
dt
y
, and hence
r, is seen to be the average rate of growth per individual. It is a remarkable fact that, in a
wide range of growth processes in nature and economics, this quantity is constant during the
initial stages of such growth.
From (3) we see that
ln
_
y
y(0)
_
= rt, (4)
where from now on ln x means log
e
x.
(4) may be used as a rough check to test that we have exponential growth. We simply plot
the logarithm of
y
y(0)
against t. If this is a straight line then we have exponential growth.
This exponential growth is hard to deal with intuitively and in practice it is easier to use
the doubling time t
d
dened as the interval of time required for y to double. Then if we
4
consider the population at any time t we will double our population at time t + t
d
. Hence
y(t + t
d
)
y(t)
= 2.
Substituting in our exponential expression (1) we see that this reduces to
exp(rt
d
) = 2, (that is the y(t) cancels out and the expression is thus independent of t).
From this we have
t
d
= ln
2
r
.
If r is measured as a % then this is equivalent to the following - the doubling time is roughly
equal to 70 divided by the interest in % or
td
70
r%
.
This hyperbolic (or reciprocal) relationship is considered to be easier to grasp than the expo-
nential. Thus for example with r = 7% the doubling time is 10yrs, and with r = 10% the
doubling time is 7 yrs.
Other Properties and History
e has all sorts of other interesting properties. Perhaps the most beautiful relationship is
Eulers formula which also involves and i the square root of 1. This is
e
i
= 1.
Part of the beauty of this relationship lies in its unexpected simplicity. Here it is seen that
i brings together e and and unity, three of the fundamental constants of nature, together
with the two operations minus and square root.
We will show in detail in Appendix A that Napiers logarithms are somewhat different to
natural logarithms which use base e. However through formula AA there is a relationship but
Napier was quite unaware of it so his work only implied the use of e.
5
It was not until Bernoulli (Johann I, 1704) that the symbol e and the name exponential
were rst used. This arose because it was necessary to have a way of representing e
x
. In fact
the term Naperian antilogarithm" or e was not used until 1784.
The calculation of logarithms by innite series was carried out by Gregory, Wallis and
Halley (of Comet fame) beginning in the 1670s. It was not until much later (1742) that
William Jones (1675-1749) used exponents to calculate logarithms and allowed present day
methods tobe developed. However it shouldbe emphasisedthat Euler was the rst tointroduce
this concept, and the symbol e, for the base of natural logarithms in an unpublished paper of
1728 (he quoted e to 27 places).
The Layout of this Note
The formal part of this note begins with 2 denitions rstly of a function which is the area
under the curve
1
x
and then a denition for the inverse of this function.
This is followed by a set of propositions (P1, P2 etc) about the logarithmic function
which are developed axiomatically from the denition. Then we present some properties
of the exponential function e
x
(Q1, Q2, Q3) developed axiomatically from its denition as
an inverse of the exponential. Then we show that e
x
is equal to its derivative (Q4). This is
followed by some properties of the functions a
x
and x
a
(Q5, Q6).
Q7 gives a derivation for e which we used in the interest rate problem. Then we derive an
expression for e
x
as a series (Q8). From this we can show that e is irrational (see Appendix
B). Then we give an independent proof of Q7 based on a note by Barnes which has the merit
that it is coupled to Eulers constant. At this stage we digress into Harmonic numbers.
Finally in contrast to the above axiomatic development we offer some alternative proofs
A1, A2, ... A6. These have been selected simply on the basis of their elegance - and by
elegance is meant unexpected simplicity.
6
An Axiomatic Approach
We begin with the denition
LN(x) =
_
x
1
t
1
dt. def 1
It is seen that LN(x) is simply a function of x at this stage. We hope to show that it does
have the properties of a logarithm (see P4 below).
We also dene the inverse of LN(x) as E(y)
LN(x) = y
Hence
E(y) = LN
1
(y) = x. def 2
We hope to show that this inverse function is actually e
y
(see Q3).
Properties of the logarithmic function ln x
We now show that we can derive what we normally consider to be the properties of the
logarithmic function from the above denition, def 1. We present these as a series of lemmas
(minor theorems) or propositions - (it follows that a dillema is a theoremproduced by a stupid
person).
It is not until P4 that we are convinced that LN(x) is in fact ln x.
P1 LN(1) = 0, directly from def 1
P2
dLN(x)
dx
=
1
x
.
This follows directly from def 1 by differentiation under the integral.
P3 The additive relation holds, LN(ab) = LN(a) + LN(b).
7
Proof Let z = LN(ax) then from the chain rule
dz
dx
=
_
1
ax
_
a =
1
x
.
However
d LN(x)
dx
=
1
x
fromP1. Hence as the functions LN(ax) and LN(x) have the same
derivative they must differ only by a constant, that is
LN(ax) = LN(x) + c.
In particular this relationship holds for x = 1 and thus c = LN(a), hence
LN(ax) = LN(x) + LN(a). With x = b we have proved our lemma.
P3a From P3 it follows directly that
_
xy
1
t
1
dt =
_
x
1
t
1
dt +
_
y
1
t
1
dt,
a result which is rather messy to prove directly by integration (is t
1
the only function that
has this property ?)
P4 LN(a
n
) = nLN(a).
This follows by continued application of P3. By substituting
x
2
+ y
2
. Hence the result follows. Thus
the graph of e
x
is simply the graph of ln x reected about the line y = x. Thus for example
we see that the slope of e
x
at x = 0 is the same as the slope of ln x at x = 1, namely unity.
Q5 The functions a
x
and x
a
.
Let ln e
f(x)
= g(x), then by (Q3b) e
f(x)
= e
g(x)
and thus f(x) = g(x). Similarly if
h(x) = e
ln f(x)
, taking logs we have ln h(x) = ln f(x), and again f(x) = h(x). Combining
these we see that
f(x) = ln e
f(x)
= e
ln f(x)
.
From this we see that we can always write a as power of e or as a logarithm
a = e
ln a
= ln e
a
.
We can now write our functions very neatly as
y = a
x
= e
xln a
and x
a
= e
a ln x
.
10
Also more directly from the series for e (Q8)
y = e
xln 2
1 + x ln 2
Q6
d(a
x
)
dx
= a
x
ln a and
d(x
a
)
dx
= ax
a1
.
Q7 e
z
= lim
n
_
1 +
z
n
_
n
If f(x) = ln x, then
f
(x) =
1
x
= lim
h0
f(x + h) f(x)
h
,
= lim
h0
1
h
{ln(x + h) ln(x)},
= lim
h0
1
h
ln(1 +
h
x
).
With
1
x
= z
lim
h0
1
h
ln(1 + hz) = z.
Consequently
e
z
= lim
h0
exp[ln(1 + hz)
1/h
] = lim
h0
(1 + hz)
1/h
= lim
n
(1 +
z
n
)
n
, if h =
1
n
.
In particular with z = 1 we have independently shown that
e = lim
n
_
1 +
1
n
_
n
Q8 e
x
= 1 + x +
x
2
2!
+
x
3
3!
+ . . . =
x
j
j!
, for j 0.
This follows directly from the Taylor series expansion about x = 0 for f(x) = e
x
. From
this we have for x = 1
e = 1 + 1 +
1
2!
+
1
3!
+ . . . =
1
j!
, for j 0,
11
which has reasonably fast convergence and so may be used to calculate e
An example Suppose we seek an approximation to y = 2
x
for small x. Fromthe binomial
theorem
y = (1 + 1)
x
= 1 + x +
x(x1)
2!
+
x(x1)(x2)
3!
+
x(x1)(x2)(x4)
4!
+ . . .
= 1 + x(1
1
2
+
1
3
1
4
+ . . .) + x
2
(
1
2
. . .)
1 + x ln 2.
Q9 With Q8 as our denition for e we can now give the
Proof that e is irrational. Assume the opposite that is e =
p
q
where p and q are integers.
Then
p
q
= 1 + 1 +
1
2!
+ . . . +
1
q!
+
1
(q+1)!
+ . . .
(q 1)!p =
_
q! + q! +
q!
2!
+
q!
3!
+ . . .
q!
q!
_
+
1
(q+1)
+
1
(q+1)(q+2)
+ . . .
But the expressions on the left hand side and in brackets are obviously integers, whereas
the remainder
1
(q + 1)
+
1
(q + 1)(q + 2)
+ . . . <
1
(q + 1)
+
1
(q + 1)
2
+ . . . =
1
(q+1)
1
1
(q+1)
=
1
q
1.
Thus for our original assumption to be valid we require an integer to be equal to an integer
plus a number less than one. As this is absurd our original asssumption is invalid and e must
be irrational.
Theorem M To show that
e = lim
n
_
1 +
1
n
_
n
.
This presentation is based on Barnes, C. (1984). As a by-product it establishes the
expression for Eulers constant , which is dened below. Alternative proofs are given in Q7,
and A1.
Proof From P6
_ 1
n
1
(n+1)
ln xdx =
1
n(n + 1)
_
ln
(n + 1)
n
n + 1
1
_
. (M1)
12
Also we know from the mean value theorem that there exists a c
n
somewhere in the range
of integration (that is
1
(1+n)
< cn <
1
n
) where
_ 1
n
1
(n+1)
ln xdx =
_
1
n
1
(1 + n)
_
ln c
n
. (M2)
From (M1) and (M2)
ln
_
(n + 1)
n
n
n+1
_
1 = ln c
n
,
from which
ln
_
_
(n+1)
n
n
n+1
c
n
_
_
= 1.
Hence
e =
(n+1)
n
n
n+1
c
n
,
or
e =
_
1 +
1
n
_
n
1/nc
n
. (M3)
Now let a
n
=
1
(nc
n
)
. Then because of the range of c
n
we have
1 < a
n
< 1 +
1
n
. (M4)
Hence we have the remarkable result that for all positive n
e = a
n
_
1 +
1
n
_
n
. (M5)
As a
n
> 1 it follows immediately from (M5) that for example
e >
_
1 +
1
2
_
2
= 2
1
4
.
Also because a
n
is in the above sandwich (M4), in the limit as n , a
n
= 1 and thus
our theorem is proved.
Denition We dene Eulers constant as
= lim
n
__
1 +
1
2
+
1
3
+ . . . +
1
n
_
ln(n + 1)
_
.
13
Corollary M
<
2
6
= 1.64..
Proof
From (M5)
e = a
n
_
1 +
1
n
_
n
=
a
n
(n + 1)
n
n
n
.
Thus taking logarithms
1 = ln a
n
+ n{ln(n + 1) ln n},
or
1
n
= ln a
1/n
n
+ ln(n + 1) ln n.
From this by substituting successive values 1, 2, 3, ... for n, we obtain the equations
1 = ln a
1
+ ln 2 ln 1
1
2
= ln a
1/2
2
+ ln 3 ln 2
1
3
= ln a
1/3
3
+ ln 4 ln 3
..........
1
n
= ln a
1/n
n
+ ln(n + 1) ln n.
Summing
1 +
1
2
+
1
3
+ . . . +
1
n
ln(n + 1) = ln a
1
a
1/2
2
a
1/3
3
. . . a
1/n
n
. (M6)
We note that the continued product on the right hand side is increasing in n. Also from
(M4) a
n
< 1 +
1
n
, hence
ln a
1
a
1/2
2
a
1/3
3
. . . a
1/n
n
= ln a
1
+
1
2
ln a
2
+ . . .
1
n
ln a
n
,
14
< ln
_
1 +
1
1
_
+
1
2
ln
_
1 +
1
2
_
+ . . . +
1
n
ln
_
1 +
1
n
_
.
Furthermore from (M4), a
n
> 1. Thus from (M5), 1 <
_
1 +
1
n
_
n < e. Hence ln 1 <
ln
_
1 +
1
n
_
n
< ln e, or 0 < ln
_
1 +
1
n
_
< 1/n. Substituting this into the inequality
ln a
1
a
1/2
2
a
1/3
3
. . . a
1/n
n
< 1 +
1
2
2
+
1
3
2
+ . . .
1
n
2
.
But it can be easily shown that the innite series 1 +
1
2
2
+
1
3
2
+ . . . is convergent. In fact
it has the value
2
6
= 1.644 (see Knuth - would you believe that if we select two integers at
random the probability that they dont have a common factor is
6
2
). Thus as the continued
product is bounded, the l.h.s. of (M6) is bounded and hence < 1.64 . . .
Note - It can be shown numerically that = 0.5772 . . ..
Eulers constant features widely in many diverse areas of mathematics such as number
theory and modied Bessel functions. Also Knuth (1973) (with a rather more introductory
account of this in Graham, Knuth and Patashnik) gives an interesting treatment of the above
series for the summation of the rst n reciprocals which he denes as
H
n
= 1 +
1
2
+
1
3
+
1
4
+ . . . +
1
n
=
1
k
, k n,
where H stands for harmonic number. Hence the innite series H
(r) is convergent. H
(r) is the corresponding innite series is the Rieman Zeta function and is usually
represented by (r). Now if we attempt to solve (r) = H
(r)
1
n
2
_
1 +
1
n
_
1
n
2
,
= 1,
and this shows that y = e as required.
Notice that we have assumed P2, that is that we know the derivative of a logarithm.
A2 Let f be a function such that for every pair of reals x, y
f(x + y) = f(x)f(y). (M).
We will say that f satises the multiplicative property.
Then we require to prove the following:
if f(x) is continuous and is not identically zero, and f(1) = a > 0, and furthermore f(x)
satises the multiplicative property then
f(x) = a
x
for all x.
17
Proof Now f(x) = f{(x 1) + 1} = f(x 1)f(1).
Hence if f(1) = 0 then f(x) = 0 for all x. Hence we assume the opposite.
Suppose f(1) = a > 0. Then again
f(x) = f(x 1)f(1) = f(x 1)a = f(x 2)a
2
= f{x (x 1)}a
x1
= f(1)a
x1
= a
x
.
Note that this shows that (M) holds only when x is integer as we subtracted integers from
x to obtain 1. To prove this for x rational say
m
p
where m and p are integers we rstly write
1 =
1
p
+
1
p
+ . . .
1
p
, that is p terms.
Thus applying (M) repeatedly
f(1) =
_
f
_
1
p
__
1/p
.
Also
f
_
m
p
_
= f
_
1
p
+
1
p
+ . . .
1
p
_
, where we now have m terms
=
_
f
_
1
p
__
m
= {f(1)}
m/p
= a
m/p
or f(x) = a
x
, as required for x rational. It follows that as the rationals are everywhere dense
that this result also holds for real functions provided that f(x) is continuous.
A3 We wish to show that the solution to
dy
dx
= y where y(0) = 1,
has the multiplicative property of A2.
That is we claim that y(x + v) = y(x)y(v).
Proof Let z(x) = y(a + x). Then z
(x) = z(x).
Now
d
dx
y(x)
z(x)
=
zy
y
z
2
=
zyzy
z
2
= 0.
18
Hence
y(x)
z(x)
= K, a constant, or
y(x) = Kz(x) = K y(a + x).
Hence as y(0) = 1, with x = 0 we have K =
1
y(a)
. From this
y(a + x) = y(a) y(x),
as required.
We also know independently that y(x) has a unique solution and hence the multiplicative
property is unique.
A4 Assuming A1, e = lim
t
_
1 +
1
t
_
t
, we wish to show that
e
x
= lim
n
_
1 +
x
n
_
n
.
We note rstly that t may be real, that is it does not have to range over the integers for A1
to be valid.
Now
_
1 +
x
n
_
n
=
_
_
1 +
x
n
_
n/x
_
x
=
_
_
1 +
1
t
_
t
_
x
with t =
n
x
. The result then follows by taking the limit.
A5 Let f(x) = b
x
Then from Q6
f
(x) = b
x
ln b and thus f
(0) = ln b
But from rst principles
f
(x) = lim
h0
b
x+h
b
x
h
.
Thus with x = 0
ln b = lim
h0
b
h
1
h
.
19
The relationship A5 automatically resolves the apparent conict with the expression
_
b
a
x
u
dx =
b
u+1
a
u+1
(u + 1)
,
when u 1, say 1+h, and with a = 1, it is seen fromthe above limit that this approaches
ln b as required when h 0
A6 Using Q8 to dene e
x
we wish to establish the multiplicative property.
Now
e
u
e
v
=
m=0
u
m
m!
m=0
v
m
m!
But it is known that (the Cauchy convolution form for the product of two series)
n=0
a
n
n=0
b
n
=
n=0
c
n
where c
n
=
n
k=0
a
nk
b
k
.
Hence
e
u
e
v
=
n=0
n
k=0
u
nk
(n k)!
v
k
k!
,
and by the binomial expansion
(u + v)
n
=
n
k=0
n!u
nk
v
k
(n k)!k!
,
from which
e
u
e
v
=
n=0
(u + v)
n
n!
= e
u+v
.
References
Barnes, C. Eulers Constant and e", Am Math Monthly Aug 84, p428.
Coolidge, J.L. The Mathematics of Great Amateurs", (Oxford, 1949).
20
Eves, H. Great Moments in Mathematics", (Mathematical Association of America).
Knuth, D.E. The Art of Computer Programming", (Section 1.2.7 Harmonic Numbers p73)
vol 1 (2nd edition) 1973.
Graham, Knuth, D. and Patashnik, O. Concrete Maths" - (Addison Wesley - 2nd Edition
1994) - Concrete is short for continuous plus discrete).
Maor, Eli. e: The Story of a Number" (Princeton University Press, 1994).
Struik, D.J. A Source Book in Mathematics 1200-1800", (Harvard, 1969).
APPENDIX A
Napier (1550-1617), a Scottish Laird, was one of the great amateurs of mathematics.
He began to develop logarithms in 1594 and worked on this task until his death. Amongst
other things he was obsessed with demonstrating that the then reigning pope was the anti
christ. In fact Napier believed that his now forgotten work on the Apocalypse of St John
was his greatest service to mankind and this belief was bolstered by the fact that his book on
this subject went through 20 editions. Newton had a similar obsession with the Apocalypse
which might indicate (at least to a statistician) that if we started off budding young geniuses
on this study it could lead to great mathematical breakthroughs.
Briggs was the rst occupant of a chair of mathematics in Britain. This was in Geometry at
Gresham College in London. Later he was the rst to occupy the Savillian chair of geometry
at Oxford. He extended Napiers work to base 10 and went to visit him after hearing of his
wonderful discovery and communicating with him. On the interminable journey by coach
to Edinburgh Briggs recorded some of his thoughts in his diary - how tall a forehead must
his nobleman possess in order to house the brains that have discovered so remarkable an
invention and so on. His arrival in Edinburgh was very much delayed and Napier confessed
to a friend Ah John, the Professor will not come".
21
At that very moment a knock was heard at the gate and the Professor was ushered into
the noblemans presence. For almost a quarter of an hour each man beheld the other without
speaking a word. Then Briggs said My Lord, I have undertaken this long journey purposely
to see your person, and to learn by what engine of wit or ingenuity you came rst to think of
this most excellent help in astronomy. But, my lord being by you found out, I wonder nobody
found it out before when now known it appears so easy."
In understanding the importance the work of Napier it should be realised that the math-
ematics of his day was very decient in both concepts and notation. Not only was there
no concept of logarithm but the concept of powers was only emerging and for example the
symbol x
n
as a shorthand for x x . . . had not been introduced. It is of interest to note that
only because of the authority of Napiers work did the decimal notation became generally
accepted, thus = 2.083 instead of the clumsy 2
83
1000
.
The description of the technique he used for calculating his logarithms" and tables of such
logarithms" are best described in his second book which was rst published in Latin in 1619.
This was translated from the Latin into English and published by Blackwood, Edinburgh in
1889, with the title The Construction of the Wonderful Canon of Logarithms".
Curiously enough although the term logarithm occurs in the title the text uses only the
term articial number. As we will show shortly his logarithm is not what we now know as
the natural logarithm. Detailed accounts of his work are contained in Struik, D.J. (1969) and
Coolidge, J.L. (1949). Unfortunately his presentation is almost unintelligible to a present
day mathematician and although his work is one of the great moments in mathematics it is in
effect an ugly dead end, similar to the dead end taken in the pursuit of the Ptolemaic system
of epicycles as a description of planetary motion.
The basis for his idea was probably seen in the formula
sin Asin B =
1
2
{cos(A B) cos(A + B)},
which was certainly known in Napiers time.
22
Now Stifel had earlier in Arithmetica Integra" (1543) recorded the rst observation that
the terms of the geometric progression 1, r, r
2
, . . . correspond to the arithmetic progression
0, 1, 2, . . .. Stifel also noted that multiplication of two terms in this g.p. yields a term
whose exponent is in effect the sum of the corresponding two terms in the a.p. Similarly
he implied that division was coupled with subtraction of exponents and he was even able
to manipulate negative and fractional exponents. All this was recorded without any special
mathematical notation. It appears that Napier did know of this work but, incredibly, failed to
see the relevance of exponents on his logarithms. If he had picked it up it would most surely
have accelerated the development of mathematics.
Because of his interest in spherical trigonometry he based his work on the logarithms of
sines however let us warn you now that his use of the terms radius and sine simply referred to
intervals on a straight line (AZ and DZ in the following). In effect he said the following. Lay
out two lines AZ and A
, C
, D
, E
corresponds to
his logarithm of BZ and so on.
Now he has a point initially at A moving along AZ, with a velocity proportional to its
distance from Z and this corresponds to a point moving along A
, A
, . . . to be 1,2,3,
. . ., however any arbitrary increasing a.p. would be just as satisfactory.
Translating the above into present day calculus, which was of course unavailable to Napier,
we consider the motionof twoparticles the rst whichwe will call P, travellingalongAB. The
second particle P
. Let B = a(= 10
7
),
PZ = y and let AP