%-------------------------------------------------------------------------------
% This program will take the file called 'data' gotten from the simulation run
% and produce a surface plot of the potentials along with the position of the
% cells shown along the surface.
%-------------------------------------------------------------------------------
a = load('data');
for i = 1:40
for j = 1:40
b(41-i,j) = a(((j-1)*40+i), 3);
end
end
surf(b);
xlabel('x');
ylabel('y');
zlabel('potential');
% hold graph to place all cells
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,val,'k*') where x,y are
% the indeces of the location and val is a constant location.
for i = 1:cellnum
x = a(1601+i,1)+1;
y = 39-a(1601+i,2)+1;
plot3(x,y,b(y,x),'mh');
end
% release graph
hold