Twin Primes Demystified

Twin Prime Distribution Algorithms and Symmetries

Creative Commons License

This work is licensed by its author, Gary William Croft, under a Creative Commons Attribution-ShareAlike 3.0 Unported License.

Return to www.primesdemystified.com Home Page

"This sieve is the source of two efficient prime number factorization algorithms and completely deconstructs the twin prime sequence." – Philip Gervasse Jackson, citing the Prime Spiral Sieve in his Simplicity Instinct: Why Prime Numbers are Elusive!

"It was through my search for meaning that I discovered how important the digital root space is, which in return rewarded me with a look at how elegant and beautiful the world really is." –Talal Ghannam, Ph.D., The Mystery of Numbers Revealed Through Their Digital Root



Twin Prime Distribution Channels

When the set of natural numbers not divisible by 2, 3 or 5 is spiraled within a modulo 30 factorization wheel (see graphic, below), where 1 = 12°, all twin prime candidates (n, n+2), and twin prime pairs (p, p+2), with the exception of (3,5) and (5,7), are arrayed along three distinct "distribution channels," i.e., potential twin prime pairs are separated by 24° and are parsed equally along three sets of paired radii: [12° ↔ 348°] ... [132° ↔ 156°] ... [204° ↔ 228°] ... in intervals of 30:

Twin Prime Candidate Distribution Channels

Below we translate the three twin prime candidate distribution channels into the language of modular and digit sum arithmetic. We define each in modulo 30 terms, and show their angles and digital root sequences, the latter repeating every three rotations around the sieve. Also indicated are the digital root sums of the n, n+2 dyads, which – for reasons we demonstrate below – must always be either 3, 6 or 9 while their digital root products must be 8 in all three cases. [Note: On this site digital root (also known as repeated digit sum, iterated sum-of-digits, or modulus 9 function) refers to the single digit (1-9) that results when all digit sums are summed in turn, e.g., 47 = 4+7 = 11 = 1+1 = 2; the digital root of 47.]:

    Numbers ≡ {11,13} modulo 30 (at 132° & 156°) sequence as:

(11,13) = digital roots {2,4} = digital root 6
(41,43) = digital roots {5,7} = digital root 3
(71,73) = digital roots {8,1} = digital root 9
            {digital roots repeat...}

    Numbers ≡ {17,19} modulo 30 (at 204° & 228°) sequence as:

(17,19) = digital roots {8,1} = digital root 9
(47,49) = digital roots {2,4} = digital root 6
(77,79) = digital roots {5,7} = digital root 3
            {digital roots repeat...}

    Numbers ≡ {29, 1} modulo 30 (at 348° & 12°) sequence as:

(29,31) = digital roots {2,4} = digital root 6
(59,61) = digital roots {5,7} = digital root 3
(89,91) = digital roots {8,1} = digital root 9
            {digital roots repeat...}

You'll note that all of the n, n+2 twin prime candidates listed above sequence their digital roots in one of three variations: {2,4}, {5,7} or {8,1}. We offer a proof why this is so for all twin prime pairs, starting with (11,13), in our Proof by Construction for the Digital Root Sequencing of Twin Primes.

Earlier we drew your attention to the fact that the sums of the digital root dyads must equal 3, 6 or 9 (which, of course, follows from the proof we just noted) because these represent a triad {3,6,9} which – with {1,4,7} and {2,5,8} – are the vertices of rotating equilateral triangles (the "Trinity of Triangles") that together form tensor matrices pointing to {9/3} star polygons (more about that below). We will also reveal the 'magic square' reflecting {3,6,9} hidden in plain sight betwixt and between the twin prime distribution channels.

[Note: In our 'Foundations' section, we stated that the set of natural numbers not divisible by 2, 3 or 5 can be expressed in the (modulo 30) form: 30n+1, 30n+7, 30n+11, 30n+13, 30n+17, 30n+19, 30n+23, 30n+29, once factor-sieved leaving all primes >5. Sub-sorting (to use a spreadsheet term) the prime pair candidates in the twin prime distribution channels can be expressed as (30n+11, 30n+13), (30n+17, 30n+19), (30n+29, 30n+31) ... which narrows the factor-sieving field from 8/30 to 6/30, or 20% of all natural numbers. We've documented this sequence on the Online Encyclopedia of Integer Sequences (OEIS Sequence Number A230462) where, incorporating a suggestion from Omar E. Pol, we've defined it as Numbers congruent to {1, 11, 13, 17, 19, 29} mod 30.]

When we consolidate the 6 radii used to construct the twin prime distribution channels into a single number line, we can completely isolate twin prime candidates, pairing them in tandem in a deterministic sequence, viz.: 11 {+2+4+2+10+2+10} {repeat ... ∞}. For example, six of the first seven such pairings are twin primes: [11, 13] [17, 19] [29, 31] [41, 43] [47, 49 (7 x 7)] [59, 61] [71, 73]. Below is a matrix showing the distribution of twin prime sets (both candidates and confirmed prime pairs) from [11, 13] to [569, 571]:

Twin Prime Distribution Matrix

The first three twin prime sets in this matrix configuration are an anagram, in angular terms (in a 30-sectioned circle such as encompasses the Prime Spiral Sieve, where 1 = 12°), to the first three primes, thus: [11, 13] [17, 19] and [29, 31] respectively translate to:

11 + 13 = 24 ... 24° = 2
17 + 19 = 36 ... 36° = 3
29 + 31 = 60 ... 60° = 5

