arc(center,radius,t1,t2)

function []= arc(center,radius,t1,t2)
if(radius~=0)
NOP=round(200*radius*abs(t1-t2)/(2*pi));
else
NOP=1;
end
THETA=linspace(t1,t2,NOP);
RHO=ones(1,NOP)*radius;
[X,Y] = pol2cart(THETA,RHO);
X=X+center(1);
Y=Y+center(2);
plot(X,Y,'k');
axis square;