MATLAB Function Reference    

Magic square



M = magic(n) returns an n-by-n matrix constructed from the integers 1 through n^2 with equal row and column sums. The order n must be a scalar greater than or equal to 3.


A magic square, scaled by its magic sum, is doubly stochastic.


The magic square of order 3 is

This is called a magic square because the sum of the elements in each column is the same.

And the sum of the elements in each row, obtained by transposing twice, is the same.

This is also a special magic square because the diagonal elements have the same sum.

The value of the characteristic sum for a magic square of order n is

which, when n = 3, is 15.


There are three different algorithms:

To make this apparent, type:


If you supply n less than 3, magic returns either a nonmagic square, or else the degenerate magic squares 1 and [].

See Also

ones, rand

 luinc mat2str