
n = 40;   % plots an n by n grid of quadrilateral cells
[xc,yc] = meshgrid(linspace(-1,1,n+1), linspace(-1,1,n+1));
[xp,yp] = mapc2p(xc,yc);
pcolor(xp,yp,0*yp);   
colormap([1 1 1]);   
daspect([1 1 1])

