function h=plotcircs(x,y,D)

th=0:.01:2*pi;

N=length(x);

plot(x(1)+D(1)/2*sin(th),y(1)+D(1)/2*cos(th)); 
hold on;
for n=2:N
  plot(x(n)+D(n)/2*sin(th),y(n)+D(n)/2*cos(th)); 
end
  
axis('equal');  
hold off;