function [] = circle(a,r)
n=1000;
for k=1:n
    th=2*pi*(k-1)/n;
    x(k)=real(a)+r*cos(th);
    y(k)=imag(a)+r*sin(th);
end
%plot(x,y);
patch(x,y,'red')

end