Sequences of Pythagorean Multiples (Part XVI)

Use of the Diophantine Equation x2 ± D(∑ yj2) = z2 for Generating Sequences

In Part XIII an alternative method was used to generate a couple of partial integral Pythagorean triples, while Part XIV was used to generate all integral Pythagorean triples. In this section the method of Part XIV will be used to generate a sequence of partial integral triples where the z value is a non-perfect square and using the method of Part XIII where n2 is added instead of being subtracted.

We first produce all triples where all the aj are equal and where these generating numbers aj increase consecutively from 1 to k. This is shown in Table I where for example all the ajs are equal to 3. Since the entire equation for z is too big to fit in the cell, it has been set equal to Root and this variable placed within the cell. The last root is some variable, Nk, whose square root is again a non-perfect square.


Root =
 _______________________________ √D2(a12 + a22 + a32 + ... + ak2 )2 + 6n2(a12 + a22 + a32 + ... + ak2) + n4
Table I (D = 1, n= 1)
x y1y2y3...ykz
a1a2a3...ak D(a12 + a22 + a32 + ... + ak2) + n2 2na12na22na3 Root
300(9) + 1 = 10600136
330(9 + 9) + 1 = 19660433
333(9 + 9 + 9) + 1 = 28666892
333...3(9 + 9 + 9 + ...) + 1 = 9k + 1666...6Nk

When the calculations are done for k = 1 to 16, the following results are obtained and listed in Tables IIa and IIb, where the δy2 values are used to avoid non-integral values of y. In addition, the δ1 between adjacent x and z values is 9 which may also be obtained from the arithmetic progression formula using 10 as the initial number and 9 as the common difference:

ak = 10 + 9(k − 1)     (A1)

On inspection the x column shows that these numbers are the identical numbers from the z column of Tables IIa and IIB of Part XIV, while inspection of the y column shows, that only the bold/italized values give y integral values of 6, 12, 18 and 24 and where their common difference is 6. We can extract these particular italized triples from these two tables and generate a sequence of integral Pythagorean triples (Table III). Alternatively, the xs can be obtained by multiplying (aj)2 (in this case 9) with multiples of m2 = 1,4,9,16,... then adding 1 as exhibited in Table I.

Table IIa (δy2 =36)
δ1x y zδ
106136
9
1912433
9
2818892
9
37241513
9
46302296
9
55363241
9
64424348
9
73485617
9
Table IIb (δy2=36)
δ1x y zδ
82547048
9
91608641
9
1006610396
9
1097212313
9
1187814392
9
12784416633
9
1369019036
9
1459621601

On inspection of Table III we see that the common difference between adjacent x values is δ1 and that the difference between adjacent δ1s is δ1δ1 = 18. The common difference between adjacent z2 values, however, is not 18 but is obtained thru four levels of δ differences. In addition, we use the z2 to avoid working with square roots. Moreover, when we get to the end of Table II we can continue generating new x and z2 values by using the various delta (δ) differences.

Table III contains additional triples obtained using δ1 and δ1δ1 columns. In addition, the first row (1,0,1) was added by continuing backwards from the third row.

Table III (δy=6)
δ1δ1δ1x y z2δδδ δδδδδδδ
18101
9135
181061361242
2713772916
183712151341581944
4555354860
188218704890181944
63145536804
181452421601158221944
81303758748
182263051976245701944
995494510692
1832536106921352621944
1179020712636
1844242197128478981944
13513810514580
1857748335233624781944
15320058316524

The x values may be obtained using equation (A2):

bk+1 = bk + (ai + (k − 1)d)    (A2)

where ai is the initial δ1 value, d is the common difference δδ and bk is the variable x. The bk+1, thus, obtained in the first step becomes the bk of the next step. On the other hand, the values for y are obtained using the regular arithmetic progression formula (A1).

ai + (k − 1)d     (A1)

The z values, on the other hand, must be obtained directly from the regular z2 = x2 + y2. By using the delta (δ) differences in the last four columns, the fourth level common difference δδδδ was found to be 1944.

Alternatively, to calculate the values of x using equation (A2), we place the initial numbers for bk of 10 in equation (A2) and perform the calculations for the various bks in Tables IIIx. The bk values are in perfect agreement with those of Table III.

Table IVx (d = δδ = 18)
bk+1 = bk + (ai + (k − 1)d)
b1 = 10 + (27 + 0(18)) = 37
b2 = 37 + (27 + 1(18)) = 82
b3 = 82 + (27 + 2(18)) = 145
b4 = 145 + (27 + 3(18)) = 226
b5 = 226 + (27 + 4(18)) = 325
b6 = 325 + (27 + 5(18)) = 442
b7 = 442 + (27 + 6(18)) = 577
b8 = 577 + (27 + 7(18)) = 730

Alternatively we can calculate the δδδδ values using the light green or key lime green cell values of Table V.

Table V (δy=6)
δ1δ1δ1x y z2δδδ δδδδδδδ
18101
9135
181061361242
2713772916
183712151341581944
4555354860
188218704890181944
63145536804
181452421601158221944
81303758748
182263051976245701944
995494510692
1832536106921352621944
1179020712636
1844242197128478981944
13513810514580
1857748335233624781944
15320058316524
73054535816790021944

We can take any z2row r1, subtractthe cell below, z2row r+1, and then diagonally, subtract then add, subtracting then adding until the δδδδ is reached.

δδδδ = z2row rz2row r+1 + δδδ + δδδ
δδδδ = 1 − 136 + 1377 − 4158 + 4860 = 1944
δδδδ = 1513 − 7048 + 14553 − 15822 + 8748 = 1944

To calculate a z2 starting from a δδδδ we add values diagonally down until the z2 column is reached then dropping down one cell gives a new value. For example, adding up the diagonal light blue elements of the table we eventually reach a new z2 after dropping down one cell:

z2row r+1 = δδδδ + δδδ + δδ + δ + z2row r
z2row r+1 = 1944 + 12636 + 47898 + 138105 + 335233 = 535816

This concludes Part XVI. Go to Part XVII.
Go back to Part XV.
Go back to homepage.


Copyright © 2025 by Eddie N Gutierrez. E-Mail: enaguti1949@gmail.com