function A = mapping(n)



% To create a Morton Scan order matrix



if n == 2

   A = [1 2; 3 4];

else

   B = mapping(n/2);

   A = [B B+(n/2)^2; B+(n/2)^2*2 B+(n/2)^2*3];

end