
% plot quarter circle
n = 20;  
[xc,yc] = meshgrid(linspace(0,1,n+1), linspace(0,1,n+1));
[xp,yp] = mapc2p(xc,yc);
pcolor(xp,yp,0*yp);   
colormap([1 1 1]);   
daspect([1 1 1])
axis off

