N = 20; 
xc = 0;
yc = 0;
R = 10;
[randXY] = generateRandomPupilSamplingPoints(N,xc,yc,R)  ;
[randXY2] = generateRandomPupilSamplingPoints(N) ; 
n=1;
cg1 = (circle_grid ( n, R, [xc,yc], N))';
size(cg1)

n=2;
cg2 = (circle_grid ( n, R, [xc,yc], N))';
size(cg2) 


n=3;
cg3 = (circle_grid ( n, R, [xc,yc], N))';
size(cg3)



