function wam=diskwam(n)
r=cos( pi/(2*n).*[0:2*n]);
X=cos( pi/(2*n).*[0:2*n]')*r;
Y=sin( pi/(2*n).*[0:2*n]')*r;
wam=[X(:),Y(:)];