funcprot(0);
function circle(xc,yc)
    clf();
    ap=gca();
    ap.data_bounds=[-2,0;4,5];
    A=linspace(0,2*%pi,100)
    X=2*cos(A);
    Y=2*sin(A)+2;
    c=[0]
    cy=[2]
    plot(c,cy,'m+')
    plot(X,Y);
    
    a=linspace(0,2*%pi,100)
    x1=cos(a)+xc;
    y1=sin(a)+yc;
    plot(x1,y1);
endfunction

theta=linspace(0,2*%pi,100)
xa=cos(theta);
ya=sin(theta)+2;

for j=1:5
for i=1:100
    drawlater();
    circle(xa(i),ya(i));
    drawnow();
end
end

