magic function

MATLAB magic function