Magic square


In recreational mathematics and combinatorial design, a magic square is a square grid filled with distinct positive integers in the range such that each cell contains a different integer and the sum of the integers in each row, column and diagonal is equal. The sum is called the magic constant or magic sum of the magic square. A square grid with cells on each side is said to have order n.
The mathematical study of magic squares typically deals with its construction, classification, and enumeration. Although completely general methods for producing all the magic squares of all orders do not exist, historically three general techniques have been discovered: by bordering method, by making composite magic squares, and by adding two preliminary squares. There are also more specific strategies like the continuous enumeration method that reproduces specific patterns. The magic squares are generally classified according to their order n as: odd if n is odd, evenly even if n = 4k, oddly even if n = 4k + 2. This classification is based on different techniques required to construct odd, evenly even, and oddly even squares. Beside this, depending on further properties, magic squares are also classified as associative magic squares, pandiagonal magic squares, most-perfect magic squares, and so on. More challengingly, attempts have also been made to classify all the magic squares of a given order as transformations of a smaller set of squares. Except for n ≤ 5, the enumeration of higher order magic squares is still an open challenge. The enumeration of most-perfect magic squares of any order was only accomplished in the late 20th century.
In regard to magic sum, the problem of magic squares only requires the sum of each row, column and diagonal to be equal, it does not require the sum to be a particular value. Thus, although magic squares may contain negative integers, they are just variations by adding or multiplying a negative number to every positive integer in the original square. Magic squares composed of integers are also called normal magic squares, in the sense that there are non-normal magic squares whose integers are not restricted to. However, in some places, "magic squares" is used as a general term to cover both the normal and non-normal ones, especially when non-normal ones are under discussion. Moreover, the term "magic squares" is sometimes also used to refer to various types of word squares.
Magic squares have a long history, dating back to at least 190 BCE in China. At various times they have acquired magical or mythical significance, and have appeared as symbols in works of art. In modern times they have been generalized a number of ways, including using extra or different constraints, multiplying instead of adding cells, using alternate shapes or more than two dimensions, and replacing numbers with shapes and addition with geometric operations.
's Melancholia I includes an order 4 square with magic sum 34

History

The third order magic square was known to Chinese mathematicians as early as 190 BCE, and explicitly given by the first century of the common era. The first datable instance of the fourth-order magic square occur in 587 CE in India. Specimens of magic squares of order 3 to 9 appear in an encyclopedia from Baghdad, the Encyclopedia of the Brethren of Purity. By the end of 12th century, the general methods for constructing magic squares were well established. Around this time, some of these squares were increasingly used in conjunction with magic letters, as in Shams Al-ma'arif, for occult purposes. In India, all the fourth-order pandiagonal magic squares were enumerated by Narayana in 1356. Magic squares were made known to Europe through translation of Arabic sources as occult objects during the Renaissance, and the general theory had to be re-discovered independent of prior developments in China, India, and Middle East. Also notable are the ancient cultures with a tradition of mathematics and numerology that did not discover the magic squares: Greeks, Babylonians, Egyptians, and Pre-Columbian Americans.

China and Japan

While ancient references to the pattern of even and odd numbers in the 3×3 magic square appears in the I Ching, the first unequivocal instance of this magic square appears in the chapter called Mingtang of a 1st-century book Da Dai Liji, which purported to describe ancient Chinese rites of the Zhou dynasty.
These numbers also occur in a possibly earlier mathematical text called Shushu jiyi, said to be written in 190 BCE. This is the earliest appearance of a magic square on record; and it was mainly used for divination and astrology. The 3×3 magic square was referred to as the "Nine Halls" by earlier Chinese mathematicians. The identification of the 3×3 magic square to the legendary Luoshu chart was only made in the 12th century, after which it was referred to as the Luoshu square. The oldest surviving Chinese treatise on the systematic methods for constructing larger magic squares is Yang Hui's Xugu zheqi suanfa written in 1275. The contents of Yang Hui's treatise were collected from older works, both native and foreign; and he only explains the construction of third and fourth-order magic squares, while merely passing on the finished diagrams of larger squares. He gives a magic square of order 3, two squares for each order of 4 to 8, one of order nine, and one semi-magic square of order 10. He also gives six magic circles of varying complexity.
492
357
816

216133
115810
79126
144115

12316421
151471811
24171392
20819126
53102225

132218271120
314369292
122114231625
303532347
172610191524
835281633

468162029749
340353618412
4412332319386
28261125392422
5373127171345
4891514321047
143343021424

613264577660
125455916505113
2046471724424321
3727264033313036
2935343225393828
4422234148181945
5214154956101153
559588163624

317613368118297411
224058274563203856
674497295465247
307512327714347916
213957234159254361
663486855070752
358017287310337815
264462193755244260
718536414669651

The above magic squares of orders 3 to 9 are taken from Yang Hui's treatise, in which the Luo Shu principle is clearly evident. The order 5 square is a bordered magic square, with central 3×3 square formed according to Luo Shu principle. The order 9 square is a composite magic square, in which the nine 3×3 sub squares are also magic. After Yang Hui, magic squares frequently occur in Chinese mathematics such as in Ding Yidong's Dayan suoyin, Cheng Dawei's Suanfa tongzong, Fang Zhongtong's Shuduyan which contains magic circles, cubes and spheres, Zhang Chao's Xinzhai zazu, who published China's first magic square of order ten, and lastly Bao Qishou's Binaishanfang ji, who gave various three dimensional magic configurations. However, despite being the first to discover the magic squares and getting a head start by several centuries, the Chinese development of the magic squares are much inferior compared to the Indian, Middle Eastern, or European developments. The high point of Chinese mathematics that deals with the magic squares seems to be contained in the work of Yang Hui; but even as a collection of older methods, this work is much more primitive, lacking general methods for constructing magic squares of any order, compared to a similar collection written around the same time by the Byzantine scholar Manuel Moschopoulos. This is possibly because of the Chinese scholars' enthralment with the Lo Shu principle, which they tried to adapt to solve higher squares; and after Yang Hui and the fall of Yuan dynasty, their systematic purging of the foreign influences in Chinese mathematics.
The Japanese interest in magic square began after the dissemination of Yang Hui's Suanfa and Cheng Dawei's Suanfa tongzong in the 17th century, and almost all the wasans devoted their time to its study. In the 1660 edition of Ketsugi-sho, Isomura Kittoku gave both odd and even ordered bordered magic squares as well as magic circles; while the 1684 edition of the same book contained a large section on magic squares, demonstrating that he had a general method for constructing bordered magic squares. In Jinko-ki by Muramatsu Kudayu Mosei, both magic squares and magic circles are displayed. The largest square Mosei constructs is of 19th order. Various magic squares and magic circles were also published by Nozawa Teicho in Dokai-sho, Sato Seiko in Kongenki, and Hosino Sanenobu in Ko-ko-gen Sho. One of Seki Takakazu's Seven Books is devoted completely to magic squares and circles. This is the first Japanese book to give a general treatment of magic squares in which the algorithms for constructing odd, singly even and doubly even bordered magic squares are clearly described. In 1694 and 1695, Yueki Ando gave different methods to create the magic squares and displayed squares of order 3 to 30. A fourth-order magic cube was constructed by Yoshizane Tanaka in Rakusho-kikan. The study of magic squares was continued by Seki's pupils, notably by Katahiro Takebe, whose squares were displayed in the fourth volume of Ichigen Kappo by Shukei Irie, Yoshisuke Matsunaga in Hojin-Shin-jutsu, Yoshihiro Kurushima in Kyushi Iko who rediscovered a method to produce the odd squares given by Agrippa, and Naonobu Ajima. Thus by the beginning of the 18th century, the Japanese mathematicians were in possession of methods to construct magic squares of arbitrary order. After this, attempts at enumerating the magic squares was initiated by Yamaji.

India

The 3×3 magic square first appears in India in Gargasamhita by Garga, who recommends its use to pacify the nine planets. The oldest version of this text dates from 100 CE; however passage on planets could not have been written earlier than 400 CE. The first datable instance of 3×3 magic square in India occur in a medical text Siddhayog by Vrnda, which was prescribed to women in labor in order to have easy delivery.
The oldest datable fourth order magic square in the world is found in an encyclopaedic work written by Varahamihira around 587 CE called Brhat Samhita. The magic square is constructed for the purpose of making perfumes using 4 substances selected from 16 different substances. Each cell of the square represents a particular ingredient, while the number in the cell represents the proportion of the associated ingredient, such that the mixture of any four combination of ingredients along the columns, rows, diagonals, and so on, gives the total volume of the mixture to be 18. Although the book is mostly about divination, the magic square is given as a matter of combinatorial design, and no magical properties are attributed to it.
2358
5823
4176
7641

103138
516211
49714
156121

The square of Varahamihira as given above has sum of 18. Here the numbers 1 to 8 appear twice in the square. It is a pan-diagonal magic square. It is also an instance of most perfect magic square. Four different magic squares can be obtained by adding 8 to one of the two sets of 1 to 8 sequence. The sequence is selected such that the number 8 is added exactly twice in each row, each column and each of the main diagonals. One of the possible magic squares shown in the right side. This magic square is remarkable in that it is a 90 degree rotation of a magic square that appears in the 13th century Islamic world as one of the most popular magic squares.
The construction of 4th-order magic square is detailed in a work titled Kaksaputa, composed by the alchemist Nagarjuna around 10th century CE. All of the squares given by Nagarjuna are 4×4 magic squares, and one of them is called Nagarjuniya after him. Nagarjuna gave a method of constructing 4×4 magic square using a primary skeleton square, given an odd or even magic sum. Incidentally, the special Nagarjuniya square cannot be constructed from the method he expounds. The Nagarjuniya square is given below, and has the sum total of 100.
30161836
10442224
32142034
2826406

7146
2853
5328
4671

