Advertisement
Matthen

Julia Set

Jun 4th, 2011
282
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.37 KB | None | 0 0
  1. m = 50; Fract[z_] :=
  2. Module[{c = -0.8 + 0.156 I},
  3. Length[NestWhileList[#^2 + c &, z, Abs[#] < 10 &, 1, m]]];
  4. Colour[n_, z_] := {0.5 + 0.5 Tanh[(n - 30)/20],
  5. 0.5 + 0.3 Tanh[(n - 30)/20] + 0.1 (n/m),
  6. 0.2};
  7. \[Delta] = 0.01; range =
  8. Table[x + I y, {y, -2, 2, \[Delta]}, {x, -2, 2, \[Delta]}];
  9. data = Map[Colour[Fract[#], #] &, range, {2}];
  10. Image[data]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement