% Run 10 simulations of the irreversible seeded accretion model.
% 1000 corners
% 10000 edge tiles

% For each,pick a random corner for each edge tile
edges = 10000;
corners = 1000;
simulations = 10;
arm_threshold = 5;

seeds = random('unid',corners,simulations,edges);
arm_sizes = zeros(corners,simulations);

for i=1:corners
  for j=1:simulations
    arm_sizes(i,j) = histc(seeds(j,:),i);
  end;
end;

%max_arm_size = max(max(arm_sizes));
%min_arm_size = min(min(arm_sizes));
min_arm_size = 0;
max_arm_size = 60;
bins = 21;

size(histc(arm_sizes(:,1),linspace(min_arm_size, ...
						  max_arm_size,bins)))

for i=1:simulations
  arm_sizes_hist(:,i) = histc(arm_sizes(:,i),linspace(min_arm_size, ...
						  max_arm_size,bins));
end;

arm_sizes_mean = mean(arm_sizes_hist,2);
arm_sizes_std = std(arm_sizes_hist,0,2);

figure(4);
hold off;
bar(linspace(min_arm_size,max_arm_size,bins),arm_sizes_mean,'y');
hold on;
errorbar(linspace(min_arm_size,max_arm_size,bins),arm_sizes_mean, ...
	 arm_sizes_std,'.');



axis([-2,62,0,800]);
xlabel('Corner Assembly Length, in Number of Tiles','FontSize',20);
ylabel('Number of Assemblies','FontSize',20);
title(['Seeded Irreversible Accretion 10:1'],'FontSize',20);
set(gca,'FontSize',20);
hold off;
saveas(gcf, 'irreversible_seeded_length_distribution', 'fig');
print -depsc2 irreversible_seeded_length_distribution.eps