The Nagarjuniya square is a pan-diagonal magic square. The Nagarjuniya square is made up of two arithmetic progressions starting from 6 and 16 with eight terms each, with a common difference between successive terms as 4. When these two progressions are reduced to the normal progression of 1 to 8, we obtain the adjacent square.
Around 12th-century, a 4×4 magic square was inscribed on the wall of Parshvanath temple in Khajuraho, India. Several Jain hymns teach how to make magic squares, although they are undatable.
As far as is known, the first systematic study of magic squares in India was conducted by Thakkar Pheru, a Jain scholar, in his Ganitasara Kaumudi. This work contains a small section on magic squares which consists of nine verses. Here he gives a square of order four, and alludes to its rearrangement; classifies magic squares into three according to its order; gives a square of order six; and prescribes one method each for constructing even and odd squares. For the even squares, Pheru divides the square into component squares of order four, and puts the numbers into cells according to the pattern of a standard square of order four. For odd squares, Pheru gives the method using horse move or knight's move. Although algorithmically different, it gives the same square as the De la Loubere's method.
The next comprehensive work on magic squares was taken up by Narayana Pandit, who in the fourteenth chapter of his Ganita Kaumudi gives general methods for their construction, along with the principles governing such constructions. It consists of 55 verses for rules and 17 verses for examples. Narayana gives a method to construct all the pan-magic squares of fourth order using knight's move; enumerates the number of pan-diagonal magic squares of order four, 384, including every variation made by rotation and reflection; three general methods for squares having any order and constant sum when a standard square of the same order is known; two methods each for constructing evenly even, oddly even, and odd squares when the sum is given. While Narayana describes one older method for each species of square, he claims the method of superposition for evenly even and odd squares and a method of interchange for oddly even squares to be his own invention. The superposition method was later re-discovered by De la Hire in Europe. In the last section, he conceives of other figures, such as circles, rectangles, and hexagons, in which the numbers may be arranged to possess properties similar to those of magic squares. Below are some of the magic squares constructed by Narayana:
816
357
492

114415
811510
132163
12796

161473023
241710831
322518114
528261912
136292220

135433326
2511928830
241418161722
132319212015
12262710297
362343531

3526171625344
4637211236455
5741322314566
6152433425167
2635445362711
1346556473122
2415667514233

605344374132029
314193059544338
585542392151831
116173257564140
615245365122128
611222762514635
635047347102326
89242564494833

The order 8 square is interesting in itself since it is an instance of the most-perfect magic square. Incidentally, Narayana states that the purpose of studying magic squares is to construct yantra, to destroy the ego of bad mathematicians, and for the pleasure of good mathematicians. The subject of magic squares is referred to as bhadraganita and Narayana states that it was first taught to men by god Shiva.

Middle East: Persia, Arabia, North Africa, Muslim Iberia

Although the early history of magic squares in Persia and Arabia is not known, it has been suggested that they were known in pre-Islamic times. It is clear, however, that the study of magic squares was common in medieval Islam, and it was thought to have begun after the introduction of chess into the region. The first datable appearance of magic square of order 3 occur in Jābir ibn Hayyān's Kitab al-mawazin al-Saghir where the magic square and its related numerology is associated with alchemy. While it is known that treatises on magic squares were written in the 9th century, the earliest extant treaties we have date from the 10th-century: one by Abu'l-Wafa al-Buzjani and another by Ali b. Ahmad al-Antaki. These early treatise were purely mathematical, and the Arabic designation for magic squares used is wafq al-a'dad, which translates as harmonious disposition of the numbers. By the end of 10th century, the two treatises by Buzjani and Antaki makes it clear that the Middle Eastern mathematicians had understood how to construct bordered squares of any order as well as simple magic squares of small orders which were used to make composite magic squares. A specimen of magic squares of orders 3 to 9 devised by Middle Eastern mathematicians appear in an encyclopedia from Baghdad, the Rasa'il Ikhwan al-Safa. The squares of order 3 to 7 from Rasa'il are given below:
276
951
438

414151
97612
511108
162313

21341225
15176198
102413216
18720911
11422235

11223252318
25167301320
276353643
1031123334
14198292615
24172891221

47118964549
4372017163546
2182621283248
4319272523317
38362229241412
40153033341310
13942414453

The 11th century saw the finding of several ways to construct simple magic squares for odd and evenly-even orders; the more difficult case of evenly-odd case was solved by Ibn al-Haytham with k even, and completely by the beginning of 12th century, if not already in the latter half of the 11th century. Around the same time, pandiagonal squares were being constructed. Treaties on magic squares were numerous in the 11th and 12th century. These later developments tended to be improvements on or simplifications of existing methods. From the 13th century on wards, magic squares were increasingly put to occult purposes. However, much of these later texts written for occult purposes merely depict certain magic squares and mention their attributes, without describing their principle of construction, with only some authors keeping the general theory alive. One such occultist was the Egyptian Ahmad al-Buni, who gave general methods on constructing bordered magic squares; some others were the 17th century Egyptian Shabramallisi and the 18th century Nigerian al-Kishnawi.
The magic square of order three was described as a child-bearing charm since its first literary appearances in the alchemical works of Jābir ibn Hayyān and al-Ghazālī and it was preserved in the tradition of the planetary tables. The earliest occurrence of the association of seven magic squares to the virtues of the seven heavenly bodies appear in Andalusian scholar Ibn Zarkali's Kitāb tadbīrāt al-kawākib. A century later, the Egyptian scholar Ahmad al-Buni attributed mystical properties to magic squares in his highly influential book Shams al-Ma'arif, which also describes their construction. This tradition about a series of magic squares from order three to nine, which are associated with the seven planets, survives in Greek, Arabic, and Latin versions. There are also references to the use of magic squares in astrological calculations, a practice that seems to have originated with the Arabs.

Latin Europe

Unlike in Persia and Arabia, we have better documentation of how the magic squares were transmitted to Europe. Around 1315, influenced by Arab sources, the Greek Byzantine scholar Manuel Moschopoulos wrote a mathematical treatise on the subject of magic squares, leaving out the mysticism of his Middle Eastern predecessors, where he gave two methods for odd squares and two methods for evenly even squares. Moschopoulos was essentially unknown to the Latin Europe until the late 17th century, when Philippe de la Hire rediscovered his treatise in the Royal Library of Paris. However, he was not the first European to have written on magic squares; and the magic squares were disseminated to rest of Europe through Spain and Italy as occult objects. The early occult treaties that displayed the squares did not describe how they were constructed. Thus the entire theory had to be rediscovered.
Magic squares had first appeared in Europe in Kitāb tadbīrāt al-kawākib written by Ibn Zarkali of Toledo, Al-Andalus, as planetary squares by 11th century. The magic square of three was discussed in numerological manner in early 12th century by Jewish scholar Abraham ibn Ezra of Toledo, which influenced later Kabbalists. Ibn Zarkali's work was translated as Libro de Astromagia in the 1280s, due to Alfonso X of Castille. In the Alfonsine text, magic squares of different orders are assigned to the respective planets, as in the Islamic literature; unfortunately, of all the squares discussed, the Mars magic square of order five is the only square exhibited in the manuscript.
Magic squares surface again in Florence, Italy in the 14th century. A 6×6 and a 9×9 square are exhibited in a manuscript of the Trattato d'Abbaco by Paolo Dagomari. It is interesting to observe that Paolo Dagomari, like Pacioli after him, refers to the squares as a useful basis for inventing mathematical questions and games, and does not mention any magical use. Incidentally, though, he also refers to them as being respectively the Sun's and the Moon's squares, and mentions that they enter astrological calculations that are not better specified. As said, the same point of view seems to motivate the fellow Florentine Luca Pacioli, who describes 3×3 to 9×9 squares in his work De Viribus Quantitatis by the end of 15th century.

Europe after 15th century

The planetary squares had disseminated into northern Europe by the end of 15th century. For instance, the Cracow manuscript of Picatrix from Poland displays magic squares of orders 3 to 9. The same set of squares as in the Cracow manuscript later appears in the writings of Paracelsus in Archidoxa Magica, although in highly garbled form. In 1514 Albrecht Dürer immortalized a 4×4 square in his famous engraving Melencolia I. Paracelsus' contemporary Heinrich Cornelius Agrippa von Nettesheim published his famous three volume book De occulta philosophia in 1531, where he devoted Chapter 22 of Book II to the planetary squares shown below. The same set of squares given by Agrippa reappear in 1539 in Practica Arithmetice by Girolamo Cardano. The tradition of planetary squares was continued into the 17th century by Athanasius Kircher in Oedipi Aegyptici. In Germany, mathematical treaties concerning magic squares were written in 1544 by Michael Stifel in Arithmetica Integra, who rediscovered the bordered squares, and Adam Riese, who rediscovered the continuous numbering method to construct odd ordered squares published by Agrippa. However, due to the religious upheavals of that time, these work were unknown to the rest of Europe.
492
357
816

414151
97612
511108
162313

11247203
41225816
17513219
101811422
23619215

632334351
7112728830
191416152324
182022211713
25291092612
365334231

2247164110354
5234817421129
3062449183612
1331725431937
3814321264420
213983322745
461540934328

858595462631
4915145253111056
4123224445191848
3234352928383925
4026273736303133
1747462021434224
955541213515016
642361606757

37782970216213545
63879307122631446
47739803172235515
16488408132642456
57174994173336525
26581850142743466
67275910512437535
36681960115234476
77286920611253445

In 1624 France, Claude Gaspard Bachet described the "diamond method" for constructing Agrippa's odd ordered squares in his book Problèmes Plaisants. Blaise Pascal, Bernard Frenicle de Bessy and Pierre Fermat are also known to have constructed concentric bordered magic squares, while an early account of their method was given by Antoine Arnauld in his Nouveaux éléments de géométrie. In the two treatise Des quarrez magiques and Table générale des quarrez magiques de quatre de côté, published posthumously in 1693, twenty years after his death, Bernard Frenicle de Bessy demonstrated that there were exactly 880 distinct magic squares of order four and also gave methods to find magic squares of any even order. In 1691, Simon de la Loubère described the Indian continuous method of constructing odd ordered magic squares in his book Du Royaume de Siam, which he had learned while returning from a diplomatic mission to Siam, which was faster than Bachet's method. In an attempt to explain its working, de la Loubere used the primary numbers and root numbers, and rediscovered the method of adding two preliminary squares. This method was further investigated by Abbe Poignard in Traité des quarrés sublimes, by Philippe de La Hire in Mémoires de l’Académie des Sciences for the Royal Academy, and by Joseph Sauveur in Construction des quarrés magiques. Concentric bordered squares were also studied by De la Hire in 1705, while Sauveur introduced magic cubes and lettered squares, which was taken up later by Euler in 1776, who is often credited for devising them. In 1750 d'Ons-le-Bray rediscovered the method of constructing doubly even and singly even squares using bordering technique; while in 1767 Benjamin Franklin published a semi-magic square which had the properties of eponymous Franklin square. By this time the earlier mysticism attached to the magic squares had completely vanished, and the subject was treated as a part of recreational mathematics.
In the 19th century, Bernard Violle gave a comprehensive treatment of magic squares in his three volume Traité complet des carrés magiques, which also described magic cubes, parallelograms, parallelopipeds, and circles. Pandiagonal squares were extensively studied by Andrew Hollingworth Frost, who learned it while in the town of Nasik, India, in a series of articles: On the knight's path, On the General Properties of Nasik Squares, On the General Properties of Nasik Cubes, On the construction of Nasik Squares of any order. He showed that it is impossible to have normal singly-even pandiagonal magic square. Frederick A.P. Barnard constructed inlaid magic squares and other three dimensional magic figures like magic spheres and magic cylinders in Theory of magic squares and of magic cubes. In 1897, Emroy McClintock published On the most perfect form of magic squares, coining the words pandiagonal square and most perfect square, which had previously been referred to as perfect, or diabolic, or Nasik.

