% 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 = 1000;
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));
max_arm_size = 8;
min_arm_size = 0;
bins = 9;

size(histc(arm_sizes(:,1),0:max_arm_size))

for i=1:simulations
  arm_sizes_hist(:,i) = histc(arm_sizes(:,i),0:max_arm_size);
end;

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

figure(1);
hold off;
bar(linspace(min_arm_size,max_arm_size,bins),arm_sizes_mean,'y');
hold on;
axis([-1,9,0,1000]);
errorbar(linspace(min_arm_size,max_arm_size,bins),arm_sizes_mean, ...
	 arm_sizes_std,'.');
hold off;



xlabel('Corner Assembly Length, in Number of Tiles','FontSize',20);
ylabel('Number of Assemblies','FontSize',20);
title(['Seeded Irreversible Accretion 1:1'],'FontSize',20);
set(gca,'FontSize',20);
saveas(gcf, 'irreversible_seeded_length_distribution_one_to_one', 'fig');
print -depsc2 irreversible_seeded_length_distribution_one_to_one.eps