Ed Pegg Jr
This picture starts with the greatest common denominator of the coordinates, then applies a Fourier transform. The complete code for this picture is as follows: With[{n = 512}, ListDensityPlot[ Abs[ Fourier[ Table[ GCD[i, j], {i,n}, {j, n}]]], ColorFunction -> Hue, Mesh -> False, Frame -> False]]

© Wolfram Research, Inc.