%Animated spinning circles
for i=1:30*60
   clf;
   t=[0:0.1:2*pi+0.1];
   for j = 1:12
       for k = 1:5
        R = (k/5)*(1 + sin((i+j)/30*2*pi))/2;   
        theta = 2*pi*(0.5*i/30 +j/12);    
        xpos = 3*cos(theta);
        ypos = 3*sin(theta); 
        x1 = xpos + R*cos(t);
        y1 = ypos + R*sin(t);
        plot(x1, y1)
        hold on;
        axis equal;
        axis([-5 5 -5 5]);    
       end;
   end; 
   pause(1/30);
end