function circle(xc,yc,r,s)
% draw a circle of center (xc,yc)and radius r
theta = linspace(0,2*pi,150);
x = xc + r* cos(theta); y = yc +r*sin(theta);
plot(x,y,s)
end