What was being plotted at WWDC 2003?

At the WWDC 2003 keynote, Theodore Gray took the stage, along with Phill Schiller, and demoed Mathematica 5, comparing performance on G5 and Xeon processors.

In the presentation, a 40-step fractal battle is conducted and I'd like to recreate this on my new Macbook. So I was wondering what function is being used to make the DensityPlots (if that's what they are)?

• You don't, by chance, have any more pixels, do you? Commented Nov 21, 2019 at 20:44
• I looked but the best YouTube video is only 480p and the text is unreadable
– M.R.
Commented Nov 22, 2019 at 0:41
• Looks like it computing some high order polynomial roots
– M.R.
Commented Nov 22, 2019 at 2:43
• @jasonb ok I added a gif of the screen capture...
– M.R.
Commented Dec 3, 2019 at 7:23
• From Michael Trott: "I darkly remember making something for Theo. But I don't think I have it anymore." Commented Dec 3, 2019 at 15:34

Yes, that would be ListDensityPlot.

Some example code:

SqrtFractalSetup[{xmin_, xmax_}, {ymin_, ymax_}, n_] :=
With[{xminn = N[xmin], yminn = N[ymin],
xstep = N[(xmax - xmin)/n], ystep = N[(ymax - ymin)/n]},
re0 = Table[r, {r, xminn, xmax, xstep}, {i, yminn, ymax, ystep}];
im0 = Table[i, {r, xminn, xmax, xstep}, {i, yminn, ymax, ystep}]];

SqrtFractalDraw[c_, {xmin_, xmax_}, {ymin_, ymax_}, steps_] :=
Module[{quadrant}, {re, im} = {re0, im0};
Do[{re, im} = {(c - Sqrt[Abs[re]] + Sqrt[Abs[im]])^2,
(1 - Sqrt[Abs[im]] - Sqrt[Abs[re]])^2}; , {steps}];