
% some potential here, just need to plot orbits of circles under action of deck group

 clear all ,close all, 
hold on
 
center=[0,0];
radius=1;  
THETA=linspace(0,2*pi,2000);
RHO=ones(1,2000)*radius;
[X,Y] = pol2cart(THETA,RHO);
X=X+center(1);
Y=Y+center(2); 
fill(X,Y,'r'),   axis square , grid on
 
axis([-1 1 -1 1]) 

    

k = .01; 
r = 1;

[xx,yy] = meshgrid(-r:k:r,-r:k:r);


x = 0:pi/4:2*pi;
P = sin(x);
Q = cos(x);

for j = 1:length(P)
for i = 1:length(P)

p = P(i);
q = Q(i);

h = (xx-p).^2 + (yy-q).^2;

h = h + repmat([p;q],1,length(h));

contour(xx,yy,h,.25/(max(1,5*(j-1))));
plot(p,q,'*','MarkerSize',10)

end
end

