Lsystem CylindricalPhyllotaxis;



  define N = 100;             · number of elements
  define h = 0.02;            · vertical distance between adjacent elements
  define r0 = 8;              · radius at the base
  define rf = 8;              · radius at the top
  define di = (r0-rf)/N;      · variation of radius

  colorspace = hsb;
  step = 0.3;
  angle = 1.0*fibangle;
  baseelement = {
    shape = quad;
  };

  useralphabet = {A,B};
  lifespan = N+1;

  axiom = @(40,1,0.9)A(r0);
  developmentrules = {
    p1: A(r) : true -> [&(90)f(sqrt(r))B]f(h)/A(r-di);
    p2: B    : true -> ^(90)F;
  };