

%---------------------------------------

% function to plot unit circles in the plane.
% this is called by dla2 to draw the dla disks
% for small aggregates.  The disks are colored
% red or black, depending on the value in 'list'

function plotcirc(z,list)

    figure;
    hold on;
    axis equal;
    axis off;
    t=[0:.1:2.1*pi];
    rad=1;
    for k=1:length(z)
       x=real(z(k));
       y=imag(z(k));
       if list(k)==0
          plot([x+rad*cos(t)],[y+rad*sin(t)],'k');
       else
          plot([x+rad*cos(t)],[y+rad*sin(t)],'r');
       end
    end % for k

return;