[You'll note that the sum of the intervals for each row of twin prime candidates = 30 = 2 x 3 x 5; and interestingly the product of the sums of the first three twin prime candidate sets = 24 x 36 x 60 = 51840, while 51840/360 = 122 and 51840/432 = 120.]



Factorization Algorithms for the Twin Prime Candidates

All composite numbers in the twin prime distribution channels can be accounted for by the eight-fold algorithmic progressions shown below, where we first show factorization by multiplication, then follow with factorization by division. [Note: The algorithms shown here are a sub-set of the chordal progressions discussed at some length earlier that factored all composite numbers in the array populated by the set of natural numbers not divisible by 2, 3 and 5, as opposed to only factoring those impacting twin prime candidates, demonstrated below.]:

Twin Prime Candidate Factorization Sequences

Another way to demonstrate twin prime candidate factorization is by sequential division. As you can see below, the left-most column consists of 1 and members of the twin prime distribution channels in sequence, thus delimiting the factorization to only twin prime candidates. Essentially, we've taken the sequenced elements in the twin prime distribution channels (30n+11, 30n+13), (30n+17, 30n+19), (30n+29, 30n+31) and divided them by the set of natural numbers not divisible by 2, 3 or 5: 30n+1, 30n+7, 30n+11, 30n+13, 30n+17, 30n+19, 30n+23, 30n+29. A more concise way to express this is that we divide numbers congruent to {1, 11, 13, 17, 19, 29} mod 30 by numbers congruent to {1, 7, 11, 13, 17, 19, 23, 29} mod 30. All modulo n ≡ 0 results generated (except when a number is divided by itself) indicate factors (arrayed horizontally, hi-lited in pink) of n (arrayed vertically, hi-lited in yellow, except for twin primes hi-lited in bright blue):

Twin Prime Factorization Matrix


The Prime Spiral Sieve as Star Birthing Chamber: Rotating Equilateral Triangles, Tensor Matrices, {9/3} Star Polygons and the Digital Root of Modulo 90 Factorization Sequences in 24 Dimensions

Analyzing the factorization sequences above from the perspective of digital root multiplication, we learn that numbers congruent to {1, 7, 11, 13, 17, 19, 23, 29} modulo 30 can be parsed sequentially by digital root into three groups {2,4}, {8,1} and {5,7} (and, as we discuss at length below and give proof to here, twin primes starting with (11,13) fall into these three categories, as well). We've documented these as twin prime digital root dyad sequences on the Online Encyclopedia of Integer Sequences (and we explore the beautiful symmetries they create under the heading "Digital Root Parsing of the Twin Primes," below). Here are the OEIS sequence numbers and links:

    A232880: Twin primes with digital root 2 or 4

    A232881: Twin primes with digital root 5 or 7

    A232882: Twin primes with digital root 8 or 1

This digital root parsing of our domain leads us to the 6 x 6 = 36 fundamental operations of modulo 30 factorization (see matrix, below). To illustrate the result of this parsing, we show the digital root distribution of the first 36 elements of our domain, including the number 1, and the 29 prime numbers and 6 composite numbers (displaying both their natural and digital root factors) encountered within this range:

Digital Root Dyad Factorization Matrix

Here's another view of the digital root multiplication matrix, this time sorting (horizonally) by digital root products, showing their 6 x 6 = 36 combinations:

The 36 fundamental digital root multiplication dyads.

Can you find this 36 digit palindrome weaving through the matrix above like two opposing square waves?: 1,2,4,5,7,8,4,8,7,2,1,5,7,5,1,8,4,2,2,4,8,1,5,7,5,1,2,7,8,4,8,7,5,4,2,1

When we stretch our modulo 30 / 8-dimension horizon times 3, and look at the above factorization sequences from a modulo 90 / 24-dimension digital root perspective, we make a profound discovery (illustrated in the three graphics that follow): Straight out of group theory comes a triad of rotating equilateral triangles, [{1,4,7}, {2,5,8}, {3,6,9}], interacting to form 3 x 3 tensor matrices (including magic squares) that transform into 24 groups of 24 each {9/3} star polygons – each star's unique sequence of rotations digitally rooted in a matrix that acts as a correspondent between sets of geometric coordinates. These digital root sequences and their associated geometries are infinitely repeating (modulo 90), and the symmetries they create are simultaneously reflectional, rotational and translational. You could say that the Prime Spiral Sieve, when telescoped x 3 to modulo 90, is a star birthing chamber. These symmetry groups are digitally rooted in numbers congruent to {1, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 77, 79, 83, 89} modulo 90, and form the basis of the 24 x 24 "Magic Mirror Matrix" (see also section below titled "Cirque de Primes Juxtaposed with the Magic Mirror Matrix").

These triangles seem reminiscent of–and are perhaps isomorphic to–the quark triads that form protons and neutrons. In both cases we have a fundamental dyad of interacting triads. With regard to our modulo 90 triangulations, {1,4,7} and {2,5,8}, each has six possible rotational symmetry operations referenceable to the values of vertices 1 thru 9 of a {9/3} star polygon; thus we get: {1,4,7}, {1,7,4}, {4,1,7}, {4,7,1}, {7,1,4}, {7,4,1} vs. {2,5,8}, {2,8,5}, {5,2,8}, {5,8,2}, {8,2,5}, {8,5,2}. These triangles interact combinatorially via digital root dyad sums to create {3,6,9}'s six possible operations: {3,6,9}, {3,9,6}, {6,3,9}, {6,9,3}, {9,3,6} and {9,6,3}:

{3,6,9} equilateral triangle's six possible rotations

As we've noted previously, the {3,6,9} triangles provide fundamental structure to our domain in the form of repeating modulo 30, 60, 90 {3,6,9} cycles ... These progressions appear to be isomorphic to the geometric hierarchy of star polygons constructed from equilateral triangles such that (Schäfli notation in braces): {3/1}, or simply {3}, is to modulo 30 as {6/2} is to modulo 60 as {9/3} is to modulo 90. In much the same way that a Buddhist might say "It's the interaction between Yin and Yang that moves fate forward," we're arguing that it's the interaction between {1,4,7} and {2,5,8} that moves structure foward shaped by {3,6,9}'s six combinatorial symmetries. Writing about the importance of group theory developed by Évariste Galois, well known mathematician Marcus du Sautoy captures this idea in his book Symmetry: A Journey into the Patterns of Nature: "[Galois] discovered that it is the interactions between the symmetries in a group that encapsulate the essential qualities of the symmetry of an object."

We stated earlier that this is a 24-dimension model. Not only do each of the 24 modulo 90 factorization sequences have a repetition cycle of 24 (shown vertically on the graphic below), but there is also a hidden dimension generated by successive digital root dyads (starting with n x 7 coupled with n x 11) for each factorization sequence that is also period-24 given that (and there's a proof for this): Any initial dyad of nonzero single-digit numbers enters a periodic cycle of length-24, except for eight pairs that enter cycles of length-8, {(3,3), (3,6), (3,9), (6,3), (6,6), (6,9), (9,3), (9,6)}, and one dyad, (9,9), which is periodic of length-1. (Source: Adapted from a comment posted on the Online Encyclopedia of Integer Sequences: OEIS #A030132: Digital Root of Fibonacci(n)). For example, the first two terms of the Fibonacci digital root sequence just referenced, (1,1), enter a periodic cycle of length-24, as follows: {1, 1, 2, 3, 5, 8, 4, 3, 7, 1, 8, 9, 8, 8, 7, 6, 4, 1, 5, 6, 2, 8, 1, 9} {repeat ... }.

Below are shown the first two of the 24 sets of period-24 digital root factorization sequences i.e., the sequences of numbers congruent to {7} and {11} modulo 90 squentially multiplied times numbers congruent to {1, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 77, 79, 83, 89} modulo 90, starting with 7. (Note: The author is in the process of mapping all 24 similar but unique variations for publication here soon.):

Digital Root Factorization Sequence for numbers congruent to 7 modulo 90 Digital Root Factorization Sequence for numbers congruent to 11 modulo 90

Scanning the factorization products in the two examples above, we can quickly determine which of them impact twin prime candidates, given all twin primes greater than (5,7) are either {11,13}, {17,19} or {1, 29} modulo 30. For example, for the {7} modulo 90 factorization configurations pictured above, we see that exactly 18 of the 24 products (or 75%), index (identify) composite twin prime candidates (shown in brackets): 47 [49], [77], 79, 89 [91], [119] 121, 131 [133], [161] 163, 257 [259], [287] 289, 299 [301], [329] 331, 341 [343], [371] 373, 467 [469], [497] 499, 509 [511], [539] 541, 551 [553], and [581] 583 (Note how the lower vs. higher valued positioning of products alternates with each pair of twin prime candidates.). The remaining 6 products, of course, all render "isolated" prime candidates composite. We go into greater depth regarding digital root dyad cycles specifically related to twin primes under the headings "Digital Root Factorization Dyads of the Twin Prime Distribution Channels" and "Twin Prime Magic Squares," toward the bottom of this page.



Additive Compilation of the Twin Prime Candidates

From what we learned above about segregating twin prime candidates, we can demonstrate that they compile additively in perfect progression, completing an infinite sequence of circles (multiples of 30 and 360):

.
Twin Prime Candidate Additive Compilation

We've multiplied the twin prime candidate pair sums by 15° to show the equivalent arc of rotation of each pair in a modulo 90/24-element per rotation spiral, where 1 = 15° and 24 x 15 = 360°. Thus, when you divide the products by 360°, you get circular equivalents, e.g., 11+13=24; 24x15=360; 360/360=1 (thus, the 1st pair in the first thirt is equivalent to one 360° rotation).



Cirque de Primes

A thirt, in case you're wondering, is a useful unit of measure when discussing intervals in natural numbers not divisible by 2, 3 or 5. A thirt, equivalent to one rotation around the Prime Spiral Sieve is like a mile marker on the prime number highway. If we take the Modulo 30 Prime Spiral Sieve and expand it to Modulo 360, we see that there are 12 thirts in one complete circle, or 'cirque' as we've dubbed it. Each thirt consists of 8 elements. The first two thirts are {1,7,11,13,17,19,23,29} and {31,37,41,43,47,49,53,59}. The last element in thirt1, 29, and the 1st element in thirt2, 31, are rooted in -1 and 1, respectively: where -1+30=29 and 1+30=31. These are twin primes, followed by an infinite sequence in intervals of 30 of twin prime candidates, that serve as the center stripes of the prime number highway, and thus: {-1,+1} +30 {29,31} +30 {59,61} +30 {89,91} +30 {119,121} +30 {149,151} +30 {179,181} +30 {209,211} +30 {239,241} +30 {269,271} +30 {299,301} +30 {329,331} +30 {359,361} completes cirque1, albeit, paradoxically, 361 is the start of cirque2. You can say that {359,361} straddle the 1st and 2nd cirques exactly at 360°, just as {719,721} straddle cirque2 and cirque3 at 720 (or 360x2) ... This is congruent with the fact that all repeating patterns at the digital root level, including the digital roots of Fibonnaci numbers and Lucas numbers indexed by our domain as well as the digital roots of the set constituting our domain, achieve in toto end-to-end symmetrical resolution at 360° intervals, which, of course, is aesthetically satisfying.

Here's a graphic that illustrates these relationships:

Cirque de Primes Illustrating Pattern Resolution at 360°

All prime numbers among the index numbers and the Prime Root Set have been hi-lited in blue. Of the 24 index numbers that are prime, 22 point to elements of the Prime Root Set that are also prime; and the two that aren't, 161 (7x23) & 221 (13x17), are both semi-primes with factors that sum to 30 (and, interestingly, the 30th index # points to 109, a prime number that–like 89–has a decimal expansion based on the Fibonnaci sequence).

You can see from the graphic above that the digital roots of the Fibonacci numbers indexed to our domain (Numbers congruent to {1,7,11,13,17,19,23,29} modulo 30) repeat palindromically every 32 digits (or 4 thirts) consisting of 16 pairs of bilateral 9 sums. The digital root sequence of our domain, on the other hand, repeats every 24 digits (or 3 thirts) and possesses 12 pairs of bilateral 9 sums. The entire Prime Root sequence end-to-end covering 360° has 48 pairs of bilateral 9 sums. And finally, the Prime Root elements themselves within the Cirque, consisting of 96 elements, has 48 pairs of bilateral sums totaling 360. Essentially, the prime number highway consists of infinitely telescoping circles ... Also note, the digital roots of the Prime Root Set as well as the digital roots of Fibonnaci numbers and Lucas numbers (the latter not shown above) indexed to it all sum to 432 (48x9) in 360° cycles. The sequence involving Fibonacci digital roots repeats every 120°, and has been documented by the author on the Online Encyclopedia of Integer Sequences (Sequence #A227896): Digital root of Fibonacci numbers indexed by natural numbers not divisible by 2, 3 or 5 (A007775). (As an aside, we believe the genesis of the Fibonacci sequence is deepened when rooted in -1,1 rather than 1,1, thus creating identity symmetry: -1,1,0,1,1,2,3,5,8 ... Reverse the first two terms (-1 & 1) to generate negative Fibonacci series: 1,-1,0,-1,-1,-2,-3,-5,-8 ...)

[Note: In some circles (no pun intended) the number 432 (as in 432Hz) has been declared the "frequency of the cosmos" to which the "Music of the Spheres" is tuned. A growing number of musicians and musicologists support tuning musical instruments to A=432Hz (or Pythagorean A) rather than the standard concert pitch of A=440Hz. They claim that 432Hz, with its 12-scale octave overtones compared to 440Hz's 8 octave overtones, is harmonically richer and cleaner in sound (The author, who plays American World Music on a flamenco guitar, will attest to this.). It's clear, upon analyzing the cirque de primes, that the digital roots sum to 432 in 360° cycles encompassing 12 thirts. Common sense, if not mathematical rigor, would suggest that there is some correspondence to 432hz's 12-scale octave overtones (and in fact, the Prime Spiral Sieve is featured prominently on this Latin American website's page titled "Fundamental Theorem of Arithmetic Applied Musically to the Zodiacal Circle"). Go to YouTube.com, search on "432," and you'll be impressed by the number of videos featuring this number.]

We get another cirque perspective when we sum the intervals between the twin prime candidate pairs, as matrixed below:

Summing the Intervals Between Twin Prime Candidates

Another way to demonstrate the progression of a Cirque is with the simple function stair-stepping in increments of 152 shown below. If you add 2 to each result you create the twin prime distribution channel rooted in n ≡ 29 (modulo 30) and n ≡ 1 (modulo 30), the 1st prime pair in this sequence being [29,31]. In a profound sense, ignoring the fact that 1 is not a prime number, you could say that the twin prime sequence is effectively grounded and zeroed at {-1,1}:

02 -  12 = -1
152 - 142 = 29
302 - 292 = 59
452 - 442 = 89
  602 - 592 = 119
  752 - 742 = 149
  902 - 892 = 179
  1052 - 1042 = 209
  1202 - 1192 = 239
  1352 - 1342 = 269
  1502 - 1492 = 299
  1652 - 1642 = 329
  1802 - 1792 = 359



Cirque de Primes Juxtaposed with the Magic Mirror Matrix

Now lets zoom in to examine the first 90° partition (encompassing 3 thirts) of our Cirque de Primes, that starts with 1 and ends with 89. By juxtapositioning the 'Magic Mirror Matrix' against this range (see graphic, below), beautiful hidden symmetries are revealed at the digital root level when matrix factorized. This symmetry repeats every 90 degrees, repeatedly incrementing all 24 elements +90 (jumping from 1 thru 89 to 91 thru 179 then +90...n) while the digital roots never change. And thus, this 'geoarithmetric' object, telescoping to infinity, represents the factorization sequencing that accounts for all composite numbers extant in the set of natural numbers not divisible by 2, 3 or 5, leaving all prime numbers greater than 5. For an in depth discussion of the 'Magic Mirror Matrix,' including stunning facts about the number 89 (which, among other things, has the Fibonnaci sequence embedded in the decimal expansion of its reciprocal), click here: Magic Mirror Matrix.


Magic Mirror Matrix Juxapositioned with 3 Thirts of One Cirque


Digital Root Parsing of the Twin Primes

Parsing Twin Primes by their digital roots (in context with natural numbers not divisible by 2, 3 or 5) reveals extremely interesting patterns. Study the matrix below, and you'll discover that it's a beautiful construct of perfectly balanced, multi-directional, 'logic tight' symmetries:

Twin primes parsed by digital roots

Note 1: As indicated, the pink squares marked with x's denote elements with terminating digits = td(5) (their values neither shown nor added into the row sums, above). Although they are not members of our defined set, they are nonetheless interesting. For one thing, the first six td(5) elements sum to 360, as follows: 25 + 35 + 55 + 65 + 85 + 95 = 360. For another, they are always adjacent to modulo 30 "Isolated Primes," discussed below.

Note 2: The first six Isolated Primes, starting with 23, also compile to 360, as follows: 23 + 37 + 53 + 67 + 83 + 97 = 360.

Note 3: Seven is the only number arrayed along the two Isolated Prime radii that isn't 'isolated,' given it is twinned with number 5. Were we to include 7 in our summing of isolated primes, above, the result would be 367, which is itself prime and the sum of seven prime numbers.

Note 4: The first 12 elements that are not td(5), consisting of 11 consecutive prime numbers plus 72, also sum to 360, as follows: 11 + 13 + 17 + 19 + 23 + 29 + 31 + 37 + 41 + 43 + 47 + 49 = 360.

Note 5: Every row possesses "sum symmetry" involving three identical sums. The first row, for example, gives us three sums (working outside-in) totaling 36: 11 + 25; 13 + 23; and 17 + 19 = 36 (and given 36 is 360/10, it comes as no surprise that the 10th row has three sums totaling 360). The 2nd row has three sums totaling 72; the 3rd row has three sums totaling 108, etc. Each successive row thus increments +36, and the row sums accumulate in overlapping cycles of 360 and 432 (the curious relationship between these two numbers discussed elsewhere on this site).



The 900 Block

When we described the Digital Root Parsing Matrix above as 'logic tight,' we were not overstating it. This is evidenced by the fact that, when elements with a terminating digit of 5 are included, every element in the first 49 rows of the matrix can be coupled with (added to) another element in the matrix to total 900 (with the paired elements spatially located in "equal and opposite" positions). [Note: Please keep in mind that it is completely arbitrary of us to fucus on the '900 Block' given there are an infinite number of blocks in an infinite number of incremented configurations within the bounds of the infinitely expandable matrix (e.g., from 47 in Row 3 to 313 in Row 17 forms a '360 Block' and 11 in Row 1 to 1789 in Row 99 (not shown above) would create an '1800 Block').] Obviously, each pair totaling 900 must be comprised of elements producing digital root sums that decompose to dr(9). We've dubbed the 48 dyads consisting of two prime numbers totaling 900 in this block '900 Block Primes' (listed below):

900 block prime number dyads

The Digital Root Parsing Matrix is structured such that any rectangular configuration with 'n' number of columns and 'n' number of rows will produce lateral dyadic sum symmetry, determined by the three dr(9) sum variations: {2,7}, {4,5} and {8,1} (Note: To calculate the bilateral sums for a row, multiply the row's index number times 36.). Here are two examples involving complete (6-element) rows: Take the 17 rows from 2 thru 18 and you'll produce 51 combinations summing to 360; or take the 5 rows from 19 thru 23, which produces 15 combinations summing to 756 (To calculate the number of combinations for any and all rectangles configured from complete rows, multiply the given number of rows times 3.). Now lets carve out a rectangle that doesn't utilize complete rows: we create a 4-column/6-row rectangle with the four corners = 11, 19, 101 and 109 (all prime), and find that all cross-lateral sums equal 120.

digital root parsing matrix with 147 dyads summing to 900

What does this object tell us about twin primes? First, we see that twin primes fall into three categories when parsed by their digital roots: {2,4}, {8,1} and {5,7} (2x4, 8x1 and 5x7 all equal dr(8), while 2+4=dr(6); 8+1=dr(9) and 5+7=dr(3): the three legs of the {3,6,9} triangulation which is a member of the Trinity of Triangles, discussed toward the bottom of this page. Next we learn that every twin prime candidate set wihin the '900 Block' has a twin prime candidate pair 'in opposition' creating two sums of 900. Of those falling into this category within the block above, our favorite consists of four sets of twin primes, thus:

Four sets of twin primes with four dyads intertwined summing to 900


Digital Root Factorization Dyads of the Twin Prime Distribution Channels

As we demonstrated above, employing the logic of digit sum arithmetic, there are a finite number of factorization combinations that correlate to each of the six digital roots, for a total of 6 x 6 = 36 dyads.

The diagram below shows the first 26 twin prime candidates (13 pairs) from the Digital Root Parsing Matrix including the 6 composite numbers falling within this range (from {11,13} to {131,133}). For each composite number we show both its digital root and natural factors.

Twin Prime Digital Root Factorizaton Dyads

Examining the six digital root dyad factorization groups shown above, we discover that the six sets of multiplicands in every case consist of the triangles {1,4,7} and {2,5,8} (and their rotating variations) interlineated, as we show below ({1,4,7} triangles are in bold to help delineate them):

    Digital Root Product = 2: {2, 1, 5, 4, 8, 7}

    Digital Root Product = 4: {4, 2, 1, 8, 7, 5}

    Digital Root Product = 8: {8, 4, 2, 7, 5, 1}

    Digital Root Product = 1: {1, 5, 7, 2, 4, 8}

    Digital Root Product = 5: {5, 7, 8, 1, 2, 4}

    Digital Root Product = 7: {7, 8, 4, 5, 1, 2}

If you think of the above as a numerical block, all rows and columns sum to 27, and working from the outside in (both vertically and horizontally) we find perfect symmetry involving bilateral 9 sums. This digital root factorization object displays the six possible rotational variations for both {1,4,7} and {2,5,8} triangulations and–speaking of symmetry–does so palindromically. Thus:

    {1,4,7}, {4,1,7}, {4,7,1}, {1,7,4}, {7,1,4}, {7,4,1} ... ← palindrome

    {2,5,8}, {2,8,5}, {8,2,5}, {5,2,8}, {5,8,2}, {8,5,2} ... ← palindrome

If we string together the six digital root multiplicand sets to create a 36-digit number, we produce this beautiful palindrome:

    215487421875842751157248578124784512.

The significance of these triangles, as members of rotational symmetry groups that can be configured to create a triad of magic squares, is discussed at length toward the bottom of this page.

As we've discussed previously, every pair of twin prime candidates has one element each from the {1,4,7} and {2,5,8} alternating triangulations (and as indicated, all squares distribute to {1,4,7}–never to {2,5,8}).

Most notably, perhaps, digital root parsing reveals the clustering of twin prime candidates into "six packs," i.e., groupings into three prospective prime pairs demarcated by the isolated prime radii. [Note: The only twin prime set in our domain not in a six-pack {17, 19} sums to 36, which, as 1/10th of 360, seems poetic. If we were to decrement {17, 19} in two -18 steps, the result would be 17 + 19 - 1 + 1 - 17 - 19 = 0, which makes logical sense (from zero to 180 being consistent with the set incrementation, shown below).] The first eleven six packs follow (and you'll find the first three of these in the graphic above). And as we alluded, the sets sequence in +180 increments, thus:

    11+13 + 29+31 + 47+49 = 180
    41+43 + 59+61 + 77+79 = 360
    71+73 + 89+91 + 107+109 = 540
    101+103 + 119+121 + 137+139 = 720
    131+133 + 149+151 + 167+169 = 900
    161+163 + 179+181 + 197+199 = 1080
    191+193 + 209+211 + 227+229 = 1260
    221+223 + 239+241 + 257+259 = 1440
    251+253 + 269+271 + 287+289 = 1620
    281+283 + 299+301 + 317+319 = 1800
    311+313 + 329+331 + 347+349 = 1980

If you think of the above as an 11-row/6-column matrix (excluding the row totals), you'll discover that its 66 elements form 33 dyads that sum to 360 (e.g, 11+349 = 360). Each column stair-steps in increments of 30 and each row, working from the outside in, has three identical sums (e.g., the first row has three dyads that sum to 60 each). Also, we find that the repeating digital root sequences reflect the Trinity of Triangles {3,6,9} {1,4,7} {2,5,8}:

    2+4+2+4+2+4: dr(2) + dr(4) = dr(6)
    5+7+5+7+5+7: dr(5) + dr(7) = dr(3)
    8+1+8+1+8+1: dr(8) + dr(1) = dr(9).

Pointing our discussion of twin prime candidate sequencing triangles back to the Prime Spiral Sieve, we see that three rotations around the Sieve produce six twin prime candidate-related digital root triangles:

    11 (2), 13 (4) ... 17 (8), 19 (1) ... 29 (2), 31 (4)

    41 (5), 43 (7) ... 47 (2), 49 (4) ... 59 (5), 61 (7)

    71 (8), 73 (1) ... 77 (5), 79 (7) ... 89 (8), 91 (1)



Twin Prime Domino Theory

Switching metaphors, below we represent our repeating 'six-pack' clusters as dominoes. The first row of dominoes reveals patterns shared by all three of the 3-pair twin prime candidate configurations, namely the terminating digit and modulo 30 values array identically across the three dyads parsed by digital roots. The second row of dominoes reveals the repeating patterns by digital root and then by digital root with the terminating digit values excluded. [Note: Every category below possesses 'equal and opposite' sum symmetry. For example, the modulo 30 configuration has 9 lateral sums of 30 each.]

Twin Prime Fundamental Patternization

It is evident studying the above that the difference sequence between elements in like positions in the 3-pair clusters is 90. We know from digit sum arithmetic that when you increment +90, both the original elements and the +90 incremented sums will possess the same terminating digits, digital roots, n(modulo 30) results, as well as the same digital roots when terminating digit values are excluded. For example, take the twin prime candidate pairs {89, 91} +90 = {179, 181}:

    89: dr(8); td(9); 29(mod 30); dr(n-td) = 8
    89+90 = 179
    179: dr(8); td(9); 29(mod 30); dr(n-td) = 8

    91: dr(1); td(1); 1(mod 30); dr(n-td) = 9
    91+90 = 181
    181: dr(1); td(1); 1(mod 30); dr(n-td) = 9

The dominoes also tell us that:

The lower-valued (left-hand) member of a prime pair will always possess one of three terminating digit ↔ modulo 30 combinations:
    td(1) ↔ 11(mod 30)
    td(9) ↔ 29(mod 30)
    td(7) ↔ 17(mod 30)

The higher-valued (right-hand) member of a prime pair will always possess one of three terminating digit ↔ modulo 30 combinations:
    td(3) ↔ 13(mod 30)
    td(1) ↔ 1(mod 30)
    td(9) ↔ 19(mod 30)

And lastly, twin prime sets that are: n ≡ 29(mod 30) + 2 = n ≡ 1(mod 30) decompose to: dr(n-td) = 2, 3, 5, 6, 8 or 9, while all other twin primes sets [namely n ≡ 11(mod 30) + 2 = n ≡ 13(mod 30) and n ≡ 17(mod 30) + 2 = n ≡ 19(mod 30)] decompose to: dr(n-td) = 1, 4 or 7.



The Six Digital Root Multiplication Matrices

When we convert the Twin Prime Distribution Channels and/or the paired radii whose angles sum to 360° into digital root multiplication matrices, the results are interesting–to say the least. At the digital root level there are six interrelated matrices of special interest to us: three represent the twin prime distribution channels; the fourth captures the two 'isolated' radii of the Prime Spiral Sieve, i.e. where n ≡ 7 (modulo 30), at 84°, and n ≡ 23 (modulo 30), at 276° (isolated, because, with the exception of 7, which is twinned with 5, the closest proximity they can be to other primes is plus or minus 4, and therefore they can never be twin primes); the fifth and sixth matrices complete the set of four derived from radii whose angles = 360° and in modulo 30 terms sum to 30.

Here's a diagram showing the derivation of the six matrices:

Diagram showing origin of the six digital root matrices

Each of the six matrices draws 24 elements from the first 96 elements of our domain (starting with 1 and ending with 359; you'll find a deep analysis of these 96 elements under the heading 'Cirque de Primes,' above ...). When matrix multiplied, each sequence produces sixteen identical 36-element squares comprised of six 6-element sequences. The matrix derived from the twin prime distribution channel where n ≡ 1 (modulo 30), at 12°, and n ≡ 29 (modulo 30), at 348°, is pictured directly below. It has sequences identical to those generated when you take a digital root multiplication matrix (aka Vedic Square) and remove all 3's, 6's and 9's (This "Imaginary Square," as we dubbed it, is discussed in more detail under prime factorization).

Twin Prime Distribution Channel A_12 and 348 degrees

The matrices for twin prime distribution channels B & C follow:

Twin Prime Distribution Channel B_132 and 156 degrees
Twin Prime Distribution Channel C_204 and 228 degrees

Our next matrix is derived from the Isolated Prime Distribution Radii:

Isolated Prime Distribution Channel D_84 and 276 degrees

And lastly we present the third and fourth matrices derived from paired diagonals whose angles sum to 360°:

Prime Spiral Sieve Radii_132 and 228 degrees
Prime Spiral Sieve Radii_156 and 204 degrees

Each of the digital root multiplication matrices produced by the six channels consists of what are known in mathematics as 'Orthogonal Latin Squares' (defined in Wikipedia as "an n x n array filled with n different symbols, each occurring exactly once in each row and exactly once in each column" ... in our case every row and column of the repeating 6x6 matrices possesses the six elements: 1, 2, 4, 5, 7, 8 in some order). Also, the sum of the multiplicative digital roots = 108 x 24 = 2592 = 432 x 6.

[Note: Channels A, D, E and F combined represent the set of natural numbers not divisible by 2, 3 and 5, the first 24 elements of which form the basis of the Magic Mirror Matrix.]

The graphic below illustrates the transformative relationships between the matrices employing their primary building blocks (one of the sixteen identical 6 x 6 (36 element) Latin Squares that constitute each matrix):

Prime Spiral Sieve Digital Root Multiplication Matrix Transforms

Focusing on just the twin prime distribution channels, we see the relationships shown below [and, directly above, we show that two of the channels (B & C) transform bi-directionally by rotating 180° around one of their principal (lower-left to upper-right) diagonal axes]:

Twin Prime Channel Digital Root Multiplication Relationships


Fibonacci Digital Root Multiplication Matrix and Parsing Symmetries

In our discussion under Cirque de Primes, above, we noted that the digital root of Fibonacci numbers indexed by natural numbers not divisible by 2, 3 or 5 form a 32-beat palindromic sequence that repeats every 120° and sums to 432 every 360°. When matrix multiplied this sequence forms a pyramidal mandala of reflectional symmetries that is vertically, horizontally and diagonally palindromic. If there were room here, we'd expand this matrix times 3 to encompass 360°. Its 96 palindromic rows and 96 palindromic columns would all equal the magic number 432 ... Go figure ...

    Fibonacci Digital Root Multiplication Matrix

When we index the Fibonacci digital roots matrixed above by the twin prime digital root parsing method described earlier, where we parse by digital root dyads {2,4}{1,8} and {5,7}, bilateral sum symmetries are revealed:

Fibonacci Digital Root Parsing of Twin Primes

To complete the picture (the 4th dimension of our 108 x 4 = 432 geometry), we illustrate the Fibonacci digital root sequencing for the 'Isolated Prime' candidates:

Fibonacci Digital Root Parsing of Isolated Primes



Distribution of Squared Prime Numbers via Fibonacci Digital Roots

Analysis of the modulo 120 sequences discussed above reveals that the digital root of Fibonacci numbers indexed by squared members of natural numbers not divisible by 2, 3 or 5 is always 1, while the squared numbers they're indexed by are always congruent to either 1(mod 120) or 49 (mod 120), i.e., The squares of all prime numbers >5 are congruent to {1 or 49} mod 120. This confines the distribution of squared prime numbers >5 to a sequence that is 2/120 or 1.66% of natural numbers. (The author has documented this sequence on the Online Encyclopedia of Integer Sequences [Sequence #A227863], where it's named: Numbers congruent to {1, 49} mod 120.)

This sequence is mysteriously prismatic, as prime numbers embedded in squares are split off in two angles (determined by the terminating digits of their factors), i.e., within a modulo 120 factorization wheel where 1 = 3°, the squares distribute to 3° and 147°, and are thus separated by 144° (And we have to ask whether it's coincidental that the digital roots of the Fibonacci numbers that form the 32-beat repeating palindromic sequence also sum to 144, not to mention, as we've already noted above, when the first three twin prime candidate pairs are summed and multiplied together, 24 x 36 x 60 = 51840, and 51840/360 = 144. This last fact will take on even greater significance when we diagram the 'square roots of squares,' second image, below. One more riff on 144: it's the 12th Fibonacci number following 89 and the only Fibonacci number that's a square other than Fibo1. When you divide 144 by 89, a number with magical properties we explore under 'Fibonacci Twins,' below, you produce an approximation of the golden ratio, or phi, with a difference of only +0.000056460660007.). Another fascinating feature of this array is that any even number of–not necessarily contiguous–factors drawn from any one of the 32 angles in this modulo 120 configuration distribute products to 1(mod 120) or 49 (mod 120), along with the squares. Two examples: 1) at 33°: 11 x 131 x 251 x 371 = 134187361, which is congruent to 1 (mod 120) ; ... 2) at 291°: 97 x 217 = 21049, which is congruent to 49 (mod 120). The matrix below illustrates these relationships:

Squares Congruent to 1 and 49 mod 120

Referencing the matrix above, the illustration below sieves the squares congruent to {1,49} mod 120 and takes their square roots. Twin prime candidates are emphasized to show their repeating zig-zag sequencing on the rungs of a perfectly counter-balanced lattice structure mysteriously populated by the set of natural numbers not divisible by 2, 3 or 5:

Twin Prime Candidates Squared Congruent to 1 and 49 mod 120



Fibonacci Convergence

It is well known that the ratio of any Fibo(n+2)/Fibo(n) will converge to a limit φ + 1 = φ2 (an irrational number) as n approaches infinity (φ being the symbol for phi or the "golden ratio" which = 1.6180339887498948482045868343656 ...). And so it follows that the ratios of the Fibonacci numbers indexed to the twin primes (n and n+2) and/or twin prime candidates in sequence converge accordingly (in other words, the square roots of the ratios of Fibonacci numbers indexed to the twin primes and/or twin prime candidates in sequence converge to φ).

Here are some examples (square roots not taken):

Fibo(13)/Fibo(11) = 2.617977528089887640449438202247 ...
Fibo(19)/Fibo(17) = 2.618033813400125234815278647464 ...
Fibo(31)/Fibo(29) = 2.618033988748203621343798191078 ...
Fibo(43)/Fibo(41) = 2.618033988749894831892914017992 ...
Fibo(49)/Fibo(47) = 2.618033988749894848153928976786 ...
Fibo(61)/Fibo(59) = 2.618033988749894848204586345776 ...

Here's an interesting connection to Fibo(13)/Fibo(11) = 2.617977528089887640449438202247: Fibo(13) = 233, while Fibo(11) = 89. The square root of Fibo(13) divided by Fibo(11) = 122 divided by 89 (and 89 has the basis of the Fibonacci sequence in the decimal expansion of its reciprocal). Then consider that 3 x 122 = 432, a fascinating number we discuss above under 'Cirque de Primes,' below under 'Fibonacci Twins,' as well as on our page detailing the Magic Mirror Matrix where we also examine the amazing properties of number 89.



Perfect Twins

One other twin prime related deterministic algorithm worth noting stair-steps its way up the twin prime distribution channels with a periodicity of 6. We've dubbed these "perfect twins," given that the square root of their 'perfect' square sums is the first perfect number, 6, or one of its multiples. To calculate this interesting twin pair candidate sequence, start with x = 6 (then add 6 and repeat for each successive step, as shown below):

  _________________
 (x2/2 - 1 + x2/2 +1) 

x=6: 62 = 36; 36/2 = 18; 18-1 = 17; 18+1 = 19; thus, twin pair candidates = 17 and 19; take the square root of 17+19 which = 6

add 6 (6+6 = 12)

x=12: 122 = 144; 144/2 = 72; 72-1 = 71; 72+1 = 73; thus, twin pair candidates = 71 and 73; take the square root of 71+73 which = 12

add 6 (12+6 = 18) {repeat ... n}

The table below shows the first 20 steps in the sequence. All "perfect twins" are hi-lited in gold. [Note: The table is arrayed in five columns to show the +30 and repeating modulo 30 vertical incrementation. Also note that the digital roots for all prime pair candidates and/or prime pairs in this sequence equate to 8 + 1 = 9, e.g., 17, or dr(8) + 19, or dr(1) = 36 or dr(9)].

Perfect Twin Primes

[Note: We've documented these "perfect twins" on the Online Encyclopedia of Integer Sequences (Sequence #A232878) where they're defined as Twin prime pairs which sum to perfect squares.]



Infinite Twins

Alluding to the Twin Prime Conjecture, given that all composite numbers extant in the infinite set of natural numbers not divisible by 2, 3 or 5 (which by definition contains all prime numbers ≥ 7 and their multiples and is defined as 1 {+6 +4 +2 +4 +2 +4 +6 +2} {repeat ... ∞}), can be produced algorithmically, employing eight geometrically-expanding factorization sequences, ultimately leaving the equal distribution of all prime numbers ≥ 7 ... n along eight diagonals; and given that six of these diagonals form the three twin prime distribution channels (described above), which, when combined into a single number line, reveal a divergent series: 11 {+2 +4 +2 +10 +2 +10} {repeat ... ∞}; and given that the infinity of prime numbers is proven, it's axiomatic that there is an infinity of prime pairs in the form n, n + 2. Put differently, twin primes, starting with [11, 13] are sequenced infinitely within a divergent (aka harmonic) series while paradoxically converging to φ2 when indexed to Fibonacci numbers like meshing gears.



Seeing Stars

Another interesting aspect of the three twin prime distribution channels becomes apparent when the eight radials of the Prime Spiral Sieve (or modulo 30 wheel factorization), are superimposed upon a {10/3} regular star polygon (see illustration, below). Thirty (30) is the total unit length of the ten straight line segments used to construct the {10/3} star polygon, where the straight lines connect every 3rd point of 10 equally spaced points lying on a circle's circumference. Note that three points of the polygon, 144°(12), 216°(18) and 360°(30), precisely and symmetrically split the twin prime distribution channels. And as you rotate this object, several other symmetries become apparent, e.g.: between 7 → 17; 13 → 23; and 7 → 23.

Modulo 30 Prime Spiral Sieve and Star Polygon Superimposed

Even more interesting, is the result when you superimpose a 15-point star over a Prime Spiral Sieve with its 'aperture' increased to modulo 90. As shown below, every one of the 9 twin prime distribution channels is hit dead center by points of the star:

Modulo 90 Prime Spiral Sieve and 15-Point Star Superimposed

This superimposition reveals some interesting circularities: When you sum the 9 angles that hit the twin prime distribution channels dead center you get 1800. 1800/360=5; 1800/9=200. In turn, sum the 6 angles that hit equidistant between the twin prime distribution channels and you get 1080. 1080/360=3. 1080/6=180. In toto, 1800+1080=2880; 2880/360=8.



Twin Prime Magic Squares

Before leaving the twin primes, let's take one more look at them from a digital root perspective. The illustration below shows the infinitely repeating digital root sequences for each of the 8 radii of the spiral sieve, as well the sums of the digital roots (in yellow) between each radius, including the twin prime distribution channels. It's notable that the twin pair digital sum sequences, {3,9,6}, {9,6,3} and {6,3,9} when positioned adjacent to each other in ascending order, create a "magic square" (in that all vertical, horizontal and principal diagonal sums equal 18), as pictured below:

Twin Prime Digital Root Magic Square

The {3,6,9} triangulation is a member of the so-called 'Trinity of Triangles,' along with {1,4,7} and {2,5,8}. Looking at the graphic above, you'll note that each of the three Twin Prime Distribution Channels has one leg that sequences as {1,4,7} and the other as {2,5,8}, albeit the starting digits may vary, viz. {1,4,7} vs. {4,7,1}. When you combine the three legs that sequence as {1,4,7} as well as the three legs that sequence as {2,5,8}, these also form magic squares, as illustrated below. Another magic square hidden in plain sight becomes evident when you examine the sequences you encounter in rotation: {1,4,7} ... {2,5,8} ... {3,6,9} in various sequence combinations. Taking the digits 1 thru 9 these three triangles represent, they can be configured into a magic square representing 'et alia', in this case the famous 'Lo Shu Magic Square':

Trinity of Triangles and Trinity of Magic Squares

Under "Factorization Algorithms for the Twin Prime Candidates," above, we detailed how twin primes come in three digital root configurations: {2,4}, {5,7} and {8,1} and demonstrated that these dyads, in turn, can form the basis of period-24 digital root cycles that generate tensor matrices. So with these cycles and matrices in mind, take the three twin prime digital root dyads under discussion, sequence their 24-length digital root cycles, then tier them, and you'll discover a direct means to constructing the Trinity of Magic Squares via the Trinity of Triangles, as illustrated here:

Magic Squares Derived from Twin Prime Digital Root Dyad Cycles

When you rotate either the {1,4,7} or {2,5,8} magic square around its horizontal axis, i.e. columns {A,B,C} become {C,B,A}, then add the {1,4,7} {2,5,8} magic squares together, you produce a square with nine 9's. For example, adding the first rows of each gives us: {2,8,5} + {7,1,4} = {9,9,9}.

Triangles and magic squares similar–or identical–to those shown above can be derived from the digital root sequence cycles of all three twin prime distribution channels (namely numbers congruent to {11,13}, {17,19} and {1,29} modulo 30). This is also true of dyads formed by paired radii of the Prime Spiral Sieve that sum to 30, i.e., numbers congruent to {1,29}, {7,23}, {11,19}, or {13,17} modulo 30, as well as dyads formed when {n, n + 10} are congruent to {1, 11}, {7, 17}, {13, 23} or {19, 29} modulo 30 (note their pairing by terminating digits). One example relating to twin primes: The first three candidate pairs in the twin prime distribution channel congruent to {11,13} modulo 30 (all three of which are indeed twin primes) sequence their digital roots as follows:

   {11,13} = digital roots 2 & 4
   {41,43} = digital roots 5 & 7
   {71,73} = digital roots 8 & 1.

As you can see, this is the same digital root sequence illustrated above.

It appears that the triangulations and magic squares structuring the distribution of twin primes (and as it turns out, all prime numbers) have a genesis in universal principles involving symmetry groups rotated by the 8-dimensional algorithms discussed at length on this site. You can see this universal principle at work, for example, with regard to the Fibonacci digital root sequence when coupled to a pair of dyads that follow certain incremental rules. As we illustrated above, the initializing dyad of the period-24 Fibonacci digital root sequence is {1,1, ...}. We can generate triangles and magic squares by tiering the Fibonacci digital root sequence with two pairs of terms that are + 3 or + 6 from the initial terms {1,1}. The values of the 2nd and 3rd tiers, or rows, must differ, or symmetry is lost. In other words, the first two columns should read either {1,4,7 + 1,7,4, or vice versa} but not {1,4,7 + 1,4,7, or 1,7,4, + 1,7,4}. Here's a picture that displays the {1,4,7 + 1,7,4} digital root sequencing used to construct triangles and magic squares:

Magic Squares Derived from Fibonacci digital root dyad cycle

We'll end our exploration of twin primes with two powerful magic squares. First: The Trinity of Magic Squares x 3 can be constructed into a magic square where all columns, rows and principal diagonals equal 45 (and were you to construct an 8x8 square employing 64 Trinity Magic Squares shown below, all columns, rows and principal diagonals would equal 360).

Trinity Magic Square

Next we show how the Lo Shu Magic Square and Trinity of Magic Squares can be combined in a "4-Fold Way" to create their own magic square when configured in a 4 x 4 matrix of magic squares consisting of 144 elements where all rows, columns and principal diagonals sum to 60 and the sum total of the matrix's elements = 720 = 360 x 2 (and were we to construct a 6x6 square employing 36 4-Fold Way Magic Squares shown below, all columns, rows and principal diagonals would equal 360).

Lo Shu Magic Square with Trinity of Magic Squares create Fourfold Way Magic Square

Return to www.primesdemystified.com Home Page

Your feedback welcome! Email: gwc@hemiboso.com