RECURRENT SEQUENCES IN THE EQUATION DQ 2 = R 2 + N INTRODUCTION

Similar documents
Table of Contents. 2013, Pearson Education, Inc.

FIFTH ROOTS OF FIBONACCI FRACTIONS. Christopher P. French Grinnell College, Grinnell, IA (Submitted June 2004-Final Revision September 2004)

Complex Numbers: Definition: A complex number is a number of the form: z = a + bi where a, b are real numbers and i is a symbol with the property: i

PELL S EQUATION, II KEITH CONRAD

Problem Set 5 Solutions

CONVOLUTION TREES AND PASCAL-T TRIANGLES. JOHN C. TURNER University of Waikato, Hamilton, New Zealand (Submitted December 1986) 1.

CLOSED FORM CONTINUED FRACTION EXPANSIONS OF SPECIAL QUADRATIC IRRATIONALS

USING SHANKS BABY-STEP GIANT-STEP METHOD TO SOLVE THE GENERALIZED PELL EQUATION x 2 Dy 2 = N. Copyright 2009 by John P. Robertson. 1.

PELL S EQUATION NATIONAL INSTITUTE OF TECHNOLOGY ROURKELA, ODISHA

EQUATIONS WHOSE ROOTS ARE T I E nth POWERS OF T I E ROOTS OF A GIVEN CUBIC EQUATION

CONGRUUM PROBLEM. Manju Somanath 1 and J. Kannan 2. National College, Trichy - 01, India

Algebra 2 (2006) Correlation of the ALEKS Course Algebra 2 to the California Content Standards for Algebra 2

Algebra. Table of Contents

Warm-Up. Simplify the following terms:

A Guide to Arithmetic

P.5 Solving Equations

ZEROS OF POLYNOMIAL FUNCTIONS ALL I HAVE TO KNOW ABOUT POLYNOMIAL FUNCTIONS

E.J. Barbeau. Polynomials. With 36 Illustrations. Springer

THESIS. Presented in Partial Fulfillment of the Requirements for the Degree Master of Science in the Graduate School of The Ohio State University

Properties of Radicals

Cryptography. Number Theory with AN INTRODUCTION TO. James S. Kraft. Lawrence C. Washington. CRC Press

ALGEBRA I FORM I. Textbook: Algebra, Second Edition;Prentice Hall,2002

Section 1.3 Review of Complex Numbers

Neverending Fractions

Pell's Equation. Luke Kanczes

Algebra I Unit Report Summary

Sums of Squares. Bianca Homberg and Minna Liu

Math 1302 Notes 2. How many solutions? What type of solution in the real number system? What kind of equation is it?

Solving x 2 Dy 2 = N in integers, where D > 0 is not a perfect square. Keith Matthews

COMPOSITIONS AND FIBONACCI NUMBERS

198 VOLUME 46/47, NUMBER 3

Complex Numbers. Essential Question What are the subsets of the set of complex numbers? Integers. Whole Numbers. Natural Numbers

Solving Quadratic Equations Review

Power series and Taylor series

DIVISIBILITY BY FIBONACCI AND LUCAS SQUARES

12x + 18y = 50. 2x + v = 12. (x, v) = (6 + k, 2k), k Z.

Sect Complex Numbers

Part II. Number Theory. Year

Perform the following operations. 1) (2x + 3) + (4x 5) 2) 2(x + 3) 3) 2x (x 4) 4) (2x + 3)(3x 5) 5) (x 4)(x 2 3x + 5)

Pell s Equation Claire Larkin

MATH98 Intermediate Algebra Practice Test Form A

The Diophantine equation x n = Dy 2 + 1

CONTINUED FRACTIONS, PELL S EQUATION, AND TRANSCENDENTAL NUMBERS

Exploring Number Theory via Diophantine Equations

Constructions with ruler and compass.

EXERCISES. a b = a + b l aq b = ab - (a + b) + 2. a b = a + b + 1 n0i) = oii + ii + fi. A. Examples of Rings. C. Ring of 2 x 2 Matrices

Notes on Systems of Linear Congruences

In memoriam Péter Kiss

8. Dirichlet s Theorem and Farey Fractions

Advanced Counting Techniques. Chapter 8

Introduction to Number Theory

QUOTIENTS OF FIBONACCI NUMBERS

