MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

3d Mandelbulb Douady Julia

  • To: mathgroup at smc.vnet.net
  • Subject: [mg105858] 3d Mandelbulb Douady Julia
  • From: Roger Bagula <rlbagula at sbcglobal.net>
  • Date: Tue, 22 Dec 2009 04:04:21 -0500 (EST)
  • Reply-to: rlbagulatftn at yahoo.com

Resolution isn't much...
http://www.flickr.com/photos/fractalmusic/4203226808/
Clear[x, y, z, pc, n]
n = 30; norm[x_] := x.x; pc = {-0.122561, 0.744862, 0}; m0 = 3;
(* By Roger Bagula and Paul Nylander 21 Dec 2009©*)
Square[{x_, y_, z_}] := If[x == y == 0, {-z^2, 0, 0}, Module[{a = 1 -
z^2/(x^2 + y^2)}, {(x^2 - y^2)a, 2 x y a, -2 z Sqrt[x^2 + y^2]}]];
Julia3D[p0_] := Module[{p = p0, i = 0}, While[i <24 && norm[p] < 4, p =
N[Square[p] + pc]; i++]; i];

image1 = Table[z = 1.5; While[z >= -1.5 && Julia3D[{x, y, z}] < 24, z -=
m0/n]; If[z <= -1.51, {EdgeForm[], Cuboid[{0, 0,
0}, {m0/n, m0/n, m0/n}]}, {EdgeForm[], Cuboid[{x, y, z}, {x + m0/n,
y + m0/n, z + m0/n}]}], {y, -1.5, 1.5, m0/n}, {x, -1.5, 1.5, m0/n}];
image2 = Table[z = -1.5; While[z <= 1.5 && Julia3D[{
x, y, z}] < 24, z += m0/n]; If[z â?¥ 1.491.5, {EdgeForm[],
Cuboid[{0, 0, 0}, {m0/n,
m0/n, m0/n}]}, {EdgeForm[], Cuboid[{x, y, z}, {x + m0/n, y + m0/n, z +
m0/n}]}], {y, -1.5, 1.5, m0/n}, {x, -1.5, 1.5, m0/n}];
image3 = Table[y = 1.5; While[y >= -1.5 && Julia3D[{x, y, z}] < 24,
y -= m0/n]; If[y <= -1.51, {EdgeForm[],
Cuboid[{0, 0, 0}, {m0/n, m0/n, m0/n}]}, {
EdgeForm[], Cuboid[{x, y, z}, {x + m0/n, y +
m0/n, z + m0/n}]}], {z, -1.5, 1.5, m0/n}, {x, -1.5, 1.5, m0/n}];
image4 = Table[y = -1.5; While[y <= 1.5 && Julia3D[{x, y, z}] < 24,
y += m0/n]; If[y >= 1.49, {EdgeForm[], Cuboid[{0, 0, 0}, {m0/n,
m0/n, m0/n}]}, {
EdgeForm[], Cuboid[{x, y, z}, {
x + m0/n, y +
m0/n, z + m0/n}]}], {z, -1.5, 1.5, m0/n}, {x, -1.5, 1.5, m0/n}];

image0 = Join[image1, image2, image3, image4];

gr1 = Show[Graphics3D[image0, Boxed -> False, PlotRange -> All]]

The Mandelbulb by this method:
http://www.flickr.com/photos/fractalmusic/4198510584/in/photostream/
Respectfully, Roger L. Bagula
11759 Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 :
http://www.google.com/profiles/Roger.Bagula
alternative email: roger.bagula at gmail.com







  • Prev by Date: Re: intersection of two circles
  • Next by Date: Re: Marbles.nb [was Re: CellChangeTimes?]
  • Previous by thread: Export[..., ImageResolution, ... ImageSize]
  • Next by thread: problems with parameter lumping using ReplaceAll