Some famous magic squares

Luo Shu magic square

Legends dating from as early as 650 BCE tell the story of the Lo Shu or "scroll of the river Lo". According to the legend, there was at one time in ancient China a huge flood. While the great king Yu was trying to channel the water out to sea, a turtle emerged from it with a curious pattern on its shell: a 3×3 grid in which circular dots of numbers were arranged, such that the sum of the numbers in each row, column and diagonal was the same: 15. According to the legend, thereafter people were able to use this pattern in a certain way to control the river and protect themselves from floods. The Lo Shu Square, as the magic square on the turtle shell is called, is the unique normal magic square of order three in which 1 is at the bottom and 2 is in the upper right corner. Every normal magic square of order three is obtained from the Lo Shu by rotation or reflection.

Magic square in Parshavnath temple

There is a well-known 12th-century 4×4 magic square inscribed on the wall of the Parshvanath temple in Khajuraho, India.
712114
213811
163105
96154

This is known as the Chautisa Yantra since its magic sum is 34. It is one of the three 4×4 pandiagonal magic squares and is also an instance of the most-perfect magic square. The study of this square led to the appreciation of pandiagonal squares by European mathematicians in the late 19th century. Pandiagonal squares were referred to as Nasik squares or Jain squares in older English literature.

Albrecht Dürer's magic square

The order four magic square Albrecht Dürer immortalized in his 1514 engraving Melencolia I, referred to above, is believed to be the first seen in European art. The square associated with Jupiter appears as a talisman used to drive away melancholy. It is very similar to Yang Hui's square, which was created in China about 250 years before Dürer's time. The sum 34 can be found in the rows, columns, diagonals, each of the quadrants, the center four squares, and the corner squares. This sum can also be found in the four outer numbers clockwise from the corners and likewise the four counter-clockwise, the two sets of four symmetrical numbers, the sum of the middle two entries of the two outer columns and rows, and in four kite or cross shaped quartets. The two numbers in the middle of the bottom row give the date of the engraving: 1514. The numbers 1 and 4 at either side of the date correspond respectively to the letters "A" and "D," which are the initials of the artist.
163213
510118
96712
415141

Dürer's magic square can also be extended to a magic cube.

Sagrada Família magic square

The Passion façade of the Sagrada Família church in Barcelona, conceptualized by Antoni Gaudí and designed by sculptor Josep Subirachs, features a 4×4 magic square: The magic constant of the square is 33, the age of Jesus at the time of the Passion. Structurally, it is very similar to the Melancholia magic square, but it has had the numbers in four of the cells reduced by 1.
114144
11769
810105
132315

While having the same pattern of summation, this is not a normal magic square because two numbers are duplicated and two are absent, failing the rule that an n x n magic square must contain each of the positive integers 1 through n 2. Such squares are not generally mathematically interesting and may be described, in a slightly deprecative sense, as trivial. Lee Sallows has pointed out that, due to Subirachs's ignorance of magic square theory, the renowned sculptor made a needless blunder, and supports this assertion by giving several examples of non-trivial 4 x 4 magic squares showing the desired magic constant of 33.
Similarly to Dürer's magic square, the Sagrada Familia's magic square can also be extended to a magic cube.

Parker's magic square

Parker's square is an attempt to create the 3x3 magic square of squares. Parker's Square is not a magic square: it uses some numbers more than once, and the diagonal -- sums to, not as for all the other rows, columns, or diagonal. Parker's Square became a "mascot for people who give it a go, but ultimately fall short". It is also a metaphor for something that is almost right, but is a little off.

Properties of magic squares

Magic constant

The constant that is the sum of any row, or column, or diagonal is called the magic constant or magic sum, M. Every normal magic square has a constant dependent on the order, calculated by the formula. This can be demonstrated by noting that the sum of is. Since the sum of each row is, the sum of rows is, which when divided by the order yields the magic constant. For normal magic squares of orders n = 3, 4, 5, 6, 7, and 8, the magic constants are, respectively: 15, 34, 65, 111, 175, and 260.
For example, a normal 6x6 square will always equate to 111 for each row, column, or diagonal.

Magic square of order 1 is trivial

The 1×1 magic square, with only one cell containing the number 1, is called trivial, because it is typically not under consideration when discussing magic squares; but it is indeed a magic square by definition, if we regard a single cell as a square of order one.

Magic square of order 2 cannot be constructed

Normal magic squares of all sizes can be constructed except 2×2.

Number of magic squares of a given order

Excluding rotations and reflections, there is exactly one 3×3 magic square, there are exactly 880 4×4 magic squares, and there are exactly 275,305,224 5×5 magic squares. For the 6×6 case, there are estimated to be approximately 1.8 × 1019 squares.

Center of mass

If we think of the numbers in the magic square as masses located in various cells, then the center of mass of a magic square coincides with its geometric center.

Moment of inertia

The moment of inertia of a magic square has been defined as the sum over all cells of the number in the cell times the squared distance from the center of the cell to the center of the square; here the unit of measurement is the width of one cell. Then all magic squares of a given order have the same moment of inertia as each other. For the order-3 case the moment of inertia is always 60, while for the order-4 case the moment of inertia is always 340. In general, for the n×n case the moment of inertia is

Birkhoff–von Neumann decomposition

Dividing each number of the magic square by the magic constant will yield a doubly stochastic matrix, whose row sums and column sums equal to unity. However, unlike the doubly stochastic matrix, the diagonal sums of such matrices will also equal to unity. Thus, such matrices constitute a subset of doubly stochastic matrix. The Birkhoff–von Neumann theorem states that for any doubly stochastic matrix, there exists real numbers, where and permutation matrices such that
This representation may not be unique in general. By Marcus-Ree theorem, however, there need not be more than terms in any decomposition. Clearly, this decomposition carries over to magic squares as well, since we can recover a magic square from a doubly stochastic matrix by multiplying it by the magic constant.

Classification of magic squares

While the classification of magic squares can be done in many ways, some useful categories are given below. An n×n square array of integers 1, 2,..., n2 is called:

For any magic square

103138
516211
49714
156121

71449
121156
138103
211516

816
357
492

618
753
294

276
951
438

438
951
276

294
753
618

492
357
816

834
159
672

672
159
834

115144
12679
810115
133216

163213
96712
510118
415141

115144
12679
810115
133216

12679
115144
133216
810115

61297
151414
313162
108511

115144
12679
810115
133216

810115
133216
115144
12679

115810
216133
144115
79126

115144
12679
810115
133216

115810
216133
144115
79126

17241815
23571416
46132022
101219213
11182529

213191012
29251118
20221346
81511724
14167235

For associative magic squares

163213
96712
510118
415141

163213
510118
96712
415141

115144
12679
810115
133216

12679
115144
133216
810115

For pan-diagonal magic squares

103138
516211
49714
156121

516211
49714
156121
103138

211516
71449
121156
138103

For bordered magic squares

12316421
151471811
24171392
20819126
53102225

11623421
151471811
24171392
20819126
51032225

12316421
151471811
24171392
20819126
53102225

21112625
41471822
161713910
23819123
11524205

For composite magic squares

There are many ways to construct magic squares, but the standard way is to follow certain configurations/formulas which generate regular patterns. After a magic square has been constructed using some standard formula, the transformations described in the previous section can be applied to yield further magic squares.
Magic squares exist for all values of n, with only one exception: it is impossible to construct a magic square of order 2. Magic squares can be classified into three types according to their order: odd, doubly even and singly even. Odd and doubly even magic squares are easy to generate; the construction of singly even magic squares is more difficult but several methods exist, including the LUX method for magic squares and the Strachey method for magic squares.
Group theory was also used for constructing new magic squares of a given order from one of them.
The numbers of different n×n magic squares for n from 1 to 5, not counting rotations and reflections are: 1, 0, 1, 880, 275305224. The number for n = 6 has been estimated to be
Cross-referenced to the above sequence, a new classification enumerates the magic tori that display these magic squares. The numbers of magic tori of order n from 1 to 5, are: 1, 0, 1, 255, 251449712.

Method for constructing a magic square of order 3

In the 19th century, Édouard Lucas devised the general formula for order 3 magic squares. Consider the following table made up of positive integers a, b and c:
c − bc + c − a
c − cc +
c + ac − c + b

These nine numbers will be distinct positive integers forming a magic square so long as 0 < a < b < c − a and b ≠ 2a. Moreover, every 3×3 square of distinct positive integers is of this form.
In 1997 Lee Sallows discovered that leaving aside rotations and reflections, then every distinct parallelogram drawn on the Argand diagram defines a unique 3×3 magic square, and vice versa, a result that had never previously been noted.

A method for constructing a magic square of odd order

A method for constructing magic squares of odd order was published by the French diplomat de la Loubère in his book, A new historical relation of the kingdom of Siam, in the chapter entitled The problem of the magical square according to the Indians. The method operates as follows:
The method prescribes starting in the central column of the first row with the number 1. After that, the fundamental movement for filling the squares is diagonally up and right, one step at a time. If a filled square is encountered, one moves vertically down one square instead, then continues as before. When an "up and to the right" move would leave the square, it is wrapped around to the last row or first column, respectively.
1

