A NEW METHOD FOR GENERATING MAGIC SQUARES OF SQUARES

RECURSION METHODS TO GENERATE NEW INTEGER SEQUENCES (Part VIF)

Picture of a square

Introduction

The tables of partial imaginary tuples listed in   Part IB,   Part IC,     Part IIIC,   Part ID   and Part IE are useful, when squared, as right diagonals for magic square of squares. I will show that new recursive formulas can be used to generate new integer sequences from the tuple numbers a, b and c using without employing R as part of the multiplier as was done previously in Part VF.

Recursive Formulas for a, b and c Integer Sequences

This section will use the values of a. b and c from line 8 of tables V thru XXIV. It has been found that recursive formulas for the sequences produce linear sequences without resorting to the use of R as part of the multiplier as was shown in Part VF. The indeces x(0) and x(1) are initially given and from these two values the sequences are built up in a recursive fashion. Note that x is a placeholder for a, b or c. The differences (δs) between values of both sets of a and c sequences and the average of the as and cs shown below are identical. Only the a(n) and c(n) values are similar but not quite identical. In addition, since the b values are identical in both tables only one sequence is obtained.

The recursive formulas uses a factor of 6 to multiply to the first number in each formula line. It appears that this factor is obtained from the ratio b(1)/b(0) since the same number pops up in each of the other lines in the tables V/VII or VI/VIII. A second constant ka or kc is added or subtracted to the recursive formulas of the as or the cs sequences but none is added/subtracted to the b sequences. In addition, if the values of either the a or c sequences are averaged, no constant is required and the formulas are structually similar to those of b (see Average Sequence Ia + IIa or Average Sequence Ic + IIc).

The reason the a or c are separated out into separate tables is that they are part of well formed tuples giving the requisite correct square values for the diagonals. If we average out the a or c values, on the other hand, non well formed tuples are generated, i.e., at least one number is a non integer square and, thus, is not suitable for use as diagonals in squares of magic squares. In other words, the equation:

c2 = a2 + 2b2

is not satisfied when the averages of a and/or c are used.

Sequence of a(n) Values

Sequence Ia

From tables V and VII, line 8 the a values of 97 and 433 are initialized followed by generation of a new sequence using the following recursive equations:

a(0) = 97
a(1) = 433
a(2) = a(1)*6 − a(0) − ka
a(3) = a(2)*6 − a(1) − ka
a(4) = a(3)*6 − a(2) − ka
a(5) = a(4)*6 − a(3) − ka
.
.
a(n+1) = a(n)*6 - a(n-1) − ka

