Instead of an IFS this is a regular iterative calculation
of the Julia in 25 parts.
Clear[nz, f0, g, arraym, a0, b0, sd, p]
a0 = 1 + 5/9;
sd = Sqrt[7];
b0 = 4/9;
f0[r_] := p*(r^3 + (a0)*r^2 - (r/(a0) - 1))/sd + (1 - p)*(b0*r^3 + r^2)
3D Julia of Mandelbrot like cubic Bezier
(*Julia with SQRT(x^2 + y^2) limited measure*)
(*by R. L. BAGULA 19 July 2007 © *)
numberOfz2ToEscape[z_] := Block[
{escapeCount, nz = N[z], nzold = 0},
For[
escapeCount = 0,
(Sqrt[Re[nz]^2 + Im[nz]^2] < 16) && (escapeCount <
255) && (Abs[nz - nzold] > 10^(-3)),
nzold = nz;
nz = f0[nz];
++escapeCount
];
escapeCount
]
FractalPureM[{{ReMin_, ReMax_, ReSteps_},
{ImMin_, ImMax_, ImSteps_}}] :=
Table[
numberOfz2ToEscape[x + y I],
{y, ImMin, ImMax, (ImMax - ImMin)/ImSteps},
{x, ReMin, ReMax, (ReMax - ReMin)/ReSteps}
]
p = n/26;
Table[ListDensityPlot[FractalPureM[{{-3.5, 1.5, 100}, {-2.5, 2.5, 100}}],
Mesh -> False,
AspectRatio -> Automatic,
ColorFunction -> (Hue[2#] &)];, {n, 1, 25}]
Respectfully, Roger L. Bagula
11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814
:http://www.geocities.com/rlbagulatftn/Index.html
alternative email: rlbagula@[EMAIL PROTECTED]
|
