ࡱ > n p m ` / bjbj . yb 1 1 1 1 V2 T 2 3 3 3 3 3 5 5 5 IS KS KS KS KS KS KS $ 8U h W \ oS Q 5 5 5 5 5 oS 3 3 S S S S 5 3 3 IS S 5 IS S S S 3 2 g|< 1 E . S IS S 0 T S W O W S W S , 5 5 S 5 5 5 5 5 oS oS S 5 5 5 T 5 5 5 5 D+ 0 0 Xnumbers.dll tutorial
General description
The current, industry-wide standard precision for personal computers is 32 bits, which is equivalent to approximately eight decimals. Double precision, standard in Excel and the default in VBA, has a precision of 64 bits, or about 15 decimals. For most problems this is completely satisfactory, but there are some mathematical operations, such as the multiplication or inversion of large matrices, that require better precision and/or special approaches to keep truncation errors from accumulating in the answer. For such problems, as may be encountered in least squares fitting to multi-term polynomials, or in multivariate applications, extended precision may be needed. It may also be useful as a diagnostic tool in program development and testing, to pinpoint precisely where numerical accuracy can be lost. And as scientists increasingly learn how to use Excel for non-trivial data analysis tasks, the need for extended precision will only increase.
Xnumbers.dll is an ActiveX dynamic link library, compatible with Windows 2000, NT, and XP, which allows Excel macros to operate in extended precision, of up to 200 decimals. It was created by Leonardo Volpi and his students, and is freely downloadable from http://digilander.libero.it/foxes. It handles numbers internally as character strings, divides them into packets of up to six decimals, and then manipulates these packets to perform elementary mathematical operations such as addition, subtraction, multiplication, and division. More complicated functions then use these elementary functions as building blocks. In the end, the strings are converted back to numbers. The method used is similar to that in the Fortran77 program MPFUN and in its recent successor, ARPREC, which allows arbitrary precision in Fortran90 and C++. For a description of MPFUN see D. H. Bailey, ACM Trans. Math. Software 19 (1993) 288-319; for ARPREC, see http://crd.lbl.gov/~dhbailey/mpdist/.
One of the beauties of Xnumbers.dll is that it allows the user to modify existing macros while keeping its user interface untouched. In other words, all high-precision operations are performed in the background, and the software user need not see a difference in either input or output. (High-precision results can of course be displayed in full when so desired, as in the code example of Test 2b shown below). Xnumbers is free, it is powerful, and it is fun, but it does slow down the computation considerably. Extended-precision macros are therefore not recommended for routine use, and they are best labeled as such, e.g., with the prefix x, as we have done with the macros in the xMacroBundle. .
Installation
You will need to have an unzip routine handy, because the download is a zip file, and an Adobe Acrobat reader, because the documentation comes as a pdf (portable document format) file. Both of these can be freely downloaded from the web, and should already be on your desktop as standard tools.
Here is how you can install the dynamic link library. (A specific location is selected here for your convenience, but you can of course place it somewhere else, in which case you should make the corresponding changes in what follows.) Go to the website http:// digilander.libero.it/foxes, select Downloads, and download the xnumberdll.zip to your desktop. Unzip the file, and store its three components, readme, xnumbers.dll, and Xnumbers_DLL_Ref, on the desktop. The readme file contains nothing important: read it and either discard or store somewhere. The reference file, Xnumbers_DLL_Ref, is very worthwhile and has many useful details. You may want to print it out, but please do that later, because it is a 115-page tome. For now we will focus on xnumber.dll.
On your desktop, click on My Computer ( Local Disk (C:) ( Program Files ( Microsoft Office ( Office (Office11 in Excel2003) ( Library. Then press File ( New ( Folder, and change the name of the so generated folder from New Folder to Xnumbers, then open that file (which, as a new file, should be empty). If necessary, reduce the space taken up by this folder and/or move it, so that you can find xnumbers.dll on the desktop. Click on xnumbers.dll, copy it into the open folder, and close the folder. Now, xnumbers.dll is stored safely, and presumably in good company, next to Solver.
In order to make sure that Excel can find it, open Excel, go to the Visual Basic Editor with Tools ( Macro ( Visual Basic Editor or, faster, with Alt(F11 (on the Mac: Opt(F11), and then, in the VBEditor menu bar, select Tools ( References. This will display the dialog box for the VBA Project. Click on Browse, select Files of type: Type Libraries (*.olb; *.tlb; *.dll), in the Add Reference dialog box go to My Computer ( Local Disk (C:) ( Program Files ( Microsoft Office ( Office ( Library ( Xnumbers, open the file, and type Xnumbers.dll in the Filre name window if it is empty. Click Open it, which will bring you back to the References VBA Project dialog box. Now use the Priority Up button to ease xnumbers up the list so that it is contiguous with the other tick-marked items. Click OK. From now on, Excel will know where to find xnumbers.dll.
Using Xnumbers in Excel VBA
Before the extended-precision instructions and functions (here abbreviated as Xfunctions, Xinstructions, etc.) can be used in any particular function, macro, or subroutine, the class of Xnumbers must first be invoked (in computerese: instantated) with the line
Dim MP As New Xnumbers
and, at the end of the code, deactivated with
Set MP = Nothing
Thereafter, individual Xfunctions can be called between those two statements with the prefix MP. followed by the specific Xinstruction, as in
Sub Test1a()
Dim a, b, c
Dim MP As New Xnumbers
a = "1234567890.0987654321"
b = MP.xInv(a)
c = MP.xExp(a, 40)
Debug.Print "a = ", a
Debug.Print "1 / a = ", b
Debug.Print "e ^ a = ", c
Set MP = Nothing
End Sub
which will display the values of a, b = 1/a, and c = ea in the Immediate Window. (You can open the Immediate Window in the VBEditor, with Ctrl(G.) For 1/a in the above example you will obtain answers to the default value of 30 decimal places. The numerical precision of each result can be specified individually in its instruction, as illustrated for c, which is computed and displayed to 40 decimals. In this example, the Immediate Window will show
a = 1234567890.0987654321
1 / a = 8.1000000730620006590111459442E-10
e ^ a = 1.60263297808435399238741617564E+536166022
If a had been defined without double quotes in the macro, i.e., as
a = 1234567890.0987654321
then the result for x in the Immediate Window would have read
a = 1234567890.09877
but those for 1/a and ea would not have been changed.
For more than occasional use of extended precision, it is often more convenient to use a collective statement, as in
Sub Test1b()
Dim a, b, c
Dim MP As New Xnumbers
With MP
a = "1234567890.0987654321"
b = .xInv(a)
c = .xExp(a, 40)
Debug.Print "a = ", a
Debug.Print "1 / a = ", b
Debug.Print "e ^ a = ", c
End With
Set MP = Nothing
End Sub
In each macro, subroutine, or function, the number of decimals used can also be set globally, as with
DgtMax = 150
Upon completion of the macro, its results can be written onto the spreadsheet, where they will usually be converted back to double precision, and in their display may be limited further by the chosen cell format. Therefore, all calculations that are critical to obtaining high computational accuracy should be performed within the macro, or inside the functions and subroutines it calls, before returning the final result to the spreadsheet. From the user point of view, this is the simplest approach, because in this way the final output doesnt even show that any special processing was used, although it may generate complaints about the apparent sluggishness of the calculation.
Below is a macro that takes a number from the spreadsheet, raises it to the tenth power, and returns the answer to the spreadsheet. To use it, place a number in a spreadsheet cell, or highlight any other cell containing a number, and call the macro. The answer will appear in the spreadsheet cell immediately below it, as well as in two message boxes.
Sub Test2a()
Dim MP As New Xnumbers
Dim a
Dim b As Double
CellValue = Selection.Value
a = MP.xPow(CellValue, 10, 40)
Selection.Offset(1, 0).Select
Selection.Value = a
MsgBox "x ^ 10 = " & a
b = a
Selection.Offset(1, 0).Select
Selection.Value = b
MsgBox "x ^ 10 = " & b
Selection.Offset(-2, 0).Select
Set MP = Nothing
End Sub
The message boxes display both a and b; for an input value of 12.34567, the message boxes will show
a = 82252032957.44822844133332922210978090071
and
b = 82252032957.4482
respectively, which only differ in that b is dimensioned As Double, which will override the extended precision. Note that the answers shown on the spreadsheet, like that for b in the message box, will have been reduced to double precision, because upon converting the string back to a number, it will be stored in 64 bits. But neither the message box nor the Immediate Window has that constraint, and both can therefore display the extended-precision result in full.
Still, you may want to find the extended-precision results on the spreadsheet rather than having to copy them from a message box or VBA debugging tool. You can prevent conversion into double precision by concatenating an apostrophe with the numerical value, which will make Excel treat the result as a text string:
Sub Test2b()
Dim MP As New Xnumbers
Dim a, c
Dim b As Double
CellValue = Selection.Value
a = MP.xPow(CellValue, 10, 40)
Selection.Offset(1, 0).Select
Selection.Value = a
MsgBox "a = " & a
b = a
Selection.Offset(1, 0).Select
Selection.Value = b
MsgBox "a ^ 10 = " & b
c = "'" & a
MsgBox "a ^ 10 = " & c
Selection.Offset(1, 0).Select
Selection.Value = c
Selection.Offset(-3, 0).Select
Set MP = Nothing
End Sub
The column with the receiving cell should be wide enough to accommodate the extended-precision result. Alternatively, you can left-align that receiving cell, and make sure that cells to its right are empty, so that the extended-precision result can overflow into those empty cells.
Now that you have the tools to write macros using Xnumbers, you will want to know the available Xinstructions. Some of the most useful ones are listed below. Apart from the usual set of elementary instructions, and an ingenious function evaluator, Xnumbers has an extensive set of instructions for manipulating vectors and matrices, and for computations involving complex numbers. In fact, so many of the latter functions are absent from VBA that Xnumbers also contains their standard, double-precision equivalents, distinguishable by the prefix s.
Function Instruction Operation
items within straight brackets are optional
the prefix MP is implied
Elementary single-argument instructions
Negation .xNeg(a) a
Absolute value .xAbs(a) (a(
Inversion .xInv(a) 1/a
Elementary two-argument operations
Addition .xAdd(a,b [,Dgt Max]) a+b
Subtraction .xSub(a,b [,Dgt Max]) ab
or: .xAdd(A1, .xNeg(A2)); do NOT use xAdd(A1, A2)
Multiplication .xMult(a,b [,Dgt Max]) ab
Division .xDiv(a,b [,Dgt Max]) a/b
or: .xMult(A1, .xInv(A2))
Integer division .xDivint(a,b [,Dgt Max]) int (a/b)
b ( 0
Integer modulus .xMod(a,b [,Dgt Max]) b mod (a/b)
Integer power .xPow(a,n [,Dgt_Max]) a n
if n is a real number, int(n) is used for the exponent
Arbitrary power .xExp_Base(a, b [,Dgt Max]) ab
a > 0
Square root .xSqr(a [,Dgt Max]) (a
a > 0
Integer root .xRoot(a,n [,Dgt Max]) a1/n
if n is non-integer, int(n) is substituted for n
xroot(100,9,45) ( 1.66810053720005875359979114908865584747919268
Multi-argument operations
V is a vector , i . e . , a n a r r a y o f n u m b e r s v i i n a c o n t i g u o u s r o w o r c o l u m n
S u m m a t i o n . x S u m ( V [ , D g t M a x ] ) S v i
P r o d u c t . x P r o d ( V [ , D g t M a x ] ) P v i
T r i g o n o m e t r i c f u n c t i o n s
a l l a r g u m e n t a r e i n r a d i a n s
S i n e . x S i n ( a [ , D g t M a x ] ) s i n a
C o s i n e . x C o s ( a [ , D g t M a x ] ) c o s a
T a n gent .xTan(a [,Dgt Max]) tan a
Inverse Sine .xAsin(a [,Dgt Max]) arcsin a
Inverse Cosine .xAcos(a [,Dgt Max]) arccos a
Inverse Tangent .xAtan(a [,Dgt Max]) arctan a
Exponential and logarithmic functions
Exponential .xExp(a [,Dgt Max]) ea
Natural Logarithm .xLn(a [,Dgt Max]) ln a
Decimal Logarithm .xLog(a [,Dgt Max]) log a
Hyperbolic Sine .xSinh(a [,Dgt Max]) sinh a
sinh a = ( ex ex)/2
Hyperbolic Cosine .xCosh(a [,Dgt Max]) cosh a
cosh a = ( ex + ex)/2
Hyperbolic Tangent .xTanh(a [,Dgt Max]) tanh a
tanh a = (ex ex)/( ex + ex)
Inverse Hyperbolic Sine .xAsinh(a [,Dgt Max]) arsinh a
arsinh a = ln[a+((a2 +1)]
Inverse Hyperbolic Cosine .xAcosh(a [,Dgt Max]) arcosh a
arcosh a = ln[a+((a2 1)], a > 1
Inverse Hyperbolic Tangent .xAtanh(a [,Dgt Max]) artanh a
artanh a = (1/2) ln[(1+a)/(1a)]
Statistical functions
Factorial .xFact(n [,Dgt Max]) n!
for DgtMax < 100, returns result in scientific notation
Binomial coefficient .xComb(n, m [,Dgt Max]) EMBED Equation.3
Some useful co n s t a n t s
T h e b r a c k e t s a r e r e q u i r e d , e v e n i f e m p t y .
D g t M a x c a n g o u p t o 4 1 5 f o r p, p/ 2 , p/ 4 , 2 p, e , l n ( 2 ) , l n ( 1 0 ) , a n d g
p . x P i ( [ D g t M a x ] ) p
a l s o : . x P i 2 ( ) f o r p/ 2 , . x P i 4 ( ) f o r p/ 4 , a n d . x 2 P i ( ) f o r 2 p
e . x E ( [ D g t M a x ] ) e
a l s o : . x L n 2 ( ) ( 0.693147180559945309417232121458
a n d : . x L n 1 0 ( ) ( 2.30258509299404568401799145467
g . x E u ( [ D g t M a x ] ) g
E u l e r s c o n s t a n t , 0 . 5 7 7 2 1 5 6 6 4 9 0 1 5 3 2 8 6 0 6 0 6 5 1 2 0 9 0 0 8 &
M o r e s i n g l e - a r g u m e n t i n s t r u c t i o n s
I n t e g e r p a r t . x I n t ( a )
c o m p u t e s t h e l a r g e s t i n t e g e r s m a l l e r t h a n a, e.g.,
xInt(2.14) = 2, xInt(2.99) = 2, xInt(2.14) = 3
Decimal part .xDec(a)
computes the decimal part of a, e.g.,
xDec(2.14) = 0.14, xDec(2.14) = 0.14
Round .xRound(a, n)
rounds a to n decimal places, e.g.,
xRound(2.14) = 0.14, xDec(2.14) = 0.14
Comparison .xComp(a [,b])
.xComp(a, b) = 1 for a > b, = 0 for a = b, = 1 for a < b
b = 0 when omitted, i.e., .xComp(a) returns the sign of a
Do NOT include logic operators inside the argument, as then the logic is applied before the extended precision. Therefore, replace If ab Then by If .xComp(a,b) > 0 Then.
Formating instructions
Format .xFormat(a [,Digit_Sep])
formats a in comma-separated groups of Digit_Sep
default: Digit_Sep = 6. For a = 1234567.89012345,
.xFormat(a) = 1,234567.890123,45 and
.xFormat(a,3) = 1,234,567.890,123,45
Unformat .xUnformat(a)
Removes formatting commas from a
Matrix and vector operations
We distinguish here between vectors V of dimension m, a positive integer, and matrices M (or A, B, C, etc.) of dimension m(n, where n is another positive integer. A matrix is called square when m = n. Individual matrix elements are denoted by aij.
Absolute value .xMatMod(M [,DgtMax]) ((M((
Absolute value or modulus = EMBED Equation.3
Determinant .xMatDet(M [,DgtMax]) (M(
M must be square. Returns ? when M is singular.
Addition .xMatAdd(M1, M2 [,DgtMax]) M1 + M2
M1 and M2 must have the same dimension m(n
Subtraction .xMat Sub(M1, M2 [,DgtMax]) M1 M2
M1 and M2 must have the same dimension m(n
Scalar product .xProdScal(V1, V2 [,DgtMax]) V1 ( V2
V1 and V2 must have the same dimension m
Vector product .xProdVect(V1, V2 [,DgtMax]) V1 ( V2
V1 and V2 must have the same dimension m
Multiplication .xMat Mult(M1, M2 [,DgtMax]) M1 M2
when M1 is m(p , M2 must be p(n
Integral power .xMat Pow(M, n [,DgtMax]) Mn
M must be square, m(m, and n a positive integer
Inversion .xMat Inv(M [,DgtMax]) M1
M must be square, m(m
Special matrix functions
LU decomposition .xMat _LU(M [,Pivot] [,DgtMax])
M must be square, m(m. Returns the Lower and Upper triangular matrices that satisfy M = L U or, when Pivot is true, M = P L U where P is the permutation matrix.
LL decomposition .xMat _LL(M [,DgtMax])
Cholesky decomposition. M must be square, m(m.
Linear equation solver .xSysLin(A, B [,DgtMax])
Solves A X = B to yield X = A1AX = A1B
A must be square, m(m, and B must be m(1 or m(n
Similarity transform .xMat BAB(A, B [,DgtMax])
Computes C = B1 A B. A and B must be square, m(m.
Gauss-Jordan .xGaussJordan(M, n, m, Det, Algo, DgtMax)
performs Gauss-Jordan matrix reduction.
Special least squares functions
Least squares fit .xRegLin _Coeff(y, x [,DgtMax] [,Intercept])
Finds least squares fit of a vector y to a vector or matrix x for both polynomial or multivariate fits, without uncertainty estimates. Optional Intercept forces the fit through the spe-cified y-value Intercept at x = 0. Passes all NIST linear least squares test sets with the default setting DgtMax = 30.
Least squares result .xRegLin _Eval(coeff, x [,DgtMax])
Computes the values of ycalc based on the coefficients found (e.g., with .xRegLin_ Coeff) for the specified set of x-values
Functions operating on complex numbers
For functions using complex arithmetic, the prefix MP.Cplx is required. Complex numbers, here denoted by z = a + jb, where j = ((1), must be defined in terms of their separate, real and imaginary components, a and b. Note that many of these instructions have the same form as those for real numbers, but are distinguishable by their prefix Cplx. For their double-precision equivalents, in the instruction just change the prefix x to s.
Negation .xNeg(z [,DgtMax]) z
= (a + jb) = a jb
Absolute value .xAbs(z [,DgtMax]) (z(
= EMBED Equation.3 = EMBED Equation.3
Inversion .xInv(z [,DgtMax]) 1/z
= EMBED Equation.3 = EMBED Equation.3
Conjugate .xCon(z [,DgtMax]) z*
= a jb
Zero check .xIsZero(z, tiny [,Dgt Max]) z ( 0
(z( = EMBED Equation.3 ( tiny
Addition .xAdd(z1,z2 [,DgtMax]) z1+z2
= (a1 + a2) + j (b1 + b2)
Subtraction .xSub(z1,z2 [,DgtMax]) z1z2
= (a1 a2) + j (b1 b2)
Multiplication .xMult(z1,z2 [,DgtMax]) z1 z2
= (a1a2b1b2) + j (a1b2+a2b1)
Division .xDiv(z1,z2 [,DgtMax]) z1/z2
= EMBED Equation.3
Integer power .xPow(z,n [,DgtMax]) zn
= EMBED Equation.3
Integer root .xRoot(z,n [,k] [,DgtMax]) z1/n
= EMBED Equation.3 , k = 0, 1, , n1
Square root .xSqr(z [,k] [,DgtMax]) z
= EMBED Equation.3 , k = 0, 1
Exponential .xExp(z [,DgtMax]) ez
=ea cos(a) + jeb cos(b)
Natural logarithm .xLn(z [,DgtMax]) ln z
= ln(a+jb) + arctan(a/b)
Decimal logarithm .xLog(z [,DgtMax]) log z
= ln(z) / ln(10)
Convert to polar .xPolar(z [,DgtMax])
z = EMBED Equation.3
Convert to rectangular .xRect(z [,DgtMax])
z = EMBED Equation.3
Evaluating an algebraic function
There are some practical restrictions with the above list. First, it may lack an operation you may want to perform. Secondly, while the elementary instructions can be nested to evaluate rather complicated algebraic expressions, this approach quickly becomes rather unwieldy. In order to calculate, e.g., the root d = [b+((b24ac)]/2a of the quadratic equation y = ax2 + bx + c one would need to write the instruction root = .xDiv(.xAdd(.xNot(b), .xSqr(.xSub(.xPow(b, 2), .xMult(4, .xMult(a, c))))), .xMult(2,a)), or, broken up, as discr = .xSub(.xPow (b, 2), .xMult(4, .xMult(a, c))), root = .xDiv(.xAdd(.xNot(b), .xSqr(discr)), .xMult(a, 2)), still quite clumsy. Therefore, a convenient function evaluator has been made available, in which f is a function in regular VBA format, and x is either a single number or a vector.
Function evaluation xeval(f[,x] [,Dgt Max] [,Angle]) f (x)
For trigonometric functions: Angle = Rad (default), Deg, or Grad (= Degrees ( 10/9)
x_eval("(1.2345^6.789)") ( 4.17958381726967068826150545493
x_eval("(sqr(5/2)+7^3.21)") ( 517.717800424940440699345256799
x_eval ("(sin(22.5)+3*cos(22.5))") ( 1.66810053720005875359979114908865584747919268
The following example (with the value of a placed inside quotation marks in order to preserve all of its more than fifteen digits) involves a single variable:
Sub Test3a()
Dim MP As New Xnumbers
Dim a, b
With MP
a = "1.234567890123456789"
b = .x_eval("(a ^3.456)", a)
Debug.Print "a = ", a
Debug.Print "a^3.456 = ", b
End With
Set MP = Nothing
End Sub
and yields, in the Immediate Window,
a = 1.234567890123456789
a^3.456 = 2.07145623088457274932940268363
Expressions containing multiple variables must enter the latter in a single matrix, with their names specified as labels in the first matrix column, and their numerical values in the second matrix column, as in
Sub Test3b()
Dim MP As New Xnumbers
Dim x(1 To 2, 1 To 2) 'dimension matrix
With MP
x(1, 1) = "a": x(1, 2) = .xE 'specify matrix elements
x(2, 1) = "b": x(2, 2) = .xPi
b = .x_eval("(a ^ b)", x)
Debug.Print "e = ", x(1, 2)
Debug.Print "Pi= ", x(2, 2)
Debug.Print "e^Pi = ", b
End With
Set MP = Nothing
End Sub
which produces the following output in the Immediate Window:
e = 2.71828182845904523536028747135
Pi= 3.14159265358979323846264338327
e^Pi = 23.1406926327793214925637779471
And here, finally, is a macro to compute the root of a quadratic equation with the usual (but, as we will see, sometimes not quite satisfactory) standard expression:
Sub Test4()
Dim a, b, c, root1, root2
Dim y(1 To 3, 1 To 2)
Dim MP As New Xnumbers
With MP
a = 1
b = 100000
c = 0.000001
y(1, 1) = "a": y(1, 2) = a
y(2, 1) = "b": y(2, 2) = b
y(3, 1) = "c": y(3, 2) = c
root1 = (-b + Sqr(b ^ 2 - 4 * a * c)) / (2 * a)
root2 = .x_eval("((-b + Sqr(b ^ 2 - 4 * a * c)) _
/ (2 * a))", y)
Debug.Print "root1 = ", root1
Debug.Print "root2 = ", root2
End With
Set MP = Nothing
End Sub
with the following results displayed in the Immediate Window:
root1 = 7.27595761418343E-12
root2 = 0.00000000001000000000000005
The difference between these two answers shows that, when 4ac b2, double precision has difficulties evaluating this rather common formula, because it must compute the difference between two nearly equal numbers, b and the square root of b2 4ac. In this case, the quadruple precision of Xnumbers is still sufficient to get the answer right. A better way, of course, would be to use an alternative formula that does not involve that difference, such as d = 2c/[b((b24ac)]. But such an obvious choice might not be available for more complicated problems, in which case extended precision would be a best first defense, before one finds out how to improve the algorithm used.
Summary
Here are the rules you must follow when writing new macros that use Xnumbers, or when modifying old ones.
(1) Give the macro a distinct name, such as one starting with an x.
(2) If you modify an already existing macro, keep a copy of the unmodified original. High-precision macros should not be used as routine updates of general-purpose macros, since they are too slow, even when set to the standard double precision of 15 decimals. This is so because much of the extra execution time is up-front, regardless of the precision specified.
(3) Before the dimensions are listed, insert the lines
Dim MP As New Xnumbers
With MP
where MP stands for multi-precision. The default (quadruple) precision yields 30 decimals, and need not be specified. If you want a different precision, specify that here instead, as with
DgtMax = 50
or, if you so prefer, with an input box that lets you specify the value of DigitsMax from the spreadsheet.
(4) Just before the end of the subroutine, insert the lines
End With
Set MP = Nothing
(5) In the dimension statements, remove the dimension specifiers As Single and As Double from any parameter you want to use in extended precision, because these dimension statements will overrule Xnumbers and keep the parameter in their dimension-specified precisions. By leaving their dimensions unspecified (but of course included in the list, otherwise Option Explicit would complain), these parameters will be assumed to be dimensioned As Variant, which is how Xnumbers can use them.
(6) Then modify all precision-critical instructions to fit the Xnumber format. For complicated expressions, either break up the expression into several simpler parts, or use the xEval function.
By comparing the macros in the xMacroBundle with their parents in the MacroBundle you will see that this need not involve very much work, because you need not rethink the structure of your code, merely some of its instructions. For more specifics, including additional functions and more details, see the Xnumbers_DLL_Ref manual.
*****
In summary, with xnumbers.dll and a few relatively minor modifications of your macros, subroutines, and functions, you can have high-precision computations while the structure of t * * + u v 7
8
; <
- .
% ' V o f g q r k l v 纳 ht B*ph ht 6] ht
ht CJ H*ht 6CJ NH ] ht 6CJ ] ht CJ NH
ht CJ ht 6B*]ph ht 6B*CJd ]ph B * v J f n U
$
x a$ $$
P x 1$ `a$gd!%
$x `a$gd!%
$x `a$ $$
pP@ x `a$ $
P x $$
P x a$ / % & . / P
$ % & ' A B G H O P U V W X ] ^ _ ` E F ˻˻˱˱˻˩˩ h!% CJ NH j h!% CJ H* j h!% CJ
h!% >*CJ
h!% CJ
h!% NH j h!% h!%
ht NH ht ht B*ph ht 6B*]ph BF ? @ I J f n o
# & U z * 2 S ͻ ht B*CJ OJ QJ ^J ph ht CJ OJ QJ ^J
ht NH ht ht 5CJ OJ QJ \^J ht B*CJ ph ht CJ NH
ht CJ ht 6B*]ph h!%
h!% >*CJ h!% CJ NH j h!% CJ
h!% CJ 3U b z * 2 ! :! m! $$
@ `x ^`a$ $$
@ `^`a$
$
x a$ $$
@ `x ^`a$ $$
@ `^`a$ $
`^`a$ ` ^` S T V W \ ] c d h i } ~ p! q! ! ! ! ! $" 4" 5" ;" <" [" " " " #
# *# 2# 7# C# # # # $ $ )$ p$ q$ $ $ 0% ټĭĭĭĭٝټټ ht 5CJ OJ QJ \^J ht B*CJ OJ QJ ^J ph ht 5CJ OJ QJ \^J ht CJ NH ht CJ OJ QJ ^J ht 6CJ ]
ht CJ j ht H*
ht NH ht 6H*]ht ht 6]8m! ! ! " $" [" " " " #
# &# 3# D# ^# x# # # # $$
@ `^`a$ x $$
@ `x ^`a$ $$
pP@ x x a$ $
a$ $$
pP@ x `a$ # # $ )$ & 9( F( ]( c( s( ( ( ( ( ( ( ) /) F) e) v) ~) `x ^` `^` $$
pP@ x a$ $$
pP@ x `a$ $$
@ `x ^`a$ 0% 1% % % & & }' ~' ' ' ' ' 9( i( r( }) ~) ) ) ) ) ) ) ) ) * * U* V* * * * * 0+ 1+ + + + + , , , , , =- H. T. . . "/ #/ / / 0 0 I1 J1 f1 g1 1 1 1 1 2 ,2 ľľľľľľľĴľĬ ht B*ph ht 56\]
ht NH ht ht 6CJ ]
ht CJ ht 5CJ OJ QJ \^J ht CJ OJ QJ ^J
ht CJ ht CJ NH @~) ) * * -* , =- J- a- j- z- - - - - - - . 1. H. `^`
$x `a$ $$
pP@ x a$ $$
pP@ a$ $$
pP@ `a$ $$
pP@ x `a$ H. T. k. . . . . . / 2 52 d2 2 2 2 2 2 3 $
pP8@ x
8@ x $
P8 x
8@ x
8@ $
pP@ l $
x `a$ `x ^` `^` ,2 -2 52 \2 ]2 2 2 2 2 2 2 2 2 2 2 2 2 3 &3 '3 33 43 53 63 L3 M3 Y3 Z3 [3 \3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 )4 *4 +4 ,4 .4 14 24 34 44 84 Q4 R4 _4 `4 f4 g4 h4 i4 k4 4 4 »鱪ҵ j hI hI 6] hI
ht CJ ht 6] j |ht ht 6CJ ]
ht CJ ht 6B*]ph
ht NH ht ht B*ph ht B*NH ph @3 73 ]3 3 3 3 3 .4 84 k4 4 4 4 5 ,5 65 _5 $
pP8@ gdI
8@ $
pP8@ l x
8@ gdI $
pP8@ x
@ x
$
pP8@ $
pP8@ x gdI 4 4 4 4 4 4 4 4 4 4 4 4 4 4 5 5 5 (5 )5 *5 /5 05 65 M5 N5 Z5 [5 ]5 ^5 e5 f5 {5 |5 5 5 5 5 5 5 5 5 86 :6 <6 >6 6 6 6 6 6 6 Ӯڦڀ ht OJ QJ hI 6H*]ht 6B*]ph j ht B*ph ht B*ph ht H* hI 6] j hI hI ht 6H*]
ht CJ ht 6] ht ht 6CJ H*]ht 6CJ ] ht 6CJ ]
ht CJ 2_5 5 5 5 |6 6 7 87 p7 7 7 !8 M8 {8 8 8 8 9 L9 $
P8 x
$
pP8@ $
pP8@ x
@ x
$
P8
8@ `x ^`
8@ 6 6 6 6 6 6 6 7 7 7 7 7 87 7 7 7 7 7 7 7 7
8 8 8 8 68 78 J8 K8 d8 e8 x8 y8 8 8 8 8 8 8 8 8 8 8 9
9 9 9 89 99 I9 J9 d9 e9 v9 w9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 ht H*ht 6H*]ht 6B*]ph ht OJ QJ
ht CJ hI ht ht 6H*] ht 6]ML9 y9 9 9 9
: ,: d: : : : ; ;; Q; s; ; ; <
$
P8
$
pP8@ $
P8 x
x
$
pP8
@ x
$
pP8@ 9 9 : : : : : : : : : : : !: #: %: &: (: ): *: M: N: a: b: m: n: t: u: v: w: x: y: z: : : : : : : : : : : : : : : : ; ; ; ; !; "; &; '; (; ); 1; 2; 7; 8; ;; Q; c; d; p; q; ; ; ; ; ;
ht NH ht 6B*]ph ht CJ EH ht CJ EH j ht ht H*ht 6H*] ht 6] ht
ht CJ H; ; ; ; ; ; ; < < < < < < < < < < < < < < < = = = = >= @= L= N= b= d= p= r= = = = = = = = = = = = = = 4> T> V> Z> \> ^> > > > > > > > 8? |? ht OJ QJ ht B*ph j ht B*ph
ht CJ ht 6] ht 6OJ QJ ] ht 6B*]ph j ht EHU jcE
ht CJ UVaJ ht j ht U =< f< < = = = 6> > > 8? |? ?
@ =@ T@ }@ @ @ @ A
@ x
@ $
P8 x
8@
8@ x
$
pP8@
@ x
@ |? @ @ s@ t@ @ @ @ @ !A "A )A 2A 3A 5A 6A @A AA DA EA OA PA SA TA _A `A cA dA fA hA iA A A A A A A B
B B NB OB B B B B B B 0C 1C QC RC yC zC C C C C #D $D +D ,D iD jD kD lD tD uD D D D j ht ht 6CJ ]
ht CJ ht 6B*]ph h
ht NH ht 6B* ]ph ht B* ph
ht CJ ht 6] ht FA )A fA A OB }B B B B C FC nC C C C C $
pP8@ l < $
pP8@ x
@ x
$
pP8@ $
P8 x
@ @ ^@
@ $
p8@ D D D D D D D E E E E E 3E 4E GE HE IE JE KE nE oE pE qE rE E E E E E E E E E E E E E E E E F F F F F F F -F .F 2F 3F 4F 7F 8F 9F ?F @F ^F _F ǻֲಝht B*CJ ph j ht CJ
ht CJ H* ht H*ht 6CJ ] j ht CJ EHU jCE
ht CJ UVaJ j ht CJ U
ht CJ j |ht ht 6H*]ht ht 6]:C D E KE rE E E F 4F bF F F F $G WG zG G G H H 6H $
pP8@ l < x $
pP8@ x $
pP8@ l x
$
pP8@
$
@ x a$ _F `F aF bF ~F F F F F F F F F F F F F F F F F F F F F F F F F F F F F F F F F G G "G #G $G @G AG DG EG RG SG UG VG WG `G aG eG fG gG hG lG mG vG wG xG zG G G G G G ݳݫ ht 6H*] j ht
ht CJ H*
ht CJ ht B*CJ ph j ht ht H*ht ht 6B*CJ ]ph ht 6CJ ] j ht CJ CG G G G G G G G H H H H H H 6H jH |H }H ~H H H H I 7I KI LI dI eI fI gI jI I I I I I I I I I I I I I I I I I I ,J =J ?J ]J ^J _J `J cJ J J J J J J K 仯 ht 6B*CJ ]ph ht B*CJ NH ph ht B*CJ ph
ht CJ H*ht CJ NH ht 6B*]ph ht H*ht
ht CJ ht 6CJ ] j ht CJ ?6H jH I 7I jI I I I ,J cJ J J J (K YL L M
@ @ ^@ $
@ @ x ^@ a$ $
P8 x $
pP8@ x $$
pP8@ @ x ^@ a$
$
pP8@ K K LK MK dK eK K K K K K K L L L L L L L M M M 5M M M M M M M M M M M M M M M N N
N N )N *N N N O
O O O O O O O "O $O IO JO KO LO QO RO eO fO gO hO iO jO j' ht EHU jNE
ht CJ UVaJ j ht U j |ht j ht ht 6B*]ph ht 6H*]
ht NH ht 6] ht
ht CJ BM 5M N O &O MO O O O O P 4P ^P P P P P Q
@ l x
@ x $
pP8@
@ x
$
pP8@
$
@ x a$
$
P8 jO }O ~O O O O O O O O O O O O O O O O O O O O O P P !P "P .P /P 0P 1P :P ;P U IU JU YU ZU ]U _U `U cU ht H* j ht ht 6B*]ph jZ ht EHU j E
ht CJ UVaJ jD ht EHU j E
ht CJ UVaJ j ht U
ht NH ht 6] ht
ht CJ <`S S S S S +W fW W W 3X X &Y 3Y JY SY ~ ~ ~
8@ l `^`
8@ l x ^
8@ l ^
$
pP8@ $
x `a$ $
pP8@ l x
@ x $
pP8@
@ x cU eU hU iU V V V V V V W W JW KW `W bW cW dW fW W W W W W W W W W W W W W W X X X X 3X VX WX X X X &Y zY Y Y PZ $[ H[ q[ y[ [ ĶĶīĶĶīĕĶ!ht 5B*OJ QJ \^J ph ht 6B*]ph j ht B*ph j ht B*ph ht B*OJ QJ ^J ph ht B*ph j ht B*ph ht B*CJ ph ht B*ph
ht CJ
ht NH ht ht 6]4SY [Y vY Y Y Y Y Y Y Z &Z PZ $[ 1[ H[ q[ y[ [ [ [
\
8@ l `^`
8@ dl `^`
8@ l `^`
8@ l x x
8@ l `^` [ [ [ f\ ] ] ^ =^ ^ ^ _ e_ %` '` (` )` +` ,` 9` :` ` ` ` ` ` ` ` ` ` a a a a a a a a a a a a a a a a b b b c c c c ͶججͶجب؈
ht NH ht 6B*CJ( ]ph ht H* j ht ht 6] ht ht B*NH ph ht B*H*ph ht 6B*CJ ]ph ht 6B*]ph ht B*ph !ht 5B*OJ QJ \^J ph ht B*OJ QJ ^J ph 3
\ )\ B\ K\ \\ f\ \ \ \ ] ] ] ] ^ ^ ^ %^ 0^ =^ X^ s^ ^ ^
8@ l x x `
8@ l `x ^`
8@ l x x
8@ l `^` ^ ^ _ &_ D_ M_ ^_ f_ _ _ _ b b c Ic d d $x a$ $x x a$ $
8@ l x x `a$ $
8@ l `x ^`a$ $
8@ l `x ^`a$
8@ l x x
8@ l `^` c d d d
e ae be e e e e $f %f :f ;f f f f f f f f f =g >g Dg Eg g g 'h (h Sh ]h h h mi ni Ij Jj j j Jk U V / Uht CJ OJ QJ ^J ht OJ QJ ^J
ht NH ht 3d e
e e e Df f f f h Fi j j . /
8pPl
@
$
@ a$ $x x a$ $a$ $x a$
$`x ^`a$ $`^`a$ hese routines, as well as your data input and output formats, remain unchanged. This can provide a very useful back-up for any custom macros, as it allows you to check the accuracy of your code. It can also help you identify areas of code where it is worthwhile to pay extra attention to the algorithm used. Some of the macros of the MacroBundle have already spawned x-versions which you can find in the xMacroBundle. You are of course welcome to try them out. Please let me know by e-mail, at rdelevie@ bowdoin.edu, if you encounter any problems with them.
, 1h / =!"#$% D d
J
C A ? " 2 N41T>#5` D `! N41T>#5`& 8> xڕ?Ka繳SH,BNB
R= ԐNm
SC4-p2TM!]?Sj={b )d1K$$Tu+bRnfΐdcߏ&@Qg1'^Uw eO|1Gȣ=OEѹ4Ǥ7j b{"$K-wv ^)=>)+Ush7j9&9KP|GA|A9.%djL1ub @= S$ʖZr2#̍o){H?꼂l D d
J
C A ? " 2 AC@CdSLG `! AC@CdSLG @ 0= xڥS=H@~ԦKu]A lVhŴB AW{ g .
[E%P %{w =0bDa?ad4u. M3F?'Ȑѩ@QCk4'v필x ua)ADlT G2o%X/Hqn|Qjg<
bsGY\(3E3LwH[P{I Dž3 K@&Ñ$>Rp~r=Yrˠn% D d
J
C A ? " 2 {F*폏|) c k `![ {F*폏|) (+ h ) xcdd`` $d@9`,&FF(` Ti A ? d- P=7T
obIFHeA*CPD 32Bab ԘMO_-CۺnuGvEε+[V ul0y_>rx}ǐYv:
%%օk ?D+}1G\X E8.p0MfdogM+ D d
J
C A ? " 2 mxp@ $ i `!a mxp@ $ & ` P d / xcdd``f 2
ĜL0##0KQ*
Wä"d3H1)f Y;X3PT
obIFHeA*/& *fe- f Y, wfjQ9B&br(b5~=Wͷ zzf7zH܌{ fj Wv Mމ?CemU D d
J
C A ? " 2
bsUEa] a `!
bsUEa] ~ h Y xڕRJ@$M
Q *} Qzy
-ĽH#!R)|q'+ 0o/
6,\GrD!eYPLZj4!`.a2I
+֘_5@F =(st\6j&]Bt?hfhn=ڨk9GG'VBOGq4? >MO_-Cۺ
! " # $ % & ' ( ) * + , - . / 0 1 2 3 4 5 6 7 8 9 : ; < = > ? @ A B C D E F G H I J K L M N O P Q R S T U V W X Y Z [ \ ] ^ _ ` a b c d f g h i j k l o r s t v u w x y { z | ~ } R o o t E n t r y A F j|<q D a t a
e W o r d D o c u m e n t @ . O b j e c t P o o l C ' ]^|<j|< _ 1 1 6 7 7 4 6 0 3 7 F ]^|<_|< O l e
C o m p O b j f O b j I n f o
" % ( + , - . 1 4 5 8 ; > A D G J M O P Q R S T V W X Y [
F Microsoft Equation 3.0 DS Equation Equation.3 9q *f
n m ( ) = n ! m ! ( n "m ) ! E q u a t i o n N a t i v e _ 1 1 6 7 8 0 3 2 9 3 ; F _|<_|< O l e
C o m p O b j
f
F Microsoft Equation 3.0 DS Equation Equation.3 9q *l h
a i , j j = 1 n
" i = 1 m
"
F Microsoft Equation 3.0 DS Eq O b j I n f o
E q u a t i o n N a t i v e _ 1 1 6 7 7 6 3 2 7 8 F _|< ka|< O l e
C o m p O b j f O b j I n f o E q u a t i o n N a t i v e D _ 1 1 6 7 7 6 3 3 2 5 " F b|<b|< uation Equation.3 9q *( <
a + j b
F Microsoft Equation 3.0 DS Equation Equation.3 9q O l e
C o m p O b j f O b j I n f o E q u a t i o n N a t i v e T *8 1
a 2
+ b 2
F Microsoft Equation 3.0 DS Equation Equation.3 9q *& p
1 a + j b _ 1 1 6 7 7 9 2 5 7 5 6 F ka|<b|< O l e
C o m p O b j f O b j I n f o E q u a t i o n N a t i v e B _ 1 1 6 7 7 9 2 5 8 5 F b|<b|< O l e
C o m p O b j ! f
F Microsoft Equation 3.0 DS Equation Equation.3 9q *Q F
"a b + j b a 2
+ b 2
F Microsoft Equation 3.0 DS Eq O b j I n f o ! # E q u a t i o n N a t i v e $ m _ 1 1 6 7 7 6 6 4 2 8 $ F b|<`xd|< O l e
& C o m p O b j # % ' f O b j I n f o & ) E q u a t i o n N a t i v e * _ 1 1 6 7 7 6 6 5 4 7 1 ) F `xd|<`xd|< uation Equation.3 9q * =
( a 1
a 2
+ b 2
b 2
) + j ( "a 1
b 2
+ a 2
b 1
) a 2 2
+ b 2 2 O l e
/ C o m p O b j ( * 0 f O b j I n f o + 2 E q u a t i o n N a t i v e 3
F Microsoft Equation 3.0 DS Equation Equation.3 9q *
8 <
a 2
+ b 2
e x p [ n a r c t a n ( a / b ) ]
F Microsoft Equation 3.0 DS Eq_ 1 1 6 7 7 9 1 9 8 7 . F `xd|< e|< O l e
6 C o m p O b j - / 7 f O b j I n f o 0 9 uation Equation.3 9q *' h
a + j b n
F Microsoft Equation 3.0 DS Equation Equation.3 9q E q u a t i o n N a t i v e : C _ 1 1 6 7 7 9 2 2 9 0 , 3 F e|< e|< O l e
< C o m p O b j 2 4 = f O b j I n f o 5 ? E q u a t i o n N a t i v e @ C _ 1 1 6 7 7 9 3 4 2 0 8 F e|<g|< O l e
B *' 8
a + j b
F Microsoft Equation 3.0 DS Equation Equation.3 9q *' 8
e
C o m p O b j 7 9 C f O b j I n f o : E E q u a t i o n N a t i v e F C _ 1 1 6 7 7 9 3 6 1 2 = F g|<g|< O l e
H C o m p O b j < > I f O b j I n f o ? K E q u a t i o n N a t i v e L i
F Microsoft Equation 3.0 DS Equation Equation.3 9q *M p
c o s ( ) + j s i n ( ) Oh +'0 nuGvEε+[V ul0y_>rx}ǐYv:
%%օk ?D+}1G\X E8.p0MfdogM D d
J
C A ? " 2 ^&Z0r `! ^&Z0r. ` P xڝ=HP]b?jP!*(C `n BPqru(C $`| iAMy./y Z7:c;BƼ2e=\EqV 1W> \Dx(LJJB|
UFZ;eMB;L|j5\+/Ena5/m&rFPKGDjQuAm\ n `! $ _ukV>\R xڝRJP>礿ii-E!
E n(bIT
-ZHuG7 DpT\ ''87 rs{ %`
3KSb!/
NL9yʉ Lɠ11Y?up&&X bQjׄ7NW b3)Cg@e,x.gAO|B@kY桾aV2pu?=fhT<$j@g |3 ZjfLjyNh;1CW"Q=
ŹD ,IZ>+5~1Zr$̖w~m9 ~ u{} ǩ?}uӇ6ڂ%o??K ߔ ,y'Ao+n0GN<^;) u bifZ()wô?@+d0 #Xƥ[">{ ǓgT^ib? D d
J
C A
? "
2 X7Ғw-kGo } I `!u X7Ғw-kGo* h C xmJPƿsM tE E`cBp)9(,?)Z=!g P$QB1UykFVqU pVjLLcK)|2Uv41%O6'r%Vxh5gK=°u|bhJrWRe8{Q4q$X}XۆѣFsvj4Hp]wW&+гSI^>&u1/^?Ua Vo 5=
.GW D d
Xh J
C A ? " 2 x [GZ%DF T `!L [GZ%DF @ | xcdd`` @c112BYL%bpu 13*ߍ/b/I~׀ dBV>0!V\28' F&&\A z> 1. a. . . . . . !/ M/ {/ / / / 0 L0 y0 0 0 0
1 ,1 d1 1 1 1 2 ;2 Q2 s2 2 2 3 33 y3 3 3 3 4 O4 f4 4 4 4
5 =5 T5 }5 5 5 5 6 )6 f6 6 O7 }7 7 7 7 8 F8 n8 8 8 8 8 9 : K: r: : : ; 4; b; ; ; ; $< W< z< < < = = 6= j= > 7> j> > > > ,? c? ? ? ? (@ YA A B 5B C D &D MD D D D D E 4E ^E E E E E F 9F `F }F F F F G EG kG G G G G "H 7H `H H H H H +L fL L L 3M M &N 3N JN SN [N vN N N N N N N O &O PO $P 1P HP qP yP P P P
Q )Q BQ KQ \Q fQ Q Q Q R R R R S S S %S 0S =S XS sS S S S T &T DT MT ^T fT T T T W W X IX Y Y Z
Z Z Z D[ [ [ [ ] F^ _ _ xb {b 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 @ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 J {b ۊ 0 DlZي 0 ي 0 A B G H ۊ 0 D F S 0% ,2 4 6 9 ; |? D _F G K jO P yQ R cU [ c / 6 9 : < ? B D F H I K M O P R T U W X Z \ _ U m! # ~) H. 3 _5 L9 < A C 6H M Q `S SY
\ ^ d / 7 ; = > @ A C E G J L N Q S V Y [ ] ^ ` / 8 2 2 2 3: G: I: QD eD gD iD }D D D D D D D D ?E SE UE dF xF zF F F F F G G IG ]G _G gH {H }H H H H yb ::::::::::::: l , b $ =W;e]vL٧
@ 0 (
B
S ? yb f
r
?
E
G
S
Y
i
J C N Q \ @ L
+ / 8 < D O ^ i x " s |
! ! *! F! V! z% % % % % % % % % % % % % & & ,& k& {& & & & & P( ]( * * * * * * * * * * * * + "+ #+ &+ )+ 0+ 3+ 6+ D+ H+ I+ L+ O+ V+ d+ h+ n+ r+ + + + + + + + + + + + + + + + + + + + + , , , , , !, $, ', -, 8, I, M, N, Q, T, [, z, ~, , , , , , , , , , , , , - 6- D- I- J- M- P- W- w- z- - - - - . . I. M. N. O. R. Y. k. p. q. r. u. |. . . . . . . . . . . / / / / // 4/ 9/ @/ C/ I/ ]/ b/ g/ n/ q/ w/ / / / / / / / / / / 0
0 0 0 0 0 20 60 ;0 B0 ]0 b0 g0 n0 q0 u0 {0 0 0 0 0 0 0 0 0 0 0 0 0 0 E1 K1 P1 W1 Z1 `1 f1 l1 q1 s1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 +2 -2 \2 a2 f2 m2 y2 2 2 2 2 2 2 2 33 93 c3 e3 j3 l3 ~3 3 3 3 3 3 3 3 T4 W4 Z4 a4 i4 p4 4 4 4 4 5 5 5 5 ,5 05 K5 O5 Q5 Q5 5 5 5 5 5 5 5 5 5 5 5 5 6 6 ,6 16 6 6 7 7 7 7 7 7 7 7 8 8 8 &8 I8 P8 q8 x8 8 8 8 8 9 9 9 : : : 3: J: Y: `: e: k: : : : : ; ; #; ); s; |; ; ; ; ; ; ; 5< >< H< N< < < < < < < < < < < I= Q= _= e= > '> ,> 2> :> B> > > > > ? ? "? (? r? ~? ? ? ? ? ? ? ? @ @ @ @ @ LA RA oA |A }A A A A A A A A hB oB B B C C C C C D D D "D $D 6D :D ?D ED QD hD iD D D D D D D D D D D D D D E E E E $E +E ?E VE hE lE uE {E E E E E E E F F CF GF PF VF dF {F F F F F F F F F F F F F F F F F F G *G .G 8G >G IG `G xG |G G G G G G G G G G G G G G G G G G G G G G G H H H H 'H )H /H 1H IH OH TH ZH gH ~H H H H H H H bI bI cJ eJ J J J J J J J J J J J J J J J J J J K K
K K K K !K &K +K 0K BK FK HK LK NK RK XK \K ]K bK gK lK ?L DL ML TL L L L L L L L M 3M 9M IM LM {N N N N N N P P P P P P P P
Q Q )Q 4Q 6Q :Q Q Q S S S S S S T T &T 1T xZ Z Z Z "[ +[ [ [ ] ] 6^ ;^ e^ q^ w_ _ a a 8b @b Bb Mb {b J e v x ' / Y a 2 7 : ; $ ' + 0 8 = = E M O i k ! F! W! " " -" 9" A% I% Q% S% p% r% % % % % % & k& |& & & 8* =* g* j* * * * * * * + #+ D+ I+ _+ a+ + + + + + + , , 1, 2, I, N, z, , , , , , - - - - /- 0- D- J- b- d- - - . . I. N. k. q. . . . . . . / / // 5/ ]/ c/ / / / / 0 0 20 70 ]0 c0 {0 0 0 0 0 0 0 0 1 1 E1 L1 f1 l1 1 1 1 1 1 2 2 2 \2 b2 u2 x2 2 2 c3 f3 ~3 3 3 3 3 3 3 3 4 !4 T4 X4 4 4 4 4 5 5 K5 P5 V5 ^5 5 5 5 5 5 5 5 5 6 6 ,6 26 6 6 97 ?7 d7 j7 7 7 7 7 7 7 8 8 I8 Q8 q8 y8 8 8 9 : Y: a: : : : : ; ; s; }; ; ; 5< ?< Z< ^< < < < < O= R= %> (> :> Q> > > > > ? ? E? J? r? ? ? ? @ @ j@ |@ xA }A C C D D 6D ;D D D D D E E hE mE E E E F CF HF F F F F G G *G /G xG }G G G G G G G H
H 'H *H IH PH H H ?L EL L L L L =M AM *N 2N :N '> *> *> :> R> g> i> > > ? ? ? ? ,? c? ? ? ? ? QD hD iD D D D D D ?E VE dF {F F F F G IG `G gH ~H H H bI bI J J J J $P $P P P W W
Z
Z 0Z 3Z TZ [Z ]Z Z Z Z Z Z Z Z Z Z 7^ 8^ >_ r_ d` d` xb {b J O {b t !% B I @ J J hJ J J L . / 3 5 J` yb p @ p 6 p p @ p < p @ p @ U n k n o w n G : T i m e s N e w R o m a n 5 S y m b o l 3& : A r i a l ; W i n g d i n g s ?5 : C o u r i e r N e w " 1 h FFX} S 2 S 2 ! 24 d Gb Gb 2Q HP ) ? I 2 F u n c t i o n R o b e r t d e L e v i e I T