HP Calculator
HP Calculator
HP Calculator
In the Spotlight
Issue 25
» 30th Anniversary Edition of the HP 12c October 2011
Performance never goes out of style.
Celebrate 30 years of this one of a kind Welcome to the twenty-
business calculator with a limited edition fifth edition of the HP
collector's special. Solve newsletter. Learn
calculation concepts, get
advice to help you
Your articles succeed in the office or
the classroom, and be
the first to find out about
new HP calculating
solutions and special
offers.
» Download the
PDF version
» The HP 15c Limited Edition » The HP-12C, 30 Years and
of newsletter articles.
calculator Counting
One of nine machines in the Richard J. Nelson and Gene Wright
Voyager Series of HP » Contact the editor
Read this extensive review of
calculators, HP announces a how the HP-12C calculator From the Editor
re-release of their legendary has withstood the test of
HP 15c calculator. professionals and time
longer than any other
calculator.
HP respects your privacy. If you'd like to discontinue receiving e-mails from HP regarding special offers and information,
please unsubscribe here.
For more information regarding HP's privacy policy or to obtain contact information, please visit our privacy statement or write
to us at: HP Privacy Mailbox, 11445 Compaq Center Drive W., Mailstop 040307, Houston, Texas 77070, ATTN: HP Privacy
Mailbox.
© 2011 Hewlett-Packard Development Company, L.P. The information contained herein is subject to change without notice.
30th Anniversary Edition of the HP 12c
HP Solve #25 page 3
From the Editor
Article – Next
Return to Top
HU U
HP Solve # 25 Page 4 Page 1 of 1
PERFORMANCE NEVER GOES OUT OF Style
HP is celebrating the 30th Anniversary of its iconic gold trimmed financial calculator, the HP 12c.
With its sophisticated look and feel, this calculator is still in its heyday, withstanding the test of time
for professionals in the world of banking, finance, and real estate.
On May 4, 2011 the Wall Street Journal referred to the HP 12c as the “Cult Calculator” asserting that
for finance professionals, “Ignorance of the 12c can flash more warning signs than a scuffed pair of
shoes.” This article became HP’s inspiration for the 12c 30th Anniversary ad campaign which will
run in numerous financial publications word-wide. The ad portrays a well-dressed, sophisticated
business man who has successfully advanced his career from a broker fighting his way through the
bull pen on Wall Street to a talented financier in the boardroom of America’s finest corporations—all
with the assistance of his trusted companion—the 12c. The WSJ article states: “the 12c is just as
commonplace for financial analysts as papers and cell phones” demonstrating that it has become more
than just a tool of the trade; it’s a fundamental staple.
Exclusive to the 30th Anniversary, HP is offering this exemplary calculator in a limited edition
collector’s special. Each 12c comes in a leather gift box, is embossed with a “30th Anniversary
Edition” emblem, and is stamped with a unique production number, making it the perfect gift for
business professionals, students, and calculator connoisseurs.
The 12c “30th Anniversary Edition” is available for purchase September 1, 2011 in retail stores and
online at: www.hp.com/go/voyager
*Peterson, Kristina. “A Cult Calculator Endures.” Wall Street Journal, May 4, 2011: Print
http://online.wsj.com/article/SB10001424052748703841904576257440326458056.html
This decision devastated 15c loyalists. Since then, they have beaten their chests loudly for HP to
bring it back, creating websites such as, “Bring Back the HP15c”.
Their passion for the dearly departed 15c was so strong that over 15,000 people in 175 countries
have signed a petition—“HP, bring back the 15c.”
With demand on this petition exceeding 84,000 units, bidding wars on Ebay reaching over $500,
and tireless forums of fans lamenting the loss of “The long gone but not forgotten” 15c, HP
relented—“It’s time to resurrect the legendary 15c!”
The 15c PR campaign will commence on September 1, 2011 and be featured in numerous
engineering and scientific publications world-wide. The ad is a play on the Voyager Calculator
Series and the Voyager Space Mission— both symbols of the technology revolution during the
early 1980s. The Voyager Space Shuttle discovers the 15c on the outer edges of the Milky Way,
floating further and further away from planet earth and its devotees. The Golden Record,
intended to relay messages from human life to extra terrestrial life, communicates the sentiments
for the return of the 15c.
Today, fans everywhere can celebrate the return of the HP 15c with this limited-edition
collector’s special. Each 15c comes in a leather gift box, is embossed with a “Limited Edition”
emblem, and is stamped with a unique production number, making it the perfect gift for 15c
aficionados!
The 15c Limited-Edition Calculator will be available for purchase September 1, 2011 in retail
stores and online at: www.hp.com/go/voyager.
Introduction
HP has a long list of calculator firsts(1) and the 1981 HP-12C
business, finance, and real estate calculator significantly
contributes to the list. The 12C has been continuously made
for 30+ years. Made in five countries(2), and often in two
countries at the same time, HP-12Cs made today look pretty
much the same as they did 30+ years ago. The HP-12C
calculator has withstood the test of professionals and time
Fig. 1 – HP-12C 30 Anniversary Edition. longer than any other calculator (and possibly longer than
th
The financial HP-12C was introduced September 16th 1981 along with a scientific HP-11C. These two
machines were eventually joined by six other similar looking calculators as part of the Voyager series of
HP calculators. See Fig. 2 below. Also known as the “Slim Line(3) ” calculators they were quite different
from other calculators in HP’s line up because of their smaller size and thinness. Some sellers/reviewers/
websites also refer to the voyager calculators as the “10 Series.”
The HP-12C integrates the five Time Value of Money, TVM, variables into one solver that provides the
easy and very accurate solution of n, i, PV, PMT, and FV using a consistent approach to solving complex
compound interest problems.
The HP-12C has a clear readable display: The popularity of the HP-12C was debated at a user’s
conference in London in 2002. The intent of the one hour discussion was to identify the most important
features of the HP-12C that makes it so popular. Six qualities of excellence were given to start the
discussion. Everyone had an opinion, but after an hour one quality was mentioned over all others; a clear
readable display(4). The HP-12C uses a 10 digit 7 segment Liquid Crystal Display, LCD, with 0.2 inch
high digits.
The HP-12C has a clean uncluttered keyboard: The HP-12C is a tool to be consulted in the “heat of the
deal” for the business community. It must not appear intimidating to either the customer or the user. The
gold display bezel and the gold frame around the keys give the calculator a timeless quality look. The
gold is in contrast to the textured dark brown case. Of course the classical HP calculator keyboard key
feel and tactile feedback click is also an important part of the keyboard. Two shift keys, yellow and blue,
are used to provide access to the 120+ functions of the calculator.
The HP-12C uses a horizontal format: Most of the calculators HP makes use a vertical format. The
Voyager series of calculators, however, broke from the tall thick vertical cases of other models to use a
thin horizontal format. This was quickly accepted for two reasons. First, the intended user was often in
the office and at a desk. Second, the light weight, 4.1(5) oz., makes it easy to hold in the hand. Holding a
machine in your two hands and operating the keyboard with your thumbs (pre-cell phone) is natural for a
12C if used while standing.
The HP-12C has a powerful business function set: The HP-12C has a 2,465 year (1582 – 4046) calendar
built in with a nice set of date functions. It also has an excellent mix of business and statistical functions
with the normal [no trigonometry(7)] math functions. Specific applications such as IRR and NPV are also
included. See Appendix A for the details.
The HP-12C uses RPN: Most business calculator users have no idea what RPN means, but they easily
learn the beauty of simply pressing a key to execute a function. You want to add two numbers you press
The HP-12C is programmable: See article on the importance of calculator programmability elsewhere in
this issue. As with RPN most business users of the HP-12C are not very concerned with its program-
mability. When they have a number of similar problems to solve, however, they will discover how RPN
programmability may be easily learned to “get the work done.” The programming capability is limited to
99 program steps and two logic compares (X≤Y? & X=0?). Many of the nice programming features of
more capable calculators are missing: absolute value, flags, increment or decrement for looping, indirect
addressing, labels, sign, subroutines, etc. This feature is further illustrated with the program examples in
appendix C.
The HP-12C is accurate: Most people who use an HP calculator expect that the answer is correct. The
situation is that this is not the case in the real world, and HP calculators are renowned for their accuracy.
There are two issues of accuracy involved in the HP-12C. First is the normal issue of making calculations
accurate rounded to ± 1 count in the 10th place. The second accuracy issue is that of the applications
programs, especially the TVM. Solutions are iteratively solved and the process may be error prone. The
algorithms used are critical and the HP-12C uses the best. The critical feature (and characteristic) of the
HP-12C is that a problem solved on the first (1981) calculator will get the same answer on the most recent
(2011) calculator. This is quite remarkable in the calculator world – and another first.
While the insides have changed and the outside looks much the same there are subtle mechanical (case)
changes that have been made over the years. HHC 2011 attendees will be asked to bring their HP-12Cs to
the conference so that we may photograph them.
The HP-12C has been in continuous production for 32 years. Research was made and while this is far
from being 100% accurate it illustrates how things have changed over three decades. The HP-12C
production was in the countries listed in Table 1 below. This data was obtained from interviews with
retired HP personnel, serial numbers from friends and the Internet, and Don O’Rourke of
http://internationalcalculator.com/ The serial numbers of HP-12Cs made in China have changed in the
last few years and serial number changes and inconsistencies have been observed by the HP User
Community. One conclusion that may be made is that HP had a significant overlap of countries whenever
a country change was made to insure a continuous supply.
In addition to the technology the competition has changed greatly over the last 30 years and this has
influenced the price(10) of the HP-12C. The current HP website price is $70, less than half of the original
introductory price of $150.
Pressing [FV] should display one year’s accretion, but different financial calculators display different
amounts. Here is the result if you use an HP-27, HP-92, HP-37, HP-38, and HP-12C = $ 331,667.0067.
As an aside, the new HP 10bII calculates $ 331,667.006691.
Less accurate, non-HP, calculators may get answers like: $ 293,539.16035, $ 334,858.18373, or $
331,559.383549.
Even with the price increases in 1989 and 1996 the HP-
12C popularity did not change. HP suggested list price history for the HP-12C.
(11) William Kahan of the University of California Berkeley consulted with HP on the HP-12C and other
calculators (HP-15C & HP-34C) to insure algorithmic accuracy. He has a long and productive career in
Fig. D1 – Handy sized Owner’s Fig. D2 – Solutions handbook Rev. F Fig. D3 – Real Estate Applications
Handbook Rev. G 11/85. 7/87. Handbook Rev. B 3/84.
Early third party support booklets are represented below. Chary is still active in HP-12C support using
videos. D3 & D4 - http://hp-12c-calculator.com/?page_id=60
Gene’s business math book is at: http://www.amazon.com/Quantitative-Analysis-Business-Statistics-Calculator/dp/1888840382/
Fig. D3 – Chary Software Program- Fig. D4 – Chary Software Program- Fig. D6 - Renaissance Publications.
ming Hints book Rev. A 12/85 ming Hints book Rev. A 8/86 Problem examples book © 1985
Return to Top
Previous Article – Next
U UH HU U
Calculator Programmability – How Important is it in 2011?
Richard J. Nelson
Introduction
When the non-programmable scientific HP-35A calculator first appeared many HP
enthusiasts wrote down keystroke lists(1) on index cards for the most efficient method of
solving problems, especially those that required iterating to a solution. Two years later
the programmable HP-65A appeared and tens of thousands of HP calculator users
became programmers. This was in the mid 70’s and the HP-65A was the first personal
computer that gave programmability to everyone, not just to college students or
employees of large corporations that could afford room sized computers.
A program may also “manage” a data base and serve as a ‘look up” mechanism for information such as
prices, specifications, or any numeric data in the form of tables or equations.
Programming Advantages
What are the advantages of programmability for the calculator user?
1. Programming provides a quick simple method of solving the same problem over and
over again. A single keystroke could replace dozens of keystrokes. An example is
calculating the volume of a dozen grapefruit from a quality sampling.
Fig. 3 - Languages 2. A program could solve more complex problems that the user couldn’t solve
otherwise. The program is written by an expert and used by anyone. An example is a complex
statistical calculation such as mean and standard deviation.
3. A program could contain data that isn’t conveniently looked up or remembered. Making English
Metric conversions is a simple example of this program usage e.g. converting grains to grams or troy
ounces to grams.
HP 15c LE programming
The HP 15c RPN program instruction set has a respectable collection
of programming features. The large number of applications
programs, especially those that may be adapted from many other
Fig. 5 – HP-15c+ Program Mode. RPN programmable calculators is a substantial advantage
the exceptional speed of the HP 15c LE provides an additional advantage of this programmability.
More complex problems are often solved on a personal computer, which tends to be
turned on for extended periods, such as a work week of five days. Comparing the
Fig. 7 – Tools. calculator solution with a computer solution doesn’t yield a clear advantage in turn-on
convenience (instant on vs. boot time). Still, a single key press for a programmed solution is nearly
impossible to improve and the usefulness of programmability is easily justified in the midrange to the
high end models for most of the dozen reasons listed above.
The programming method, however, must be as simple as possible so the average user may easily master
it. The HP 15c LE with its FOCAL like keystroke programming is easy to learn. Even if the student is
less interested in programming their own machine he or she will still need those programs suitable for the
course at hand. Being able to share programs via the Internet and down load them via a USB port will
always be a justification to include the programming feature of any serious scientific or graphing
calculator.
The large number of features, memory size, and super speed of the HP
15c makes its programmability an important feature when compared to
other calculators - for the 12 programming reasons listed above.
To answer the title question, “Calculator Programmability – How Important is it Today?” I would have to
say that it is as important as it ever was. A more modern consideration is having an I/O capable of down
loading the programs regardless of the language of the operating system and even if most users are not
writing programs themselves.
Fig. A1 – Owner’s Handbook (initial1982 version) and Advanced Functions handbook (later 1984 version).
What fraction is π equal to? Most people would reply “22/7” and of course they would be wrong, since
π is an irrational number which means that it cannot be expressed as a ratio of integers. But 22/7 is
close to π . Is there a fraction that is closer?
The answer originally flowed from the following line of brute-force observation: “ π ≅ 3 . Well, not quite
1 1
3, but more precisely π ≅ 3 . (Hence 22/7). Well, not quite 7, but 7 , which means that
7 15
1
π ≅ 3+ . (Hence the better approximation 333/106).” And so on. Following this line of reasoning,
1
7+
15
closer and closer approximations to π can be generated. The trick is to pick the correct adjustments.
There is a method for doing this that avoids brute-force searching. It is called “The Continued-Fraction
Algorithm.” Here it is.
This will generate the following list for π : { 3, 7, 15, 1, 293 }. Recognize the first three numbers? They
are the denominators printed in red in the large equation for π above. Try using the first four numbers
from this list instead of just the first three; it yields this larger expression:
1 355
π ≅ 3+ =
1 113
7+
1
15 +
1
Hence a better approximation for π is 355/113. Now use all five numbers from the list:
1
π ≅ 3+
1
7+
1
15 +
1
1+
293
This huge expression of nested reciprocals (called a “continued fraction”) simplifies to 104348/33215,
which is therefore a fraction that is even closer to π than 355/113. It is so close, in fact, that it evaluates
on the HP-15 to the same value as the π key.
This procedure can be used to quickly generate a fraction that approximates any decimal number to any
desired number of digits. Let’s do it by hand one more time to appreciate the method.
1 16 1 39 1 257
1+ =1 + = 1+ =
1 9 1 22 1 145
1+ 1+ 1+
1 1 1
3+ 3+ 3+
2 1 1
2+ 2+
1 1
1+ 1+
1 1
1+
6
The fractions generated by this method are called the “convergents” to the number being approximated.
The infinite set of convergents for π begins thus:
Since the continued fraction algorithm is mathematically well-defined, it can be programmed on any
programmable HP calculator. Appendix A of this article is an HP-15C program that implements it.
Appendix B is an HP 35s program. The HP-15C program only needs 45 steps to convert any decimal
number into a fraction with any desired accuracy. Just input the decimal number, set the display FIX to
control the accuracy, and GSB A. The numerator appears. Press X<>Y to see the denominator.
Optionally press ÷ RCL – 0 to see the difference between the input decimal and the output fraction.
5 +1
The input that makes the program work the hardest is the golden ratio, , because its partial
2
quotients are { 1, 1, 1, 1, 1, 1, 1, … }. The HP-15C program says that the fraction that best approximates
6765
the golden ratio to 7 digits is . What are its other convergents? It has more than any other number!
4181
If you write the first several down, you will see the Fibonacci Series hidden in it. Many other Number
Theory goodies await discovery along the path of the Continued Fraction Algorithm. Happy exploring!
Instructions:
(1) Use a FIX setting as desired.
(2) Place decimal number in X.
(3) GSB A ➯ see numerator.
(4) X<>Y ➯ see denominator.
Optional: ÷ RCL - 0 ➯ “error” (difference between input and output)
Example:
Convert π to a fraction accurate to 5 decimal places.
(1) FIX 5
(2) π
(3) GSB A ➯ 355
(4) X<>Y ➯ 113
Answer: 355/113
Notes:
(a) The program uses registers 0 through 2, so make sure that enough register memory is allocated.
(b) As always, turning off USER mode makes it easier to key in the program.
(c) When the program finishes, register 0 contains your original input. Registers 1 and 2 are used for
temporary storage.
(d) The Golden ratio, 1.618033989, is the most challenging (takes longest time).
Note: Program execution is the same as for the HP-15C program except it is launched by XEQ A ENTER
instead of GSB A.
No attempt was made to optimize the program by using line number addressing instead of labels.
Introduction
The famous French writer Voltaire seems to have captured the driving force behind progress when he said
“The better is the enemy of the good.” Recently, Hewlett-Packard launched the HP 15c Limited Edition
(which I will nickname the HP-15C+ in this article) due to popular demand by users of vintage HP
calculators. The vintage HP-15C was the second calculator to support the Solve and Integrate features,
after the HP-34C, and the first one to support both complex math and matrix operations. The HP-15C has
earned many fans that went online and posted a petition for HP to re-launch that machine. After many
years of waiting, I am happy to report that HP has listened and indeed responded to these very loyal HP-
15C fans. The HP 15c Limited Edition is here and comes with the same features of the original HP-15C,
with added speed. I mean this little puppy is fast! How fast, you may ask? This article answers this
question by comparing the speeds of the old HP-15C and the new HP 15c.
LBL B
DSE 0
GTO B
RTN
Label A initializes the loop counter stored in register 0. The initial loop counter value depends on the test
and the machine version. Label B performs the test to time an empty loop. I use the DSE command
because I can count down from numbers greater than 1000. If I use the ISG command instead, I am
limited to a maximum of 1000 iterations. On the HP15C+ I use tests programs that loop10000 times.
This high loop iteration justifies selecting the DSE command to control the loop iterations.
Here is a sample code that calculates the sine of a 0.6 radians, stored in register 9:
LBL A
RAD # switch to radians
.
6
STO 9 # store 0.6 in memory register 9
EEX
2
STO 0 # initialize the loop counter for the HP-15C
RTN
LBL B
RCL 9 # recall argument value
SIN # calculate the sine
DSE 0
GTO B # end of the test loop
RTN
Testing the Σ+ operation is a bit different. The test repeats the process of adding 100 random values for Y
and X in the statistical registers. The test uses the following code:
LBL A
EEX
2
STO
0 # initialize the loop counter for the HP-15C
RTN
LBL B
CLEAR Σ # clear the statistical summations
EEX
2
STO 1 # set counter for the inner loop
LBL 0 # start inner loop
RAN # # generate random Y
RAN # # generate random X
Σ+ # add to statistical summations
DSE 1 # test end of inner loop
GTO 1 # end of inner loop
DSE 0 # test end of outer loop
GTO B # end of outer loop
RTN
• Calculating e^π using a Taylor expansion polynomial. The computations add the values of
different Taylor expansion terms to obtain the exponential value. When the program rounds the
value of a term to zero (based on the current number of displayed digits), the iteration stops. I
tested calculating e^ π in FIX 7 and FIX 4 display modes.
• Calculating the integral of 1/x from 1 to 100.
• Calculating the root of the function e^x-3*x^2 between 3 and 4.
• Converting decimals to fractions. This test is based on a program written by Joseph Horn . The
test uses the value of the golden ratio, since it requires more calculations than many other common
numbers with fractions, like π or e. The test program uses the FIX 9 display mode.
The code in this section is aimed at the new HP-15C+. The code for the vintage HP-15C differs only in
the values stored in the loop counter (in register 0), which is set to 10.
The code for calculating e^π using Taylor approximation is:
LBL A
FIX 7 # can set to other values
Π # π
STO 1
EEX
3
STO 0 # initialize the loop counter for the HP-15C+
RTN
LBL B
0
STO 2 # counter
STO 3 # sum
1
STO 4 # x^counter
LBL 1 # start of inner loop
RCL 4
RCL 2
X!
/
STO+ 3
RND # round term
LBL B
RCL 4
GSB A # call Joe Horn’s program in LBL A
DSE 3
GTO B
RTN
LBL C initializes the test. Enter the value for the number with a fractional part and then press the key C.
In the case of the golden ratio, enter 1.6180339887 and then press the key C. Label B starts the test by
repeatedly recalling the value with a fractional part from register 4 and then invoking Joe’s code in label
A as a subroutine. The test code uses register 3 as a loop counter since Joe’s code already uses registers
0, 1, and 2.
Results
Table 1 shows the speed comparison summary. The table has the following columns:
1. The Test column specifies the test conducted.
2. The HP-15C column lists the speed for the various operation sets for the HP-15C.
3. The HP-15C+ column lists the speed for the various operation sets for the HP-15C+.
4. The Speed Ratio column displays the ratio for the speeds of the HP-15C+ divided by the speeds
of the HP-15C.
Tables 2 and 3 show the detailed timings for the HP-15C and HP-15C+, respectively.
Looking at the above results you can see that speed ratios of 150 and higher appear in register operations,
rolling down the stack, addition, subtraction, the INT function, the FRAC function, the program that
calculates e^ π (using the Taylor approximation), and the decimal to fraction program. The two programs
in the list show high speed ratios, perhaps, because they do not use functions that exhibit lower speed
ratios (like the trigonometric, hyperbolic, and logarithmic functions.) This conclusion agrees with the
relatively lower speed ratio of the statistical summation test. The listing for that test shows that it uses
two RAN# commands (which show speed ratios of about 113). The speed ratio for the statistical
summation test is very close to the speed ratio of the RAN# command. This observation tells me that the
RAN# is the slowest part of the statistical summation test—akin to a bottleneck, if you like.
The HP-15C Limited Edition shows speed ratios, compared to the vintage
HP-15C, in the range of 102 to 189. Even at the lower end of this range, the
new HP-15C offers very impressive speed!
Introduction
The dynamic range of many scientific calculators (entry level) is specified as -9.999999999 x 10-99 for the
smallest number to 9.999999999 x 1099 for the largest number if it is a ten digit machine. This was the
range of the first scientific calculator the 1972 HP-35A. How large or small is the number represented by
10±99?
As a number 1099 is 1 followed by 99 zeros. 1099 = 1,000,000,000,000,000,000,000,000,000,000,000,000,
000,000,000,000,000,000,000,000,000,000,000,000,000,000,000,000,000,000,000,000,000.
Add another zero and the number has a special name, googol (not Googel). Ten raised to the power of
googol is called a googolplex. This number is impossible to write out.
In the early 1970’s scientific calculators were often called electronic slide rules and this 10±99 range was
simply astounding. It was so unbelievable that the literature of the time often explained that all known
values in the universe would easily be covered by numbers within the 10±99range. Certainly all known
physical constants were not known beyond six or seven significant figures so the HP-35A was what
would be called today “The Killer Ap” of its day. Let’s explore a few real world “numbers” just to see
how meaningful 10±99 really is.
An answer might be 2.75, 0.000275, or 2.75 x 1013. Two problems had to be solved, one for the
significant digits and one the power of ten of the answer. The power of ten part was done in your head.
The student quickly learned to express all problems with the numbers in scientific notation. Let’s use the
example of the problem of multiplying 321 x 57,943.
Mentally the problem is expressed as 3 x 102 multiplied by 6 x 104. 3 x 6 = 18. You add exponents if
they are multiplied so the answer is 18 x 106 or in normal form 1.8 x 107. The correct answer from my
calculator is 18,599,703 or 1.8 x 107 as an approximation. Teachers often lament that students all too
often just mindlessly accept the answer they see in the calculator display without their having any real
world sense for the number. This is the reason you should always “sketch out the problem” just as a
technician or engineer must always sketch out a circuit he or she is wiring no matter how simple it may be
– even a battery, switch, resistor, and LED series circuit. To do otherwise is just plain sloppy work (prone
to human error) no matter how brilliant you may be.
A light year is the distance light travels in one year. Since we need to calculate the answer in cubic inches
The distance light travels per year in mm = 3.72461748223 x 1017 in/lt yr. x 25.4 mm/in =
9.4605284052 x 1018 mm/lt yr.
The Universe volume in cubic millimeters = 4*π(9.4605284052 x 1018)3/3 = 3.54678441413 x 1057 mm3.
6.0221367 x 1023 * 5.55555555556 x 10-2 = 3.3456315 x 1022 molecules per cc. The number of atoms is
3 times larger because each molecule has three atoms, so there are 1.00368945 x 1023 atoms of water per
cubic centimeter, cc, of water.
The answer from example 1 calculated the universe as 3.54678441413 x 1057 mm3. This value converted
to cm3 is divided by 1,000 mm3/cm3 which = 3.54678441413 x 1054 cm3 in the Universe.
The number of water molecules to fill the Universe is 3.54678441413 x 1054 cm3 x 1.00368945 x 1023
molecules /cm3.= 3.55987009789 x 1077 H2O molecules.
Only 1077? That number is still quite small compared to 1099.
14 billion is 1.4 x 1010 yrs x 3.1556,9259747 x 107sec./yr. = 4.41796963646 x 1017 = Universe age in
seconds.
If the big bang theory theorizes the early events in one nanosecond increments how old is the Universe in
nanoseconds (light travels about one foot in one nanosecond)? There are 109 nanoseconds in one second.
The age of the Universe in nanoseconds is approximately 4.41796963646 x 1026 nanoseconds.
Only 1026? That number is almost tiny compared to 1099.
Can the reader think of any physical reality object that requires its value to be expressed as a numeric
value that cannot be calculated by an HP scientific or graphing calculator?
Obversations and Conclusions
How Large is 1099?
Our current generation of calculators have an even wider dynamic range with the upper limit of 10±499.
As mentioned with our piano key example it is easy to generate very large or very small numbers
mathematically. The new question is, “How Large is 10499?”
It is certainly larger than any reality numbers you may think about. Any factorial calculation greater than
253 (253! = 5.17346099264 x 10499) will cause overflow, but when it comes to “pure” mathematics a
finite state machine like a scientific calculator is essentially useless. The important question to ask is:
What are realistic dynamic range limits for a scientific or graphing calculator?
Epilog
Modern graphing calculators have morphed into symbolic math machines. The student asks: Is the hand
held calculator adequate?
“In the second presentation I asked if Einstein could have done some General Relativity work on an HP50g, even if
this were to be very difficult. Say he was shipwrecked on a deserted island with only an HP50g for company, could
he have used it? Fundamentally, the Computer Algebra System on an HP50g is just about powerful enough to do
some symbolic tensor calculus, but I doubt that Einstein or anyone else would happily do so. Instead, I suggested
that if any General Relativity were to be done on an HP50g, the programming would have to be done first in some
other language than RPL. If Einstein were to use an HP50g he could do so, but he would have to learn to program
in an appropriate language, as well as learn to use Tensor Calculus, which took him some effort. He would need
some other computer to write the programs first, then translate them into code to run on the HP50g, and transfer it
to the HP50g. So he could just about have used an HP50g on his desert island, but with considerable difficulty.”
“This means that if E wants to work with tensors on an HP calculator, then he or she has to get someone to write
the programs for them. This was discussed on comp.sys.hp48 two years ago, but at that time no-one had done it
yet. I just wonder if anyone will ever have the time or the wish to do so. And that is the answer to “Should E use
one?” They won’t know until they try! So, they should find someone willing to write the necessary programs for an
HP50g! I’m not volunteering unless I get a substantial university grant to do it. If I do, I’ll let you know next year.
Otherwise I might go back to programming on my HP-41 and give a talk on that.
____________________________________________________________________________________
Notes
(1) See HP Solve issue 21 for an article, Calculating before Calculators, that describes the various pre-calculator
methods of performing numerical calculations.
http://h20331.www2.hp.com/hpsub/downloads/Newsletters_HP_Calculator_eNL_12_Dec_2010_v1.pdf
(2) This number is from: http://en.wikipedia.org/wiki/Observable_universe .
(3) See http://www-users.cs.york.ac.uk/susan/cyc/g/graham.htm for an explanation of Graham’s number and other
large numbers and how they are named and expressed.
(4) All of the numbers given in this article are silly in terms of the number of digits provided. They are given as
shown in your 12 digit calculator display as a means of verifying that you have pressed the correct keys and
you have the correct answer.
HP News
HP hosted an HP
August 30, 2011, NYC @ Harry’s Cafe & Steak.
th
12c 30 Anniversary Celebration for the financial
press. This event also served to announce the re-
release of the popular HP-15C now known as the
HP 15c Limited Edition.
There was even a business card drawing for one of each of the two new products. Radiris Diaz of
CuteGeek.com won the HP 12c 30th Anniversary Edition and Stacy Cowley of CNNMoney.com won the
HP 15c Limited Edition.
The press event was planned in advance of the storm, but everything dried out just in time for everyone to
have a great evening sharing stories of how the finance world has been impacted by the over 30 years of
HP-12C usage.
Community News
As I mentioned in the last issue, the HP user Community is participating in the “designing” of two new
calculators. If you think that it is an easy task to design a calculator you soon realize that the devil is in
the details. One of the most challenging aspects is accuracy - a hallmark of HP calculators. Another
aspect is the documentation.
Suppose you are King and you make all the decisions. You have the factory standing by to start the new
model on a specific date. If you don’t release the new design the factory will (instead) make another
model (HP-12C’s are always needed ) and the new model production will have to be moved to another
date – many months later. Now you have to make an important change at the last minute. What do you
do? Do you delay the machine and miss an important sales period such as Christmas or back to school?
In some situations the change could ripple through several teams (to change the hardware and/or the
documentation, etc.) to extend the introduction a full year. The market radically changes in one year.
Another topic being presented at HHC 2011 is a new HP calculator restoration book by Geoff Quickfall.
HP Solve readers will remember that Geoff was interviewed in HP Solve issue 12, and the book was
S03 – The HP-12C, 30 Years and Counting by Gene Wright and me. This is a short review of the very
popular business calculator with new information related to the manufacturing history and the MSRP
pricing over the 32 years of its being manufacturered.
S05 – Converting Decimal Numbers to Fractions by Jospeh Horn. Joseph explains the continued
Fraction Algorithm and provides an RPN HP-15C program that quickly converts any decimal to a fraction
to the accuracy set by the FIX display value. The program listing provides program line numbers for
easier conversion to other RPN calculators.
S06 – Limited Edition HP-15C Execution Times by Namir Shammas. Namir examined the new
incarnation of the 15C LE and compared the execution times for many of the instructions, programs, and
built-in applications with those of the original 15C. The results were surprising because the expectation
was that it would be faster by some factor. This turned out to be very untrue. The speed advantage
actually ranged from 102 to 180 times faster. A total of 44 speed ratios were measured, the total is
5,659.06 and the average is 128.615. Saying that that the HP-15C LI is about 128 times faster is a good
average value to use.
S09 – #8 in the Math Review Series: Mathematical Constants – π This constant is probably one of
the most well-known of all mathematical constants. Here is an overview of this interesting number
with lots of links for further astonishment and study.
That is it for this issue. I hope you enjoy it. If not, tell me!
Also tell me what you liked, and what you would like to read about.
X < > Y,
Richard
Email me at: hpsolve@hp.com
16 commands, 88.0 Bytes, # 670Bh. Timing: 1987⇒ 2.131987, 3.131987, 11.131987 in 0.148_sec.
Observations:
1. Occurrences per century. Depending on how you wish to define a century - I called one
traditional and the other popular - the total number of occurrences of F13 per century vary. The
popular belief of the recent new millennium starting on January 1st 2,000 picks up the cycle “out of
sync” so that there are not 172 Friday The 13ths for each century.
2. The calendar repeats every 4 centuries. Roger Hill participated in a contest to print calendars and
he investigated them. I asked him for his thoughts. His comment is: "One thing to note is that, for
the presently-used Gregorian calendar, you only have to look at a 400-year interval, say 1600-1999,
because the calendar repeats every 400 years." The total number of Friday The 13ths per cycle is
688.
3. Having three Friday The 13ths in one year occurs 15% of the time. Roger observes: "I did some
research on Friday the 13ths (or F13's for short) a long time ago and (I think) figured out the
probability of having one, two or three in a year chosen at random. I don't remember what I got
now, but with a little luck I might be able to find it."
"Another thing is that, in any year, the 7 months from May through November (inclusive) always
start on 7 different days, so there will always be one and only one F13 in that 7-month interval. For
non-leap years, the ones that have three F13's are the ones in which February 13 and March 13 are
on a Friday, which happens in about 1 year out of 7. The criterion is different for a leap year, but
again it happens in about 1 leap year out of 7. Hence the 14 or 15 per century that you found." I did
not pick up that bit of information from the tables.
4. Having one or two Friday The 13ths in one year is about equal probability. This is obvious
from the tables. Roger notes: "An interesting fact about the Gregorian calendar is that, while you
would expect a year chosen at random to be equally likely to start on any day of the week, over the
Function accuracy, however, is only one part of the accuracy “specification.” If the calculator has built in
applications such as complex numbers, solve, integrate, or matrices the accuracy specification is even
less. In some cases surprisingly, perhaps three digits, less. The average user will seldom consult the
User’s Guide for an accuracy statement – if there is one – to know how reliable, how accurate, all those
digits are. After all, the calculator designers should know what they are doing.
Customers needing to buy a calculator won’t spend very much time delving into the details of calculator
accuracy other than to hear about the reputation of a particular manufacturer.
The calculator programmers that microcode the internal math routines spend their time trying to reduce
the number of instructions and to optimize them to run as fast as possible. They cannot be expected to be
numerical analysts in terms of mathematical accuracy. They just implement the algorithms they are
given.
In order to determine the accuracy of the number in the display, with its limited fixed number of digits,
you must know the correct answer to several digits greater than that produced by the function. Casual
calculator users just want to know how the various models compare, and they have devised various “tests”
to compare machines.
Comparison testing
One such scientific model test has been proposed by Mike Sebastian. He has a website dedicated to what
he calls Calculator Forensics(1). Below is the “forensics” data for HP calculators. The results in Tables 1
and 2 are for his single forensics trigonometry function test. Some of the data is current data not found in
Mikes 2007 table(2). The tables also add a Rating column and an Error column for comparison purposes.
The models are arranged in increasing accuracy order.
Here is the test problem. In degrees mode calculate arcsin (arccos (arctan (tan (cos (sin (9) ) ) ) ) )
algebraically, and 9 SIN COS TAN ARCTAN ARCCOS ARCSIN in RPN. A rating is suggested for
comparison purposes with the higher the value the better. The idea is that the function-inverse is
supposed to return the input, in this case 9 degrees. Supposedly the more accurate calculator (overall?) is
the one that is closest to 9. All calculators in the tables use BCD arithmetic except as noted + in table 1.
* Simply the number of zeros or nines following the decimal point. The higher the better. The rating
is one less than the exponent of the error.
+ Uses 80 bit microprocessor binary arithmetic which provides at least 23 decimal digits.
* Simply the number of zeros or nines following the decimal point. The higher the better.
The rating is one less than the exponent of the error.
Previous Article
U UH
Return to Top
#8 in Fundamentals of Applied Math Series
Mathematical Constants – π
Richard J. Nelson
Then there is the π joke. (Forwarded to me by Roger Hill Hill, PPC ROM Manual contributor, HP
calculator author and enthusiast, meet Roger at HHC 2011).
Question: How can you tell if a person is a mathematician, a physicist, an engineer, or a salesman?
Answer: Ask him what π is.
The mathematician will say, "It is the ratio of the circumference of a circle to the diameter".
The physicist will say, "It's 3.14".
The engineer will say, "It's about 3".
The salesman will say, "It's normally 3, but I can get it to you for 2.6".
Discovering π
π has a very long colorful and interesting history. Table 1 is a brief portion of the Chronological Table
taken from A History of π by Petr Beckmann, see page 196 shown as reference 1 of Appendix B.
Table 1 – The Early History of π
Date π Status
Circa 2000 BC Babylonians use π = 3 1/8.
Circa 2000 BC Egyptians use π = (16/9)2 = 3.1605.
th
12 century BC Chinese use π = 3.
Circa 500 BC I Kings vii, 23 implies π = 3.
Circa 440 BC Hippocrates of Chios squares the lune.
Circa 434 BC Anaxagoras attempts to square the circle.
Circa 430 BC Antiphon enunciates the principle of exhaustion.
Circa 420 BC Hippias discovers the quadratix.
Circa 335 BC Dinostratos uses the quadratix to square the circle.
The relationship of the circle circumference to its diameter has been known for about 4,000 years, but the
first use of the π symbol was by William Jones in 1706 because it was the first letter of the Greek word
for "perimeter" or "periphery.” The great mathematician Leonhard Euler used it in 1737 and it has been
accepted since. π is a lower case greek letter pronunced “pie” that is not capitalized if it begins a
sentence because the upper case greek letter, Π, is used for designating the product of a sequence. If C is
the circumference of the circle and D is the diameter of the circle, π is defined as shown below.
Measuring π
You may determine a desired number of digits of π in one of three ways. You may measure it, you may
calculate it, or you may look it up on the internet. I really wanted to make an attempt to measure π, but
the daily summer 107 degrees here in the Sonoran desert discouraged me. Instead, let’s perform a thought
measurement. Here are the materials I was going to use.
1. Three metal stakes (12” nails).
2. 210 feet of white polypropylene cord.
3. 110 feet of bailing wire.
4. 100 foot measuring tape.
5. Bicycle (front wheel and fork). I would have to borrow one.
6. Masking tape.
7. Carpenters marking pencil and ball point pen.
8. Large flat level smooth layout area at least 220 x 110 feet. (half of a 200 foot diameter circle.)
Here is the measurement procedure. The goal is to make all measurements to a resolution of 1/8 inches.
A. Remove the front fork from the bicycle.
B. Wrap a piece of masking tape around the tire and wheel and mark a line down the outside
perpendicular to the ground.
C. Attach one end of the wire to the frame near the wheel axle.
D. Move all of your materials to the layout area and pound in two stakes about 220 feet apart and
stretch the white polypropylene cord to serve as the diameter of the circle.
E. Pound in the third stake in the center of the white cord to mark the center of the circle.
F. Position the wheel 100 feet from the center (an assistant is most helpful here) and make a loop to
be placed over the center stake.
G. Holding the wheel with the wire tight (it should not stretch) rotate and align the mark on the rim
with the white cord.
H. Keeping the wire tight walk the wire around the semicircle keeping the wheel perpendicular to the
surface counting the number of times the tape comes in contact with the surface. When you are
You must account for the diameter of the center stake in your measurements.
Let’s assume the radius (wire length) is 100 feet or 1200 inches. The diameter is 2400 inches. After
doubling the thought wheel measured total length (measured circumference x revolutions + extra partial
revolution) the circumference is 7,539 3/4hs (7,539.75 inches). I realize that making these measurements
to a resolution of 1/8th inches takes great skill and is very difficult. Making the division we have π =
3.1415625.
The theoretical circumference with a 100.00000 feet diameter is 7,539.82236862 inches. Let’s assume an
error of ¼ inch too high for the circumference and ¼ inch too low for the diameter for a maximum error
on the high side. The calculated theoretical value for π is 7,540.07236862 divided by 2,400.25 =
3.14136959426. The difference (high error) with 12 digits of π is 0.00022305933. This means that we
can easily depend on three accurate decimal digits or 3.141 using our measurement method. If our work
is really practiced and careful we might hope for 3.1415.
Calculating π
Most students don’t even think of the numerical value of π because it is a keystroke
constant of most HP scientific and graphing calcualtors. This most popular and very
easily entered irrational and transcendental number is found on the HP 35s shifted COS
key as shown in Fig. 2.
All scientific and graphing calcualtors have π internally stored to a value much greater
Fig. 2 - HP35s π than what the display will show. This is required to perform the trig functions to their
specified accuracy. The number of digits of π stored in the HP-48/49/50 calculators, as a constant, is 31.
Here is a technique suggested by Joseph K. Horn to obtain the digits of π that are used by your HP
calculator for the sin function. Using RAD mode type in as many digits of π as you wish up to the limit
of the display. Do not round, use the truncated value. Suppose I only remember that π is 3.14159 and I
have an HP-15C handy. After typing in the 6 digits I remember (let’s pretend that the π key is broken) I
press the SIN key. The display will show: 0.000002654. Use FIX 9. Now multiply the number by
10000 to see 0.026535900. That is all the useful digits we can get, but we now have 12 accurate rounded
digits of π.
For 10 digit calculators the guard digits are usually two. For 12 digit calculators there are typically three
guard digits.
Now let’s use the “SIN π Digits” technique using the HP 35s. Keying in the maximum number (12) of
truncated digits, 3.14159265358, and taking the SIN we see 0.000000000001. Multiplying this by 1012
we will see 9.79323846264. We now have 24 accurate digits of π if we append these digits to our 12
input digits. This is shown in Fig. 1. No matter what you do you can’t get more digits. Checking the
normal π key for guard digits (the only other technique we have to squeeze out extra digits of a
calculation) all we get are 12 digits, 3.14159265359 which is the correctly rounded value of π to 12 digits.
Exactly what this means is unclear. We know that the Tangent function is the most sensitive to errors at
the extremes and that π is used internally by HP for the Tangent function to 31 digits in the later (1990)
RPL calculators. Below is a table of the current machines with a few historical models for reference.
π has been calculated to more than 5 trillion digits and memorized to 67,890 digits(5) by Lu Chao a 24-
year old chinese graduate student as recognized by the Guinness book of records. It took 24 hours and 4
minutes to speak them – without error! See Exploring π below to learn how this is done.
Calculating the decimal digits of π is a popular HP calculator programming challenge and various authors
have written programs to calculate the first 1,000 digits of π for the calculators of the 70’s and 80’s.
Often the various machines were compared for speed by using the number of hours it would take to
produce the first 1,000 π digits. This number was used because of memory restrictions. RPL calculators
of the 90’s and later have math libraries that calculate hugh numbers of the digits of π.
Calculating an approximation for π may be also be done by converting the decimal value desired to a
fraction using Joseph’s program found elsewhere in this issue. Using his program running on the HP 35s
results in the values shown in Table 3. The input is 3.14159265359. The accuracy is set by the FIX mode
which is one less than the Digits column.
If you use the program you will have two choices for the 12 digits. You may use a truncted input with the
twelfth digit 8 or you may use a rounded input with the twelfth digit 9. Joseph says the the rounded value
must be used because that is the mathematical basis for the continued fractions algorithm. If you use
3.14159265358 as the input several fractions will change and the error will be off by one digit.
Mathematicians have been calculating (approximating) π by ratios since the beginning and Beckmann
provides a table of these values as shown in Table 4.
12D Divided
N D Value Error
3 1 3.00000000000 -0.14159265359
22 7 3.14285714286 0.00126448927
333 106 3.14150943396 -0.00008321963
355 113 3.14159292035 0.00000026676
103,993 33,102 3.14158265301 -0.00000000058
104,348 33,215 3.14159265392 0.00000000033
208,341 66,317 3.14159265347 0.00000000012
312,689 99,532 3.14159265362 0.00000000003
833,719 26,5381 3.14159265358 0.00000000001
1,146,408 364,913 3.14159265359 0
4,272,943 1,360,120 3.14159265359 0
1,725,033 5,419,351 3.14159265359 0
Exploring π
The Internet provides a wealth of information about π. For a list of formulas that use π see note (6). You
may download millions of digits if you are interested and a text editor (Word) may be used to explore the
digits as you wish. Appendix A shows 10,000 digits. If you copy the desired digits into Word you may
use the word counting, find, and the find and replace functions to discover a few of the properties of this
irrational number. Highlight the digits before the ones of interest and use the word count function to
determine the location in the π string.
One question to ask is the distribution of the digits. If the digits 0 through 9 were evenly distributed there
would be 1,000 of each in 10,000 digits. I searched the first 10,000 digits for each digit and replaced it
with the same digit. Word returns the number of replacements made which is the count of each digit.
Table 5 shows the results. This is a trivial observation and Bailey(2) (1988) found that the first 30 million
digits are very uniformly distributed. I wonder how many hours it would take to use Word to make a
similar table for 30 million digits.
The highest digit count = 1,046, the lowest = 968, and the mean (of the occurances) is 1000. The standard
deviation, σx, (population) = 30.3743312684.
Are the digits of π random(7)? Speculation is that they are. Should the digits of π be used as an encoding
key? Probably not because they are so “common.”
How may you memorize long strings that contain the digits of π? Usually this is done using a “poem” or
specifically, a piem. Here is one devised by Sir James Jeans.
How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics.
The number of letters in each word is a digit of π. A more refined and very long (thousands of words)
piem starts:
Poe3, E1.
Near4 a1 Raven5
Midnights9 so2 dreary6, tired5 and3 weary5,
Silently8 pondering9 volumes7 extolling9 all3 by2-now3 obsolete8 lore4.
During6 my2 rather6 long4 nap3 – the3 weirdest8 tap3!
An2 ominous7 vibrating9 sound5 disturbing0 my2 chamber's8 antedoor8.
"This4", I1 whispered9 quietly7, "I1 ignore6".
This verse provides the first 42 digits of π. - Mike Keith, First verse of Near a Raven (with respects to
Edgar Allan Poe).
Randomly π
The French scientist Georges Buffon (1707 – 1788) conducted what he called his needle experiment to
calculate π which caused quite a stir in his day. Here is how it is described. If we have a uniform grid of
parallel lines, unit distance apart and if we drop a needle of length k< 1 on the grid, the probability that
the needle falls across a line is 2k/π.
The temptation to repeat this experiment is strong and in 1901 Lazzerine is reported to have found that π
= 355/113 = 3.1415929 after 34,080 tosses. This correct six decimal digit value had been previously
determined by Chinese mathematician Zu Chongzhi (429 – 501). It is obvious that “working backwards”
you could stop the number of tosses to achieve the desired accuracy and the 34,080 number seems to
indicate this.
The connection of randomness to π is strong. Note (8) provides a more modern approach to the fun of
using random numbers to predict π. These methods, however, are very very inefficient.
Memorizing the value of π is a popular pastime for young students and one of the easiest ways to
remember π is to think of the first three odd numbers, 1, 3 & 5. Write them down twice, 113355 and
divide the last three by the first three to get 3.141592^9205 on a 12 digit calculator. You used six digits to
get the answer and your result is correct to a truncated six decimal digits.
The review provides several resources and links to websites that are especially useful. Note (1) shows a
list of mathematical constants that will or have been discussed in this math review series. The next one
just may be the most interesting number you could put into a calculator display.
An appendix (A) includes a single page of the first 10,000 digits of π with a few of my observations of
“interesting” sequences including the Feynman point, note (3), which occurs at decimal digit 762. Word
(or other word processor) may be used to search through and discover various series of digits from a down
loaded file of digits from the Internet. One website, note (4), is especially designed to search 200 million
digits of π very very quickly. Reading through the techniques and methods used for this project provides
an interesting perspective on the fun aspect of π,
____________________________________________________________________________________
Notes for Mathematical Constants – π
(1). A selected short list of the first ten of the more common 47 mathematical constants from Wikipedia
is shown below. # 1 – 3, & 5 – 8 have been discussed in previous math review articles. #10 will be
discussed in Issue 26.
First # of Known
# Symbol Value Name Field N
Described Digits
1 c. 7th–5th
0 =0 Zero Gen R N/A
century BC
2 1 =1 One, Unity Gen R N/A
3 Gen,
i = Imaginary unit A 16th century N/A
Ana
4 ≈ 3.14159 26535 89793
Pi, Archimedes' constant or Gen, 5,000,000,
π 23846 26433 83279 T by c. 2000 BC
Ludolph's number Ana 000,000[30]
50288
5 ≈ 2.71828 18284 59045 Napier's constant, or Euler's
Gen, 100,000,
e 23536 02874 71352 number, base of Natural T 1618
Ana 000,000
66249 logarithm
6 ≈ 1.41421 35623 73095
Pythagoras' constant, square 137,438,
√2 04880 16887 24209 Gen A by c. 800 BC
root of 2 953,444
69807
(2). See a much larger table of the distribution of π digits at: Weisstein, Eric W. "Pi Digits." From MathWorld--A
Wolfram Web Resource. http://mathworld.wolfram.com/PiDigits.html This source discusses many interesting
patterns that occur in the decimal digits of π. A list of references is included.
(3). From the (2) source. “The sequence 9999998 occurs at decimal 762 (which is sometimes called the Feynman
point; Wells 1986, p. 51). This is the largest value of any seven digits in the first million decimals.
(6). A very nice list of formulas involving π may be found at: http://en.wikipedia.org/wiki/List_of_formulae_involving_%CF%80
(7). Here is a 2005 research study announcement that addresses using the digits of π as a source of random
numbers. The ScienceDaily article also provides some insights as to the methods used to determine
randomness. "Our work showed no correlations or patterns in pi's number set – in short, pi is indeed a good
source of randomness," Fischbach said. "However, there were times when pi's performance was outdone by the
RNGs."
"These tests are simple to reproduce with a desktop computer. All you need is time," he said. "It took us almost a
year of work to crunch these numbers. We have included the program we used in the paper if anyone would like to
try doing the analysis with a larger number set. I hope someone will because pi shows up in security systems,
(8). Here is a link for the results of a more modern approach to using random numbers to calculate π. The method
uses Python for the calculations.
http://www.stealthcopter.com/blog/2009/09/python-calculating-pi-using-random-numbers/
Fig. B1 – 1971 hardbound Golem Press Edition. Fig. B2 – 1971 softbound Golem Press Edition.