|
Guest
|
http://www.geocities.com/rlbagulatftn/eyeofra_ifs.gif
The Eye of Ra fractal
by doing an affine inside the
kiss ellipse with reduced the number of transforms in
Mathematica:
Clear[f, dlst, pt, cr, ptlst, x, y]
RandomSeed[];
dlst = Table[ Random[Integer, {1, 2}], {n, 100000}];
f[1, {x_, y_}] := N[ {2*x*y/(x^2 + y^2) , (y^2 - x^2)/(y^2 + x^2)}];
f[2, {x_, y_}] := N[ {(2*((x - y)/
Sqrt[2]) - (x + y)/Sqrt[2])/(
2.83), (2*((x - y)/Sqrt[2]) + (x + y)/Sqrt[2])/(2.83)}];
pt = {0.5, 0.75};
cr[n_] := If[n - 2 == 0, RGBColor[
0, 0, 1], If[n - 3 == 0, RGBColor[0,
1, 0], If[n - 1 == 0, RGBColor[1, 0, 0], RGBColor[0, 0, 0]]]]
ptlst = Table[{cr[dlst[[j]]], Point[pt = f[dlst[[j]], Sequence[pt]]]},
{j, Length[dlst]}];
Show[Graphics[Join[{PointSize[.001]}, ptlst]],
AspectRatio -> Automatic, PlotRange -> All] |
|
|
