close all;

n = 30;   % 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);

% We want the figure centered in a square, not a rectangle.
set(gcf,'paperposition',[0 0 6 6]);

% we want to be able to fill up the figure box, if necessary.
ha = axes('position',[0 0 1 1]);

pcolor(ha,xp,yp,0*yp);
colormap([1 1 1]);
daspect([1 1 1]);

% Finally, zoom to make the plot fill more of the figure box.
set(gca,'cameraviewanglemode','manual');
set(gca,'cameraviewangle',7);

axis off;
