%-------------------------------------------------------------------------------
% This program will take the file called 'data' gotten from the simulation run
% and produce a contour plot of the potentials along with the position of the
% cells.
%-------------------------------------------------------------------------------

a = load('data');
for i = 1:40
  for j = 1:40
    b(41-i,j) = a(((j-1)*40+i), 3);
  end
end
contour(b);
xlabel('x');
ylabel('y');

hold

% the next line in 'a' (1601) has the data for how many cells need to be depicted
cellnum = a(1601, 3);
 
% place cells by using plot3(x,y,min,'k*') where x,y are the indeces of the location
% and val is a constant location.  
for i = 1:cellnum
  plot(a(1601+i,1)+1,40-a(1601+i,2),'k*');
end

hold