function h = circle(x,y,r)
hold on
time = 0:pi/50:2*pi;
xunit = r * cos(time) + x;
yunit = r * sin(time) + y;
h = plot(xunit, yunit);
hold off