For a nice introduction to roots and their fractal nature, see this post from 2011 on John Baez’s Azimuth blog or his 2023 UCR page here. Some of my code here is sourced from a great repository of interesting math visualizations from many others including the author of the pywonderland GitHub of Zhao Liang.
The SVG images below can be clicked to zoom in on another tab, but here is a Mathematica Notebook with Tooltip info on the selected roots.
The roots are complex valued (Real x-axis and Imaginary y-axis) with normalized Abs@binValue defining the color from 0 to 1. The ColorRules are as follows:
The analysis below provides more detail on the Roots-of-Unity Circle (RUC) at the center of the light blue area. These are interesting in that there are 4 quadrants each with 12 roots that all have 1/2< Abs@binVal <1 at the center of holes.
There are 6 roots that sit in holes surrounded by zeros (i.e. 2π/{12,10,8,6,5,4}) and 6 roots in holes with values <1/2 (i.e. 2π/{24, 14, 7, 11/2, 29/6 (75°), 437/90 (77°)}). The only other binVal roots >1/2 are on the 1 pixel wide x-axis with 5936 of 8000 (~75%) from the x-y grid of 8000×8000 (64 million) bins which sum over each of the 224 (~16 million) roots of polynomials up to with -/+1 coefficients. The rest of the x-axis roots are 2020 bins equal to 0 and only 44 bins between 0 and 1/2).
The upper-right quadrant of the SVG image has the RUC roots ringed in colors assigned by their order in each quadrant. To the right is an annotated band at the center of the black ring with the detail (below) with the angle at which the root is located. Click the image below to see this SVG and zoom further.
Below are 100 detail areas of “special” algorithmically searched roots with a 400×400 ArrayPlot and an 11×11 zoom of that Abs@binVal matrix with its ArrayPlot. This includes the 48 RUC roots and 12 x-axis roots with binVal>1/2, and 40 off-axis roots with binVal<1/2. This off-axis list includes one set of 4 quadrant roots that is a proper hole (in the 4 yellow rings located just above/below the bottom/top apex in the full image or 10-14th from the bottom in the detail list below). The other 9×4 roots are what I call “pseudo-holes” that have no enclosed high central binVal.
FYI – The color scheme on these detail zoom MatrixPlots are default to highlight numeric differences, so they don’t match the ColorRules in the large 8000×8000 MatrixPlot.
