
% circlegrid.m
% Produces figure 1(a) of 
%  Wave Propagation Software, Computational Science, and Reproducible Research
%  by R. J. LeVeque, http://www.amath.washington.edu/~rjl/pubs/icm06

n = 30;
[xc,yc] = meshgrid(linspace(-1,1,n+1),linspace(-1,1,n+1));
[xp,yp] = mapc2p(xc,yc);

plot(xp,yp,'k')
hold on
plot(xp',yp','k')
hold off
axis square
axis([-1 1 -1 1])
axis off

%set(gcf,'paperunits','inches','paperposition',[0 0 3 3])
%print -djpeg90 circle
