The Pellian equation
In addition, the method of converting a quadratic surd √D into continued fractions (pages 261-262) are also methods that can be used to generate these least solutions. On the other hand, this method is now available on this website as a computer program making it extremely easy to generate the convergents.
As to the topic of this article, it concerns the sequence Sn generated from the the equation 4(n)(n+3) + 5 which was found in the OEIS database with the Sloane number A078371 as the equation (2n + 5)(2n + 1). Note the equations though dissimilar are out of phase by one, i.e., n = 1 in the former is the same as n = 0 in the latter. The comment section describes the numbers in the sequence as quote "generic form of D in the (nontrivially) solvable Pell equation
These D values entered into the Pell calculator program gives initial values of Qn = +4 at n = 3 and 5 which leads to x and y values at positions n = 2 and 4. The Sloane comment that Qn is +4 is confusing and most important it makes no mention that the Qn = +1 for the numbers in this series all fall initially at
Table I shows the first 13 D numbers in the sequence Sn.
D | 21 | 45 | 77 | 117 | 165 | 221 | 285 | 357 | 437 | 525 | 621 | 725 | 837 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
x | 55 | 161 | 351 | 649 | 1079 | 1665 | 2431 | 3401 | 4599 | 6049 | 7775 | 9801 | 12151 |
y | 12 | 24 | 40 | 60 | 84 | 112 | 144 | 180 | 220 | 264 | 312 | 364 | 420 |
The Pell equations are set up according to the following formats with the x and ys taken from Table I. The first equation with n1 and n2 gives the alternate expression for x + y√D where the square term is divided by ½n1 and n2 is a triangular number nt equal to (n2 + n)/2.
Table II gives the formulas for all 13 Pell equations according to the format listed above where both formulas on the right are equivalent. Entry one with
Pell Equation | Equal Expressions |
---|---|
x2 − 5y2 = 1 | R5 = (2 + 1√5)2 = 9 + 4√5 |
x2 − 21y2 = 1 | R21 = (14 + 3√21)2 ∕7 = 55 + 12√21 |
x2 − 45y2 = 1 | R45 = (40 + 6√45)2 ∕20 = 161 + 24√45 |
x2 − 77y2 = 1 | R77 = (88 + 10√77)2 ∕44 = 351 + 40√77 |
x2 − 117y2 = 1 | R117 = (162 + 15√117)2 ∕81 = 649 + 60√117 |
x2 − 165y2 = 1 | R165 = (270 + 21√165)2 ∕135 = 1079 + 84√165 |
x2 − 221y2 = 1 | R221 = (416 + 28√221)2 ∕208 = |
x2 − 285y2 = 1 | R285 = (608 + 36√285)2 ∕304 = 2431 + 144√285 |
x2 − 357y2 = 1 | R357 = (850 + 45√357)2 ∕425 = 3401 + 180√357 |
x2 − 437y2 = 1 | R437 = (1150 + 55√437)2 ∕575 = 4599 + 220√437 |
x2 − 525y2 = 1 | R525 = (1512 + 66√525)2 ∕756 = 6049 + 264√525 |
x2 − 621y2 = 1 | R621 = (1944 + 78√621)2 ∕972 = 7775 + 312√621 |
x2 − 725y2 = 1 | R725 = (2450 + 91√725)2 ∕1225 = 9801 + 364√725 |
x2 − 837y2 = 1 | R837 = (3038 + 105√837)2 ∕1519 = 12151 + 420√837 |
This concludes Part XVII.
To go to Part XVIII where the triangular numbers are used in the calculation of
n1.
Go back to homepage.
Copyright © 2021 by Eddie N Gutierrez. E-Mail: enaguti1949@gmail.com