Algebra. Mathematics Help Sheet. The University of Sydney Business School

Chapter 4 ARITHMETIC AND GEOMETRIC PROGRESSIONS 2, 5, 8, 11, 14,..., 101

Day 6: 6.4 Solving Polynomial Equations Warm Up: Factor. 1. x 2-2x x 2-9x x 2 + 6x + 5

1 Basic Combinatorics

Course Name: MAT 135 Spring 2017 Master Course Code: N/A. ALEKS Course: Intermediate Algebra Instructor: Master Templates

Diophantine quadruples and Fibonacci numbers

Diophantine equations

Necessary and Sufficient Conditions for the Central Norm to Equal 2 h in the Simple Continued Fraction Expansion of 2 h c for Any Odd Non-Square c > 1

ON THE POSSIBLE QUANTITIES OF FIBONACCI NUMBERS THAT OCCUR IN SOME TYPES OF INTERVALS

2005 Palm Harbor February Invitational Algebra II Answer Key

SOME REMARKS ON THE DIOPHANTINE EQUATION x3+y3+z3 = x+y+z

Skills Practice Skills Practice for Lesson 4.1

Theorems About Roots of Polynomial Equations. Rational Root Theorem

PYTHAGOREAN TRIPLES KEITH CONRAD

MATH Fundamental Concepts of Algebra

Mathematics (MEI) Advanced Subsidiary GCE Core 1 (4751) May 2010

THE GOLDEN RATJO AMD A GREEK CRJSJS*

Solving Quadratic Equations by Formula

ADVANCED PROBLEMS AND SOLUTIONS

Section 4.1: Sequences and Series

MATH98 Intermediate Algebra Practice Test Form B

At the end of this section, you should be able to solve equations that are convertible to equations in linear or quadratic forms:

THE CONGRUENT NUMBER PROBLEM

Generalization of Fibonacci sequence

#A5 INTEGERS 17 (2017) THE 2-ADIC ORDER OF SOME GENERALIZED FIBONACCI NUMBERS

Advanced Elementary Algebra

LEGENDRE S THEOREM, LEGRANGE S DESCENT

Algebra I. Book 2. Powered by...

Pure Mathematics P1

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0

QUADRATIC EQUATIONS. + 6 = 0 This is a quadratic equation written in standard form. x x = 0 (standard form with c=0). 2 = 9

Unit 5 Evaluation. Multiple-Choice. Evaluation 05 Second Year Algebra 1 (MTHH ) Name I.D. Number

Gaussian integers. 1 = a 2 + b 2 = c 2 + d 2.

SECOND-ORDER RECURRENCES. Lawrence Somer Department of Mathematics, Catholic University of America, Washington, D.C

Lesson 3.5 Exercises, pages

Extended Binet s formula for the class of generalized Fibonacci sequences

Chapter 9 Factorisation and Discrete Logarithms Using a Factor Base

AROUND (1 + 2) n. (1 + x) 5 = 1 + 5x + 10x x 3 + 5x 4 + x 5.

MATHEMATICS AND GAMES* Nagayoshi lwahori University of Tokyo

AMBIGUOUS FORMS AND IDEALS IN QUADRATIC ORDERS. Copyright 2009 Please direct comments, corrections, or questions to

ADVANCED PROBLEMS AND SOLUTIONS

Divisibility in the Fibonacci Numbers. Stefan Erickson Colorado College January 27, 2006

MCR3U Unit 7 Lesson Notes

#A11 INTEGERS 12 (2012) FIBONACCI VARIATIONS OF A CONJECTURE OF POLIGNAC

Notes on Continued Fractions for Math 4400

Chapter Five Notes N P U2C5

RATIONAL POINTS ON CURVES. Contents

ON SUMS AND RECIPROCAL SUM OF GENERALIZED FIBONACCI NUMBERS BISHNU PADA MANDAL. Master of Science in Mathematics

Transcription:

