function H=circle(center,radius)


N=1000;
THETA=linspace(0,2*pi,N);
RHO=ones(1,N)*radius;
[X,Y] = pol2cart(THETA,RHO);
X=X+center(1);
Y=Y+center(2);
H=plot(X,Y);
axis square;