clear all
clf

circle(0,0,1,'b')
hold on

circle(1,3,1.5,'g')

teta=0:pi/5:2*pi;

cx=cos(teta);
cy=sin(teta);

for i=1:length(cx)
   circle(cx(i),cy(i),3,'m')
   plot([0,cx(i)],[0,cy(i)],'k*',[0,cx(i)],[0,cy(i)],'k')
   pause(0.2)
end