RECURRENT SEQUENCES IN THE EQUATION DQ 2 = R 2 + N EDGAR I. EMERSON Rt. 2, Box 415, Boulder, Colorado INTRODUCTION The recreational exploration of numbers by the amateur can lead to discovery, or to a different way of looking at problems, because he often does not know the conventional approaches* Sometimes, as a form of amusement, I picked a quadratic expression at rom, set it equal to a square then tried to solve the resulting equation in positive integers. Whenever I was able to solve the problem I noticed that recurrency was evident. One of the most satisfying results came from the solution of 5x 2 =t 6y + 1 = y 2 where the recurrent relationships involved Fibonacci Lucas sequences. However, the method reported [ l j for this solution is not general. An improvement in the method resulted from exploring the Pell Lagrange* equations* As experimental data accumulated I was able to make some conjectures when I d i s - cussed the results with my friend, Professor Burton W, Jones, he urged me to try to prove them. For his encouragement, l a m grateful. The following are some of these conjectures: a) For any recurrent equation such as U - = cu + U - or U - = cu - U -, c constant even, there exists at least one Pell equation such n n-1 ^ that the sequence of X f s of Y f s follow the given recurrent law. b) In a Pell equation if DY = X\ + 1 then the recurrent law for the sequence of X's or Y f s is U n + 1 = cu n + U n _ 1 if DY = X*- 1 then the governing law is U - = cu - U -. c) In Lagrange equations having the same D as a Pell equation, there exists a recurrent law common to both. communication.) (Proof to be offered in another *The Lagrange equations Dy 2 communication. = x 2 ± N, N > 1 will be discussed in another 231

232 RECURRENT SEQUENCES IN THE PELL EQUATIONS [Oct. Since a method of developing the sequence of one of the variables, in a Pell equation, independent of the other is so easy since the proof justifying such treatment uses only elementary algebra, without the use of continued fractions or convergents, I thought that the method might be of interest. As will be demonstrated, problems, relating to the Pell equations which seem difficult, are solved in an almost trivial fashion by means of the theorems to be developed here. (Before continuing the reader is invited to try solving problems 1-5.) PART 1 - THE PELL EQUATIONS DY* = X* -(-l) n DY* = X* - 1 For a given D > 1 not a square the complete* Pell equations are either of the forms (1) DY^ = X^ - (-l) n or (2) DY^ - X 2 n - 1, n = 0, 1, 2, 3,.-. While both of these equations have the trivial solution X 0 = 1, YQ = 0, the key to the general solution is in finding X l9 Yj, either by inspection or otherwise. How this may be done by convergents is explained by Burton W. Jones p ], C. D. Olds [3], R. Kortum G. M c N e i l ^ others in books on number theory. The least positive, non-trivial solution '(X^Yj) is variously called the minimal or fundamental or generating solution. Once this solution is found, the general solution is given by (3) X n + Y n ^ = ( x i + YiA/D) n. *The equation DY* = X* - (-l) n, n = 0,1,2,3, #, is complete. However, it is commonly treated as two equations, e. g., DY, = X, - 1 DY, = X, - + 1, k = 0,1,2,3,-. Unless otherwise stated, we will assume that for the given D, the Pell equation is complete we are dealing with all possible solutions.

1969] RECURRENT SEQUENCES IN THE PELL EQUATIONS 233 The sum of the rational terms in the binomial expansion of (Xj + Yi^O)* 1 is X n the sum of the irrational terms is Y VD* That equation (3) gives all of the possible solutions was first shown by Robert D. Carmichael later explained in his book Diophantine Analysis [5 ]. When the minimal solution (X l9 Yi) is substituted in equations (1) (2) we have respectively the minimal equations (4) DY 2 = X 2 + 1* (5) DY 2! = X 2 i - 1. In either case, irrespective of the sign preceding 1, the general solution is given by the single equation (3). PROOF OF THREE THEOREMS ON RECURRENCY IN THE PELL EQUATIONS Theorem 1. In the integer solution of a Pell equation, the sequence of X f s is recurrent as is the sequence of Y f s according to the recurrent law, U - = cu ± U 1, c = 2Xj. The + sign is used if the minimal equation is DY 2 = Xi + 1 the - sign is used if DY 2 = X 2-1. To prove this theorem we combine the minimal equations (4) (5) so that (6) DY 2 = Xi dt 1. Then for reference we prepare, from the general solution (3), the following set of equations: (7a) (Xi + YiVD) 1 1 " 1 = X n _ 1 + Y n 4 ^ *If the minimal equation for a certain D is DY 2 = X 2 + 1 then there are s o - lutions for DY 2 = X 2 ± 1. If the minimal equation is DY 2 = X 2-1 then the only solutions are for DY 2 = X 2-1. Thus DY 2 = X 2 + 1 is not solvable for all D f s nor does it have a trivial solution.

