function [x,y]=arcM(xC,yC,r,theta1,theta2,N)

% xC and yC are center of circle
% r is radius of circle


ang=linspace(theta1,theta2,N);%0:0.01:2*pi;
xp=r*cosd(ang);
yp=r*sind(ang);

x=xp+xC;
y=yp+yC;


end