1
2

1
3
2

1
3
42

1
35
42

16
35
42

16
357
42

816
357
42

816
357
492

Starting from other squares rather than the central column of the first row is possible, but then only the row and column sums will be identical and result in a magic sum, whereas the diagonal sums will differ. The result will thus be a semimagic square and not a true magic square. Moving in directions other than north east can also result in magic squares.
816
357
492

17241815
23571416
46132022
101219213
11182529

47586980112233445
57687991122334446
67788102132435456
77718203142535566
61719304152636576
16272940516264755
26283950617274415
36384960717331425
37485970812132435

A method of constructing a magic square of doubly even order

means that n is an even multiple of an even integer; or 4p, where p is an integer.
Generic pattern
All the numbers are written in order from left to right across each row in turn, starting from the top left hand corner. Numbers are then either retained in the same place or interchanged with their diametrically opposite numbers in a certain regular pattern. In the magic square of order four, the numbers in the four central squares and one square at each corner are retained in the same place and the others are interchanged with their diametrically opposite numbers.
A construction of a magic square of order 4
Starting from top left, go left to right through each row of the square, counting each cell from 1 to 16 and filling the cells along the diagonals with its corresponding number. Once the bottom right cell is reached, continue by going right to left, starting from the bottom right of the table through each row, and fill in the non-diagonal cells counting up from 1 to 16 with its corresponding number. As shown below:
14
67
1011
1316

115144
12679
810115
133216

An extension of the above example for Orders 8 and 12
First generate a pattern table, where a '1' indicates selecting from the square where the numbers are written in order 1 to n2, and a '0' indicates selecting from the square where the numbers are written in reverse order n2 to 1. For M = 4, the pattern table is as shown below. When we shade the unaltered cells, we get a criss-cross pattern.
1234
5678
9101112
13141516

16151413
1211109
8765
4321

1001
0110
0110
1001

115144
12679
810115
133216

The patterns are a) there are equal number of '1's and '0's in each row and column; b) each row and each column are "palindromic"; c) the left- and right-halves are mirror images; and d) the top- and bottom-halves are mirror images. The pattern table can be denoted using hexadecimals as for simplicity. The simplest method of generating the required pattern for higher ordered doubly even squares is to copy the generic pattern for the fourth-order square in each four-by-four sub-squares.
For M = 8, possible choices for the pattern are ; ; .
10011001
01100110
01100110
10011001
10011001
01100110
01100110
10011001

1458
10111415
18192223
25282932
33363740
42434647
50515455
57606164

163624559588
5610115352141549
4818194544222341
2539382829353432
3331303637272640
2442432120464717
165051131254559
577660613264

For M = 12, the pattern table yields a magic square It is possible to count the number of choices one has based on the pattern table, taking rotational symmetries into account.

Method of superposition

The earliest discovery of the superposition method was made by the Indian mathematician Narayana in the 14th century. The same method was later re-discovered and studied in early 18th century Europe by de la Loubere, Poignard, de La Hire, and Sauveur; and the method is usually referred to as de la Hire's method. Although Euler's work on magic square was unoriginal, he famously conjectured the impossibility of constructing the evenly odd ordered mutually orthogonal Graeco-Latin squares. This conjecture was disproved in the mid 20th century. For clarity of exposition, we have distinguished two important variations of this method.

Euler's method

This method consists in constructing two preliminary squares, which when added together gives the magic square. As a running example, we will consider a 3×3 magic square. We can uniquely label each number of the 3×3 natural square by a pair of numbers as
123
456
789

αaαbαc
βaβbβc
γaγbγc

where every pair of Greek and Latin alphabets, e.g. αa, are meant to be added together, i.e. αa = α + a. Here, = and =. The numbers 0, 3, and 6 are referred to as the root numbers while the numbers 1, 2, and 3 are referred to as the primary numbers. An important general constraint here is
Thus, the original square can now be split into two simpler squares:
ααα
βββ
γγγ

abc
abc
abc

The lettered squares are referred to as Greek square or Latin square if they are filled with Greek or Latin letters, respectively. A magic square can be constructed by ensuring that the Greek and Latin squares are magic squares too. The converse of this statement is also often, but not always, true: A magic square can be decomposed into a Greek and a Latin square, which are themselves magic squares. Thus the method is useful for both synthesis as well as analysis of a magic square. Lastly, by examining the pattern in which the numbers are laid out in the finished square, it is often possible to come up with a faster algorithm to construct higher order squares that replicate the given pattern, without the necessity of creating the preliminary Greek and Latin squares.
During the construction of the 3×3 magic square, the Greek and Latin squares with just three unique terms are much easier to deal with than the original square with nine different terms. The row sum and the column sum of the Greek square will be the same, α + β + γ, if
This can be achieved by cyclic permutation of α, β, and γ. Satisfaction of these two conditions ensures that the resulting square is a semi-magic square; and such Greek and Latin squares are said to be mutually orthogonal to each other. For a given order n, there are at most n - 1 squares in a set of mutually orthogonal squares, not counting the variations due to permutation of the symbols. This upper bound is exact when n is a prime number.
In order to construct a magic square, we should also ensure that the diagonals sum to magic constant. For this, we have a third condition:
The mutually orthogonal Greek and Latin squares that satisfy the first part of the third condition are said to be mutually orthogonal doubly diagonal Graeco-Latin squares.
Odd squares: For the 3×3 odd square, since α, β, and γ are in arithmetic progression, their sum is equal to the product of the square's order and the middle term, i.e. α + β + γ = 3 β. Thus, the diagonal sums will be equal if we have βs in the main diagonal and α, β, γ in the skew diagonal. Similarly, for the Latin square. The resulting Greek and Latin squares and their combination will be as below. The Latin square is just a 90 degree anti-clockwise rotation of the Greek square with the corresponding letters interchanged. Substituting the values of the Greek and Latin letters will give the 3×3 magic square.
βαγ
γβα
αγβ

cab
abc
bca

βcαaγb
γaβbαc
αbγcβa

618
753
294

For the odd squares, this method explains why the Siamese method and its variants work. This basic method can be used to construct odd ordered magic squares of higher orders. To summarise:
A peculiarity of the construction method given above for the odd magic squares is that the middle number /2 will always appear at the center cell of the magic square. Since there are ! ways to arrange the skew diagonal terms, we can obtain ! Greek squares this way; same with the Latin squares. Also, since each Greek square can be paired with ! Latin squares, and since for each of Greek square the middle term may be arbitrarily placed in the main diagonal or the skew diagonal, we can construct a total of 2 × ! × ! magic squares using this method. For n = 3, 5, and 7, this will give 8, 1152, and 1,036,800 different magic squares, respectively. Dividing by 8 to neglect equivalent squares due to rotation and reflections, we obtain 1, 144, and 129,600 essentially different magic squares, respectively.
As another example, the construction of 5×5 magic square is given. Numbers are directly written in place of alphabets. The numbered squares are referred to as primary square or root square if they are filled with primary numbers or root numbers, respectively. The numbers are placed about the skew diagonal in the root square such that the middle column of the resulting root square has 0, 5, 10, 15, 20. The primary square is obtained by rotating the root square counter-clockwise by 90 degrees, and replacing the numbers. The resulting square is an associative magic square, in which every pair of numbers symmetrically opposite to the center sum up to the same value, 26. For e.g., 16+10, 3+23, 6+20, etc. In the finished square, 1 is placed at center cell of bottom row, and successive numbers are placed via elongated knight's move, or equivalently, bishop's move. When a collision occurs, the break move is to move one cell up. All the odd numbers occur inside the central diamond formed by 1, 5, 25 and 21, while the even numbers are placed at the corners. The occurrence of the even numbers can be deduced by copying the square to the adjacent sides. The even numbers from four adjacent squares will form a cross.
105
1020
10
010
1510

1005
10200
010
010
15010

10152005
51015200
05101520
20051015
15200510

21543
15432
54321
43215
32154

12162548--
61519232--
59131721--
2437112024-
1822110141822
481216
26-

A variation of the above example, where the skew diagonal sequence is taken in different order, is given below. The resulting magic square is the flipped version of the famous Agrippa's Mars magic square. It is an associative magic square and is the same as that produced by Moschopoulos's method. Here the resulting square starts with 1 placed in the cell which is to the right of the centre cell, and proceeds as De la Loubere's method, with downwards-right move. When a collision occurs, the break move is to shift two cells to the right.
1020
1015
10
510
010

10015520
20100155
52010015
15520100
01552010

14253
42531
25314
53142
31425

114171023
24125186
72513119
20821142
31692215

In the previous examples, for the Greek square, the second row can be obtained from the first row by circularly shifting it to the right by one cell. Similarly, the third row is a circularly shifted version of the second row by one cell to the right; and so on. Likewise, the rows of the Latin square is circularly shifted to the left by one cell. The row shifts for the Greek and Latin squares are in mutually opposite direction. It is possible to circularly shift the rows by more than one cell to create the Greek and Latin square.
This essentially re-creates the knight's move. All the letters will appear in both the diagonals, ensuring correct diagonal sum. Since there are n! permutations of the Greek letters by which we can create the first row of the Greek square, there are thus n! Greek squares that can be created by shifting the first row in one direction. Likewise, there are n! such Latin squares created by shifting the first row in the opposite direction. Since a Greek square can be combined with any Latin square with opposite row shifts, there are n! × n! such combinations. Lastly, since the Greek square can be created by shifting the rows either to the left or to the right, there are a total of 2 × n! × n! magic squares that can be formed by this method. For n = 5 and 7, since they are prime numbers, this method creates 28,800 and 50,803,200 pandiagonal magic squares. Dividing by 8 to neglect equivalent squares due to rotation and reflections, we obtain 3,600 and 6,350,400 equivalent squares. Further dividing by n2 to neglect equivalent panmagic squares due to cyclic shifting of rows or columns, we obtain 144 and 129,600 essentially different panmagic squares. For order 5 squares, these are the only panmagic square there are. The condition that the square's order not be divisible by 3 means that we cannot construct squares of orders 9, 15, 21, 27, and so on, by this method.
In the example below, the square has been constructed such that 1 is at the center cell. In the finished square, the numbers can be continuously enumerated by the knight's move. When collision occurs, the break move is to move one cell up, one cell left.The resulting square is a pandiagonal magic square. This square also has a further diabolical property that any five cells in quincunx pattern formed by any odd sub-square, including wrap around, sum to the magic constant, 65. For e.g., 13+7+1+20+24, 23+1+9+15+17, 13+21+10+19+2 etc. Also the four corners of any 5×5 square and the central cell, as well as the middle cells of each side together with the central cell, including wrap around, give the magic sum: 13+10+19+22+1 and 20+24+12+8+1. Lastly the four rhomboids that form elongated crosses also give the magic sum: 23+1+9+24+8, 15+1+17+20+12, 14+1+18+13+19, 7+1+25+22+10.
10152005
05101520
15200510
51015200
20051015