234 RECURRENT SEQUENCES IN THE P E L L EQUATIONS [Oct. (7b) (Xi + YiN/D) n = X +Y N/D (7c) (Xi + Yt N/D ) n + 1 = X n + 1 + Y n + 1 N/D. When X 2! + 2XiY A N/D is added to both sides of DY\ = x\ ± 1 we obtain x } + 2XtYi N/D + D Y I = 2X 2 i + 2XiYt N/D ± 1 o r (8) (Xj + Y i V D ) 2 = 2Xi(Xi + Y W D ) ± 1. n 1 Multiplying both sides of this equation by (X A + Yt \/D) ~ we derive (9) (Xi + Yi *T5 ) n + 1 = 2Xi(Xi + Yi N/D f dt (X t + Y t N/D J 11 " 1. Now when the appropriate substitutions a r e made in this equation from s e t (7) we get X _ + Y ^ VD - 2Xj(X l + Y N / D ) ± (X - + Y, *1T5) n+1 n+1 n n n - 1 n - 1 rearranging this equation we have ( 1 0 ) X n + 1 + Y n + 1 ^ = { 2 X i X n * X n - 1 > + < 2 X i Y n * V l W ' After equating the rational then the irrational t e r m s in (10) we finally derive (11) X '= 2XjX + X, n+1 x n - n - 1 ( 1 2 ) Y n + 1 = 2 X i Y n ± Y n - l Thus the proof of T h e o r e m 1 is complete equations (11) (12) a r e the equations of the Theorem. *The equations of the Theorem s e e m s i m i l a r to expressions found for the convergents of continued fractions. F o r instance, the n u m e r a t o r of the n convergent is p = a p 1 + p 0. This equation s e e m s s i m i l a r to X = cx - n n n x n & n n x ± X 0 but in the equations of T h e o r e m 1, + o r - signs a r e used w h e r e a s n u in the convergent expression only the + sign a p p e a r s.

1969] RECURRENT SEQUENCES IN THE PELL EQUATIONS 235 As a consequence of Theorem 1 we have Theorem 2. For every recurrent equation, U,- = cu +U., or U n+1 = c U n " U n 1 w n e r e c i s e v e n * there exists at least one Pell equation for which the sequence of either variable is governed by the given recurrent law. To prove this theorem we note from Theorem 1 that c = 2X^ whence Xi = c/2. When this value of Xj is substituted in the minimal equations DY 2 i = Xj ± 1 we have DYi = m ± i ( «) ' * Except for a trivial case 9 (t) " ± i *. therefore we can let ( * ) 2 ± 1 = D whence Yi = 1 thus we have proved Theorem 2, If (s) f * > contains a square factor >1 there may be other solutions as demonstrated by problem 1. In equation (1), DY 2 = X 2 - (-l) n, we notice that when n = 2k then (13) DY k = X k - 1 when n = 2k + 1 then (14) D Y k+l = X k + 1 + 1 > k = O. 1 ' 2 ' 3 ' ' "

236 RECURRENT SEQUENCES IN THE PELL EQUATIONS [Oct. In order to study the sequence of every other term in a Pell equation we have Theorem 3. The sequence of every other X or Y in a Pell equation in recurrent. If the recurrent law for the Pell equation is U ^ = c U +U - ^ n+1 n n-1 then the sequence of every other X or Y is U n + 3 = (o» + 2)U n + 1 - U n _ x if the recurrent law is U - = cu - U 1 then the sequence of every other X or Y is governed by U n + 3 " (c2 " 2 ) U n " U n - 1 ' sign. We prove the two parts of Theorem 3 together using the ambiguous ± U n + 1 = c U n ± U n-1 then U n + 2 = c U n + l ± U n U n + 3 - c U n + 2 * U n + 1 But U n + 2 - c U n + l * U n therefore U n + 3 = c ( c U n + l ± U n ) ± U n + l or U a + 3 - c 2 U n + l * c U n * U n + 1

1969] RECURRENT SEQUENCES IN THE PELL EQUATIONS 237 U n + 3 = < c 2 * 1 ) U n + l ± c U n But ±cu = ztu _,_- - U - n n+1 n-1 therefore U n+3 = ( o ^ D U ^ i U ^ - U ^ or ( 1 5 ) U n + 3 = < c2 ± 2 > U n + l - U n-1 ' With the derivation of equation (15) we have proved Theorem 3. For convenience we let c 2 ± 2 = c 2 then the equations of Theorem 3 become ( 1 6 ) U k+1 = C2U k " U? k - 1 ' U l = U ' U f l = U 2 or Ui = U l s W 2 = U 3. The method of proof for Theorem 3 demonstrates that the properties of the sequences of X f s or of Y f s in the Pell equations are simply the properties to be expected from considerations of the recurrent equations U n+1 c U n ± U n - r EXAMPLES 2 2 Example 1, When D = 2 the minimal solution is 2Y* = X 4 + 1, Yi = 1 9 Xj = 1. From Theorem 1 we know that we must use the recurrent equation with the + sign that the constant c = 2Xj = 2. of X f s develops from X n + 1 = 2X n +X _ r X 0 = 1, X 4 = 1. Thus, the sequence

238 RECURRENT SEQUENCES IN THE P E L L EQUATIONS [Oct. X 2 = 2Xi + X 0 = 2 1 + 1 = 3 X 3 = 2X 2 + Xj = 2 3 + 2 = 7 X 4 = 2X 3 + X 2 = 2 7 + 3 = 17, etc. Thus X = 1, 1, 3, 7, 17, 4 1, 9 9, -. Similarly for Y we have Y n + 1 = 2 Y + Y ^ Y 0 = 0, Y* = 1. Y 2 = 2Yi + Y 0 = 2 1 + 0 = 2 Y 3 = 2Y 2 + Yi = 2 2 + 1 = 5 Y 4 = 2Y 3 + Y 2 = 2 5 + 2 = 12, etc., Y = 0, 1, 2, 5, 12, 29, 70, ". Example 2. F o r D = 3 the minimal solution is Xj = 2, Y* = 1 2 2 the minimal equation is 3Yi = X A - 1, whence the r e c u r r e n t law for D = 3 is U n + l = c U n- U n-l' c = 2X1 = 2-2 = 4. Then X 2 = 4X t - X 0 = 4 2-1 = 7 X 3 = 4X 2 - Xj = 4 7-2 = 26 X 4 = 4X 3 - X 2 = 4 26-7 = 99, etc., for the Y f s Y? = 4Yi - Y 0 = 4 1-0. = 4

1969] RECURRENT SEQUENCES IN THE P E L L EQUATIONS 239 Y 8 ; = 4Y 2 - Yi = 4 4-1 = 15 Y 4 = 4Y 3 - Y 2 = 5 15-4 = 56, etc., X = 1, 2, 7, 26, 9 9, Y = 0, 1, 4, 15, 5 6, -. PROBLEMS The following problems illustrate the use of the theorems developed h e r e. Without knowledge of these t h e o r e m s, I believe the problems might be difficult to solve. P r o b l e m 1. The numbers 2024 32257 a r e consecutive values of one of the variables in a Pell equation. What a r e the corresponding values of the other v a r i a b l e? (There a r e two solutions.) Problem 2. F o r 8Y 2 = X 2-1 we have X = 1, 3, 17, 9 9, Y = 0, 1, 6, 35, e o e U _ = 6U - U -. n+1 n n - 1 Find another Pell equation(s) for which this r e c u r r e n t law holds. P r o b l e m 3. Prove that v - n Xi Y ± Y - n n " 1 Yi

240 RECURRENT SEQUENCES IN THE PELL EQUATIONS [Oct. v n X * X n * Y t D \ - l Use the + sign if DY* = x\ + 1 the - sign if DYi = X* - 1. Notice that in this problem the recurrent sequence of one variable is developed in terms of constants the other variable. Problem 4. In a Pell equation where D = a 2-1, a > l, prove that X ± X 1 = 0 m.od (Xj i 1) using corresponding signs on each side of the congruence. Problem 5. In Pell equations if DYj = xf + 1, prove: Ex. = X _s±i s i, i + X - Xj - 1 j=l 3=1 c = 2Xi Note that if c = 1 the X? s are Lucas numbers the Y f s are Fibonacci numbers then we have the summation equations for the Lucas Fibonacci 2 2 sequences. If DYj = Xj - 1, show that the comparable summations are E x = ^ I ' V * 1 " 1 j c - 2 ir% Y ^ - Y -Yi L Y 3 = n + 1 c, n 2 \ 0 = 2 X ^ 1. 5 = 1 Problem 6, In each of the following equations find recurrent sequences of rational x f s such that y is integral. The ambiguous sign is used to avoid negative roots.

1969] RECURRENT SERIES IN THE PELL EQUATIONS 241 a) 3x 2 ± 4x + 1 = y 2 b) 3x 2 ± 5x + 2 = y 2 c) 2x 2 ± 6x + 5 = y 2 d) 6x 2 ± 5x + 1 = y 2 EPILOGUE In this part of the paper, some terms notations are introduced which were found to be convenient, a) In the Pell equation, DY 2 = X 2 - (-l) n, n = 0, l, 2, 3 /, e, we notice that as n increases, 1 is alternately subtracted added to the X 2 term. Thus the equation is referred to as an alternating equation. For the equation DY^ = X^ - 1, n = 0, l, 2, 3, - - -, l i s always subtracted from X 2 is referred to as non-alternating. The term alternating Pell equation im- 2 2 plies the minimal equation DYj = X t + 1 the recurrent law U - = cu + U -, whereas the term non-altemating Pell equation implies DYf = x\ - 1 the recurrent law U - = cu - U -. I n this connection it is interesting to note that in recurrent equations where the n f s are negative, the neighboring terms in the sequence developed from U - = U - - cu have opposite signs thus the signs in the sequence alternate. If U - = cu - U - n < 1, the neighboring terms of the sequence have the same signs the sequence is non-alternating. The use of non-positive n f s in the equations of Theorem 1 leads to the conjugate solutions of the Pell equations. b) In the recurrent equation U - = cu +U.,,c > 1 is associated with n+1 n n- 1 - the + sign preceding the U 1 term in the equation U - =.. cu fl - U - c > l is associated with the - sign preceding the U - term. A convenient notation for these recurrent equations is c + c~. For example 6 + implies U n+1 ± 1 = 6 U n +U n-1-4r implies ^ U n+1 _,, = 4U n - U n-1 -. Since c + or c" indicates the manner in which the recurrent sequence is developed they are called the indicator, I, of the sequence. If a /3 are the first two terms of a sequence, then the development of the sequence is completely determined by the indicator the first two terms as l(a,0). For example, if I = 3 +, a = 2, 0 = 3 then 3 +, (2 f 3) d e - fines the sequence implies U - = 3U +U -, U 0 = 2, Uj = 3.

