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