31425
42531
53142
14253
25314

131624210
47151821
20231912
61417253
22581119

We can also combine the Greek and Latin squares constructed by different methods. In the example below, the primary square is made using knight's move. We have re-created the magic square obtained by De la Loubere's method. As before, we can form 8 × ! × n! magic squares by this combination. For n = 5 and 7, this will create 23,040 and 29,030,400 magic squares. After dividing by 8 in order to neglect equivalent squares due to rotation and reflection, we get 2,880 and 3,628,800 squares.
15200510
20051015
05101520
51015200
10152005

24135
35241
41352
52413
13524

17241815
23571416
46132022
101219213
11182529

For order 5 squares, these three methods give a complete census of the number of magic squares that can be constructed by the method of superposition. Neglecting the rotation and reflections, the total number of magic squares of order 5 produced by the superposition method is 144 + 3,600 + 2,880 = 6,624.
Even squares: We can also construct even ordered squares in this fashion. Since there is no middle term among the Greek and Latin alphabets for even ordered squares, in addition to the first two constraint, for the diagonal sums to yield the magic constant, all the letters in the alphabet should appear in the main diagonal and in the skew diagonal.
An example of a 4×4 square is given below. For the given diagonal and skew diagonal in the Greek square, the rest of the cells can be filled using the condition that each letter appear only once in a row and a column.
αδ
δα
γβ
βγ

αβγδ
γδαβ
δγβα
βαδγ

abcd
dcba
badc
cdab

αaβbγcδd
γdδcαbβa
δbγaβdαc
βcαdδaγb

161116
121525
14983
741310

Using these two Graeco-Latin squares, we can construct 2 × 4! × 4! = 1,152 magic squares. Dividing by 8 to eliminate equivalent squares due to rotation and reflections, we get 144 essentially different magic squares of order 4. These are the only magic squares constructible by the Euler method, since there are only two mutually orthogonal doubly diagonal Graeco-Latin squares of order 4.
Similarly, an 8×8 magic square can be constructed as below. Here the order of appearance of the numbers is not important; however the quadrants imitate the layout pattern of the 4×4 Graeco-Latin squares.
08162432404856
24168056484032
48563240162408
40325648802416
56484032241680
32404856081624
80241640325648
16240848563240

12345678
34127856
56781234
78563412
43218765
21436587
87654321
65872143

110192837465564
27209263564538
536239481726312
474061541142518
605142333223145
344152596132431
167302144355849
222981550573643

Euler's method has given rise to the study of Graeco-Latin squares. Euler's method for constructing magic squares is valid for any order except 2 and 6.
Variations: Magic squares constructed from mutually orthogonal doubly diagonal Graeco-Latin squares are interesting in themselves since the magic property emerges from the relative position of the alphabets in the square, and not due to any arithmetic property of the value assigned to them. This means that we can assign any value to the alphabets of such squares and still obtain a magic square. This is the basis for constructing squares that display some information in the square and for creating "reversible squares". For example, we can display the number π ≈ 3·141592 at the bottom row of a 4×4 magic square using the Graeco-Latin square given above by assigning = and =. We will obtain the following non-normal magic square with the magic sum 124:
1029319
9418120
1790413
3141592

Narayana-De la Hire's method for even orders

Narayana-De la Hire's method for odd square is the same as that of Euler's. However, for even squares, we drop the second requirement that each Greek and Latin letter appear only once in a given row or column. This allows us to take advantage of the fact that the sum of an arithmetic progression with an even number of terms is equal to the sum of two opposite symmetric terms multiplied by half the total number of terms. Thus, when constructing the Greek or Latin squares,
As a running example, if we take a 4×4 square, where the Greek and Latin terms have the values = and =, respectively, then we have α + β + γ + δ = 2 = 2. Similarly, a + b + c + d = 2 = 2. This means that the complementary pair α and δ can appear twice in a column and still give the desired magic sum. Thus, we can construct:
In the example given below, the main diagonal is filled with sequence ordered as α, β, γ, δ, while the skew diagonal filled in the same order. The remaining cells are then filled column wise such that the complementary letters appears only once within a row, but twice within a column. In the first column, since α appears on the 1st and 4th row, the remaining cells are filled with its complementary term δ. Similarly, the empty cells in the 2nd column are filled with γ; in 3rd column β; and 4th column α. Each Greek letter appears only once along the rows, but twice along the columns. As such, the row sums are α + β + γ + δ while the column sums are either 2 or 2. Likewise for the Latin square, which is obtained by flipping the Greek square along the main diagonal and interchanging the corresponding letters.
αδ
βγ
βγ
αδ

αγβδ
δβγα
δβγα
αγβδ

adda
cbbc
bccb
daad

αaγdβdδa
δcβbγbαc
δbβcγcαb
αdγaβaδd

112813
156103
147112
49516

The above example explains why the "criss-cross" method for doubly even magic square works. Another possible 4×4 magic square, which is also pan-diagonal as well as most-perfect, is constructed below using the same rule. However, the diagonal sequence is chosen such that all four letters α, β, γ, δ appear inside the central 2×2 sub-square. Remaining cells are filled column wise such that each letter appears only once within a row. In the 1st column, the empty cells need to be filled with one of the letters selected from the complementary pair α and δ. Given the 1st column, the entry in the 2nd row can only be δ since α is already there in the 2nd row; while, in the 3rd row the entry can only be α since δ is already present in the 3rd row. We proceed similarly until all cells are filled. The Latin square given below has been obtained by flipping the Greek square along the main diagonal and replacing the Greek alphabets with corresponding Latin alphabets.
αγ
γα
βδ
δβ

αβδγ
δγαβ
αβδγ
δγαβ

adad
bcbc
dada
cbcb

αaβdδaγd
δbγcαbβc
αdβaδdγa
δcγbαcβb

181312
141127
45169
151036

We can use this approach to construct singly even magic squares as well. However, we have to be more careful in this case since the criteria of pairing the Greek and Latin alphabets uniquely is not automatically satisfied. Violation of this condition leads to some missing numbers in the final square, while duplicating others. Thus, here is an important proviso:
Below is a construction of a 6×6 magic square, where the numbers are directly given, rather than the alphabets. The second square is constructed by flipping the first square along the main diagonal. Here in the first column of the root square the 3rd cell is paired with its complement in the 4th cells. Thus, in the primary square, the numbers in the 1st and 6th cell of the 3rd row are same. Likewise, with other columns and rows. In this example the flipped version of the root square satisfies this proviso.
0241812630
3061218240
0241218630
3024121860
3061812240
0618122430

161661
525522
433344
344433
252255
616116

13019181231
3581723262
42715211034
3328162293
32112014295
6724132536

As another example of a 6×6 magic square constructed this way is given below. Here the diagonal entries are arranged differently. The primary square is constructed by flipping the root square about the main diagonal. In the second square the proviso for singly even square is not satisfied, leading to a non-normal magic square where the numbers 3, 13, 24, and 34 are duplicated while missing the numbers 4, 18, 19, and 33.
6301218024
2401218306
2401812306
6301812024
2401812306
6301218024

255252
611616
334443
443334
166161
522525

8351720526
30113243112
2732216349
10342115328
2562413367
11321423229

The last condition is a bit arbitrary and may not always need to be invoked, as in this example:
6301224180
6018241230
2401261830
6301824120
2430126180
2401861230

225255
611661
343434
552522
434343
166116

8321726235
12119301831
27415102134
11352029142
2833169223
2562471336

As one more example, we have generated an 8×8 magic square. Unlike the criss-cross pattern of the earlier section for evenly even square, here we have a checkered pattern for the altered and unaltered cells. Also, in each quadrant the odd and even numbers appear in alternating columns.
04816322440856
56840243216480
04816322440856
56840243216480
56840243216480
04816322440856
56840243216480
04816322440856

18188181
72722727
36366363
54544545
45455454
63633636
27277272
81811818

156174032411657
631047263423507
354193830431459
611245283621425
601344293720534
651223527461162
581542313918552
84724332548964

Variations: A number of variations of the basic idea are possible: a complementary pair can appear n/2 times or less in a column. That is, a column of a Greek square can be constructed using more than one complementary pair. This method allows us to imbue the magic square with far richer properties. The idea can also be extended to the diagonals too. An example of an 8×8 magic square is given below. In the finished square each of four quadrants are pan-magic squares as well, each quadrant with same magic constant 130.
04856816324024
56804840241632
04856816324024
56804840241632
48085632162440
85648024403216
48085632162440
85648024403216

18187272
72721818
81812727
27278181
36365454
54543636
63634545
45456363

156571623344726
631075041321740
84964918394231
581525548252433
516116237202944
136053427463522
543145936212845
126152530433819

Method of borders

Bordering method for order 3

In this method, the objective is to wrap a border around a smaller magic square which serves as a core. Consider the 3×3 square for example. Subtracting the middle number 5 from each number 1, 2,..., 9, we obtain 0, ± 1, ± 2, ± 3, and ± 4, which we will, for lack of better words, following Harry S. White, refer to as bone numbers. The magic constant of a magic square, which we will refer to as the skeleton square, made by these bone numbers will be zero since adding all the rows of a magic square will give nM = Σ k = 0; thus M = 0.
It is not difficult to argue that the middle number should be placed at the center cell: let x be the number placed in the middle cell, then the sum of the middle column, middle row, and the two diagonals give Σ k + 3 x = 4 M. Since Σ k = 3 M, we have x = M / 3. Here M = 0, so x = 0.
Putting the middle number 0 in the center cell, we want to construct a border such that the resulting square is magic. Let the border be given by:
uav
b*0b
v*a*u*

