function [X,Y,C]=CircleAM(n,c)
if nargin < 2, c=2; end
nu=pi/(2*acos(1/sqrt(c)));
m=ceil(nu*n);
A=cos(pi.*[0:2:2*m]'./(2*m))*exp(1i*pi/m.*[0:m-1]);
X=real(A);
Y=imag(A);
C=c;