Fractal tree demo

[Nelson-numerical-software]

%==============================================================================
function demo_fractal_tree()
  fractaltree(10, 25)
end
%==============================================================================
function fractaltree(r, angle)
  totalIterations=r+1;
  f = gcf();
  c1=[0 0 0];
  len1=3*1.2^totalIterations;
  w1=r^0.6;
  plot([0 0],[-len1,0],'LineWidth',w1,'Color',c1)
  xlim([-2*len1,2*len1]);
  ylim([-2*len1,2*len1]);
  axis off
  hold on
  drawbranches(90,[0,0],totalIterations-1,angle,totalIterations);
end
%==============================================================================
function drawbranches(initAngle,pt,iterations,angle,totalIterations)
  len=1.2^(iterations);
  x1=pt(1);
  y1=pt(2);
  ang1=initAngle+angle;
  ang2=initAngle-angle;
  y2=len*sind(ang1)+y1;
  y3=len*sind(ang2)+y1;
  x2=len*cosd(ang1)+x1;
  x3=len*cosd(ang2)+x1;
  w=iterations^0.6;
  p1=plot([x1,x2],[y1,y2],'LineWidth',w);
  p2=plot([x1,x3],[y1,y3],'LineWidth',w);
  if iterations-1>0
    drawbranches(ang1,[x2,y2],iterations-1,angle,totalIterations);
    drawbranches(ang2,[x3,y3],iterations-1,angle,totalIterations);
  end
end
%==============================================================================

fractal_tree.mp4