Since the sum of each row, column, and diagonals must be a constant, we have
Now, if we have chosen a, b, u, and v, then we have a* = - a, b* = - b, u* = - u, and v* = - v. This means that if we assign a given number to a variable, say a = 1, then its complement will be assigned to a*, i.e. a* = - 1. Thus out of eight unknown variables, it is sufficient to specify the value of only four variables. We will consider a, b, u, and v as independent variables, while a*, b*, u*, and v* as dependent variables. This allows us to consider a bone number ± x as a single number regardless of sign because its assignment to a given variable, say a, will automatically imply that the same number of opposite sign will be shared with its complement a*, and two independent variables, say a and b, cannot be assigned the same bone number. But how should we choose a, b, u, and v? We have the sum of the top row and the sum of the right column as
Since 0 is an even number, there are only two ways that the sum of three integers will yield an even number: 1) if all three were even, or 2) if two were odd and one was even. Since in our choice of numbers we only have two even non-zero number, the first statement is false. Hence, it must be the case that the second statement is true: that two of the numbers are odd and one even.
The only way that both the above two equations can satisfy this parity condition simultaneously, and still be consistent with the set of numbers we have, is when u and v are odd. For on the contrary, if we had assumed u and a to be odd and v to be even in the first equation, then u* = - u will be odd in the second equation, making b odd as well, in order to satisfy the parity condition. But this requires three odd numbers, contradicting the fact that we only have two odd numbers which we can use. This proves that the odd bone numbers occupy the corners cells. When converted to normal numbers by adding 5, this implies that the corners of a 3×3 magic square are all occupied by even numbers.
Thus, taking u = 1 and v = 3, we have a = - 4 and b = - 2. Hence, the finished skeleton square will be as in the left. Adding 5 to each number, we get the finished magic square.
1-43
20-2
-34-1

618
753
294

Similar argument can be used to construct larger squares. Since there does not exist a 2×2 magic square around which we can wrap a border to construct a 4×4 magic square, the next smallest order for which we can construct bordered square is the order 5.

Bordering method for order 5

Consider the fifth-order square. For this, we have a 3×3 magic core, around which we will wrap a magic border. The bone numbers to be used will be ± 5, ± 6, ± 7, ± 8, ± 9, ± 10, ± 11, and ± 12. Disregarding the signs, we have 8 bone numbers, 4 of which are even and 4 of which are odd. In general, for a square of any order n, there will be 4 border cells, which are to be filled using 2 bone numbers. Let the magic border be given as
uabcv
d*d
e*e
f*f
v*a*b*c*u*

As before, we should
It is sufficient to determine the numbers u, v, a, b, c, d, e, f to describe the magic border. As before, we have the two constraint equations for the top row and right column:
Multiple solutions are possible. The standard procedure is to
There are 28 ways of choosing two numbers from the set of 8 bone numbers for the corner cells u and v. However, not all pairs are admissible. Among the 28 pairs, 16 pairs are made of an even and an odd number, 6 pairs have both as even numbers, while 6 pairs have them both as odd numbers.
We can prove that the corner cells u and v cannot have an even and an odd number. This is because if this were so, then the sums u + v and v + u* will be odd, and since 0 is an even number, the sums a + b + c and d + e + f should be odd as well. The only way that the sum of three integers will result in an odd number is when 1) two of them are even and one is odd, or 2) when all three are odd. Since the corner cells are assumed to be odd and even, neither of these two statements are compatible with the fact that we only have 3 even and 3 odd bone numbers at our disposal. This proves that u and v cannot have different parity. This eliminates 16 possibilities.
Using similar type reasoning we can also draw some conclusions about the sets and. If u and v are both even, then both the sets should have two odd numbers and one even number. If u and v are both odd, then one of the sets should have three even numbers while the other set should have one even number and two odd numbers.
As a running example, consider the case when both u and v are even. The 6 possible pairs are:,,,,, and two of them are odd and one is even, or 2) when all three are even. The fact that the two corner cells are even means that we have only 2 even numbers at our disposal. Thus, the second statement is not compatible with this fact. Hence, it must be the case that the first statement is true: two of the three numbers should be odd, while one be even.
Now let a, b, d, e be odd numbers while c and f be even numbers. Given the odd bone numbers at our disposal: ± 5, ± 7, ± 9, and ± 11, their differences range from D = while their sums range from S =. It is also useful to have a table of their sum and differences for later reference. Now, given the corner cells, we can check its admissibility by checking if the sums u + v + c and v + u* + f fall within the set D or S. The admissibility of the corner numbers is a necessary but not a sufficient condition for the solution to exist.
For example, if we consider the pair =, then u + v = 20 and v + u* = 6; and we will have ± 6 and ± 10 even bone numbers at our disposal. Taking c = ± 6, we have the sum u + v + c to be 26 and 14, depending on the sign of ± 6 taken, both of which do not fall within the sets D or S. Likewise, taking c = ± 10, we have the sum u + v + c to be 30 and 10, both of which again do not fall within the sets D or S. Thus, the pair is not admissible. By similar process of reasoning, we can also rule out the pair.
As another example, if we consider the pair =, then u + v = 22 and v + u* = 2; and we will have ± 6 and ± 8 even bone numbers at our disposal. Taking c = ± 6, we have the sum u + v + c to be 28 and 16. While 28 does not fall within the sets D or S, 16 falls in set S. By inspection, we find that if =, then a + b = -16; and it will satisfy the first constraint equation. Also, taking f = ± 8, we have the sum v + u* + f to be 10 and -6. While 10 does not fall within the sets D or S, -6 falls in set D. Since -7 and -9 have already been assigned to a and b, clearly = so that d + e = 6; and it will satisfy the second constraint equation.
Likewise, taking c = ± 8, we have the sum u + v + c to be 30 and 14. While 30 does not fall within the sets D or S, 14 falls in set S. By inspection, we find that if =, then a + b = -14. Also, taking f = ± 6, we have the sum v + u* + f to be 8 and -4. While 8 does not fall within the sets D or S, -4 falls in set D. Clearly, = so that d + e = 4, and the second constraint equation will be satisfied.
Hence the corner pair = is admissible; and it admits two solutions: = and =. The finished skeleton squares are given below. The magic square is obtained by adding 13 to each cells.
10-7-9-612
5-5
-1111
8-8
-12796-10

2364725
188
224
215
12022193

10-5-9-812
7-7
-1111
6-6
-12598-10

2384525
206
224
197
11822213

Using similar process of reasoning, we can construct the following table for the values of u, v, a, b, c, d, e, f expressed as bone numbers as given below. There are only 6 possible choices for the corner cells, which leads to 10 possible border solutions.
u, va, b, cd, e, f
12, 10-6, -7, -9-11, 5, 8
12, 10-5, -8, -9-11, 6, 7
11, 56, -10, -12-9, 7, 8
10, 65, -9, -12-11, 7, 8
10, 67, -11, -12-9, 5, 8
9, 75, -10, -11-12, 6, 8
9, 76, -10, -12-11, 5, 8
8, 67, -10, -11-12, 5, 9
8, 69, -11, -12-10, 5, 7
7, 59, -10, -11-12, 6, 8

Given this group of 10 borders, we can construct 10×8×2 = 2880 essentially different bordered magic squares. Here the bone numbers ± 5,..., ± 12 were consecutive. More bordered squares can be constructed if the numbers are not consecutive. If non-consecutive bone numbers were also used, then there are a total of 605 magic borders. Thus, the total number of order 5 essentially different bordered magic squares is 174,240. See history. It is worth noting that the number of fifth-order magic squares constructible via the bordering method is about 26 times larger than via the superposition method.

Continuous enumeration methods

Exhaustive enumeration of all the borders of a magic square of a given order, as done previously, is very tedious. As such a structured solution is often desirable, which allows us to construct a border for a square of any order. Below we give three algorithms for constructing border for odd, doubly even, and singly even squares. These continuous enumeration algorithms were discovered in 10th century by Arab scholars; and their earliest surviving exposition comes from the two treatises by al-Buzjani and al-Antaki, although they themselves were not the discoverers. Since then many more such algorithms have been discovered.
Odd ordered squares: The following is the algorithm given by al-Buzjani to construct a border for odd squares. A peculiarity of this method is that for order n square, the two adjacent corners are numbers n - 1 and n + 1.
Starting from the cell above the lower left corner, we put the numbers alternately in left column and bottom row until we arrive at the middle cell. The next number is written in the middle cell of the bottom row just reached, after which we fill the cell in the upper left corner, then the middle cell of the right column, then the upper right corner. After this, starting from the cell above middle cell of the right column already filled, we resume the alternate placement of the numbers in the right column and the top row. Once half of the border cells are filled, the other half are filled by numbers complementary to opposite cells. The subsequent inner borders is filled in the same manner, until the square of order 3 is filled.
Below is an example for 9th-order square.
88078767512141610
672264626126282415
695532525136342713
715747384540352511
73594943413933239
51929423744536377
31748303146506579
15818202156546081
72246770686674

Doubly even order: The following is the method given by al-Antaki. Consider an empty border of order n = 4k with k ≥ 3. The peculiarity of this algorithm is that the adjacent corner cells are occupied by numbers n and n - 1.
Starting at the upper left corner cell, we put the successive numbers by groups of four, the first one next to the corner, the second and the third on the bottom, and the fourth at the top, and so on until there remains in the top row six empty cells. We then write the next two numbers above and the next four below. We then fill the upper corners, first left then right. We place the next number below the upper right corner in the right column, the next number on the other side in the left column. We then resume placing groups of four consecutive numbers in the two columns as before. Once half of the border cells are filled, the other half are filled by numbers complementary to opposite cells.
The example below gives the border for order 16 square.
15125525445251250891024624524424316
24017
18239
19238
23720
23621
22235
23234
23324
23225
26231
27230
22928
22829
30227
241256232532526724924824711121314242

For order 8 square, we just begin directly with the six cells.
712626160598
569
1055
1154
5312
5213
1451
576463345658