where the general formula for the sequences is the last line of the recursive formulas with n ≥ 1 and ka = [(aI(n) − aII(n)]*2. The aI(n) − aII(n) being the difference between two numbers having the same offset taken from each of the sequences of Table V and Table VI, viz., aI(0) − aII(0) = 97 − 47 = 50. In addition, every aI(n) − aII(n) has the same value.

a(2) = 433*6 − (97 + 100) = 2401
a(3) = 2401*6 − (433 + 100) = 13873
a(4) = 13873*6 − (2401 + 100) = 80737
a(5) = 80737*6 − (13873 + 100) = 470449
.
.

The sequence up to 10 values is shown as follows and can be compared to the originals from tables V thru XXIV. Note that the first row corresponds to the differences (δs) between the numbers (a(n)) in the second row with the sequence continuing after the arrow:

δ 336 1968 11472  66864  389712 2271408 
a(n)97 433 2401 13873  80737 470449 2741857
13238736 77161008 449727312 
 15980593 93141601 542868913

Sequence IIa

From tables VI and VIII, line 8 the a values of 47 and 383 are initialized followed by generation of a new sequence :

a(0) = 47
a(1) = 383
a(2) = a(1)*6 − a(0) + ka
a(3) = a(2)*6 − a(1) + ka
a(4) = a(3)*6 − a(2) + ka
a(5) = a(4)*6 − a(3) + ka
.
.
a(n+1) = a(n)*6 − a(n-1) + ka

where the general formula for the sequences is the last line of the recursive formulas with n ≥ 1 and ka = [(aI(n) − aII(n)]*2. The aI(n) − aII(n) being the difference between two numbers having the same offset taken from each of the sequences of Table V and Table VI, viz., aI(0) − aII(0) = 97 − 47 = 50. In addition, every aI(n) − aII(n) has the same value.

a(2) = 383*6 − (47 − 100) = 2351
a(3) = 2351*6 − (383 − 100) = 13823
a(4) = 13823*6 − (2351 − 100) = 80687
a(5) = 80687*6 − (13823 − 100) = 470399
.
.

The sequence up to 10 values is shown as follows and can be compared to the originals from tables V thru XXIV. Note that the first row corresponds to the differences (δs) between the numbers (a(n)) in the second row with the sequence continuing after the arrow:

δ 336 1968 11472  66864  389712 2271408 
a(n)47 383 2351 13823  80687 470399 2741807
13238736 77161008 449727312 
 15980543 93141551 542868863

Sequence Ia + IIa Average

A new sequence is also obtained from the averages of sequence Ia and IIa where the initialization values are 72 and 408 by followed by generation of a new sequence:

a(0) = 72
a(1) = 408
a(2) = a(1)*6 − a(0)
a(3) = a(2)*6 − a(1)
a(4) = a(3)*6 − a(2)
a(5) = a(4)*6 − a(3)
.
.
a(n+1) = a(n)*6 − a(n-1)

where the general formula for the sequences is the last line of the recursive formulas with n ≥ 1.

a(2) = 408*6 − 72 = 2376
a(3) = 2376*6 − 408 = 13848
a(4) = 13848*6 − 2376 = 80712
a(5) = 80712*6 − 13848 = 470424
.
.

The sequence up to 10 values is shown as follows and can be compared to the originals from tables V thru XXIV. Note that the first row corresponds to the differences (δs) between the numbers (a(n)) in the second row with the sequence continuing after the arrow:

δ 336 1968 11472  66864  389712 2271408 
a(n)72 408 2376 13848  80712 470424 2741832
13238736 77161008 449727312 
 15980568 93141576 542868888

Sequence of b(n) Values

Sequence Ib

From tables V and VII, line 8 the b values of 14 and 84 are initialized followed by generation of a new sequence using the following recursive formulas:

b(0) = 14
b(1) = 84
b(2) = b(1)*6 − b(0)
b(3) = b(2)*6 − b(1)
b(4) = b(3)*6 − b(2)
b(5) = b(4)*6 − b(3)
.
.
b(n+1) = b(n)*6 − b(n-1)

where the general formula for the sequences is the last line of the recursive formulas with n ≥ 1.

b(2) = 84*6 − 14 = 490
b(3) = 490*6 − 84 = 2856
b(4) = 2856*6 − 490 = 16646
b(5) = 16646*6 − 2856 = 97020
.
.

The sequence up to 10 values is shown as follows and can be compared to the originals from tables V thru XXIV. Note that the first row corresponds to the differences (δs) between the numbers (b(n)) in the second row with the sequence continuing after the arrow:

δ 70 406 2336  13790 80374 468454 
b(n)14 84 490 2856  16646 97020 565474
2730350 15913646 92751526 
 3295824 19209470 111960996

Sequence of c(n) Values

Sequence Ic

From tables V and VII, line 8 the c values of 99 and 449 are initialized followed by generation of a new sequence using the following recursive formulas:

c(0) = 99
c(1) = 449
c(2) = c(1)*6 − (c(0) + kc)
c(3) = c(2)*6 − (c(1) + kc)
c(4) = c(3)*6 − (c(2) + kc)
c(5) = c(4)*6 − (c(3) + kc)
.
.
c(n+1) = c(n) − c(n-1) − kc

and where the general formula for the sequences is the last line of the recursive formulas with n ≥ 1 and kc = [(cI(n) − cII(n)]*2. The cI(n) − cII(n) being the difference between two numbers having the same offset taken from each of the sequences of Table V and Table VI, viz., cI(0) − cII(0) = 99 − 51 = 48. In addition, every cI(n) − cII(n) has the same value.

c(2) = 449*6 − (99 + 96) = 2499
c(3) = 2499*6 − (449 + 96) = 14449
c(4) = 14449*6 − (2499 + 96) = 84099
c(5) = 84099*6 − (14449 + 96) = 490049
.
.

The sequence up to 10 values is shown as follows and can be compared to the originals from tables V thru XXIV. Note that the first row corresponds to the differences (δs) between the numbers (c(n)) in the second row with the sequence continuing after the arrow:

δ 350 2050 11950  69650  405950 2366050 
c(n)99 449 2499 14449  84099 490049 2856099
13790350 80376050 468465950 
 16646449 97022499 565488449

Sequence IIc

From table VI and VIII, line 8 the c values of 51 and 401 are initialized followed by generation of a new sequence using the following recursive formulas:

c(0) = 51
c(1) = 401
c(2) = c(1)*6 − (c(0) − kc)
c(3) = c(2)*6 − (c(1) − kc)
c(4) = c(3)*6 − (c(2) − kc)
c(5) = c(4)*6 − (c(3) − kc)
.
.
c(n+1) = c(n) − c(n-1) + kc

and where the general formula for the sequences is the last line of the recursive formulas with n ≥ 1 and kc = [(cI(n) − cII(n)]*2. The cI(n) − cII(n) being the difference between two numbers having the same offset taken from each of the sequences of Table V and Table VI, viz., cI(0) − cII(0) = 99 − 51 = 48. In addition, every cI(n) − cII(n) has the same value.

c(2) = 401*6 − (51 − 96) = 2451
c(3) = 2451*6 − (401 − 96) = 14401
c(4) = 14401*6 − (2451 − 96) = 84051
c(5) = 84051*6 − (14401 − 96) = 490001
.
.

The sequence up to 10 values is shown as follows and can be compared to the originals from tables V thru XXIV. Note that the first row corresponds to the differences (δs) between the numbers (c(n)) in the second row with the sequence continuing after the arrow:

δ 350 2050 11950  69650  405950 2366050 
c(n)51 401 2451 14401  84051 490001 2856051
13790350 80376050 468465950 
 16646401 97022451 565488401

Sequence Ic + IIc Average

A new sequence is also obtained from the averages of sequence Ic and IIc where the initialization values are 75 and 425 by followed by generation of a new sequence:

c(0) = 75
c(1) = 425
c(2) = c(1)*6 − c(0))
c(3) = c(2)*6 − c(1)
c(4) = c(3)*6 − c(2)
c(5) = c(4)*6 − c(3)
.
.
c(n+1) = c(n) − c(n-1)

and where the general formula for the sequences is the last line of the recursive formulas with n ≥ 1.

c(2) = 425*6 − 75 = 2475
c(3) = 2475*6 − 425 = 14425
c(4) = 14425*6 − 2475 = 84075
c(5) = 84075*6 − 14425 = 490025
.
.

The sequence up to 10 values is shown as follows and can be compared to the originals from tables V thru XXIV. Note that the first row corresponds to the differences (δs) between the numbers (c(n)) in the second row with the sequence continuing after the arrow:

δ 350 2050 11950  69650  405950 2366050 
c(n)75 425 2475 14425  84075 490025 2856075
13790350 80376050 468465950 
 16646425 97022475 565488425

This concludes Part VIF using two recursive progressions to generate sequences identical to those generated in tables V thru XXIV. Finally to see what started it all go to Part IG which lists the new set of magic square of seven squares.

Go back to homepage.


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