Welcome to Fractal Forums

#version 3.7;

#global_settings {
  assumed_gamma 1.0
}

camera {
  location  <0, 2, 10>
  direction <0, 0, -2>
  up        <0, 1, 0>
  right     <1, 0, 0>
  sky       <0,1,0>
  look_at   <0, 1, 0>
}

// spline path through Mandelbrot set
// left number is spline parameter, i.e. tree height
// right vector is point in complex plane
#declare mandelPath = spline {
  natural_spline
  0.0 <  0.0,  0.0>
  1.0 < -2.0,  0.0>
}


#declare numPts = 4000;  // number of points to be rendered
#declare treeHeight = 2.0;  // height of treetop

#declare iterSkip = 100;  // do not display first 100 iterations

#declare curPt = 0;

union {  // group all spheres together

#while (curPt < numPts)
  #declare curT = curPt/(numPts - 1);     // current spline parameter
  #declare curHeight = treeHeight*curT;   // current point height
  #declare curC = mandelPath(curT);       // current constant C

  #declare Zr = 0;
  #declare Zi = 0;
  #declare Zr2 = 0;
  #declare Zi2 = 0;
  #declare iter = 0;  

  #declare numIters = iterSkip + curPt;

  #while (((Zr2 + Zi2) <= 4) & (iter < numIters))
    #declare iter = iter + 1;
    #declare Zi = 2*Zr*Zi + curC.v;
    #declare Zr = Zr2 - Zi2 + curC.u;
    #declare Zr2 = Zr*Zr;
    #declare Zi2 = Zi*Zi;
  #end

  #if (iter = numIters)
    sphere { <Zr, curHeight, Zi>, 0.025 }
  #end

  #declare curPt = curPt + 1;
#end

texture {
  pigment {rgb <0.5, 0.5, 0.5>}
  finish {
    ambient 0.0
    diffuse 0.8
    specular 0.35
    roughness 0.04
  } 
}

}  // close union around all spheres

// ground box
box {<-1, -1, -1>, <1, 0, 1>
  texture {
    pigment {rgb <0.5, 0.5, 0.5>}
    finish {
      ambient 0.0
      diffuse 0.8
      specular 0.35
      roughness 0.04
    } 
  }
}

light_source { <18, 15, 8> color rgb <.6, .7, .8>}
light_source { <-18, -10, -8> color rgb <.8, .7, .6>}

background { color rgb <0.3, 0.5, 0.8> }