07
Apr
06

### Magnet Fractals

These fractals were originally designed for predicting magnetic phase-transitions.

Here is some Mathematica code:
```Mandelbrot[zc_] := Length[FixedPointList[f[#, zc] &, 0, 100, SameTest -> (Abs[#] > 2 &)]]; Magnet[] := DensityPlot[Log[Mandelbrot[xc + I yc]], {xc, -1, 3}, {yc, -2, 2}, PlotPoints -> 275, Mesh -> False, Frame -> False, ColorFunction -> (If[# != 1, Hue[#], RGBColor[0, 0, 0]] &)];```

``` (* Magnet 1: runtime: 5 minutes *) f[z_, zc_] := (z^2 + zc - 1.0)^2 / (2 z + zc - 2)^2; Magnet[]; ```

```(* Magnet 2: runtime: 15 minutes *) f[z_, zc_] := (z^3 + 3 (zc - 1) z + (zc - 1) (zc - 2))^2 / (3 z^2 + 3 (zc - 2) z + (zc - 1) (zc - 2) + 1.0)^2; Magnet[]; ```
POV-Ray has a built-in function for these fractals:
```//Magnet 1: runtime: 0.5 second camera{orthographic location <1.25,0,-2> look_at <1.25,0,0> angle 90} plane{z,0 pigment{magnet 1 mandel 50 interior 1,1 color_map{[0 rgb 0][1/6 rgb <0,0,1>][1/3 rgb 1]}} finish{ambient 1}}```

``` ```

```//Magnet 2: runtime: 0 seconds camera{orthographic location <1,0,-2> look_at <1,0,0> angle 90} plane{z,0 pigment{magnet 2 mandel 50 color_map{[0 rgb 0][1/6 rgb <1,0,0>][1/3 rgb 1][1 rgb 1][1 rgb 0]}} finish{ambient 1}}```

## Welcome !

You will find here some of my favorite hobbies and interests, especially science and art.

I hope you enjoy it!

Subscribe to the RSS feed to stay informed when I publish something new here.

I would love to hear from you! Please feel free to send me an email : bugman123-at-gmail-dot-com