Singly even order: For singly even order, we have the algorithm given by al-Antaki. Here the corner cells are occupied by n and n - 1. Below is an example of 10th-order square.
Start by placing 1 at the bottom row next to the left corner cell, then place 2 in the top row. After this, place 3 at the bottom row and turn around the border in anti-clockwise direction placing the next numbers, until n - 2 is reached on the right column. The next two numbers are placed in the upper corners. Then, the next two numbers are placed on the left column, then we resume the cyclic placement of the numbers until half of all the border cells are filled. Once half of the border cells are filled, the other half are filled by numbers complementary to opposite cells.
910029859488158410
8318
1685
8714
1289
1190
938
695
974
91199396713861792

Method of composition

For squares of order ''m × n'' where ''m'', ''n'' > 2

This is a method reminiscent of the Kronecker product of two matrices, that builds an nm × nm magic square from an n × n magic square and an m × m magic square. The "product" of two magic squares creates a magic square of higher order than the two multiplicands. Let the two magic squares be of orders m and n. The final square will be of order m × n. Divide the square of order m × n into m × m sub-squares, such that there are a total of n2 such sub-squares. In the square of order n, reduce by 1 the value of all the numbers. Multiply these reduced values by m2, and place the results in the corresponding sub-squares of the m × n whole square. The squares of order m are added n2 times to the sub-squares of the final square. The peculiarity of this construction method is that each magic subsquare will have different magic sums. The square made of such magic sums from each magic subsquare will again be a magic square. The smallest composite magic square of order 9, composed of two order 3 squares is given below.
816
357
492

636363000454545
636363000454545
636363000454545
181818363636545454
181818363636545454
181818363636545454
272727727272999
272727727272999
272727727272999

716469816534651
666870357485052
677265492495447
261924443742625560
212325394143575961
222720404538586356
352833807378171015
303234757779121416
313629768174131811

Since each of the 3×3 sub-squares can be independently rotated and reflected into 8 different squares, from this single 9×9 composite square we can derive 89 = 134,217,728 essentially different 9×9 composite squares. Plenty more composite magic squares can also be derived if we select non-consecutive numbers in the magic sub-squares, like in Yang Hui's version of the 9×9 composite magic square. The next smallest composite magic squares of order 12, composed of magic squares of order 3 and 4 are given below.
294
753
618

114118
127213
69163
154510

294119126121929994657267
753124122120979593706866
618123118125969198696471
101108103566358111813110117112
106104102615957161412115113111
105100107605562151017114109116
475449748176137144139202722
525048797775142140138252321
514653787380141136143241926
128135130293631384540839085
133131129343230434139888684
132127134332835423744878289

1730272412914213913649625956
2823182914013513014160555061
2225321913413714413154576451
3120212614313213313863525358
971101071046578757233464340
108103981097671667744393445
102105112997073806738414835
1111001011067968697447363742
81949188114118113126123120
92878293127213124119114125
8689968369163118121128115
95848590154510127116117122

For the base squares, there is only one essentially different 3rd order square, while there 880 essentially different 4th-order squares that we can choose from. Each pairing can produce two different composite squares. Since each magic sub-squares in each composite square can be expressed in 8 different forms due to rotations and reflections, there can be 1×880×89 + 880×1×816 ≈ 2.476×1017 essentially different 12×12 composite magic squares created this way, with consecutive numbers in each sub-square. In general, if there are cm and cn essentially different magic squares of order m and n, then we can form cm × cn × composite squares of order mn, provided m ≠ n. If m = n, then we can form 2 × 8m2 composite squares of order m2.

For squares of doubly even order

When the squares are of doubly even order, we can construct a composite magic square in a manner more elegant than the above process, in the sense that every magic subsquare will have the same magic constant. Let n be the order of the main square and m the order of the equal subsquares. The subsquares are filled one by one, in any order, with a continuous sequence of m2/2 smaller numbers together with their complements to n2 + 1. Each subsquare as a whole will yield the same magic sum. The advantage of this type of composite square is that each subsquare is filled in the same way and their arrangement is arbitrary. Thus, the knowledge of a single construction of even order will suffice to fill the whole square. Furthermore, if the subsquares are filled in the natural sequence, then the resulting square will be pandiagonal. The magic sum of the subsquares is related to the magic sum of the whole square by where k = n/m.
In the examples below, we have divided the order 12 square into nine subsquares of order 4 filled each with eight smaller numbers and, in the corresponding bishop's cells, their complements to n2 + 1 = 145. Each subsquare is pandiagonal with magic constant 290; while the whole square on the left is also pandiagonal with magic constant 870.
114213989134131161712612324
1407214113215101331242318125
6137144314129136112212112819
1434513813512131301272021122
25118115323311010740411029948
116312611710839341091004742101
30113120273810511235469710443
11928291141113637106103444598
499491565786836465787572
925550938463588576716677
548996516281885970738067
955253908760618279686974

69747968291141192861828760
75726578115322511883645786
66777671261171163158858463
80677073120273011388596281
2112212720539095521313013512
123241712691564994131169134
1812512423509392551013313215
1281922121965154891361114129
459810344513814343710611136
994841102139811421074033110
4210110047214114073410910839
104434697144361371123538105

In another example below, we have divided the order 12 square into four order 6 squares. Each of the order 6 squares are filled with eighteen small numbers and their complements using bordering technique given by al-Antaki. If we remove the shaded borders of the order 6 subsquares and form an order 8 square, then this order 8 square is again a magic square. In its full generality, we can take any m2/2 smaller numbers together with their complements to n2 + 1 to fill the subsquares, not necessarily in continuous sequence.
608288569059241181242012623
646974796881283311011532117
837572657862119111362911426
846677767161120301131123525
588067707387221163134109123
866357895585122272112519121
613614221445421001063810841
101512813314135465192975099
137129181113281019354479644
138121311301771024895945343
413413161271414098495291105
140931431139104453910737103

Medjig-method for squares of even order 2''n'', where ''n'' > 2

In this method a magic square is "multiplied" with a medjig square to create a larger magic square. The namesake of this method derives from mathematical game called medjig created by Willem Barink in 2006, although the method itself is much older. An early instance of a magic square constructed using this method occurs in Yang Hui's text for order 6 magic square. The LUX method to construct singly even magic squares is a special case of the medjig method, where only 3 out of 24 patterns are used to construct the medjig square.
The pieces of the medjig puzzle are 2×2 squares on which the numbers 0, 1, 2 and 3 are placed. There are three basic patterns by which the numbers 0, 1, 2 and 3 can be placed in a 2×2 square, where 0 is at the top left corner:
01
23

01
32

02
31

Each pattern can be reflected and rotated to obtain 8 equivalent patterns, giving us a total of 3×8 = 24 patterns. The aim of the puzzle is to take n2 medjig pieces and arrange them in an n × n medjig square in such a way that each row, column, along with the two long diagonals, formed by the medjig square sums to 3n, the magic constant of the medjig square. An n × n medjig square can create a 2n × 2n magic square where n > 2.
Given an n×n medjig square and an n×n magic square base, a magic square of order 2n×2n can be constructed as follows:
Assuming that we have an initial magic square base, the challenge lies in constructing a medjig square. For reference, the sums of each medjig piece along the rows, columns and diagonals, denoted in italics, are:
101-
523-
3243

101-
532-
4332

202-
431-
5331

Doubly even squares: The smallest even ordered medjig square is of order 2 with magic constant 6. While it is possible to construct a 2×2 medjig square, we cannot construct a 4×4 magic square from it since 2×2 magic squares required to "multiply" it does not exist. Nevertheless, it is worth constructing these 2×2 medjig squares. There exists 96 such 2×2 medjig squares. In the examples below, each 2×2 medjig square is made by combining different orientations of a single medjig piece.
0132
2310
3201
1023

0123
3210
1032
2301

0231
3102
0231
3102

We can use the 2×2 medjig squares to construct larger even ordered medjig squares. One possible approach is to simply combine the 2×2 medjig squares together. Another possibility is to wrap a smaller medjig square core with a medjig border. The pieces of a 2×2 medjig square can form the corner pieces of the border. An example of an 8×8 magic square is constructed below by combining four copies of the left most 2×2 medjig square given above:
114415
811510
132163
12796

01320132
23102310
32013201
10231023
01320132
23102310
32013201
10231023

11762464206347
3349301436523115
5640112753371026
24843592154258
1329503416325135
45611824864193
60447235741622
281239552593854

The next example is constructed by bordering a 2×2 medjig square core.
114415
811510
132163
12796

01012332
23321010
03023103
12310212
21023121
30310230
32012301
10321023

117143036526347
334962462043115
856114353211058
244059275372642
452923464323519
611350181648513
60447234157622
281255392593854

Singly even squares: Medjig square of order 1 does not exist. As such, the smallest odd ordered medjig square is of order 3, with magic constant 9. There are only 7 ways of partitioning the integer 9, our magic constant, into three parts. If these three parts correspond to three of the medjig pieces in a row, column or diagonal, then the relevant partitions for us are
Once a 3×3 medjig square has been constructed, we can convert it into a 6×6 magic square. For example:
816
357
492

230202
103131
311220
020331
322002
013113

2635119624
17828103315
30121423257
3215323416
3122279220
41336181129

There are 1,740,800 such 3×3 medjig squares. An easy approach to construct higher order odd medjig square is by wrapping an smaller odd ordered medjig square with a medjig border, just as with even ordered medjig squares. Approaches such as the LUX method can also be used. In the example below, a 5×5 medjig square is created by wrapping a medjig border around a 3×3 medjig square given previously:
17241815
23571416
46132022
101219213
11182529

0131013132
2302230210
3023020212
1210313130
0231122013
1302033102
3032200212
1201311330
3213101301
1020322023

1742994912683339065
6792247451768584015
9823558075714644166
4873305823289399116
4548131386370204797
2979656138895452272
85108762691921712853
3560123794444696783
86614393100252777934
3611681895755225984

Solving partially completed magic squares

Similar to the Sudoku and KenKen puzzles, solving partially completed magic squares has become a popular mathematical puzzle. Puzzle-solving centers on analyzing the initial given values and possible values of the empty squares. One or more solutions arise as the participant uses logic and permutation group theory to rule out all unsuitable number combinations.

Variations of the magic square

Extra constraints

Certain extra restrictions can be imposed on magic squares.
If raising each number to the nth power yields another magic square, the result is a bimagic, a trimagic, or, in general, a multimagic square.
A magic square in which the number of letters in the name of each number in the square generates another magic square is called an alphamagic square.
There are magic squares consisting entirely of primes. Rudolf Ondrejka discovered the following 3×3 magic square of primes, in this case nine Chen primes:
178971
113595
4729101

The Green–Tao theorem implies that there are arbitrarily large magic squares consisting of primes.
The following "reversible magic square" has a magic constant of 264 both upside down and right way up:
96118968
88699116
61861899
19986681

When the extra constraint is to display some date, especially a birth date, then such magic squares are called birthday magic square. An early instance of such birthday magic square was created by Srinivasa Ramanujan. He created a 4×4 square in which he entered his date of birth in DD-MM-CC-YY format in the top row and the magic happened with additions and subtractions of numbers in squares. Not only do the rows, columns, and diagonals add up to the same number, but the four corners, the four middle squares, the first and last rows two middle numbers, and the first and last columns two middle numbers all add up to the sum of 139.

Different constraints

Sometimes the rules for magic squares are relaxed, so that only the rows and columns but not necessarily the diagonals sum to the magic constant.
In heterosquares and antimagic squares, the 2n + 2 sums must all be different.

Multiplicative magic squares

Instead of adding the numbers in each row, column and diagonal, one can apply some other operation. For example, a multiplicative magic square has a constant product of numbers. A multiplicative magic square can be derived from an additive magic square by raising 2 to the power of each element, because the logarithm of the product of 2 numbers is the sum of logarithm of each. Alternatively, if any 3 numbers in a line are 2a, 2b and 2c, their product is 2a+b+c, which is constant if a+b+c is constant, as they would be if a, b and c were taken from ordinary magic square. For example, the original Lo-Shu magic square becomes:
165124
832128
256264

Other examples of multiplicative magic squares include:
2912
3661
3418

162056
402823
145244
128710

2750668413232
2452340547011
569204436656
5572911163630
4244560771226
1022483954863
787818403360

Multiplicative magic squares of complex numbers

Still using Ali Skalli's non iterative method, it is possible to produce an infinity of multiplicative magic squares of complex numbers belonging to set. On the example below, the real and imaginary parts are integer numbers, but they can also belong to the entire set of real numbers.
The product is: −352,507,340,640 − 400,599,719,520 i.
21+14i−70+30i−93−9i−105−217i16+50i4−14i14−8i
63−35i28+114i−14i2+6i3−11i211+357i−123−87i
31−15i13−13i−103+69i−261−213i49−49i−46+2i−6+2i
102−84i−28−14i43+247i−10−2i5+9i31−27i−77+91i
−22−6i7+7i8+14i50+20i−525−492i−28−42i−73+17i
54+68i138−165i−56−98i−63+35i4−8i2−4i70−53i
24+22i−46−16i6−4i17+20i110+160i84−189i42−14i

Additive-multiplicative magic and semimagic squares

Additive-multiplicative magic squares and semimagic squares satisfy properties of both ordinary and multiplicative magic squares and semimagic squares, respectively.

It is unknown if any additive-multiplicative magic squares smaller than 8×8 exist, but it has been proven that no 3×3 or 4×4 additive-multiplicative magic squares and no 3×3 additive-multiplicative semimagic squares exist.

Geometric magic squares

Magic squares may be constructed which contain geometric shapes instead of numbers. Such squares, known as geometric magic squares, were invented and named by Lee Sallows in 2001.
In the example shown the shapes appearing are two dimensional. It was Sallows' discovery that all magic squares are geometric, the numbers that appear in numerical magic squares can be interpreted as a shorthand notation which indicates the lengths of straight line segments that are the geometric 'shapes' occurring in the square. That is, numerical magic squares are that special case of a geometric magic square using one dimensional shapes.

Area magic squares

In 2017, following initial ideas of and , the first linear area magic square was constructed by Walter Trump.

Other magic shapes

Other two dimensional shapes than squares can be considered. The general case is to consider a design with N parts to be magic if the N parts are labeled with the numbers 1 through N and a number of identical sub-designs give the same sum. Examples include magic circles, magic rectangles, magic triangles magic stars, magic hexagons, magic diamonds. Going up in dimension results in magic spheres, magic cylinders, magic cubes, magic parallelepiped, magic solids, and other magic hypercubes.
Possible magic shapes are constrained by the number of equal-sized, equal-sum subsets of the chosen set of labels. For example, if one proposes to form a magic shape labeling the parts with, the sub-designs will have to be labeled with and.

Related problems

Over the years, many mathematicians, including Euler, Cayley and Benjamin Franklin have worked on magic squares, and discovered fascinating relations.

''n''-Queens problem

In 1992, Demirörs, Rafraf, and Tanik published a method for converting some magic squares into n-queens solutions, and vice versa.

Enumeration of magic squares

As mentioned above, the set of normal squares of order three constitutes a single equivalence class-all equivalent to the Lo Shu square. Thus there is basically just one normal magic square of order 3. But the number of distinct normal magic squares rapidly increases for higher orders. There are 880 distinct magic squares of order 4 and 275,305,224 of order 5. These squares are respectively displayed on 255 magic tori of order 4, and 251,449,712 of order 5. The number of magic tori and distinct normal squares is not yet known for any higher order.
Algorithms tend to only generate magic squares of a certain type or classification, making counting all possible magic squares quite difficult. Traditional counting methods have proven unsuccessful, statistical analysis using the Monte Carlo method has been applied. The basic principle applied to magic squares is to randomly generate n × n matrices of elements 1 to n2 and check if the result is a magic square. The probability that a randomly generated matrix of numbers is a magic square is then used to approximate the number of magic squares.
More intricate versions of the Monte Carlo method, such as the exchange Monte Carlo, and Monte Carlo backtracking have produced even more accurate estimations. Using these methods it has been shown that the probability of magic squares decreases rapidly as n increases. Using fitting functions give the curves seen to the right.

Magic squares in occultism

Magic squares of order 3 through 9, assigned to the seven planets, and described as means to attract the influence of planets and their angels during magical practices, can be found in several manuscripts all around Europe starting at least since the 15th century. Among the best known, the Liber de Angelis, a magical handbook written around 1440, is included in Cambridge Univ. Lib. MS Dd.xi.45. The text of the Liber de Angelis is very close to that of De septem quadraturis planetarum seu quadrati magici, another handbook of planetary image magic contained in the Codex 793 of the Biblioteka Jagiellońska. The magical operations involve engraving the appropriate square on a plate made with the metal assigned to the corresponding planet, as well as performing a variety of rituals. For instance, the 3×3 square, that belongs to Saturn, has to be inscribed on a lead plate. It will, in particular, help women during a difficult childbirth.
In about 1510 Heinrich Cornelius Agrippa wrote De Occulta Philosophia, drawing on the Hermetic and magical works of Marsilio Ficino and Pico della Mirandola. In its 1531 edition, he expounded on the magical virtues of the seven magical squares of orders 3 to 9, each associated with one of the astrological planets, much in the same way as the older texts did. This book was very influential throughout Europe until the counter-reformation, and Agrippa's magic squares, sometimes called kameas, continue to be used within modern ceremonial magic in much the same way as he first prescribed.
of Venus, drawn on the magic square of Venus. Each Hebrew letter provides a numerical value, giving the vertices of the sigil.
The most common use for these kameas is to provide a pattern upon which to construct the sigils of spirits, angels or demons; the letters of the entity's name are converted into numbers, and lines are traced through the pattern that these successive numbers make on the kamea.
In a magical context, the term magic square is also applied to a variety of word squares or number squares found in magical grimoires, including some that do not follow any obvious pattern, and even those with differing numbers of rows and columns. They are generally intended for use as talismans. For instance the following squares are: The Sator square, one of the most famous magic squares found in a number of grimoires including the Key of Solomon; a square "to overcome envy", from The Book of Power; and two squares from The Book of the Sacred Magic of Abramelin the Mage, the first to cause the illusion of a superb palace to appear, and the second to be worn on the head of a child during an angelic invocation:
SATOR
AREPO
TENET
OPERA
ROTAS

666848938
811544839
111383839
273774447

HESEB
EQAL
S
EG
B

ADAM
DARA
ARAD
MADA

Magic squares in popular culture

  • In Goethe's Faust, the witch's spell used to make a youth elixir for Faust, the , has been interpreted as a construction of a magic square.
  • The English composer Peter Maxwell Davies has used magic squares to structure many of his compositions. For example, his 1975 Ave Maris Stella uses the 9×9 magic square of Moon while his 1977 A Mirror of Whitening Light uses the 8×8 magic square of Mercury to create the entire set of notes and durations for the piece. His other works that employ magic squares include The Lighthouse, Resurrection, Strathclyde Concerto No. 3 for Horn and Trumpet, as well as many of his symphonies. According to Davies' own account:

A magic square in a musical composition is not a block of numbers – it is a generating principle, to be learned and known intimately, perceived inwardly as a multi-dimensional projection into that vast area of the internal ear – the space/time crucible – where music is conceived.... Projected onto the page, a magic square is a dead, black conglomeration of digits; tune in, and one hears a powerful, orbiting dynamo of musical images, glowing with numen and lumen.

  • Magic squares, including Benjamin Franklin's, appear as clues to the mystery in Katherine Neville's novels The Eight and The Fire.
  • Dürer's magic square and his Melencolia I both also played large roles in Dan Brown's 2009 novel, The Lost Symbol.
  • In the 2011 Korean television drama Deep Rooted Tree, King Sejong is shown attempting to construct a 33×33 magic square using lunch boxes. He ultimately discovers the "pyramid method" and completes the magic square with the help of an army of court attendants. This inspires him to create a more just form of government ruled by reason and words rather than military might.
  • On October 9, 2014 the post office of Macao in the People's Republic of China issued a series of stamps based on magic squares. The figure below shows the stamps featuring the nine magic squares chosen to be in this collection.
  • The metallic artifact at the center of The X-Files episode "Biogenesis" is alleged by Chuck Burks to be a magic square.
  • Mathematician Matt Parker attempted to create a 3x3 magic square using square numbers in a YouTube video on the Numberphile channel. His failed attempt is known as the Parker Square.
  • The first season Stargate Atlantis episode "Brotherhood" involves completing a magic square as part of a puzzle guarding a powerful Ancient artefact.
OWIKI.org. Text is available under the Creative Commons Attribution-ShareAlike License.