242 RECURRENT SEQUENCES IN THE PELL EQUATIONS Oct. 1969 Throughout my notes I have used this notation because of its convenience brevity. Since each of the Pell equations, (1) (2) have a unique recurrent law for a given D then it follows that they have a unique indicator but a given indicator does not necessarily determine a Pell equation uniquely. c) If a sequence is determined by I, (#,j8) a,j3 have a common factor, f, then all terms of the sequence contain this factor. Let a = fa t js = fft then I, (<*,/3) = I, tfori.fft) = 1, fiflufo). The n term of the sequence can be developed from I, {a i9^) to the n term which is then multiplied by f by this procedure we can use smaller numbers. d) Applying these concepts to the Pell equations we have for the general recurrent solution X = I, Y = I, (1, Xt) Yi(0,l) where I = 2X t if DY? = x\ + 1 I = 2Xj if DYi = X^ - 1. We see that in general for any Pell equation Y s 0 mod Y 4. REFERENCES 1. Edgar Emerson, "On the Integer Solution of 5x 2 ± 6x + 1 = y 2," Fibonacci Quarterly, Vol. 4, No. 1, 1966, pp. 63-69. 2. Burton W. Jones, The Theory of Numbers, Holt, Rinehart Winston, New York, 1961, pp. 82-104. 3. C. D. Olds, Continued Fractions, New Mathematical Library 9, New York, Rom House, Inc., 1963, pp. 61-119. 4. R. Korum, G. McNeil, "A Table of Periodic Continued Fractions from 2-10,000," Lockheed Missiles Space Division, Lockheed Corporation, Sunnyvale, California, pp. I-VII. 5. Robert W. Carmichael, Diophantine Analysis, Dover Publications, Inc., New York, 1959, pp. 26-33.